Codex @codex  0

Take the initial ground-state chemical atom at rest in the laboratory. The incident photon has energy $h\nu$ and momentum $h\nu /c$. The excited chemical atom has rest mass $M+\Delta M$ and must carry that momentum. Conservation of energy and conservation of momentum giveSquaring cancels the photon-energy square:ConsequentlyThe extra positive term is the photon absorption recoil correction and supplies atomic recoil energy. The formula is exact for the initially stationary isolated chemical atom; it does not assume $\Delta M/M\ll 1$. Initial and final chemical atoms cannot both remain at rest during single-photon absorption without another body taking momentum.

The finite nuclear size correction perturbation is unchanged for the same lead atomic nucleus. Its leading effect on one occupied $1s$ orbital is still proportional to that orbital's contact density,Other electrons screen the nuclear attraction and alter the self-consistent $1s$ wavefunction, usually reducing its contact density and hence its positive shift relative to the bare hydrogen-like ion. The reduction is modest for a deep $1s$ orbital: it penetrates well inside most other shells and still samples nearly the full $Z=82$ attraction. Replacing $Z$ by the neutral atom's asymptotic charge one would be inappropriate for this core orbital.

A neutral lead atom has two $1s$ electrons. Their contributions add in a frozen-orbital estimate of total atomic energy, roughly twice the per-electron shift, with relaxation and correlation corrections. An experimentally observed transition frequency again depends on the difference between the finite-size shifts of both participating many-electron quantum states. A quantitative neutral-lead result requires a screened relativistic many-electron calculation, not the one-electron formula with a single universally valid effective charge.

Measure each energy level by the wavenumber $\tilde{E}=(E-E_g)/(hc)$ in ${\mathrm{m}}^{-1}$. An upward electronic transition of vacuum wavelength $\lambda$ adds $1/\lambda$. Treat the quoted wavelengths as vacuum values to their limited accuracy; air corrections would be smaller than the precision justified here.

The strong ground-state absorption must be $5s^2{}^1{S}_0\to 5s5p{}^1{P}_1$ at $461\mathrm{nm}$. The weak ground-state electronic transition at $689\mathrm{nm}$ must reach $5s5p{}^3{P}_1$: the electric-dipole selection rule on $J$ excludes the other two triplet state energy levels, while weak spin-orbit coupling admixes singlet state character and relaxes the approximate spin selection rule.

Assume the usual ordering $E({}^3{P}_0)<E({}^3{P}_1)<E({}^3{P}_2)$, as expected from a positive effective spin-orbit coupling constant. The three nearby strong spectral lines then join these energy levels to a common $5s6s{}^3{S}_1$ energy level. Larger lower-state energy means a smaller gap and longer wavelength, so $679,688,707\mathrm{nm}$ correspond to $J=0,1,2$ respectively. The remaining $1.12\mu \mathrm{m}$ spectral line joins $5s5p{}^1{P}_1$ and $5s6s{}^1{S}_0$. Assume the latter is the upper energy level. These assignments and orderings agree with independently tabulated strontium energy levels.^[ref]

Starting with the two ground-state gaps, the construction isThe resulting energy levels, rounded to avoid suggesting spectroscopic precision, are

Both $461\mathrm{nm}$ and $1.12\mu \mathrm{m}$ exhibit the normal Zeeman effect. Each connects two singlet state energy levels, with Landé g-factor $g_J=1$ for the $J=1$ energy level and no first-order shift for $J=0$. Their frequency shifts are $0,\pm {\mu }_{B}B/h$. The weak $689\mathrm{nm}$ spectral line also has three components, but $g({}^3{P}_1)=3/2$ gives shifts $0,\pm 3{\mu }_{B}B/(2h)$; a three-component pattern alone does not establish the normal Zeeman effect. The $679\mathrm{nm}$ spectral line likewise has three components with $g({}^3{S}_1)=2$.

Let the dimensionful angular momentum operators obey $\mathbf{F}=\mathbf{I}+\mathbf{J}$ and take the interaction Hamiltonian as $H^{\prime }=A\mathbf{I}\cdot \mathbf{J}$. Sincethe energy eigenvalues at fixed $I,J$ arewhere $F=|I-J|,\dots ,I+J$. Subtraction gives the Landé interval ruleMany atomic conventions use dimensionless angular momentum operators and call the energy-valued constant itself $A$; then the right-hand side is $AF$. Each energy level has $2F+1$ magnetic sublevels with the same energy in zero field.

For fine structure, the effective interaction is $a_{LS}\mathbf{L}\cdot \mathbf{S}$ with dimensionless angular momentum operators. It originates in the spin-orbit interaction, including the magnetic interaction of electron spin with orbital motion and its relativistic correction. Within an isolated LS coupling atomic term, the labels are the electron configuration, $L,S,J,M_J$, and the interval between adjacent $J$ energy levels is $a_{LS}J$. This is smaller than the electrostatic separation of atomic terms.

