Rabi frequency

Consider two energy eigenstates of a quantum system with Hamiltonian ${\displaystyle {\hat {H}}_{0}}$ (for example, this could be the Hamiltonian of a particle in a ${\displaystyle {\frac {1}{r}}}$ potential, like the Hydrogen atom or the Alkali atoms):

$${\displaystyle {\begin{aligned}\psi _{1}(\mathbf {r} ,t)&=e^{-i\omega _{1}t}|1\rangle \\\psi _{2}(\mathbf {r} ,t)&=e^{-i\omega _{2}t}|2\rangle \end{aligned}}}$$

We want to consider the time dependent Hamiltonian

$${\displaystyle {\hat {\mathcal {H}}}={\hat {H}}_{0}+{\hat {V}}(t)}$$

where ${\displaystyle {\hat {V}}(t)=e\mathbf {r} \cdot \mathbf {E} _{0}\cos(\omega t)}$ is the potential of the electromagnetic field. Treating the potential as a perturbation, we can expect the eigenstates of the perturbed Hamiltonian to be some mixture of the eigenstates of the original Hamiltonian with time dependent coefficients:

$${\displaystyle \Psi (\mathbf {r} ,t)=c_{1}(t)e^{-i\omega _{1}t}|1\rangle +c_{2}(t)e^{-i\omega _{2}t}|2\rangle }$$

Plugging this into the time dependent Schrödinger equation

$${\displaystyle i\hbar {\frac {\partial \Psi (\mathbf {r} ,t)}{\partial t}}={\hat {\mathcal {H}}}\Psi (\mathbf {r} ,t)}$$

taking the inner product with each of ${\displaystyle e^{i\omega _{1}t}\langle 1|}$ and ${\displaystyle e^{i\omega _{2}t}\langle 2|}$, and using the orthogonality condition of eigenstates ${\displaystyle \langle i|j\rangle =\delta _{i,j}}$, we arrive at two equations in the coefficients ${\displaystyle c_{1}(t)}$ and ${\displaystyle c_{2}(t)}$:

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)&={\frac {c_{2}(t)\cos(\omega t)}{\hbar }}(\langle 1|e\mathbf {r} \cdot \mathbf {E} _{0}|1\rangle +e^{i\omega _{0}t}\langle 1|e\mathbf {r} \cdot \mathbf {E} _{0}|2\rangle )\\i{\dot {c}}_{2}(t)&={\frac {c_{1}(t)\cos(\omega t)}{\hbar }}(e^{-i\omega _{0}t}\langle 2|e\mathbf {r} \cdot \mathbf {E} _{0}|1\rangle +\langle 2|e\mathbf {r} \cdot \mathbf {E} _{0}|2\rangle )\end{aligned}}}$$

where ${\displaystyle \omega _{0}=\omega _{1}-\omega _{2}}$. The two terms in parentheses are dipole matrix elements dotted into the polarization vector of the electromagnetic field. In considering the spherically symmetric spatial eigenfunctions ${\displaystyle |i\rangle }$ of the Hydrogen atom potential, the diagonal matrix elements go to zero, leaving us with

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)&={\frac {c_{2}(t)\cos(\omega t)}{\hbar }}\langle 1|e\mathbf {r} \cdot \mathbf {E} _{0}|2\rangle e^{i\omega _{0}t}\\i{\dot {c}}_{2}(t)&={\frac {c_{1}(t)\cos(\omega t)}{\hbar }}\langle 2|e\mathbf {r} \cdot \mathbf {E} _{0}|1\rangle e^{-i\omega _{0}t}\end{aligned}}}$$

or

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)&=\Omega c_{2}(t)\cos(\omega t)e^{i\omega _{0}t}\\i{\dot {c}}_{2}(t)&=\Omega ^{*}c_{1}(t)\cos(\omega t)e^{-i\omega _{0}t}\end{aligned}}}$$

Here ${\displaystyle \Omega$ :=\Omega _{1,2}} , where ${\displaystyle \Omega _{i,j}={\frac {\langle i|e\mathbf {r} \cdot \mathbf {E} _{0}|j\rangle }{\hbar }}}$ is the Rabi Frequency.

