Skip to content

Effective Hamiltonians and Scale Separation

An effective Hamiltonian is a reduced generator designed to reproduce specified predictions of a more detailed quantum model within a stated regime and accuracy. It may act on a low-energy subspace, on slow degrees of freedom, in a rotating frame, or only at stroboscopic times. Its value is not that it is smaller. Its value is that the discarded physics has been encoded systematically in corrected energies, induced couplings, dressed observables, or controlled error terms.

The word “effective” is incomplete unless four questions have answers:

  1. What is retained? A subspace, band, manifold, slow coordinate, or set of resonant states.
  2. What is reproduced? Energies, transition amplitudes, scattering poles, slow dynamics, or one-period evolution.
  3. In which regime? Below an energy cutoff, away from resonances, for weak coupling, at high drive frequency, or over a specified time window.
  4. To what accuracy? An order in a dimensionless ratio, a norm estimate, a leakage probability, or a benchmark against the full model.

This chapter organizes the main constructions. Projection Methods owns exact P/Q block reduction and state reconstruction, Feshbach Projection Formalism owns the continuum and resonance extension of energy-dependent elimination, Schrieffer–Wolff Transformation owns perturbative unitary block diagonalization, Rotating-Wave Approximation owns near-resonant time averaging, and Magnus Expansion together with Floquet-Magnus Expansion owns exponential and periodic effective generators.

Folded Effective Hamiltonians owns the Q-box derivative and folded-diagram route from an energy-dependent model-space equation to one energy-independent interaction.

Adiabatic Elimination owns the dynamical slow–fast reduction, memory and derivative corrections, the detuned Λ system, and trace-preserving effective operators for decaying fast sectors.

Born–Oppenheimer Approximation as Scale Separation owns the heavy–light mass expansion, exact electronic-channel equations, potential-energy surfaces, and diagnostics for nonadiabatic failure.

Average Hamiltonian Theory owns periodic control cycles, toggling-frame interval averages, pulse-order corrections, and selective interaction engineering.

A full Hilbert space split into retained P and eliminated Q sectors, followed by four controlled routes to a reduced effective Hamiltonian.

Scale separation can be spectral, dynamical, or temporal. Projection with a resolvent, unitary block diagonalization, adiabatic elimination, and time averaging answer different questions, but each must carry the retained sector, expansion parameter, transformed observables, and error window into the reduced model.

Let Hfull\mathcal H_{\mathrm{full}} be the Hilbert space of the detailed model and HP\mathcal H_P the retained model space. An encoding isometry

E:HP⟶Hfull\mathcal E: \mathcal H_P \longrightarrow \mathcal H_{\mathrm{full}}

identifies reduced states with physical states. A useful effective Hamiltonian satisfies a relation of the form

E†e−iHt/ℏE≃e−iHefft/ℏ\mathcal E^\dagger e^{-iHt/\hbar} \mathcal E \simeq e^{-iH_{\mathrm{eff}}t/\hbar}

for the states, observables, and times of interest. The symbol ≃\simeq hides the entire approximation problem. It may mean equality of selected eigenvalues through order ηk\eta^k, agreement of matrix elements up to an error ε\varepsilon, or equality only at integer multiples of a drive period.

No single reduced operator is required to reproduce every property of the full theory. Common targets include:

  • the low-lying spectrum and level splittings;
  • transition amplitudes among a selected group of states;
  • projected Green functions and resonance poles;
  • slow evolution after fast transients have decayed;
  • dynamics in a rotating frame near a chosen resonance;
  • one-period or stroboscopic evolution in a driven system.

The target determines the construction. A Hamiltonian optimized for low-energy eigenvalues need not reproduce short-time leakage. A stroboscopic Floquet Hamiltonian does not by itself describe micromotion inside each period. An energy-dependent optical Hamiltonian can reproduce a projected scattering resolvent without being a conventional Hermitian operator with a fixed spectrum.

Simply replacing HH by its upper-left block can miss the leading physics. States outside the retained space may be only virtually occupied, yet excursions through them can shift energies and mediate interactions. Effective theory keeps those consequences while avoiding explicit evolution of the remote states.

