Next Article in Journal
Casimir Interaction of Chern–Simons Layers on Substrates via Vacuum Stress Tensor
Next Article in Special Issue
Ultrafast Resonant Photon Emission from a Molecule Driven by a Strong Coherent Field in Terms of Complex Spectral Analysis
Previous Article in Journal
A Theory of Best Choice Selection through Objective Arguments Grounded in Linear Response Theory Concepts
Previous Article in Special Issue
Electric Octupole-Dependent Contributions to Optical Binding Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation

by
Benedetto Militello
1,2 and
Anna Napoli
1,2,*
1
Dipartimento di Fisica e Chimica—Emilio Segrè, Università degli Studi di Palermo, Via Archirafi 36, 90123 Palermo, Italy
2
Istituto Nazionale di Fisica Nucleare, Sezione di Catania, Via Santa Sofia 64, 95123 Catania, Italy
*
Author to whom correspondence should be addressed.
Physics 2024, 6(2), 483-495; https://doi.org/10.3390/physics6020032
Submission received: 22 December 2023 / Revised: 31 January 2024 / Accepted: 5 February 2024 / Published: 28 March 2024

Abstract

:
The Stimulated Raman Adiabatic Passage (STIRAP) on a three-state system interacting with a spin bath is considered, focusing on the efficiency of the population transfer. Our analysis is based on the perturbation treatment of the interaction term evaluated beyond the Rotating Wave Approximation, thus focusing on the limit of weak system–bath coupling. The analytical expression of the correction to the efficiency and the consequent numerical analysis show that, in most of the cases, the effects of the environment are negligible, confirming the robustness of the population transfer.

1. Introduction

A quantum system ruled by a slowly varying Hamiltonian undergoes a dynamics known as adiabatic following of the eigenstates, based on the adiabatic theorem, according to which the populations of the instantaneous eigenstates of the Hamiltonian do not change [1,2] (details are given in Appendix A). This physical behavior is the key ingredient of many protocols aimed at controlling a quantum system [3,4,5,6,7]. Stimulated Raman Adiabatic Passage (STIRAP) [8,9,10,11,12] represents an essential and known example of an adiabatic process. The STIRAP technique was introduced to realize a complete population transfer from one state to another one by exploiting a Raman scheme involving suitable pulses with time-dependent amplitudes which couple each of the two previously mentioned states with an auxiliary one. The pulses have to be set in such a way that they should ensure the validity of the adiabatic approximation. Moreover, their specific time-dependence must be such that an eigenstate of the Hamiltonian coincides with the initial state of the system at the beginning of the process and with the target state at the end of the application of the pulses. Therefore, differently to what one might think at first glance, the total population transfer from the initial state to the target state is not due to radiative two-photon processes but is plainly due to an adiabatic following of an Hamiltonian eigenstate. This is particularly evident when the so called counter-intuitive sequence is considered, where the pulse which couples the auxiliary state and the target one precedes the pulse which couples the initial state to the auxiliary one. Indeed, in such a case it is not reasonable to interpret the whole process as a (possibly virtual) photon absorption concomitant to the initial auxiliary state transition followed by a (possibly virtual) photon emission concomitant to the auxiliary target state transition. Moreover, it is worth observing that, not only does the counter-intuitive sequence work well, but it usually works better than the so called intuitive sequence where the coupling between the auxiliary state and the initial one precedes the other coupling. A detailed analysis of the counter-intuitive sequence is given in Section 2.2.
The STIRAP technique is still extensively investigated [13,14,15,16,17,18] and has been exploited in different physical contexts ranging from cold gases [19,20,21] to condensed matter [22,23,24,25,26,27,28], plasmonic systems [29,30], superconducting devices [31,32,33], trapped ions [34,35] and optomechanical systems [36]. Recently, in order to improve the original technique by shortening the population transfer process, modifications to the original scheme including shortcuts to adiabaticity have been proposed [37,38,39,40,41,42]. However, in this case, a more complicated apparatus is required, which constitutes a disadvantage with respect to the original scheme. Indeed, the implementation of shortcuts requires the activation of additional interactions which somehow compensate for any deviation from the adiabatic following of the eigenstates of the Hamiltonian, even in cases where the Hamiltonian does not change slowly. In this way, the original time-dependent Hamiltonian of the system, H ( t ) , does not rule the dynamics anymore, since new terms, δ H ( t ) , are also considered, and the total Hamiltonian, H ( t ) + δ H ( t ) , induces a unitary evolution which coincides with an adiabatic following of the eigenstates of H ( t ) only. The structure of δ H ( t ) is generally complicated and the relevant pulses need to be very precise.
As a general feature, a quantum system is subjected to the effects of noise either rising from the interaction with external systems, such as, for example, the constituents of the environment, or related to imperfections in the apparatus. Therefore, on the one hand, the fidelity of the population transfer with respect to uncertainty or fluctuations in the amplitudes and phases of the pulses has been analyzed [43,44]. On the other hand, the interaction with the quantized electromagnetic field has also been considered, for example, by exploiting effective non-Hermitian Hamiltonians, which is limited to cases where the states involved in the STIRAP scheme decay towards states not involved in the procedure [45]. Beyond such a specific scenario, a more general approach based on the theory of open quantum systems is required [46,47]. Along this line, master equations have been proposed to study STIRAP-manipulated systems [48,49,50,51]. Moreover, under suitable assumptions fitting the Davies–Spohn theory for open quantum systems ruled by time-dependent Hamiltonians [52], time-dependent master equations in the Lindblad form can be obtained from a microscopic model of the interactions between the atomic system and the quantized field [53,54].
Since in many cases the system to be manipulated is close to other atomic systems and interacts with them, it might be the case that the interaction with the electromagnetic field is not the main source of quantum noise or, at least, not the only significant source. Nitrogen vacancies in diamonds [24,27,28] or rare-earth doped crystals [22,23] are two typical examples of such a scenario. The manipulation of spin defects in magnetic materials through adiabatic following-based techniques has recently been studied in the presence of their interaction with the surrounding spins [55]. Moreover, very recently, the effect of STIRAP processes on a system interacting with a spin bath has been theoretically analyzed under the assumptions that allow for the Rotating Wave Approximation (RWA) in the interaction between the three-state system and the spins in the environment [56]. This study has been developed through an evaluation of the unitary dynamics, despite the fact that master equations can also be derived for the spin environments [57,58,59].
In this paper, we extend the previous study in ref. [56]; we still exploit the unitary evolution of the universe but overcome the RWA in the system–environment interaction, which makes the physical model more realistic. In fact, while the RWA implies the conservation of the total number of excitations (a feature which has been extensively used in the previous analysis), it somehow excludes a variety of possible transitions. On the contrary, in this study, we take into account all the terms of the system–bath interaction, thus introducing processes which can be responsible for a reduction in the efficiency of population transfer. The paper is structured as follows. In Section 2, we describe the physical system and the relevant Hamiltonian model, and also provide a brief sketch of the STIRAP technique in the ideal case and beyond, also introducing the theoretical analysis based on the perturbation approach, specifically for the truncation of the Dyson series after two changes to picture. In Section 3, we give the explicit form of the correction to the efficiency of the population transfer according to our theory and then we show predictions based on numerical calculations. Finally, in Section 4, we present an extensive discussion on the results. Two Appendixes complete the presentation: in the Appendix A, the adiabatic approximation is recalled while in the Appendix B, we provide details about the matrix elements involved in the perturbation treatment.

