Three-Level System: Electromagnetically Induced Transparency

Consider an atomic -system: a three-level system with two same-parity lower states |1 and |2, and an upper state |3 with opposite parity. A near-resonant light field consisting of two phase-coherent components of frequency a and b is applied. When the difference frequency ab=a-b is equal to the frequency splitting 12 of the two lower levels (the Raman resonance condition), the system is pumped into a coherent superposition of the two lower states which no longer absorbs the bichromatic field. When one frequency component of the light field is much stronger than the other, it is referred to as a "coupling" or "drive" field, and the reduction of the absorption of the weaker probe field is called "electromagnetically induced transparency". Other names used to refer to this situation are "dark resonance" and "coherent population trapping".
This loads the package.
In[1]:=
Click for copyable input
We define an atomic system consisting of two even-parity lower states and an odd-parity upper state. We apply a light field with components at frequencies a (near resonant with the |1→|3 transition) and b (near resonant with the |2→|3 transition).
Define the atomic system.
In[2]:=
Click for copyable input
Specify that there will be no explicit time dependence in the density matrix (steady-state condition).
In[3]:=
Click for copyable input
Out[3]=
Here is what the density matrix looks like.
In[4]:=
Click for copyable input
Out[4]//MatrixForm=
Define two optical fields with frequencies a and b.
In[5]:=
Click for copyable input
Out[5]=
In[6]:=
Click for copyable input
Out[6]=
The Hamiltonian for the system subject to the optical fields. Each field is assumed to interact with only one transition—the other terms are set to zero.
In[7]:=
Click for copyable input
Out[7]//MatrixForm=
The level diagram for the system.
In[8]:=
Click for copyable input
Out[8]=
Apply the rotating-wave approximation to the Hamiltonian.
In[9]:=
Click for copyable input
Out[9]//MatrixForm=
IntrinsicRelaxation and TransitRelaxation supply the relaxation matrices.
In[10]:=
Click for copyable input
Out[10]//MatrixForm=
OpticalRepopulation and TransitRepopulation supply the repopulation matrices.
In[11]:=
Click for copyable input
Out[11]//MatrixForm=
Here are the evolution equations.
In[12]:=
Click for copyable input
Out[12]//TableForm=
Observables supplies the absorption and phase shift experienced by the light in terms of the density-matrix elements.
In[13]:=
Click for copyable input
Out[13]=
Solve the system analytically. (We temporarily set the Method option of RowReduce to speed up Solve in Mathematica 9 and later.)
In[14]:=
Click for copyable input
Here is the probe absorption spectrum as a function of probe detuning.
In[15]:=
Click for copyable input
Out[15]=
The increased transmission at the Raman resonance is associated with coherence induced between the |1 and |2 states. Here is the coherence as a function of probe frequency.
In[16]:=
Click for copyable input
Out[16]=