Next Article in Journal
Shannon Entropy in Confined He-Like Ions within a Density Functional Formalism
Next Article in Special Issue
How Does the Planck Scale Affect Qubits?
Previous Article in Journal
Heterogeneous Network Architecture for Integration of AI and Quantum Optics by Means of Multiple-Valued Logic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems

1
SUNY Polytechnic Institute, Albany, NY 12203, USA
2
Air Force Research Laboratory, Information Directorate, Rome, NY 13441, USA
*
Author to whom correspondence should be addressed.
Quantum Rep. 2020, 2(1), 166-188; https://doi.org/10.3390/quantum2010011
Submission received: 16 January 2020 / Revised: 23 January 2020 / Accepted: 26 January 2020 / Published: 5 February 2020
(This article belongs to the Special Issue Relevance of Information Geometry in Quantum Information Science)

Abstract

:
We present an information geometric analysis of off-resonance effects on classes of exactly solvable generalized semi-classical Rabi systems. Specifically, we consider population transfer performed by four distinct off-resonant driving schemes specified by su 2 ; time-dependent Hamiltonian models. For each scheme, we study the consequences of a departure from the on-resonance condition in terms of both geodesic paths and geodesic speeds on the corresponding manifold of transition probability vectors. In particular, we analyze the robustness of each driving scheme against off-resonance effects. Moreover, we report on a possible tradeoff between speed and robustness in the driving schemes being investigated. Finally, we discuss the emergence of a different relative ranking in terms of performance among the various driving schemes when transitioning from on-resonant to off-resonant scenarios.

1. Introduction

It is well-known that the resonance phenomenon concerns the amplification that occurs when the frequency of an externally applied field (for instance, a magnetic field) is equal to a natural characteristic frequency of the system on which it acts. When a field is applied at the resonant frequency of another system, the system oscillates at a higher amplitude than when the field is applied at a non-resonant frequency. In an actual experimental laboratory, off-resonance effects can occur for a number of reasons. For instance, a common source of off-resonance effects is the presence of unforeseen magnetic field gradients that are not part of the externally applied magnetic field. These gradients can cause signal loss and, as a consequence, signals at the “wrong” resonant frequency (that is, off-resonance). Furthermore, when the intensity of an applied magnetic field gradient is changed in a laboratory, decaying eddy currents can produce time-varying off-resonant effects [1]. Additionally, the presence of concomitant gradients can create off-resonance effects in the form of x y -spatial variations in the Larmor frequency when a magnetic field gradient in the z-direction is applied [2]. For the interested reader, we discuss two illustrative examples concerning the departure from the on-resonance condition in the context of classical and quantum mechanics in Appendix A.
The resonance phenomenon happens to play a key role in problems where one seeks to drive an initial source state towards a final target state. For instance, typical problems of this kind can be found when controlling population transfer in quantum systems [3,4,5] or when searching for a target state [6,7]. In the framework of quantum control theory, it is commonly believed that an off-resonant driving scheme cannot help achieving population transfer with high fidelity between two quantum states. However, there do exist theoretical investigations relying on the phenomenon of superoscillations [8] where population transfer occurs by using only off-resonant driving fields [9]. In particular, during resonant superoscillations, the quantum system is temporarily excited and behaves as if truly driven at resonance [9]. Furthermore, the existence of suboptimal transfer efficiencies, caused by the sequential use of two lasers in order to physically realize a transition from highly excited Rydberg states to ground states [10,11], leads to question the possibility of transition from a ground state to the Rydberg state by means of lasers far-off-resonant with respect to the transition frequency. Indeed, motivated by such possibility, driving schemes specified by periodic square-wells and Gaussian pulses in order to implement coherent transfer with far-off-resonant driving fields have been proposed in reference [12].
It is also well-known that on-resonance driving is not a good control strategy for a quantum system if, from an experimental standpoint, one allows the Rabi frequency (that is, the frequency that describes the coupling between the driving field and the two-level system) to be higher than the Larmor frequency (that is, the characteristic transition frequency of the two-level system in the absence of the external driving field) [13]. This is called the strong driving regime. In such a scenario, the so-called rotating wave approximation (RWA, [14]) is not applicable and the system is driven with an external field whose amplitude is greater than or equal to the energy splitting between the system’s states. For an analysis on the deviation from sinusoidal Rabi oscillations by studying the time-evolution dynamics of a strongly driven dressed electron spin in silicon, we refer to reference [15]. Furthermore, for an explicit manifestation of strong sensitivity to the initial phase of the driving field in the dynamics of a general semi-classical Rabi model in regimes of arbitrary strong driving, we refer to reference [16].
Based upon our considerations, it is clear that off-resonance phenomena exhibit both experimental (for instance, limited control of spurious electromagnetic signals and strong driving regimes) and theoretical (for instance, sensitivity to initial conditions and emergence of more complex dynamical scenarios) interest. In the framework of analog quantum searching, we have recently presented a detailed investigation concerning the physical connection between quantum search Hamiltonians and exactly solvable time-dependent two-level quantum systems in reference [7]. In references [17,18], instead, we analyzed the possibility of modifying the original Farhi-Gutmann Hamiltonian quantum search algorithm [19] in order to speed up the procedure for producing a suitably distributed unknown normalized quantum mechanical state provided only a nearly optimal fidelity is sought. In reference [20], we presented an information geometric characterization of the oscillatory or monotonic behavior of statistically parametrized squared probability amplitudes originating from special functional forms of the Fisher information function: constant, exponential decay, and power-law decay. Furthermore, for each case, we computed both the computational speed and the availability loss of the corresponding physical processes by exploiting a convenient Riemannian geometrization of useful thermodynamical concepts. Finally, building upon our works presented in references [7,17,18,20], we presented in reference [21] an information geometric analysis of geodesic speeds and entropy production rates in geodesic motion on manifolds of parametrized quantum states. These pure states emerged as outputs of suitable su 2 ; time-dependent Hamiltonian evolutions used to describe distinct types of analog quantum search schemes viewed as driving strategies. In particular, by evaluating the geodesic speed and the total entropy production along the optimum transfer paths in a number of physical scenarios of interest in analog quantum search problems (for instance, constant, oscillatory, power-law decay, and exponential decay of the driving magnetic fields), we showed in an explicit quantitative manner that to a faster transfer there corresponds necessarily a higher entropy production rate. Thus, we concluded that lower entropic efficiency values do appear to accompany higher entropic geodesic speed values in quantum transfer processes. Our information geometric analysis in reference [21] was limited to the case in which the on-resonance driving condition was satisfied for all quantum driving strategies. However, based upon our previous considerations, it would be of theoretical interest to investigate from an information geometric perspective the effects of off-resonance effects on these driving strategies. In particular, we would like to address the following questions:
(i)
What are the main effects of deviations from the on-resonance condition in these su 2 ; time-dependent Hamiltonian evolutions?
(ii)
Do off-resonance effects modify (with respect to the on-resonance scenario) the relative ranking in terms of performance (quantified in terms of geodesic speed and/or minimum transfer time) among the driving schemes being considered?
(iii)
Are there driving schemes that are especially robust against deviations from the on-resonance condition and that, in addition, are capable of reaching sufficiently high fidelity values?
The rest of the paper is organized as follows. In Section 2, we present our driving strategies in terms of classes of exactly solvable generalized semi-classical Rabi systems. Specifically, we consider the population transfer performed by four distinct off-resonant driving schemes specified by su 2 ; time-dependent Hamiltonian models. In Section 3, we propose our information geometric analysis of off-resonance effects on such classes of systems. For each scheme, we study the consequences of a departure from the on-resonance condition in terms of both geodesic paths and geodesic speeds on the corresponding manifold of transition probability vectors. In particular, we analyze the robustness of each driving scheme against off-resonance effects. Moreover, we report on a possible tradeoff between speed and robustness in the driving schemes being investigated. Finally, we discuss the emergence of a different relative ranking in terms of performance among the various driving schemes when transitioning from the on-resonant to the off-resonant scenarios. Our conclusive remarks appear in Section 4. Auxiliary illustrative examples together with some technical details are placed in Appendices Appendix A and Appendix B.

2. The su ( 2 ; ) Quantum Driving

