I was going through Chapter 9 of Schwartz's QFT book and one of the results bothers me.
Suppose we have a complex scalar field theory, and we want to find the propagator associated to the complex scalar field $\phi$. In principle, there are two ways of going about finding propagators, we could take the time-ordered product and manually calculate it, like the book did previously for a real scalar field. In which case I understand why the propagator must be $\langle 0|T\{\phi^*(x)\phi(y)\}|0\rangle$, and why this coincides with the real scalar field propagator, which, in momentum space, is given by
$$ \frac{i}{p^2-m^2+i\epsilon}. $$ This propagates both $\phi$ and $\phi^*$, you cannot disentangle them.
Alternatively, we could do what the book does for massless and massive spin 1 particles, which is taking the momentum space equations of motion, inverting them, and manually adding $+i\epsilon$ to make it time ordered. This is where the issue lies, I understand that the equations of motion for both $\phi$ and $\phi^*$ are Klein Gordon equations, justifying the fact that the propagator must be $\sim1/(p^2-m^2+i\epsilon)$. My issue is that, through this method, I cannot see why they must be entangled, since we get independent equations for both $\phi$ and $\phi^*$.
EDIT: I think I understand now, the fact that the propagator obtained from each equation of motion is the same, means that at the level of free theories, you can't distinguish one from the other (clearly once interactions are turned on they couple with $A_\mu$ with opposite charges). Which in turn means you can't separate propagating one from propagating the other. If this is correct, then a follow-up would be what if I have multiple complex scalar fields? Then they'd be distinguished by their mass in the propagator, correct?