In Peskin and Schroeder's book p 83, the authors said
At any fixed time $t_0$, we can of course expand $\phi$ in terms of ladder operators: $$\phi(t_0, \textbf{x})=\displaystyle\int\frac{d^3p}{(2\pi)^3}\frac{1}{\sqrt{2E_{\mathbf{p}} }}(a_{\mathbf{p}} e^{i\textbf{p.x}}+a_{\mathbf{p}} ^\dagger e^{-i\textbf{p.x}}) (1) $$ ... $e^{iH_0 (t-t_0)} \phi (t_0, \mathbf{x}) e^{-iH_0 (t-t_0)} \equiv \phi_I(t,\mathbf{x}).$...
Since we can diagonize $H_0$, it is easy to construct $\phi_I$ explicitly:
$\phi_I(t, \mathbf{x} ) = \int \frac{d^3p}{ (2\pi)^3} \frac{1}{ \sqrt{2 E_{\mathbf{p}}} } \left( a_{\mathbf{p}} e^{-ipx} + a^{\dagger}_{\mathbf{p}} e^{ipx} \right) |_{x^0 = t - t_0} $ (4.15)
To derive (4.15), I thought about $$ e^{i H_0 (t-t_0)} a_{\mathbf{p}} e^{-i H_0 (t-t_0)} $$ use BCH expansion in terms of $e^{-B} A e^B$ to work out (4.15), similar to the free field theories.
As answered in Interacting Fields in QFT $a_{\mathbf{p}}$ is time-indenpendent for free theory (I thought about $(\partial^2 + m^2) \phi(x) =0 $ in Eq. (2) in the aforementioned link for free theories) but time-dependent on interacting theories. Will this difference affect commutation relation $[a^{int}_{\mathbf{p}},H_0]$? More specifically, $$ [a^{int}_{\mathbf{p}},H_0 = \int \frac{d^3p }{(2\pi)^3} \omega_p a_{\mathbf{p}}^{free, \dagger} a_{\mathbf{p}}^{free} ] =? $$