In this section, we introduce our driving strategies in terms of classes of exactly solvable generalized semi-classical Rabi systems [22,23].
The quantum evolution we take into consideration is specified by means of an Hamiltonian operator H su 2 ; written as the most general linear superposition of the three traceless and anti-Hermitian generators i σ x , i σ y , i σ z of su 2 ; , the Lie algebra of the special unitary group SU 2 ; ,
H su 2 ; t = def a t i σ x + b t i σ y + c t i σ z .
In Equation (1), a t , b t , and c t are time-dependent complex coefficients while σ = def σ x , σ y , σ z denotes the Pauli vector operator. In particular, by putting a t = def i ω x t , b t = def i ω y t , and c t = def i Ω t , Equation (1) yields
H su 2 ; t = def ω x t σ x + ω y t σ y + Ω t σ z .
In the framework of su 2 ; Hamiltonian models, ω t = def ω x t i ω y t = ω H t e i ϕ ω t and Ω t are the so-called complex transverse field and real longitudinal field, respectively. Obviously, ω H t denotes the modulus of ω t . In what follows, we assume that longitudinal fields Ω t are oriented along the z-axis while transverse fields ω t lie in the x y -plane. Taking into consideration the quantum evolution of an electron (or, more generally, a spin- 1 / 2 particle) in an external time-dependent magnetic field B t , the Hamiltonian H su 2 ; t in Equation (2) can be rewritten as
H su 2 ; t = def μ · B t ,
where μ = def e ħ / 2 m c σ denotes the magnetic moment of the electron with μ Bohr = def e ħ / ( 2 m c ) 9.27 × 10 21 cgs being the so-called Bohr magneton. The quantity m 9.11 × 10 31 MKSA is the mass of an electron while e 1.60 × 10 19 MKSA denotes the absolute value of the electric charge of an electron. Moreover, c 3 × 10 8 MKSA and ħ 1.05 × 10 34 MKSA are the speed of light and the reduced Planck constant, respectively. The magnetic field B t in Equation (1) can be recast as
B t = def B t + B t ,
with B t = def B x t x ^ + B y t y ^ and B t = def B z t z ^ . Equating the Hamiltonians in Equations (2) and (3) and making use of the magnetic field decomposition in Equation (4), the connection between the set of field intensities ω H t , Ω H t and the set of magnetic field intensities B t , B t becomes clear. More specifically, we observe that B t ω H t and B t Ω H t = def Ω t . The exact connection between the components B x t , B y t , B z t and ω x t , ω y t , Ω t is expressed by
B x t = 2 m c e ħ ω x t , B y t = 2 m c e ħ ω y t , and B z t = 2 m c e ħ Ω t .
Moreover, in terms of field intensities, we get
B t = 2 m c e ħ ω H t , and B t = 2 m c e ħ Ω H t .
Studying the quantum evolution of an electron governed by the Hamiltonian in Equation (3) by means of exact analytical expressions of complex probability amplitudes and/or real transition probabilities from an initial source state to a final target state can be quite challenging. In particular, the canonical matrix representations of H su 2 ; in Equation (2) and U t , the unitary evolution operator arising from H su 2 ; , with i ħ U ˙ t = H su 2 ; U t and U ˙ = def t U , are given by
H su 2 ; = def Ω t ω t ω * t Ω t , and U t = def α t β t β * t α * t ,
respectively. The unitarity of the quantum evolution demands that the complex probability amplitudes α t and β t fulfill the normalization condition, α t 2 + β t 2 = 1 . Once the unitary evolution operator U t in Equation (7) is given, the temporal evolution of a quantum source state s ,
s = def x w + 1 x 2 w ,
can be characterized by the following transformation law,
x 1 x 2 α t x + β t 1 x 2 β * t x + α * t 1 x 2 ,
with x = def w | s denoting the quantum overlap. The set of orthonormal state vectors w , w generate the two-dimensional search space of the N = 2 n -dimensional complex Hilbert space H 2 n . Therefore, using Equations (7), (8), and (9), the probability that the source state s transitions into the target state w under U t becomes
P s w t = def w | U t | s 2 = α t 2 x 2 + β t 2 1 x 2 + α t β * t + α * t β t x 1 x 2 .
It is evident from Equation (10) that in order to calculate the exact analytical expression of transition probabilities, one must have the exact analytical formula of the evolution operator U t expressed in terms of the complex probability amplitudes α t and β t . Ideally, having specified the fields ω t and Ω t by means of physical arguments, one would solve the coupled system of first order ordinary differential equations with time-dependent coefficients emerging from the relation i ħ U ˙ t = H su 2 ; U t with U 0 = I ,
i ħ α ˙ t = Ω t α t ω t β * t , and i ħ β ˙ t = ω t α * t + Ω t β t ,
where α 0 = 1 and β 0 = 0 . However, it is often the case that this general approach does not yield exact analytical solutions. It is recognized that it is rather difficult to find exact analytical solutions in time-dependent two-level quantum systems specified by equations as the ones in Equation (11). A very powerful technique for investigating these types of quantum evolutions is the so-called rotating coordinates technique. This method, originally proposed by Rabi, Ramsey, and Schwinger in magnetic resonance problems, can be presented as a three-step technique: First, recast the original problem from a stationary to a rotating frame of reference by performing a suitable change of coordinates; second, find the solution to the simplified problem viewed in a rotating frame of reference by means of rotating coordinates; third, solve the original problem by an inverse transformation from the rotating to the stationary frame. Usually, the rotating frame of reference rotates about the axis specified by the magnetic field B with the angular frequency characterized by the magnetic field B . In quantum terms, assume that in a static frame of reference the quantum evolution of the qubit is ruled by the Schrödinger equation,
i ħ t ψ t = H t ψ t .
Then, considering a unitary transformation T between two distinct vector bases ψ t and ψ t , the Schrödinger equation in the new frame of reference can be written as
i ħ t ψ t = H t ψ t ,
with ψ t = def T ψ t and H t = def T H t T + i ħ t T T . More specifically, assuming that the unitary transformation T = def e i σ z ϕ ω t 2 is a rotation by an angle ϕ ω t about the z-axis, the Hamiltonians in Equations (2) and (3) can be recast as
H su 2 ; t = def Ω t + ħ 2 ϕ ˙ ω t σ z + ω H t σ x ,
and,
H su 2 ; t = def e ħ 2 m c B t + ħ 2 ϕ ˙ ω t σ z + e ħ 2 m c B t σ x ,
respectively, where ϕ ˙ ω = def d ϕ ω / d t . In the so-called original Rabi scenario, ϕ ω t = def ω 0 t with the angular frequency ω 0 being a negative real constant while the magnetic field intensities B and B are constant. In particular, assuming the so-called static resonance condition,
ω 0 = e m c B ,
the Hamiltonian in Equation (15) specifies a time-independent quantum mechanical problem. Unlike the original Rabi scenario, in the so-called generalized Rabi scenario proposed by Messina and collaborators in references [22,23], B , B , and ϕ ω are arbitrary time-dependent dynamical variables. In particular, assuming the so-called generalized Rabi condition,
ϕ ˙ ω t = 2 ħ Ω t ,
the Hamiltonians in Equations (14) and (15) do not yield a time-independent quantum mechanical problem. Inspired by our research reported in reference [7] and exploiting the findings in references [22,23], we take into consideration here four quantum mechanical scenarios where the transition probability P w w t from an initial state w to a final state w , with w | w = δ w , w , σ z w = + w , and σ z w = w , can be analytically described. In all four cases, we assume to be in a physical scenario in which,
ϕ ˙ ω t = ω 0 , and Ω t = ħ 2 ω 0 ,
with ω 0 being a negative real constant. We point out that, from a pure mathematical viewpoint, more general temporal behaviors of ϕ ˙ ω t and Ω t in Equation (18) could have been selected provided Equation (17) is fulfilled. However, the choice adopted in Equation (18) seems to be more suitable from an experimental standpoint. The four scenarios can be formally distinguished by means of the temporal behavior of the intensity ω H t of the complex transverse field ω t . In the first case, we assume a constant field intensity ω H t ,
ω H 1 t = def Γ .
This first case specifies the original Rabi scenario where P w w t is given by,
P w w 1 t = sin 2 Γ ħ t .
In the remaining three cases, we take into consideration three generalized Rabi scenarios where the field intensity ω H t manifests oscillatory, power law decay, and exponential law decay behaviors,
ω H 2 t = def Γ cos λ t , ω H 3 t = def Γ 1 + λ t 2 , and ω H 4 t = def Γ e λ t ,
respectively. Note that ω H 2 t in Equation (21) is a positive quantity on a temporal scale with 0 t π / 2 λ 1 . In all three cases, it can be proven that the transition probability P w w j t is given by [23],
P w w j t = sin 2 0 t ω H j t ħ d t ,
for any j 2 , 3 , 4 . Interestingly, being on resonance, the transition probabilities In Equations (20) and (22) for all four cases depend only on the integral of the transverse field intensity ω H t .
Before beginning our information geometric analysis, we formally introduce here an adimensional parameter β 0 that quantifies the departure from the on-resonance condition,
β 0 = def ħ 2 ω H t ϕ ˙ ω t + 2 ħ Ω t .
When β 0 assumes a nonzero constant value, the generalized expression of the transition probability in Equation (22) is given by [23],
P w w j t ; β 0 = def 1 1 + β 0 2 sin 2 1 + β 0 2 0 t ω H j t ħ d t .
We emphasize that β 0 in Equation (23) is generally a time-dependent quantity. However, in what follows, we shall limit our information geometric analysis to the case in which β 0 is a constant parameter with transition probabilities given in Equation (24). For further details on the link between β 0 in Equation (23) and its analogue counterpart in the framework of Rabi’s original static off-resonance condition, we refer to Appendix A.

