# Landau-Zener transition in a two-level system coupled to a single highly-excited oscillator

## Abstract

Two-level system strongly coupled to a single resonator mode (harmonic oscillator) is a paradigmatic model in many subfields of physics. We study theoretically the Landau-Zener transition in this model. Analytical solution for the transition probability is possible when the oscillator is highly excited, i.e. at high temperatures. Then the relative change of the excitation level of the oscillator in the course of the transition is small. The physical picture of the transition in the presence of coupling to the oscillator becomes transparent in the limiting cases of slow and fast oscillator. Slow oscillator effectively renormalizes the drive velocity. As a result, the transition probability either increases or decreases depending on the oscillator phase. The net effect is, however, the suppression of the transition probability. On the contrary, fast oscillator renormalizes the matrix element of the transition rather than the drive velocity. This renormalization makes the transition probability a non-monotonic function of the coupling amplitude.

###### pacs:

73.40.Gk, 05.40.Ca, 03.65.-w, 02.50.Ey## I Introduction

Since the publication of seminal papers (1); (2), see also the review Ref. (3), the effect of environment on the dynamics of a two-level system is modeled by introducing the coupling of the levels to the infinite set of harmonic oscillators.

In the course of the Landau-Zener (LZ) transition,(4); (5); (6); (7) when the energy levels of the two-level system undergo the avoided crossing under the action of external drive, the effect of environment amounts to the loss of adiabaticity of the transition. More quantitatively, the probability for the system to stay in the ground state after the transition is diminished by the environment.(8); (9); (10); (11); (12); (13); (14); (15); (16); (17); (18); (19) The underlying reason for this is the absorption of “quanta” of the environment. This absorption leads to decoherence, which suppresses the interference of different virtual tunneling pathways.

The situation is more delicate when the environment is represented by a single oscillator.(20); (21); (22); (23); (24); (25); (26); (27); (28) In experiment, the role of such an oscillator, which is coupled to the two-level system, is played e.g. by the transmission line resonator(29) or by the optical resonator.(30)

In the absence of coupling, the amplitude of the LZ transition can be viewed as coherent superposition of many amplitudes corresponding to virtual trajectories. With coupling, each virtual transition is accompanied by the excitation of the oscillator. On the other hand, the stronger oscillator is excited, the stronger is the feedback that it exercises on the two-level system. Then it is a compound object, two-level system dressed by many oscillator quanta, that undergoes the LZ transition.

For the two-level system coupled to the environment with a continuous spectrum only the weak-coupling regime is of interest. This is because, upon increasing coupling, the interference is completely suppressed, so that the transition probability, , assumes the value . On the contrary, when the two level system is coupled to a single oscillator, there is a wide domain of parameters when the coupling is strong while is still a strong function of the drive velocity.

Nontriviality of the Landau-Zener (LZ) transition in the presence of coupling to the oscillator is highlighted by the exact result reported in Ref. (21). This result pertains strictly to zero temperature when at time the oscillators is in the ground state. It was demonstrated in Ref. (21) that if the two-level system starts in the state and ends in the state , then the oscillator remains in the ground state at . This result can be viewed as a manifestation of the “no-go” theorem(31); (32); (33) in application to the spin-boson model. Certainly, at intermediate times, the oscillator can be excited. As a consequence of the above restriction, the two-level system and the oscillator end up entangled.

Another manifestation of nontriviality of the LZ transition with coupling to a single oscillator is the dependence of on the coupling strength, . In particular, numerical results of Ref. (25) suggest that, for finite-temperature oscillator, the dependence of on is a non-monotonic curve with a minimum. In other words, upon increasing , the adiabaticity of the transition first decreases and then increases again. There is no clear physical picture explaining the emergence of this minimum. In theory, coupling to the oscillator turns a single avoided crossing, taking place in the course of the LZ transition, into a network of avoided crossings(34) corresponding to different oscillator levels. The coupling strength quantifies the “talking” between the and amplitudes pertaining to a certain oscillator level to the corresponding amplitudes for two neighboring oscillator levels.

