Skip to content

A dressed state is an eigenstate of a Hamiltonian that includes both a matter system and its coherent coupling to a field. The word “dressed” emphasizes that the interaction changes both the energies and the state vectors. The eigenstates are no longer bare atomic states, and, when the field is quantized, they are generally not product states of an atom and a definite photon number.

For the standard near-resonant two-state Hamiltonian,

HRWA=ℏ2[Δσz+Ω(cos⁡ϕ σx+sin⁡ϕ σy)],H_{\mathrm{RWA}} = \frac{\hbar}{2} \left[ \Delta\sigma_z + \Omega \left( \cos\phi\,\sigma_x + \sin\phi\,\sigma_y \right) \right],

the dressed-state separation is

ℏΩR,ΩR=Δ2+Ω2.\hbar\Omega_R, \qquad \Omega_R = \sqrt{ \Delta^2+\Omega^2 }.

The same square-root spectrum appears in one excitation manifold of the Jaynes–Cummings model, with the replacement

Ω⟷2gN.\Omega \longleftrightarrow 2g\sqrt N.

That resemblance is powerful, but it must be interpreted carefully. A semiclassical rotating-frame eigenstate, a joint atom–photon eigenstate, and a Floquet mode are related descriptions, not literally the same object.

Required background. Two-Level Atom fixes the rotating-frame Hamiltonian; Rabi Oscillations explains coherent population transfer; Tensor Products and the Oscillator Ladder Operators supply joint states and photon-number matrix elements.

The phrase “atom plus photons” is often used before the field model has been declared. Three common constructions should be separated.

Semiclassical RWA. The state space is the atomic Hilbert space, and the stationary object is a rotating-frame eigenvector. There is no extra field index.

Classical periodic drive. Floquet analysis extends the atomic space by Fourier harmonics. The stationary object is a Floquet mode, and the extra integer labels a harmonic or quasienergy zone.

Quantized single mode. The state space is the atomic Hilbert space tensored with oscillator Fock space. The stationary object is a joint atom–mode eigenvector, and the extra integer is a physical mode occupation.

For a prescribed classical drive, the source is not a quantum subsystem of the model. The rotating-frame dressed vectors live only in the matter Hilbert space. A Floquet construction introduces a Fourier index, sometimes drawn as a photon ladder, but that index is bookkeeping in an extended space.

For a quantized cavity mode, by contrast, ∣g,n⟩|g,n\rangle and ∣e,n⟩|e,n\rangle are physical tensor-product states. Their superpositions can entangle the atom with the mode, and photon-number fluctuations can produce dynamics with no deterministic classical-drive counterpart.

These descriptions agree in overlapping limits. They answer different questions outside those limits.

Use bare matter states satisfying

H0∣g⟩=Eg∣g⟩,H0∣e⟩=Ee∣e⟩,H_0|g\rangle = E_g|g\rangle, \qquad H_0|e\rangle = E_e|e\rangle,

with

ℏω0=Ee−Eg>0.\hbar\omega_0 = E_e-E_g > 0.

The drive has angular frequency ωL\omega_L, phase ϕ\phi, and atom-minus-laser detuning

Δ=ω0−ωL.\Delta = \omega_0-\omega_L.

Thus red detuning means Δ>0\Delta>0. In the ordered basis (∣e⟩,∣g⟩)(|e\rangle,|g\rangle), after removing a common quasienergy, use

HRWA=ℏ2(ΔΩe−iϕΩeiϕ−Δ),H_{\mathrm{RWA}} = \frac{\hbar}{2} \begin{pmatrix} \Delta & \Omega e^{-i\phi} \\ \Omega e^{i\phi} & -\Delta \end{pmatrix},

where Ω≥0\Omega\ge0. Frequency variables are angular frequencies; decay rates are inverse times. The two-level and rotating-wave approximations require the relevant coupling, detuning, and envelope bandwidth to be small compared with the optical sum frequency and the separation from excluded resonances.

Changing the sign convention for Δ\Delta, reversing the basis order, or absorbing ϕ\phi into a bare state changes the displayed formulas. It does not change consistently computed gaps or probabilities.

Define

ΩR=Δ2+Ω2\Omega_R = \sqrt{ \Delta^2+\Omega^2 }

and choose a mixing angle in the range 0≤θ≤π0\le\theta\le\pi through

cos⁡θ=ΔΩR,sin⁡θ=ΩΩR.\cos\theta = \frac{\Delta}{\Omega_R}, \qquad \sin\theta = \frac{\Omega}{\Omega_R}.

For a nonzero gap, this definition remains unambiguous on both sides of resonance. The normalized eigenvectors may be chosen as

∣+⟩=cos⁡θ2∣e⟩+eiϕsin⁡θ2∣g⟩,∣−⟩=−sin⁡θ2∣e⟩+eiϕcos⁡θ2∣g⟩.\begin{aligned} |+\rangle &= \cos\frac{\theta}{2}|e\rangle + e^{i\phi} \sin\frac{\theta}{2}|g\rangle, \\ |-\rangle &= - \sin\frac{\theta}{2}|e\rangle + e^{i\phi} \cos\frac{\theta}{2}|g\rangle. \end{aligned}