3. Information Geometric Analysis

In this section, we present our information geometric analysis of off-resonance effects on the chosen quantum driving strategies. For each scheme, we analyze the consequences of a departure from the on-resonance condition in terms of both geodesic paths and geodesic speeds on the corresponding manifold of transition probability vectors.

3.1. Preliminaries

In order to provide some conceptual background for our forthcoming information geometric analysis that emerges from quantum driving strategies, we begin by introducing some preliminary remarks.
Assume to consider an N = def 2 n -dimensional complex Hilbert space H 2 n together with two neighboring normalized pure states ψ θ and ψ θ where,
ψ = def l = 1 N p l e i φ l l and ψ = def l = 1 N p l + d p l e i φ l + d φ l l ,
respectively. In Equation (25), l with 1 l N denotes an orthonormal basis of H 2 n while p l = p l θ and φ l = φ l θ are real functions of a continuous real parameter θ . The distinguishability metric on this manifold of Hilbert space rays is specified by the Fubini–Study metric [24],
d s FS 2 = def cos 1 ψ | ψ 2 = g a b FS θ d θ a d θ b ,
where the Fubini–Study metric tensor components g a b FS are related to the Fisher–Rao metric tensor components F a b θ by the condition,
g a b FS θ = 1 4 F a b θ + 4 σ a b 2 θ .
The quantities F a b θ and σ a b 2 θ in Equation (27) are given by,
F a b θ = def l = 1 N 1 p l θ p l θ θ a p l θ θ b ,
and,
σ a b 2 θ = def l = 1 N φ l θ θ a φ l θ θ b p l θ k = 1 N φ l θ θ a p l θ k = 1 N φ l θ θ b p l θ ,
respectively. In what follows, we assume that the non-negative term σ a b 2 θ in Equation (29) (which denotes the variance of the phase changes) is equal to zero. This working assumption can be justified by rephasing in a suitable fashion the basis vectors used to write the state ψ θ . Specifically, the rephasing procedure demands that Im ψ θ | l l | d ψ θ = 0 , for any 1 l N . The state d ψ = def d ψ ψ | d ψ d ψ denotes the projection of d ψ orthogonal to ψ where d ψ = def ψ ψ while ψ and ψ are given in Equation (25). In summary, for a convenient choice of the basis vectors l used for the decomposition in Equation (25), g a b FS θ becomes proportional to F a b θ as evident from Equation (27). We refer to reference [24] for further details. For the sake of completeness, we emphasize that if the basis l satisfies the above mentioned conditions, the basis vectors can always be rephased so that both l | ψ and l | d ψ are real. These latter conditions, in turn, are reminiscent of the so-called parallel transport conditions that occur in the context of Berry’s description of a geometric effect in the shape of an additional phase factor emerging after an adiabatic and cyclic transport of a quantum system [25]. Specifically, in the framework of Berry’s phase analysis generalized by Aharonov and Anandan by removing the adiabaticity constraint in the cyclic quantum evolution [26], it is possible to show that there exists one distinct curve in the Hilbert space H fulfilling the parallel transport conditions. Namely, two neighboring states ψ θ and ψ θ + d θ in H have the same phase, that is to say, ψ θ | ψ θ + d θ is real and positive. For such a curve, the dynamical phase vanishes. In summary, for special Hamiltonians that satisfy a special gauge choice (that is, the parallel transport conditions), the dynamical phase vanishes [27].
In our paper, the parameter θ is the statistical version of the elapsed time t. That is, we assume θ is a parameter that can be experimentally specified by measurement of a suitable observable that changes with time (for instance, the transverse magnetic field intensity B t ). For further details on the concept of statistical elapsed time, we refer to reference [28]. In particular, we suppose that the output ψ θ of a quantum driving Hamiltonian acting on the input defined by the normalized source state s = def ψ θ 0 can be described as,
ψ θ 0 ψ θ = def e i φ w θ p w θ w + e i φ w θ p w θ w .
In general, the normalized output state ψ θ is an element of the two-dimensional subspace of the n-qubit complex Hilbert space H 2 n accommodating the source state s and generated by the set of orthonormal state vectors w , w . The modulus squared of the probability amplitudes p w θ = def w | ψ θ 2 and p w θ = def w | ψ θ 2 stand for the success and failure probabilities of the given driving strategy, respectively. Furthermore, φ w θ and φ w θ denote the real quantum phases of the states w and w , respectively. We note that the quantum state ψ θ in Equation (30) is parameterized by means of a single continuous real parameter that emerges from the elapsed time of the driving strategy. As previously mentioned, this parameter θ mimics a statistical macrovariable used to distinguish neighboring quantum states ψ θ and ψ θ + d ψ θ along a path through the space of pure quantum states. In summary, we focus on the space of probability distributions p θ ,
ψ θ p θ = def p w θ , p w θ = w | ψ θ 2 , w | ψ θ 2 ,
where the natural Riemannian distinguishability metric between two neighboring probability distributions p θ and p θ + d θ is given by the Fisher–Rao information metric in Equation (28). For an overview of the use of information geometric techniques to analog quantum search problems, we refer to references [20,29,30,31].
In any quantum driving scheme considered in this paper, following references [22,23], we assume the success probability p 0 θ and the failure probability p 1 θ can be recast as
p 0 θ = def A sin 2 Σ θ and , p 1 θ = def 1 A sin 2 Σ θ ,
respectively. In Equation (32), the quantity θ denotes the temporal parameter while the quantity A is defined as
A = A β 0 = def 1 1 + β 0 2 ,
with β 0 being the (assumed) constant parameter that quantifies the deviation from the on-resonance condition. The quantity Σ θ , instead, specifies the type of driving being considered. The Fisher information F θ corresponding to the probability vector p 0 θ , p 1 θ is given by [32,33],
F θ = def p ˙ 0 2 θ p 0 θ + p ˙ 1 2 θ p 1 θ ,
with p ˙ = def d p / d θ . Substituting Equation (32) into Equation (34), the Fisher information F θ becomes
F θ = def 4 A 2 d Σ d θ 2 cos 2 Σ θ A + cos 2 Σ θ sin 2 Σ θ 1 A sin 2 Σ θ .
Note that when we are on-resonance, β 0 = 0 , A = 1 , and F θ reduces to
F θ = def 4 d Σ d θ 2 .
In particular, when Σ θ is linear in θ , F θ = 4 as reported in reference [20]. Finally, once the expression of the Fisher information F θ is given, the geodesic equation to integrate becomes
d 2 θ d ξ 2 + 1 2 F θ d F θ d θ d θ d ξ 2 = 0 .
For more details on how to derive Equation (37) and how to express the quantum mechanical infinitesimal Fubini–Study line element in terms of the Fisher information and/or the Fisher–Rao information metric, we refer to reference [20]. For the sake of completeness, we observe that Equation (37) can be recast as,
θ ¨ + F ˙ 2 F θ ˙ = 0 ,
where θ ˙ = def d θ / d ξ and F ˙ = def d F / d ξ . Defining w = def θ ˙ and after some straightforward algebraic manipulations, Equation (38) reduces to
d F F = 2 d w w .
Integration of Equation (39) yields a closed-form implicit relation between the dependent variable θ and the independent variable ξ ,
c θ , ξ = θ ξ 0 θ ξ F θ d θ ξ 0 ξ a 0 d ξ = 0 ,
with a 0 being a real constant of integration. Equation (40) is an implicit constraint equation between the variables θ and ξ . Despite the closed-form implicit relation in Equation (40), in what follows we shall numerically find the geodesic paths satisfying Equation (37) for the four driving strategies defined in Equations (19) and (21).