In general, the problem of the LZ transition in the presence of coupling to the oscillator contains, in addition to , three other parameters with the dimensionality of frequency: the matrix element between and levels, the inverse bare LZ transition time, and the oscillator frequency. Definitely, it is impossible to derive an analytical expression for for arbitrary relations between these parameters. In the present paper we focus on the situation when the oscillator is highly excited. Under this simplifying assumption we identify the domain of parameters where the asymptotic analytical expression for can be found. Roughly speaking, the two domains correspond to slow and fast oscillator depending on whether the LZ transition time is shorter or longer than the oscillator period. Coupling to a slow oscillator effectively renormalizes the drive velocity. In the case of a fast oscillator, the LZ transition splits, as a results of coupling, into a sequence of individual transitions even-spaced in time. The corresponding gaps are the oscillating functions of coupling strength, . Non-monotonic behavior of the survival probability with is the result of interference of partial transition amplitudes. We confirm this behavior by solving the many-level Schrödinger equation numerically.

## Ii Basic relations

We start from the standard Hamiltonian of the spin-boson model with drive

(1) |

where is the oscillator frequency, while and are, respectively, the annihilation and the creation operators of the oscillator. The drive is characterized by the rate, , of the change of the energies of the and states coupled directly by the matrix element . It is assumed that the coupling between the two-level system and the oscillator is longitudinal.

Denote with and the amplitudes to find the system in the states and with quanta excited. As follows from Eq. (1), these amplitudes satisfy the following infinite system of coupled equations

(2) |

(3) |

Our goal is to find the analytical solution of this system in the limit when the oscillator is highly excited, so that the relevant -values are big. In this limit, we can neglect the difference between and . Denote the initial state of the oscillator with . A crucial simplification is achieved if, in the course of the transition, the excitation level of the oscillator changes relatively weakly, i.e. by quanta with much smaller than . This allows to eliminate the explicit -dependence from the system Eq. (3). Upon introducing the new variables

(4) |

we get

(5) |

where

(6) |

Obviously, the partial solutions of the system Eq. (5) are the plane waves

(7) |

where is the wave vector. The amplitudes , satisfy the system

(8) |

This system describes the LZ transition within a given . The form Eq. (8) suggests the interpretation of as a phase of the classical oscillator.

Original system Eq. (3) describes the “spreading” of the initial state with over the states with . The survival probability is the probability for the system, which starts at from a single nonzero amplitude , to remain in one of the states at . After reducing the original system to the form Eq. (8) the amplitude represents a combination . The initial condition that at the -dependence of amplitude is suggests that the solutions of Eq. (8) corresponding to different have the same absolute value at . This allows to express the net LZ survival probability via the survival probabilities corresponding to all -values

(9) |

In the remainder of the paper we study the dependence of found from Eqs. (8), (9) in the limits of slow and fast oscillator.

## Iii Slow oscillator

To illustrate how delicate is the effect of coupling to the low-frequency oscillator on the survival probability, we plot in Figs. 1a, b the numerical solutions of the system Eq. (8) for different -values. Fig. 1a corresponds to the wave vector , while Fig. 1b corresponds to . The solutions correspond to the “subgap” frequency . The drive velocity is chosen to be , so that, without coupling, the survival probability is small, . Comparing Figs. 1a and 1 b, we conclude that the effect of coupling on is very different for these two values of . For the survival probability increases monotonically with the coupling strength, , while for , already for the minimal coupling , the value of is smaller than in the absence of coupling.

Formal explanation of this peculiar dependence of on follows from the expression of the time-dependent energy levels of the system Eq. (8). This expression reads

(10) |

In Fig. 2a we plot these levels for and . It is seen that the plots for lie on the opposite sides of the curve. Thus, for , the coupling to the oscillator effectively increases the drive velocity, and, thus, gets enhanced. For , the effective drive velocity is decreased due to the coupling to the oscillator and, correspondingly, is diminished,

Dramatic difference of the survival probabilities for different -values becomes even more dramatic upon further increase of coupling strength. This is illustrated in Fig. 3, where the curves obtained numerically are plotted for , , and . While all three curves start from , the and curves increase with , while the curve decreases with . It also follows from Fig. 3 that beyond certain -value all three curves exhibit strong oscillations.

The goal of the theory is to account for the shapes of the curves. To this end, we recall that, in the absence of coupling, the most concise way to derive expression is to perform the analytic continuation of the semiclassical solutions, , for the , amplitudes to the complex plane(35). Then emerges in the form of the following integral between the turning points on the imaginary axis

(11) |

Here and are the left and the right turning points, , which, in the absence of coupling, are simply equal to . General expression Eq. (11) suggests that the extension to the finite coupling at amounts to the modification of

