Atomic Density Matrix TutorialMore About »

The Density Matrix and the Liouville Equation

The density matrix is a generalization of the state vector that allows graceful handling of ensemble averages and relaxation effects. For a single atom with state vector |, we define the density operator by =||.
Load the package.
In[1]:=
Click for copyable input
Define an atomic system and form its ket vector.
In[2]:=
Click for copyable input
Out[2]=
In[3]:=
Click for copyable input
Out[3]//MatrixForm=
Taking the outer product of the ket with its dual gives the density matrix in terms of the state vector coefficients.
In[4]:=
Click for copyable input
Out[4]//MatrixForm=
Evidently, the diagonal matrix elements of the density operator are the probabilities of the corresponding states (here the states 1 and 2). The off-diagonal matrix elements are known as coherences.
We can retrace the example given in "The State Vector and the Schrödinger Equation" in terms of the density-matrix elements rather than the state-vector amplitudes. The equation governing the evolution of the density matrix, , known as the Liouville equation, can be found directly from the Schrödinger equation.
The function DensityMatrix forms the density matrix for an atomic system.
In[5]:=
Click for copyable input
Out[5]//MatrixForm=
The Hamiltonian describing the interaction with an electric field. We normalize the electric field to frequency units by dividing by the reduced dipole matrix element between states 1 and 2.
In[6]:=
Click for copyable input
Out[6]//MatrixForm=
The Liouville equation.
In[7]:=
Click for copyable input
Out[7]//TableForm=
Initial conditions.
In[8]:=
Click for copyable input
Out[8]=
List of the density-matrix variables.
In[9]:=
Click for copyable input
Out[9]=
Solve the Liouville equation.
In[10]:=
Click for copyable input
Out[10]=
The diagonal density-matrix elements are the probabilities of finding each state.
In[11]:=
Click for copyable input
Out[11]=
If instead of a single atom, we are interested in averages over an ensemble of n atoms, each with state vector |i, we simply work with the averaged density matrix: =|ii|.
Everything then follows identically as above. Note that averaging the state vectors, rather than the density matrices, would not produce satisfactory results. Also, the averaged density matrix can represent states of the ensemble that can not correspond to a single atomic state vector. For example, the density matrix represents atoms that are equally likely to be in state 1 or 2, but have no well-defined phase between the amplitudes to be in these states, whereas a state vector always has a particular phase between its components.