Hidden facts in Landau-Zener transitions revealed by the Riccati Equation
Pith reviewed 2026-05-23 04:14 UTC · model grok-4.3
The pith
The Riccati equation bridges the two Landau-Zener probability amplitudes and shows why the Markov approximation gives the exact asymptotic result for one but not the other.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The dynamics of the two probability amplitudes in the elementary Landau-Zener problem can be expressed in terms of the solution of the corresponding Riccati differential equation. This solution provides the bridge between the two amplitudes. Neglecting the nonlinearity in the Riccati equation is equivalent to the Markov approximation, which yields the exact asymptotic expression for one of the probability amplitudes. The Riccati equation itself identifies the origin of the failure of the Markov approximation to provide the correct asymptotic expression for the other amplitude. The approach relies on approximate yet analytical solutions of the Riccati equation in different time regimes.
What carries the argument
The Riccati differential equation obtained from the Landau-Zener amplitudes; its solution connects the two amplitudes while its nonlinearity controls the validity of the Markov approximation.
If this is right
- The Markov approximation supplies the exact asymptotic probability for one of the two final states.
- The retained nonlinear term in the Riccati equation is the precise reason the Markov approximation fails for the second amplitude.
- Separate analytic approximations to the Riccati solution in successive time regimes suffice to track the effect of nonlinearity on both final amplitudes.
- Once one amplitude is known, the bridge property of the Riccati solution determines the other.
Where Pith is reading between the lines
- The same Riccati reformulation could be applied to other two-level systems with linear time dependence to test whether the Markov-Markov distinction persists.
- Direct numerical integration of the Riccati equation itself would provide a quantitative check on the accuracy of the piecewise analytic approximations used for the asymptotics.
- The explicit identification of nonlinearity as the source of discrepancy suggests a systematic route for constructing improved, non-Markovian approximations in related driven two-level models.
Load-bearing premise
The approximate analytical solutions of the Riccati equation constructed separately in different time regimes correctly capture the impact of its nonlinear term on the asymptotic amplitudes.
What would settle it
A numerical solution of the original time-dependent Schrödinger equation for the Landau-Zener Hamiltonian whose long-time amplitudes deviate from the values obtained by applying the Markov approximation to the Riccati equation.
Figures
read the original abstract
We express the dynamics of the two probability amplitudes in the elementary Landau-Zener problem in terms of the solution of the corresponding Riccati differential equation and identify three key features: (i) The solution of the Riccati equation provides the bridge between the two probability amplitudes. (ii) Neglecting the non-linearity in the Riccati equation is equivalent to the Markov approximation which yields the exact asymptotic expression for one of the probability amplitudes, and (iii) the Riccati equation identifies the origin of the failure of the Markov approximation not being able to provide us in general with the correct asymptotic expression of the other probability amplitude. Our approach relies on approximate yet analytical solutions of the Riccati equation in different time regimes, highlighting the impact of its non-linear nature on the time evolution of the system.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reformulates the Landau-Zener (LZ) two-level dynamics in terms of the solution to the associated Riccati differential equation. It asserts three features: (i) the Riccati solution bridges the two probability amplitudes, (ii) dropping the nonlinear term is equivalent to the Markov approximation and recovers the exact asymptotic expression for one amplitude, and (iii) the retained nonlinearity explains why the Markov approximation fails to give the correct asymptotic for the second amplitude. These conclusions are reached by constructing approximate analytical solutions of the Riccati equation separately in distinct time regimes and matching them.
Significance. If the central claims are rigorously established, the work supplies a compact mathematical link between the two LZ amplitudes and a transparent origin for the selective success of the Markov approximation. This perspective could be useful for analyzing other time-dependent two-level problems where nonlinearity in an auxiliary equation controls asymptotic accuracy. The significance is currently tempered by the absence of explicit error bounds or direct numerical/exact comparisons that would confirm the piecewise approximations transmit the nonlinear contribution without uncontrolled artifacts.
major comments (2)
- [Abstract] Abstract (final sentence) and the associated derivations: the claim that the regime-specific approximate solutions of the Riccati equation correctly transmit the effect of the nonlinear term onto the asymptotic amplitudes is load-bearing for all three listed features. No independent error estimate, matching-condition analysis, or direct comparison to the known exact LZ asymptotics (e.g., the standard exponential form for the transition probability) is supplied to substantiate that the piecewise construction does not introduce or mask nonlinear contributions.
- [Abstract] Abstract (item ii): the asserted equivalence between the linear (Markov) truncation of the Riccati equation and the exact asymptotic expression for one amplitude requires an explicit derivation showing that the neglected term vanishes identically in the relevant limit; without this step the identification of the 'origin of the failure' for the second amplitude does not follow.
minor comments (1)
- [Abstract] Abstract, last sentence: the phrasing 'the failure of the Markov approximation not being able to provide us in general with the correct asymptotic expression' is grammatically awkward and should be revised for clarity.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. We address each major comment below and indicate the revisions we will incorporate.
read point-by-point responses
-
Referee: [Abstract] Abstract (final sentence) and the associated derivations: the claim that the regime-specific approximate solutions of the Riccati equation correctly transmit the effect of the nonlinear term onto the asymptotic amplitudes is load-bearing for all three listed features. No independent error estimate, matching-condition analysis, or direct comparison to the known exact LZ asymptotics is supplied to substantiate that the piecewise construction does not introduce or mask nonlinear contributions.
Authors: We acknowledge that the manuscript would be strengthened by explicit validation of the piecewise construction. The matching is performed in overlapping time intervals where both approximations remain valid, which is intended to preserve continuity and carry the nonlinear contribution forward. To address the concern directly, we will add a matching-condition analysis together with direct comparisons of the resulting asymptotic amplitudes against the known exact Landau-Zener expressions and numerical integrations of the Riccati equation. revision: yes
-
Referee: [Abstract] Abstract (item ii): the asserted equivalence between the linear (Markov) truncation of the Riccati equation and the exact asymptotic expression for one amplitude requires an explicit derivation showing that the neglected term vanishes identically in the relevant limit; without this step the identification of the 'origin of the failure' for the second amplitude does not follow.
Authors: We agree that an explicit demonstration of the vanishing is necessary for rigor. In the appropriate large-time limit the nonlinear term is suppressed by the asymptotic form of the solution for the amplitude in question, recovering the exact Markov result. We will insert a dedicated derivation of this limit in the revised manuscript to make the equivalence and the contrasting behavior of the second amplitude fully transparent. revision: yes
Circularity Check
Derivation self-contained from standard LZ setup; no circular reductions
full rationale
The paper starts from the elementary Landau-Zener Hamiltonian, expresses the amplitudes via the Riccati equation, and constructs approximate analytical solutions in separate time regimes to identify the listed features. No self-citations, parameter fittings, ansatzes smuggled via prior work, or renamings of known results appear as load-bearing steps. The equivalence of linear truncation to Markov approximation and the attribution of failure to nonlinearity follow directly from the equation structure and regime-specific approximations without reducing to the target asymptotics by definition. This is the common case of an independent mathematical reformulation.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The Landau-Zener problem is governed by a pair of coupled first-order linear ODEs for the probability amplitudes.
- ad hoc to paper Approximate analytical solutions of the Riccati equation can be constructed separately in distinct time regimes and then matched.
Reference graph
Works this paper leans on
-
[1]
(33), of η into the ansatz, eq
Probability amplitude a We start our analysis by considering the probability amplitude a and substitute this representation, eq. (33), of η into the ansatz, eq. (26), and arrive at the identity a(τ) = eiφ(τ)A(τ) (34) where the amplitude A(τ) ≡ e− Re[H(τ)] = exp − τZ −τ0 dτ ′η(R)(τ ′) (35) is given by the negative exponential of the integral of the ...
-
[2]
(3) as follows b = i˙a + ϵτ a, (37) and substitute eq
Connection between a and b by η In order to find an expression for b we rearrange eq. (3) as follows b = i˙a + ϵτ a, (37) and substitute eq. (28) into this expression. We arrive at the representation b = −iηa (38) for the probability amplitude for b in terms of the complex- valued function η and the probability amplitude a. Hence, η given by the Riccati d...
-
[3]
(31), representing two real-valued functions
Summary In summary, we have shown that the two complex- valued probability amplitudes, a and b, corresponding to four real-valued functions, are governed by the complex- valued Riccati differential equation, eq. (31), representing two real-valued functions. For the specific initial condition, eq. (32), the time-integral of the real part η(R) and the phase...
-
[4]
(38), for b, we have used the differential equation, eq
Exponential representation In the derivation of the expression, eq. (38), for b, we have used the differential equation, eq. (3), along with the ansatz, eq. (26), for a. Equivalently, we can formally integrate the differential equation, eq. (4), for b subjected to the initial condition, eq. (6), to arrive at the formula b(τ) = −ie−iϵτ 2/2 τZ −τ0 dτ ′eiϵτ ...
-
[5]
(55) due to the presence of η in the exponent of eq
Approximations of the non-linear integral equation It is interesting that η2 appears in eq. (55) due to the presence of η in the exponent of eq. (54). When we neglect η in this exponent, that is the non-linearity in the Riccati equation, eq. (55), we arrive at the linear differential equation ˙ηM = −2iϵτ ηM + 1 (56) with the solution ηM(τ) ≡ e−iϵτ 2 τZ −τ...
-
[6]
(54), into the form η = i b a (60) following from eq
A different integral equation for η In order to cast η given by the non-linear integral equation, eq. (54), into the form η = i b a (60) following from eq. (38), we take the term a−1(τ) out of the integral which yields the expression η(τ) = e−iϵτ 2 exp τZ −τ0 dτ ′η(τ ′) × τZ −τ0 dτ ′eiϵτ ′2 exp − τ ′ Z −τ0 dτ ′′η(τ ′′) . (61) Next, we decom...
-
[7]
An integral identity Indeed, we now use the non-linear integral equation, eq. (62), to establish the identity |I(τ)| = p 1 − A2(τ) (64) 9 where we have introduced the integral I(τ) ≡ τZ −τ0 dτ ′eiλ(τ ′)A(τ ′). (65) From the derivative d dτ |I|2 = d dτ (II ∗) = ˙II ∗ + I ˙I ∗ (66) of the absolute value squared of I, we find with the help of eq. (65) the id...
-
[8]
(62), for η and the representation eq
Summary We are now in the position to make the connection between the non-linear integral equation, eq. (62), for η and the representation eq. (47). Indeed, the identity, eq. (64), leads to the representation e−iλ(τ) τZ −τ0 dτ ′eiλ(τ ′)A(τ ′) = p 1 − A(τ)2eiϕη(τ) (73) where the phase ϕη(τ) is determined by the argument of the left-hand side of this equati...
-
[9]
(49), for b into the differential equation, eq
Schr¨ odinger picture When we insert the formal solution, eq. (49), for b into the differential equation, eq. (3), we obtain the differential equation ˙a(τ) = iϵτ a(τ) − e−iϵτ 2/2 τZ −τ0 dτ ′eiϵτ ′2/2a(τ ′) (76) of Lippmann-Schwinger type [28] for the probability am- plitude a. Up to this point, our derivation remains exact, and no approximation has been ...
-
[10]
Interaction picture Next, we compare and contrast this expression with the one derived from the Markov approximation in the interaction picture, which was previously obtained in Ref. [24]. From eq. (13), we find with the initial condition, eq. (15), the expression for ˜b which when substituted into eq. (12) yields the integral differential equation ˙˜a(τ)...
-
[11]
(32), at a finite but large negative time τ = −τ0
Emergence of oscillations We start our discussion with the case of a vanishing initial condition, eq. (32), at a finite but large negative time τ = −τ0. Given that η is initially zero, we expect it to remain small during the early-time dynamics. This assumption allows us to linearize the Riccati equation by Figure 7. Oscillations in the functions η and ηM...
-
[12]
No oscillations Deeper insight into the oscillations emerges when we take the first derivative in time of of ηM given by eq. (98). The first term, proportional to 1 /τ turns into 1/τ 2 which we neglect. However, due to the quadratic chirp the oscillatory term creates an additional term linear in τ, and as a result, the amplitude is now of the order τ /τ0,...
-
[13]
Connection formula From the definition, eq. (79), of ηM we find the expres- sion ηM(−τ) = e−iϵτ 2 −τZ −τ0 dτ ′eiϵτ ′2 (113) which with the new integration variable ¯τ ≡ − τ ′ takes the form ηM(−τ) = e−iϵτ 2 τ0Z τ d¯τ eiϵ¯τ 2 . (114) Next, we subtract and add an appropriate integral to 16 obtain the relation ηM(−τ) = e−iϵτ 2 − τZ −τ0 d¯τ eiϵ¯τ 2 + τ0Z...
-
[14]
(117), of the Markov solution, eq
Stueckelberg oscillations Two contributions appear in the connection formula, eq. (117), of the Markov solution, eq. (80): (i) For large positive times the term −ηM(−|τ |) corresponds to ηM at large negative times, and (ii) phase factor quadratic in time. We now show that this phase factor is the origin of the Stueckelberg oscillations. In this derivation...
-
[15]
Exact Landau-Zener result We are now in the position to rederive the exact Landau- Zener result [1–5, 24, 35], eq. (22). For this purpose, we set in eq. (121) τ = τ0 and use the initial condition, eq. (82), which yields aM(τ = τ0) = exp −1 2 |F(τ0)|2 . (122) Next, we perform the limit τ0 → ∞ which with the integral relation F(τ0 → ∞) = ∞Z −∞ dτ eiϵτ 2 = r...
-
[16]
(49) the Markov approximation, eq
Markov approximation for b To bring out this fact most clearly, we perform in eq. (49) the Markov approximation, eq. (77), which yields bM(τ) = −ie−iϵτ 2 τZ −τ0 dτ ′eiϵτ ′2 aM(τ), (125) and with the definition, eq. (79), of ηM, we arrive at the expression bM = −iηM aM (126) for the probability amplitude bM in Markov approxima- tion. We note that we also a...
-
[17]
(49) the probability amplitude a by the exponential ansatz, eq
No factorization of the Markov approximation The origin of the failure of the Markov approximation to predict correctly the asymptotic value of b stands out most clearly when we express in eq. (49) the probability amplitude a by the exponential ansatz, eq. (26), which yields b(τ) = −ie−iϵτ 2 0 /2 × e−iϵτ 2/2 τZ −τ0 dτ ′eiϵτ ′2 exp − τ ′ Z −τ0 dτ ′′η(τ ...
-
[18]
(A2) We emphasize, that z differs from the corresponding quan- tity in Ref
Expression for probability amplitude a We start our discussion by recalling the linear second order differential equation d2 dz2 a(z) − z2 4 + i 2ϵ + 1 2 a(z) = 0 (A1) for the probability amplitudes a with the complex scaling z ≡ √ 2ϵτei π 4 . (A2) We emphasize, that z differs from the corresponding quan- tity in Ref. [35] by the factor √ϵ. According to R...
-
[19]
(3) which in the complex scaling, eq
Expression for probability amplitude b In order to derive an expression for b we recall eq. (3) which in the complex scaling, eq. (A2), reads i d dz a(z) = i z 2 a(z) − e−i π 4 √ 2ϵ b(z), (A26) and cast it into the form b(z) = √ 2ϵe−i π 4 d dz a(z) − z 2 a(z) . (A27) Together with the expressions, eq. (A3) and eq. (A8), for a and its derivative, we find b...
-
[20]
Expression for solution η of the Riccati equation In section III, we have presented an ansatz, eq. (26), for the probability amplitude a introducing the function η which leads us directly to the non-linear differential equation of first order, eq. (31), for η of the Riccati type. We now derive a representation of η in terms of the parabolic cylinder funct...
-
[21]
Asymptotics of a, b and η We are now in the position to evaluate the expressions for a, b and η in the limit of τ0 → ∞ using the appropriate asymptotic expansions of the parabolic cylinder functions. For this purpose, we first collect the corresponding ex- pressions, and then apply them to the formulae for a, b and η derived in the preceding sections of t...
-
[22]
For this purpose, we cast the equations, eq
Instantaneous eigenstates We start by recalling the expressions of the instanta- neous eigenvalues and eigenstates. For this purpose, we cast the equations, eq. (3) and eq. (4), for the probability amplitudes a and b into the matrix form i d dτ a(τ) b(τ) = −ϵτ 1 1 ϵτ a(τ) b(τ) (B1) which leads us to the instantaneous eigenvalues θ±(τ) ≡ ±θ(τ) ≡ ± q (ϵτ)2 ...
-
[23]
(2), using the instantaneous eigenstates, eq
Diagonalization Next, we diagonalize the Hamiltonian, eq. (2), using the instantaneous eigenstates, eq. (B3). Needless to say, this technique is only correct under appropriate conditions. In particular, we show it is valid for large negative times which represents the domain of interest in the present discussion. For the diagonalization we now find the tr...
-
[24]
(B15) only when we can neglect the time derivative of G governing the last term
Solution for the diagonalized system We emphasize that we can solve eq. (B15) only when we can neglect the time derivative of G governing the last term. Since in the limit of large negative times the matrix G(−|τ |) ∼= 1 − 1 2ϵ|τ | 1 2ϵ|τ | 1 ! (B18) is almost constant we can indeed neglect it. Therefore, we solve the equation i d dτ ˜a ˜b = D ˜a ˜b (B19)...
-
[25]
(B30) Hence, the difference of v+ and v− is given by v+ − v− = 2ϵ|τ | + 1 2ϵ|τ |
Asymptotic behavior for large negative times Since we are interested in an approximate expression for a in the neighborhood of τ = −τ0 we expand θ and v± for large negative times, and find the expressions θ±(−|τ |) ∼= ± ϵ|τ | + 1 2ϵ|τ | (B28) as well as v+(−|τ |) ∼= 2ϵ|τ | (B29) and v−(−|τ |) ∼= − 1 2ϵ|τ | . (B30) Hence, the difference of v+ and v− is giv...
- [26]
-
[27]
L. D. Landau, Sov. Phys. 2, 46 (1932)
work page 1932
-
[28]
L. D. Landau, Collected Papers of L.D. Landau , 63 (1965)
work page 1965
- [29]
-
[30]
E. C. G. Stueckelberg, Helvetica Physica Acta 5, 369 (1932)
work page 1932
-
[31]
D. L. Hill and J. A. Wheeler, Phys. Rev. 89, 1102 (1953)
work page 1953
-
[32]
E. Peik, M. Ben Dahan, I. Bouchoule, Y. Castin, and C. Salomon, Phys. Rev. A 55, 2989 (1997)
work page 1997
-
[33]
V. V. Konotop, P. G. Kevrekidis, and M. Salerno, Phys. Rev. A 72, 023611 (2005)
work page 2005
-
[34]
P. Clad´ e, S. Guellati-Kh´ elifa, F. Nez, and F. Biraben, Phys. Rev. Lett. 102, 240402 (2009)
work page 2009
- [35]
-
[36]
D. Heim, W. Schleich, P. Alsing, J. Dahl, and S. Varro, Physics Letters A 377, 1822 (2013)
work page 2013
- [37]
-
[38]
M. Gebbe, J.-N. Siemß, M. Gersemann, H. M¨ untinga, S. Herrmann, C. L¨ ammerzahl, H. Ahlers, N. Gaaloul, C. Schubert, K. Hammerer, S. Abend, and E. M. Rasel, Nature Communications 12, 2544 (2021)
work page 2021
-
[39]
P. O. Kofman, O. V. Ivakhnenko, S. N. Shevchenko, and F. Nori, Scientific Reports 13, 5053 (2023)
work page 2023
-
[40]
P. O. Kofman, S. N. Shevchenko, and F. Nori, Phys. Rev. A 109, 022409 (2024)
work page 2024
-
[41]
B. Konrad and M. Efremov, Quantum Science and Tech- nology 10, 015025 (2024)
work page 2024
-
[42]
V. M. Akulin and W. P. Schleich, Phys. Rev. A 46, 4110 (1992)
work page 1992
- [43]
-
[44]
N. V. Vitanov, Phys. Rev. A 59, 988 (1999)
work page 1999
-
[45]
J. Liu, L. Fu, B.-Y. Ou, S.-G. Chen, D.-I. Choi, B. Wu, and Q. Niu, Phys. Rev. A 66, 023404 (2002)
work page 2002
-
[46]
M. Wubs, K. Saito, S. Kohler, Y. Kayanuma, and P. H¨ anggi, New Journal of Physics7, 218 (2005)
work page 2005
-
[47]
M. Wubs, K. Saito, S. Kohler, P. H¨ anggi, and Y. Kayanuma, Phys. Rev. Lett. 97, 200404 (2006)
work page 2006
-
[48]
A. G. Rojo, Matrix exponential solution of the Landau- Zener problem (2010), arXiv:1004.2914 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[49]
E. P. Glasbrenner and W. P. Schleich, Journal of Physics B: Atomic, Molecular and Optical Physics 56, 104001 (2023)
work page 2023
- [50]
- [51]
-
[52]
M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- bridge University Press, 1997)
work page 1997
-
[53]
V. Paulisch, H. Rui, H. K. Ng, and B.-G. Englert, The European Physical Journal Plus 129, 12 (2014)
work page 2014
-
[54]
V. V. Zheleznyakov, V. V. Kocharovski˘ ı, and V. V. Kocharovski˘ ı, Soviet Physics Uspekhi26, 877 (1983)
work page 1983
-
[55]
Gardas, Journal of Mathematical Physics 51, 062103 (2010)
B. Gardas, Journal of Mathematical Physics 51, 062103 (2010)
work page 2010
-
[56]
S. Bittanti, A. Laub, and J. Willems, The Riccati Equa- tion, Communications and Control Engineering (Springer Berlin Heidelberg, 2012)
work page 2012
-
[57]
Schuch, Journal of Physics: Conference Series 538, 012019 (2014)
D. Schuch, Journal of Physics: Conference Series 538, 012019 (2014)
work page 2014
- [58]
-
[59]
M. Abramowitz and I. Stegun, Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965)
work page 1965
-
[60]
O. V. Ivakhnenko, S. N. Shevchenko, and F. Nori, Physics Reports 995, 1 (2023)
work page 2023
-
[61]
L. A. J. Guttieres, M. D. Petrovi´ c, and J. K. Freericks, American Journal of Physics 91, 885 (2023)
work page 2023
-
[62]
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nel- son, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Polat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perk- told, R. Cimrman, I. Henr...
work page 2020
-
[63]
J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, SIAM Review 59, 65 (2017)
work page 2017
-
[64]
A. Kazantsev, G. Surdutovich, and V. Yakovlev,Mechani- cal Action Of Light On Atoms (World Scientific Publishing Company, 1990)
work page 1990
- [65]
-
[66]
E. P. Glasbrenner, Y. Gerdes, S. Varr´ o, and W. P. Schle- ich, in Proceedings of the Second International Work- shop on Quantum Nonstationary Systems , edited by A. Dodonov and C. C. H. Ribeiro (LF Editorial, 2024) 1st ed., Chap. 17, pp. 273–297
work page 2024
-
[67]
J. S. Briggs and A. Eisfeld, Phys. Rev. A 88, 062104 28 (2013)
work page 2013
-
[68]
I. Gradshteyn and I. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 2014)
work page 2014
-
[69]
K.-A. Suominen, B. Garraway, and S. Stenholm, Optics Communications 82, 260 (1991)
work page 1991
-
[70]
V. Kocharovsky, V. Kocharovsky, E. Derishev, S. Litvak, I. Shereshevsky, and S. Tasaki, Computers and Mathe- matics with Applications 34, 727 (1997)
work page 1997
- [71]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.