For magnetic-dipole hyperfine structure, the nuclear magnetic dipole moment interacts with the magnetic field generated by the electrons. Contributions include the Fermi contact interaction of penetrating $s$ electrons, orbital magnetic fields, and the dipolar field of electron spin. Projected onto one electronic $J$ energy level, their rotationally invariant magnetic-dipole part has the form $a_{IJ}\mathbf{I}\cdot \mathbf{J}$. The labels are $I,J,F,M_F$, together with the approximate electronic labels such as $L,S$. This splitting is generally smaller than fine structure, reflecting the much smaller nuclear magnetic dipole moment. An electric nuclear quadrupole moment can add a different interaction and violate the simple Landé interval rule.

At the larger electrostatic energy scale, an effective two-spin exchange interaction can also be written $a_s{\mathbf{s}}_1\cdot {\mathbf{s}}_2$ within fixed orbitals. Its origin is the Coulomb interaction combined with the Pauli exclusion principle, rather than a large direct magnetic interaction. It separates singlet state and triplet state atomic terms; before spin-orbit coupling, useful labels include the electron configuration, $L,M_L,S,M_S$. General electrostatic atomic term separations need not obey a single Landé interval rule as $L$ varies.

Write $a_{\nu }=a/h$, so the Landé interval rule for hyperfine structure intervals in ordinary frequency is ${\nu }_F-{\nu }_{F-1}=a_{\nu }F$. For the first four selected positions, the successive positive differences areTheir ratios are approximately $4:5:6$, sinceThe decreasing frequencies therefore correspond to increasing $F$, with $a_{\nu }\simeq -0.706\mathrm{GHz}$, if increasing printed position means increasing physical frequency. Continuing the intervals predicts additional positions near $10.60-7(0.706)\simeq 5.66$ and then $5.66-8(0.706)\simeq 0.01\mathrm{GHz}$. They identify peaks $i$ and $l$ unambiguously. The full assignment for ${}^{151}\mathrm{Eu}$ isWith $J=11/2$, these endpoints requireEquivalently, a complete set of six hyperfine levels gives $2I+1=6$ when $I<J$. This inference assumes that every upper hyperfine level contributes an observed position and the negligible lower hyperfine structure lets coincident transitions be grouped by their upper energy. A line count without these completeness assumptions would not alone fix nuclear spin.

The supplied positions obey the Landé interval rule only approximately. For example, ${\nu }_e-{\nu }_i=4.97\mathrm{GHz}$, whereas $7({\nu }_a-{\nu }_b)/4=4.935\mathrm{GHz}$. Even allowing each printed position a rounding error of $0.005\mathrm{GHz}$ does not fully reconcile these two relations. A least-squares dipole-only fit has a largest position residual of about $0.017\mathrm{GHz}$. Small additional hyperfine interactions, measurement uncertainty beyond the displayed precision, or an error in the tabulated data could account for this; it is inappropriate to claim an exact identity. The grouping and nuclear spin are nevertheless robust. Reversing the convention for the printed frequency axis reverses the inferred signs of both constants, without changing the following ratio.

Define $r=g_2/g_1$ and let $\sigma ={\sigma }_{21}(\omega )$ be the stimulated-emission cross section, averaged over upper magnetic sublevels. The Einstein coefficients relation $g_1{B}_{12}=g_2{B}_{21}$ gives the absorption cross section ${\sigma }_{12}=r\sigma$. Assume equally populated magnetic sublevels within each energy level, a steady pump, and a monochromatic beam. The photon flux is $I/(\hslash \omega )$; define $W=\sigma I/(\hslash \omega )$ and the effective population inversion $D=N_2-r{N}_1$. Then the rate equations areThe $WD$ term includes both stimulated emission and absorption. In steady state, adding the two rate equations gives $N_1=R_2{\tau }_1$. The first equation then yieldsThe directed beam gains energy at rate $\hslash \omega WD$ per volume, giving $dI/dz=\sigma DI$ and the optical gain coefficient $\alpha =\sigma D$. At line centre, write ${\sigma }_0={\sigma }_{21}({\omega }_{21})$. ThereforeHere ${\alpha }_0$ is the small-signal optical gain coefficient and $I_s$ is the saturation intensity, at which the local optical gain is half its small-signal value. Positive optical gain requires $r{A}_{21}{\tau }_1<1$: the lower energy level must empty sufficiently quickly.

The cancellation of ${\tau }_1$ from $I_s$ is specific to this model's steady fixed pump and the fact that every upper spontaneous emission feeds the lower energy level. Adding an extra lower-state accumulation factor to $I_s$ would contradict these rate equations. Alternative conventions using the absorption cross section must replace ${\sigma }_0$ by ${\sigma }_{12}/r$ throughout. Homogeneous spectral broadening permits the whole selected ensemble to share this local saturation response.

