Skip to content

Born Approximation

The Born approximation is the weak-coupling step in many open-system derivations. It says that, to the perturbative order being kept, the environment can be treated as remaining near a fixed reference state and system-environment correlations do not need to be promoted to independent dynamical variables.

In a common interaction-picture form,

ρSEI(t)≈ρSI(t)⊗ρE,\rho_{SE}^{I}(t) \approx \rho_S^{I}(t)\otimes\rho_E,

inside the second-order memory kernel. Here SS denotes the system and EE the environment. The approximation is not the Born rule for measurement probabilities. It is also not the claim that the exact state never becomes entangled with the environment.

The point is more delicate: weak coupling generates correlations, but their feedback on the reduced state is approximated by bath correlation functions computed in a stationary reference state.

For what is meant by “bath” and “reference state” in this context, see Baths, Reservoirs, and Environments.

Start from a microscopic Hamiltonian, as in System–Bath Hamiltonians:

H=HS+HE+λHI,HI=∑αSα⊗Bα.H = H_S+H_E+\lambda H_I, \qquad H_I = \sum_\alpha S_\alpha\otimes B_\alpha.

The parameter λ\lambda may be a physical coupling constant or a bookkeeping parameter for perturbation theory. The bath reference state ρE\rho_E is usually assumed stationary:

[HE,ρE]=0.[H_E,\rho_E]=0.

It is also convenient to center the bath operators:

Tr⁡E(BαρE)=0.\operatorname{Tr}_E(B_\alpha\rho_E)=0.

If this condition does not hold, write

Bα=δBα+Tr⁡E(BαρE)IEB_\alpha = \delta B_\alpha + \operatorname{Tr}_E(B_\alpha\rho_E)I_E

and absorb the mean term into the system Hamiltonian. Otherwise a first-order coherent shift can be mistaken for a dissipative effect.

In the interaction picture with respect to HS+HEH_S+H_E,

HI(t)=∑αSα(t)⊗Bα(t),H_I(t) = \sum_\alpha S_\alpha(t)\otimes B_\alpha(t),

where

Sα(t)=eiHSt/ℏSαe−iHSt/ℏ,Bα(t)=eiHEt/ℏBαe−iHEt/ℏ.S_\alpha(t) = e^{iH_St/\hbar} S_\alpha e^{-iH_St/\hbar}, \qquad B_\alpha(t) = e^{iH_Et/\hbar} B_\alpha e^{-iH_Et/\hbar}.

The exact total state obeys

ddtρSEI(t)=−iλℏ[HI(t),ρSEI(t)].\frac{d}{dt}\rho_{SE}^{I}(t) = -\frac{i\lambda}{\hbar} \left[ H_I(t), \rho_{SE}^{I}(t) \right].

The Born approximation enters only after this exact equation has been expanded perturbatively.

Second-Order Equation Before the Approximation

Section titled “Second-Order Equation Before the Approximation”

Integrating the interaction-picture equation once gives

ρSEI(t)=ρSEI(0)−iλℏ∫0tds [HI(s),ρSEI(s)].\rho_{SE}^{I}(t) = \rho_{SE}^{I}(0) - \frac{i\lambda}{\hbar} \int_0^t ds\, \left[ H_I(s), \rho_{SE}^{I}(s) \right].

Substitute this expression back into the equation for ρ˙SI(t)\dot\rho_S^{I}(t) and trace over the environment. The result is exact at this stage:

ddtρSI(t)=−iλℏTr⁡E[HI(t),ρSEI(0)]−λ2ℏ2∫0tds Tr⁡E[HI(t),[HI(s),ρSEI(s)]].\begin{aligned} \frac{d}{dt}\rho_S^{I}(t) =& -\frac{i\lambda}{\hbar} \operatorname{Tr}_E \left[ H_I(t), \rho_{SE}^{I}(0) \right] \\ & -\frac{\lambda^2}{\hbar^2} \int_0^t ds\, \operatorname{Tr}_E \left[ H_I(t), \left[ H_I(s), \rho_{SE}^{I}(s) \right] \right]. \end{aligned}