2. Physical Model and Methods

2.1. Hamiltonian Model

The physical system we are focusing on consists of a three-state system subjected to two coherent fields and interacting with the surrounding environment, consisting of a spin bath. The three-state system has two ground states, g 1 and g 2 , and an excited state, e , and its free dynamics are governed by the Hamiltonian H A as given just below. The action of the STIRAP pulses coupling each of the two ground states with the excited one (see Figure 1) is described by H S . In addition, the free spin-bath (an ensemble of two-state systems) is described by H B while the system–bath interaction is described by H A B , which associates spin flips with atomic transitions between the excited state and each of the two ground states. Therefore, the total Hamiltonian is given by ( = 1 ):
H = H A + H S + H B + H A B ,
with:
H A = ν e e
H S = m = 1 2 Ω m ( t ) cos ( ν t ) ( g m e + e g m ) ,
H B = k = 1 L ω 2 σ z ( k ) ,
H A B = m = 1 2 g m e + e g m k = 1 L η k ( m ) σ x ( k ) ,
where ν is the energy gap between the free excited states and the two grounds and ν is the frequency of the two pulses, whereas the Ω m ’s describe the profiles of the pulses; the natural frequency of the spins of the bath is ω , while the quantity η k ( m ) is the spin–atom coupling constant related to the k-th spin and to the atomic transitions e g m . Finally σ α ( k ) are the Pauli operators associated to the k-th spin.
Introducing the Hermitian operator, G 1 = ν e e + H B , and the relevant unitary operator U 1 , such that i U ˙ 1 = G 1 U 1 , in the new picture given by | ψ ˜ ( t ) = U 1 ( t ) ψ ( t ) and A ˜ ( t ) = U 1 ( t ) A ( t ) U 1 ( t ) one finds that the generator of the time evolution is:
H ˜ ( t ) G 1 = H A + H ˜ S ( 0 ) ( t ) + H ˜ S ( ω ) ( t ) + H ˜ A B ( t ) ,
with
H A = ( ν ν ) e e ,
H ˜ S ( 0 ) ( t ) = m = 1 2 Ω m ( t ) ( g m e + e g m ) ,
H ˜ S ( ω ) ( t ) = m = 1 2 Ω m ( t ) g m e e i 2 ν t + h . c . ,
H ˜ A B ( t ) = m = 1 2 ( g m e e i ν t + e g m e i ν t ) k η k ( m ) ( σ ( k ) e i ω t + σ + ( k ) e i ω t ) ,
with ⊗ denoting the tensor product, “h.c.” standing for the “Hermitian conjugate” term, Δ = ν ν being the detuning between the atomic frequency and the field frequency. It is worth mentioning that, in view of the further treatments, we have split the transformed pulse Hamiltonian H ˜ S into two contributions, H ˜ S = H ˜ S ( 0 ) + H ˜ S ( ω ) , where the first corresponds to the so called rotating terms (characterized by the absence of fast oscillations) while the second is related to the counter-rotating terms (rapidly oscillating). For our treatment of the system–bath interaction terms in H ˜ A B ( t ) , this separation of rotating and counter-rotating terms is not necessary.

2.2. Ideal STIRAP

Let us first sum up the basic form of the ideal STIRAP, which corresponds to the absence of interaction with the environment (referring to our model, this condition is accomplished assuming η k ( m ) = 0 , m , k ). Moreover, assuming high atomic and field frequencies, one is legitimated to neglect the counter-rotating terms in the STIRAP Hamiltonian, so that the system can be assumed to be approximately ruled only by H A + H ˜ S ( 0 ) in this new picture usually, which is referred to as the rotating frame. The operator
H A + H ˜ S ( 0 ) ( t ) = Δ e e + m = 1 2 Ω m ( t ) ( g m e + e g m ) ,
can be diagonalized at every instant of time, and its instantaneous eigenstates are:
+ ( t ) = sin φ ( t ) sin θ ( t ) g 1 + cos φ ( t ) e + sin φ ( t ) cos θ ( t ) g 2 ,
0 ( t ) = cos θ ( t ) g 1 sin θ ( t ) g 2 ,
( t ) = cos φ ( t ) sin θ ( t ) g 1 sin φ ( t ) e + cos φ ( t ) cos θ ( t ) g 2 ,
where
tan θ ( t ) = Ω 1 ( t ) Ω 2 ( t ) ,
tan 2 φ ( t ) = 2 Ω L ( t ) Δ ,
Ω L ( t ) = Ω 1 2 ( t ) + Ω 2 2 ( t ) .
Such eigenstates correspond to the following eigenvalues:
E 0 = 0 , E ± = 1 2 Δ ± Δ 2 + 4 Ω L 2 ( t )
The standard STIRAP process aimed at transferring population from the state g 1 to g 2 is realized through a so-called counterintuitive sequence, where the pulse Ω 2 precedes Ω 1 . Accordingly, at the initial time θ = 0 and φ = 0 , so that + = e , 0 = g 1 and = g 2 , while in the final time θ = π / 2 and φ = 0 , so that + = e , 0 = g 2 and = g 1 . When the pulse profiles are slowly varying functions, the hypotheses of the adiabatic theorem are satisfied and the population of each eigenstate is preserved during the evolution. Consequently, in particular, the population of the state g 1 , initially coinciding with 0 , is totally transferred to the state g 2 , which equals 0 in the final time.

2.3. Time Evolution and Efficiency