3.2. Geodesic Paths

In what follows, after specifying Σ θ in Equation (32) and F θ in Equation (35), we shall focus on numerically finding from Equation (37) the geodesic paths θ = θ ξ that correspond to the various off-resonant driving schemes.

3.2.1. Constant Behavior

In the first case, we assume Σ θ is defined as
Σ θ = def B θ .
Therefore, from Equations (32) and (41), the probabilities p 0 θ and p 1 θ are given by
p 0 θ = def A sin 2 B θ and , p 1 θ = def 1 A sin 2 B θ ,
respectively. In Equation (42), the quantity A is defined in Equation (33) while B is given by
B = B β 0 , Γ , ħ = def 1 + β 0 2 Γ ħ ,
respectively. In Figure 1, we plot the success probability p 0 θ versus the temporal parameter θ in the case of constant driving with the assumption that β 0 0 , 1 / 4 , 1 / 2 , 1 . From Figure 1, we note that the maximum value of the success probability decreases with respect to increasing values of the parameter β 0 . Furthermore, the periodic oscillatory behavior of the success probability is also affected by nonzero values of β 0 . In particular, to higher values of β 0 there correspond smaller (higher) values of the period (frequency) of oscillation. In this first case, putting Equation (42) into Equation (35), the Fisher information becomes
F constant off resonance θ = def 4 AB 2 1 + cos 2 B θ 2 A 1 cos 2 B θ .
For the sake of completeness, note that when β 0 = 0 , A = 1 , and F θ in Equation (44) reduces to
F constant on resonance θ = 4 B 2 .
Finally, employing Equations (37) and (44), the geodesic equation to integrate becomes
d 2 θ d ξ 2 1 A B tan B θ 1 A sin 2 B θ d θ d ξ 2 = 0 .
The output arising from the numerical integration of Equation (46) is plotted in Figure 2.

3.2.2. Oscillatory Behavior

In this second scenario, we assume Σ θ is given by
Σ θ = def B λ sin λ θ .
Therefore, from Equations (32) and (47), the probabilities p 0 θ and p 1 θ become
p 0 θ = def A sin 2 B λ sin λ θ and , p 1 θ = def 1 A sin 2 B λ sin λ θ ,
respectively. In this case, employing Equation (47), the Fisher information in Equation (35) becomes
F oscillatory off resonance θ = def d A sin 2 B λ sin λ θ d θ 2 A sin 2 B λ sin λ θ + d 1 A sin 2 B λ sin λ θ d x 2 1 A sin 2 B λ sin λ θ ,
that is, after some straightforward but tedious algebra,
F oscillatory off resonance θ = 8 AB 2 cos 2 λ θ cos 2 B λ sin λ θ 2 A 1 cos 2 B λ sin λ θ .
Observe that when β 0 = 0 , A = 1 , and, recalling that 1 + cos 2 x = 2 cos 2 x , F θ in Equation (50) reduces to
F oscillatory on resonance θ = 4 B 2 cos 2 ( λ θ ) .
Finally, employing Equations (37) and (50), the geodesic equation to integrate in this second scenario can be formally written as
d 2 θ d ξ 2 + 1 2 F o θ d F o θ d θ d θ d ξ 2 = 0 ,
with F o θ = def F oscillatory off resonance θ as in Equation (50). Using a symbolic mathematical software (Mathematica, for instance) together with the analytical knowledge of what happens in the limiting case of β 0 = 0 [21], Equation (52) can be finally recast as
d 2 θ d ξ 2 λ tan λ θ + 1 A B cos λ θ tan B λ sin λ θ 1 A sin 2 B λ sin λ θ d θ d ξ 2 = 0 .
The output emerging from numerically integrating Equation (53) is plotted in Figure 2.

3.2.3. Power Law Decay

In this third scenario, we assume Σ θ is defined as
Σ θ = def B λ 1 1 1 + λ θ .
Therefore, the probabilities p 0 θ and p 1 θ in Equation (32) can be rewritten as
p 0 θ = def A sin 2 B λ 1 1 1 + λ θ and , p 1 θ = def 1 A sin 2 B λ 1 1 1 + λ θ ,
respectively. In this case, the Fisher information in Equation (35) can be formally recast as
F power law decay off resonance θ = def d A sin 2 B λ 1 1 1 + λ θ d θ 2 A sin 2 B λ 1 1 1 + λ θ + d 1 A sin 2 B λ 1 1 1 + λ θ d x 2 1 A sin 2 B λ 1 1 1 + λ θ ,
that is, after some tedious algebra,
F power law decay off resonance θ = 4 AB 2 1 + λ θ 4 1 + cos 2 B θ 1 + λ θ 2 A 1 cos 2 B θ 1 + λ θ .
Once again, for the sake of completeness, we point out that when β 0 = 0 , A = 1 , and F θ in Equation (57) reduces to
F power law decay on resonance θ = 4 B 2 1 + λ θ 4 .
Finally, making use of Equations (37) and (57), the geodesic equation to integrate in this third scenario can be formally written as
d 2 θ d ξ 2 + 1 2 F pld θ d F pld θ d θ d θ d ξ 2 = 0 ,
with F pld θ = def F power law decay off resonance θ as in Equation (57). Employing a symbolic mathematical software (Mathematica, for instance) together with the analytical knowledge of what occurs in the on-resonance case where β 0 = 0 [21], Equation (59) can be finally recast as
d 2 θ d ξ 2 2 λ 1 + λ θ + ( 1 A ) B tan B θ 1 + λ θ 1 + λ θ 2 1 A sin 2 B θ 1 + λ θ d θ d ξ 2 = 0 .
In Figure 2, we plot the output that arises from the numerical integration of Equation (60).

3.2.4. Exponential Decay

In this last scenario, we assume
Σ θ = def B λ 1 e λ θ .
Therefore, the success and failure probabilities probabilities p 0 θ and p 1 θ are given by
p 0 θ = def A sin 2 B λ 1 e λ θ and , p 1 θ = def 1 A sin 2 B λ 1 e λ θ ,
respectively. In this case, the Fisher information in Equation (35) is formally given by
F exponential decay off resonance θ = def d A sin 2 B λ 1 e λ θ d θ 2 A sin 2 B λ 1 e λ θ + d 1 A sin 2 B λ 1 e λ θ d x 2 1 A sin 2 B λ 1 e λ θ ,
that is,
F exponential decay off resonance θ = def 4 AB 2 e 2 λ θ 1 + cos 2 B λ 1 e λ θ 2 A 1 cos 2 B λ 1 e λ θ
Observe that when β 0 = 0 , A = 1 , F θ in Equation (64) reduces to
F exponential decay on resonance θ = 4 B 2 e 2 λ θ .
Lastly, using Equations (37) and (64), the geodesic equation to integrate in this third scenario can be formally written as
d 2 θ d ξ 2 + 1 2 F exp θ d F exp θ d θ d θ d ξ 2 = 0 ,
with F exp θ = def F exponential decay off resonance θ as in Equation (64). Making use of a symbolic mathematical software (Mathematica, for instance) together with the analytical knowledge of what happens in the limiting case of β 0 = 0 [21], Equation (66) can be finally recast as
d 2 θ d ξ 2 λ + 1 A B e λ θ tan B λ 1 e λ θ 1 A sin 2 B λ 1 e λ θ d θ d ξ 2 = 0 .
In Figure 2, we report the numerical plots of the geodesic paths θ versus the affine parameter ξ for the four driving strategies with geodesic equations in Equations (46), (53), (60) and (67). In the LHS of Figure 2, we consider the on-resonance scenario with β 0 = 0 . In the RHS of Figure 2, instead, we consider the off-resonance scenario with β 0 = 1 / 2 . In both sides, we set Γ / ħ = 1 , λ = 2 / π , θ 0 = 0 , and θ ˙ 0 = 1 . The most noticeable feature in Figure 2 is the deviation from the straight line behavior of the geodesic path that emerges from the constant driving scheme when the on-resonance condition is satisfied. In particular, a first visual comparison of the corresponding plots in the two cases (on-resonance v.s. off-resonance) seems to suggest that the constant driving scheme might be the least robust to departures from the on-resonance condition. Indeed, this preliminary remark will be quantitatively confirmed in our forthcoming information geometric analysis.

