pith. sign in

arxiv: 2502.04033 · v1 · pith:5FKR4MUDnew · submitted 2025-02-06 · 🪐 quant-ph

Hidden facts in Landau-Zener transitions revealed by the Riccati Equation

Pith reviewed 2026-05-23 04:14 UTC · model grok-4.3

classification 🪐 quant-ph
keywords Landau-Zener transitionRiccati equationMarkov approximationprobability amplitudesasymptotic behaviorquantum dynamicstwo-level systemnonlinear differential equation
0
0 comments X

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.

In the Landau-Zener problem the time evolution is carried by two probability amplitudes. Their dynamics can be rewritten exactly in terms of the solution of a Riccati differential equation. Dropping the nonlinear term in that equation is identical to the Markov approximation and produces the precise long-time value of one amplitude. The same Riccati equation locates the source of the Markov approximation's failure to recover the correct asymptotic value of the second amplitude. These relations are obtained from approximate analytic solutions of the Riccati equation constructed in successive time intervals.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2502.04033 by Eric P. Glasbrenner, S\'andor Varr\'o, Wolfgang P. Schleich, Yannik Gerdes.

Figure 1
Figure 1. Figure 1: Time dependence of the probability amplitudes [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Time dependence of the probability amplitudes [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Numerical solution η of the Riccati differential equation, eq. (31), subjected to the initial condition eq. (32), visualized as a trajectory in the complex plane. The trajectory starts at τ = −τ0 from the origin marked by the grey cross and approaches rapidly without circling the origin an oscillatory regime. Indeed, the dark grey dot indicates the time τ = 0. The grey unit circle serves as a guide for the… view at source ↗
Figure 5
Figure 5. Figure 5: Emergence of the time dependence of the two probability amplitudes [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison between the exact numerical solution [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Oscillations in the early time evolution of the proba [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparison between the numerical evaluation of [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison between the numerical solution [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard two-level LZ Hamiltonian and the validity of piecewise analytical approximations to the resulting Riccati equation; no free parameters or new entities are introduced.

axioms (2)
  • domain assumption The Landau-Zener problem is governed by a pair of coupled first-order linear ODEs for the probability amplitudes.
    Standard starting point in quantum mechanics; invoked implicitly when the Riccati reformulation is introduced.
  • ad hoc to paper Approximate analytical solutions of the Riccati equation can be constructed separately in distinct time regimes and then matched.
    Explicitly stated in the final sentence of the abstract as the method used to highlight nonlinearity effects.

pith-pipeline@v0.9.0 · 5681 in / 1316 out tokens · 46254 ms · 2026-05-23T04:14:32.961606+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

71 extracted references · 71 canonical work pages · 1 internal anchor

  1. [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. [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. [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. [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. [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. [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. [7]

    (62), to establish the identity |I(τ)| = p 1 − A2(τ) (64) 9 where we have introduced the integral I(τ) ≡ τZ −τ0 dτ ′eiλ(τ ′)A(τ ′)

    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. [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. [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. [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. [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. [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. [13]

    (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

    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. [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. [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. [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. [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. [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. [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. [20]

    (26), for the probability amplitude a introducing the function η which leads us directly to the non-linear differential equation of first order, eq

    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. [21]

    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 this appendix

    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. [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. [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. [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. [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. [26]

    Zener, Proc

    C. Zener, Proc. R. Soc. London A 137, 696 (1932)

  27. [27]

    L. D. Landau, Sov. Phys. 2, 46 (1932)

  28. [28]

    L. D. Landau, Collected Papers of L.D. Landau , 63 (1965)

  29. [29]

    Majorana, Il Nuovo Cimento 9, 43 (1932)

    E. Majorana, Il Nuovo Cimento 9, 43 (1932)

  30. [30]

    E. C. G. Stueckelberg, Helvetica Physica Acta 5, 369 (1932)

  31. [31]

    D. L. Hill and J. A. Wheeler, Phys. Rev. 89, 1102 (1953)

  32. [32]

    E. Peik, M. Ben Dahan, I. Bouchoule, Y. Castin, and C. Salomon, Phys. Rev. A 55, 2989 (1997)

  33. [33]

    V. V. Konotop, P. G. Kevrekidis, and M. Salerno, Phys. Rev. A 72, 023611 (2005)

  34. [34]

    Clad´ e, S

    P. Clad´ e, S. Guellati-Kh´ elifa, F. Nez, and F. Biraben, Phys. Rev. Lett. 102, 240402 (2009)

  35. [35]

    Shevchenko, S

    S. Shevchenko, S. Ashhab, and F. Nori, Physics Reports 492, 1 (2010)

  36. [36]

    D. Heim, W. Schleich, P. Alsing, J. Dahl, and S. Varro, Physics Letters A 377, 1822 (2013)

  37. [37]

    Pagel, W

    Z. Pagel, W. Zhong, R. H. Parker, C. T. Olund, N. Y. Yao, and H. M¨ uller, Phys. Rev. A102, 053312 (2020)

  38. [38]

    Gebbe, J.-N

    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)

  39. [39]

    P. O. Kofman, O. V. Ivakhnenko, S. N. Shevchenko, and F. Nori, Scientific Reports 13, 5053 (2023)

  40. [40]

    P. O. Kofman, S. N. Shevchenko, and F. Nori, Phys. Rev. A 109, 022409 (2024)

  41. [41]

    Konrad and M

    B. Konrad and M. Efremov, Quantum Science and Tech- nology 10, 015025 (2024)

  42. [42]

    V. M. Akulin and W. P. Schleich, Phys. Rev. A 46, 4110 (1992)

  43. [43]

    Suominen and B

    K.-A. Suominen and B. M. Garraway, Phys. Rev. A 45, 374 (1992)

  44. [44]

    N. V. Vitanov, Phys. Rev. A 59, 988 (1999)

  45. [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)

  46. [46]

    M. Wubs, K. Saito, S. Kohler, Y. Kayanuma, and P. H¨ anggi, New Journal of Physics7, 218 (2005)

  47. [47]

    M. Wubs, K. Saito, S. Kohler, P. H¨ anggi, and Y. Kayanuma, Phys. Rev. Lett. 97, 200404 (2006)

  48. [48]

    A. G. Rojo, Matrix exponential solution of the Landau- Zener problem (2010), arXiv:1004.2914 [quant-ph]

  49. [49]

    E. P. Glasbrenner and W. P. Schleich, Journal of Physics B: Atomic, Molecular and Optical Physics 56, 104001 (2023)

  50. [50]

    J. R. F. Lima and G. Burkard, Superadiabatic Landau- Zener transitions (2024), arXiv:2408.03173 [quant-ph]

  51. [51]

    Brewer, E

    W. Brewer, E. Fick, and G. Sauermann, The Quantum Statistics of Dynamic Processes , Springer Series in Solid- State Sciences (Springer Berlin Heidelberg, 1990)

  52. [52]

    M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- bridge University Press, 1997)

  53. [53]

    Paulisch, H

    V. Paulisch, H. Rui, H. K. Ng, and B.-G. Englert, The European Physical Journal Plus 129, 12 (2014)

  54. [54]

    V. V. Zheleznyakov, V. V. Kocharovski˘ ı, and V. V. Kocharovski˘ ı, Soviet Physics Uspekhi26, 877 (1983)

  55. [55]

    Gardas, Journal of Mathematical Physics 51, 062103 (2010)

    B. Gardas, Journal of Mathematical Physics 51, 062103 (2010)

  56. [56]

    Bittanti, A

    S. Bittanti, A. Laub, and J. Willems, The Riccati Equa- tion, Communications and Control Engineering (Springer Berlin Heidelberg, 2012)

  57. [57]

    Schuch, Journal of Physics: Conference Series 538, 012019 (2014)

    D. Schuch, Journal of Physics: Conference Series 538, 012019 (2014)

  58. [58]

    Ndiaye, Appl

    M. Ndiaye, Appl. Math. (Irvine) 13, 774 (2022)

  59. [59]

    Abramowitz and I

    M. Abramowitz and I. Stegun, Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables (Dover Publications, 1965)

  60. [60]

    O. V. Ivakhnenko, S. N. Shevchenko, and F. Nori, Physics Reports 995, 1 (2023)

  61. [61]

    L. A. J. Guttieres, M. D. Petrovi´ c, and J. K. Freericks, American Journal of Physics 91, 885 (2023)

  62. [62]

    Virtanen, R

    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...

  63. [63]

    Bezanson, A

    J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, SIAM Review 59, 65 (2017)

  64. [64]

    Kazantsev, G

    A. Kazantsev, G. Surdutovich, and V. Yakovlev,Mechani- cal Action Of Light On Atoms (World Scientific Publishing Company, 1990)

  65. [65]

    Kayanuma, Phys

    Y. Kayanuma, Phys. Rev. A 55, R2495 (1997)

  66. [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

  67. [67]

    J. S. Briggs and A. Eisfeld, Phys. Rev. A 88, 062104 28 (2013)

  68. [68]

    Gradshteyn and I

    I. Gradshteyn and I. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 2014)

  69. [69]

    Suominen, B

    K.-A. Suominen, B. Garraway, and S. Stenholm, Optics Communications 82, 260 (1991)

  70. [70]

    Kocharovsky, V

    V. Kocharovsky, V. Kocharovsky, E. Derishev, S. Litvak, I. Shereshevsky, and S. Tasaki, Computers and Mathe- matics with Applications 34, 727 (1997)

  71. [71]

    Steck, Quantum and Atom Optics (2024)

    D. Steck, Quantum and Atom Optics (2024)