In order to better evaluate the effects of the remaining terms of the Hamiltonian, we perform a new change of picture, applying the transformation ψ ¯ = U 2 ψ ˜ with U 2 such that i U ˙ 2 = ( H A + H ˜ S ( 0 ) ( t ) ) U 2 .
The operator U 2 is well approximated by the unitary evolution describing the adiabatic following of the instantaneous eigenstates of the Hamiltonian H A + H ˜ S ( 0 ) ( t ) , given by U 2 ( t , t 0 ) m e i α m ( t , t 0 ) ϕ m ( t ) ϕ m ( t 0 ) , with ϕ m ( t ) being the instantaneous eigenstates of H A + H ˜ S ( 0 ) ( t ) , and α m ( t , t 0 ) = t 0 t E m ( w ) d w +   i t 0 t ( w ϕ m ( w ) ) ϕ m ( w ) d w being the relevant phase factors. It is worth mentioning that all the geometric phases are zero, which comes from the feature that ( w ϕ m ( w ) ) ϕ m ( w ) , with w / w , is always an imaginary number, while the coefficients of the eigenstates are all real, which implies that all such terms are zero. Applying the transformation, one obtains:
H ¯ ( t ) G ¯ 2 ( t ) = H ¯ S ( ω ) ( t ) + H ¯ A B ( t ) ,
with
H ¯ S ( ω ) ( t ) = m = 1 2 Ω m ( t ) g ˚ m ( t ) e ˚ ( t ) e i 2 ν t + h . c . ,
H ¯ A B ( t ) = m = 1 2 ( g ˚ m ( t ) e ˚ ( t ) e i ν t + e ˚ ( t ) g ˚ m ( t ) e i ν t ) k η k ( m ) ( σ ( k ) e i ω t + σ + ( k ) e i ω t ) ,
with g ˚ m ( t ) = U 2 ( t , t 0 ) g m and e ˚ ( t ) = U 2 ( t , t 0 ) e .
In the new picture, the generator of the time evolution is the following:
H ¯ G = H ¯ S ( ω ) ( t ) + H ¯ A B ( t ) ,
and the relevant approximated dynamics can be evaluated, to the second order, by truncation of the iterated formal solution:
T ( t , t 0 ) 1 1 i t 0 t H ¯ G ( w ) d w t 0 t d w t 0 w d w H ¯ G ( w ) H ¯ G ( w ) T 2 ( t , t 0 ) ,
which is essentially the truncation of the Dyson series (expressed without the chronological ordering operator) to the second order.
Moreover, since we have moved to this new picture by removing the adiabatic evolution operator responsible for a perfect population transfer from g 1 to g 2 , remaining in the state g 1 is equivalent to undergoing a perfect transition from g 1 to g 2 in the Schrödinger picture. Therefore, in the new picture the efficiency of the population transfer process through STIRAP pulses is given by the survival probability of the initial state of the three-state system:
P ( t ) = tr [ T ( t , t 0 ) ψ ( 0 ) ψ ( 0 ) ρ B ( 0 ) T ( t ) ψ ( 0 ) ψ ( 0 ) 1 1 B ] ,
where ρ B ( 0 ) is the density operator describing the initial configuration of the bath, 1 1 B is the identity operator of the bath and T ( t , t 0 ) is possibly replaced by T 2 ( t , t 0 ) and “tr” stands for the trace operation. In the case considered here, ψ ( 0 ) = g 1 . In order to better understand Equation (24), consider that the complete time evolution of the initial state g 1 g 1 ρ B ( 0 ) in the Schrödinger picture is given by: ρ A B ( t ) = U 1 ( t , t 0 ) U 2 ( t , t 0 ) T ( t , t 0 ) g 1 g 1 ρ B ( 0 ) T ( t , t 0 ) U 2 ( t , t 0 ) U 1 ( t , t 0 ) . Now, since the target state is g 2 , irrespective of the state of the bath, we need to evaluate P ( t ) = tr [ ρ A B ( t ) g 2 g 2 1 1 B ] , which, after performing two cyclic permutations inside the trace functional leads to the following expression: tr [ T ( t , t 0 ) g 1 g 1 ρ B ( 0 ) T ( t , t 0 ) U 2 ( t , t 0 ) U 1 ( t , t 0 ) g 2 g 2 1 1 B U 2 ( t , t 0 ) U 1 ( t , t 0 ) ] , which is equivalent to Equation (24) once it is considered that U 2 ( t , t 0 ) g 1 = g 2 , and conversely, U 2 ( t , t 0 ) g 2 = g 1 , so that U 2 ( t , t 0 ) U 1 ( t , t 0 ) g 2 g 2 1 1 B U 1 ( t , t 0 ) U 2 ( t , t 0 ) = g 1 g 1 1 1 B .

3. Results