(12) |

The equation for , becomes transcendental. Still, for a given set of parameters, the dependence determined by Eqs. (11) and (12) can be obtained by the numerical integration. The results are shown in Fig. 3. They agree very well with found from numerical solution of the system (8).

At weak coupling, the term amounts to the modification of the drive velocity. For the effective velocity obtained by expansion at small we find

(13) |

Upon substituting into the LZ survival probability we get the results which agree perfectly with the result obtained above using the semiclassical , Eq. (11). This agreement could be expected only in “perturbative” regime , but, for numerical reasons, this agreement holds up to the maximal value of . For this maximum value at turns to zero. Thus, we use this simplified procedure for arbitrary . In doing this, it is very important to take into account that for different from the transition point is shifted from to some . The combination should be replaced by . Fig. 3c illustrates that theoretical dependence is in good agreement with numerical solution of the system Eq. (8) at .

Summarizing, we write the expression for survival probability in the limit of a slow oscillator in the form

(14) |

where is determined by the condition

(15) |

The final step is averaging of Eq. (14) over . This averaging can be performed analytically when the renormalization of the velocity due to the coupling to the oscillator is weak. Expanding the denominator, we get

(16) |

where is the modified Bessel function. We note that, while deriving this result, we assumed that is much smaller than , the argument of can be bigger than . This is because this argument contains an additional big factor . In other words, Eq. (16) captures strong enhancement of the survival probability caused by the coupling to the oscillator.

In Fig. 3d we compare the results of three approaches to the calculation of the evolution of with coupling strength. The first result, shown with blue curve, is purely numerical. Namely, the dependence was obtained for each -value and then averaged over numerically. The second result (red line) is semi-analytical, obtained from Eqs. (14), (15), and, finally, the analytical result Eq. (16) (green line). As could be expected, Eq. (16) captures the behavior only for small couplings. The semi-analytical descriptions works well until . The origin of the discrepancy between this description and the numerics is that is the result of the averaging of rapidly growing and rapidly decaying contributions.

In the limit of large the curves in Fig. 3 start to oscillate. The oscillations survive the averaging over . The origin of these oscillations becomes clear from Fig. 2b. LZ transition at small evolves into a sequence of individual well-resolved LZ transitions upon increasing . The net number of transitions, , grows linearly with coupling. Passage of these transitions may result in constructive or destructive interference depending on the phase accumulated between the subsequent transitions. The situation is fully analogous to the Landau-Zener-Stückelberg interferometry.(36)

## Iv Fast oscillator

From Fig. 2c we realize that the description based on time-dependent energy levels is inadequate for the fast oscillator. This is also clear from physical arguments, since the “local” velocity is much bigger than the drive velocity. The role of the oscillator at large is to renormalize not the drive velocity but rather the matrix element, , between the levels. To see this, we make the following substitution in the system Eq. (8)

(17) |

after which it acquires the form

(18) |

In the form Eq. (IV), the inter-level matrix element oscillates with time. If the time of the LZ transition, , is much longer than the period of oscillations, the gap oscillates many times in the course of the transition. This suggests that the oscillating factor can be replaced(37) by its average , where is the Bessel function. This replacement immediately leads to the survival probability

(19) |

This result suggests that a very small at zero coupling increases rapidly with coupling, reaches , when the Bessel function passes through zero, and then drops down.

To check numerically the validity of averaging over the oscillator period, in Fig. 4 we show the time dependence calculated by numerical solution of Eq. (IV) for particular value , when the Bessel function turns to zero. We see that, for this coupling, does not change with time in a certain domain around suggesting that the effective gap is zero. However, unlike what Eq. (19) predicts, the value of is not one in this domain. The reason is that the system Eq. (IV) encodes a number of individual transitions which take place around the times . This becomes apparent if we make the following substitution in the system Eq. (IV)

(20) |

Upon this substitution, Eq. (IV) assumes the form

(21) |

(22) |

The same argument as above suggests that, as a result of being fast, the exponent in the right-hand side can be averaged over the period, , Then the system Eq. (21) will describe a regular LZ transition taking place around with a gap reduced by , where is the first-order Bessel function. The corresponding survival probability reads

(23) |

Naturally, a similar transition taking place around

is described
by the same . Note also, that, in addition to ,
the averaged matrix element is multiplied by .
The physical meaning of the moments is
transparent. The energy separation between and state
changes with time as . At this separation becomes
equal to the oscillator quantum.