They satisfy

HRWA∣±⟩=ε±∣±⟩,ε±=±ℏΩR2.H_{\mathrm{RWA}}|\pm\rangle = \varepsilon_\pm|\pm\rangle, \qquad \varepsilon_\pm = \pm \frac{\hbar\Omega_R}{2}.

An overall phase on either eigenvector is arbitrary. The relative phase inside a dressed state is not arbitrary once the phases of the bare basis and drive have been fixed.

The excited-state weights are

∣⟨e∣+⟩∣2=12(1+ΔΩR),∣⟨e∣−⟩∣2=12(1−ΔΩR).\begin{aligned} |\langle e|+\rangle|^2 &= \frac12 \left( 1+\frac{\Delta}{\Omega_R} \right), \\ |\langle e|-\rangle|^2 &= \frac12 \left( 1-\frac{\Delta}{\Omega_R} \right). \end{aligned}

The corresponding ground-state weights are obtained by subtraction from one. Equivalently,

⟨±∣σ∣±⟩=±1ΩR(Ωcos⁡ϕ,Ωsin⁡ϕ,Δ).\langle\pm|\boldsymbol\sigma|\pm\rangle = \pm \frac1{\Omega_R} \left( \Omega\cos\phi, \Omega\sin\phi, \Delta \right).

Each dressed state points along or opposite the effective field on the Bloch sphere.

On resonance,

Δ=0,θ=π2,\Delta=0, \qquad \theta=\frac{\pi}{2},

so the states are equal-weight superpositions:

∣±⟩=12(∣e⟩±eiϕ∣g⟩),|\pm\rangle = \frac1{\sqrt2} \left( |e\rangle \pm e^{i\phi}|g\rangle \right),

up to an overall phase and the sign convention used to label the branches. The gap is ℏΩ\hbar\Omega.

For large positive detuning,

Δ≫Ω,\Delta\gg\Omega,

the upper state becomes excited-like and the lower state becomes ground-like. For large negative detuning, their bare character is exchanged. A branch label such as “ground-like” is therefore local in parameter space; the fixed labels ++ and −- denote upper and lower quasienergies within this selected RWA block; they are not a global ordering of all equivalent Floquet-zone representatives.

The bare states can be expanded in the dressed basis:

∣e⟩=cos⁡θ2∣+⟩−sin⁡θ2∣−⟩,eiϕ∣g⟩=sin⁡θ2∣+⟩+cos⁡θ2∣−⟩.\begin{aligned} |e\rangle &= \cos\frac{\theta}{2}|+\rangle - \sin\frac{\theta}{2}|-\rangle, \\ e^{i\phi}|g\rangle &= \sin\frac{\theta}{2}|+\rangle + \cos\frac{\theta}{2}|-\rangle. \end{aligned}

This inverse relation explains resonant Rabi flopping. A suddenly prepared bare state is generally a coherent superposition of two stationary dressed states. Their relative phase evolves at ΩR\Omega_R, and transforming back to the bare basis produces population oscillation.

The dressed-state picture and the Bloch-sphere picture are therefore not competing explanations. One emphasizes stationary eigenvectors and phase accumulation; the other emphasizes precession about the same effective Hamiltonian vector.

With Ω=0\Omega=0, the rotating-frame bare quasienergies are

εe(0)=+ℏΔ2,εg(0)=−ℏΔ2.\varepsilon_e^{(0)} = +\frac{\hbar\Delta}{2}, \qquad \varepsilon_g^{(0)} = -\frac{\hbar\Delta}{2}.

They cross at Δ=0\Delta=0. A nonzero coupling replaces this crossing by

ε±=±ℏ2Δ2+Ω2,\varepsilon_\pm = \pm \frac{\hbar}{2} \sqrt{ \Delta^2+\Omega^2 },

whose minimum separation is

ε+−ε−∣Δ=0=ℏΩ.\left. \varepsilon_+-\varepsilon_- \right|_{\Delta=0} = \hbar\Omega.

An avoided crossing records two facts at once:

  • the eigenvalues repel because the states are coupled;
  • the eigenvectors exchange their bare-state character across resonance.

The second fact is easy to miss if only the energy curves are plotted. Hellmann–Feynman differentiation exposes it:

∂ε±∂Δ=ℏ2⟨±∣σz∣±⟩=±ℏΔ2ΩR.\frac{\partial\varepsilon_\pm}{\partial\Delta} = \frac{\hbar}{2} \langle\pm|\sigma_z|\pm\rangle = \pm \frac{\hbar\Delta}{2\Omega_R}.

Far from resonance, a branch slope approaches the slope of one bare state. At resonance both slopes vanish because each branch contains equal bare weight.

The avoided-crossing plot and exact far-detuned expansion are developed for light shifts on AC Stark Shift. The important lesson here is about basis identity: following an eigenvalue continuously is not the same operation as holding a bare label fixed.