We now focus on the zero-temperature bath, which means assuming that the spin bath is initially in its ground state: { } = k k . Therefore, the complete initial state is g 1 { } . Since we are considering the approximation T ( t , t 0 ) T 2 ( t , t 0 ) according to Equation (23), the only transitions considered in our calculations are those involving zero, one or two spin flips in the bath. This implies the following form for the efficiency:
P ( t ) | { } g 1 T 2 ( t , t 0 ) g 1 { } | 2 + Σ l | { } l { } g 1 T 2 ( t , t 0 ) g 1 { } | 2 + Σ j l | { } j { } l { } g 1 T 2 ( t , t 0 ) g 1 { } | 2 ,
where { } l { } is the bath state with all spins in the state, except for the l-th spin, which is the state, while { } j { } l { } has only the l-th and j-th spins in the state, all the others being in the state . The overlaps involving only one spin flip turn out to be zero (see Appendix B for details). The overlaps involving two spin flips involve only second-order terms and, once their squared modulus is evaluated, such terms give rise to fourth-order contributions. To make the calculation consistent with the truncation of the Dyson series to the second order, only terms up to the second order are to be kept in the probability, which gives the following expression:
P ( t ) 1 2 J η k ( m ) ( t ) ,
where denotes the real part, and its argument is the following integral:
J η k ( m ) ( t ) = t 0 t d w t 0 w d w [ e i w w E + ( s ) d s cos φ ( w ) cos φ ( w ) + e i w w E ( s ) d s sin φ ( w ) sin φ ( w ) ] e i ( ν + ω ) ( w w ) × k η k ( 1 ) cos θ ( w ) η k ( 2 ) sin θ ( w ) η k ( 1 ) cos θ ( w ) η k ( 2 ) sin θ ( w ) .
We now assume that the quantities η k ( m ) , first defined after Equation (5) without any constraint, are related in such a way that the ratios η k ( m ) / η 1 ( m ) do not depend on m, so that one can introduce λ k η k ( m ) / η 1 ( m ) , define the following quantities:
η η 1 ( 1 ) , r η η 1 ( 2 ) / η 1 ( 1 ) , Λ k λ k 2 ,
and recast the probability in the following form:
P ( t ) 1 2 η 2 Λ [ J ( t ) ] ,
with
J ( t ) = t 0 t d w t 0 w d w e i w w E + ( s ) d s cos φ ( w ) cos φ ( w ) + e i w w E ( s ) d s sin φ ( w ) sin φ ( w ) × e i ( ν + ω ) ( w w ) cos θ ( w ) r η sin θ ( w ) cos θ ( w ) r η sin θ ( w ) .
On the basis of Equation (29), the survival probability of the initial state in the interaction picture, which corresponds to the efficiency of the population transfer in the Schrödinger picture, turns out to differ from unity by a term proportional to the real part of the the integral J ( t ) (having the dimensions of the square of time), on which we focus in what follows. Every specific correction should take into account the specific value of the square of the quantity η Λ (having the dimension of a frequency), which somehow is a cumulative measure of the coupling strength between the three-level system and the whole environment. It is also interesting to observe that, in the RWA, the counterpart of J ( t ) would be zero. This can be straightforwardly proven by recalculating the relevant matrix elements. Moreover, it is already clearly visible from the expression of J ( t ) , where all the terms contain rapidly oscillating factors coming from the fact that all the contributions come out as matrix elements of the counter-rotating terms.
It is worth mentioning that the assumption η k ( m ) / η 1 ( m ) , independent from the index m, is not as restrictive as one could think, since the coupling strengths are supposed to be proportional to a function of the spin distance from the central three-state system, and this proportionality function is supposed to be the same independently of the specific transitions involved, whether e g 1 or e g 2 . Concerning the shape of the pulses, they are usually taken as Gaussian, with the peaks occurring at different times. In particular, we assume two pulses centered at ± τ and with a width τ / 2 :
Ω 1 = Ω 0 e ( t / τ + 1 ) 2 ,
Ω 1 = Ω 0 e ( t / τ 1 ) 2 .
The process is supposed to start at time t 0 = T (with T > τ ) and finish at t = T . In Figure 2 and Figure 3 the numerically calculated quantity 2 { J ( T ) } is reported as a function of different parameters, which allows us to evaluate the efficiency of the population transfer in several regimes: the smaller the quantity 2 { J ( T ) } , the more efficient the population transfer, according to Equation (29). In all the plots we have assumed that the parameters of the pulses satisfy the following conditions, which guarantee an optimal transfer in the ideal case: Ω 0 τ = 10 , T / τ = 5 , ν τ = 10 .
In Figure 2, 2 { J ( T ) } is shown as a function of ω and Δ for two values of the ratio r η , in particular, for r η = 1 (Figure 2a) and r η = 1 (Figure 2b). In both plots it is clearly visible that the value of 2 { J ( T ) } is always small, never exceeding the value 0.03 , and that, for high values of the frequency ω , the corrections become smaller and smaller. For any fixed value of Δ , one can see that, when varying the value of ω , the quantity 2 { J ( T ) } exhibits an oscillatory behavior, which a posteriori can be related to the presence of the phase factor exp [ ( ω + ν ) ( w w ) ] in the integrand. Actually, in spite of the presence of other functions, the mentioned phase factor is indeed the main rapidly changing factor, whereas trigonometric functions of ϕ and θ are smoothly changing and the phase factor associated with the integral of E + changes rapidly only in the region between the two peaks. Thus, both the presence of oscillations and the vanishing of J ( T ) for higher values of ν are traceable back to the oscillatory character of the counter-rotating terms. By comparing Figure 2a and Figure 2b, it emerges that, whether the value of the ratio parameter r η is 1 or 1 , the behavior is quite similar, though significant differences are present, particularly in the region corresponding to small values of Δ .
In Figure 3, 2 { J ( T ) } is shown as a function of ω and r η for two values of the detuning Δ , in particular, for Δ / ν = 0.01 (Figure 3a) and Δ / ν = 0.05 (Figure 3b). One can see that the values of the quantity in a particular parameter region are much higher than in Figure 2, reaching the value of 0.8 . These high values can be justified by the fact that increasing the value of r η implies increasing the amplitude of the coupling term between the three-state system and the environment. In particular, a high value of r η means dealing with a high value of η 1 ( 2 ) (and consequently all η k ( 2 ) ), which in turn implies a higher value of J ( T ) .
In order to fix the ideas about the meaning of our results, let us consider the case where η Λ τ = 1 means η Λ / Ω 0 1 , thus the perturbation treatment is kept valid. In such a case, in the greener region of Figure 3, the correction would be about 0.8 , implying a quite low efficiency, of about 0.2 . On the contrary, in the most problematic region of Figure 2, consisting of the yellow parts corresponding to values close to 0.03 , the efficiency would be about 0.97 .

4. Discussion