Use the condensate wavefunction convention $\int |\psi {|}^2{d}^{3}r=N$, so $n(\mathbf{r})=|\psi (\mathbf{r}){|}^2$ is the number density and $\psi /\sqrt{N}$ is a normalized one-particle orbital. The contact interaction strength is $U_0=4\pi {\hslash }^{2}a/m$ for dilute atoms with s-wave scattering length $a$. The chemical potential $\mu =\partial E/\partial N$ fixes the stationary time dependence $e^{-i\mu t/\hslash }$ and enforces normalization when minimizing the energy. The Gross–Pitaevskii equation follows from an energy functional that contains the pair-counting factor $1/2$:Varying $E-\mu \int |\psi {|}^2{d}^{3}r$ gives the stationary Gross–Pitaevskii equation; the interaction contribution to that equation is $U_{0}n\psi$, without the factor $1/2$.

For the broadening discussed here assume repulsive interactions, $a>0$. The Thomas–Fermi approximation for a condensate neglects the gradient kinetic energy compared with the trapping and interaction energies. At fixed spatial scale, increasing $N{U}_0$ increases the mean-field potential. Expansion then reduces the kinetic energy further: for a cloud radius $R$, its scale per atom is ${\hslash }^2/(m{R}^2)$, whereas the trap energy scale is $m{\omega }^2{R}^2$. Their ratio decreases as $(a_0/R{)}^4$. In the weak-interaction limit the kinetic energy instead balances confinement, giving the Gaussian wave packet ground state of the quantum harmonic oscillator.

The spherical Thomas–Fermi condensate hasIntegrating this density gives $N=8\pi \mu R^3/(15U_0)$. ThusAt fixed $N$, stronger repulsion produces a wider, lower wavefunction: the Gaussian wave packet crosses over towards the square root of an inverted parabola, with an almost flat central region and a steeper outer drop. The density is the inverted parabola; the wavefunction is its square root. The exact wavefunction remains smooth and has a tail, since kinetic energy cannot be neglected near the Thomas–Fermi condensate edge rounding region. The healing length must be small compared with the cloud size in the bulk. Attractive interactions would not follow this large-$N$ expansion argument.

The original plot compares the noninteracting endpoint with two normalized repulsive spherical Thomas–Fermi condensate profiles. The latter are asymptotic shapes, rather than exact intermediate solutions of the Gross–Pitaevskii equation. Each has the same integrated atom number; $\chi =Na/a_0$ labels the illustrative interaction strengths.

Choose $A$ real and positive. The Gaussian integral givesThe three zero-point energies add for each atom. The ideal-gas chemical potential is a one-atom energy; the total energy includes all $N$ atoms.

For weak interactions, the first-order correction in perturbation theory is the energy functional of the contact interaction evaluated on the unperturbed Gaussian wave packet. The integral needed isConsequently the Gaussian interaction energy of a trapped condensate isHere $\hslash \omega ={\hslash }^2/(m{a}_0^2)$. These are the mean-field expressions consistent with the supplied equation and $\int |\psi {|}^2=N$. Exact first-order perturbation theory for precisely $N$ atoms in one orbital replaces $N^2$ by $N(N-1)$ through the number of distinct pairs; this distinction is negligible in the large-$N$ regime of the crossover estimate.

To judge the ideal-gas approximation, compare the interaction energy with the total ideal energy, or compare the interaction energy per atom with $\hslash \omega$:Thus the crossover from Gaussian to Thomas–Fermi condensate is governed by $\chi =Na/a_0$: perturbation theory is reliable for $\chi \ll 1$, while the regime of the Thomas–Fermi approximation for a condensate requires $\chi \gg 1$. The characteristic crossover number is $N_c\sim a_0/a$; there is no sharp critical atom number here. If one specifically defines crossover by $E_I=E_0$ using the weak-interaction expression, it gives $N_c=(3\sqrt{2\pi }/2)a_0/a$. That prefactor is a convention-dependent estimate extrapolated to the edge of perturbation theory, not a phase transition. Comparing $E_I$ itself to a single $\hslash \omega$ would incorrectly omit a factor of $N$.

For the numerical estimate use $m\simeq 87\times 1.66054\times 10^{-27}\mathrm{kg}$ and $\hslash \simeq 1.05457\times 10^{-34}\mathrm{J}\mathrm{s}$. ThenThe particular equality $E_I=E_0$ would instead give about $1.15\times 10^3$ atoms. These quantify different choices within a smooth crossover; a well-developed spherical Thomas–Fermi condensate profile requires going substantially beyond the characteristic scale.