Suppose the system begins in ∣g⟩|g\rangle and a constant drive is switched on much faster than 1/ΩR1/\Omega_R. The state does not instantaneously rotate into an eigenvector. Immediately after the switch it is still ∣g⟩|g\rangle, which is a superposition of ∣+⟩|+\rangle and ∣−⟩|-\rangle. Subsequent interference produces Rabi oscillations.

This is a quench into the dressed basis, not preparation of one dressed state.

If Ω(t)\Omega(t) and Δ(t)\Delta(t) change slowly while the dressed gap stays open, an initial instantaneous eigenstate follows its branch up to dynamical and geometric phases. For the phase-fixed two-state model,

θ˙=ΔΩ˙−ΩΔ˙ΩR2,\dot\theta = \frac{ \Delta\dot\Omega-\Omega\dot\Delta }{ \Omega_R^2 },

and a useful local condition is

∣θ˙∣≪ΩR.|\dot\theta| \ll \Omega_R.

The off-diagonal coupling in the instantaneous dressed basis is proportional to θ˙/2\dot\theta/2. For a varying phase the magnitude is instead 12θ˙2+sin⁡2θ ϕ˙2\tfrac12\sqrt{\dot\theta^2+\sin^2\theta\,\dot\phi^2}. Compare it with ΩR\Omega_R. These are local diagnostics: the full ramp history and duration also matter, because small transitions can accumulate resonantly.

At exact resonance, turning on from Ω=0\Omega=0 begins at a degeneracy. There is then no unique resonant dressed branch for the ordinary nondegenerate adiabatic theorem to select. A finite initial detuning, an explicit symmetry-breaking term, or a deliberately prepared superposition is needed to define the outcome.

For a linear sweep

Δ(t)=vt,Ω(t)=Ω,\Delta(t) = vt, \qquad \Omega(t) = \Omega,

the ideal infinite-time Landau–Zener probability to remain on the same diabatic branch is

Pdiab=exp⁡(−πΩ22∣v∣).P_{\mathrm{diab}} = \exp \left( - \frac{\pi\Omega^2}{2|v|} \right).

For finite-time corrections and the derivation, see Landau–Zener Transition.

An abrupt switch-off leaves the state unchanged during the ideal sudden step. A subsequent bare-state measurement samples its bare-basis probabilities; the switch itself is not a projection measurement. An adiabatic switch-off maps each dressed branch onto the bare state it approaches, provided the chosen path ends at nonzero detuning so the gap does not close. These protocols answer different experimental questions:

  • sudden release measures the instantaneous bare-state composition;
  • adiabatic release measures branch population;
  • a phase-controlled ramp can convert dressed coherence into a bare-state population signal.

Statements such as “the atom is half excited” are incomplete unless the switch-off and measurement protocol are specified.

The Jaynes–Cummings manifold gives a concrete example of quantized field dressing. For the derivation from the quantum Rabi interaction, exact propagator, number-statistics dependence, and collapse and revival, see the Jaynes–Cummings Model.

The near-resonant two-level, single-mode Hamiltonian is

HJC=ℏωca†a+ℏω02σz+ℏg(aσ++a†σ−),H_{\rm JC}=\hbar\omega_c a^\dagger a+ \frac{\hbar\omega_0}{2}\sigma_z+ \hbar g(a\sigma_++a^\dagger\sigma_-),

where a phase convention makes g≥0g\ge0. The conserved excitation number N=a†a+σ+σ−\mathcal N=a^\dagger a+\sigma_+\sigma_- selects, for N≥1N\ge1, the pair ∣e,N−1⟩|e,N-1\rangle, ∣g,N⟩|g,N\rangle. Notice that NN is total excitation, not the photon number of both components. With Δc=ω0−ωc\Delta_c=\omega_0-\omega_c, this pair has matrix

HN=ℏωc(N−12)I+ℏ2(Δc2gN2gN−Δc).H_N=\hbar\omega_c(N-\tfrac12)\mathbb I+ \frac{\hbar}{2} \begin{pmatrix}\Delta_c&2g\sqrt N\\2g\sqrt N&-\Delta_c\end{pmatrix}.

Thus the mixing-angle construction above applies with Ω↦2gN\Omega\mapsto2g\sqrt N and Δ↦Δc\Delta\mapsto\Delta_c:

ΩN=Δc2+4g2N,cos⁡θN=Δc/ΩN,∣N,+⟩=cos⁡(θN/2)∣e,N−1⟩+sin⁡(θN/2)∣g,N⟩,∣N,−⟩=−sin⁡(θN/2)∣e,N−1⟩+cos⁡(θN/2)∣g,N⟩,EN,±=ℏωc(N−12)±ℏΩN/2.\begin{aligned} \Omega_N&=\sqrt{\Delta_c^2+4g^2N},\qquad \cos\theta_N=\Delta_c/\Omega_N,\\ |N,+\rangle&=\cos(\theta_N/2)|e,N-1\rangle +\sin(\theta_N/2)|g,N\rangle,\\ |N,-\rangle&=-\sin(\theta_N/2)|e,N-1\rangle +\cos(\theta_N/2)|g,N\rangle,\\ E_{N,\pm}&=\hbar\omega_c(N-\tfrac12)\pm\hbar\Omega_N/2. \end{aligned}