For an initially factorized state,

ρSEI(0)=ρS(0)⊗ρE,\rho_{SE}^{I}(0) = \rho_S(0)\otimes\rho_E,

with centered bath operators, the first-order term vanishes. Without centering, it produces a Hamiltonian correction.

The second-order term is still exact in form because it contains the true ρSEI(s)\rho_{SE}^{I}(s). It is not closed on ρSI\rho_S^{I}.

The Born approximation replaces the total state inside the second-order term by

ρSEI(s)≈ρSI(s)⊗ρE.\rho_{SE}^{I}(s) \approx \rho_S^{I}(s)\otimes\rho_E.

This gives the second-order Born equation

ddtρSI(t)=−λ2ℏ2∫0tds Tr⁡E[HI(t),[HI(s),ρSI(s)⊗ρE]],\frac{d}{dt}\rho_S^{I}(t) = -\frac{\lambda^2}{\hbar^2} \int_0^t ds\, \operatorname{Tr}_E \left[ H_I(t), \left[ H_I(s), \rho_S^{I}(s)\otimes\rho_E \right] \right],

after the centered first-order term has been removed.

This equation is not yet Markovian. The earlier state ρSI(s)\rho_S^{I}(s) still appears under the integral. The Born approximation is the weak-coupling factorization step; the Markov Approximation is the later memory-shortening step.

For a stationary bath, define

Cαβ(t−s)=Tr⁡E[Bα(t)Bβ(s)ρE].C_{\alpha\beta}(t-s) = \operatorname{Tr}_E \left[ B_\alpha(t)B_\beta(s)\rho_E \right].

Then the Born equation can be written as

ddtρSI(t)=−λ2ℏ2∑α,β∫0tds (Cαβ(t−s)[Sα(t),Sβ(s)ρSI(s)]+Cβα(s−t)[ρSI(s)Sβ(s),Sα(t)]).\begin{aligned} \frac{d}{dt}\rho_S^{I}(t) = -\frac{\lambda^2}{\hbar^2} \sum_{\alpha,\beta} \int_0^t ds\, \big( & C_{\alpha\beta}(t-s) \left[ S_\alpha(t), S_\beta(s)\rho_S^{I}(s) \right] \\ & + C_{\beta\alpha}(s-t) \left[ \rho_S^{I}(s)S_\beta(s), S_\alpha(t) \right] \big). \end{aligned}

This expression shows what the approximation keeps. The bath affects the system through two-point correlation functions computed in ρE\rho_E. It does not require the bath to be featureless. Long-lived, oscillatory, or structured correlations can still appear here; they become a problem only if later approximations discard memory that is physically important.

The Born approximation neglects dynamical variables for the correlated part of the total state. Write schematically

ρSE(t)=ρS(t)⊗ρE+χSE(t),\rho_{SE}(t) = \rho_S(t)\otimes\rho_E + \chi_{SE}(t),

where

Tr⁡EχSE(t)=0.\operatorname{Tr}_E\chi_{SE}(t)=0.

The approximation does not assert χSE(t)=0\chi_{SE}(t)=0 exactly. Rather, χSE\chi_{SE} is generated perturbatively and its leading influence is represented by bath correlation functions. Corrections involving the bath’s changed state, higher cumulants, repeated backaction, and nonperturbative dressing are neglected.

Typical neglected effects include:

  • bath depletion or heating;
  • coherent reabsorption from a narrow mode;
  • persistent system–bath entanglement;
  • initial correlations from prior equilibration;
  • strong-coupling renormalization of the system Hamiltonian;
  • memory variables that cannot be summarized by a fixed ρE\rho_E.

The approximation is therefore a statement about a controlled reduced description, not about the absence of correlations in the exact total state.

