Rabi cycle explained

See also: Rabi problem. In physics, the Rabi cycle (or Rabi flop) is the cyclic behaviour of a two-level quantum system in the presence of an oscillatory driving field. A great variety of physical processes belonging to the areas of quantum computing, condensed matter, atomic and molecular physics, and nuclear and particle physics can be conveniently studied in terms of two-level quantum mechanical systems, and exhibit Rabi flopping when coupled to an optical driving field. The effect is important in quantum optics, magnetic resonance, and quantum computing, and is named after Isidor Isaac Rabi.

A two-level system is one that has two possible energy levels. One level is a ground state with lower energy, and the other is an excited state with higher energy. If the energy levels are not degenerate (i.e. don't have equal energies), the system can absorb or emit a quantum of energy and transition from the ground state to the excited state or vice versa. When an atom (or some other two-level system) is illuminated by a coherent beam of photons, it will cyclically absorb photons and emit them by stimulated emission. One such cycle is called a Rabi cycle, and the inverse of its duration is the Rabi frequency of the system. The effect can be modeled using the Jaynes–Cummings model and the Bloch vector formalism.

Mathematical description of spin flopping

One example of Rabi flopping is the spin flipping within a quantum system containing a spin-1/2 particle and an oscillating magnetic field. We split the magnetic field into a constant 'environment' field, and the oscillating part, so that our field looks like\mathbf = \mathbf_ + \mathbf_ = B_0 \mathbf + B_1 (\cos(\omega t)\mathbf + \sin(\omega t)\mathbf)where

B0

and

B1

are the strengths of the environment and the oscillating fields respectively, and

\omega

is the frequency at which the oscillating field oscillates. We can then write a Hamiltonian describing this field, yieldingH = -\vec \cdot \mathbf = \omega_0 S_z + \omega_1(\cos(\omega t)S_x + \sin(\omega t)S_y)where

\omega0=\gammaeB0/2m

,

\omega1=\gammaeB1/2m

, and

Sx,Sy,Sz

are the spin operators. The frequency

\omega1

is known as the Rabi frequency. We can substitute in their matrix forms to find the matrix representing the Hamiltonian:\beginH &= \omega_0 \frac \begin 1 & 0 \\ 0 & -1 \end + \omega_1 \left(\cos(\omega t)\frac \begin 0 & 1 \\ 1 & 0 \end + \sin(\omega t)\frac \begin 0 & -i \\ i & 0 \end\right) \\&= \begin \omega_0 & \omega_1 e^ \\ \omega_1 e^ & -\omega_0\end\endwhere we have used

\cos(\omegat)+i\sin(\omegat)=ei\omega

. This Hamiltonian is a function of time, meaning we cannot use the standard prescription of Schrödinger time evolution in quantum mechanics, where the time evolution operator is

U(t)=e-iHt/\hbar

, because this formula assume that the Hamiltonian is constant with respect to time.

The main strategy in solving this problem is to transform the Hamiltonian so that the time dependence is gone, solve the problem in this transformed frame, and then transform the results back to normal. This can be done by shifting the reference frame that we work in to match the rotating magnetic field. If we rotate along with the magnetic field, then from our point of view, the magnetic field is not rotating and appears constant. Therefore, in the rotating reference frame, both the magnetic field and the Hamiltonian are constant with respect to time.

We denote our spin-1/2 particle state to be

|\psi(t)\rangle=c+(t)|+\rangle+c-(t)|-\rangle

in the stationary reference frame, where

|+\rangle

and

|-\rangle

are spin up and spin down states respectively, and
2
|c
+(t)|

+

2
|c
-(t)|

=1

. We can transform this state to the rotating reference frame by using a rotation operatorR_z(\theta) = \begin e^ & 0 \\ 0 & e^ \endwhich rotates the state counterclockwise around the positive z-axis in state space, which may be visualized as a Bloch sphere. At a time

t

and a frequency

\omega

, the magnetic field will have precessed around by an angle

\omegat

. To transform

|\psi(t)\rangle

into the rotating reference frame, note that the stationary x and y-axes rotate clockwise from the point of view of the rotating reference frame. Because the

Rz(\theta)

operator rotates counterclockwise, we must negate the angle to produce the correct state in the rotating reference frame. Thus, the state becomes|\tilde(t)\rangle = R_z(-\omega t)|\psi(t)\rangle = c_+(t)e^|+\rangle + c_-(t)e^|-\rangleWe may rewrite the amplitudes so that|\tilde(t)\rangle = \alpha_+(t)|+\rangle + \alpha_-(t)|-\rangleThe time dependent Schrödinger equation in the stationary reference frame isi\hbar \frac|\psi(t)\rangle = H(t)|\psi(t)\rangleExpanding this using the matrix forms of the Hamiltonian and the state yieldsi\hbar\begin\frac \\ \frac\end = \frac \begin \omega_0 & \omega_1 e^ \\ \omega_1 e^ & -\omega_0 \end \begin c_+(t) \\ c_-(t) \endApplying the matrix and separating the components of the vector allows us to write two coupled differential equations as follows\begin& i\hbar \frac = \frac c_+(t) + \frac e^ c_-(t) \\& i\hbar \frac = \frac e^ c_+(t) - \frac c_-(t)\endTo transform this into the rotating reference frame, we may use the fact that

c+(t)=\alpha+(t)e-i

and

c-(t)=\alpha-(t)ei

to write the following:\begin& i\hbar \frac = - \frac \alpha_+(t) + \frac \alpha_-(t) \\& i\hbar \frac = \frac \alpha_+(t) + \frac \alpha_-(t)\endwhere

\Delta\omega=\omega-\omega0

. Now define\tilde = \frac \begin -\Delta\omega & \omega_1 \\ \omega_1 & \Delta\omega \endWe now write these two new coupled differential equations back into the form of the Schrödinger equation:i\hbar \frac|\tilde(t)\rangle = \tilde|\tilde(t)\rangleIn some sense, this is a transformed Schrödinger equation in the rotating reference frame. Crucially, the Hamiltonian does not vary with respect to time, meaning in this reference frame, we can use the familiar solution to Schrödinger time evolution:|\tilde(t)\rangle = \tilde(t)|\tilde(0)\rangle = e^ |\tilde(0)\rangleThis transformed problem is equivalent to that of Larmor precession of a spin state, so we have solved the essence of Rabi flopping. The probability that a particle starting in the spin up state flips to the spin down state can be stated asP_ = \frac \sin^2\left(\fract\right)where

\Omega=\sqrt{\Delta\omega2+

2
\omega
1}
is the generalized Rabi Frequency. Something important to notice is that

P+

will not reach 1 unless

\Delta\omega=0

. In other words, the frequency of the rotating magnetic field

\omega

must match the environmental field's Larmor frequency

\omega0

in order for the spin to fully flip; they must achieve resonance. When resonance (i.e.

\omega=\omega0

) is achieved,

\Omega=\omega1

.

Within the rotating reference frame, when resonance is achieved, it is as if there is no environmental magnetic field, and the oscillating magnetic field looks constant. Thus both mathematically (as we have derived) and physically, the problem reduces to the precession of a spin state under a constant magnetic field (Larmor precession).

To transform the solved state back to the stationary reference frame, we reuse the rotation operator with the opposite angle

|\psi(t)\rangle=Rz(\omegat)|\tilde{\psi}(t)\rangle

, thus yielding a full solution to the problem.

Applications

The Rabi effect is important in quantum optics, magnetic resonance and quantum computing.

Quantum optics

See main article: Jaynes–Cummings model. Rabi flopping may be used to describe a two-level atom with an excited state and a ground state in an electromagnetic field with frequency tuned to the excitation energy. Using the spin-flipping formula but applying it to this system yields

2
|c
b(t)|

\propto\sin2(\omegat/2),

where

\omega

is the Rabi frequency.

Quantum computing

Any two-state quantum system can be used to model a qubit. Rabi flopping provides a physical way to allow for spin flips in a qubit system. At resonance, the transition probability is given by P_(t)=\sin^2\left(\frac\right).To go from state

|0\rang

to state

|1\rang

it is sufficient to adjust the time

t

during which the rotating field acts such that

\omega1t/2=\pi/2

or

t=\pi/\omega1

. This is called a

\pi

pulse. If a time intermediate between 0 and

\pi/\omega1

is chosen, we obtain a superposition of

|0\rang

and

|1\rang

. In particular for
t=\pi
2\omega1
, we have a
\pi
2
pulse, which acts as: |0\rang \to \frac
1\rang
The equations are essentially identical in the case of a two level atom in the field of a laser when the generally well satisfied rotating wave approximation is made, where

\hbar\omega0

is the energy difference between the two atomic levels,

\omega

is the frequency of laser wave and Rabi frequency

\omega1

is proportional to the product of the transition electric dipole moment of atom

\vec{d}

and electric field

\vec{E}

of the laser wave that is

\omega1\propto\hbar\vec{d}\vec{E}

. On a quantum computer, these oscillations are obtained by exposing qubits to periodic electric or magnetic fields during suitably adjusted time intervals.[1]

See also

References

Notes and References

  1. A Short Introduction to Quantum Information and Quantum Computation by Michel Le Bellac,