Skip to content

Baker–Campbell–Hausdorff

The Baker–Campbell–Hausdorff formula, abbreviated BCH, replaces a product of noncommuting exponentials by one exponential:

eAeB=eZ,Z=log⁡(eAeB).e^A e^B=e^Z, \qquad Z=\log(e^Ae^B).

The exponent ZZ is built from AA, BB, and their nested commutators. BCH is therefore the right tool when the order of two finite transformations matters. The closely related Hadamard lemma evaluates a conjugation eABe−Ae^A B e^{-A}.

This page is a formula and diagnostic reference. The extended physical development is at Baker–Campbell–Hausdorff Formula.

Define the adjoint action

ad⁡A(B)≡[A,B],ad⁡A 0(B)=B.\operatorname{ad}_A(B)\equiv[A,B], \qquad \operatorname{ad}_A^{\,0}(B)=B.

The principal formulas are:

CaseResult
General BCH serieslog⁡(eAeB)=A+B+12[A,B]+⋯\log(e^Ae^B)=A+B+\frac12[A,B]+\cdots
Commuting operators[A,B]=0 ⟹ eAeB=eA+B[A,B]=0\ \Longrightarrow\ e^Ae^B=e^{A+B}
Central commutator[A,[A,B]]=[B,[A,B]]=0 ⟹ eAeB=eA+B+12[A,B][A,[A,B]]=[B,[A,B]]=0\ \Longrightarrow\ e^Ae^B=e^{A+B+\frac12[A,B]}
Central reorderingeAeB=e[A,B]eBeAe^Ae^B=e^{[A,B]}e^Be^A
Hadamard lemmaeABe−A=ead⁡ABe^A B e^{-A}=e^{\operatorname{ad}_A}B
Leading splitting correctionetAetB=exp⁡ ⁣(t(A+B)+t22[A,B]+O(t3))e^{tA}e^{tB}=\exp\!\left(t(A+B)+\frac{t^2}{2}[A,B]+O(t^3)\right)

The central-commutator rows require [A,B][A,B] to commute with both AA and BB. They do not follow merely from [A,B][A,B] being easy to calculate.

Through commutator degree four,

Z=A+B+12[A,B]+112[A,[A,B]]+112[B,[B,A]]−124[B,[A,[A,B]]]+⋯ .\begin{aligned} Z &= A+B+\frac12[A,B] \\ &\quad +\frac1{12}[A,[A,B]] +\frac1{12}[B,[B,A]] \\ &\quad -\frac1{24}[B,[A,[A,B]]] +\cdots . \end{aligned}

Here degree counts occurrences of AA and BB, not the number of bracket symbols. Thus [A,B][A,B] has degree two and [A,[A,B]][A,[A,B]] has degree three. Equivalent-looking expressions can be related by antisymmetry and the Jacobi identity. Fix one ordering convention before comparing signs.

A reliable way to organize approximations is to introduce a bookkeeping parameter:

log⁡(etAetB)=t(A+B)+t22[A,B]+t312([A,[A,B]]+[B,[B,A]])−t424[B,[A,[A,B]]]+O(t5).\begin{aligned} \log(e^{tA}e^{tB}) &= t(A+B) +\frac{t^2}{2}[A,B] \\ &\quad +\frac{t^3}{12} \left( [A,[A,B]] +[B,[B,A]] \right) \\ &\quad -\frac{t^4}{24} [B,[A,[A,B]]] +O(t^5). \end{aligned}

This form makes the order of a discarded term explicit. An expansion in tt is controlled only when the operator scales and the relevant time or step size make the neglected terms small.

For scalars, the power series for eAeBe^Ae^B can be freely reordered and collected into eA+Be^{A+B}. Operator products cannot generally be reordered. The first mismatch is

eAeB=I+(A+B)+12(A2+B2)+AB+⋯ ,\begin{aligned} e^Ae^B &= I+(A+B) \\ &\quad +\frac12(A^2+B^2)+AB+\cdots , \end{aligned}

whereas

eA+B=I+(A+B)+12(A2+AB+BA+B2)+⋯ .\begin{aligned} e^{A+B} &= I+(A+B) \\ &\quad +\frac12(A^2+AB+BA+B^2)+\cdots . \end{aligned}

Their quadratic difference is

12(AB−BA)=12[A,B].\frac12(AB-BA)=\frac12[A,B].

BCH recursively packages every higher ordering mismatch into nested commutators. If AA and BB belong to a Lie algebra closed under commutators, then ZZ remains in that Lie algebra wherever the logarithm and series are well defined.

If

[A,B]=0,[A,B]=0,

all nested commutators vanish and

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