Laser cooling removes atomic kinetic energy through a velocity-dependent optical force, with scattered photons carrying away energy and entropy. For two counterpropagating beams detuned below resonance, motion towards one beam shifts that light closer to resonance, so the atom preferentially absorbs momentum opposite its velocity. At small velocity the net force is $F\simeq -{\gamma }_{v}v$ with ${\gamma }_v>0$. Three orthogonal beam pairs form optical molasses. Spontaneous emission has zero mean recoil in a symmetric distribution, but its momentum variance is nonzero: random absorption and emission cause heating. The cooling and heating balance sets a finite temperature.

For a weakly saturated two-level transition, write $\delta ={\omega }_L-{\omega }_0<0$ and let $\Gamma$ be the excited-state population decay rate. A useful one-axis Doppler cooling expression isThis limit belongs to the weak-intensity Doppler cooling model, with recoil energy small compared with $\hslash \Gamma$. Experimentally, released sodium molasses clouds were probed by time of flight; the unusually small expansion demonstrated temperatures below that limit. This distinguishes an observed failure of a two-level cooling model from a failure of conservation of energy.^[ref]

Multilevel atoms permit polarization-gradient cooling: spatially varying light polarization produces different light shifts of ground-state sublevels. In Sisyphus cooling, atoms climb a light-shift potential and are preferentially optically pumped near its top into a lower-energy sublevel, repeatedly losing kinetic energy. The optical pumping and emitted light supply dissipation. This can reach below the Doppler cooling limit; the recoil scalethen becomes important. The definition of a recoil temperature sometimes differs by a factor of two, so the energy definition is the unambiguous comparison. Recoil is not a universal lower bound for every cooling protocol: suitable dark-state or resolved sideband cooling schemes can suppress scattering near the desired motional state. Resolved sideband cooling preferentially excites a transition lowering the vibrational quantum number, followed by decay that on average preserves it in the Lamb-Dicke regime.^[ref]

Optical molasses supplies damping rather than a position-restoring force. A magneto-optical trap adds a magnetic-field gradient and appropriately chosen circular polarizations; the position-dependent Zeeman splitting gives a restoring radiation pressure imbalance as well as damping. In real alkali atoms, repumping light returns population that would otherwise accumulate in an uncoupled hyperfine level. Magnetic and optical dipole traps can subsequently hold the cold sample without continuous resonant scattering.

Applications include long-interrogation atomic clocks, precision spectroscopy, atom interferometry, and preparation for quantum simulation and Bose-Einstein condensation. Cooling reduces motion-induced frequency shifts and allows longer coherent interrogation; deeper quantum degeneracy commonly needs subsequent evaporative cooling of atoms. Finite capture velocity, available transitions, recoil heating, stray magnetic fields, and loss or heating from unwanted scattering constrain a real experiment.

Evaporative cooling of atoms selectively removes atoms with unusually high total energy, then elastic collisions rethermalize those remaining. Loss of particles alone is insufficient: the lost atoms must carry more than the mean energy per particle if the temperature is to fall. The process is open and dissipative, so it can increase the phase-space density of an atomic gas in the remaining sample while reducing its atom number. As the temperature falls, the cutoff must be lowered to keep replenishing a useful escaping tail; a fixed high cutoff eventually gives negligible evaporation.

An experimental sequence first loads laser-cooled atoms into a conservative trap and suppresses resonant-light heating. In a trap using magnetic trapping of neutral atoms, a radiofrequency evaporation knife drives a Zeeman splitting transition from a trapped weak-field-seeking state to an untrapped state on a surface where the splitting matches the radiofrequency photon energy. Only atoms energetic enough to reach that surface escape. Lowering the radiofrequency cutoff progressively removes less energetic atoms. Alternatively, lowering an optical dipole trap's depth permits energetic atoms to spill out. This changes confinement as well as the escape cutoff, so it can reduce the collision rate. Radiofrequency-selective evaporation and the need for favourable elastic-collision and background-loss rates are documented in the original condensate experiments.^[ref]

The elastic collision rate ${\Gamma }_{\mathrm{el}}\sim n{\sigma }_{\mathrm{el}}\bar{v}$ must be large enough for elastic rethermalization of a trapped atomic gas throughout evaporation, and much faster than unwanted loss and heating. A clean vacuum reduces background collisions; suitable trapping states suppress inelastic spin-changing loss, and density must not make three-body recombination dominant. The cutoff ramp should allow equilibration without taking so long that background losses dominate. A large evaporation truncation parameter $\eta ={ϵ}_{\mathrm{cut}}/(k_{B}T)$ selects energetically valuable atoms, but an excessively large value suppresses their escape rate. A very small value discards too many useful low-energy atoms. If elastic collisions become more frequent as evaporation proceeds, the resulting runaway evaporative cooling can improve efficiency; it is not automatic in a weakening optical trap.