For this reason,

Heff≠PHPH_{\mathrm{eff}} \ne PHP

in general. The projected Hamiltonian PHPPHP is only the zeroth approximation. Corrections record how the eliminated sector reacts to the retained one.

Choose orthogonal projectors

P2=P,Q2=Q,P+Q=I,PQ=0.\begin{aligned} P^2=P, &\qquad Q^2=Q, \\ P+Q=I, &\qquad PQ=0. \end{aligned}

In the decomposition

H=PH⊕QH,\mathcal H = P\mathcal H \oplus Q\mathcal H,

the Hamiltonian has block form

H=(HPPHPQHQPHQQ),H = \begin{pmatrix} H_{PP} & H_{PQ}\\ H_{QP} & H_{QQ} \end{pmatrix},

with HPP=PHPH_{PP}=PHP and similarly for the other blocks.

Suppose the retained states occupy an energy window separated from the unwanted states by a gap or detuning Δ\Delta, while the off-diagonal coupling has characteristic size vv. The natural expansion parameter is

η=vΔ.\eta = \frac{v}{\Delta}.

When η≪1\eta\ll1, leakage amplitudes are often of order η\eta, leakage probabilities of order η2\eta^2, and virtual energy shifts of order

δE∼v2Δ.\delta E \sim \frac{v^2}{\Delta}.

These scalings explain why “unoccupied” levels can matter. Their real population is small, but the phase accumulated from virtual excursions can be measurable over long times.

The retained space should contain every state that mixes strongly on the scale being resolved. If two levels have detuning δ\delta comparable with their coupling vv, placing one in PP and one in QQ creates a denominator of order δ\delta and destroys the expansion. Both belong in the model space and should be diagonalized together.

This is the central rule of Quasi-Degenerate Perturbation Theory:

Model-space rule: treat nearby states exactly; eliminate remote states perturbatively.

The boundary of PP is therefore physical and resolution dependent. A two-level model can be excellent for microwave dynamics and inadequate for a faster pulse that resolves leakage levels.

For an exact eigenstate,

H∣Ψ⟩=E∣Ψ⟩,H\lvert\Psi\rangle = E\lvert\Psi\rangle,

write

∣ψP⟩=P∣Ψ⟩,∣ψQ⟩=Q∣Ψ⟩.\lvert\psi_P\rangle = P\lvert\Psi\rangle, \qquad \lvert\psi_Q\rangle = Q\lvert\Psi\rangle.

Projection gives

HPP∣ψP⟩+HPQ∣ψQ⟩=E∣ψP⟩,HQP∣ψP⟩+HQQ∣ψQ⟩=E∣ψQ⟩.\begin{aligned} H_{PP}\lvert\psi_P\rangle + H_{PQ}\lvert\psi_Q\rangle &= E\lvert\psi_P\rangle, \\ H_{QP}\lvert\psi_P\rangle + H_{QQ}\lvert\psi_Q\rangle &= E\lvert\psi_Q\rangle. \end{aligned}

Whenever the required inverse exists with the appropriate boundary condition,

∣ψQ⟩=(E−HQQ)−1HQP∣ψP⟩.\lvert\psi_Q\rangle = \left( E-H_{QQ} \right)^{-1} H_{QP}\lvert\psi_P\rangle.

Substitution yields the exact projected equation

Heff(E)∣ψP⟩=E∣ψP⟩,H_{\mathrm{eff}}(E) \lvert\psi_P\rangle = E\lvert\psi_P\rangle,

with

Heff(E)=HPP+HPQ(E−HQQ)−1HQP.\begin{aligned} H_{\mathrm{eff}}(E) &= H_{PP} \\ &\quad+ H_{PQ} \left( E-H_{QQ} \right)^{-1} H_{QP}. \end{aligned}

The second term is a self-energy. It contains virtual propagation through QQ. This formula is exact but nonlinear in EE, and singularities of the QQ-space resolvent can invalidate a naive expansion.