At resonance the states have equal bare weights and their separation is 2ℏgN2\hbar g\sqrt N. The unpaired state ∣g,0⟩|g,0\rangle has energy −ℏω0/2-\hbar\omega_0/2; there is no ∣e,−1⟩|e,-1\rangle. The first doublet is split by 2ℏg2\hbar g even with no photon initially: ∣e,0⟩|e,0\rangle can exchange its excitation with ∣g,1⟩|g,1\rangle. This is the single-emitter vacuum Rabi splitting.

Bare atom–photon manifolds reorganized into Jaynes–Cummings dressed doublets

At atom–cavity resonance, the interaction couples ∣e,N−1⟩|e,N-1\rangle and ∣g,N⟩|g,N\rangle within each fixed-excitation manifold. Their symmetric and antisymmetric combinations form a doublet separated by 2ℏgN2\hbar g\sqrt N. The vacuum state ∣g,0⟩|g,0\rangle remains unpaired.

As Δc\Delta_c is swept through zero, every N≥1N\ge1 manifold has an avoided crossing. Its minimum gap is

ΔENmin⁡=2ℏgN.\Delta E_N^{\min} = 2\hbar g\sqrt N.

The bare content changes with detuning:

∣⟨e,N−1∣N,+⟩∣2=12(1+ΔcΩN).|\langle e,N-1|N,+\rangle|^2 = \frac12 \left( 1+\frac{\Delta_c}{\Omega_N} \right).

Thus the same spectral curve contains both a coupling measurement and a state-composition map. Measuring only the separation does not by itself determine which bare component a branch contains; the sign of detuning and the branch convention are also required.

For a physical quantized mode, the resonant state

∣N,+⟩=12(∣e,N−1⟩+∣g,N⟩)|N,+\rangle = \frac1{\sqrt2} \left( |e,N-1\rangle + |g,N\rangle \right)

is not an atomic superposition multiplied by a field state. The two photon states are orthogonal, so the atom and mode are entangled.

Tracing over the field gives

ρatom(N,+)=cos⁡2θN2∣e⟩⟨e∣+sin⁡2θN2∣g⟩⟨g∣.\rho_{\mathrm{atom}}^{(N,+)} = \cos^2\frac{\theta_N}{2} |e\rangle\langle e| + \sin^2\frac{\theta_N}{2} |g\rangle\langle g|.

There is no off-diagonal atomic coherence in this reduced state, even though the joint state is pure. The entanglement entropy is

SN=−pNln⁡pN−(1−pN)ln⁡(1−pN),S_N = - p_N\ln p_N - (1-p_N)\ln(1-p_N),

where

pN=cos⁡2θN2.p_N = \cos^2\frac{\theta_N}{2}.

At resonance, pN=1/2p_N=1/2 and SN=ln⁡2S_N=\ln2. Far from resonance, one component dominates and the entropy tends to zero.

This is a sharp difference from a semiclassical dressed state

ce∣e⟩+cg∣g⟩,c_e|e\rangle+c_g|g\rangle,

which is a pure state of the atom because the classical source is not part of the Hilbert space. Replacing a quantized field by a classical amplitude can therefore reproduce mean coherent dynamics while discarding atom–field entanglement.

Let the mode begin in a coherent state ∣α⟩|\alpha\rangle with

nˉ=∣α∣2.\bar n = |\alpha|^2.

The Jaynes–Cummings interaction samples many excitation manifolds. Near the center of a large photon-number distribution,

2gN≃2gnˉ,2g\sqrt N \simeq 2g\sqrt{\bar n},

so the corresponding classical Rabi-frequency coefficient is

Ω≃2g∣α∣.\Omega \simeq 2g|\alpha|.

This is a local, large-occupation correspondence, not an operator identity. A coherent state has photon-number variance

Var⁡(n)=nˉ,\operatorname{Var}(n) = \bar n,

and therefore contains a spread of Rabi frequencies: 2gn+12g\sqrt{n+1} for an initially excited atom and 2gn2g\sqrt n for an initially ground-state atom. At short times and large nˉ\bar n, the fractional spread is small and a classical drive is accurate. At longer times, the components dephase and can later rephase, producing collapse and revival of the inversion. That photon-number-resolved behavior is absent from a deterministic single-frequency classical field.

The precise large-field limit must also specify what is held fixed. Taking ∣α∣→∞|\alpha|\to\infty at fixed gg makes the drive diverge. A controlled semiclassical scaling keeps g∣α∣g|\alpha| finite while suppressing relative field fluctuations.

For a classical periodic Hamiltonian,

H(t+T)=H(t),T=2πωL,H(t+T) = H(t), \qquad T = \frac{2\pi}{\omega_L},

