Demkov–Osherov Model (DOM) Overview

• Demkov–Osherov Model (DOM) is an exactly solvable multilevel Landau–Zener problem where a single linearly changing diabatic level crosses a band of constant-energy states.
• It employs a star-topology coupling that factorizes transition probabilities into independent two-level Landau–Zener processes, highlighting its integrable structure.
• The model underpins advanced studies in nonadiabatic scattering, integrable systems, and extensions to curved diabatic levels in multivariable Painlevé-II analyses.

The Demkov–Osherov model (DOM), often written as the DO model, is an exactly solvable multilevel Landau–Zener problem in which a single diabatic level with linear time dependence crosses a band of parallel or time-independent diabatic levels, with constant couplings only between the slanted level and each band state. In the standard reduction used in the literature, the Hamiltonian is obtained from a generic multilevel Landau–Zener form $H_{LZ}(t)=\hat A+\hat B\,t$ when $\hat B$ has one eigenvalue distinct from an $n$-fold degenerate subspace, after subtracting a trivial diagonal term and rescaling time. The model is a paradigmatic exactly solvable system of nonadiabatic transitions, and later work uses it both as a reference point for nonlinear generalizations with curved levels and as an exact local scattering problem in integrable constructions based on Knizhnik–Zamolodchikov equations and multivariable Painlevé-II systems (Sedrakyan et al., 2021, Lin et al., 2014, Sinitsyn, 23 Mar 2026).

1. Canonical Hamiltonian and reduction from multilevel Landau–Zener dynamics

A standard starting point is the generic multilevel Landau–Zener Hamiltonian

$H_{LZ}(t)=\hat A+\hat B\,t.$

If $\hat B$ has one eigenvalue $b_1$ and $n$ degenerate eigenvalues $b_2$, one diagonalizes $\hat B$, uses the degeneracy of the $n$-dimensional subspace to diagonalize the corresponding minor of $\hat B$0, subtracts the diagonal piece $\hat B$1, and rescales time by $\hat B$2. This yields the Demkov–Osherov Hamiltonian

$\hat B$3

or, in bra–ket form,

$\hat B$4

Only the level $\hat B$5 has a linearly time-dependent diabatic energy. The levels $\hat B$6, $\hat B$7, have constant diabatic energies $\hat B$8, and there is no direct coupling between $\hat B$9 and $n$0 for $n$1. In this basis, the model is precisely one moving level interacting with a band of static levels (Sedrakyan et al., 2021).

The two-level Landau–Zener problem is recovered for $n$2. In that sense, the DOM is the minimal multilevel extension that preserves exact solvability while retaining a highly structured coupling graph: a star topology centered on the sweeping level $n$3 (Sedrakyan et al., 2021).

2. Scattering problem, exact structure, and factorization

The Schrödinger equation is considered in the scattering picture,

$n$4

with evolution from $n$5 to $n$6, where the diabatic levels are asymptotically well separated. Transition probabilities are then defined in the diabatic basis $n$7 (Sedrakyan et al., 2021).

An algebraic characterization of the exact solution uses the eigenstate ansatz

$n$8

with consistency conditions

$n$9

and

$H_{LZ}(t)=\hat A+\hat B\,t.$0

These relations show that the parameters of the Hamiltonian admit a nontrivial integrable parametrization in terms of $H_{LZ}(t)=\hat A+\hat B\,t.$1 (Sedrakyan et al., 2021).

In the standard DOM, transition probabilities factorize in a Landau–Zener-like product form. For the three-level case with one sweeping level and two parallel levels, the exact DOM scattering amplitudes coincide with the result of the independent crossing approximation: one multiplies the individual two-level Landau–Zener scattering matrices in chronological order and inserts adiabatic phase accumulation between crossings. At each isolated two-level crossing with coupling $H_{LZ}(t)=\hat A+\hat B\,t.$2 and relative slope $H_{LZ}(t)=\hat A+\hat B\,t.$3, the corresponding Landau–Zener probability is

$H_{LZ}(t)=\hat A+\hat B\,t.$4

In the DOM setting used for asymptotic analysis of multivariable Painlevé-II, this factorized construction is not merely heuristic; it coincides with the exact DOM solution (Sinitsyn, 23 Mar 2026).

3. Curved-level generalization and the DOM limit

A later exactly solvable model replaces the band of parallel levels by a Coulomb band and can be regarded as a generalization of the DOM. In the original time variable $H_{LZ}(t)=\hat A+\hat B\,t.$5, the diabatic Hamiltonian is

$H_{LZ}(t)=\hat A+\hat B\,t.$6

with all other matrix elements zero. Thus the sweeping level remains linear,

$H_{LZ}(t)=\hat A+\hat B\,t.$7

while the band states become curved,

$H_{LZ}(t)=\hat A+\hat B\,t.$8

The DOM is recovered in the specific limit

$H_{LZ}(t)=\hat A+\hat B\,t.$9

for $\hat B$0. In that limit the Coulomb band behaves locally as a band of nearly parallel levels, and the Demkov–Osherov solution is recovered (Lin et al., 2014).

After the change of variables

$\hat B$1

the system reduces to

$\hat B$2

Elimination of the band amplitudes yields an $\hat B$3-th order ordinary differential equation of Meijer $\hat B$4-type for $\hat B$5. The exact survival probability of the initially populated level $\hat B$6 is

$\hat B$7

where $\hat B$8 are the real roots of the characteristic polynomial

$\hat B$9