There is no universal single inequality called “the Born condition.” In practice, one checks a family of small parameters and time scales.

A useful schematic requirement is that the rate induced by the bath be slow compared with microscopic bath dynamics. If τE\tau_E is a bath correlation time and Γ\Gamma is a relaxation or dephasing rate generated at order λ2\lambda^2, one wants

ΓτE≪1.\Gamma\tau_E\ll1.

This condition says that the system changes little over one bath memory time. It does not by itself justify every later step, but it is often part of the same hierarchy.

Other checks are equally important:

  • the dimensionless coupling is small compared with the relevant system and bath scales;
  • the bath is large or mixing enough that its reference state is not appreciably depleted;
  • bath correlations decay or average in a way compatible with the prediction being made;
  • the time interval of interest does not accumulate strong backaction on the environment;
  • the preparation is compatible with a factorized or weakly correlated initial state.

The Approximation Checklist gives these checks in a practical audit format.

Weak-coupling master-equation derivations often combine three steps:

StepMain assumptionTypical mathematical move
Bornweak system–bath correlationsreplace ρSEI(s)\rho_{SE}^{I}(s) by ρSI(s)⊗ρE\rho_S^{I}(s)\otimes\rho_E in the second-order kernel
Markovshort bath memory relative to resolved system evolutionreplace ρSI(s)\rho_S^{I}(s) by ρSI(t)\rho_S^{I}(t) and extend a memory integral
Secularwell-separated Bohr-frequency sectorsaverage rapidly rotating cross terms

These steps are logically separate. A Born equation may remain non-Markovian. A Markovian Redfield Equation may fail to be completely positive. A secular GKSL equation may be too coarse near nearly degenerate transitions.

The Lindblad–GKSL Equation is the canonical form reached after additional positivity-preserving assumptions. Complete positivity is discussed separately in Completely Positive Maps.

An atom weakly coupled to a broad continuum of electromagnetic modes is often a good Born-approximation setting. The field is large, the vacuum or thermal reference state is not significantly depleted by one atom, and the induced decay rate is small compared with optical frequencies. Further Markov and secular approximations are still separate claims.

For a qubit with

HI=σx⊗B,H_I = \sigma_x\otimes B,

the coupling contains raising and lowering components relative to HSH_S. In weak coupling, bath spectra near the qubit frequency determine relaxation rates, while low-frequency components can contribute to dephasing. The resulting channel may resemble amplitude damping or dephasing in suitable limits.

A two-level atom strongly or resonantly coupled to a single cavity mode is a poor bath model unless the cavity itself is damped by a larger reservoir. Energy can flow into the mode and back to the atom. Treating the mode as a fixed environment can miss Rabi oscillations and memory.

A finite spin bath can store information about the system for long times. The bath state may not relax to a reference state, and repeated backaction can matter. A Born approximation may still be useful for short times or weak inhomogeneous coupling, but it is not automatic.

Confusing the Born approximation with the Born rule

Section titled “Confusing the Born approximation with the Born rule”

The Born rule assigns probabilities such as

p(a)=Tr⁡(Paρ).p(a)=\operatorname{Tr}(P_a\rho).

The Born approximation in open systems is a perturbative factorization assumption about a system and an environment. The shared name is historical, not conceptual.

Saying the system never entangles with the bath

Section titled “Saying the system never entangles with the bath”

Weak coupling still creates entanglement and classical correlations. The approximation says those correlations do not need independent dynamical variables at the order being retained.

If Tr⁡E(BαρE)≠0\operatorname{Tr}_E(B_\alpha\rho_E)\ne0, there is a coherent mean force. Centering the bath operator or absorbing the mean into HSH_S should happen before identifying dissipative rates.

Born closure alone does not guarantee a GKSL generator. Positivity and complete positivity usually require additional Markov, secular, coarse-graining, or exact channel arguments.

If the system and bath were equilibrated together before the experiment, the initial state is generally correlated. Applying a factorized weak-coupling derivation to such a state can produce wrong short-time behavior or an incorrect assignment map.