Floquet states have the form

∣ψα(t)⟩=e−iεαt/ℏ∣uα(t)⟩,∣uα(t+T)⟩=∣uα(t)⟩.\begin{gathered} |\psi_\alpha(t)\rangle = e^{-i\varepsilon_\alpha t/\hbar} |u_\alpha(t)\rangle, \\ |u_\alpha(t+T)\rangle = |u_\alpha(t)\rangle. \end{gathered}

The quasienergy is defined modulo ℏωL\hbar\omega_L:

εα∼εα+mℏωL,m∈Z.\varepsilon_\alpha \sim \varepsilon_\alpha+m\hbar\omega_L, \qquad m\in\mathbb Z.

Expanding ∣uα(t)⟩|u_\alpha(t)\rangle in Fourier harmonics produces an extended matrix with repeated bare levels shifted by mℏωLm\hbar\omega_L. Near a one-photon resonance, retaining the nearly degenerate pair gives the same 2×22\times2 RWA matrix that defines the semiclassical dressed states.

This establishes a useful dictionary:

Quantized single mode Classical Floquet drive
∣g,n⟩\lvert g,n\rangle, ∣e,n−1⟩\lvert e,n-1\rangle bare state with Fourier index
physical photon occupation harmonic bookkeeping index
joint energy eigenvalue quasienergy modulo ℏωL\hbar\omega_L
atom–field entanglement possible no source subsystem in the model
number fluctuations are physical drive amplitude prescribed

The analogy is strongest for a highly occupied, nearly monochromatic mode and near-resonant dynamics over times for which back-action and photon-number dispersion are negligible.

Outside the RWA, a Floquet calculation can retain many harmonics and reveal multiphoton avoided crossings and the Bloch–Siegert shift. A fully quantized Rabi model can retain counter-rotating terms. These extensions are again related, but their state spaces and observables should not be conflated.

Eigenvectors are not observed directly. A measurement accesses transition frequencies, line strengths, time evolution, or correlations.

Suppose the control does not dress a third state ∣r⟩|r\rangle, and a weak probe couples it only to the bare state ∣e⟩|e\rangle. If

⟨r∣D∣g⟩=0,⟨r∣D∣e⟩=dre,\langle r|D|g\rangle = 0, \qquad \langle r|D|e\rangle = d_{re},

then the dressed-state matrix elements are

⟨r∣D∣+⟩=drecos⁡θ2,⟨r∣D∣−⟩=−dresin⁡θ2.\begin{aligned} \langle r|D|+\rangle &= d_{re} \cos\frac{\theta}{2}, \\ \langle r|D|-\rangle &= - d_{re} \sin\frac{\theta}{2}. \end{aligned}

The squared dipole matrix elements therefore have the ratio

S+:S−=cos⁡2θ2:sin⁡2θ2.S_+ : S_- = \cos^2\frac{\theta}{2} : \sin^2\frac{\theta}{2}.

These give relative integrated absorption strengths from a common populated ∣r⟩|r\rangle when other frequency prefactors are effectively equal. For emission from the dressed branches, their populations also enter. Peak heights further depend on linewidths. At resonance the two dipole strengths are equal; far away one vanishes.

The observed peaks are separated by dressed energy differences only when the relevant linewidths, probe back-action, populations, and interference paths are included consistently. For the full three-level spectrum, see Autler–Townes Splitting.

If a bare state is prepared suddenly, the two dressed components acquire a relative phase

Δφ(t)=ε+−ε−ℏt=ΩRt.\Delta\varphi(t) = \frac{ \varepsilon_+-\varepsilon_- }{ \hbar } t = \Omega_Rt.

Bare-state population oscillations are an interferometric measurement of that phase. The frequency reveals the dressed gap, while the oscillation contrast reveals the mixing angle and state preparation.

Spontaneous emission acts through an operator such as σ−\sigma_-. Written in the dressed basis, this operator has several matrix elements between neighboring atom–photon manifolds. Consequently, a strongly driven two-level emitter can radiate at more than one dressed transition frequency. In the resonant, well-resolved regime this underlies the Mollow triplet.

The radiative reservoir, linewidths, and fluorescence spectrum are not contained in the Hermitian two-state eigenproblem alone. They require the open-system treatment developed in Spontaneous Emission and Optical Bloch Equations.

In cavity QED, the quantized dressed doublets can be spectroscopically resolved because the atom repeatedly exchanges an excitation with a long-lived resonator mode. In the one-excitation sector, a useful non-Hermitian pole model is

Heffℏ=(ω0−iγ/2ggωc−iκ/2),\frac{H_{\mathrm{eff}}}{\hbar} = \begin{pmatrix} \omega_0-i\gamma/2 & g \\ g & \omega_c-i\kappa/2 \end{pmatrix},

where γ\gamma is the atomic population-decay rate into noncavity channels and κ\kappa is the cavity energy-decay rate in this convention. Pure dephasing is absent; the atomic optical coherence therefore decays at γ/2\gamma/2. The complex poles are