3.3. Geodesic Speeds

In what follows, after specifying Σ θ in Equation (32) and F θ in Equation (35), we shall focus on finding the geodesic speeds v,
v = def 1 2 F θ d θ d ξ ,
corresponding to the various off-resonant driving schemes. For some technical details on the proof of the fact that geodesic paths have constant speed, we refer to Appendix B.

3.3.1. Constant Case

In this first case, using Equations (33), (43) and (44), the on-resonance geodesic speed v from Equation (68) is given by
v ON = Γ ħ θ ˙ .
In the off-resonance regime, instead, the geodesic speed v from Equation (68) can be recast as
v OFF = Γ ħ 1 + cos 2 Γ ħ 1 + β 0 2 θ 2 1 1 + β 0 2 1 cos 2 Γ ħ 1 + β 0 2 θ 1 2 θ ˙ .
Furthermore, setting Γ / ħ = 1 , θ ˙ ξ 0 = θ ˙ 0 = 1 , and θ ξ 0 = θ 0 , v ON and v OFF in Equations (69) and (70) become
v ON constant = def 1 ,
and,
v OFF constant β 0 , θ 0 = def 1 + cos 2 1 + β 0 2 θ 0 2 1 1 + β 0 2 1 cos 2 1 + β 0 2 θ 0 1 2 ,
respectively. From Equations (71) and (72), we note that the geodesic speed becomes sensitive to the initial condition θ 0 when departing from the on-resonance scenario in the case of a driving strategy with constant complex transverse field intensity. Finally, employing Equations (72) and (71), we define the robustness coefficient as the ratio
r constant β 0 , θ 0 = def v OFF constant β 0 , θ 0 v ON constant θ 0 ,
that is,
r constant β 0 , θ 0 = 1 + cos 2 1 + β 0 2 θ 0 2 1 1 + β 0 2 1 cos 2 1 + β 0 2 θ 0 1 2 .
In Figure 3, we plot the analytical (thin solid) and numerical (filled circle) values of the geodesic speed v OFF in Equation (70) versus the initial condition θ 0 in the case of constant driving. In the plot, we set β 0 = 1 / 2 , Γ / ħ = 1 , and θ ˙ 0 = 1 . The plot in Figure 3 clearly illustrates the emergence of the sensitivity on the initial condition θ 0 of the geodesic speed v OFF for the constant driving scheme when the on-resonance condition is violated. For the sake of completeness, we point out that despite the fact that this sensitivity is present in the remaining three driving strategies even when the on-resonance condition is fulfilled, off-resonance effects substantially enhance this sensitivity as evident from our information geometric analysis of geodesic speeds presented here.

3.3.2. Oscillatory Behavior

In this second case, making use of Equations (33), (43) and (50), the on-resonance geodesic speed v from Equation (68) becomes
v ON = Γ ħ cos λ θ θ ˙ .
In the off-resonance regime, instead, the geodesic speed v OFF is given by
v OFF = Γ ħ cos λ θ 2 cos 2 1 λ Γ ħ 1 + β 0 2 sin λ θ 2 1 1 + β 0 2 1 cos 2 λ Γ ħ 1 + β 0 2 sin λ θ 1 2 θ ˙ .
Moreover, letting Γ / ħ = 1 , λ = 2 / π , θ ˙ ξ 0 = θ ˙ 0 = 1 , and θ ξ 0 = θ 0 , v ON and v OFF in Equations (75) and (76) become
v ON oscillatory = def cos 2 π θ 0 ,
and,
v OFF oscillatory β 0 , θ 0 = def cos 2 π θ 0 2 cos 2 π 2 1 + β 0 2 sin 2 π θ 0 2 1 1 + β 0 2 1 cos π 1 + β 0 2 sin 2 π θ 0 1 2 ,
respectively. Finally, using Equations (77) and (78), we introduce the robustness coefficient as the ratio
r oscillatory β 0 , θ 0 = def v OFF oscillatory β 0 , θ 0 v ON oscillatory θ 0 ,
that is,
r oscillatory β 0 , θ 0 = 2 cos 2 π 2 1 + β 0 2 sin 2 π θ 0 2 1 1 + β 0 2 1 cos π 1 + β 0 2 sin 2 π θ 0 1 2 .
From Equations (77) and (78), we observe that v OFF oscillatory β 0 , θ 0 v ON oscillatory and 0 r oscillatory β 0 , θ 0 < 1 .

3.3.3. Power Law Decay

In this third case, Equations (33), (43) and (57) lead to an expression of the on-resonance geodesic speed v from Equation (68) given by
v ON = Γ ħ 1 1 + λ θ 2 θ ˙ .
In the off-resonance regime where β 0 0 , instead, the geodesic speed v OFF from Equation (68) becomes
v OFF = Γ ħ 1 1 + λ θ 2 1 + cos 2 Γ ħ 1 + β 0 2 θ 1 + λ θ 2 1 1 + β 0 2 1 cos 2 Γ ħ 1 + β 0 2 θ 1 + λ θ 1 2 θ ˙ .
In addition, assuming Γ / ħ = 1 , λ = 2 / π , θ ˙ ξ 0 = θ ˙ 0 = 1 , and θ ξ 0 = θ 0 , v ON and v OFF in Equations (81) and (82) become
v ON power law decay = def 1 1 + 2 π θ 0 2 ,
and,
v OFF power law decay β 0 , θ 0 = def 1 1 + 2 π θ 0 2 1 + cos 2 1 + β 0 2 θ 0 1 + 2 π θ 0 2 1 1 + β 0 2 1 cos 2 1 + β 0 2 θ 0 1 + 2 π θ 0 1 2 ,
respectively. Lastly, employing Equations (83) and (84), we define the robustness coefficient as the ratio
r power law decay β 0 , θ 0 = def v OFF power law decay β 0 , θ 0 v ON power law decay θ 0 ,
that is,
r power law decay β 0 , θ 0 = 1 + cos 2 1 + β 0 2 θ 0 1 + 2 π θ 0 2 1 1 + β 0 2 1 cos 2 1 + β 0 2 θ 0 1 + 2 π θ 0 1 2 .
From Equations (83) and (84), we observe that v OFF power law decay β 0 , θ 0 v ON power law decay and 0 r power law decay β 0 , θ 0 < 1 .

3.3.4. Exponential Decay

