Skip to content

Baker–Campbell–Hausdorff Formula

The Baker–Campbell–Hausdorff formula, often abbreviated BCH, describes how exponentials of noncommuting operators combine. In quantum dynamics it answers two recurring questions:

  1. What is the effective exponent of a product eAeBe^Ae^B?
  2. How does an operator change under conjugation eABe−Ae^A B e^{-A}?

Both answers are controlled by nested commutators. That is why BCH appears in Heisenberg evolution, interaction-picture transformations, coherent-state displacement operators, spin rotations, product-formula errors, and Lie-algebraic approximation schemes.

For the compact formula-card version, see Baker–Campbell–Hausdorff. For the general finite-dimensional construction of operator exponentials, see Matrix Functions and Exponentials.

If AA and BB commute, then

eA+B=eAeB.e^{A+B} = e^Ae^B.

This identity is safe because the power series can be rearranged like scalar products. If [A,B]≠0[A,B]\ne0, the rearrangement fails. The product eAeBe^Ae^B is still meaningful, but its logarithm is not usually A+BA+B.

The first correction is already enough to explain many quantum effects:

log⁡(eAeB)=A+B+12[A,B]+⋯ .\log(e^Ae^B) = A+B+\frac12[A,B]+\cdots .

The ellipsis matters. Unless nested commutators vanish or are controlled, replacing the product by only the first few terms is an approximation.

The BCH series begins

log⁡(eAeB)=A+B+12[A,B]+112[A,[A,B]]+112[B,[B,A]]+⋯ .\begin{aligned} \log(e^Ae^B) &= A+B+\frac12[A,B] \\ &\quad +\frac1{12}[A,[A,B]] +\frac1{12}[B,[B,A]] +\cdots . \end{aligned}

Equivalently,

eAeB=exp⁡(A+B+12[A,B]+⋯ ).e^Ae^B = \exp \left( A+B+\frac12[A,B]+\cdots \right).

This is an identity under appropriate convergence assumptions, and it is also used formally in Lie-algebraic calculations. In quantum mechanics, formal use is common, but it should not hide domain questions for unbounded operators.

A particularly useful special case occurs when

[A,B][A,B]

commutes with both AA and BB. Then all higher nested commutators vanish, and

eAeB=exp⁡(A+B+12[A,B]).e^Ae^B = \exp \left( A+B+\frac12[A,B] \right).

The same assumption gives the reordering identity

eAeB=eBeAe[A,B].e^Ae^B = e^Be^A e^{[A,B]}.

This central-commutator case is the algebra behind many Weyl, displacement, and canonical-commutator manipulations. It is also a common place to make mistakes: the truncation is not valid merely because [A,B][A,B] is simple-looking; it must commute with AA and BB.

The conjugation identity is often the more useful dynamics formula. Define

ad⁡A(B)=[A,B].\operatorname{ad}_A(B) = [A,B].

Then

eABe−A=ead⁡AB,e^A B e^{-A} = e^{\operatorname{ad}_A}B,

meaning

eABe−A=B+[A,B]+12![A,[A,B]]+13![A,[A,[A,B]]]+⋯ .e^A B e^{-A} = B+[A,B] +\frac{1}{2!}[A,[A,B]] +\frac{1}{3!}[A,[A,[A,B]]] +\cdots .

This identity is sometimes called the Hadamard lemma. It is the local engine behind operator picture transformations: conjugation by an exponential becomes a series of commutators.

For a time-independent Hamiltonian,

AH(t)=eiHt/ℏAe−iHt/ℏ.A_H(t) = e^{iHt/\hbar}A e^{-iHt/\hbar}.

Use the similarity formula with

Agen=iHtℏ.A_{\rm gen} = \frac{iHt}{\hbar}.

Then

AH(t)=A+itℏ[H,A]+12!(itℏ)2[H,[H,A]]+⋯ .\begin{aligned} A_H(t) &= A + \frac{it}{\hbar}[H,A] \\ &\quad + \frac{1}{2!} \left( \frac{it}{\hbar} \right)^2 [H,[H,A]] +\cdots . \end{aligned}

Differentiating this expansion at t=0t=0 gives the Heisenberg equation

A˙H=iℏ[H,AH].\dot A_H = \frac{i}{\hbar}[H,A_H].

The nested commutators describe the finite-time orbit of AA under Hamiltonian conjugation. If they close on a finite-dimensional operator subspace, the dynamics can often be solved by finite-dimensional linear algebra.

BCH also explains why splitting a Hamiltonian into noncommuting pieces produces finite-step errors. Let

X=−iℏHAΔt,Y=−iℏHBΔt.X = -\frac{i}{\hbar}H_A\Delta t, \qquad Y = -\frac{i}{\hbar}H_B\Delta t.

Then

eXeY=exp⁡[X+Y+12[X,Y]+O(Δt3)].e^Xe^Y = \exp \left[ X+Y+\frac12[X,Y]+O(\Delta t^3) \right].

Since

[X,Y]=−Δt2ℏ2[HA,HB],[X,Y] = -\frac{\Delta t^2}{\hbar^2}[H_A,H_B],

one step contains the correction

−Δt22ℏ2[HA,HB].-\frac{\Delta t^2}{2\hbar^2}[H_A,H_B].

This is the local source of the first-order Trotter error. The convergence and higher-order splitting discussion belongs to Trotter Product Formula.

For one oscillator mode, define

D(α)=exp⁡(αa†−α∗a).D(\alpha) = \exp \left( \alpha a^\dagger-\alpha^*a \right).

Let

K=αa†−α∗a.K = \alpha a^\dagger-\alpha^*a.

