The Rotating-Wave Approximation

The Hamiltonian for a light field interacting with an atomic transition has matrix elements that oscillate at the optical frequency. However, by a suitable change of basis and neglecting terms corresponding to very far-off-resonant optical coupling, the time dependence can often be eliminated from the Hamiltonian, greatly simplifying calculations. In this tutorial, we discuss this approximation, known as the rotating-wave approximation.
This loads the package.
In[1]:=
Click for copyable input
As a simple example, we consider two atomic states coupled by a near-resonant monochromatic light field.
A two-level atomic system.
In[2]:=
Click for copyable input
Out[2]=
The Hamiltonian for this system interacting with a monochromatic optical field with frequency and induced Rabi frequency R.
In[3]:=
Click for copyable input
Out[3]//MatrixForm=
If we write the matrix elements of the light-atom interaction Hamiltonian in terms of exponentials, rather than trigonometric functions, each off-diagonal matrix element is composed of two terms.
The Hamiltonian expressed in terms of exponential functions.
In[4]:=
Click for copyable input
Out[4]//MatrixForm=
We can show that, for each of these matrix elements, one of the terms induces a near-resonant coupling, while the other induces a far-off-resonant coupling.
Consider one of the matrix elements written above:
    ,
where 0 is the lower state with energy zero, and 0 is the upper state with energy 0. For a near-resonant transition, the energy of the state H0 should be close to that of 0, i.e., zero.
First we consider the term t|0. Applying the energy operator (with =0), we have
    |0 = - t|0 + t|0
            = - t|0 + t(E|0)
            = - t|0 + t0|0
            = (0-) t|0
