Laser Cooling of Atoms
In this problem you are asked to consider the mechanism of atom cooling with the help of laser radiation. Investigations in this field led to considerable progress in the understanding of the properties of quantum gases of cold atoms, and were awarded Nobel prizes in 1997 and 2001 .
Theoretical Introduction
Consider a simple two-level model of the atom, with ground state energy E_{g} and excited state energy E_{e}. Energy difference is E_{e}-E_{g}=\hbar\omega_{0}\,,, the angular frequency of used laser is w, and the laser detuning is δ =w-\omega_{0}\ll\omega_{0}. Assume that all atom velocities satisfy v « c, where c is the light speed. You can always restrict yourself to first nontrivial orders in small parameters \frac{\upsilon}{c} and \frac{\delta}{w_{0}}.Natural width of the excited state E_{e} due to spontaneous decay is \gamma« w_{0}: for an atom in an excited state, the probability to return to a ground state per unit time equals \gamma. When an atom returns to a ground state, it emits a photon of a frequency close to w_{0} in a random direction.
It can be shown in quantum mechanics, that when an atom is subject to low-intensity laser radiation, the probability to excite the atom per unit time depends on the frequency of radiation in the reference frame of the atom, w_{a} , according to
\gamma_{p}=s_{0}\,\frac{\frac{\gamma}{2}}{1+\frac{4(\omega_{a}-\omega_{0})^{2}}{\gamma^{2}}}\ll\gamma,where S_{0} « 1 is a parameter, which depends on the properties of atoms and laser intensity.
In this problem properties of the gas of sodium atoms are investigated neglecting the interactions between the atoms. The laser intensity is small enough, so that the number of atoms in the excited state is always much smaller than number of atoms in the ground state. You can also neglect the effects of the gravitation, which are compensated in real experiments by an additional magnetic field.
Numerical values:
Planck constant \hbar=1.05 \cdot 10^{-34}~{\mathrm{Js}}
Boltzmann constant k_{B}=1.38\cdot 10^{-23}~ \mathsf{J K}^{-1}
Mass of sodium atom m=3.81~ \cdot ~10^{-26}\,\mathrm{kg}
Frequency of used transition \omega_{0}=2\pi~\cdot~5.08~\cdot ~10^{14}~\mathrm{Hz}
Excited state linewidth \gamma=2\pi~\cdot~9.80~\cdot~10^{6}~\operatorname{Hz}
Concentration of the atoms n=10^{14} ~cm^{-3}
Questions
Questions
(a) Suppose the atom is moving in the positive x direction with the velocity v_{x} ,and the laser radiation with frequency ω is propagating in the negative x direction. What is the frequency of radiation in the reference frame of the atom?
(b) Suppose the atom is moving in the positive x direction with the velocity v_{x} ,and two identical laser beams shine along x direction from different sides. Laser frequencies are ω, and intensity parameters are s_{0} .
Find the expression for the average force F ( v_{x}) acting on an atom. For small v_{x} this force can be written as F(v_{x}) =-βv_{x} . Find the expression for β. What is the sign of \delta= ω-ω_o' if the absolute value of the velocity of the atom decreases? Assume that momentum of an atom is much larger than the momentum of a photon.
In what follows we will assume that the atom velocity is small enough so that one can use the linear expression for the average force.
(c) If one uses 6 lasers along x, y and z axes in positive and negative directions, then for β> 0 the dissipative force acts on the atoms, and their average energy decreases. This means that the temperature of the gas, which is defined through the average energy, decreases. Using the concentration of the atoms given above, estimate numerically the temperature T_{Q}, for which one cannot consider atoms as point-like objects because of quantum effects.
In what follows we will assume that the temperature is much larger than To and six lasers along x, y and z directions are used, as was explained in part (c).
In part (b) you calculated the average force acting on the atom. However, because of the quantum nature of photons, in each absorption or emission process the momentum of the atom changes by some discrete value and in random direction, due to the recoil processes.
(d) Determine numerically the square value of the change of the momentum of the atom, (Δp)² , as the result of one absorption or emission event.
(e) Because of the recoil effect, average temperature of the gas after long time does not become an absolute zero, but reaches some finite value. The evolution of the momentum of the atom can be represented as a random walk in the momentum space with an average step{\sqrt{(\Delta\rho)^{2}}}, and a cooling due to the dissipative force. The steady-state temperature is determined by the combined effect of these two different processes. Show that the steady state temperature {T}_{d} is of the form: T_{d}\ =\ {\frac{\hbar\gamma\left(x+{\frac{1}{x}}\right)}{4k_{\mathrm{B}}}}. Determine x.
Assume that {T}_{d} is much larger than {\frac{(\Delta\rho)^{2}}{4k_{B}m}}.
Note: If vectors P_{1}\,,\;P_{2}\,,\;\cdots,\;P_{n} are mutually statistically uncorrelated, mean square value of their sum is
((P_{1}+P_{2}+\cdots+P_{n})^{2}\,\rangle=\langle P_{1}{}^{2}\,\rangle+\langle P_{2}{}^{2}\rangle\,+\cdot\cdot\cdot+\langle P_{n}{}^{2}\,\rangle.(f) Find numerically the minimal possible value of the temperature due to recoil effect. For what ratio \frac{\delta}{\gamma} is it achieved?
(a) \omega\!\left(1+{\frac{v_{x}}{c}}\right)\!, this is classic Doppler effect.
(b) Absolute value of the momentum, transferred during each absorption, equals
{\frac{\hbar\omega_{0}}{c}}. (1)
The momentum of the emitted photon is uniformly distributed over different directions, and after averaging gives a contribution which is much smaller than {\frac{\hbar\omega_{0}}{c}}. . The average force is nonzero, since for atoms moving towards c right frequency of right laser gets larger (due to Doppler effect discussed in part (A)), while frequency of left laser goes down. Since number of scattered photons depends on the frequency in the reference frame of the atom, there is a net nonzero force. It equals
F(v_{x})=F_{\mathrm{+}}+F_{-}=\frac{\hbar\omega_{0}}{c}.\frac{s_{0}\gamma}{2}.\left[\frac{1}{1+\frac{4\left(\delta+\frac{\omega_{0}V_{x}}{c}\right)^{2}}{\gamma^{2}}}-\frac{1}{1+\frac{4\left(\delta-\frac{\omega_{0}V_{x}}{c}\right)^{2}}{\gamma^{2}}}\right].
For \frac{v_{x}}{c}\ll\frac{\delta}{\omega_{0}},
\beta=-\,\frac{8\hbar\omega_{0}^{2}\;\delta s_{0}}{\gamma c^{2}\biggl[1+4\biggl(\frac{\delta}{\gamma}\biggr)^{2}\biggr]^{2}}.
For β > 0, one needs
\delta\lt 0.(c) Characteristic de-Broglie wavelength at temperature T equals \lambda ={\frac{\hbar}{\sqrt{m k_{B}T}}}. To consider the atoms as point-like objects one needs this distance to be much smaller than characteristic inter particle separation n^{- {\frac{1}{3}}}. From the condition that these two lengths are of the same order of magnitude we get
T_{Q}={\frac{\hbar^{2}n^{\frac{2}{3}}}{k_{B}m}}\approx10^{-6}~\mathrm{K}.
(\mathrm{d})\,\,\,\langle\Delta p^{2}\rangle\,=\,\frac{\hbar^{2}\omega_{0}^{2}}{c^{2}}\approx\,10^{-54}~\mathrm{kg}^{2}{\mathrm{m}}^{2}/{\mathrm{s}}^{2}\,– this is the mean square recoil momentum of a photon.
(e) Assume that the steady state value of the average square of the momentum of atom equals P_{0}^2. In steady state regime this quantity does not change with time, and temperature is obtained according to 3k_{B}\,\frac{T_{d}}{2}=\frac{P_{0}^{2}}{2m}.Let the momentum at some point of time in steady state regime be P_{0} . Let us consider the value of the momentum after some time t. During this time the atom will participate in N=6\gamma_{p}t\gg1 absorption-emission processes (6 comes from the number of lasers). For each absorption-emission event the atom gets two recoil kicks, each with a mean square value (Δp²) calculated in part d) (one kick is during absorption and one is during emission). The directions of these kicks are uncorrelated for different events, so this leads to an increase of the mean square of the momentum by 2N (Δp² ).
On the other hand, atoms are cooled because of the dissipative force, and the change of the mean square of the momentum because of this process is -{\frac{2\beta P_{0}^{2}t}{m}}. For steady state solution these two processes compensate each other, so we obtain:
P_{0}^{2}=12\langle\Delta p^{2}\rangle\gamma_{p}\,\frac{m}{2\beta}=\frac{3\hbar m\gamma\Bigl(\frac{2\,\mid\delta\,\mid}{\gamma}+\frac{\gamma}{2\,\mid\delta\,\mid}\Bigr)}{4}.Thus the temperature
T_{d}={\frac{\hbar\gamma\left({\frac{2\mid\delta\mid}{\gamma}}+{\frac{\gamma}{2\mid\delta\mid}}\right)}{4k_{B}}}.(f) The minimum is achieved for \delta\,=-\,{\frac{\gamma}{2}}, and equals \frac{\hbar \gamma}{2k_{B}} = 2.4~\cdot~10^{-4}~K.