For a continuum with outgoing boundary conditions, the resolvent can acquire an imaginary part. The projected operator then becomes non-Hermitian and describes shifts together with resonance widths. That is a feature of the boundary-value problem, not a loss of unitarity in the full closed theory.

The exact block reduction, Schur-complement identity, state reconstruction, and projected normalization are developed in Projection Methods. Continuum boundary conditions, resonance interpretation, and scattering applications belong to Feshbach Projection Formalism.

An alternative is to seek a unitary transformation that decouples PP and QQ to a chosen order. Let S†=−SS^\dagger=-S and define

H~=eSHe−S.\widetilde H = e^SHe^{-S}.

The generator is chosen so that

PH~Q=O(ηk+1)P\widetilde H Q = O(\eta^{k+1})

after working through order ηk\eta^k. The effective Hamiltonian is then

Heff=PH~P.H_{\mathrm{eff}} = P\widetilde H P.

Expanding the exponential produces nested commutators:

H~=H+[S,H]+12[S,[S,H]]+⋯ .\widetilde H = H + [S,H] + \frac12[S,[S,H]] + \cdots.

Unlike the exact resolvent form, a Schrieffer–Wolff Hamiltonian can be made energy independent and Hermitian order by order for an isolated Hermitian problem. It is especially useful when one wants a reusable low-energy operator rather than a nonlinear eigenvalue equation.

Different block-diagonalization conventions can produce different-looking Hamiltonians related by a unitary transformation inside PP. Their matrix entries are not individually observable. Spectra and consistently transformed matrix elements are.

Scale separation can be dynamical rather than purely spectral. Consider slow amplitudes csc_s coupled to a fast, far-detuned amplitude cfc_f:

ic˙s=hsscs+vsfcf,ic˙f=Δcf+vfscs.\begin{aligned} i\dot c_s &= h_{ss}c_s + v_{sf}c_f, \\ i\dot c_f &= \Delta c_f + v_{fs}c_s. \end{aligned}

If ∣Δ∣\lvert\Delta\rvert is the dominant frequency and csc_s changes little during a time ∣Δ∣−1\lvert\Delta\rvert^{-1}, the fast amplitude follows approximately:

cf≃−vfsΔcs.c_f \simeq - \frac{v_{fs}}{\Delta} c_s.

Substitution gives

ic˙s≃(hss−vsfvfsΔ)cs.i\dot c_s \simeq \left( h_{ss} - \frac{v_{sf}v_{fs}}{\Delta} \right)c_s.

This is adiabatic elimination. The correction is the same virtual second-order process seen in projector and unitary languages. What differs is the organizing assumption: a fast amplitude is slaved to slow motion.

Setting c˙f=0\dot c_f=0 is only the first term of an expansion. Adiabatic Elimination derives the exact transient and memory kernel, local derivative corrections, validity conditions, and the open-system effective jumps. Near resonance, Δ\Delta is no longer fast and the eliminated state must be restored.

The Born–Oppenheimer approximation is a richer version of slow–fast separation. Electronic states are solved at fixed nuclear coordinates, generating potential-energy surfaces, while derivative couplings quantify the failure of perfect separation. A level crossing or small electronic gap can make those nonadiabatic couplings large.

Two long-lived states can couple through a far-detuned excited state even when that excited state remains weakly populated. Eliminating the excited amplitude generates diagonal ac Stark shifts and an off-diagonal Raman coupling with the scale

ΩR∼Ω1∗Ω2Δ.\Omega_{\mathrm R} \sim \frac{\Omega_1^*\Omega_2}{\Delta}.

The induced coupling and light shifts arise at the same perturbative order, while the excited-state population is suppressed by a factor of order ∣Ω/Δ∣2\lvert\Omega/\Delta\rvert^2. The full convention-sensitive derivation, differential light shift, bright and dark states, and spontaneous-scattering extension belong to Adiabatic Elimination.

