Landau–Zener formula

This article may be too technical for most readers to understand. Please help improve it to make it understandable to non-experts, without removing the technical details. (September 2018) (Learn how and when to remove this message)

The Landau–Zener formula is an analytic solution to the equations of motion governing the transition dynamics of a two-state quantum system, with a time-dependent Hamiltonian varying such that the energy separation of the two states is a linear function of time. The formula, giving the probability of a diabatic (not adiabatic) transition between the two energy states, was published separately by Lev Landau,^[1] Clarence Zener,^[2] Ernst Stueckelberg,^[3] and Ettore Majorana,^[4] in 1932.

If the system starts, in the infinite past, in the lower energy eigenstate, we wish to calculate the probability of finding the system in the upper energy eigenstate in the infinite future (a so-called Landau–Zener transition). For infinitely slow variation of the energy difference (that is, a Landau–Zener velocity of zero), the adiabatic theorem tells us that no such transition will take place, as the system will always be in an instantaneous eigenstate of the Hamiltonian at that moment in time. At non-zero velocities, transitions occur with probability as described by the Landau–Zener formula.

Such transitions occur between states of the entire system; hence any description of the system must include all external influences, including collisions and external electric and magnetic fields. In order that the equations of motion for the system might be solved analytically, a set of simplifications are made, known collectively as the Landau–Zener approximation. The simplifications are as follows:

1. The perturbation parameter in the Hamiltonian is a known, linear function of time
2. The energy separation of the diabatic states varies linearly with time
3. The coupling in the diabatic Hamiltonian matrix is independent of time

The first simplification makes this a semi-classical treatment. In the case of an atom in a magnetic field, the field strength becomes a classical variable which can be precisely measured during the transition. This requirement is quite restrictive as a linear change will not, in general, be the optimal profile to achieve the desired transition probability.

The second simplification allows us to make the substitution

${\displaystyle \Delta E=E_{2}(t)-E_{1}(t)\equiv \alpha t,\,}$

where ${\displaystyle E_{1}(t)}$ and ${\displaystyle E_{2}(t)}$ are the energies of the two states at time t, given by the diagonal elements of the Hamiltonian matrix, and ${\displaystyle \alpha }$ is a constant. For the case of an atom in a magnetic field this corresponds to a linear change in magnetic field. For a linear Zeeman shift this follows directly from point 1.

The final simplification requires that the time–dependent perturbation does not couple the diabatic states; rather, the coupling must be due to a static deviation from a ${\displaystyle 1/r}$ Coulomb potential, commonly described by a quantum defect.

The details of Zener's solution are somewhat opaque, relying on a set of substitutions to put the equation of motion into the form of the Weber equation^[5] and using the known solution. A more transparent solution is provided by Curt Wittig^[6] using contour integration.

The key figure of merit in this approach is the Landau–Zener velocity:

${\displaystyle v_{\rm {LZ}}={{\frac {\partial }{\partial t}}|E_{2}-E_{1}| \over {\frac {\partial }{\partial q}}|E_{2}-E_{1}|}\approx {\frac {dq}{dt}},}$

where q is the perturbation variable (electric or magnetic field, molecular bond-length, or any other perturbation to the system), and ${\displaystyle E_{1}}$ and ${\displaystyle E_{2}}$ are the energies of the two diabatic (crossing) states. A large ${\displaystyle v_{\rm {LZ}}}$ results in a large diabatic transition probability and vice versa.

Using the Landau–Zener formula the probability, ${\displaystyle P_{\rm {D}}}$, of a diabatic transition is given by

${\displaystyle {\begin{aligned}P_{\rm {D}}&=e^{-2\pi \Gamma }\\\Gamma &={a^{2}/\hbar \over \left|{\frac {\partial }{\partial t}}(E_{2}-E_{1})\right|}={a^{2}/\hbar \over \left|{\frac {dq}{dt}}{\frac {\partial }{\partial q}}(E_{2}-E_{1})\right|}\\&={a^{2} \over \hbar |\alpha |}\end{aligned}}}$

The quantity ${\displaystyle a}$ is the off-diagonal element of the two-level system's Hamiltonian coupling the bases, and as such it is half the distance between the two unperturbed eigenenergies at the avoided crossing, when ${\displaystyle E_{1}=E_{2}}$.

The simplest generalization of the two-state Landau–Zener model is a multistate system with a Hamiltonian of the form