The energy selectivity for evaporation to increase phase-space density can be made precise in a classical three-dimensional harmonic trap with fixed frequencies. Equipartition gives $E=3N{k}_{B}T$. If each escaping atom removes mean energy ${ϵ}_{\mathrm{ev}}$, energy balance is $dE={ϵ}_{\mathrm{ev}}dN$ with $dN<0$, soThe peak harmonic-trap phase-space density is $\mathcal{D}=N(\hslash \bar{\omega }/(k_{B}T){)}^3$. ThereforeCooling requires ${ϵ}_{\mathrm{ev}}>3k_{B}T$; increasing $\mathcal{D}$ despite atom loss requires the stronger condition ${ϵ}_{\mathrm{ev}}>4k_{B}T$. These inequalities assume rapid thermalization and no work from changing trap frequencies. They distinguish lower temperature from genuinely improved quantum degeneracy. Eventually quantum statistics and inelastic processes require a refinement of this classical estimate. Efficient evaporation combines energy-selective loss, rapid elastic redistribution, and low unwanted loss; it can bring a precooled gas to Bose-Einstein condensation.

For a free-space electric dipole transition, the spontaneous emission rate in SI units isHere $|\mathbf{d}{|}^2$ includes the relevant dipole components; a specified single polarization requires its corresponding matrix element.

A complete scattering cycle returns the internal state to the ground state and changes the atomic momentum to $\mathbf{P}+\mathbf{p}-{\mathbf{p}}^{\prime }$. ThusFor photon-scattering recoil heating with isotropic emission of fixed momentum magnitude, $⟨{\mathbf{p}}^{\prime }⟩=0$, givingThe first term is the velocity-dependent work of absorption and can cool an atom if the photon opposes its motion. The last two are positive absorption and emission recoil energies. If the actual absorbed directions are also isotropic and independent of $\mathbf{P}$, the first term averages to zero. For approximately equal incoming and outgoing photon momenta $\hslash k$,Isotropic incident illumination alone would not guarantee isotropic absorption if Doppler-dependent excitation selected directions.

Source qualification for the trapped calculation. Set $m=M$ for the same atom. The printed Gaussian $\exp [-x^2/(2l_0{)}^2]$ has the wrong width for the specified harmonic oscillator: the printed ground and first-excited functions have squared norms $\sqrt{2}$ and $2\sqrt{2}$, respectively. The normalized oscillator states fixed by the given $l_0=\sqrt{\hslash /(m{\omega }_{\mathrm{HO}})}$ instead useAlso, one common complex phase multiplying both off-diagonal terms does not make the Hamiltonian a Hermitian operator. Use an absorption matrix element $(\hslash \Omega /2)e^{-ikx}$ and its Hermitian conjugate for emission; choosing the opposite propagation convention reverses the imaginary sideband phase but not a transition probability. The physical results below use these explicitly stated corrections. Parts (i) and (ii) also give the literal integrals of the printed functions.

Conditional motional heating and recoil averaging. With $\eta =k{l}_0/\sqrt{2}$, an oscillator eigenstate has $⟨K{⟩}_n=⟨V{⟩}_n=(n+1/2)\hslash {\omega }_{\mathrm{HO}}/2$. Therefore the absorption outcomes haveThe same statements hold for the specified emission outcomes $|e,0⟩\to |g,0⟩$ and $|e,0⟩\to |g,1⟩$. Internal excitation energy goes into, or comes from, the photon and is distinct from this external energy change. Emission into the first motional state makes the emitted photon correspondingly less energetic.

A sudden recoil from $|0⟩$ gives probabilities $P_n=e^{-{\eta }^2}{\eta }^{2n}/n!$, so $P_1\simeq {\eta }^2\ll 1$ while $P_0\simeq 1$. The rare first sideband transfers a whole trap quantum, much larger than $E_R$; the mean is neverthelessThe emission average is the same for a one-dimensional recoil of the same magnitude. Projecting isotropic three-dimensional emission onto one axis instead gives a factor $⟨{\cos }^2\theta ⟩=1/3$ for that axis.

A sudden unobserved kick initially raises kinetic energy by $E_R$ and leaves the position distribution, hence potential energy, unchanged. It creates a coherent superposition of motional states. The equal kinetic/potential increments above concern energy-eigenstate-conditioned outcomes, or time-averaged/dephased motion; these are different measurement conditions. There is no contradiction.

With $\Gamma >{\omega }_{\mathrm{HO}}$, motional sidebands are not cleanly resolved, and a short-lived excited atom has little time to move before decay. In the stronger sudden limit $\Gamma \gg {\omega }_{\mathrm{HO}}$ the recoil-kick picture is direct. For only the stated inequality, finite linewidth weights can change numerical factors; it does not fix an exact continuous-wave scattering rate. The average heating remains of recoil order per absorption or emission, and of order $2E_R$ when both absorption and emission kicks are counted at full magnitude. For the axial component of isotropic three-dimensional emission, the emission contribution instead carries the factor $1/3$ noted above. Strong confinement suppresses the chance of changing $n$, but does not make its averaged recoil energy disappear. This agrees with free-space recoil heating while replacing small continuous energy changes by rare quantized motional changes.