The same leading virtual process can be organized as degenerate perturbation theory, a resolvent excursion, Schrieffer–Wolff block diagonalization, or adiabatic elimination. Their distinctions matter in higher-order normalization, observable reconstruction, time dependence, and dissipation.

A time-dependent problem introduces another kind of scale separation: fast oscillations can average away while slow resonant dynamics survives.

For a unitary frame transformation R(t)R(t),

HR(t)=R†(t)H(t)R(t)−iℏR†(t)R˙(t).H_R(t) = R^\dagger(t)H(t)R(t) - i\hbar R^\dagger(t)\dot R(t).

This formula is exact. An approximation begins only when terms in HR(t)H_R(t) are discarded or averaged. Confusing the frame change with the approximation hides the control parameter.

In a rotating frame, resonant terms can become slowly varying while counter-rotating terms oscillate near twice the carrier frequency. The rotating-wave approximation drops those fast terms when the coupling and detuning are small compared with the fast frequency and when the observation window does not resolve the omitted micromotion.

The leading neglected effect is often not zero. It can include the Bloch–Siegert shift, leakage, or corrections to the rotation axis. Rotating-Wave Approximation develops those conditions quantitatively.

For general time dependence, the Magnus Expansion writes evolution as

U(t,0)=exp⁡[Ω(t)],U(t,0) = \exp\left[ \Omega(t) \right],

where Ω(t)\Omega(t) is built from time integrals and nested commutators. Over an interval TT, one may define

Heff(T)=iℏTΩ(T).H_{\mathrm{eff}}^{(T)} = \frac{i\hbar}{T} \Omega(T).

For a cyclic pulse sequence, Average Hamiltonian Theory first transforms the internal Hamiltonian into the control toggling frame. Its zeroth-order weighted average selects the retained interaction, while ordered commutators quantify pulse-order and finite-cycle corrections.

For periodic driving, Floquet–Magnus Expansion applies the Magnus series directly to a one-period propagator at a chosen drive phase:

U(T,0)=e−iHFT/ℏ.U(T,0) = e^{-iH_FT/\hbar}.

The Floquet Hamiltonian HFH_F controls stroboscopic evolution. A separate micromotion operator is needed inside each period. High-Frequency Expansions owns the broader inverse-frequency framework: Floquet-space block diagonalization, the phase-independent van Vleck Hamiltonian, kick-operator dressing, resonant denominators, and local many-body prethermal regimes. Its simplest expansion parameter is schematically

ηΩ∼JℏΩ,\eta_\Omega \sim \frac{J}{\hbar\Omega},

where JJ is a local energy or coupling scale. Resonances and many-body heating can limit the useful time window even when early terms are small.

An effective Hamiltonian alone is not a complete reduction. Suppose the transformed state is

∣Ψ~⟩=eS∣Ψ⟩.\lvert\widetilde\Psi\rangle = e^S\lvert\Psi\rangle.

If ∣ψeff⟩=P∣Ψ~⟩\lvert\psi_{\mathrm{eff}}\rangle=P\lvert\widetilde\Psi\rangle, then reconstructing the physical state requires

∣Ψ⟩≃e−S∣ψeff⟩,\lvert\Psi\rangle \simeq e^{-S}\lvert\psi_{\mathrm{eff}}\rangle,

to the same perturbative order. An observable OO must be transformed consistently:

Oeff=PeSOe−SP.O_{\mathrm{eff}} = P e^S O e^{-S}P.

Using HeffH_{\mathrm{eff}} with the bare projected observable POPPOP can give the correct energy spectrum and incorrect transition strengths. The missing terms represent the small admixture of eliminated states.

The same warning appears in other languages:

  • Feshbach elimination reconstructs ∣ψQ⟩\lvert\psi_Q\rangle with the resolvent.
  • Adiabatic elimination includes derivative and initial-slip corrections.
  • Floquet theory requires micromotion dressing for observables away from stroboscopic times.
  • Open-system elimination generally modifies both the Hamiltonian and the effective jump operators.

Every reduction should report more than the formal Hamiltonian.

Identify a dimensionless ratio such as

