Suppose we have two subsystems $A$ and $B$, and the hamilitonian $H$ for the system separates into a sum $H_A(t) + H_B(t)$ of two time dependent hamiltonians $H_A(t)$ and $H_B(t)$ that act only on subsystems $A$ and $B$ respectively.
Suppose the collective system begins in a non-entangled (product) state $\psi_o \equiv |a\rangle \otimes |b\rangle$ and the system evolves under the action of $H$ for some period of time so that the state is now some $\psi'$.
Will $\psi'$ also be a product state? If not wouldn't this be surprising that correlations would develop in non-interacting systems?
A short time $dt$ later the Schrodinger equation tells us that $$ |a\rangle \otimes |b\rangle \to |a\rangle \otimes |b\rangle + \frac{dt}{i\hbar} \left( H_A | a \rangle \otimes | b \rangle + | a \rangle \otimes H_B | b \rangle \right) $$ so that to first order in $dt$ we could write $$ |a\rangle \otimes |b\rangle \to \left( 1 + \frac{dt}{i\hbar} H_A \right) |a\rangle \otimes \left( 1 + \frac{dt}{i\hbar} H_B \right) |b\rangle $$ so that it appears that the state remains in a product state. I do not know how to show this for finite $dt$.
Could we add a term $-\frac{\Delta t^2}{\hbar^2}$ in the Dyson series, i.e. \begin{align} U(t_o,t) & \equiv \Pi_{n=0}^N \left(1 + \frac{\Delta t}{i\hbar}\left( H_A(t_o + n \Delta t ) + H_B(t_o + n \Delta t) \right) \right) \\ & \to \Pi_{n=0}^N \left[ \left(1 + \frac{\Delta t}{i\hbar} H_A(t_o + n \Delta t ) \right)\left(1 + \frac{\Delta t}{i\hbar} H_B(t_o + n \Delta t ) \right) \right] \equiv \tilde{U}(t_o,t) \end{align} where $\Delta t \equiv \frac{t-t_o}{N}$ so that $U=\tilde{U}$ in the limit $N \to \infty$?