The AC Stark Effect

A formula for the energy shifts due to an optical field is obtained for a two-level atomic system.
This loads the package.
In[1]:=
Click for copyable input
We define an atomic system consisting of two states (a ground state labeled 1 and upper state labeled 2). This is a "toy" system that neglects angular momentum (J and M are not defined). We apply a light field detuned from resonance by a frequency .
Define the atomic system.
In[2]:=
Click for copyable input
Define the optical field with frequency . This is an x polarized field by default, but for a toy system the polarization does not matter.
In[3]:=
Click for copyable input
Out[3]=
The Hamiltonian for the system subject to the optical field. For simplicity of notation, we define the "reduced" Rabi frequency R=||d1||E0.
In[4]:=
Click for copyable input
Out[4]//MatrixForm=
Assuming that 0>>R, we can apply the rotating-wave approximation to the Hamiltonian. We define the detuning =-0.
In[5]:=
Click for copyable input
Out[5]//MatrixForm=
The energy eigenvalues are found by diagonalizing the Hamiltonian.
The energy eigenvalues in the presence of the optical field.
In[6]:=
Click for copyable input
Out[6]=
Plot the eigenvalues as a function of Rabi frequency.
In[7]:=
Click for copyable input
Out[7]=
To find the Stark shifts, we subtract the eigenenergies for zero applied field. This gives one shift that is always positive, and one that is negative. (This is because the eigenvalue ordering changes as goes through zero.) To obtain the desired physical behavior, in which the energy shifts flip sign as goes through zero, we multiply by (equivalent to Sign[]).
Calculate the energy shifts.
In[8]:=
Click for copyable input
Out[8]=
A plot of the energy shift for the excited state as a function of detuning.
In[9]:=
Click for copyable input
Out[9]=
To lowest order, the shifts are quadratic.
In[10]:=
Click for copyable input
Out[10]=
Note that these solutions lose validity as →0. This is because we have neglected the natural line width of the upper state. When the detuning is on the order of or smaller than the natural line width, the formula must be modified. To analyze the general case, we can model the effect of the upper-state line width by making the upper-state energy complex.
Add an imaginary part to the upper-state energy.
In[11]:=
Click for copyable input
Out[11]//MatrixForm=
Find the eigenvalues of the modified Hamiltonian.
In[12]:=
Click for copyable input
Out[12]=
The energies are the real parts of the eigenvalues of the modified Hamiltonian. In order to find compact expressions for the real parts, it helps to do a little manipulation by hand.
Add /2 to each eigenvalue and take the real part.
In[13]:=
Click for copyable input
Out[13]=
The energies are now equal in magnitude and opposite in sign. Square each and simplify.
In[14]:=
Click for copyable input
Out[14]=
Reverse the procedure by taking the square root, restoring the opposite signs, and subtracting /2.
In[15]:=
Click for copyable input
Out[15]=
Again, we calculate the energy shifts by subtracting the zero-field energies and multiplying by Sign[].
In[16]:=
Click for copyable input
Out[16]=
The shift now turns around and goes through zero as the field goes on resonance.
In[17]:=
Click for copyable input
Out[17]=
A plot of the upper-state shift as a function of Rabi frequency and detuning.
In[18]:=
Click for copyable input
Out[18]=
To lowest order in the Rabi frequency the shifts are given by dispersive Lorentzian functions.
The series expansion.
In[19]:=
Click for copyable input
Out[19]=