Rabi cycle

This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (December 2015) (Learn how and when to remove this message)

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. These two levels are a ground state with lower energy and an excited state with higher energy. If the energy levels are not degenerate (i.e. not having equal energies), the system can absorb a quantum of energy and transition from the ground state to the "excited" state. When an atom (or some other two-level system) is illuminated by a coherent beam of photons, it will cyclically absorb photons and re-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.

Rabi flopping refers to 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$${\displaystyle \mathbf {B} =\mathbf {B} _{env}+\mathbf {B} _{osc}=B_{0}\mathbf {k} +B_{1}(\cos(\omega t)\mathbf {i} +\sin(\omega t)\mathbf {j} )}$$where ${\displaystyle B_{0}}$ and ${\displaystyle B_{1}}$ are the strengths of the environment and the oscillating fields respectively, and ${\displaystyle \omega }$ is the frequency at which the oscillating field oscillates. We can then write a Hamiltonian describing this field, yielding$${\displaystyle H=-{\vec {\mu }}\cdot \mathbf {B} =\omega _{0}S_{z}+\omega _{1}(\cos(\omega t)S_{x}+\sin(\omega t)S_{y})}$$where ${\displaystyle \omega _{0}=\gamma eB_{0}/2m}$, ${\displaystyle \omega _{1}=\gamma eB_{1}/2m}$, and ${\displaystyle S_{x},S_{y},S_{z}}$ are the spin operators. The frequency ${\displaystyle \omega _{1}}$is known as the Rabi frequency. We can substitute in their matrix forms to find the matrix representing the Hamiltonian:$${\displaystyle {\begin{aligned}H&=\omega _{0}{\frac {\hbar }{2}}{\begin{bmatrix}1&0\\0&-1\end{bmatrix}}+\omega _{1}\left(\cos(\omega t){\frac {\hbar }{2}}{\begin{bmatrix}0&1\\1&0\end{bmatrix}}+\sin(\omega t){\frac {\hbar }{2}}{\begin{bmatrix}0&-i\\i&0\end{bmatrix}}\right)\\&={\begin{bmatrix}\omega _{0}&\omega _{1}e^{-i\omega t}\\\omega _{1}e^{i\omega t}&-\omega _{0}\end{bmatrix}}\end{aligned}}}$$where we have used ${\displaystyle \cos(\omega t)+i\sin(\omega t)=e^{i\omega t}}$. 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 ${\displaystyle 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 independence 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 ${\displaystyle |\psi (t)\rangle =c_{+}(t)|+\rangle +c_{-}(t)|-\rangle }$ in the stationary reference frame, where ${\displaystyle |+\rangle }$ and ${\displaystyle |-\rangle }$ are spin up and spin down states respectively, and ${\displaystyle |c_{+}(t)|^{2}+|c_{-}(t)|^{2}=1}$. We can transform this state to the rotating reference frame by using a rotation operator$${\displaystyle R_{z}(\theta )={\begin{bmatrix}e^{-i\theta /2}&0\\0&e^{i\theta /2}\end{bmatrix}}}$$which rotates the state counterclockwise around the positive z-axis in state space, which may be visualized as a Bloch sphere. At a time ${\displaystyle t}$ and a frequency ${\displaystyle \omega }$, the magnetic field will have precessed around by an angle ${\displaystyle \omega t}$. To transform ${\displaystyle |\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 ${\displaystyle R_{z}(\theta )}$ operator rotates counterclockwise, we must negate the angle to produce the correct state in the rotating reference frame. Thus, the state becomes$${\displaystyle |{\tilde {\psi }}(t)\rangle =R_{z}(-\omega t)|\psi (t)\rangle =c_{+}(t)e^{i\omega t/2}|+\rangle +c_{-}(t)e^{-i\omega t/2}|-\rangle }$$We may rewrite the amplitudes so that$${\displaystyle |{\tilde {\psi }}(t)\rangle =\alpha _{+}(t)|+\rangle +\alpha _{-}(t)|-\rangle }$$The time dependent Schrödinger equation in the stationary reference frame is$${\displaystyle i\hbar {\frac {d}{dt}}|\psi (t)\rangle =H(t)|\psi (t)\rangle }$$Expanding this using the matrix forms of the Hamiltonian and the state yields$${\displaystyle i\hbar {\begin{bmatrix}{\frac {dc_{+}}{dt}}\\{\frac {dc_{-}}{dt}}\end{bmatrix}}={\frac {\hbar }{2}}{\begin{bmatrix}\omega _{0}&\omega _{1}e^{-i\omega t}\\\omega _{1}e^{i\omega t}&-\omega _{0}\end{bmatrix}}{\begin{bmatrix}c_{+}(t)\\c_{-}(t)\end{bmatrix}}}$$Applying the matrix and separating the components of the vector allows us to write two coupled differential equations as follows$${\displaystyle {\begin{aligned}&i\hbar {\frac {dc_{+}}{dt}}={\frac {\hbar \omega _{0}}{2}}c_{+}(t)+{\frac {\hbar \omega _{1}}{2}}e^{-i\omega t}c_{-}(t)\\&i\hbar {\frac {dc_{-}}{dt}}={\frac {\hbar \omega _{1}}{2}}e^{i\omega t}c_{+}(t)-{\frac {\hbar \omega _{0}}{2}}c_{-}(t)\end{aligned}}}$$To transform this into the rotating reference frame, we may use the fact that ${\displaystyle c_{+}(t)=\alpha _{+}(t)e^{-i\omega t/2}}$ and ${\displaystyle c_{-}(t)=\alpha _{-}(t)e^{i\omega t/2}}$ to write the following:$${\displaystyle {\begin{aligned}&i\hbar {\frac {d\alpha _{+}}{dt}}=-{\frac {\hbar \Delta \omega }{2}}\alpha _{+}(t)+{\frac {\hbar \omega _{1}}{2}}\alpha _{-}(t)\\&i\hbar {\frac {d\alpha _{-}}{dt}}={\frac {\hbar \omega _{1}}{2}}\alpha _{+}(t)+{\frac {\hbar \Delta \omega }{2}}\alpha _{-}(t)\end{aligned}}}$$where ${\displaystyle \Delta \omega =\omega -\omega _{0}}$. Now define$${\displaystyle {\tilde {H}}={\frac {\hbar }{2}}{\begin{bmatrix}-\Delta \omega &\omega _{1}\\\omega _{1}&\Delta \omega \end{bmatrix}}}$$We now write these two new coupled differential equations back into the form of the Schrödinger equation:$${\displaystyle i\hbar {\frac {d}{dt}}|{\tilde {\psi }}(t)\rangle ={\tilde {H}}|{\tilde {\psi }}(t)\rangle }$$In 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:$${\displaystyle |{\tilde {\psi }}(t)\rangle ={\tilde {U}}(t)|{\tilde {\psi }}(0)\rangle =e^{-i{\tilde {H}}t/\hbar }|{\tilde {\psi }}(0)\rangle }$$This 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 as$${\displaystyle P_{+\to -}={\frac {\omega _{1}^{2}}{\Omega ^{2}}}\sin ^{2}\left({\frac {\Omega }{2}}t\right)}$$where ${\displaystyle \Omega }$ is the generalized Rabi Frequency. Something important to notice is that ${\displaystyle P_{+\to -}}$ will not reach 1 unless ${\displaystyle \Delta \omega =0}$. In other words, the frequency of the rotating magnetic field ${\displaystyle \omega }$ must match the environmental field's Larmor frequency ${\displaystyle \omega _{0}}$ in order for the spin to fully flip; they must achieve resonance. Note that when resonance is achieved, ${\displaystyle \Omega =\omega _{1}}$.

To transform the solved state back to the stationary reference frame, we reuse the rotation operator with the opposite angle${\displaystyle |\psi (t)\rangle =R_{z}(\omega t)|{\tilde {\psi }}(t)\rangle }$, thus yielding a full solution to the problem.

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

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

${\displaystyle |c_{b}(t)|^{2}\propto \sin ^{2}(\omega t/2),}$

where ${\displaystyle \omega }$ is the Rabi frequency.

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 $${\displaystyle P_{0\to 1}(t)=\sin ^{2}\left({\frac {\omega _{1}t}{2}}\right).}$$To go from state ${\displaystyle |0\rangle }$ to state ${\displaystyle |1\rangle }$ it is sufficient to adjust the time ${\displaystyle t}$ during which the rotating field acts such that ${\displaystyle \omega _{1}t/2=\pi /2}$ or ${\displaystyle t=\pi /\omega _{1}}$. This is called a ${\displaystyle \pi }$ pulse. If a time intermediate between 0 and ${\displaystyle \pi /\omega _{1}}$ is chosen, we obtain a superposition of ${\displaystyle |0\rangle }$ and ${\displaystyle |1\rangle }$. In particular for ${\displaystyle t={\frac {\pi }{2\omega _{1}}}}$, we have a ${\displaystyle {\frac {\pi }{2}}}$ pulse, which acts as: $${\displaystyle |0\rangle \to {\frac {|0\rangle +i|1\rangle }{\sqrt {2}}}}$$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 ${\displaystyle \hbar \omega _{0}}$ is the energy difference between the two atomic levels, ${\displaystyle \omega }$ is the frequency of laser wave and Rabi frequency ${\displaystyle \omega _{1}}$ is proportional to the product of the transition electric dipole moment of atom ${\displaystyle {\vec {d}}}$ and electric field ${\displaystyle {\vec {E}}}$ of the laser wave that is ${\displaystyle \omega _{1}\propto \hbar \ {\vec {d}}\cdot {\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]

• Atomic coherence
• Bloch sphere
• Laser pumping
• Optical pumping
• Rabi problem
• Vacuum Rabi oscillation
• Neutral particle oscillation

• Quantum Mechanics Volume 1 by C. Cohen-Tannoudji, Bernard Diu, Frank Laloe, ISBN 9780471164333
• A Short Introduction to Quantum Information and Quantum Computation by Michel Le Bellac, ISBN 978-0521860567
• The Feynman Lectures on Physics, Volume III
• Modern Approach To Quantum Mechanics by John S Townsend, ISBN 9788130913148