vΔ,τfastτslow,JℏΩ.\frac{v}{\Delta}, \qquad \frac{\tau_{\mathrm{fast}}}{\tau_{\mathrm{slow}}}, \qquad \frac{J}{\hbar\Omega}.

If several gaps or frequencies exist, the smallest relevant denominator controls the worst case.

Monitor

L(t)=⟨Ψ(t)∣Q∣Ψ(t)⟩.L(t) = \langle\Psi(t)\rvert Q \lvert\Psi(t)\rangle.

Small leakage supports a subspace description, but it is not sufficient by itself. Virtual phase shifts can accumulate while L(t)L(t) remains small.

An operator error δH=Hexact−Heff\delta H=H_{\mathrm{exact}}-H_{\mathrm{eff}} can produce an evolution error that grows with time. Duhamel’s identity gives, for bounded Hermitian generators,

∥e−iHexactt/ℏ−e−iHefft/ℏ∥≤∣t∣ℏ∥δH∥.\left\| e^{-iH_{\mathrm{exact}}t/\hbar} - e^{-iH_{\mathrm{eff}}t/\hbar} \right\| \le \frac{\lvert t\rvert}{\hbar} \lVert\delta H\rVert.

Thus a small energy error does not justify arbitrary evolution times. If ∥δH∥=O(ηk+1)\lVert\delta H\rVert=O(\eta^{k+1}), the naive useful time can shrink as the desired phase accuracy becomes more demanding.

Inspect denominators explicitly. A drive, threshold, avoided crossing, or continuum pole can make a nominally remote sector resonant. Near such a point:

  • enlarge the retained space;
  • change to a resonant frame;
  • retain the relevant channel exactly;
  • or abandon the expansion in favor of direct numerical evolution.

For a finite model, compare low-energy eigenvalues and eigenvectors with direct diagonalization. For dynamics, compare projected populations, phases, and observables over the intended time window. Vary the cutoff, expansion order, and retained space. A single matching eigenvalue is not enough.

Effective Hamiltonians are not unique. If WW is unitary inside the retained space,

Heff′=WHeffW†,Oeff′=WOeffW†,\begin{aligned} H_{\mathrm{eff}}' &= W H_{\mathrm{eff}}W^\dagger, \\ O_{\mathrm{eff}}' &= W O_{\mathrm{eff}}W^\dagger, \end{aligned}

then all consistently computed predictions agree. Time-dependent frame choices, Floquet gauges, and different perturbative generators create similar representational freedom.

The reliable objects are matched observables, not isolated coefficients. A practical matching procedure is:

  1. choose the retained degrees of freedom and operator basis;
  2. identify the expansion parameter and symmetries;
  3. calculate selected full-model amplitudes or energies;
  4. choose effective coefficients so the reduced model reproduces them;
  5. estimate the first omitted operators and verify against new observables.

This is the quantum-mechanical precursor of effective field theory. The QFT Bridge: EFT and Effective Hamiltonians explains what changes when infinitely many modes, locality, renormalization, and power counting enter.

Physical separationNatural methodTypical outputPrincipal warning
Chosen subspace coupled to remote statesProjection or Feshbach methodEnergy-dependent self-energyResolvent poles and boundary conditions
Several model-space states need one reusable interactionFolded effective HamiltonianEnergy-independent model-space operatorQ-box poles, root selection, and induced many-body terms
Low-energy and high-energy blocks with weak mixingSchrieffer–Wolff transformationHermitian energy-independent low-energy HamiltonianTransform states and observables
Slow amplitudes coupled to fast amplitudesAdiabatic eliminationAlgebraic slow-sector generatorInitial transients and derivative corrections
Heavy and light coordinatesBorn–Oppenheimer methodPotential-energy surfacesSmall gaps and nonadiabatic couplings
Near-resonant and fast counter-rotating termsRotating-wave approximationSlowly varying rotating-frame HamiltonianBloch–Siegert shifts and strong-drive failure
Rapid cyclic pulse controlAverage Hamiltonian theoryToggling-frame cycle HamiltonianCommutator corrections and finite pulses
General time-ordered evolutionMagnus expansionSingle exponential generatorConvergence and long-time accumulation
One-period logarithm at a chosen drive phaseFloquet–Magnus expansionStroboscopic HamiltonianBranch, phase, and micromotion conventions
Off-resonant rapid periodic driveHigh-frequency expansionsEffective Hamiltonian plus periodic kickResonances, asymptotic truncation, and heating
Weakly occupied decaying sectorEffective open-system operatorsEffective Hamiltonian and jump operatorsEliminating only the Hamiltonian is incomplete