Pairwise commutation is essential when extending this statement to more than two exponentials.

Suppose

[A,[A,B]]=0,[B,[A,B]]=0.[A,[A,B]]=0, \qquad [B,[A,B]]=0.

Equivalently, [A,B][A,B] is central within the algebra generated by AA and BB. Then the BCH series stops after its first commutator:

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

Reversing AA and BB changes the commutator sign:

eBeA=exp⁡ ⁣(A+B−12[A,B]).e^Be^A = \exp\!\left( A+B-\frac12[A,B] \right).

Because the commutator is central,

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

This identity is the source of the phase or scalar factor that appears when Weyl operators and oscillator displacement factors are reordered.

More generally, BCH terminates if sufficiently long nested commutators vanish. For example, if every commutator of degree three vanishes, only A+B+12[A,B]A+B+\frac12[A,B] remains. Closure on a finite-dimensional operator space does not by itself imply termination: the infinite series may instead sum to a finite combination of basis operators, as it does for rotations.

Conjugation by an exponential is generated by repeated commutation:

eABe−A=ead⁡AB=B+[A,B]+12![A,[A,B]]+13![A,[A,[A,B]]]+⋯ .\begin{aligned} e^A B e^{-A} &= e^{\operatorname{ad}_A}B \\ &= B+[A,B] +\frac1{2!}[A,[A,B]] \\ &\quad +\frac1{3!}[A,[A,[A,B]]] +\cdots . \end{aligned}

The opposite ordering changes the sign:

e−ABeA=e−ad⁡AB=B−[A,B]+12![A,[A,B]]−⋯ .e^{-A}Be^A = e^{-\operatorname{ad}_A}B = B-[A,B] +\frac1{2!}[A,[A,B]] -\cdots .

To derive the series, define

F(s)=esABe−sA.F(s)=e^{sA}Be^{-sA}.

Differentiation gives

dFds=[A,F(s)]=ad⁡AF(s),F(0)=B.\frac{dF}{ds} = [A,F(s)] = \operatorname{ad}_A F(s), \qquad F(0)=B.

The formal solution is F(s)=esad⁡ABF(s)=e^{s\operatorname{ad}_A}B; setting s=1s=1 produces the formula. For unbounded operators, this differentiation assumes a domain on which the products and derivatives are legitimate.

Calculate nested commutators in order:

C1=[A,B],C2=[A,C1],C3=[A,C2],…C_1=[A,B], \qquad C_2=[A,C_1], \qquad C_3=[A,C_2], \quad\ldots

If Cn+1=0C_{n+1}=0, Hadamard’s series stops at Cn/n!C_n/n!. If instead the CnC_n close on a small operator basis, convert the commutator action into a matrix on that basis and exponentiate that smaller matrix.

Let [a,a†]=I[a,a^\dagger]=I and define

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

Take

A=αa†,B=−α∗a.A=\alpha a^\dagger, \qquad B=-\alpha^*a.

Their commutator is central:

[A,B]=∣α∣2I.[A,B] = |\alpha|^2 I.

The central BCH formula gives

eαa†e−α∗a=e∣α∣2/2D(α),e^{\alpha a^\dagger}e^{-\alpha^*a} = e^{|\alpha|^2/2}D(\alpha),

or the normally ordered factorization

D(α)=e−∣α∣2/2eαa†e−α∗a.D(\alpha) = e^{-|\alpha|^2/2} e^{\alpha a^\dagger} e^{-\alpha^*a}.

Reversing the factors yields

D(α)=e∣α∣2/2e−α∗aeαa†.D(\alpha) = e^{|\alpha|^2/2} e^{-\alpha^*a} e^{\alpha a^\dagger}.

The scalar prefactor is not optional. Dropping it changes normalization and composition phases.

Hadamard’s lemma also gives

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

because the first commutator is αI\alpha I and all further nested commutators vanish.

For the canonical pair [X,P]=iℏI[X,P]=i\hbar I, the translation operator

T(a)=e−iaP/ℏT(a)=e^{-iaP/\hbar}

acts by conjugation as

T†(a)XT(a)=X+aI.T^\dagger(a)XT(a)=X+aI.

Indeed, with A=iaP/ℏA=iaP/\hbar,

[A,X]=aI,[A,X]=aI,

and the Hadamard series truncates immediately. The result is a useful sign check: T(a)T(a) translates states to the right, while T†XTT^\dagger X T shifts the position observable by +a+a in the transformed state.

For

U(a)=e−iaP/ℏ,V(b)=eibX/ℏ,U(a)=e^{-iaP/\hbar}, \qquad V(b)=e^{ibX/\hbar},