In this fourth and last case, making use of Equations (33), (43) and (64), the on-resonance geodesic speed v from Equation (68) becomes
v ON = Γ ħ e λ θ θ ˙ .
In the off-resonance regime, instead, the geodesic speed v OFF is given by
v OFF = Γ ħ e λ θ 1 + cos 2 λ Γ ħ 1 + β 0 2 1 e λ θ 2 1 1 + β 0 2 1 cos 2 λ Γ ħ 1 + β 0 2 1 e λ θ 1 2 θ ˙ .
Assuming Γ / ħ = 1 , λ = 2 / ħ , θ ˙ ξ 0 = θ ˙ 0 = 1 , and θ ξ 0 = θ 0 , v ON and v OFF in Equations (87) and (88) become
v ON exponential decay = def e 2 π θ 0 ,
and,
v OFF exponential decay β 0 , θ 0 = def e 2 π θ 0 1 + cos π 1 + β 0 2 1 e 2 π θ 0 2 1 1 + β 0 2 1 cos π 1 + β 0 2 1 e 2 π θ 0 1 2 ,
respectively. Finally, by means of Equations (89) and (90), we define the robustness coefficient as the ratio
r exponential decay β 0 , θ 0 = def v OFF exponential decay β 0 , θ 0 v ON exponential decay θ 0 ,
that is,
r exponential decay β 0 , θ 0 = 1 + cos π 1 + β 0 2 1 e 2 π θ 0 2 1 1 + β 0 2 1 cos π 1 + β 0 2 1 e 2 π θ 0 1 2 .
In Figure 4, going from left to right, we plot: (L) Plot of the geodesic speeds v ON in Equations (69), (75), (81) and (87) versus the initial condition θ 0 with β 0 = 0 ; (C) Plot of the geodesic speeds v OFF in Equations (70), (76), (82) and (88) versus the initial condition θ 0 with β 0 = 1 / 2 ; (R) Plot of the robustness coefficients r in Equations (73), (79), (85) and (91) versus the affine parameter ξ . Plots that correspond to the constant, oscillatory, exponential decay, and power law decay driving strategies appear in dotted, dashed, thin solid, and thick solid lines, respectively. In all plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 . The plots in Figure 4 help exhibiting the effects of the off-resonance condition on the speed of geodesic paths that correspond to the variety of driving strategies being considered in this paper. In particular, we observe that to faster (slower) driving strategies, there seem to correspond smaller (larger) robustness coefficients. In other words, slower driving strategies appear to be more robust against departures from the on-resonance condition. The ranking of the various driving strategies is quite straightforward in the on-resonance condition. However, when departing from this condition, off-resonance effects lead to a richer set of dynamical scenarios. This, in turn, makes the comparison of the performance of our chosen driving schemes more delicate. For example, despite our illustrative depiction in Figure 4, it is possible to uncover two-dimensional parametric regions P β 0 , θ 0 specified by the parameters β 0 and θ 0 where the constant driving scheme is not only the least robust but also the slowest one. Indeed, we plot in Figure 5 a two-dimensional parametric region P β 0 , θ 0 ,
P β 0 , θ 0 = def β 0 , θ 0 R + × R + : v OFF strategy i v OFF strategy j and , r strategy i r strategy j ,
with i j 1 , 2 , 3 , 4 and where the oscillatory (light grey), the exponential decay (grey), and the power law decay (black) driving strategies outperform the constant driving scheme in terms of both geodesic speed and robustness. More specifically, in Figure 5 we have
P power law decay β 0 , θ 0 P exponential decay β 0 , θ 0 P oscillatory β 0 , θ 0 ,
and, in addition,
r power law decay β 0 , θ 0 r exponential decay β 0 , θ 0 r oscillatory β 0 , θ 0 .
In all region plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 . In addition, the maximum success probability of each driving scheme is assumed to be greater than 25 / 26 , that is, 0 β 0 0.20 in Figure 5. The main take-home message from Figure 5 is that in the off-resonance scenario it is possible to uncover two-dimensional parametric regions P β 0 , θ 0 where the best on-resonance driving scheme (that is, the constant one) can be outperformed both in terms of speed and robustness.

4. Conclusions

In this paper, we presented an information geometric analysis of off-resonance effects on classes of exactly solvable generalized semi-classical Rabi systems. Specifically, we considered population transfer performed by four distinct off-resonant driving schemes specified by su 2 ; time-dependent Hamiltonian models. For each scheme, we studied the consequences of a departure from the on-resonance condition in terms of both geodesic paths and geodesic speeds on the corresponding manifold of transition probability vectors. In particular, we analyzed the robustness of each driving scheme against off-resonance effects. Moreover, we reported on a possible tradeoff between speed and robustness in the driving schemes being investigated. Finally, we discussed the emergence of a different relative ranking in terms of performance among the various driving schemes when transitioning from the on-resonant to the off-resonant scenarios.
Our main findings can be outlined as follows.
(1)
In the presence of off-resonance effects, the success probability of the various quantum strategies is affected both in terms of amplitude and periodicity. In particular, focusing on the constant driving case, the success probability p 0 θ in Equation (42) is dampened by the Lorentzian-like factor 1 / ( 1 + β 0 2 ) . Moreover, the periodicity of the oscillations of this probability changes from 2 π to 2 π / 1 + β 0 2 1 2 . In summary, oscillations become smaller in amplitude but higher in frequency. These facts are clearly visible in Figure 1.
(2)
Departing from the on-resonance condition (that is, β 0 = 0 ), we observe a change in the geodesic paths on the underlying manifolds. The presence of off-resonance effects (that is, β 0 0 ) generates geodesic paths with a more complex structure. In particular, unlike what happens when the on-resonance condition is satisfied, numerical integration of the geodesic equations is required. In particular, the numerical plots of the geodesic paths θ versus the affine parameter ξ for the four driving strategies with geodesic equations in Equations (46), (53), (60) and (67) appear in Figure 2.
(3)
Each and every numerical value of the geodesic speeds corresponding to the various quantum driving strategies considered in this paper becomes smaller when departing from the on-resonance condition (see Figure 3, for instance). In general, quantum driving strategies characterized by a high geodesic speed appear to be very sensitive to off-resonance effects. Instead, strategies yielding geodesic paths with smaller geodesic speed values seem to be more robust against departures from the on-resonance condition. These observations can be understood from Figure 4.
(4)
In the off-resonance regime, there emerges a sensitive dependence of the geodesic speeds on the initial conditions. This, in turn, can cause a change in the ranking of the various quantum driving strategies. In particular, the strategy specified by a constant Fisher information is no longer the absolute best strategy in terms of speed. Indeed, it is possible to find two-dimensional (2D) parametric regions where this strategy is being outperformed by all the remaining strategies both in terms of speed and robustness. In general, these 2D regions are more extended for the less robust strategies (see Equations (94) and (95)). For instance, we can identify 2D regions where the speed-based ranking becomes: (1) Oscillatory strategy; (2) exponential-law decay strategy; (3) power-law decay strategy. Instead, considering the very same 2D regions, the robustness-based ranking is given by: (1) Power-law decay strategy; (2) exponential-law decay strategy; (3) oscillatory strategy. In these 2D regions, the constant strategy is the worst, both in terms of speed and robustness. These findings are illustrated in Figure 5.
In conclusion, our information geometric analysis suggests that speed and robustness appear to be competing features when departing from the on-resonance condition. Ideally, a quantum driving strategy should be fast, robust, and thermodynamically efficient [34,35,36]. Since speed and thermodynamic efficiency quantified in terms of minimum entropy production paths seem to be conflicting properties as well [20,21], it appears reasonable to think there might be a connection between robustness and thermodynamic efficiency. We leave the exploration of this conjectured link to future scientific efforts.

Author Contributions

Conceptualization, C.C.; Numeri- cal computations, S.G. and C.C.; Writing-original draft, C.C.; Writing-review and editing, C.C., S.G., and P.M.A.; Supervision, P.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

C.C. acknowledges support under AFRL Visiting Faculty Fellowship Program, AFRL Contract No. FA8750-15-3-6003.

Acknowledgments

C.C. is grateful to the United States Air Force Research Laboratory (AFRL) Summer Faculty Fellowship Program for providing support for this work. Any opinions, findings and conclusions or recommendations expressed in this manuscript are those of the authors and do not necessarily reflect the views of AFRL.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Violating a Resonance Condition

In this Appendix, we provide two illustrative examples concerning the violation of a resonance condition.

Appendix A.1. A Classical Scenario

In the framework of classical mechanics, one can imagine violating a resonance condition in a number of manners. For instance, consider a mass-spring system in the presence of damping and sinusoidal forcing term,
m x ¨ + λ x ˙ + k x = f 0 cos γ t .
In Equation (A1), x ˙ = def d x / d t . Furthermore, m, λ , k, and f 0 denote the mass, the damping coefficient, the spring constant, and the amplitude of the forcing term, respectively. The classical resonance curve C γ ,
C γ = def k m γ 2 2 + b 2 γ 2 1 / 2 ,
for this specific physical system is proportional to the amplitude A γ = def f 0 C γ of the steady-state solution to Equation (A1). Moreover, the resonance condition γ = γ * with
γ * = γ * m , k , λ = def k m λ 2 2 m 2 1 / 2 ,
is specified by imposing the maximum C ¯ = def C γ * of the resonance curve C γ in Equation (A2). Therefore, we clearly note from Equation (A3) that departures from the resonance condition γ = γ * can happen by varying the stiffness (k), the mass (m), and/or the damping ( λ ).

Appendix A.2. A Quantum Scenario