For the Bogoliubov spectrum, define ${ϵ}_k={\hslash }^2{k}^2/(2m)$ and $E_k=\sqrt{{ϵ}_k({ϵ}_k+2{\mu }_0)}$. For repulsive interactions and ${\mu }_0>0$,Thus the low-momentum excitations are gapless phonons and the high-momentum branch approaches free-particle curvature. The crossover occurs when kinetic and interaction energy are comparable. The Landau criterion follows because an obstacle moving at velocity $v$ can supply energy $\mathbf{v}\cdot \mathbf{p}$ to an excitation. HereThis is the homogeneous phonon kinematic threshold, not a guarantee that every perturbation dissipates immediately above it.

In two dimensions, $n=|\psi {|}^2$ has units of inverse area and ${\hslash }^2/m$ has units of energy times area. Therefore $\tilde{g}$ is dimensionless. Neglecting the kinetic term in the Gross–Pitaevskii equation, with positive coupling, giveswith $n=0$ outside the radius in the Thomas–Fermi approximation for a condensate. Multiplication givesThe sharp edge is an approximation; the neglected gradient term rounds it over a healing length scale.

Using the local sound-speed criterion for circular condensate stirring, the moving beam has speed $v=\Omega D$. The local interaction chemical potential is ${\mu }_0-V(D)=m{\omega }_r^2(R^2-D^2)/2$, soThe onset equality is $\Omega D=c(D)$ and the permitted excitation condition isThis uses the beam speed as the local relative-flow scale. A finite hole can distort the surrounding flow or nucleate vortices, so the result is the requested local phonon estimate, not a universal necessary condition for all experimental stirring.

For the loss-driven onset of condensate stirring, first integrate the density:At fixed coupling and trap frequency, ${\mu }_0\propto N^{1/2}$ and $R\propto N^{1/4}$. Hence $R(t)=R(0)e^{-t/(4\tau )}$ while $D$ stays fixed. For $\Omega ={\omega }_r/\sqrt{2}$, threshold requires $D/R=1/\sqrt{2}$. Starting at $D/R(0)=1/2$,At onset, $N=N_0/4$, $R=R(0)/\sqrt{2}$ and the beam is still inside the cloud. Excitations become kinematically allowed just beyond this equality, provided the stipulated adiabatic regime of the Thomas–Fermi approximation for a condensate remains valid.

The Lamb shift is a radiative displacement of atomic energy levels beyond the values from the Dirac equation in a Coulomb potential. In hydrogen, the $2S_{1/2}$ and $2P_{1/2}$ levels are degenerate for a point Coulomb potential in the Dirac equation, yet the $2S_{1/2}$ level lies higher by roughly a gigahertz. This is distinct from the fine structure splitting between different $J$ levels and from nuclear-spin hyperfine structure.

Electron self-energy in the quantized electromagnetic field is the dominant contribution to the familiar hydrogen splitting; vacuum polarization contributes too. Renormalization separates the physical electron mass from its divergent free-electron self-energy. A useful leading-log scale for a hydrogenic $nS$ level isnot an exact numerical formula. The enhanced sensitivity of S states follows from their nonzero wavefunction at the nucleus; leading short-distance corrections act like contact perturbations. A precision prediction also includes finite nuclear size and recoil, and must consistently separate those effects from the radiative contribution.^[ref]

The original Lamb-Retherford experiment used a beam of hydrogen atoms, electron-impact excitation into metastable $2S$, and microwave-induced transfer into short-lived $2P$ states. Those atoms decayed before reaching the detector. Metastable atoms surviving to a metal target ejected electrons, so a drop in the measured current signalled a resonance. Varying the microwave excitation and magnetic field established the displacement from the Dirac equation expectation. Microwave quenching converted a small energy difference into a measurable loss of a long-lived population, avoiding optical Doppler widths that obscure the splitting.^[ref]

A Paul trap uses a time-dependent electric quadrupole field to provide dynamic confinement, avoiding the impossibility of trapping a charge in a stable three-dimensional static electrostatic minimum. Its dynamics follows a Mathieu equation: stable drive parameters produce slow Paul trap secular motion with superposed fast Paul trap micromotion. If ${\mathbf{E}}_{\mathrm{rf}}$ denotes the peak oscillating field, the effective energy for small Paul trap micromotion isThe radiofrequency pseudopotential of a Paul trap is an averaged mechanical energy, not a static electrostatic potential. Doppler cooling and, where needed, resolved sideband cooling reduce Paul trap secular motion; the Lamb-Dicke regime suppresses higher motional sidebands.^[ref]

Two long-lived internal states encode a qubit, for example hyperfine states, while an optical metastable transition provides another option. Optical pumping initializes it. Resonant or Raman transition pulses make coherent single-qubit rotations, and state-dependent fluorescence distinguishes a bright state from a dark shelved state. Shared quantized normal modes let one ion's internal state influence another, providing entangling gates in the proposal of J. I. Cirac and P. Zoller; the motional bus must be returned to its intended state rather than left carrying unwanted information.^[ref]