the exponent commutator is central, and reordering gives the Weyl relation

U(a)V(b)=e−iab/ℏV(b)U(a).U(a)V(b) = e^{-iab/\hbar}V(b)U(a).

This exponentiated relation is often better behaved than manipulating the unbounded operators XX and PP directly.

For a small parameter tt,

etAetB=exp⁡ ⁣(t(A+B)+t22[A,B]+O(t3)).e^{tA}e^{tB} = \exp\!\left( t(A+B)+\frac{t^2}{2}[A,B]+O(t^3) \right).

Thus replacing et(A+B)e^{t(A+B)} by etAetBe^{tA}e^{tB} has a local discrepancy whose leading exponent is quadratic in tt and proportional to [A,B][A,B]. For a Hamiltonian split H=HA+HBH=H_A+H_B, set

A=−iℏHA,B=−iℏHB.A=-\frac{i}{\hbar}H_A, \qquad B=-\frac{i}{\hbar}H_B.

Then

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

The symmetric product

etA/2etBetA/2e^{tA/2}e^{tB}e^{tA/2}

cancels the degree-two commutator correction; its exponent differs from t(A+B)t(A+B) first at order t3t^3. BCH diagnoses these local error structures. Norm convergence over many steps belongs to the Trotter Product Formula.

Let

A=iασx,B=iβσy,A=i\alpha\sigma_x, \qquad B=i\beta\sigma_y,

with real α\alpha and β\beta. Since

[σx,σy]=2iσz,[\sigma_x,\sigma_y]=2i\sigma_z,

the first commutator is

[A,B]=−2iαβσz.[A,B] = -2i\alpha\beta\sigma_z.

Therefore, for small angles,

log⁡(eiασxeiβσy)=iασx+iβσy−iαβσz+O(α2β,αβ2).\log(e^{i\alpha\sigma_x}e^{i\beta\sigma_y}) = i\alpha\sigma_x +i\beta\sigma_y -i\alpha\beta\sigma_z +O(\alpha^2\beta,\alpha\beta^2).

The generated σz\sigma_z term is a direct signature of noncommuting rotations. Its sign reverses when the two exponentials are reversed.

BCH has three common mathematical interpretations:

  1. Formal series. Treat AA and BB as noncommuting symbols and organize terms by degree. This is enough for perturbative algebra but makes no analytic convergence claim.
  2. Local matrix or bounded-operator identity. Near the identity, eAeBe^Ae^B has a consistent logarithm branch and the convergent BCH series represents it.
  3. Lie-theoretic coordinate formula. Locally, multiplying two Lie-group elements corresponds to a nonlinear composition law in the Lie algebra.

Even for finite matrices, a global logarithm is branch dependent, so the symbol log⁡(eAeB)\log(e^Ae^B) must be interpreted locally or with a specified branch. For unbounded quantum operators, A+BA+B, [A,B][A,B], and their nested products may have different domains. The existence of the unitary exponentials does not by itself validate termwise BCH manipulation on the full Hilbert space.

When AA and BB are anti-Hermitian finite matrices, eAe^A and eBe^B are unitary. Every displayed BCH term is then anti-Hermitian, so any consistent finite-order truncation preserves the expected generator type, although the truncation need not reproduce the exact product.

  1. State the ordering to be combined: eAeBe^Ae^B is not interchangeable with eBeAe^Be^A.
  2. Compute [A,B][A,B] and then the nested commutators needed at the desired order.
  3. Test whether [A,B][A,B] is central or whether a nested-commutator chain terminates.
  4. Introduce a small parameter when using a truncation, and count total degree consistently.
  5. Check Hermiticity or anti-Hermiticity, dimensions, and the commuting limit.
  6. For unbounded operators, identify a common invariant domain or use a rigorous exponentiated relation.

