{'fraction_non_alphanumeric': 0.03735640255034688, 'fraction_numerical': 0.012639384321545937, 'mean_word_length': 4.600755073147711, 'pattern_counts': {'":': 0, '<': 15, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Double Occupancy Errors in Quantum Computing Operations: Corrections to Adiabaticity

We study the corrections to adiabatic dynamics of two coupled quantum dot spin-qubits, each dot singly occupied with an electron, in the context of a quantum computing operation. Tunneling causes double occupancy at the conclusion of an operation and constitutes a processing error. We model the gate operation with an effective two-level system, where non-adiabatic transitions correspond to double occupancy. The model is integrable and possesses three independent parameters. We confirm the accuracy of Dykhne's formula, a nonperturbative estimate of transitions, and discuss physically intuitive conditions for its validity. Our semiclassical results are in excellent agreement with numerical simulations of the exact time evolution. A similar approach applies to two-level systems in different contexts.


Introduction
Quantum information processing is an active and fascinating direction of research with participation from various fields of physics and neighboring scientific disciplines [1]. This extraordinary interest has generated a fairly vast amount of theoretical and experimental studies. Possible experimental realizations of quantum information processing are presently being investigated. Among the different approaches, those in a solid state setting are attractive, because they offer the potential of scalability -the integration of a large number of quantum gates into a quantum computer once the individual gates and qubits are established. With that in mind, several proposals for using electron and/or nuclear spins in solid state systems have been put forward in recent years [2,3,4,5,6,7]. Specifically, in Ref. [2] it was proposed to use the spin of electrons residing in semiconductor quantum dots as qubits [8,9,10,11,12,13,14,15]. In this paper we revisit the quantum dynamics of gate operations between qubits of this type. Such two-qubit operations are performed by varying the amplitude of electron tunneling between the dots via external electric potentials. In a generic scenario, the tunneling amplitude between the dots is zero (or, more precisely, exponentially small) before and after the gate operation, while it is finite and appreciable during such a process. Thus, the typical time dependence of the tunneling amplitude is a pulse roughly characterized by its duration, amplitude, and ramp time (see Figure 1.) During such a pulse, the tunneling amplitude is finite and essentially constant, and both electrons can explore the total two quantum dot system. Therefore, their indistinguishable fermionic character is of relevance [10,16,17]. In particular, in such gate operations entanglement-like quantum correlations arise which require a description different from the usual entanglement between distinguishable parties (Alice, Bob, ...) in bipartite (or multipartite) systems. In such a case the proper statistics of the indistinguishable particles has to be taken into account [10,16,17].
Another important aspect of having a finite (as opposed to infinitely high) tunneling barrier between the dots is that it necessarily leads to (partially) doubly occupied states in the two-electron wave function, i.e. contributions to the wave function where both electrons are on the same dot (having different spin) occur with finite amplitude. Doubly occupied states which arise as the result of a measurement after the gate operation, destroy the information in those qubits and lead to errors in the information processing. Therefore, it is desirable to reduce the probability of such errors, i.e. the occurrence of doubly occupied states, in the resulting two-electron state after the gate operation, while it is necessarily finite during the operation [9,10]. If the error probability can be sufficiently reduced, error events can be tolerable and handled with quantum error correction schemes. An effective way of guaranteeing error suppression is to maintain nearly adiabatic time evolution. Doubly occupied states then correspond to corrections to adiabatic evolution, which are often called "non-adiabatic transitions". Numerical simulations [10] have shown that the adiabatic region, in terms of the pulse parameters such as ramp time and amplitude, is rather large. On a heuristic level, this numerical result is plausible on the basis of the classic papers on adiabatic quantum motion in two-level systems by Landau [18], Zener [19], Stueckelberg [20], and Rosen and Zener [21]. For an overview see Ref. [22].
In this work we study the quantum dynamics of the two-qubit gate operations described above and use Dykhne's semiclassical result to estimate the probability of non-adiabatic transition [23]. The applicability of Dykhne's formula is analyzed from the standpoint of the theory of semiclassical approximations. These semiclassical estimates are found to be in excellent agreement with numerical simulations of the exact time evolution. Moreover, in a certain limit our model is integrable, allowing us to explicitly calculate and interpret the corrections to Dykhne's formula.
This paper is organized as follows. Section 2 reviews the derivation [10] of an effective two-level model. In section 3, we present our main result -the asymptotic estimate of double occupancy, which in section 4 is compared with an integrable model and a numerical integration of the Schrödinger equation. In the Appendix, we construct the scattering matrix for the integrable model, which has three independent parameters. 3

Mapping to an Effective Two-level System
For the purpose of studying double occupancy it is practical to examine the dynamics of the quantum gate operation in a subspace spanned by singly and doubly occupied states. Following Ref. [10] with only minor changes of notation, we now detail how to reduce the description of a system of two coupled quantum dot spin qubits to an effective two-level Hamiltonian. The system is described by a Hamiltonian of the form H = T + C, where C denotes the Coulomb repulsion between the electrons, and T = i=1,2 h i is the single-particle part with The single-particle Hamiltonian h i describes electron dynamics confined to the xyplane in a perpendicular magnetic field B. The effective mass m is a material dependent parameter. The coupling of the dots (which includes tunneling) is modeled by a quartic potential which separates into two harmonic wells of frequency ω 0 (one for each dot) in the limit a ≫ a 0 , where a is half the distance between the dots and a 0 = h/mω 0 is the effective Bohr radius of a dot. Following Burkard et al. [8] we employ the Hund-Mulliken method of molecular orbits to describe the low-lying spectrum of our system. This approach concentrates on the lowest orbital states in each dot and is an extension of the Heitler-London method [8]. The Hund-Mulliken approach accounts for the fact that both electrons can, in the presence of a finite tunneling amplitude, explore the entire system of the two dots, and therefore adequately includes the possibility of doubly occupied states. In the usual symmetric gauge A = B(−y, x, 0)/2 the Fock-Darwin ground state of a single dot with harmonic confinement centered around r = (±a, 0, 0) reads where l B = hc/eB is the magnetic length, and the frequency ω is given by ω 2 = ω 2 0 + (ω L /2) 2 where ω L = eB/mc is the usual Larmor frequency. From these nonorthogonal single-particle states we construct the orthonormalized states |A and |B 4 with wave functions with S being the overlap between the states (3) and g = (1 − √ 1 − S 2 )/S. For appropriate values of system parameters such as the interdot distance and the external magnetic field, the overlap S becomes exponentially small [8]. In this limit an electron in one of the states |A , |B is predominantly localized around r = (±a, 0, 0). In the following we consider this case and use these states as basis states to define qubits, i.e. qubits are realized by the spin state of an electron in either orbital |A , or orbital |B .
An appropriate basis set for the six-dimensional two-particle Hilbert space is given (using standard notation) by the three spin singlets and the triplet multiplet, As the Hamiltonian conserves spin, the three triplet states are degenerate eigenstates (typically we can ignore possible Zeeman splittings [8]) and have the eigenvalue, where we have defined and the expectation value of Coulomb energy An important further observation is that, as a consequence of inversion symmetry along the axis connecting the dots, the Hamiltonian does not have any non-zero matrix elements between the singlet state |S 3 and other states. Hence, |S 3 is, independently of the system parameters, an eigenstate. The eigenvalues of the triplet and |S 3 , however, do depend on system parameters. The Hamiltonian acting on the remaining space spanned by |S 1 and |S 2 can be written as where and The nontrivial part of (15) is a simple Hubbard Hamiltonian on two sites and can be identified as the Hamiltonian of a pseudospin-1 2 object in a pseudomagnetic field having a component U H in theẑ-direction and U H t H in thex-direction of pseudospin space. [Note that this pseudospin is not related to the spin degree of freedom which constitutes the qubit!] The basis states themselves are eigenstates only in the case of vanishing tunneling amplitude t H where |S 1 is the ground state and |S 2 is a higher lying state due of the Coulomb (Hubbard) energy. In all other cases, the ground state has an admixture of doubly occupied states contained in |S 2 . The energy gap between the triplet and the singlet ground state is A key challenge for state-of-the-art quantum information processing is the construction of systems composed of two coupled quantum dots which can be coupled to perform swap operations U SW , i.e. unitary two-qubit operations which interchange the spin states (qubits) of the electrons on the two dots. By combining the "square root" U 1/2 SW of such a swap with other isolated-qubit manipulations one can construct a quantum XOR gate. A quantum XOR gate, along with isolated-qubit operations, has been shown to be sufficient for the implementation of any quantum algorithm [25]. Hence a practical and reliable realization of a swap gate would be an important step toward the fabrication of a solid state quantum computer. A swap operation in the present system is a unitary transformation which turns a state having the qubits in different states, say, into a state where the contents of the qubits are interchanged, These two states are eigenstates in the case V + = V − and t H = 0 for which the singlet-triplet splitting vanishes.
As discussed in references [2,8,10], swapping may be achieved by the action of a gate that lowers the potential barrier between the quantum dots. This leads to exponentially larger values for both V + −V − and t H . It is adequate for our purposes to consider a model where V + = V − (consistent with the above limit of small overlap S), and the singlet-triplet splitting results entirely from t H . If the duration and amplitude of a tunneling pulse is adjusted appropriately, the relative dynamical phase between the singlet and the triplet state accumulates a shift of π, and the swapping operation between states (19) and (20) is performed. However, during the operation the state |S 1 is coupled to |S 2 , and they evolve according to (15). Double occupancy errors are thus generically introduced. The reduction of the dynamics to the time evolution of a two-level system relies on the fact that the system has inversion symmetry along thex-axis in real space connecting the dots. This symmetry can be broken if odd powers of the particle coordinates x i are added to the Hamiltonian (1), for example the potential of a homogeneous electric field. The breaking of inversion symmetry introduces additional matrix elements between |S 3 and the other two singlets leading to an effective three level Hamiltonian. However, as it was shown in Ref. [10], this more inclusive Hamiltonian has qualitatively the same properties concerning non-adiabatic dynamics as the two-level system on which we shall concentrate in the following.
So far we have not considered a possible Zeeman coupling to the electron spin. This would not change the situation essentially since all states involved in the swapping process (|T 0 , |S 1 , |S 2 , and possibly |S 3 ) have the total spin quantum number S z = 0.

Analysis of Non-adiabatic Transitions
In this section we use Dykhne's formula for non-adiabatic transitions to derive an asymptotic expression for the probability of final double occupancy, given physically motivated properties of the two-qubit operation.
As described in the previous section, the modulation of the tunneling barrier during the swapping process induces a coupling between the singly occupied qubit state |S 1 and the doubly occupied state |S 2 . Their dynamics are governed by the effective Hamiltonian in the |S 1,2 basis. The terms omitted from (15) do not contribute to transitions, because the identity operator in the |S 1,2 basis commutes with the remainder of the Hamiltonian. The large energy offset U H between singly and doubly occupied states, primarily due to the Coulomb repulsion, is perturbed only by an exponentially small additive quantity (proportional to the overlap, S) during the swapping operation and is hereafter assumed to be a constant. Our specification of the pulse ( Figure 1) with dimensionless strength δ is considered to realistically reflect the tunneling amplitude that would arise from a modulation of the gate potential [10]. The exponential dependence of the ramping near t = ±T /2 has its origin in the exponential sensitivity of the coupling to the gate voltage and in turn the exponential decay of the single-particle wavefunctions (3) in the inter-dot region [26]. The pulse mimics a step of duration T and magnitude δU H /2, whose ramping on and off have a characteristic time τ . The perturbation of the instantaneous eigenvalues by the pulse is shown in Figure 2. The Schrödinger equation is Our task is to find the component of double occupancy in the final state, S 2 |ψ(∞) , given that the prepared state is purely singly occupied, | S 1 |ψ(−∞) | = 1. Our model involves three dimensionless scales, assigned for our purposes as follows: δ, λ ≡ U H τ /2h, and η = T /τ . Presently, the case of interest is The first of these conditions reflects the adiabaticity of the problem. The second requires that the ramping on and ramping off of the pulse be temporally well-separated and distinct events.    (22) is diagonal. By force of the adiabatic theorem, the probability of transition among eigenstates vanishes in the limit τ → ∞, where the ramping on and off of the pulse is adiabatic. In the zeroth order of adiabatic perturbation theory, there are no transitions, and the leading behavior of the general solution is simply the dynamical phase of each component eigenstate where |ξ 1,2 (t) are the instantaneous eigenstates (given explicitly in (54)) of Hamiltonian (22) corresponding to eigenvalues respectively. In general, the zeroth order approximation should also include a factor representing Berry's phase. However, for a real symmetric Hamiltonian such as (22), Berry's phase is irrelevant, because the Hamiltonian has an inherent planarity. In pseudospin one-half notation, H eff = H · σ, the time evolution of the pseudomagnetic field H = H(t) is in a plane. If the azimuthal axis (north pole) is chosen to lie within that plane, the solid angle subtended by the pseudomagnetic field vanishes identically. Although Berry's phase is out of consideration, there are interesting circumstances where Berry's phase is relevant to transitions. It can correct the transition amplitude [31] and produce topological selection rules for spin tunneling [32,33]. Our problem is one of a class initiated by the work of Landau, Zener, and Stueckelberg (LSZ) [18,19,20]. However, we emphasize that for our model (with the pulse specified as (23)) the linearization of Hamiltonian matrix elements near the times where adiabaticity is most severely violated is not applicable and leads to an incorrect result. As we will see the shape of the pulse is important.

Application of Dykhne's Formula
Returning to our model, we observe that if the time interval t ∈ (−∞, ∞) is divided into two domains t < 0 and t > 0, and in the limit η ≡ T τ ≫ 1, the pulse (23) is approximated by In each domain the pulse behaves as a step, and the dynamics are integrable (see Section 4.) We will focus first on the interval t < 0, where the probability of transition to a doubly occupied state P < may be estimated with Dykhne's formula [23] where the approximation (28) is used implicitly for the instantaneous eigenenergies ∓ε(t) defined above by (27). The turning point t = t 1 , given explicitly below, is a complex root of the function ε = ε(t); in other words, it is an intersection of the energy surfaces (curves) of the two instantaneous ("frozen") eigenstates. Our model is the patching together of two domains of time, and transitions that occur during t < 0 and t > 0 interfere. The expression for the probability of transition during the time evolution from t = −∞ to t = ∞ is where the contours C a,b are shown in Figure 6, and according to the sign of the integration variable, sgn(ℜ z), one or the other of the approximations (28) is used. The turning points t = t 1,2 appearing in the limits of integration of (30), are chosen as the two roots of ε = ε(t) that are closest to and above the real time axis (see Figure 6) They are nonreal because the Hamiltonian (22) is nondegenerate for real times. Equation (31) follows from (30), because the symmetry of the pulse implies ℑ(t 1 ) = ℑ(t 2 ) and P < = P > . The oscillatory first factor of (31) is the interference of the dynamical phase of each term of (30). The magnitude of P is dominated by the second factor, P < , whose exponent is given by the following integral Substituting this result in (29) we have From (31), we have our main result, an asymptotic estimate for the probability of final double-occupancy which is shown as a function of δ in Figure 3. The probability P is characteristically nonperturbative in the adiabatic limit τ → ∞ with U H fixed, or equivalently λ → ∞. Hence, the dimensionless quantity associated with the exponential suppression is λ and has been called the "adiabaticity parameter." For η ≡ T τ ≫ 1, the approximation (28) allows us to estimate the argument of the prefactor of (35) to exponential accuracy, The oscillation with respect to the duration of the pulse T is reminiscent of a similar factor in the Rosen-Zener model. The phenomenon of pulsed perturbations that return the full amplitude/occupation to the initial state have been studied in the context of atom-laser interactions [27,28,29,30]. In Figures 3, 4, and 5, we compare our semiclassical estimate (35) with results from numerical simulations of the exact quantum mechanical time evolution, following Ref. [10]. Both results are in excellent agreement and differ only at very small δ, i.e. for weak pulses. Of course, the non-adiabatic transition probability vanishes in this limit, whereas the semiclassical approximation breaks down (see Section 3.3). This regime is beyond the exponential accuracy of Dykhne's formula. The integrability of our model allows us to make precise statements about the form and magnitude of the corrections to Dykhne's formula (see Section 4). For example, in the limit λ ≫ 1 and δλ ≪ 1 we have from the expansion (63) that P < ∼ (2πδλ) 2 e −4πλ , while in the same limit the result of Dykhne's formula (35) gives only the exponential factor e −4πλ without information about the prefactor. This explains a trend among Figures 3, 4, and 5, namely the increasing range, in terms of δ, of validity of Dykhne's formula with increasing λ. The value for the adiabaticity parameter λ = 2, represented in Figure 3, corresponds to a ramp time τ = 4h/U H , which was identified in Ref. [10] as a practical lower bound to ensure sufficient adiabatic behavior in a gate operation between two quantumdot spin-qubits. It is interesting that Dykhne's formula remains accurate for smaller  values of λ, in particular λ = 1 as seen in Figure 5. The reason is that the results (35) and (61) have an incidental factor of 2π in the exponent, giving in practical terms the requirement for exponential suppression 2πλ ≫ 1. The expressions (34) and (35), along with Figures 3, 4, 5 comprise our main results. For the remainder of this section we will address the justification and limitations of these results.

Origin of Dykhne's Formula
Dykhne derived a concise expression for non-adiabatic transitions from a local analysis of the Schrödinger equation in the vicinity of the turning point [23]. Dykhne's formula can be viewed as a semiclassical approximation, and an elegant interpretation and proof was given by Hwang and Pechukas [24] (see also [34].) We will briefly discuss  the key elements and scope of the proof. Their method was to study the solution of the Schrödinger equation in the complex plane of the independent variable, time.
According to the adiabatic theorem, the projection of the solution onto any eigenstate other than the initial eigenstate approaches zero in the adiabatic limit. One might suppose that weak statement is all the adiabatic theorem can tell us about transitions; however, it does not exhaust its capacities. The reason lies in the following: a basis of eigenstates |ξ 1,2 (t) , when extended into the complex time plane, is multi-valued.
In particular, as a basis state is analytically continued across a branch cut of the function ε = ε(t), its long time asymptotics are discontinuously changed. In accord with our above two-level problem, we uniquely specify the basis by its asymptotics The multi-valued nature is not manifest on the real time axis, because owing to the nondegeneracy of the spectrum ±ε(t), the branch points are nonreal. We can choose a single-valued basis ξ 1,2 (t) , which makes reference to |ξ 1,2 (t) but has fixed  asymptotics, by defining rules for continuing the basis states across branch cuts. Equivalently, this new basis is said to be defined over a Riemann surface with two sheets (copies of the complex time plane), one corresponding to each of the two branches of the function [ε(t) 2 ] 1/2 . Crossing a branch cut means passing to the other sheet of the Riemann surface. We assign the following relations among the eigenstates ξ 1,2 (t) = |ξ 1,2 (t) t ∈ sheet 1 and ξ 1,2 (t) = e i α 1,2 |ξ 2,1 (t) t ∈ sheet 2 (38) where α 1,2 are phase conventions that are chosen to maintain continuity of the basis ξ 1,2 across the branch cut (no physical quantity will depend on α 1,2 .) Given |ψ(−∞) = ξ 1 (−∞) , the conclusion of the adiabatic theorem may be restated on a Riemann surface as Figure 6: The analytic structure of the function ε = ε(t) shown only in a segment of the upper half plane. The contour C a is associated with transitions that occur due to the ramping on of the pulse, while contour C b is associated with the ramping off. By Cauchy's theorem, an integral on the contour C a is equal to the integral on the contour C ′ a . Bold lines represent branch cuts, dots represent branch points, and poles are denoted with a ×.
where τ is the characteristic time scale for variation of H eff (t). The only exception to (39) is for times within τ O(λ −2/3 ) of a turning point, for there the semiclassical criterion (45) is invalid. As remarked above, the zeroth order approximation (26) of the solution as τ → ∞ is the dynamical phase. The zeroth order approximation may be extended into the complex plane by evaluating the dynamical phase on a contour C. Continuing with the above example, a state that is purely singly occupied at t = −∞ is for complex time given by where C is a contour from z = −∞ to z = t. The amplitude of transition is readily obtained as the projection of the solution, as t → ∞ on the second Riemann sheet (see Figure 7), onto the doubly occupied state |S 2 , i.e.
where the contour C crosses the branch cut emanating from the branch point that is closest to the real axis. Dykhne's formula is simply the square modulus of this amplitude.
In the adiabatic regime, in contrast to the perturbative regime, the leading contribution to transitions comes from the zeroth order term of perturbation theory instead of the first order term. By retaining only the zeroth order term, it appears that we have neglected completely the coupling t H between states. However, the coupling enters implicitly in the multivalued function [ε(t) 2 ] 1/2 and influences the location of the turning points -the complex roots of ε(t). Transition amplitudes are obtained by carefully considering the different branches of this function. In the next section, we consider the validity of keeping only the zeroth order term.

Validity and Accuracy of Dykhne's Formula
The theory of semiclassical approximations, especially WKB analysis, provides a foundation from which to evaluate the validity of Dykhne's formula. The calculation of non-adiabatic transitions is closely related to the semiclassical approximation [24], because the semiclassical limith → 0 can be mathematically equivalent to the adiabatic limit τ → ∞. An essential element of the proof by Hwang and Pechukas is the existence of a complex time contour that 1) connects the two sheets of the Riemann surface and 2) on which the zeroth order approximation of adiabatic perturbation theory is the correct leading behavior of the solution in the adiabatic limit. These are sufficient conditions for Dykhne's formula to give the correct asymptotic form of the transition probability in the adiabatic limit λ → ∞. Having established the existence of such a contour, one can calculate a more precise value for the prefactor of Dykhne's formula by applying time dependent perturbation theory along the contour. We expect Dykhne's formula to breakdown when the contour ceases to exist. At the limit of its range of validity, the higher order terms become comparable to the zeroth order term. Introducing the unitary transformation U that diagonalizes the Hamiltonian, i.e. U † HU = εσ 3 , we can write the Schrödinger equation in the basis of instantaneous eigenstates with the off-diagonal perturbationâ(t) = −U † i ∂ t U. A dominancy balance among the terms gives the condition for the accuracy of the zeroth order approximation For our model of the dynamics, λ ∝ τ is the largest scale and |â| ∼ 1. The condition (44) must be maintained at all points on the contour. Applying (44) on the real axis, where |ε(x)| ∼ U H , gives the adiabaticity condition λ ≫ 1. Additionally, in order to connect two Riemann sheets, the contour must pass between two turning points (see Figure 6), where |ε(x)| ∼ δU H , giving the condition δλ ≫ 1.
Beginning instead from an intuitive approach, we can evaluate the adiabaticity of the dynamics along a given contour. To test whether a given contour is adequate, we can exploit the analogy between quasi-adiabatic dynamics and semiclassical scattering. Recall the semiclassical criterion The analog of the de Broglie wavelength λ(x) = 2πh/p(x) in scattering problems is the period Λ(t) ≡ 2πh/ε(t). In other words, the condition (45) says that the change of period over the course of one period is small. We now require that the semiclassical criterion be obeyed everywhere along an admissible contour, i.e. one that connects the two Riemann sheets. To find an admissible contour, we must appeal to the analytic structure of the eigenenergy ε = ε(t), see Figure 6. For clarity we will focus on the time interval t < 0 and operate under the approximation (28). The singularities of ε(t) are branch points at t = − T 2 − τ ln(−1 ± iδ) and poles at t = − T 2 − τ ln(−1).
If we agree to define a branch cut connecting the nearest and next nearest branch points to the real time axis ) and respectively, then an admissible contour is one that crosses this branch cut exactly once. For the semiclassical criterion to be obeyed, the admissible contour cannot pass too close to a branch point. In essence, if δ is too small, the contour is pinched between the branch points t 1 and t 3 . Evaluating the maximum of dΛ/dt on a contour C that crosses the branch cut between t = t 1 and t = t 3 we arrive at the condition δ ≫ λ −1 . Together with the adiabatic limit λ ≫ 1, we have the following conditions on the interdependence of the physical parameters Each of these dimensionless quantities is a product of a characteristic energy and time scale. If these conditions are not satisfied, there does not exist a contour on which the motion is adiabatic. The integrability (Section 4) of our model allows us to investigate the intermediate regime λ ≫ 1 and δλ ≪ 1, where Dykhne's formula cannot be justified with the analysis of Hwang and Pechukas. with The Hamiltonians H < and H > can be obtained as a special case of by identifying ±c = a = −δλ/2, b = −λ and rescaling time x = t/τ ± T /2τ , respectively. The Schrödinger equation is exactly solvable [36] (see Appendix.) In analogy with one-dimensional scattering, the transition amplitude from a singly occupied state |S 1 to a doubly occupied state |S 2 may be viewed as an off-diagonal element of the scattering matrix S that connects the coefficients of the asymptotic final states to the asymptotic initial states. The asymptotic states are the limit as t → ±∞ of the instantaneous eigenstates |ξ 1,2 (t) of H eff corresponding to eigenvalues ∓ε(t), respectively, The leading behavior of the long time asymptotics of a general solution has the form and the scattering matrix relates the coefficients 20 Referring to the statement of our problem in Section 3, the amplitude of final double occupancy is the element S 21 of the scattering matrix, which we parametrize as The scattering matrices S < and S > associated with the Hamiltonians H < and H > may be obtained by substitution from the exact scattering matrix W (derived in the Appendix) associated with the Hamiltonian H exact . By the symmetry of the pulse, we have S > = S † < . Patching together the two domains of time evolution (49), we find the scattering matrix where the integral of the exponent has been estimated in (36), and the elements of S < are obtained from W where µ = λ √ 1 + δ 2 . Dykhne's formula (34) for P < is recovered as exactly the leading term of |(S < ) 21 | 2 in the limit λ, δλ ≫ 1, The nonperturbative corrections are typically very small. For our problem λ = 2, δ = 1/2, and the relative contribution of the second term of (62) is less than one percent. This accounts for the excellent agreement in Figures 3, 4, and 5, between the probability of double occupancy as given by the semiclassical result (35) based on Dykhne's formula and the result of a numerical integration of the Schrödinger equation. We can interpret the subleading term of (62) as the contribution from the t Re(t) Figure 8: The contour C gives a subdominant correction to Dykhne's formula.
contour C of Figure 8, which crosses the branch cut three times. The sign of the correction is negative and arises from the factor e iα 1 e iα 2 associated with matching the basis (38) across the branch cut. For our problem α 1 = 0 and α 2 = π. Similarly, it may be possible to extend Dykhne's formula and reproduce the exact result (61) by summing over all homology classes of the Riemann surface and within each homology class over all paths that give distinct nonzero values for ℑ C dt ε(t) [35]. For many physical problems this type of nonperturbative correction is dominated by perturbative corrections along the contour -those mentioned in Section 3.3. The striking absence of perturbative corrections in the limit λ, δλ ≫ 1 is a unique artifact of the integrability of our model. With a knowledge of the exact result, we can also investigate the intermediate regime λ ≫ 1 and δλ ≪ 1, where the analysis of Hwang and Pechukas (Section 3.2) cannot be used to prove Dykhne's formula. In this limit, the transition probability (61) becomes Although Dykhne's formula does not apply in this limit, because δ −1 and not λ is the largest scale, it nevertheless gives the correct controlling factor e −4πλ of (63), except for δ that are exponentially small with respect to λ. This exponential factor is resilient and remains the controlling factor for a range of parameters beyond the naive expectation based on the arguments of Section 3.3.

Conclusions
The dynamics of two coupled quantum dot spin-qubits can be mapped to an effective two-level system, where non-adiabatic transitions correspond to double occupancy. We have estimated the probability of final double occupancy with Dykhne's formula.
In the adiabatic regime, the pervasive feature of transitions is their exponential suppression by a dimensionless adiabaticity parameter, λ. Our main result (35) was expressed in terms of the dimensionless quantities λ, δ and η. An integral constraint (21) on the swapping operation gives one relation among the three dimensionless parameters. The problem is uniquely defined by specifying any two, and in a solid state setting, conservative estimates are λ ≈ 2 and δ ≈ 1/2. The probability of double occupancy P ≈ 10 −10 is sufficiently rare that the operation of a quantum gate will not be obstructed by this type of error. It is noteworthy that the probability of double occupancy (35) has nodes for However, this property is not immediately relevant to the suppression of transitions, because the oscillatory factor, sin 2 √ 1 + δ 2 λη , of (35) vanishes algebraically and for it to provide an improvement upon the exponentially small factor from Dykhne's formula, the argument would have to be tuned exponentially close to kπ. Thus, naively the errors associated with inaccuracies in satisfying the integral swapping constraint (21) will be much greater than double occupancy errors. Other important sources of error are dephasing and decoherence of the qubit states.
We have reviewed a physically motivated derivation of Dykhne's formula [24]. The theory of semiclassical approximations underlies Dykhne's formula and its validity is appropriately judged within that framework. The semiclassical estimates obtained from this approach are in excellent agreement with numerical simulations of the full quantum mechanical time evolution. The corrections to Dykhne's formula are of two types: perturbative and nonperturbative. The former appear to vanish for integrable models, and we have interpreted a nonperturbative correction as the contribution of a contour in the complex time plane [35].

Acknowledgments
AGA and RR wish to thank Harold Metcalf for helpful discussions. JS and DL acknowledge financial support from the Swiss NSF, the NCCR Nanoscience, EU RTN Spintronics, DARPA, ARO, and ONR. AGA received support from the Theory Institute for Strongly Correlated and Complex Systems at BNL as well as from NSF grant DMR-0348358.

A Integrable Dynamics
The following Hamiltonian [36], which possesses three independent parameters has integrable dynamics in the sense that the Schrödinger equation, has a solution in terms of special functions. Our aim is to construct the scattering matrix W associated with the dynamical problem. For clarity of presentation, we will consider the Hamiltonian which differs from H exact by a constant unitary transformation V Introducing the mixing angle φ defined as tan φ = b a + c tanh