In the framework of quantum mechanics, one can imagine violating a resonance condition in a number of manners. For instance, consider a two-level quantum system described by a time-dependent Hamiltonian,
H t = def H 0 + V t ,
where H 0 is the part of the Hamiltonian H that does not contain time explicitly while V t is the time-dependent sinusoidal oscillatory potential. The quantities H 0 and V t are defined as,
H 0 = def E 1 1 1 + E 2 2 2 , V t = def γ e i ω t 1 2 + γ e i ω t 2 1 ,
respectively, where γ and ω belong to R + 0 . Furthermore, 1 and 2 are two eigenstates of H 0 with corresponding eigenvalues E 1 and E 2 , respectively, with E 2 > E 1 . The quantum resonance curve Q ω ,
Q ω = def 1 + ħ γ 2 ω ω 21 2 4 1 ,
for this specific physical system described by the Hamiltonian in Equation (A4) is the amplitude squared that characterizes the transition probability P 1 2 (assuming that at t = 0 only the level 1 is populated) between the two quantum states 1 and 2 with ω 21 = def E 2 E 1 / ħ being the characteristic frequency of the system. Moreover, the resonance condition ω = ω * with
ω * = ω 21 = ω 21 E 1 , E 2 = def E 2 E 1 / ħ ,
is specified by imposing the maximum Q ¯ = def Q ω 21 of the resonance curve Q ω in Equation (A6). Therefore, we clearly note from Equation (A7) that departures from the resonance condition ω = ω * can happen by varying the larger ( E 2 ) and the smaller ( E 1 ) energy levels. We point out that considering the following correspondences
ω ω , γ e ħ B 2 m c , and ω 21 e B m c ,
it can be shown that the Hamiltonian in Equation (A4) describes a two-level quantum mechanical system represented by a spin- 1 / 2 particle immersed in an external magnetic field B t given by,
B t = def B z ^ + B cos ω t x ^ + sin ω t y ^ .
The resonance condition ω = ω * is satisfied whenever the positively assumed frequency ω of the rotating magnetic field in the x y -plane equals Larmor’s precessional frequency specified by the intensity of the uniform magnetic field B z ^ [2],
ω * = ω * e , B , m , c = def e B m c .
Finally, we emphasize that the analogue of the adimensional parameter β 0 used throughout our paper becomes in this quantum mechanical framework the quantity defined as
β 0 = def m c e B ω ω * = m c e B ω e B m c .
Observe that Equation (A11) is the static version of Equation (23). For further details on the quantum mechanical motion in two-level quantum systems described by time-dependent Hamiltonians, we refer to reference [2].

Appendix B. Speed of Geodesics

In the Appendix, we briefly show that geodesics have constant speed. For further details, we refer to reference [37].
Recall that if Θ ξ = def θ 1 ξ , . . . , θ n ξ is a curve in a n-dimensional Riemannian manifold M , the speed of Θ at any time is the length of its velocity vector Θ ˙ ξ ,
Θ ˙ ξ = def Θ ˙ ξ , Θ ˙ ξ 1 2 = g a b θ θ ˙ a θ ˙ b 1 2 ,
where g a b θ denotes a Riemannian metric on the manifold M while θ a ξ are the component functions of the curve Θ ξ . We say Θ ξ is constant speed if Θ ˙ ξ in Equation (A12) does not depend on ξ , and unit speed if the speed is identically equal to one.
It is well-known that all Riemannian geodesics are constant speed curves. Indeed, let M be a manifold with a linear connection , and let Θ be a curve on M . The acceleration of Θ is the vector field D ξ Θ ˙ along Θ (that is, the covariant derivative of Θ ˙ ξ along Θ ). A curve Θ is called a geodesic with respect to if its acceleration D ξ Θ ˙ equals zero (that is, the velocity vector field Θ ˙ is parallel along the curve Θ ). Therefore, for any geodesic Θ ξ , we have
d d ξ Θ ˙ ξ 2 = d d ξ Θ ˙ ξ , Θ ˙ ξ = D ξ Θ ˙ ξ , Θ ˙ ξ + Θ ˙ ξ , D ξ Θ ˙ ξ = 0 ,
that is,
d d ξ Θ ˙ ξ 2 = 0 .
From Equation (A14), we conclude that geodesics Θ ξ have constant speed Θ ˙ ξ .

References

  1. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons, Inc.: New York, NY, USA, 1999. [Google Scholar]
  2. Sakurai, J.J. Modern Quantum Mechanics; Addison-Wesley Publishing Company, Inc.: Boston, MA, USA, 1994. [Google Scholar]
  3. D’Alessandro, D.; Dahleh, M. Optimal control of two-level quantum systems. IEEE Trans. Autom. Control 2001, 46, 866. [Google Scholar] [CrossRef]
  4. Boscain, U.; Charlot, G.; Gauthier, J.-P.; Guerin, S.; Jauslin, H.-R. Optimal control in laser-induced population transfer for two- and three-level quantum systems. J. Math. Phys. 2002, 43, 2107. [Google Scholar] [CrossRef]
  5. Romano, R. Geometric analysis of minimum-time trajectories for a two-level quantum system. Phys. Rev. A 2014, 90, 062302. [Google Scholar] [CrossRef] [Green Version]
  6. Byrnes, T.; Forster, G.; Tessler, L. Generalized Grover’s algorithm for multiple phase inversion states. Phys. Rev. Lett. 2018, 120, 060501. [Google Scholar] [CrossRef] [Green Version]
  7. Cafaro, C.; Alsing, P.M. Continuous-time quantum search and time-dependent two-level quantum systems. Int. J. Quantum Inf. 2019, 17, 1950025. [Google Scholar] [CrossRef]
  8. Aharonov, Y.; Anandan, J.; Popescu, S.; Vaidman, L. Superpositions of time evolutions of a quantum system and a quantum time-translation machine. Phys. Rev. Lett. 1990, 64, 2965. [Google Scholar] [CrossRef] [Green Version]
  9. Kempf, A.; Prain, A. Driving quantum systems with Superoscillations. J. Math. Phys. 2017, 58, 082101. [Google Scholar] [CrossRef]
  10. Barredo, D.; Labuhn, H.; Ravets, S.; Lahaye, T.; Browaeys, A. Coherent excitation transfer in a spin chain of three Rydberg atoms. Phys. Rev. Lett. 2015, 114, 113002. [Google Scholar] [CrossRef] [Green Version]
  11. Schempp, H.; Gunter, G.; Wuster, S.; Weidemuller, M.; Whitlock, S. Correlated exciton transport in Rydberg-dressed-atom spin chains. Phys. Rev. Lett. 2015, 115, 093002. [Google Scholar] [CrossRef]
  12. Shi, Z.C.; Wang, W.; Yi, X.X. Population transfer driven by off-resonant fields. Opt. Express 2016, 24, 21971. [Google Scholar] [CrossRef] [Green Version]
  13. Hirose, M.; Cappellaro, P. Time-optimal control with finite bandwidth. Quantum Inf. Process. 2018, 17, 88. [Google Scholar] [CrossRef] [Green Version]
  14. Allen, L.; Eberly, J.H. Optical Resonance and Two-Level Atoms; Dover Publications: New York, NY, USA, 1987. [Google Scholar]
  15. Laucht, A.; Simmons, S.; Kalra, R.; Tosi, G.; Dehollain, J.P.; Muhonen, J.T.; Freer, S.; Hudson, F.E.; Itoh, K.M.; Jamieson, D.N.; et al. Breaking the rotating wave approximation for a strongly driven dressed single-electron spin. Phys. Rev. B 2016, 94, 161302. [Google Scholar] [CrossRef] [Green Version]
  16. Dai, K.; Wu, H.; Zhao, P.; Li, M.; Liu, Q.; Xue, G.; Tan, X.; Yu, H.; Yu, Y. Quantum simulation of the general semi-classical Rabi model in regimes of arbitrarily strong driving. Appl. Phys. Lett. 2017, 111, 242601. [Google Scholar] [CrossRef]
  17. Cafaro, C.; Alsing, P.M. Theoretical analysis of a nearly optimal analog quantum search. Phys. Scr. 2019, 94, 085103. [Google Scholar] [CrossRef] [Green Version]
  18. Gassner, S.; Cafaro, C.; Capozziello, S. Transition probabilities in generalized quantum search Hamiltonian evolutions. Int. J. Geom. Methods Mod. Phys. 2020. accepted. [Google Scholar] [CrossRef]
  19. Farhi, E.; Gutmann, S. Analog analogue of a digital quantum computation. Phys. Rev. A 1998, 57, 2403. [Google Scholar] [CrossRef] [Green Version]
  20. Cafaro, C.; Alsing, P.M. Decrease of Fisher information and the information geometry of evolution equations for quantum mechanical probability amplitudes. Phys. Rev. E 2018, 97, 042110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Cafaro, C.; Alsing, P.M. Information geometry aspects of minimum entropy production paths from quantum mechanical evolutions. Phys. Rev. E 2019, in press. [Google Scholar]
  22. Messina, A.; Nakazato, H. Analytically solvable Hamiltonians for quantum two-level systems and their dynamics. J. Phys. A Math. Theor. 2014, 47, 445302. [Google Scholar] [CrossRef] [Green Version]
  23. Grimaudo, R.; de Castro, A.S.M.; Nakazato, H.; Messina, A. Classes of exactly solvable generalized semi-classical Rabi systems. Ann. Phys. 2018, 2018, 1800198. [Google Scholar] [CrossRef] [Green Version]
  24. Braunstein, S.L.; Caves, C.M. Statistical distance and geometry of quantum states. Phys. Rev. Lett. 1994, 72, 3439. [Google Scholar] [CrossRef] [PubMed]
  25. Berry, M.V. Quantal phase factors accompanying adiabatic changes. Proc. R. Soc. Lond. A 1984, 392, 45. [Google Scholar]
  26. Aharonov, Y.; Anandan, J. Phase change during a cyclic quantum evolution. Phys. Rev. Lett. 1987, 58, 1593. [Google Scholar] [CrossRef] [PubMed]
  27. Resta, R. Manifestations of Berry’s phase in molecules and condensed matter. J. Phys. Condens. Matter 2000, 12, R107. [Google Scholar] [CrossRef]
  28. Braunstein, S.L.; Caves, C.M.; Milburn, G.J. Generalized uncertainty relations: Theory, examples, and Lorenz Invariance. Ann. Phys. 1996, 247, 135. [Google Scholar] [CrossRef] [Green Version]
  29. Cafaro, C.; Mancini, S. An information geometric viewpoint of algorithms in quantum computing. AIP Conf. Proc. 2012, 1443, 374. [Google Scholar]
  30. Cafaro, C.; Mancini, S. On Grover’s search algorithm from a quantum information geometry viewpoint. Physica A 2012, 391, 1610. [Google Scholar] [CrossRef] [Green Version]
  31. Cafaro, C. Geometric algebra and information geometry for quantum computational software. Physica A 2017, 470, 154. [Google Scholar] [CrossRef] [Green Version]
  32. Amari, S.; Nagaoka, H. Methods of Information Geometry; Oxford University Press: Providence, RI, USA, 2000. [Google Scholar]
  33. Felice, D.; Cafaro, C.; Mancini, S. Information geometric methods for complexity. Chaos 2018, 28, 032101. [Google Scholar] [CrossRef]
  34. Castelvecchi, D. Clash of the physical laws. Nature 2017, 543, 597. [Google Scholar] [CrossRef] [Green Version]
  35. Faist, P.; Renner, R. Fundamental work cost of quantum processes. Phys. Rev. X 2018, 8, 021011. [Google Scholar] [CrossRef] [Green Version]
  36. Campbell, S.; Deffner, S. Trade-off between speed and cost in shortcuts to adiabaticity. Phys. Rev. Lett. 2017, 118, 100601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lee, J.M. Riemannian Manifolds: An Introduction to Curvature; Springer: New York, NY, USA, 1997. [Google Scholar]