${\displaystyle H(t)=A+Bt}$,

where A and B are Hermitian NxN matrices with time-independent elements. The goal of the multistate Landau–Zener theory is to determine elements of the scattering matrix and the transition probabilities between states of this model after evolution with such a Hamiltonian from negative infinite to positive infinite time. The transition probabilities are the absolute value squared of scattering matrix elements.

There are exact formulas, called hierarchy constraints, that provide analytical expressions for special elements of the scattering matrix in any multi-state Landau–Zener model.^[7] Special cases of these relations are known as the Brundobler–Elser (BE) formula,^[8][9][10]), and the no-go theorem,.^[11][12] Discrete symmetries often lead to constraints that reduce the number of independent elements of the scattering matrix.^[13][14]

There are also integrability conditions that, when they are satisfied, lead to exact expressions for the entire scattering matrices in multistate Landau–Zener models. Numerous such completely solvable models have been identified, including:

• Demkov–Osherov model^[15] that describes a single level that crosses a band of parallel levels. A surprising fact about the solution of this model is coincidence of the exactly obtained transition probability matrix with its form obtained with a simple semiclassical independent crossing approximation. With some generalizations, this property appears in almost all solvable Landau–Zener systems with a finite number of interacting states.
• Generalized bow-tie model.^[16] The model describes coupling of two (or one in the degenerate case limit) levels to a set of otherwise noninteracting diabatic states that cross at a single point.
• Driven Tavis–Cummings model^[17] describes interaction of N spins-⁠1/2⁠ with a bosonic mode in a linearly time-dependent magnetic field. This is the richest known solved system. It has combinatorial complexity: the dimension of its state vector space is growing exponentially with the number of spins N. The transition probabilities in this model are described by the q-deformed binomial statistics.^[18] This solution has found practical applications in physics of Bose-Einstein condensates.^[19]
• Spin clusters interacting with time-dependent magnetic fields.^[20] This class of models shows relatively complex behavior of the transition probabilities due to the path interference effects in the semiclassical independent crossing approximation.
• Reducible (or composite) multistate Landau–Zener models.^[21][22] This class consists of systems that can be decoupled to subsets of other solvable and simpler models by a symmetry transformation. The notable example is an arbitrary spin Hamiltonian ${\textstyle H=gS_{x}+btS_{z}}$, where S_z and S_x are spin operators, and S>1/2; b and g are constant parameters. This is the earliest known solvable system, which was discussed by Majorana in 1932. Among the other examples there are models of a pair of degenerate level crossing,^[23] and the 1D quantum Ising chain in a linearly changing magnetic field.^[24][25]
• Landau–Zener transitions in infinite linear chains.^[26] This class contains the systems with formally infinite number of interacting states. Although most known their instances can be obtained as limits of the finite size models (such as the Tavis–Cummings model), there are also cases that do not belong to this classification. For example, there are solvable infinite chains with nonzero couplings between non-nearest states.^[27]

Applications of the Landau–Zener solution to the problems of quantum state preparation and manipulation with discrete degrees of freedom stimulated the study of noise and decoherence effects on the transition probability in a driven two-state system. Several compact analytical results have been derived to describe these effects, including the Kayanuma formula ^[28] for a strong diagonal noise, and Pokrovsky–Sinitsyn formula ^[29] for the coupling to a fast colored noise with off-diagonal components.

Using the Schwinger–Keldysh Green's function, a rather complete and comprehensive study on the effect of quantum noise in all parameter regimes were performed by Ao and Rammer in late 1980s, from weak to strong coupling, low to high temperature, slow to fast passage, etc. Concise analytical expressions were obtained in various limits, showing the rich behaviors of such problem. ^[30] The effects of nuclear spin bath and heat bath coupling on the Landau–Zener process was explored by Sinitsyn and Prokof'ev^[31] and Pokrovsky and Sun,^[32][33][34] respectively.

Exact results in multistate Landau–Zener theory (no-go theorem and BE-formula) can be applied to Landau–Zener systems which are coupled to baths composed of infinite many oscillators and/or spin baths (dissipative Landau–Zener transitions). They provide exact expressions for transition probabilities averaged over final bath states if the evolution begins from the ground state at zero temperature, see in Ref. for oscillator baths^[35] and for universal results including spin baths in Ref.^[36]

• Nonadiabatic transition state theory
• Adiabatic theorem
• Bond softening
• Bond hardening
• Froissart-Stora equation