The extension of Eqs. (19), (23) to arbitrary is straightforward:

(24) |

The transitions taking place at can be viewed as well separated if the time, , is much bigger than the individual LZ transition time , yielding the criterion . The second condition to be met is that the effective time averaging takes place during the LZ transition time, so that . The second condition can be cast in the form . If the bare survival probability is small, , then the first condition is more restrictive.

Fig. 4 offers an insight into a general scenario of the LZ transition in a two-level system coupled to a fast oscillator. The system undergoes a number of individual transitions at times characterized by survival probabilities . It is also seen from Fig. 4 that the evolution of with time depends on the wave vector, . This is a consequence of interference of partial transition amplitudes. To find the net survival probability analytically, we first assume that averaging over suppresses the interference effects completely. Then we can write the recurrent relation for , which are the successive values of after transitions

(25) |

This relation expresses the fact that the occupation of the state after the transition comes from the occupation of this state before the transition as well as from the occupation of the state which flips in the course of the transition. It is straightforward to derive from Eq. (25) the expression for after the arbitrary number of transitions

(26) |

In Fig. 5 the behavior of versus the coupling is plotted from Eq. (26) for different bare values of . Naturally, all the curves approach at very large coupling. This is because, at large coupling, the information about the initial state of the two-level system is erased. For the bare value the approach to is monotonic. However, for a wiggle emerges at . This feature evolves into a well pronounced maximum for bare .

At the position of maximum in Fig. 5 (green line) the argument, , of the Bessel functions is about . For this value and also for smaller couplings we can restrict consideration to only three LZ transitions taking place at and , since all higher Bessel functions are small. With only three transitions, we can incorporate the interference effects into the theory. To do so, we should take into account that each partial LZ transition is characterized by a scattering matrix which contains the survival probability and a phase, . The evolution of the amplitudes to find the system in and states is described by the product of scattering matrices

(27) |

From this product one can infer the following expression for the amplitude to change the level after the three transitions(38)

(28) |

Different contributions to describe partial amplitudes to change the level at one transition and to not change the level at other two transitions. The survival probability is given by .

If we assume that the phases , , and are completely uncorrelated the averaging over these phases yields

(29) |

which is nothing but the result Eq. (26) in which only the terms and are kept. Now we take into account that the transitions and are identical, set , and perform the averaging over the two phases. This gives

(30) |

We note that the interference contribution to Eq. (30) is negative. In fact, it leads to a maximum in behavior even for the bare , as illustrated in Fig. 6. Thus, it is the result Eq. (30) that should be compared to the numerical calculations. The results of these calculations are shown in Fig. 7, where probability, , calculated by solving the system Eq. (IV) and averaging over is shown for three oscillator frequencies versus the dimensionless coupling amplitude. In the domain all three curves coincide and agree with the “incoherent” and “coherent” theoretical results shown in Fig. 6. The position of the maxima also agrees with the prediction of the “coherent” theory Eq. (30). However, the peak value is higher than predicted by the theory. The possible origin of the discrepancy is that the theory Eq. (30) takes into account only and intermediate transitions. Our overall conclusion is that non-monotonic behavior of is the result of the interference of intermediate LZ transitions.

## V Discussion and concluding remarks

In the present paper we focused on the question: how the longitudinal coupling to a harmonic oscillator affects the survival probability of the Landau-Zener transition in a driven two-level system.

On general grounds, one would expect the following answer to this question. Weak coupling, by making the transition less adiabatic, increases the survival probability. At very strong coupling this probability should approach , since the memory about the initial state of the two-level system gets erased due to coupling. There is, however, the evidence that these expectations are not entirely correct. Firstly, the exact result obtained in Ref. (21) states that, for purely longitudinal coupling, the effect is identically zero for any coupling strength. Secondly, the numerical simulations of Ref. (25), which pertain to longitudinal coupling suggest that, at finite temperature, is a non-monotonic function with a maximum. The domain of parameters investigated in Ref. (25) is intermediate in all respects: the temperature, the oscillator frequency, and the LZ tunneling gap were of the same order. In this regime it is difficult to infer the underlying origin of this maximum.

