Maxwell-Bloch Equations for a Two-Level System

Requires a Wolfram Notebook System
Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.
This Demonstration shows the solution of the Maxwell–Bloch equation for a two-level system. This simplified treatment models the passage of a laser through a gas tube containing rubidium vapor, commonly studied in quantum optics experiments. Here, the "atoms" are ideal two-level quantum mechanical systems, with a ground state and an excited state. The atoms are assumed to be in their ground state at the start of the experiment. A laser beam enters the cell at and exits at
m. Although 32 meters is much longer than in a typical experiment, the parameters have been chosen to follow the example in Fig. 2 of the tutorial article by Siddons [1]. Siddons uses a completely different numerical method, so this gives us something to check our results against.
Contributed by: Zachary H. Levine (June 13)
(Quantum Optics Group, National Institute of Standards and Technology)
After work by P. Siddons [1]
Open content licensed under CC BY-NC-SA
Details
Siddons [1] gives examples of the Maxwell–Bloch equations, some of which are replicated here. The Maxwell–Bloch equations for a two-level system are given in (S9), (S26), (S43) and (S53)–(S55), where S is for Siddons. Here are the equations in slightly simplified notation:
,
,
,
,
.
The following symbols are used in these equations:
, the speed of light,
, the vacuum permittivity,
, the number density of atoms,
, the carrier frequency,
, the dipole matrix element,
, the population decay rate,
, the coherence decay rate,
, the detuning.
In some cases, [1] gives the general form for indices , which is copied here, but these are specialized to 21 in the equation. Although the equations are given with vector symbols, we assume that
and
are aligned with the
axis and treat these as scalars.
In solving the equations, we can eliminate in favor of
, leading to
.
The units are balanced in this equation. The left side has dimensions . The right side has dimensions
,
which agree with those on the left side. The density matrix is dimensionless. The other two equations have units of .
Before solving these equations, we divide them by or
to get a dimensionless time and Rabi frequency. Figure 2 in [1] uses this time unit. The spatial unit will also be converted to a dimensionless unit.
Let
,
,
,
,
,
,
.
Here, is a dimensionless distance and
is a dimensionless retarded time. The population
and coherence
are dimensionless. Also,
is the dimensionless rate of decoherence,
is the dimensionless detuning and
is a dimensionless constant whose physical meaning is how much electric field is generated by the coherence:
.
There are three constants in these equations, namely ,
and
. If the decoherence is due to spontaneous emission, then
[1, p. 12].
Siddons only solves for , leaving one dimensionless parameter
. Siddons only solves for this in a single case. In this case,
is a real number and
is purely imaginary. Of course,
is real, which is true for any
. Hence, we may rewrite the equations as a set of three differential equations in
:
.
Physically, is the Rabi frequency and is proportional to the electric field,
is the upper state population of a two-level system and
is the coherence.
The equations above are solved subject to ,
and there is an incident electric field that obeys
with
. (See initialization code for parameters and [1, p. 13] for details on how these are chosen.)
Reference
[1] P. Siddons, "Light Propagation through Atomic Vapours," Journal of Physics B: Atomic, Molecular and Optical Physics, 47(9), 2014 093001. doi:10.1088/0953-4075/47/9/093001.
Snapshots
Permanent Citation