In this paper, we have considered the effects of the interaction with a spin bath on the efficiency of a population transfer realized through a STIRAP process. Our present study is an evolution of the earlier studyreported in ref. [56], where the system–bath interaction has been considered in the RWA. Here, in order to improve the analysis, we have also considered the effects of the counter-rotating terms. Nevertheless, since the model analyzed is not exactly solvable and its numerical treatment would be challenging, we have faced the problem through a perturbation approach. In particular, second order corrections have been calculated by evaluating the Dyson series truncated to the second order contributions. The complete deviation of the efficiency of the population transfer from unity, according to Equation (29), turns out to be proportional to the square of η Λ , which plays the role of a perturbation parameter. It is important to note that η Λ = k η k ( 1 ) 2 and r η η Λ = k η k ( 2 ) 2 , which correspond to the effective strengths of the couplings with the environment involving g 1 e transitions and g 2 e transitions, respectively. Beyond this fact, according to Equation (29), the correction to the efficiency of population transfer is proportional to the real part of the integral J ( T ) defined in Equation (30).
Remarkably, in a second order perturbation treatment under the RWA, the correction to the efficiency would be zero, thus predicting a perfect population transfer. On the contrary, beyond such an approximation, deviations come up. As a matter of fact, our analysis shows that, in the second perturbation treatment, the counter-rotating terms in the interaction are the only ones really contributing to the corrections. Indeed, in the integrand of J ( T ) only contributions proportional to the phase factor exp [ ( ω + ν ) ( w w ) ] appear to be traceable back the rapidly oscillating terms in the interaction (namely, the counter-rotating terms), while no term proportional to a phase factor associated with a lower frequency ω ν (the rotating terms) is present. This implies that, in a model involving RWA, the correction is zero, which is perfectly in agreement with the results of ref. [56] where the efficiency always exhibits a plateau in the weak coupling limit and up to the weak–intermediate coupling regime.
The numerical evaluation of the integral J ( T ) reported in Figure 2 and Figure 3 shows that, in the range of the parameters analyzed the effects of the environment are mainly negligible, except for specific regions where the ratio r η assumes high values, since this implies that the coupling constants η k ( 2 ) are significant. Further numerical results (not reported in this paper because they would essentially provide uniform functions corresponding to zero values) show in a clear way that higher and higher values of ν imply smaller and smaller corrections, due to the fact that counter-rotating terms become more and more rapidly oscillating and thus ineffective. Summing up, negligibility of the effects of the environment is obtained for high values of ω (according to the plates shown in the paper) or high values of ν (according to other simulations), provided r η is not exceptionally large. Naturally, the values of η Λ are to be kept small enough to maintain the validity of our analysis based on perturbation theory.
It is worth concluding with two more comments concerning a comparison with the results reported in ref. [56]. First, since our second-order perturbation treatment is valid only in the weak coupling regime, we cannot obtain corrections in the strong coupling regime where a generalized quantum Zeno effect were predicted to occur. Second, the shape of the pulses considered here is properly Gaussian, while in ref. [56], different expressions have been considered. Nevertheless, in spite of the more complicated analytical expression, the pulses in ref. [56] are essentially Gaussian too, in the sense that they only slightly differ form the Gaussian counterparts, which makes the two situations comparable.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge financial support from the FFR2021 grant from the University of Palermo, Italy.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Adiabatic Approximation

Here we recall the proof of the adiabatic theorem [1,2]. Consider a system subjected to time-dependent Hamiltonian H ( t ) which can be diagonalized at every instant of time in such a way that H ( t ) ϕ n j ( t ) = E n ( t ) ϕ n j ( t ) , where each energy eigenspace has a possible degeneracy taken into account through the index j. Expanding the state of the system in terms of such instantaneous eigenstates,
ψ ( t ) = n j a n j ( t ) e i t 0 t E n ( s ) d s ϕ n j ( t ) ,
and inserting this expression in the Schrödinger equation leads to a set of equations for the coefficients a n j . Now, from ϕ n j ( t ) | ϕ n j ( t ) = δ n n δ j j (with δ i j the Kronecker delta) and ϕ n j ( t ) | H ( t ) | ϕ n j ( t ) = E n ( t ) δ n n δ j j one obtains: t ( ϕ n j ( t ) | ϕ n j ( t ) ) = 0 and t ( ϕ n j ( t ) | H ( t ) | ϕ n j ( t ) ) = 0 unless n = n and j = j . After some algebra, one finds that for n k the following relation occurs: ϕ n j ( t ) | ϕ ˙ k l ( t ) = ϕ n l ( t ) | H ˙ ( t ) | ϕ k j ( t ) / ( E k ( t ) E n ( t ) ) . On this basis once can eventually write down the following set of equations:
a ˙ n j ( t ) = ϕ n j ( t ) | ϕ ˙ n j ( t ) a n j ( t ) l j ϕ n j ( t ) | ϕ ˙ n l ( t ) a n l ( t ) k n l e i 0 t ( E k ( s ) E n ( s ) ) d s ϕ n l ( t ) | H ˙ ( t ) | ϕ k j ( t ) E k ( t ) E n ( t ) a k l ( t ) ,
with time derivative denoted by the dot. Under the assumption of very small matrix elements of the operator H ˙ , one can neglect the terms in the second line, and for non degenerate eigenspaces one eventually obtains a ˙ n j ϕ n j | ϕ ˙ n j a n j or, equivalently, after using χ | χ ˙ = χ ˙ | χ ,
a ˙ n j ( t ) ϕ ˙ n j ( t ) | ϕ n j ( t ) a n j ( t ) .
The complete evolution is then given by:
ψ ( t ) n j a n j ( 0 ) e t 0 t ϕ ˙ n j ( s ) | ϕ n j ( s ) d s e i t 0 t E n ( s ) d s ϕ n j ( t ) ,
where the quantity ϕ ˙ n j ( s ) | ϕ n j ( s ) , known as the geometric phase, is an imaginary number. Therefore, if the coefficients of the expansion of ϕ n j ( s ) with respect to a given basis are all real, the phase turns out to be the sum of real numbers, which is then supposed to be equal to zero. This is the case for the eigenstates of the STIRAP Hamiltonian considered above.

Appendix B. Second Order Term for the Zero-Temperature Bath

In the calculation of the corrections up to the second order we need some quantities. Since g 1 | ϕ ˚ ( t ) = g 1 U 2 ( t , t 0 ) ϕ = ϕ U 2 ( t , t 0 ) g 1 * , considered the expressions of the states of the adiabatic basis and the fact that at the initial time ( t = t 0 ) of the counter-intuitive sequence 0 = g 1 , + = e and = g 2 , one immediately finds:
g 1 | e ˚ ( t ) = 0 ,
g 1 | g ˚ 1 ( t ) = cos θ ( t ) = Ω 2 ( t ) / Ω L ( t )
g 1 | g ˚ 2 ( t ) = sin θ ( t ) = Ω 1 ( t ) / Ω L ( t )
e ˚ ( t ) | e ˚ ( t ) = e i t t E + ( s ) d s cos φ ( t ) cos φ ( t ) + e i t t E ( s ) d s sin φ ( t ) sin φ ( t ) ,
where θ ( t ) and φ ( t ) given by Equations (15) and (16).
Since g 1 | e ˚ ( t ) = 0 , all the first order terms turn out to be zero, when it is calculated any probability to find the three-state system in the state g 1 . Indeed, whatever the bath states Ψ B and Ψ B , the matrix elements Ψ B g 1 H ¯ S ( ω ) ( t ) g 1 Ψ B and Ψ B g 1 H ¯ A B ( t ) g 1 Ψ B give rise to terms proportional to the overlap g 1 | e ˚ ( t ) (or its adjoint), which is zero.
Let us now then consider the second order contributions. First of all, observe that the matrix element { } | g 1 H ¯ G ( w ) H ¯ G ( w ) g 1 | { } admits potential contributions coming only from the terms H ¯ S ( ω ) ( w ) H ¯ S ( ω ) ( w ) and H ¯ A B ( w ) H ¯ A B ( w ) . Indeed, H ¯ S ( ω ) ( w ) H ¯ A B ( w ) and H ¯ A B ( w ) H ¯ S ( ω ) ( w ) involve one spin flip coming from H ¯ A B while we are considering the matrix element between g 1 { } and itself. That said, let us then focus on the other two terms:
{ } g 1 H ¯ S ( ω ) ( w ) H ¯ S ( ω ) ( w ) g 1 { } = e 2 i ν ( w w ) e ( w ) | e ( w ) × m m Ω m ( w ) Ω m ( w ) g 1 | g ˚ m ( w ) g ˚ m ( w ) | g 1 = 0 ,
whose being zero is based on Equations (A6) and (A7), and
{ } g 1 H ¯ A B ( w ) H ¯ A B ( w ) g 1 { } = e i ( ν + ω ) ( w w ) e ( w ) | e ( w ) × g 1 | g ˚ m ( w ) g ˚ m ( w ) | g 1 k η k ( m ) η k ( m ) ,
which, after using Equations (A6)–(A8), gives the integrand in Equation (27).