In other words, the energy of t|0 is 0-, which, if 0, is near zero, the energy of the state |0. Thus, if 0-, where is the line width of the upper state, the term 0| t|0 represents a resonant transition, for which energy is conserved. On the other hand, performing a calculation similar to the above, we find that the energy of - t|0 is 0+. If 0, then this energy is much greater than the energy of the lower state (0+>> for optical frequencies). The term 0|- t|0 represents a far-off-resonant transition. This sort of transition is highly suppressed, and can be neglected.
We can visualize the situation by drawing a level diagram, and indicating the resonant and off-resonant couplings.
Level diagram for the system, showing both resonant and off-resonant couplings.
In[5]:=
Click for copyable input
Out[5]=
One arrow, pointing up from the lower level, nearly matches the separation between the lower and upper states, and represents an optical field detuned by -0 from resonance. The other arrow points down from the lower level, and couples the lower level to a point far below resonance (by an amount +0) with the excited state.
In order to see how we might remove the optical frequencies from the Hamiltonian, we first consider a physically analogous situation — that of a spin 1/2 atomic system subject to an oscillating magnetic field.
Define a second system consisting of the Zeeman sublevels of a spin-1/2 state.
In[6]:=
Click for copyable input
Out[6]=
The Hamiltonian for this system subject to an oscillating magnetic field directed along x and a static magnetic field directed along z is formally identical to the optically coupled two-level system above (allowing for an arbitrary shift in the zero reference energy by adding a constant term to the diagonal).
Find the Hamiltonian for this system interacting with a magnetic field directed along x and oscillating with frequency .
In[7]:=
Click for copyable input
Out[7]//MatrixForm=
The x-directed oscillating magnetic field can be thought of as the sum of two counter-rotating magnetic fields, one rotating clockwise around the z-axis and the other rotating counterclockwise. Each term in the off-diagonal elements of the Hamiltonian corresponds to one or the other of the rotating fields.
Physically, we expect that if we transform into a frame rotating about the z-axis at frequency , the component of the magnetic field that the frame is co-rotating with will appear static, while the counter-rotating field will appear to rotate at twice the frequency it has in the non-rotating frame. The unitary transformation required to perform this task is given by the matrix for rotations about the z axis:
Matrix for rotation by an angle - t about the z-axis.
In[8]:=
Click for copyable input
Out[8]//MatrixForm=
There is a slight complication in applying the transform to the Hamiltonian, however, due to the time dependence of the unitary transformation. We would like the transformed Hamiltonian to have the properties of a Hamiltonian, i.e., to satisfy the Schrödinger equation | in the rotating frame. The fact that the rotating frame is non-inertial, however, means that the Schrödinger equation must be modified by a correction term, analogous to the fictitious force used in classical mechanics to account for the modifications that Newton's laws undergo in an accelerating frame. Alternatively, it is more convenient to define a modified effective Hamiltonian in the rotating frame, so that the Schrödinger equation is the same in all frames. If we define a modified effective Hamiltonian in the rotating frame by , it can be easily shown, starting from the Schrödinger equation in the stationary frame, that this results in the same form for the Schrödinger equation in both the stationary and rotating frames.
Find the effective Hamiltonian in the rotating frame.
In[9]:=
Click for copyable input
Out[9]//MatrixForm=
The function EffectiveHamiltonian can be used to give the same result.
In[10]:=
Click for copyable input
Out[10]//MatrixForm=
We see that transforming the Hamiltonian to the effective rotating-frame Hamiltonian has two effects: the oscillation frequencies in the off-diagonal elements are shifted by , and the energy splitting between the diagonal terms is shifted by the same amount. This second effect is a result of the correction term in the formula for the effective Hamiltonian.
In the optical case that we considered initially, the analogous unitary transformation does not have the same physical interpretation as taking the stationary frame to a rotating frame. Nonetheless, on the strength of the analogy, we continue to refer to the transformed frame as a rotating frame.
In constructing the unitary transformation matrix for the optical case, we can think of shifting the upper state down by a frequency in order to simultaneously shift the optical frequency in the resonant coupling term to zero. The unitary transformation is then written as U= t D, where D is a diagonal matrix with the desired energy shift for each atomic state on the diagonal:
The unitary matrix for transformation to the interaction frame in the optical case.
In[11]:=
Click for copyable input
Out[11]//MatrixForm=
The unitary matrix that we find by this method is equivalent to the one used in the spin-1/2 case, up to an overall phase factor. We then find the effective Hamiltonian in the rotating frame:
Apply EffectiveHamiltonian to the stationary frame and chosen unitary matrix.
In[12]:=
Click for copyable input
Out[12]//MatrixForm=
The level diagram helps us visualize what we have accomplished. The upper state has been shifted down by the optical frequency , so that it now is only 0- above the lower state. The frequencies of the two optical fields have been shifted by the same amount, so that the frequency of one of them is zero (the two ends of the arrow are at the same energy) and the frequency of the other is 2 . Note, however, that each field has the same detuning from resonance as it did in the stationary frame: 0- below resonance for the near-resonant field, and +0 below resonance for the off-resonant field.
Level diagram for the system in the rotating frame, showing both resonant and off-resonant couplings.
In[13]:=
Click for copyable input
Out[13]=
Finally, we note that, in a typical atomic system, the far-off-resonant terms (those that oscillate at twice the optical frequency in the effective Hamiltonian) lead to excitations that are so weak that the terms can be neglected to an excellent approximation. This is the rotating-wave approximation.
The function DropFastTerms sets terms that oscillate at twice the optical frequency to zero.
In[14]:=
Click for copyable input
Out[14]//MatrixForm=
In the case considered here, we have eliminated all time dependence from the Hamiltonian. Furthermore, when written in terms of the optical detuning , there is no longer any explicit dependence on either the optical or transition frequency.
Write the optical frequency in terms of detuning from resonance.
In[15]:=
Click for copyable input
Out[15]//MatrixForm=
All succeeding calculations can be performed using this approximate effective Hamiltonian, and the results can be transformed back into the stationary frame, if necessary.
A number of functions are available that perform various steps of the rotating-wave approximation, or all of the steps together. RotatingWaveTransformMatrix can be used to automatically create an appropriate transformation matrix for going to the rotating frame. RotatingWaveTransform can create the transform matrix and apply it to the Hamiltonian, and RotatingWaveApproximation can perform the transformation and drop the fast oscillating terms, thus implementing the RWA in one step.
EffectiveHamiltonian[h,u,t]apply the unitary transformation matrix u to the Hamiltonian h, assuming time variable t
DropFastTerms[h,,t]drop terms from h that oscillate at twice the frequency
RotatingWaveTransformMatrix[sys,]automatically (for simple systems) create a unitary matrix appropriate for performing the RWA for a atomic system sys subject to an optical field of frequency
RotatingWaveTransform[sys,h,]automatically create a unitary transformation to the rotating frame and apply it to a Hamiltonian h
RotatingWaveApproximation[sys,h,]automatically transform the Hamiltonian h to the rotating frame and drop fast-oscillating terms

Making the rotating-wave approximation.

For a simple system, RotatingWaveApproximation generates a unitary transformation matrix by assuming that states with Energy→0 are lower states that should remain unshifted, while states with energies that are not explicitly zero are upper states that should be shifted down by the optical frequency . The RWA is then applied using this transformation. RotatingWaveTransformMatrix can be used to return the transform matrix itself, while RotatingWaveTransform applies the transform matrix to the Hamiltonian, without dropping fast-oscillating terms.
Automatically apply the RWA to a Hamiltonian, assuming an optical field of frequency .
In[16]:=
Click for copyable input
Out[16]//MatrixForm=
Return the automatically generated transform matrix.
In[17]:=
Click for copyable input
Out[17]//MatrixForm=
Automatically generate and apply the transform to the Hamiltonian, without dropping fast-oscillating terms.
In[18]:=
Click for copyable input
Out[18]//MatrixForm=
For more complicated systems, possibly including multiple optical fields, this simple method of generating the transformation to the rotating frame does not work. In this case, more information can be supplied to RotatingWaveApproximation, allowing an appropriate transform matrix to be generated, or the transform matrix itself can be supplied directly. See the documentation for RotatingWaveApproximation for more information.