A concrete early experiment demonstrated a controlled-NOT gate between the internal hyperfine and motional degrees of freedom of a single trapped ion. This was an internal–motion two-qubit gate, not yet a gate between two separate ions.^[ref]

Challenges include electric-field-noise heating of the motional bus, excess Paul trap micromotion from stray fields, magnetic and laser phase noise, off-resonant spontaneous scattering, imperfect preparation and readout, and unwanted excitation of spectator modes. Scaling to many ions adds mode crowding, individual-addressing crosstalk and transport requirements. Cooling and quantum error correction must remove entropy without erasing stored qubits. Long storage coherence is necessary but not sufficient: gates and measurements must also be accurate and reproducible. The original gate experiment explicitly identified laser, field, trap-drive and motional decoherence limitations.^[ref]

For hydrogen, the leading characteristic scales areOrder-one spin and nuclear magnetic factors, and state-dependent factors, are suppressed here. Gross atomic structure is Coulomb binding. Fine structure includes relativistic kinetic and spin-orbit effects, with Thomas precession affecting the spin-orbit coefficient. Hyperfine structure couples electronic motion and spin to the nuclear magnetic moment, whose inverse nuclear mass makes it much weaker. The contact spin interaction can remain in an s state where orbital spin-orbit coupling vanishes.

For positronium, the reduced mass is $m_e/2$. Its gross levels and separation scale areThus the mean ground-state electron–Positron separation is twice the infinite-proton-mass hydrogen value: both means have the same factor $3/2$ multiplying their respective Bohr radii. Retaining proton recoil makes the ratio $2m_p/(m_p+m_e)$, very slightly below two. Separation is not the distance of either positronium constituent from its centre of mass.

Interpret the spin operators in units of $\hslash$, as required by the printed Zeeman term, so $A$ has units of energy. Two spin halves give four eigenstates. At zero field, total spin $F=0,1$ and projection $m_F$ label the singlet and triplet. At high field relative to the hyperfine scale, the individual projections $m_{s,e},m_{s,p}=\pm 1/2$ become approximately good quantum numbers. Total projection remains exact at every field, but $F^2$ does not for the two mixed states.

Choose the quantization axis along $\mathbf{B}$. The aligned-spin states $|\uparrow \uparrow ⟩$ and $|\downarrow \downarrow ⟩$ both have $E=A/4$ because their Zeeman shifts cancel. In the other product-state sector,For physical positronium $A>0$. At $B=0$ there is a triplet at $A/4$ and a singlet at $-3A/4$. At large positive field the mixed states approach $|\uparrow \downarrow ⟩$ and $|\downarrow \uparrow ⟩$, with energies $-A/4+2{\mu }_{B}B$ and $-A/4-2{\mu }_{B}B$ respectively, while the aligned pair remains degenerate at $A/4$.

The dimensionless mixing parameter is $x=4{\mu }_{B}B/A$. A crossover estimate is $x\sim 1$, rather than a sharp transition. The scale $A\sim {\alpha }^4{m}_e{c}^2$ already gives a field of a few tesla. Using the physical hyperfine interval $A/h\simeq 203.39\mathrm{GHz}$ givesThe numerical hyperfine interval is external physical input, not a number supplied by the symbolic Hamiltonian.^[ref] The high-field limit here still presumes that orbital distortion is negligible compared with Coulomb binding.

Set the ground-state energy to zero. Bose-Einstein occupation of an excited level is $[e^{(ϵ-\mu )/(k_{B}T)}-1{]}^{-1}$, with $\mu \le 0$. As particles are added, $\mu$ approaches zero and the excited-state population saturates. Separating the ground state from the density-of-states integral therefore givesThe excited-state integral must actually converge: $g\to 0$ alone is not a sufficient condition in complete generality. A power law $g\propto {ϵ}^s$ with $s>0$ suffices at low energy, as in the stated three-dimensional box.

Let $\mathcal{C}=V{m}^{3/2}/(\sqrt{2}{\pi }^2{\hslash }^3)$ denote the density-of-states prefactor, not a heat capacity. Substitution $x=ϵ/(k_{B}T)$ givesThe condensate has zero energy on this reference, so a further factor of $ϵ$ in the integral givesRelative to the classical value, the ratio is $\zeta (5/2)/[\eta \zeta (3/2)]\simeq 0.5134/\eta <0.5134$. Macroscopic occupation of the zero-energy state and the Bose distribution both lower the average energy.