For conjugation, use Hadamard’s lemma directly. Computing log⁡(eAeB)\log(e^Ae^B) first is usually unnecessary.

  • Writing eA+B=eAeBe^{A+B}=e^Ae^B without verifying [A,B]=0[A,B]=0.
  • Using the central formula after checking only that [A,B]≠0[A,B]\ne0, rather than checking both second nested commutators.
  • Forgetting that reversing AA and BB changes the sign of the first commutator correction.
  • Confusing eABe−Ae^ABe^{-A} with e−ABeAe^{-A}Be^A.
  • Mixing commutator degree with the number of nested brackets.
  • Dropping the scalar or phase factor when displacement or Weyl operators are reordered.
  • Treating an asymptotic or formal truncation as an exact identity.
  • Assuming a finite-dimensional closed commutator algebra makes the BCH series terminate.
  • Confusing BCH with time ordering. BCH combines discrete exponentials; time ordering organizes evolution generated at continuously varying times.
  • Ignoring logarithm branches for finite matrices or domains for unbounded operators.
  • B. C. Hall, Lie Groups, Lie Algebras, and Representations: An Elementary Introduction, 2nd ed., Springer, 2015, Chapters 2–3.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013, Sections 14.2–14.3.
  • N. J. Higham, Functions of Matrices: Theory and Computation, SIAM, 2008, Chapters 10–11.
  • M. W. Reinsch, “A simple expression for the terms in the Baker–Campbell–Hausdorff series,” Journal of Mathematical Physics 41, 2434–2442 (2000).
  • A. Van-Brunt and M. Visser, “Special-case closed form of the Baker–Campbell–Hausdorff formula,” Journal of Physics A: Mathematical and Theoretical 48, 225207 (2015).
  • W. Magnus, “On the exponential solution of differential equations for a linear operator,” Communications on Pure and Applied Mathematics 7, 649–673 (1954).
  • M. Suzuki, “Generalized Trotter’s formula and systematic approximants of exponential operators and inner derivations with applications to many-body problems,” Communications in Mathematical Physics 51, 183–190 (1976).
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  1. Suppose [A,B]=cI[A,B]=cI. Express eAeBe^Ae^B as a single exponential and then reorder it as a multiple of eBeAe^Be^A.
Solution

Because cIcI commutes with both AA and BB, all nested commutators vanish:

eAeB=exp⁡ ⁣(A+B+c2I).e^Ae^B = \exp\!\left( A+B+\frac{c}{2}I \right).

Reversing the factors gives

eBeA=exp⁡ ⁣(A+B−c2I).e^Be^A = \exp\!\left( A+B-\frac{c}{2}I \right).

The exponents differ by cIcI, which commutes with everything, so

eAeB=eceBeA.e^Ae^B=e^c e^Be^A.
  1. Derive the normally ordered factorization of D(α)D(\alpha) from [a,a†]=I[a,a^\dagger]=I.
Solution

Set A=αa†A=\alpha a^\dagger and B=−α∗aB=-\alpha^*a. Then

[A,B]=−∣α∣2[a†,a]=∣α∣2I.\begin{aligned} [A,B] &= -|\alpha|^2[a^\dagger,a] \\ &= |\alpha|^2 I. \end{aligned}

This commutator is central, so

eAeB=exp⁡ ⁣(A+B+∣α∣22I)=e∣α∣2/2D(α).e^Ae^B = \exp\!\left( A+B+\frac{|\alpha|^2}{2}I \right) = e^{|\alpha|^2/2}D(\alpha).

Therefore

D(α)=e−∣α∣2/2eαa†e−α∗a.D(\alpha) = e^{-|\alpha|^2/2} e^{\alpha a^\dagger} e^{-\alpha^*a}.
  1. Let H=HA+HBH=H_A+H_B. Find the leading correction in the exponent of e−itHA/ℏe−itHB/ℏe^{-itH_A/\hbar}e^{-itH_B/\hbar} relative to e−it(HA+HB)/ℏe^{-it(H_A+H_B)/\hbar}.
Solution

Choose

A=−iℏHA,B=−iℏHB.A=-\frac{i}{\hbar}H_A, \qquad B=-\frac{i}{\hbar}H_B.

BCH gives

log⁡(etAetB)=t(A+B)+t22[A,B]+O(t3).\log(e^{tA}e^{tB}) = t(A+B)+\frac{t^2}{2}[A,B]+O(t^3).

Since

[A,B]=−1ℏ2[HA,HB],[A,B] = -\frac{1}{\hbar^2}[H_A,H_B],

the leading correction is

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

It vanishes when the two Hamiltonian pieces commute.

  1. Use Hadamard’s lemma to evaluate eiλP/ℏXe−iλP/ℏe^{i\lambda P/\hbar}Xe^{-i\lambda P/\hbar} when [X,P]=iℏI[X,P]=i\hbar I.
Solution

Take A=iλP/ℏA=i\lambda P/\hbar. Then

[A,X]=iλℏ[P,X]=iλℏ(−iℏI)=λI.\begin{aligned} [A,X] &= \frac{i\lambda}{\hbar}[P,X] \\ &= \frac{i\lambda}{\hbar}(-i\hbar I) = \lambda I. \end{aligned}

All further nested commutators vanish because [A,I]=0[A,I]=0. Hence

eiλP/ℏXe−iλP/ℏ=X+λI.e^{i\lambda P/\hbar} X e^{-i\lambda P/\hbar} = X+\lambda I.