Atomic Density Matrix TutorialMore About »

The State Vector and the Schrödinger Equation

The main question that the AtomicDensityMatrix package is intended to address is
Given the initial state of an atomic system subject to light and/or other external fields, and relaxation mechanisms, what is the subsequent evolution of the system and what effect does this evolution have on the interacting light fields?
The evolution of a quantum-mechanical state vector | is governed by the Schrödinger equation i t|(t)=H|(t), where H is the Hamiltonian. We will always use units in which H has units of frequency and the reduced Planck's constant is equal to one.
As a simple example, we will find the evolution of system consisting of two opposite-parity states with no Zeeman or hyperfine substructure, separated by an energy and subject to a static electric field.
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=
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[4]:=
Click for copyable input
Out[4]//MatrixForm=
The Schrödinger equation.
In[5]:=
Click for copyable input
Out[5]=
Create the initial conditions.
In[6]:=
Click for copyable input
Out[6]=
A list of the dependent variables.
In[7]:=
Click for copyable input
Out[7]=
Solve the Schrödinger equation.
In[8]:=
Click for copyable input
Out[8]=
The real and imaginary parts of the elements of the state vector.
In[9]:=
Click for copyable input
Out[9]=
Plot the trajectory of the elements of the state vector in the complex plane.
In[10]:=
Click for copyable input
Out[10]=
The probability of finding each state as a function of time.
In[11]:=
Click for copyable input
Out[11]=
Plot the probabilities.
In[12]:=
Click for copyable input
Out[12]=
The state vector and Schrödinger equation describe the evolution of a single atom (or a completely polarized atomic ensemble) when there are no relaxation processes at work. However, we are often interested in the evolution of partially polarized or unpolarized ensembles, for which the experimental observables are averages over the ensemble. We also may need to include relaxation effects. For this we use a formalism that generalizes the state vector: the density matrix.