In the numerator we have the transition dipole moment for the ${\displaystyle i\to j}$ transition, whose squared amplitude represents the strength of the interaction between the electromagnetic field and the atom, and ${\displaystyle \mathbf {E} _{0}={\hat {\epsilon }}E_{0}}$ is the vector electric field amplitude, which includes the polarization. The numerator has dimensions of energy, so dividing by ${\displaystyle \hbar }$ gives an angular frequency.

By analogy with a classical dipole, it is clear that an atom with a large dipole moment will be more susceptible to perturbation by an electric field. The dot product includes a factor of ${\displaystyle \cos \theta }$, where ${\displaystyle \theta }$ is the angle between the polarization of the light and the transition dipole moment. When they are parallel the interaction is strongest, when they are perpendicular there is no interaction at all.

If we rewrite the differential equations found above:

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)=\Omega c_{2}(t)\cos(\omega t)e^{i\omega _{0}t}&\to {\frac {\Omega c_{2}}{2}}(e^{i(\omega -\omega _{0})t}+e^{-i(\omega +\omega _{0})t})\\i{\dot {c}}_{2}(t)=\Omega ^{*}c_{1}(t)\cos(\omega t)e^{-i\omega _{0}t}&\to {\frac {\Omega ^{*}c_{1}}{2}}(e^{i(\omega +\omega _{0})t}+e^{-i(\omega -\omega _{0})t})\end{aligned}}}$$

and apply the Rotating-wave approximation, which assumes that ${\displaystyle \omega +\omega _{0}>>\omega -\omega _{0}}$, such that we can discard the high frequency oscillating terms, we have

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)&={\frac {\Omega c_{2}}{2}}e^{i\delta t}\\i{\dot {c}}_{2}(t)&={\frac {\Omega ^{*}c_{1}}{2}}e^{-i\delta t}\end{aligned}}}$$

where ${\displaystyle \delta =\omega -\omega _{0}}$ is called the detuning between the laser and the atomic frequencies.

We can solve these equations, assuming at time ${\displaystyle t=0}$ the atom is in ${\displaystyle |1\rangle }$ (i.e. ${\displaystyle c_{1}(0)=1}$) to find

$${\displaystyle |c_{2}(t)|^{2}={\frac {\Omega ^{2}\sin ^{2}{\bigg (}{\frac {{\sqrt {\Omega ^{2}+\delta ^{2}}}t}{2}}{\bigg )}}{\Omega ^{2}+\delta ^{2}}}}$$

This is the probability as a function of detuning and time of the population of state ${\displaystyle |2\rangle }$. A plot as a function of detuning and ramping the time from 0 to ${\displaystyle t={\frac {\pi }{\Omega }}}$ gives:

We see that for ${\displaystyle \delta =0}$ the population will oscillate between the two states at the Rabi frequency.

The quantity ${\displaystyle {\sqrt {\Omega ^{2}+\delta ^{2}}}}$ is commonly referred to as the "generalized Rabi frequency." For cases in which ${\displaystyle \delta \neq 0}$, Rabi flopping actually occurs at this frequency, where ${\displaystyle \delta }$ is the detuning, a measure of how far the light is off-resonance relative to the transition. For instance, examining the above animation at an offset frequency of ±1.73, one can see that during the 1/2 Rabi cycle (at resonance) shown during the animation, the oscillation instead undergoes one full cycle, thus at twice the (normal) Rabi frequency ${\displaystyle \Omega _{i,j}}$, just as predicted by this equation. Also note that as the incident light frequency shifts further from the transition frequency, the amplitude of the Rabi oscillation decreases, as is illustrated by the dashed envelope in the above plot.

Coherent Rabi oscillations may also be driven by two-photon transitions. In this case we consider a system with three atomic energy levels, ${\displaystyle |1\rangle }$, ${\displaystyle |i\rangle }$, and ${\displaystyle |2\rangle }$, where ${\displaystyle |i\rangle }$ represents a so-called intermediate state with corresponding frequency ${\displaystyle \omega _{i}}$, and an electromagnetic field with two frequency components:

$${\displaystyle {\hat {V}}(t)=e\mathbf {r} \cdot \mathbf {E} _{L1}\cos(\omega _{L1}t)+e\mathbf {r} \cdot \mathbf {E} _{L2}\cos(\omega _{L2}t)}$$

Now, ${\displaystyle \omega _{i}}$ may be much greater than both ${\displaystyle \omega _{1}}$ and ${\displaystyle \omega _{2}}$, or ${\displaystyle \omega _{2}>\omega _{i}>\omega _{1}}$, as illustrated in the figure on the right.

Two photon excitation schema. ${\displaystyle \omega _{i}>>\omega _{2}>\omega _{1}}$ is shown on the left, while ${\displaystyle \omega _{2}>\omega _{i}>\omega _{1}}$ is shown on the right. The vertical axis is the frequency (or energy) axis.

A two-photon transition is not the same as excitation from the ground to intermediate state, and then out of the intermediate state to the excited state. Instead, the atom absorbs two photons simultaneously and is promoted directly between the initial and final states. There are two necessary conditions for this two-photon process (also known as a Raman process), to be the dominant model of the light-matter interaction:

$${\displaystyle {\begin{aligned}\omega _{L2}+\omega _{L1}&=\omega _{1}-\omega _{2}\\\Delta =|\omega _{L1}-\omega _{1}|&>>0\end{aligned}}}$$

In words, the sum of the frequencies of the two photons must be on resonance with the transition between the initial and final states, and the individual frequencies of the photons must be detuned from the intermediate state to initial and final state transitions. If the latter condition is not met and ${\displaystyle \Delta \to 0}$, the dominant process will be one governed by rate equations in which the intermediate state is populated and stimulated and Spontaneous emission events from that state prevent the possibility of driving coherent oscillations between the initial and final states.

We may derive the two-photon Rabi frequency by returning to the equations

$${\displaystyle {\begin{aligned}i{\dot {c}}_{1}(t)&={\frac {\Omega _{1i}c_{2}}{2}}e^{i\Delta t}\\i{\dot {c}}_{i}(t)&={\frac {\Omega _{1i}^{*}c_{1}}{2}}e^{-i\Delta t}\end{aligned}}}$$

which now describe excitation between the ground and intermediate states. We know we have the solution

$${\displaystyle c_{i}(t)={\frac {\Omega _{1i}\sin {\bigg (}{\frac {{\tilde {\Omega }}_{1i}t}{2}}{\bigg )}}{{\tilde {\Omega }}_{1i}}}}$$

where ${\displaystyle {\tilde {\Omega }}_{1i}}$ is the generalized Rabi frequency for the transition from the initial to intermediate state. Similarly for the intermediate to final state transition we have the equations

$${\displaystyle {\begin{aligned}i{\dot {c}}_{i}(t)&={\frac {\Omega _{i2}c_{2}}{2}}e^{i\Delta t}\\i{\dot {c}}_{2}(t)&={\frac {\Omega _{i2}^{*}c_{i}}{2}}e^{-i\Delta t}\end{aligned}}}$$

Now we plug ${\displaystyle c_{i}(t)}$ into the above equation for ${\displaystyle {\dot {c}}_{2}(t)}$

$${\displaystyle i{\dot {c}}_{2}(t)={\frac {\Omega _{i2}^{*}\Omega _{1i}\sin {\bigg (}{\frac {{\tilde {\Omega }}_{1i}t}{2}}{\bigg )}}{2{\tilde {\Omega }}_{1i}}}e^{-i\Delta t}}$$

Such that, upon solving this equation, we find the coefficient to be proportional to:

$${\displaystyle c_{2}(t)\propto {\frac {\Omega _{i2}\Omega _{1i}}{2\Delta }}}$$

This is the effective or two-photon Rabi frequency.[1] It is the product of the individual Rabi frequencies for the ${\displaystyle |1\rangle \to |i\rangle }$ and ${\displaystyle |i\rangle \to |2\rangle }$ transitions, divided by the detuning from the intermediate state ${\displaystyle |i\rangle }$.

• Rabi cycle
• Vacuum Rabi oscillation
• Rabi resonance method