To establish the above physical picture, we have adopted a strong assumption that the oscillator is highly excited, so that the change of the excitation level in the course of the transition is relatively small. Under this assumption, we studied the effect of coupling on the LZ transition in two limiting cases of slow and fast oscillator. For the slow oscillator our analytical results for the survival probability are given by Eqs. (14), (15), and Eq. (16). For the fast oscillator they are given by Eq. (26)and Eq. (30). We can now quantify the validity of our main assumption. For both, the slow and the fast oscillator, the LZ transition in the presence of coupling transforms into a sequence of individual transitions.(39) This is illustrated in Fig. 2b for the slow oscillator and in Fig. 4 for the fast oscillator. The number of transitions for the slow oscillator is . For the fast oscillator this number can be estimated by equating the argument of the Bessel function in Eq. (24) to the index, i.e. . Since each individual transition is associated with different level of the oscillator, the criterion that the oscillator is highly excited can be quantified as follows: the initial excitation level should be bigger than for the slow oscillator, while for the fast oscillator it should be bigger than .

Although our analytical results essentially confirm the general expectations, we find that, in both limits, approaches with oscillations. These oscillations are the consequence of the interference of the amplitudes corresponding to different pathways through multiple LZ transitions. It is likely that non-monotonic established in Ref. (25) is the consequence of this interference.

Note finally, that throughout the paper we assumed the bare survival probability is small, so that, unlike Ref. (40), the perturbative treatment does not apply.

Acknowledgements

The work was supported by the Department of Energy, Office of Basic Energy Sciences, Grant No. DE- FG02-06ER46313.

### References