Using [a,a†]=1[a,a^\dagger]=1,

[K,a]=−α,[K,[K,a]]=0.[K,a] = -\alpha, \qquad [K,[K,a]]=0.

The similarity formula truncates:

D(α)aD†(α)=a−α,D(\alpha)aD^\dagger(\alpha) = a-\alpha,

and equivalently

D†(α)aD(α)=a+α.D^\dagger(\alpha)aD(\alpha) = a+\alpha.

This is why displacement operators create coherent states from the oscillator ground state. The coherent-state page uses this result as a construction; the operator reference card records common convention warnings in Displacement Operator.

Angular-momentum rotations are another closed commutator system. With

Uz(θ)=exp⁡(−iℏθSz),U_z(\theta) = \exp \left( -\frac{i}{\hbar}\theta S_z \right),

and

[Sz,Sx]=iℏSy,[Sz,Sy]=−iℏSx,[S_z,S_x]=i\hbar S_y, \qquad [S_z,S_y]=-i\hbar S_x,

the nested commutators alternate between SxS_x and SyS_y. One convention gives

Uz(θ)SxUz†(θ)=Sxcos⁡θ+Sysin⁡θ.U_z(\theta)S_xU_z^\dagger(\theta) = S_x\cos\theta+S_y\sin\theta.

The series has summed to an ordinary rotation. The spinor-specific half-angle formula is developed in Spin Rotations; BCH explains the operator-algebra mechanism behind the rotation of spin components.

BCH combines a finite product of exponentials into one effective exponential. Time ordering solves a different problem: it organizes a continuous product when the Hamiltonian at different times does not commute with itself.

For a time-dependent Hamiltonian, one cannot generally write

U(t,t0)=exp⁡[−iℏ∫t0tH(s) ds].U(t,t_0) = \exp \left[ -\frac{i}{\hbar} \int_{t_0}^t H(s)\,ds \right].

The correct object is a time-ordered exponential unless

[H(t1),H(t2)]=0[H(t_1),H(t_2)]=0

for the times being compared, or unless a special transformation reduces the problem. For that organization, see Time Ordering.

  • Using eA+B=eAeBe^{A+B}=e^Ae^B for noncommuting operators.
  • Applying the central-commutator truncation without checking nested commutators.
  • Treating BCH as a convergence theorem for unbounded operators without domain assumptions.
  • Confusing BCH with the Trotter product formula or a time-ordered exponential.
  • Dropping the phase generated when displacement operators are reordered.
  • Forgetting that active and passive rotation conventions change signs in conjugation formulas.
  • Keeping only the first commutator term while needing a finite-time transformation.
  • B. C. Hall, Lie Groups, Lie Algebras, and Representations, 2nd ed., Springer, 2015.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  • N. J. Higham, Functions of Matrices: Theory and Computation, SIAM, 2008.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • M. O. Scully and M. S. Zubairy, Quantum Optics, Cambridge University Press, 1997.
  1. Suppose [A,B][A,B] commutes with both AA and BB. What is eAeBe^Ae^B as a single exponential?
Solution

Under this central-commutator assumption, all nested commutators beyond [A,B][A,B] vanish in the BCH series. Therefore

eAeB=exp⁡(A+B+12[A,B]).e^Ae^B = \exp \left( A+B+\frac12[A,B] \right).
  1. Use the similarity formula to show that if [A,B]=cI[A,B]=cI, then eABe−A=B+cIe^A B e^{-A}=B+cI.
Solution

The similarity expansion is

eABe−A=B+[A,B]+12![A,[A,B]]+⋯ .e^A B e^{-A} = B+[A,B] +\frac{1}{2!}[A,[A,B]] +\cdots .

If [A,B]=cI[A,B]=cI, then [A,[A,B]]=[A,cI]=0[A,[A,B]]=[A,cI]=0. Thus the series truncates:

eABe−A=B+cI.e^A B e^{-A} = B+cI.
  1. For D(α)=exp⁡(αa†−α∗a)D(\alpha)=\exp(\alpha a^\dagger-\alpha^*a), derive D†(α)aD(α)=a+αD^\dagger(\alpha)aD(\alpha)=a+\alpha.
Solution

Let K=αa†−α∗aK=\alpha a^\dagger-\alpha^*a, so D=eKD=e^K and D†=e−KD^\dagger=e^{-K}. Since

[K,a]=α[a†,a]−α∗[a,a]=−α,[K,a] = \alpha[a^\dagger,a]-\alpha^*[a,a] = -\alpha,

one has

[−K,a]=α.[-K,a] = \alpha.

All higher nested commutators vanish because α\alpha is a scalar. Therefore

D†aD=e−KaeK=a+α.D^\dagger aD = e^{-K}ae^K = a+\alpha.
  1. Identify the leading BCH correction in one first-order split step with X=−iHAΔt/ℏX=-iH_A\Delta t/\hbar and Y=−iHBΔt/ℏY=-iH_B\Delta t/\hbar.
Solution

BCH gives

eXeY=exp⁡(X+Y+12[X,Y]+O(Δt3)).e^Xe^Y = \exp \left( X+Y+\frac12[X,Y]+O(\Delta t^3) \right).

Now

[X,Y]=[−iℏHAΔt,−iℏHBΔt]=−Δt2ℏ2[HA,HB].[X,Y] = \left[ -\frac{i}{\hbar}H_A\Delta t, -\frac{i}{\hbar}H_B\Delta t \right] = -\frac{\Delta t^2}{\hbar^2}[H_A,H_B].

The leading correction in the exponent is therefore

−Δt22ℏ2[HA,HB].-\frac{\Delta t^2}{2\hbar^2}[H_A,H_B].