Let

HI=S⊗BH_I = S\otimes B

and let the bath reference state satisfy ⟨B⟩E=Tr⁡E(BρE)\langle B\rangle_E=\operatorname{Tr}_E(B\rho_E). Show that the interaction can be rewritten as a centered interaction plus a system Hamiltonian correction.

Solution

Write

B=δB+⟨B⟩EIE,δB=B−⟨B⟩EIE.B = \delta B+\langle B\rangle_E I_E, \qquad \delta B = B-\langle B\rangle_E I_E.

Then

S⊗B=S⊗δB+⟨B⟩ES⊗IE.S\otimes B = S\otimes\delta B + \langle B\rangle_E S\otimes I_E.

The second term acts only on the system and can be absorbed into

HS↦HS+λ⟨B⟩ES.H_S \mapsto H_S+\lambda\langle B\rangle_E S.

The remaining bath operator is centered because

Tr⁡E(δBρE)=0.\operatorname{Tr}_E(\delta B\rho_E)=0.

Starting from

ddtρSEI(t)=−iλℏ[HI(t),ρSEI(t)],\frac{d}{dt}\rho_{SE}^{I}(t) = -\frac{i\lambda}{\hbar} \left[ H_I(t), \rho_{SE}^{I}(t) \right],

integrate once and identify the precise step at which the Born approximation is made.

Solution

Integration gives

ρSEI(t)=ρSEI(0)−iλℏ∫0tds [HI(s),ρSEI(s)].\rho_{SE}^{I}(t) = \rho_{SE}^{I}(0) - \frac{i\lambda}{\hbar} \int_0^t ds\, \left[ H_I(s), \rho_{SE}^{I}(s) \right].

Substitution into the equation for ρ˙SI\dot\rho_S^{I} gives a second-order expression containing

Tr⁡E[HI(t),[HI(s),ρSEI(s)]].\operatorname{Tr}_E \left[ H_I(t), \left[ H_I(s), \rho_{SE}^{I}(s) \right] \right].

The Born approximation is the replacement

ρSEI(s)≈ρSI(s)⊗ρE\rho_{SE}^{I}(s) \approx \rho_S^{I}(s)\otimes\rho_E

inside this term. Keeping ρSI(s)\rho_S^{I}(s) rather than replacing it by ρSI(t)\rho_S^{I}(t) means the equation is not yet Markovian.

Classify each statement as a Born assumption, a Markov assumption, both, or neither.

  1. The bath remains approximately in a stationary reference state.
  2. The memory integral can be extended from tt to ∞\infty.
  3. The system and bath never become entangled in the exact dynamics.
  4. The bath correlation time is much shorter than the relaxation time.
Solution
  1. Born.
  2. Markov, provided the bath correlations decay sufficiently fast.
  3. Neither. It is stronger than the Born approximation and generally false.
  4. It supports the Markov approximation and is also a useful consistency check for weak coupling, but by itself it is not the Born factorization step.

A qubit is coupled resonantly to one high-quality oscillator mode. The oscillator starts in its ground state. Is the Born approximation reliable for the oscillator treated as the environment?

Solution

Not generically. A single high-quality mode can store excitation and return it to the qubit. Its state can change appreciably, and correlations with the qubit can persist. A better split may include the oscillator in the system and treat the oscillator’s damping reservoir as the bath.

  • H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press, 2002.
  • Á. Rivas and S. F. Huelga, Open Quantum Systems: An Introduction, Springer, 2012.
  • C. W. Gardiner and P. Zoller, Quantum Noise, 3rd ed., Springer, 2004.
  • U. Weiss, Quantum Dissipative Systems, 4th ed., World Scientific, 2012.
  • R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, 2nd ed., Springer, 2007.
  • E. B. Davies, “Markovian master equations,” Communications in Mathematical Physics 39, 91–110 (1974).