- A. O. Caldeira and A. J. Leggett, “Influence of Dissipation on Quantum Tunneling in Macroscopic Systems,” Phys. Rev. Lett. 46, 211 (1981).
- S. Chakravarty and A. J. Leggett, “Dynamics of the Two-State System with Ohmic Dissipation,” Phys. Rev. Lett. 52, 5 (1984).
- A. Leggett, S. Chakravarty, A. Dorsey, M. Fisher, A. Garg, and W. Zwerger, “Dynamics of the dissipative two-state system,” Rev. Mod. Phys. 59, 1 (1987).
- L. D. Landau, “Zur theorie der energieubertragung,” Physics of the Soviet Union 2, 46 (1932).
- C. Zener, “Non-adiabatic crossing of energy levels,” Proc. R. Soc. London A 137, 696 (1932).
- E. Majorana, “Atomi orientati in campo magnetico vari- abile, Nuovo Cimento 9, 43 (1932).
- E. C. G. Stueckelberg, “Theorie der unelastischen St ¨osse zwischen atomen, Helv. Phys. Acta 5, 369 (1932).
- Y. Gefen, E. Ben-Jacob, and A. O. Caldeira, “Zener transitions in dissipative driven systems,” Phys. Rev. B 36, 2770 (1987).
- P. Ao and J. Rammer, “Quantum dynamics of a two-state system in a dissipative environment,” Phys. Rev. B 43, 5397 (1991).
- E. Shimshoni and A. Stern, “Dephasing of interference in Landau-Zener transitions,” Phys. Rev. B 47, 9523 (1993).
- P. Nalbach and M. Thorwart, “Landau-Zener Transitions in a Dissipative Environment: Numerically Exact Results,” Phys. Rev. Lett. 103, 220401 (2009).
- S. Gasparinetti, P. Solinas, and J. P. Pekola, “Geometric Landau-Zener Interferometry,” Phys. Rev. Lett. 107, 207002 (2011).
- R. S. Whitney, M. Clusel, and T. Ziman, “Temperature Can Enhance Coherent Oscillations at a Landau-Zener Transition,” Phys. Rev. Lett. 107, 210402 (2011).
- C. Xu, A. Poudel, and M. G. Vavilov, “Nonadiabatic dynamics of a slowly driven dissipative two-level system,” Phys. Rev. A 89, 052102 (2014).
- P. Nalbach, “Adiabatic-Markovian bath dynamics at avoided crossings,” Phys. Rev. A 90, 042112 (2014).
- S. Javanbakht, P. Nalbach, and M. Thorwart, “Dissipative Landau-Zener quantum dynamics with transversal and longitudinal noise,” Phys. Rev. A 91, 052103 (2015).
- M. Yamaguchi, T. Yuge, and T. Ogawa, “Markovian quantum master equation beyond adiabatic regime,” Phys. Rev. E 95, 012136 (2017).
- L. Arceci, S. Barbarino, R. Fazio, G. E. Santoro, “Dissipative Landau-Zener problem and thermally assisted Quantum Annealing,” Phys. Rev. B 96, 054301 (2017).
- R. K. Malla, E. G. Mishchenko, and M. E. Raikh, “Suppression of the Landau-Zener transition probability by a weak classical noise,” Phys. Rev. B 96, 075419 (2017)
- M. Wubs, K. Saito, S. Kohler, Y. Kayanuma, and P. Hänggi, “Landau-Zener transitions in qubits controlled by electromagnetic fields,” New J. Phys. 7, 218 (2005).
- M. Wubs, K. Saito, S. Kohler, P. Hänggi, and Y. Kayanuma, “Gauging a Quantum Heat Bath with Dissipative Landau-Zener Transitions,” Phys. Rev. Lett. 97, 200404 (2006).
- K. Saito, M. Wubs, S. Kohler, P. Hänggi, and Y. Kayanuma, “Quantum state preparation in circuit QED via Landau-Zener tunneling,” Europhys. Lett. 77, 22 (2006).
- K. Saito, M. Wubs, S. Kohler, Y. Kayanuma, and P. Hänggi, “Dissipative Landau-Zener transitions of a qubit: Bath-specific and universal behavior,” Phys. Rev. B 75, 214308 (2007).
- D. Zueco, P. Hänggi, and S. Kohler, “LandauZener tunnelling in dissipative circuit QED,” New J. Phys. 10, 115012 (2008).
- S. Ashhab, “Landau-Zener transitions in a two-level system coupled to a finite-temperature harmonic oscillator,” Phys. Rev. A 90, 062120 (2014).
- S. Ashhab, “Landau-Zener transitions in an open multilevel quantum system,” Phys. Rev. A 94, 042109 (2016).
- N. A. Sinitsyn and F. Li, “Solvable multistate model of Landau-Zener transitions in cavity QED,” Phys. Rev. A 93, 063859 (2016).
- C. Sun and N. A. Sinitsyn, “Landau-Zener extension of the Tavis-Cummings model: Structure of the solution,” Phys. Rev. A 94, 033808 (2016).
- O. Astafiev, K. Inomata, A. O. Niskanen, T. Yamamoto, Yu. A. Pashkin, Y. Nakamura, and J. S. Tsai, “Single artificial-atom lasing,” Nature (London) 449, 588 (2007).
- A. D. Cimmarusti, Z. Yan, B. D. Patterson, L. P. Corcos, L. A. Orozco, and S. Deffner, “Environment-Assisted Speed-up of the Field Evolution in Cavity Quantum Electrodynamics,” Phys. Rev. Lett. 114, 233602 (2015).
- N. A. Sinitsyn, “Counterintuitive transitions in the multistate Landau-Zener problem with linear level crossings,” J. Phys. A 37, 10691 (2004).
- A. V. Shytov, “Landau-Zener transitions in a multilevel system: An exact result,” Phys. Rev. A 70, 052708 (2004).
- A. A. Rangelov, J. Piilo, and N. V. Vitanov, “Counterintuitive transitions between crossing energy levels,” Phys. Rev. A 72, 053404 (2005).
- D. A. Garanin, R. Neb, and R. Schilling, “Effect of environmental spins on Landau-Zener transitions,” hys. Rev. B 78, 094405 (2008).
- A. M. Dykhne, “Adiabatic perturbation of discrete spectrum states,” Sov. Phys. JETP 14, 941 (1962).
- S. N. Shevchenko, S. Ashhab, and F. Nori, “Landau-Zener-Stückelberg interferometry,” Phys. Rep. 492, 1 (2010).
- R. Glenn, M. E. Limes, B. Pankovich, B. Saam, and M. E. Raikh, “Magnetic resonance in slowly modulated longitudinal field: Modified shape of the Rabi oscillations,” Phys. Rev. B 87, 165205 (2013).
- R. C. Roundy and M. E. Raikh, “Tunnel magnetoresistance in organic spin valves in the regime of multistep tunneling,” Phys. Rev. B 88, 205206 (2013).
- M. F. O’Keeffe, E. M. Chudnovsky, and D. A. Garanin, “Landau-Zener dynamics of a nanoresonator containing a tunneling spin,” Phys. Rev. B 87, 174418 (2013).
- D. A. Garanin and R. Schilling, “Many-body Landau-Zener effect at fast sweep,” Phys. Rev. B 71, 184414 (2005).