ω~±=ω0+ωc2−iγ+κ4±g2+14[Δc−iγ−κ2]2.\widetilde\omega_\pm = \frac{\omega_0+\omega_c}{2} - i\frac{\gamma+\kappa}{4} \pm \sqrt{ g^2 + \frac14 \left[ \Delta_c - i\frac{\gamma-\kappa}{2} \right]^2 }.

The real parts give pole frequencies and the imaginary parts give amplitude decay rates; measured spectral maxima need not coincide with the real parts. At exact resonance, the square root has a nonzero real part when

g>∣κ−γ∣4.g > \frac{|\kappa-\gamma|}{4}.

That algebraic pole-splitting condition is not, by itself, a universal criterion for visibly resolving two peaks. Resolvability depends on the combined widths, drive and detection port, background, and fitting model. Likewise, the cooperativity

C=4g2κγC = \frac{4g^2}{\kappa\gamma}

tests coherent coupling against dissipation but answers a different question from peak separation.

The same linear two-mode equations can describe coupled classical oscillators. A weak-excitation doublet alone therefore does not establish field quantization (Zhu and colleagues). Resolving the excitation-dependent Jaynes–Cummings ladder probes its anharmonicity (Fink and colleagues); nonclassical correlations provide complementary information.

The closed dressed-state spectrum explains the underlying level structure. Cavity QED develops the mode overlap, linewidth conventions, dissipative spectrum, Purcell limit, and measured transmission or reflection. The open-system cavity-QED map continues to conditional output records and trajectories.

The quantum Rabi Hamiltonian contains

ℏg(a+a†)(σ++σ−).\hbar g \left( a+a^\dagger \right) \left( \sigma_++\sigma_- \right).

Besides the Jaynes–Cummings exchange terms, it includes

aσ−anda†σ+.a\sigma_- \qquad\text{and}\qquad a^\dagger\sigma_+.

These terms do not conserve N\mathcal N. The independent two-dimensional manifolds disappear; the discrete parity σz(−1)a†a\sigma_z(-1)^{a^\dagger a} remains conserved in the ideal Rabi model. Perturbatively they produce shifts such as the Bloch–Siegert correction. In ultrastrong coupling, the exact eigenstates contain contributions from many bare excitation sectors and the phrase “one photon in the dressed ground state” becomes measurement-model dependent.

A real atom or molecule may have several nearby states coupled by the same field. Then the dressed Hamiltonian is larger than 2×22\times2, and dark states, multiple avoided crossings, tensor light shifts, and leakage can appear. The two-state mixing angle remains useful only after the neglected levels have been shown to be perturbative, or exactly decoupled by a symmetry or selection rule.

Near-degenerate magnetic or hyperfine manifolds are especially sensitive to polarization and selection rules. A measured splitting cannot be assigned to one Ω\Omega until the complete coupled subspace has been checked.

If the coupling varies with position,

Ω=Ω(r),\Omega = \Omega(\mathbf r),

then the dressed energies define position-dependent adiabatic potentials. Atomic motion can also drive transitions between local dressed branches. The relevant small parameter compares motional matrix elements of the position-dependent basis to the local dressed gap. Treating ε±(r)\varepsilon_\pm(\mathbf r) as ordinary scalar potentials without those derivative couplings can fail near nodes and avoided crossings. Even in a single adiabatic branch, projected kinetic energy generally contributes a geometric vector potential and a scalar correction, in addition to the local eigenvalue.

Hermitian dressed eigenstates are not generally eigenmodes of a Lindblad generator. Relaxation, dephasing, drive noise, and cavity loss change the steady state and broaden transitions. In weak dissipation, the Hamiltonian dressed basis is often a useful starting point for a secular master equation. Near degeneracies or for structured reservoirs, secularizing in that basis may erase important interference.

The phrase “dressed-state lifetime” must therefore identify:

  • the prepared branch;
  • the decay and dephasing channels;
  • whether the quoted quantity is a population lifetime, coherence time, or spectral linewidth;
  • the approximation used to transform the bath coupling operators.

Keep θ\theta fixed with a nonzero gap and vary ϕ\phi adiabatically from 0 to 2π2\pi. Find the geometric phase of ∣+⟩|+\rangle in the eigenvector convention above. Why can a time-dependent drive phase matter even though the instantaneous eigenvalues do not depend on it?

Solution

The chosen vector is periodic in ϕ\phi and

i⟨+∣∂ϕ+⟩=−sin⁡2(θ/2),γ+=−π(1−cos⁡θ)(mod2π).i\langle+|\partial_\phi+\rangle=-\sin^2(\theta/2), \qquad \gamma_+=-\pi(1-\cos\theta)\pmod{2\pi}.

It is minus half the solid angle traced by the effective-field direction. The lower branch gives the opposite phase modulo 2π2\pi. A phase variation rotates eigenvectors, producing geometric phase during adiabatic motion and potentially transitions otherwise. It cannot be treated as a sequence of unrelated static rephasings without their time-derivative term.