This formula retains a product structure reminiscent of DOM, but the exponents are renormalized by the full many-level interaction through the roots $b_1$0 (Lin et al., 2014).

Two limiting regimes are particularly informative. In the degenerate-band case $b_1$1 for all $b_1$2,

$b_1$3

so the survival probability remains finite even for arbitrarily strong couplings and many band states. In the well-separated regime $b_1$4,

$b_1$5

and

$b_1$6

which reproduces the product of effectively independent two-level Landau–Zener-like contributions familiar from DOM. By contrast, the curved-band model also shows that transition probabilities within the band generally do not saturate asymptotically; only probabilities involving the isolated level $b_1$7 are well defined at large time (Lin et al., 2014).

4. Distinction from the two-level Demkov pulse model

The DOM should be distinguished from the two-level Demkov model used in optical Bloch dynamics. The latter is a pulsed two-state problem with constant detuning and exponential envelope,

$b_1$8

where $b_1$9 is a pure dephasing rate (Vasilev et al., 2014).

Model	Time dependence	Role
Standard DOM	One level linear in time; band levels constant or parallel	Multilevel nonadiabatic scattering
Two-level Demkov model	$n$0, $n$1	Exactly solvable pulsed two-state problem
Coulomb-band generalization	One linear level and $n$2	Nonlinear-time generalization of DOM

In the two-level Demkov model with dephasing, the Bloch equations are

$n$3

with

$n$4

The exact solution is obtained by reducing the problem to a third-order differential equation for $n$5, solving it in terms of generalized hypergeometric functions $n$6, and matching across the cusp at $n$7. On resonance, $n$8, the solution simplifies further and reduces to Bessel functions (Vasilev et al., 2014).

This two-level Demkov problem is not the multilevel Demkov–Osherov model, but it functions as one of the elementary building blocks underlying the general DOM: in the full DOM, transitions are described as a sequence of independent two-level Demkov- or Landau–Zener-type crossings between one slanted level and several parallel levels. With dephasing included, the final population transfer is

$n$9

and increasing $b_2$0 leads to a monotonic suppression of the final transition probability (Vasilev et al., 2014).

5. Integrable reformulations through EKZ equations and boundary WZNW theory

The DOM and its close variants admit a reformulation in terms of extended Knizhnik–Zamolodchikov equations and boundary Wess–Zumino–Novikov–Witten theory. In that construction, the boundary WZNW action is

$b_2$1

with contour terms

$b_2$2

The resulting correlators satisfy extended KZ equations

$b_2$3

where $b_2$4 are Gaudin Hamiltonians and $b_2$5 (Sedrakyan et al., 2021).

The paper then studies an altered Demkov–Osherov model (ADO), in which two levels have the same slope in time,

$b_2$6

After Fourier transformation in time and elimination of the static levels, one obtains an EKZ-type equation

$b_2$7

together with companion equations in the parameters $b_2$8. The integrability condition is the rank-one coupling constraint

$b_2$9

which implies that the corresponding classical vectors $\hat B$0 are parallel and ensures the zero-curvature condition $\hat B$1 (Sedrakyan et al., 2021).

Under this constraint the ADO system is solved exactly. For the simplest nontrivial case $\hat B$2, the transition probability takes the explicit Landau–Zener form

$\hat B$3

This reformulation places DOM-type Hamiltonians at the interface of non-equilibrium quantum dynamics, Gaudin-type integrable systems, and boundary conformal field theory (Sedrakyan et al., 2021).

6. DOM as local scattering data in multivariable Painlevé-II theory

A later application uses the exact DOM solution as the central linear ingredient in an asymptotically exact WKB analysis of a multivariable Painlevé-II system. For the $\hat B$4-component nonlinear equations

$\hat B$5

the consistency condition of a two-time Schrödinger problem defines a Lax pair with Hamiltonians $\hat B$6 and $\hat B$7. In the worked-out $\hat B$8 case, one studies the $\hat B$9-evolution at fixed large $n$0 in a three-dimensional diabatic basis $n$1 (Sinitsyn, 23 Mar 2026).

For large negative $n$2, after the rescaling

$n$3

and expansion near $n$4, the local Hamiltonians become exact three-level DOM Hamiltonians. Near $n$5,

$n$6

and near $n$7,

$n$8

Each local Hamiltonian has one linearly sweeping level, two parallel levels, and time-independent couplings within the local variable. In this framework the exact DOM S-matrix replaces conventional complex-plane Stokes data: path-independence of the flat two-time connection allows the evolution operator to be computed in different asymptotic regions, and the DOM solution supplies the nonadiabatic matching data (Sinitsyn, 23 Mar 2026).

The resulting connection formulas between $n$9 and $\hat B$00 are written in terms of

$\hat B$01

and phase combinations $\hat B$02. Their structure contains the characteristic DOM ingredients: Landau–Zener exponentials $\hat B$03, coherent interference of multiple semiclassical paths through the three-level crossing, and phase factors involving $\hat B$04 (Sinitsyn, 23 Mar 2026).

In the physical application to unstable vacuum decay through a second-order phase transition, the Painlevé-II connection formulas determine the scaling of produced excitations. In the near-vacuum regime

$\hat B$05

the phase-averaged quantities satisfy

$\hat B$06

$\hat B$07

and

$\hat B$08

Semiclassically, $\hat B$09 is the number of Higgs-like excitations and $\hat B$10 the number of Goldstone-like excitations. In this sense, DOM is not merely a solvable crossing model but an exact local scattering mechanism that determines asymptotic data for an integrable nonlinear system (Sinitsyn, 23 Mar 2026).