These methods can overlap. The same three-level atom can be treated by projection, block diagonalization, or adiabatic elimination. Prefer the formulation whose assumptions and target observable are easiest to state and verify.

  • Elementary quantum mechanics: nearly degenerate levels, double-well doublets, spin–orbit reductions, and low-energy scattering channels.
  • Atomic, molecular, and optical physics: Raman transitions, light shifts, dark-state manifolds, rotating frames, cavity-mediated interactions, and Born–Oppenheimer surfaces.
  • Quantum information: qubit subspaces, leakage corrections, dispersive gates, and perturbative gadgets, together with pulse averaging.
  • Quantum matter: the exact Hubbard-dimer superexchange benchmark and effective masses, many-body lattice and impurity reductions, projected bands, and high-frequency Hamiltonian engineering.
  • Nuclear physics: model-space interactions, optical potentials, folded Hamiltonians, and reaction channels.
  • Open quantum systems: eliminated excited states can induce both coherent shifts and dissipative jump processes.
  • Quantum field theory: integrating out heavy modes produces symmetry-allowed effective operators whose coefficients are fixed by matching.

The common logic is scale aware rather than field specific: retain what the experiment resolves, encode what it does not, and carry an error estimate.

  • Calling PHPPHP the effective Hamiltonian without checking virtual QQ-space corrections.
  • Eliminating a state whose detuning is comparable with its coupling.
  • Quoting a small coupling without dividing by the relevant gap or frequency.
  • Matching eigenvalues while leaving states and observables undressed.
  • Treating an energy-dependent Feshbach operator as an ordinary fixed Hermitian Hamiltonian.
  • Setting a fast derivative to zero without estimating derivative and initial-transient corrections.
  • Using a rotating-wave or high-frequency Hamiltonian without specifying the frame or observation times.
  • Ignoring micromotion because stroboscopic evolution is accurate.
  • Assuming small leakage guarantees small phase error.
  • Extrapolating a finite-order Hamiltonian to arbitrarily long times.
  • Comparing coefficients from two effective Hamiltonians without checking whether they are related by a unitary change of basis.
  • Forgetting that open-system elimination changes dissipative operators as well as coherent dynamics.

Starting from the block eigenvalue equation, eliminate ∣ψQ⟩\lvert\psi_Q\rangle and derive Heff(E)H_{\mathrm{eff}}(E).

Solution

The projected equations are

HPP∣ψP⟩+HPQ∣ψQ⟩=E∣ψP⟩,HQP∣ψP⟩+HQQ∣ψQ⟩=E∣ψQ⟩.\begin{aligned} H_{PP}\lvert\psi_P\rangle + H_{PQ}\lvert\psi_Q\rangle &= E\lvert\psi_P\rangle, \\ H_{QP}\lvert\psi_P\rangle + H_{QQ}\lvert\psi_Q\rangle &= E\lvert\psi_Q\rangle. \end{aligned}

Rearranging the second equation gives

(E−HQQ)∣ψQ⟩=HQP∣ψP⟩.\left( E-H_{QQ} \right)\lvert\psi_Q\rangle = H_{QP}\lvert\psi_P\rangle.

Assuming the required resolvent exists,

∣ψQ⟩=(E−HQQ)−1HQP∣ψP⟩.\lvert\psi_Q\rangle = \left( E-H_{QQ} \right)^{-1} H_{QP}\lvert\psi_P\rangle.

Substitution into the first equation yields