References

  1. Messiah, A. Quantum Mechanics; Dover Publications, Inc.: Mineola, NY, USA, 1999; Available online: https://archive.org/details/quantummechanics0000mess/ (accessed on 2 February 2024).
  2. Griffiths, D.J.; Schroeter, D.F. Introduction to Quantum Mechanics; Cambridge University Press: New York, NY, USA, 2018. [Google Scholar] [CrossRef]
  3. Albash, T.; Lidar, D.A. Adiabatic quantum computation. Rev. Mod. Phys. 2018, 90, 015002. [Google Scholar] [CrossRef]
  4. Santos, A.C.; Sarandy, M.S. Sufficient conditions for adiabaticity in open quantum systems. Phys. Rev. A 2020, 102, 052215. [Google Scholar] [CrossRef]
  5. Il’in, N.; Aristova, A.; Lychkovskiy, O. Adiabatic theorem for closed quantum systems initialized at finite temperature. Phys. Rev. A 2021, 104, L030202. [Google Scholar] [CrossRef]
  6. Pyshkin, P.V.; Luo, D.W.; Wu, L.A. Self-protected adiabatic quantum computation. Phys. Rev. A 2022, 106, 012420. [Google Scholar] [CrossRef]
  7. Coello Pérez, E.A.; Bonitati, J.; Lee, D.; Quaglioni, S.; Wendt, K.A. Quantum state preparation by adiabatic evolution with custom gates. Phys. Rev. A 2022, 105, 032403. [Google Scholar] [CrossRef]
  8. Vitanov, N.V.; Fleischhauer, M.; Shore, B.W.; Bergmann, K. Coherent manipulation of atoms molecules by sequential laser pulses. Adv. At. Mol. Opt. Phys. 2001, 46, 55–190. [Google Scholar] [CrossRef]
  9. Vitanov, N.V.; Halfmann, T.; Shore, B.W.; Bergmann, K. Laser-induced population transfer by adiabatic passage. Ann. Rev. Phys. Chem. 2001, 52, 763–809. [Google Scholar] [CrossRef] [PubMed]
  10. Bergmann, K.; Theuer, H.; Shore, B.W. Coherent population transfer among quantum states of atoms and molecules. Rev. Mod. Phys. 1998, 70, 1003–1025. [Google Scholar] [CrossRef]
  11. Král, P.; Thanopulos, I.; Shapiro, M. Colloquium: Coherently controlled adiabatic passage. Rev. Mod. Phys. 2007, 79, 53–77. [Google Scholar] [CrossRef]
  12. Bergmann, K.; Nägerl, H.-C.; Panda, C.; Gabrielse, G.; Miloglyadov, E.; Quack, M.; Seyfang, G.; Wichmann, G.; Ospelkaus, S.; Kuhn, A.; et al. Roadmap on STIRAP applications. J. Phys. B At. Mol. Opt. Phys. 2019, 52, 202001. [Google Scholar] [CrossRef]
  13. Blekos, K.; Stefanatos, D.; Paspalakis, E. Performance of superadiabatic stimulated Raman adiabatic passage in the presence of dissipation and Ornstein-Uhlenbeck dephasing. Phys. Rev. A 2020, 102, 023715. [Google Scholar] [CrossRef]
  14. Ahmadinouri, F.; Hosseini, M.; Sarreshtedari, F. Stimulated Raman adiabatic passage: Effects of system parameters on population transfer. Chem. Phys. 2020, 539, 110960. [Google Scholar] [CrossRef]
  15. Cinins, A.; Bruvelis, M.; Bezuglov, N.N. Study of the adiabatic passage in tripod atomic systems in terms of the Riemannian geometry of the Bloch sphere. J. Phys. B At. Mol. Opt. Phys. 2022, 55, 234003. [Google Scholar] [CrossRef]
  16. Dogra, S.; Paraoanu, G.S. Perfect stimulated Raman adiabatic passage with imperfect finite-time pulses. J. Phys. B At. Mol. Opt. Phys. 2022, 55, 174001. [Google Scholar] [CrossRef]
  17. Liu, K.; Sugny, D.; Chen, X.; Guérin, S. Optimal pulse design for dissipative-stimulated Raman exact passage. Entropy 2023, 25, 790. [Google Scholar] [CrossRef] [PubMed]
  18. Genov, G.T.; Rochester, S.; Auzinsh, M.; Jelezko, F.; Budker, D. Robust two-state swap by stimulated Raman adiabatic passage. J. Phys. B At. Mol. Opt. Phys. 2023, 56, 054001. [Google Scholar] [CrossRef]
  19. Ni, K.-K.; Ospelkaus, S.; de Miranda, M.H.G.; Peér, A.; Neyenhuis, B.; Zirbel, J.J.; Kotochigova, S.; Julienne, P.S.; Jin, D.S.; Ye, J. A high phase-space-density gas of polar molecules. Science 2008, 322, 231–235. [Google Scholar] [CrossRef]
  20. Danzl, J.G.; Haller, E.; Gustavsson, M.; Mark, M.J.; Hart, R.; Bouloufa, N.; Dulieu, O.; Ritsch, H.; Nägerl, H.-C. Quantum Gas of deeply bound ground state molecules. Science 2008, 321, 1062–1066. [Google Scholar] [CrossRef] [PubMed]
  21. Danzl, J.G.; Mark, M.J.; Haller, E.; Gustavsson, M.; Hart, R.; Aldegunde, J.; Hutson, J.M.; Nägerl, H.-C. An ultracold high-density sample of rovibronic ground-state molecules in an optical lattice. Nat. Phys. 2010, 6, 265–270. [Google Scholar] [CrossRef]
  22. Klein, J.; Beil, F.; Halfmann, T. Robust Population transfer by stimulated Raman adiabatic passage in a Pr3+:Y2SiO5 Crystal. Phys. Rev. Lett. 2007, 99, 113003. [Google Scholar] [CrossRef] [PubMed]
  23. Alexander, A.L.; Lauro, R.; Louchet, A.; Chanelière, T.; Le Gouët, J.L. Stimulated Raman adiabatic passage in Tm3+:YAG. Phys. Rev. B 2008, 78, 144407. [Google Scholar] [CrossRef]
  24. Golter, D.A.; Wang, H. Optically driven Rabi oscillations and adiabatic passage of single electron spins in diamond. Phys. Rev. Lett. 2014, 112, 116403. [Google Scholar] [CrossRef]
  25. Yale, C.G.; Heremans, F.J.; Zhou, B.B.; Auer, A.; Burkard, G.; Awschalom, D.D. Optical manipulation of the Berry phase in a solid-state spin qubit. Nat. Photon. 2016, 10, 184–189. [Google Scholar] [CrossRef]
  26. Wolfowicz, G.; Heremans, F.J.; Anderson, C.P.; Kanai, S.; Seo, H.; Gali, A.; Galli, G.; Awschalom, D.D. Quantum guidelines for solid-state spin defects. Nat. Rev. Mater. 2021, 6, 906925. [Google Scholar] [CrossRef]
  27. Zhou, B.; Baksic, A.; Ribeiro, H.; Yale, C.G.; Heremans, F.J.; Jerger, P.C.; Auer, A.; Burkard, G.; Clerk, A.A.; Awschalom, D.D. Accelerated quantum control using superadiabatic dynamics in a solid-state lambda system. Nat. Phys. 2017, 13, 330–334. [Google Scholar] [CrossRef]
  28. Baksic, A.; Ribeiro, H.; Clerk, A.A. Speeding up adiabatic quantum state transfer by using dressed states. Phys. Rev. Lett. 2016, 116, 230503. [Google Scholar] [CrossRef]
  29. Varguet, H.; Rousseaux, B.; Dzsotjan, D.; Jauslin, H.R.; Guérin, S.; Colas des Francs, G. Dressed states of a quantum emitter strongly coupled to a metal nanoparticle. Opt. Lett. 2016, 41, 4480–4483. [Google Scholar] [CrossRef]
  30. Castellini, A.; Jauslin, H.R.; Rousseaux, B.; Dzsotjan, D.; Colas des Francs, G.; Messina, A.; Guérin, S. Quantum plasmonics with multi-emitters: Application to stimulated Raman adiabatic passage. Eur. Phys. J. D 2018, 72, 223. [Google Scholar] [CrossRef]
  31. Kubo, Y. Turn to the dark side. Nat. Phys. 2016, 12, 21–22. [Google Scholar] [CrossRef]
  32. Xu, H.K.; Song, C.; Liu, W.Y.; Xue, G.M.; Su, F.F.; Deng, H.; Tian, Y.; Zheng, D.N.; Han, S.; Zhong, Y.P.; et al. Coherent population transfer between uncoupled or weakly coupled states in ladder-type superconducting qutrits. Nat. Commun. 2016, 7, 11018. [Google Scholar] [CrossRef] [PubMed]
  33. Kumar, K.S.; Vepsäläinen, A.; Danilin, S.; Paraoanu, G.S. Stimulated Raman adiabatic passage in a three-level superconducting circuit. Nat. Commun. 2016, 7, 10628. [Google Scholar] [CrossRef] [PubMed]
  34. Soresen, J.L.; Moller, D.; Iversen, T.; Thomsen, J.B.; Jensen, F.; Staanum, P.; Voigt, D.; Drewsen, M. Efficient coherent internal state transfer in trapped ions using stimulated Raman adiabatic passage. New J. Phys. 2006, 8, 261. [Google Scholar] [CrossRef]
  35. Higgins, G.; Pokorny, F.; Zhang, C.; Bodart, Q.; Hennrich, M. Coherent control of a single trapped Rydberg ion. Phys. Rev. Lett. 2017, 119, 220501. [Google Scholar] [CrossRef]
  36. Fedoseev, V.; Luna, F.; Hedgepeth, I.; Löffler, W.; Bouwmeester, D. Stimulated Raman adiabatic passage in optomechanics. Phys. Rev. Lett. 2021, 126, 113601. [Google Scholar] [CrossRef]
  37. Guéry-Odelin, D.; Ruschhaupt, A.; Kiely, A.; Torrontegui, E.; Martínez-Garaot, S.; Muga, J.G. Shortcuts to adiabaticity: Concepts, methods, and applications. Rev. Mod. Phys. 2019, 91, 045001. [Google Scholar] [CrossRef]
  38. Vitanov, N.V. High-fidelity multistate stimulated Raman adiabatic passage assisted by shortcut fields. Phys. Rev. A 2020, 102, 023515. [Google Scholar] [CrossRef]
  39. Evangelakos, V.; Paspalakis, E.; Stefanatos, D. Optimal STIRAP shortcuts using the spin-to-spring mapping. Phys. Rev. A 2023, 107, 052606. [Google Scholar] [CrossRef]
  40. Stefanatos, D.; Paspalakis, E. Optimal shortcuts of stimulated Raman adiabatic passage in the presence of dissipation. Philos. Trans. R. Soc. Math. Phys. Eng. Sci. 2022, 380, 20210283. [Google Scholar] [CrossRef] [PubMed]
  41. Messikh, C.; Messikh, A. Robust stimulated Raman shortcuts to adiabatic passage with deep learning. EPL (Europhys. Lett.) 2022, 140, 48003. [Google Scholar] [CrossRef]
  42. Stefanatos, D.; Blekos, K.; Paspalakis, E. Robustness of STIRAP shortcuts under Ornstein-Uhlenbeck noise in the energy levels. Appl. Sci. 2020, 10, 1580. [Google Scholar] [CrossRef]
  43. Genov, G.T.; Vitanov, N.V. Dynamical suppression of unwanted transitions in multistate quantum systems. Phys. Rev. Lett. 2013, 110, 133002. [Google Scholar] [CrossRef] [PubMed]
  44. Yatsenko, L.P.; Shore, B.W.; Bergmann, K. Detrimental consequences of small rapid laser fluctuations on stimulated Raman adiabatic passage. Phys. Rev. A 2014, 89, 013831. [Google Scholar] [CrossRef]
  45. Vitanov, N.V.; Stenholm, S. Population transfer via a decaying state. Phys. Rev. A 1997, 56, 1463–1471. [Google Scholar] [CrossRef]
  46. Breuer, H.-P.; Petruccione, F. The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002. [Google Scholar] [CrossRef]
  47. Gardiner, C.W.; Zoller, P. Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  48. Ivanov, P.A.; Vitanov, N.V.; Bergmann, K. Spontaneous emission in stimulated Raman adiabatic passage. Phys. Rev. A 2005, 72, 053412. [Google Scholar] [CrossRef]
  49. Akram, M.J.; Saif, F. Adiabatic population transfer based on a double stimulated Raman adiabatic passage. J. Russ. Laser Res. 2014, 35, 547–554. [Google Scholar] [CrossRef]
  50. Sukharev, M.; Malinovskaya, S.A. Stimulated Raman adiabatic passage as a route to achieving optical control in plasmonics. Phys. Rev. A 2012, 86, 043406. [Google Scholar] [CrossRef]
  51. Mathisen, T.; Larson, J. Liouvillian of the open STIRAP problem. Entropy 2018, 20, 20. [Google Scholar] [CrossRef]
  52. Davies, E.B.; Spohn, H. Open quantum systems with time-dependent Hamiltonians and their linear response. J. Stat. Phys. 1978, 19, 511–523. [Google Scholar] [CrossRef]
  53. Scala, M.; Militello, B.; Messina, A.; Vitanov, N.V. Stimulated Raman adiabatic passage in an open quantum system: Master equation approach. Phys. Rev. A 2010, 81, 053847. [Google Scholar] [CrossRef]
  54. Scala, M.; Militello, B.; Messina, A.; Vitanov, N.V. Stimulated Raman adiabatic passage in a Λ system in the presence of quantum noise. Phys. Rev. A 2011, 83, 012101. [Google Scholar] [CrossRef]
  55. Onizhuk, M.; Miao, K.C.; Blanton, J.P.; Ma, H.; Anderson, C.P.; Bourassa, A.; Awschalom, D.D.; Galli, G. Probing the Coherence of Solid-State Qubits at Avoided Crossings. Phys. Rev. X Quantum 2021, 2, 010311. [Google Scholar] [CrossRef]
  56. Militello, B.; Napoli, A. Adiabatic manipulation of a system interacting with a spin bath. Entropy 2023, 15, 2028. [Google Scholar] [CrossRef]
  57. Fischer, J.; Breuer, H.-P. Correlated projection operator approach to non-Markovian dynamics in spin baths. Phys. Rev. A 2007, 76, 052119. [Google Scholar] [CrossRef]
  58. Ferraro, E.; Breuer, H.-P.; Napoli, A.; Jivulescu, M.A.; Messina, A. Non-Markovian dynamics of a single electron spin coupled to a nuclear spin bath. Phys. Rev. B 2008, 78, 064309. [Google Scholar] [CrossRef]
  59. Bhattacharya, S.; Misra, A.; Mukhopadhyay, C.; Pati, A.K. Exact master equation for a spin interacting with a spin bath: Non-Markovianity and negative entropy production rate. Phys. Rev. A 2017, 95, 012122. [Google Scholar] [CrossRef]