2. Recover Rabi oscillations from dressed phases

Section titled “2. Recover Rabi oscillations from dressed phases”

Take ϕ=0\phi=0 and prepare ∣g⟩|g\rangle suddenly at t=0t=0. Use the dressed basis to find the excited-state probability for constant Ω\Omega and Δ\Delta.

Solution

The inverse transformation gives

∣g⟩=sin⁡θ2∣+⟩+cos⁡θ2∣−⟩.|g\rangle = \sin\frac{\theta}{2}|+\rangle + \cos\frac{\theta}{2}|-\rangle.

After time tt,

∣ψ(t)⟩=sin⁡θ2e−iΩRt/2∣+⟩+cos⁡θ2e+iΩRt/2∣−⟩.|\psi(t)\rangle = \sin\frac{\theta}{2} e^{-i\Omega_Rt/2} |+\rangle + \cos\frac{\theta}{2} e^{+i\Omega_Rt/2} |-\rangle.

Projecting onto ∣e⟩|e\rangle gives

ce(t)=−2isin⁡θ2cos⁡θ2sin⁡ΩRt2=−isin⁡θsin⁡ΩRt2.\begin{aligned} c_e(t) &= - 2i \sin\frac{\theta}{2} \cos\frac{\theta}{2} \sin\frac{\Omega_Rt}{2} \\ &= - i\sin\theta \sin\frac{\Omega_Rt}{2}. \end{aligned}

Therefore

Pe(t)=Ω2ΩR2sin⁡2(ΩRt2).P_e(t) = \frac{\Omega^2}{\Omega_R^2} \sin^2 \left( \frac{\Omega_Rt}{2} \right).

The oscillation is interference between stationary dressed components. Detuning increases the phase-evolution frequency but reduces the mixing factor Ω2/ΩR2\Omega^2/\Omega_R^2.

For fixed Ω>0\Omega>0, determine the limiting bare character of ∣+⟩|+\rangle and ∣−⟩|-\rangle as Δ→+∞\Delta\to+\infty and Δ→−∞\Delta\to-\infty. Explain why “the excited dressed state” is an unsafe global label.

Solution

When Δ→+∞\Delta\to+\infty,

θ→0,\theta\to0,

so

∣+⟩→∣e⟩,∣−⟩→eiϕ∣g⟩.|+\rangle\to|e\rangle, \qquad |-\rangle\to e^{i\phi}|g\rangle.

When Δ→−∞\Delta\to-\infty,

θ→π,\theta\to\pi,

so

∣+⟩→eiϕ∣g⟩,∣−⟩→−∣e⟩.|+\rangle\to e^{i\phi}|g\rangle, \qquad |-\rangle\to-|e\rangle.

The upper and lower eigenvalue branches are continuous, but each exchanges its bare identity through the avoided crossing. A label based on bare character must therefore include the side of resonance or an explicit adiabatic path.

Let Ω\Omega be constant and sweep

Δ(t)=vt.\Delta(t) = vt.

Find ∣θ˙∣|\dot\theta|, locate its maximum, and obtain a simple local adiabaticity condition. For

Ω2π=100 kHz,∣v∣2π=1.0×1010 Hz s−1,\begin{gathered} \dfrac{\Omega}{2\pi} = 100\ \mathrm{kHz}, \\ \dfrac{|v|}{2\pi} = 1.0\times10^{10}\ \mathrm{Hz\,s^{-1}}, \end{gathered}

estimate the Landau–Zener diabatic probability. Convert all quantities consistently to angular units.

Solution

For constant Ω\Omega,

∣θ˙∣=Ω∣v∣Δ2+Ω2.|\dot\theta| = \frac{ \Omega|v| }{ \Delta^2+\Omega^2 }.

It is largest at resonance:

∣θ˙∣max⁡=∣v∣Ω.|\dot\theta|_{\max} = \frac{|v|}{\Omega}.

Comparing this with the minimum gap frequency Ω\Omega gives

∣v∣≪Ω2.|v| \ll \Omega^2.

Numerically,

Ω=2π×105 s−1,\Omega = 2\pi\times10^5\ \mathrm{s^{-1}},

and

∣v∣=2π×1010 s−2.|v| = 2\pi\times10^{10}\ \mathrm{s^{-2}}.

The exponent is

πΩ22∣v∣=π(2π×105)22(2π×1010)=π2≃9.87.\begin{aligned} \frac{\pi\Omega^2}{2|v|} &= \frac{ \pi(2\pi\times10^5)^2 }{ 2(2\pi\times10^{10}) } \\ &= \pi^2 \simeq 9.87. \end{aligned}

Thus

Pdiab≃e−9.87≃5.2×10−5.P_{\mathrm{diab}} \simeq e^{-9.87} \simeq 5.2\times10^{-5}.

The transfer is highly, though not perfectly, adiabatic in the ideal infinite-sweep model.

5. The detected channel selects different lines

Section titled “5. The detected channel selects different lines”