Figure 1. Plot of the success probability p versus the temporal parameter θ in the case of constant driving for β = 0 (dotted), β 0 = 0.25 (dashed), β 0 = 0.5 (thin solid), and β 0 = 1 (thick solid).
Figure 1. Plot of the success probability p versus the temporal parameter θ in the case of constant driving for β = 0 (dotted), β 0 = 0.25 (dashed), β 0 = 0.5 (thin solid), and β 0 = 1 (thick solid).
Quantumrep 02 00011 g001
Figure 2. Numerical plots of the geodesic paths θ versus the affine parameter ξ for the four driving strategies: Power-law (dotted line), exponential (dashed line), oscillatory (thin solid line), and constant (thick solid line). In the LHS (left-hand-side), we consider the on-resonance scenario with β 0 = 0 . In the RHS (right-hand-side), we consider the off-resonance scenario with β 0 = 1 / 2 . In both sides, we set Γ / ħ = 1 , λ = 2 / π , θ 0 = 0 , and θ ˙ 0 = 1 .
Figure 2. Numerical plots of the geodesic paths θ versus the affine parameter ξ for the four driving strategies: Power-law (dotted line), exponential (dashed line), oscillatory (thin solid line), and constant (thick solid line). In the LHS (left-hand-side), we consider the on-resonance scenario with β 0 = 0 . In the RHS (right-hand-side), we consider the off-resonance scenario with β 0 = 1 / 2 . In both sides, we set Γ / ħ = 1 , λ = 2 / π , θ 0 = 0 , and θ ˙ 0 = 1 .
Quantumrep 02 00011 g002
Figure 3. Plot of the analytical (thin solid) and numerical (filled circle) values of the geodesic speed v OFF versus the initial condition θ 0 in the case of constant driving. In the plot, we set β 0 = 1 / 2 , Γ / ħ = 1 , and θ ˙ 0 = 1 .
Figure 3. Plot of the analytical (thin solid) and numerical (filled circle) values of the geodesic speed v OFF versus the initial condition θ 0 in the case of constant driving. In the plot, we set β 0 = 1 / 2 , Γ / ħ = 1 , and θ ˙ 0 = 1 .
Quantumrep 02 00011 g003
Figure 4. From left to right: (L) Plot of the geodesic speed v ON versus the initial condition θ 0 with β 0 = 0 ; (C) Plot of the geodesic speed v OFF versus the initial condition θ 0 with β 0 = 1 / 2 ; (R) Plot of the robustness coefficient r versus the affine parameter ξ . Plots that correspond to the constant, oscillatory, exponential decay, and power law decay strategies appear in dotted, dashed, thin solid, and thick solid lines, respectively. In all plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 .
Figure 4. From left to right: (L) Plot of the geodesic speed v ON versus the initial condition θ 0 with β 0 = 0 ; (C) Plot of the geodesic speed v OFF versus the initial condition θ 0 with β 0 = 1 / 2 ; (R) Plot of the robustness coefficient r versus the affine parameter ξ . Plots that correspond to the constant, oscillatory, exponential decay, and power law decay strategies appear in dotted, dashed, thin solid, and thick solid lines, respectively. In all plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 .
Quantumrep 02 00011 g004
Figure 5. Plot of the two-dimensional parametric regions P β 0 , θ 0 where the oscillatory (light grey), the exponential decay (grey), and the power law decay (black) driving strategies outperform the constant driving scheme in terms of both geodesic speed and robustness. Note that the light grey region contains the grey region while the grey region contains the black region. Furthermore, the black region specifies a driving strategy more robust than that specified by the grey region while the grey region specifies a driving strategy more robust than that specified by the light grey region. In all region plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 . Finally, the maximum success probability of each driving scheme is assumed to be greater than 25 / 26 0.96 , that is, 0 β 0 0.20 .
Figure 5. Plot of the two-dimensional parametric regions P β 0 , θ 0 where the oscillatory (light grey), the exponential decay (grey), and the power law decay (black) driving strategies outperform the constant driving scheme in terms of both geodesic speed and robustness. Note that the light grey region contains the grey region while the grey region contains the black region. Furthermore, the black region specifies a driving strategy more robust than that specified by the grey region while the grey region specifies a driving strategy more robust than that specified by the light grey region. In all region plots, we set Γ / ħ = 1 , λ = 2 / π , and θ ˙ 0 = 1 . Finally, the maximum success probability of each driving scheme is assumed to be greater than 25 / 26 0.96 , that is, 0 β 0 0.20 .
Quantumrep 02 00011 g005

Share and Cite

MDPI and ACS Style

Cafaro, C.; Gassner, S.; Alsing, P.M. Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems. Quantum Rep. 2020, 2, 166-188. https://doi.org/10.3390/quantum2010011

AMA Style

Cafaro C, Gassner S, Alsing PM. Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems. Quantum Reports. 2020; 2(1):166-188. https://doi.org/10.3390/quantum2010011

Chicago/Turabian Style

Cafaro, Carlo, Steven Gassner, and Paul M. Alsing. 2020. "Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems" Quantum Reports 2, no. 1: 166-188. https://doi.org/10.3390/quantum2010011

APA Style

Cafaro, C., Gassner, S., & Alsing, P. M. (2020). Information Geometric Perspective on Off-Resonance Effects in Driven Two-Level Quantum Systems. Quantum Reports, 2(1), 166-188. https://doi.org/10.3390/quantum2010011

Article Metrics

Back to TopTop