Figure 1. STIRAP scheme: two lower states g 1 and g 2 are coupled to an upper state e through suitable pulses. The dashed line represents a ‘virtual level’ of energy ν . The inset shows the typical shape of the pulses.
Figure 1. STIRAP scheme: two lower states g 1 and g 2 are coupled to an upper state e through suitable pulses. The dashed line represents a ‘virtual level’ of energy ν . The inset shows the typical shape of the pulses.
Physics 06 00032 g001
Figure 2. The quantity 2 { J ( T ) } (in units of τ 2 ) as a function of ω (in units of ν and logarithmic scale) and Δ (in units of ν ) in the case r η = 1 (a) and r η = 1 (b). The other parameters are: Ω 0 τ = 10 , T / τ = 5 , ν τ = 10 . See text for details.
Figure 2. The quantity 2 { J ( T ) } (in units of τ 2 ) as a function of ω (in units of ν and logarithmic scale) and Δ (in units of ν ) in the case r η = 1 (a) and r η = 1 (b). The other parameters are: Ω 0 τ = 10 , T / τ = 5 , ν τ = 10 . See text for details.
Physics 06 00032 g002
Figure 3. The quantity 2 { J ( T ) } (in units of τ 2 ) as a function of ω (in units of ν and logarithmic scale) and r η (in logarithmic scale) for Δ = 0.01 ν (a) and Δ = 0.05 ν (b). The other parameters are: Ω 0 τ = 10 , T / τ = 5 , ν τ = 10 . See text for details.
Figure 3. The quantity 2 { J ( T ) } (in units of τ 2 ) as a function of ω (in units of ν and logarithmic scale) and r η (in logarithmic scale) for Δ = 0.01 ν (a) and Δ = 0.05 ν (b). The other parameters are: Ω 0 τ = 10 , T / τ = 5 , ν τ = 10 . See text for details.
Physics 06 00032 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Militello, B.; Napoli, A. Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation. Physics 2024, 6, 483-495. https://doi.org/10.3390/physics6020032

AMA Style

Militello B, Napoli A. Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation. Physics. 2024; 6(2):483-495. https://doi.org/10.3390/physics6020032

Chicago/Turabian Style

Militello, Benedetto, and Anna Napoli. 2024. "Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation" Physics 6, no. 2: 483-495. https://doi.org/10.3390/physics6020032

APA Style

Militello, B., & Napoli, A. (2024). Adiabatically Manipulated Systems Interacting with Spin Baths beyond the Rotating Wave Approximation. Physics, 6(2), 483-495. https://doi.org/10.3390/physics6020032

Article Metrics

Back to TopTop