I am trying to understand better the rotating wave approximation (RWA).
Consider an atom modeled as a two level system, interacting with a Laser. I have the dipole momentum operator $$\vec{D} = d \left( \vec{\epsilon_d} \sigma_{-} + \vec{\epsilon_d}^{*} \sigma_{+}\right) \, .$$
The electric field is $$\vec{E} = E_0 \left( \vec{\epsilon} e^{j(\omega_L t + \phi_L)} + \vec{\epsilon}^{*} e^{-j(\omega_L t + \phi_L)} \right) \, .$$
We have the Hamiltonian of our system $$H=\frac{\hbar \omega_q}{2} \sigma_z - \vec{D} \cdot \vec{E}$$ where $\hbar \omega_q$ is the bare energy of our two level system.
After few calculations, we can write it as \begin{align} H &= \frac{\hbar \omega_L}{2}\sigma_z + \frac{\hbar (\omega_q - \omega_L)}{2}\sigma_z \\ &- dE_0 (\vec{\epsilon_d} \cdot \vec{\epsilon} e^{j(\omega_L t + \phi_L)} \sigma_{-} + \vec{\epsilon_d} \cdot \vec{\epsilon}^{*}e^{-j(\omega_L t + \phi_L)} \sigma_{-} \\ &+ \vec{\epsilon_d}^{*} \cdot \vec{\epsilon}e^{j(\omega_L t + \phi_L)} \sigma_{+} + \vec{\epsilon_d}^{*} \cdot \vec{\epsilon}^{*}e^{-j(\omega_L t + \phi_L)} \sigma_{+}) \, . \end{align}
The usual trick to show the RWA approximation is to go in interaction picture taking $\frac{\hbar \omega_L}{2}\sigma_z$ as the non interacting part. Doing that, we end up with \begin{align} H^I &= \frac{\hbar (\omega_q - \omega_L)}{2}\sigma_z \\ &-dE_0 (\vec{\epsilon_d} \cdot \vec{\epsilon}e^{j(\phi_L)} \sigma_{-} + \vec{\epsilon_d} \cdot \vec{\epsilon}^{*}e^{-j(2*\omega_L t + \phi_L)} \sigma_{-}\\ &+ \vec{\epsilon_d}^{*} \cdot \vec{\epsilon}e^{j(2*\omega_L t + \phi_L)} \sigma_{+} + \vec{\epsilon_d}^{*} \cdot \vec{\epsilon}^{*}e^{-j(\phi_L)} \sigma_{+}) \end{align}
And at this point we say that we will neglect the fast oscillating term.
Here, we had to go to interaction picture to see what we could neglect. Isn't it possible to directly see it on the first Hamiltonian $H$?
The R.W.A approximation is equivalent to directly neglect the terms $$\vec{\epsilon_d} \cdot \vec{\epsilon}^{*}e^{-j(\omega_L t + \phi_L)} \sigma_{-} $$ and $$\vec{\epsilon_d}^{*} \cdot \vec{\epsilon}e^{j(\omega_L t + \phi_L)} \sigma_{+}$$ in the first Hamiltonian $H$.
However, they are exactly of the same order as the two others that wont be neglected. So I am a little confused.
Is there a direct argument before using the interaction picture trick to see why we neglect those terms?