At resonance choose ∣N,s⟩=(∣e,N−1⟩+s∣g,N⟩)/2|N,s\rangle=(|e,N-1\rangle+s|g,N\rangle)/\sqrt2, where s=±1s=\pm1. For N≥2N\ge2, calculate transitions to ∣N−1,t⟩|N-1,t\rangle through atomic emission σ−\sigma_- and through cavity leakage aa. Compare their relative strengths for large NN.

Solution

Apply the lowering operator before taking the final overlap:

⟨N−1,t∣σ−∣N,s⟩=t/2,⟨N−1,t∣a∣N,s⟩=N−1+stN2.\begin{aligned} \langle N-1,t|\sigma_-|N,s\rangle&=t/2,\\ \langle N-1,t|a|N,s\rangle&= \frac{\sqrt{N-1}+st\sqrt N}{2}. \end{aligned}

All four atomic-channel matrix elements have squared magnitude 1/41/4. For cavity leakage, same-branch strengths approach NN, whereas opposite-branch strengths approach 1/(16N)1/(16N). The ratio is therefore about 16N216N^2. Actual intensities also require initial populations and the channel’s coupling to the detector. The N=1N=1 transition ends in the unpaired vacuum state and must be calculated separately.

For the state ∣N,+⟩|N,+\rangle, compute the atomic purity

P=Tr⁡[(ρatom(N,+))2].\mathcal P = \operatorname{Tr} \left[ \left( \rho_{\mathrm{atom}}^{(N,+)} \right)^2 \right].

Evaluate it on resonance and in the limit ∣Δc∣≫gN|\Delta_c|\gg g\sqrt N.

Solution

The reduced atomic probabilities are

p=cos⁡2θN2,1−p=sin⁡2θN2.p = \cos^2\frac{\theta_N}{2}, \qquad 1-p = \sin^2\frac{\theta_N}{2}.

Hence

P=p2+(1−p)2.\mathcal P = p^2+(1-p)^2.

Using

2p−1=cos⁡θN=ΔcΩN,2p-1 = \cos\theta_N = \frac{\Delta_c}{\Omega_N},

one finds

P=12[1+Δc2ΩN2].\mathcal P = \frac12 \left[ 1+ \frac{\Delta_c^2}{\Omega_N^2} \right].

At resonance,

P=12,\mathcal P = \frac12,

the minimum purity for a two-level reduced state, corresponding to maximal atom–mode entanglement in this two-term Schmidt decomposition.

For ∣Δc∣≫gN|\Delta_c|\gg g\sqrt N,

Δc2ΩN2→1,\frac{\Delta_c^2}{\Omega_N^2} \to 1,

so P→1\mathcal P\to1. The dressed state approaches a bare product state.

Prepare atoms in a third state ∣r⟩|r\rangle that couples only to ∣e⟩|e\rangle. Find the relative squared dipole matrix elements for absorption to ∣+⟩|+\rangle and ∣−⟩|-\rangle for

Δ=34ΩR.\Delta = \frac34\Omega_R.

What additional information is needed before claiming that two peaks will be visible?

Solution

The condition gives

cos⁡θ=34.\cos\theta = \frac34.

Therefore

cos⁡2θ2=1+cos⁡θ2=78,\cos^2\frac{\theta}{2} = \frac{1+\cos\theta}{2} = \frac78,

and

sin⁡2θ2=18.\sin^2\frac{\theta}{2} = \frac18.

For these dipole strengths,

S+:S−=7:1.S_+:S_- = 7:1.

Two eigenvalues do not guarantee two visible peaks. One must also know the dressed separation, homogeneous and inhomogeneous linewidths, dressed-state populations, probe strength, competing matrix elements, interference paths, instrument response, and signal-to-noise ratio.

8. Separate Floquet harmonics from photons

Section titled “8. Separate Floquet harmonics from photons”

Consider a two-level atom driven by a prescribed classical periodic field. A Floquet calculation labels basis vectors by an integer mm. Decide which of the following statements are valid, and justify each answer:

  1. measuring the cavity photon number measures mm;
  2. quasienergies differing by ℏωL\hbar\omega_L represent equivalent Floquet zones;
  3. a Floquet state necessarily entangles the atom with the source;
  4. a highly occupied coherent quantized mode can approximate the semiclassical Rabi coefficient.
Solution
  1. Invalid. A prescribed classical drive has no cavity Hilbert space or photon-number operator. The integer mm labels a Fourier harmonic in extended Floquet space.
  2. Valid. Multiplying the periodic part by eimωLte^{im\omega_Lt} shifts the quasienergy by mℏωLm\hbar\omega_L without changing the physical solution.
  3. Invalid. The classical source is not a subsystem of the model, so atom–source entanglement is undefined there. The atomic state may still be a coherent superposition.
  4. Valid with qualifications. For a coherent state with large nˉ\bar n, the matrix element near the distribution center gives Ω≃2gnˉ=2g∣α∣\Omega\simeq2g\sqrt{\bar n}=2g|\alpha|. The approximation is limited by photon-number dispersion, back-action, and the timescale of interest.