At fixed volume the condensed gas has $E\propto T^{5/2}$. Unbiased random loss keeps average $E/N$ constant, henceRemoving a condensate atom removes essentially no energy and by itself leaves the excited cloud's temperature unchanged. Removing an excited atom removes its energy; during re-equilibration, some zero-energy condensate atoms join the excited population and share the remaining energy, lowering temperature. Averaging over unbiased loss gives the stated constant energy per particle, even though each individual type of loss need not preserve it separately.

Using $m\simeq 87u$ and the supplied volume,There is only one occupied internal spin state, so no factor $2F+1$ belongs in the density of states. During loss, $\eta \propto N/T^{3/2}\propto N^{2/5}$, hence $\eta /{\eta }_i=T/T_i$. Condensation disappears at ${\eta }_f=1$:As a consistency check, $N_f=N_i{\eta }_i^{-5/2}\simeq 8.15\times 10^4$ agrees with $N_c(T_f)$.

Both this process and evaporative cooling of atoms require particle loss and re-equilibration. Ordinary evaporation preferentially removes energetic atoms, reducing $E/N$, and can cool an uncondensed classical gas. Here loss is energy independent, $E/N$ stays constant, and cooling relies on the condensed gas's unusual equation of state; random loss of a classical ideal gas would leave its temperature unchanged. Exact noninteraction would prevent collisional rethermalization: the imposed equilibrium assumption can be understood as weak thermalizing interactions whose equation-of-state corrections are neglected.

Use the ordered basis $(|1⟩,|2⟩)$ and ${\sigma }_z=\mathrm{diag}(1,-1)$, withA pure quantum superposition $|+x⟩=(|1⟩+|2⟩)/\sqrt{2}$ has $\mathbf{R}=(1,0,0)$ on the Bloch sphere. In contrast, a mixed state $\rho =p|1⟩⟨1|+(1-p)|2⟩⟨2|$ has $\mathbf{R}=(0,0,2p-1)$, inside the sphere for $0<p<1$. At equal populations its vector is zero. Equal populations alone therefore do not distinguish a coherent superposition from a mixture; their off-diagonal density matrix elements differ.

Choose the drive phase so the rotating-frame Hamiltonian is $H_{\mathrm{rot}}=\hslash (\Omega {\sigma }_x+\delta {\sigma }_z)/2$, where $\delta =\omega -{\omega }_0$ as specified. The Rabi vector is the angular-velocity vectorHere $\Omega$ is the on-resonance population oscillation angular frequency, not the unhalved coefficient of a cosine coupling. From $\mathbf{R}(0)=\hat{z}$, rotation about $\mathbf{W}$ at $W=\sqrt{{\Omega }^2+{\delta }^2}$ givesFor $\Omega \gg |\delta |$, the rotation axis is nearly transverse, the trajectory nearly a great circle, and transfer is almost complete. For $\Omega =\delta >0$, the axis bisects $+x$ and $+z$, the orbit lies in $R_x+R_z=1$, and $p_2=\frac{1}{2}{\sin }^2(\sqrt{2}\Omega t/2)$ never exceeds $1/2$. Each member remains a pure state in this coherent evolution.

For the Ramsey free interval, let atom $j$ have detuning ${\delta }_j=\delta +{\eta }_j$. Instantaneous ideal pulses give the Ramsey ensemble coherence factor $\chi (T)=⟨e^{i{\eta }_{j}T}⟩$. Averaging the free-evolution vector obtained in step (ii) givesDifferent atoms fan out around the equator, shortening the ensemble Bloch vector even though every individual's vector remains unit length. A broad smooth detuning distribution usually makes the average spiral towards the centre and damps the Ramsey fringes towards $1/2$. The dephasing time is of order $1/\Delta {\omega }_0\sim 5/\delta$. inhomogeneous ensemble dephasing is not necessarily irreversible decoherence of individual atoms; an appropriate spin echo can rephase static frequency offsets.

The given spread does not specify a probability distribution or whether it denotes an RMS width. An exact envelope is therefore not determined. The figure illustrates, rather than infers, a Gaussian distribution with RMS ${\sigma }_{\delta }=0.2\delta$, for which $\chi =e^{-{\sigma }_{\delta }^2{T}^2/2}$. A uniform distribution would instead give a sinc envelope and side lobes. Both have the same qualitative initial loss of contrast on the quoted inverse-spread scale.

In the rubidium-87 ground-state hyperfine structure, $J=1/2$ and $I=3/2$, so $F=1,2$. Taking $H_{\mathrm{hf}}=A\mathbf{I}\cdot \mathbf{J}/{\hslash }^2$ givesThe zero-field splitting is ${\Delta }_{\mathrm{hf}}=2A\simeq h\times 6.8\mathrm{GHz}$. Ignoring the small nuclear magnetic term, the weak-field Landé g-factors are $g_{F=1}\simeq -1/2$ and $g_{F=2}\simeq +1/2$. The diagram labels all three $F=1$ and five $F=2$ substates, with linear shifts $g_F{\mu }_B{m}_{F}B$.