Heff(E)≡HPP+HPQRQ(E)HQP,RQ(E)≡(E−HQQ)−1,Heff(E)∣ψP⟩=E∣ψP⟩.\begin{aligned} H_{\mathrm{eff}}(E) &\equiv H_{PP} \\ &\quad+ H_{PQ} R_Q(E) H_{QP}, \\ R_Q(E) &\equiv \left( E-H_{QQ} \right)^{-1}, \\ H_{\mathrm{eff}}(E) \lvert\psi_P\rangle &= E\lvert\psi_P\rangle. \end{aligned}

The first two rows define the effective Hamiltonian and the QQ-space resolvent.

Let

H=(0vv∗Δ),∣v∣≪∣Δ∣.H = \begin{pmatrix} 0 & v\\ v^* & \Delta \end{pmatrix}, \qquad \lvert v\rvert\ll\lvert\Delta\rvert.

Take the first basis state as PP. Find the exact energy-dependent effective Hamiltonian and the low-energy eigenvalue through second order.

Solution

Here

HPP=0,HQQ=Δ,HPQ=v,HQP=v∗.\begin{aligned} H_{PP}=0, &\qquad H_{QQ}=\Delta, \\ H_{PQ}=v, &\qquad H_{QP}=v^*. \end{aligned}

Therefore

Heff(E)=∣v∣2E−Δ.H_{\mathrm{eff}}(E) = \frac{\lvert v\rvert^2}{E-\Delta}.

The eigenvalue equation is

E=∣v∣2E−Δ,E = \frac{\lvert v\rvert^2}{E-\Delta},

or

E2−ΔE−∣v∣2=0.E^2-\Delta E-\lvert v\rvert^2=0.

For the root near zero, expand the denominator using E/Δ=O(∣v/Δ∣2)E/\Delta=O(\lvert v/\Delta\rvert^2):

Elow=−∣v∣2Δ+O(∣v∣4Δ3).E_{\mathrm{low}} = - \frac{\lvert v\rvert^2}{\Delta} + O\left( \frac{\lvert v\rvert^4}{\Delta^3} \right).

The sign reverses if the eliminated state lies below rather than above the retained state.

A three-level Λ system has Rabi frequencies Ω1\Omega_1 and Ω2\Omega_2, one-photon detuning Δ\Delta, two-photon detuning δ\delta, and excited-state linewidth Γ\Gamma. State the dimensionless ratios that must be small for a coherent two-state reduction, and identify what must be added if spontaneous emission cannot be neglected.

Solution

At minimum,

∣Ω1∣∣Δ∣,∣Ω2∣∣Δ∣,∣δ∣∣Δ∣≪1.\frac{\lvert\Omega_1\rvert}{\lvert\Delta\rvert}, \quad \frac{\lvert\Omega_2\rvert}{\lvert\Delta\rvert}, \quad \frac{\lvert\delta\rvert}{\lvert\Delta\rvert} \ll 1.

Drive-envelope frequencies must also be small compared with ∣Δ∣\lvert\Delta\rvert, and the intended duration must not amplify neglected phase errors. If Γ\Gamma matters, the fast response uses the complex scale Δ−iΓ/2\Delta-i\Gamma/2, but a non-Hermitian Hamiltonian alone is incomplete. Effective Lindblad jump operators must be retained to describe scattering and preserve trace. The canonical derivation is in Adiabatic Elimination.

Let H1H_1 and H2H_2 be bounded Hermitian operators. Use Duhamel’s identity to show

∥e−iH1t/ℏ−e−iH2t/ℏ∥≤∣t∣ℏ∥H1−H2∥.\left\| e^{-iH_1t/\hbar} - e^{-iH_2t/\hbar} \right\| \le \frac{\lvert t\rvert}{\hbar} \lVert H_1-H_2\rVert.
Solution

Define

F(s)=e−iH1(t−s)/ℏe−iH2s/ℏ.F(s) = e^{-iH_1(t-s)/\hbar} e^{-iH_2s/\hbar}.

Then

