Nonlinear Magneto-Optical Rotation in an Alkali Atom

Here we examine the interaction of linearly polarized light and a longitudinal magnetic field with an atomic sample, as in the tutorial "Zeeman Structure: Nonlinear Magneto-Optical Rotation". In this case, we consider the D1 line of the 87 Rb atom, and write the input parameters in units appropriate for comparison with experimental results.
Load the package.
In[1]:=
Click for copyable input
Set options for DensityMatrix and plotting function.
In[2]:=
Click for copyable input
In[3]:=
Click for copyable input
First we set up the atomic system. We could enter the various atomic-state parameters by hand, but in this case it is convenient to pull them from a database. In this example we will focus on transitions originating from the ground-state F=2 hyperfine sublevel, and so we truncate the system by deleting the ground-state F=1 sublevels.
Choose a transition. AtomicTransition just returns the names of the states associated with the given transition, which we can use to find the associated data in the AtomicData database.
In[4]:=
Click for copyable input
Out[4]=
Get the list of quantum numbers and other information that we will use to set up the atomic system. We could also include the additional numerical atomic parameters in this list, but it can be more convenient to substitute their values in later.
In[5]:=
Click for copyable input
Out[5]=
Create the list of sublevels in the atomic system. We specify BranchingRatio[1]→1, which means that the upper J state decays with 100% probability to the lower J state. The individual F→F branching ratios are calculated automatically.
In[6]:=
Click for copyable input
DeleteStates removes the states satisfying the given criterion from the list of atomic states.
In[7]:=
Click for copyable input
Out[7]=
Here are the remaining numerical parameters for the atomic system that we will substitute in later.
In[8]:=
Click for copyable input
Out[8]=
Now we form the Hamiltonian, relaxation and repopulation matrices, and Liouville equation assuming the steady-state condition.
Form the Hamiltonian assuming an x-polarized optical field with frequency and reduced Rabi frequency R, and a z-directed magnetic field with nominal Larmor frequency L.
In[9]:=
Click for copyable input
Out[9]//MatrixForm=
Draw the level diagram for the truncated system, showing the optical and magnetic-field interactions.
In[10]:=
Click for copyable input
Out[10]=
Apply the rotating-wave approximation to the Hamiltonian, and write the optical frequency in terms of the detuning .
In[11]:=
Click for copyable input
Out[11]//MatrixForm=
Create the matrix describing relaxation due to spontaneous decay and atomic transit.
In[12]:=
Click for copyable input
Out[12]//MatrixForm=
Create the matrix describing repopulation due to spontaneous decay and atomic transit.
In[13]:=
Click for copyable input
Out[13]//MatrixForm=
Form the Liouville equation for the steady-state density matrix.
In[14]:=
Click for copyable input
Out[14]//Short=
The list of density-matrix variables.
In[15]:=
Click for copyable input
Out[15]=
Now we find the changes in light polarization parameters in terms of the atomic density-matrix elements.
The absorption, phase-shift, optical rotation, and change of ellipticity per unit path length and unit density.
In[16]:=
Click for copyable input
Out[16]=
Now we write the various input parameters in terms of experimentally relevant quantities.
ExpandDipoleRME expresses the reduced dipole matrix element in terms of the natural width and energy of the transition.
In[17]:=
Click for copyable input
Out[17]=
Convert ADM units to cgs by multiplying by , and then to SI by multiplying by .
In[18]:=
Click for copyable input
Out[18]=
Write the optical electric field (in V/m) in terms of the light intensity (in mW/cm2).
In[19]:=
Click for copyable input
Out[19]=
Write the Rabi frequency in terms of the light intensity.
In[20]:=
Click for copyable input
Out[20]=
Write the Larmor frequency in terms of the magnetic field in Gauss.
In[21]:=
Click for copyable input
Out[21]=
Write the total absorption and optical rotation in terms of the path length (in cm), and the atomic density (in cm-3).
In[22]:=
Click for copyable input
Out[22]=
Substitute all of the above relations and the atomic parameters into the Liouville equation.
In[23]:=
Click for copyable input
Out[23]//Short=
Finally, we plot the fractional electric-field amplitude absorption, and optical rotation (in radians). We can incorporate the effect of Doppler broadening by treating the Doppler shift as a shift in the light frequency. Thus the Doppler-broadened signal is given by the convolution of the Doppler-free signal and the Doppler profile.
Plot the Doppler-free absorption and rotation signals at low magnetic field and fairly high light intensity to see the nonlinear magneto-optical effect. Note that for the chosen set of parameters there is significant power broadening of the Doppler-free signal.
In[24]:=
Click for copyable input
Out[24]=
Convolve the signals with the Doppler distribution to obtain the Doppler-broadened signal.
In[25]:=
Click for copyable input
Out[25]=
Choose higher magnetic field and lower light intensity to see the linear magneto-optical effect.
In[26]:=
Click for copyable input
Out[26]=
Convolve the signal with the Doppler distribution.
In[27]:=
Click for copyable input
Out[27]=