dFds=iℏe−iH1(t−s)/ℏ(H1−H2)e−iH2s/ℏ.\frac{dF}{ds} = \frac{i}{\hbar} e^{-iH_1(t-s)/\hbar} \left( H_1-H_2 \right) e^{-iH_2s/\hbar}.

Integrating from 00 to tt gives the difference of the two propagators. Unitary invariance of the operator norm implies

∥dFds∥≤1ℏ∥H1−H2∥.\left\| \frac{dF}{ds} \right\| \le \frac{1}{\hbar} \lVert H_1-H_2\rVert.

The triangle inequality for the integral then yields the stated bound. It makes explicit why a small generator error can accumulate over long times.

Let S=O(η)S=O(\eta) be anti-Hermitian. Expand

Oeff=PeSOe−SPO_{\mathrm{eff}} = P e^S O e^{-S}P

through first order. Under what condition is POPPOP sufficient to that order?

Solution

The Baker–Campbell–Hausdorff expansion gives

eSOe−S=O+[S,O]+O(η2).e^S O e^{-S} = O+[S,O]+O(\eta^2).

Therefore

Oeff=POP+P[S,O]P+O(η2).O_{\mathrm{eff}} = POP + P[S,O]P + O(\eta^2).

The bare projection POPPOP is sufficient through first order only if

P[S,O]P=0.P[S,O]P=0.

This can occur because of symmetry or block structure, but it should be checked rather than assumed.

Three unperturbed levels have energies 00, δ\delta, and Δ\Delta, with couplings of characteristic size vv. Assume

δ∼v,v≪Δ.\delta\sim v, \qquad v\ll\Delta.

Which levels belong in PP, and what is the natural perturbative parameter?

Solution

The levels at 00 and δ\delta are separated by an amount comparable with their coupling. They mix strongly and must both be retained:

P=∣0⟩⟨0∣+∣δ⟩⟨δ∣.P = \lvert0\rangle\langle0\rvert + \lvert\delta\rangle\langle\delta\rvert.

The level at Δ\Delta is remote and belongs in QQ. Coupling to it can be expanded in

η∼vΔ.\eta \sim \frac{v}{\Delta}.

Inside PP, the detuning δ\delta and coupling vv should be diagonalized together. Expanding in v/δv/\delta would be uncontrolled because that ratio is of order unity.

  1. H. Feshbach, “Unified theory of nuclear reactions,” Annals of Physics 5, 357–390 (1958).
  2. H. Feshbach, “A unified theory of nuclear reactions. II,” Annals of Physics 19, 287–313 (1962).
  3. J. R. Schrieffer and P. A. Wolff, “Relation between the Anderson and Kondo Hamiltonians,” Physical Review 149, 491–492 (1966).
  4. S. Bravyi, D. P. DiVincenzo, and D. Loss, “Schrieffer–Wolff transformation for quantum many-body systems,” Annals of Physics 326, 2793–2826 (2011).
  5. C. Bloch, “Sur la théorie des perturbations des états liés,” Nuclear Physics 6, 329–347 (1958).
  6. C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, Atom–Photon Interactions: Basic Processes and Applications, Wiley (1992).
  7. F. Reiter and A. S. Sørensen, “Effective operator formalism for open quantum systems,” Physical Review A 85, 032111 (2012).
  8. S. Blanes, F. Casas, J. A. Oteo, and J. Ros, “The Magnus expansion and some of its applications,” Physics Reports 470, 151–238 (2009).
  9. N. Goldman and J. Dalibard, “Periodically driven quantum systems: Effective Hamiltonians and engineered gauge fields,” Physical Review X 4, 031027 (2014); see also the published erratum.
  10. A. Eckardt and E. Anisimovas, “High-frequency approximation for periodically driven quantum systems from a Floquet-space perspective,” New Journal of Physics 17, 093039 (2015).
  11. S. Teufel, Adiabatic Perturbation Theory in Quantum Dynamics, Springer (2003).
  12. P. W. Anderson, “Antiferromagnetism. Theory of superexchange interaction,” Physical Review 79, 350–356 (1950).