Avoiding leakage and errors caused by unwanted transitions in Lambda systems

Three-level Lambda systems appear in various quantum information processing platforms. In several control schemes, the excited level serves as an auxiliary state for implementing gate operations between the lower qubit states. However, extra excited levels give rise to unwanted transitions that cause leakage and other errors, degrading the gate fidelity. We focus on a coherent-population-trapping scheme for gates and design protocols that reduce the effects of the unwanted off-resonant couplings and improve the gate performance up to several orders of magnitude. For a particular setup of unwanted couplings, we find an exact solution, which leads to error-free gate operations via only a static detuning modification. In the general case, we improve gate operations by adding corrective modulations to the pulses, thereby generalizing the DRAG protocol to Lambda systems. Our techniques enable fast and high-fidelity gates and apply to a wide range of optically driven platforms, such as quantum dots, color centers, and trapped ions.


I. INTRODUCTION
Quantum information processing requires the manipulation of quantum bits (qubits) via fast gates with high fidelities. Qubits are formed when a particular two-level subspace is chosen from a larger Hilbert space of a physical system. In specific cases, energy levels outside of the qubit subspace are used as an asset for auxiliary transfer of population within the qubit subspace and for the implementation of quantum gates. An important class of such setups are Λ-type systems that occur in several optically active qubit systems such as self-assembled quantum dots [1][2][3][4][5][6], nitrogen-vacancy (NV) centers [7][8][9][10], trapped ions [11][12][13][14][15], rare-earth ions [16,17], molecular qubits [18] and even in the microwave regime of superconducting circuits [19,20]. Successful manipulation of Λ-systems is a key step to perform quantum information processing.
Various methods for the control of Λ-systems have been developed [2,3,[21][22][23][24][25][26][27][28]. However, in most platforms a bare three-level Λ-system is merely an idealization [29]. Unintended interactions with other levels cause leakage and off-resonant couplings that are detrimental to the performance of quantum gates [30,31]. The effect of these off-resonant unwanted couplings is intensified when the extra levels are close in energy to the auxiliary state. In principle, we can use longer pulses to resolve such small splittings. However, the qubit coherence time sets an upper bound to the duration of the pulse. Finding fast pulses that satisfy these opposing constraints and achieve high-fidelity gates is an open problem.
In this paper, we develop a novel framework relevant for systems with Λ-type selection rules, where the auxiliary excited state, hereafter referred to as the target state, is separated from an unwanted nearby level by a small energy splitting (Fig. 1). The key insight of the presented formalism is that the unwanted level induces two different types of errors (leakage and phase errors), and that each is more prominent in a different regime.

"
|ui |ti |0i |1i E 0 (t) E 1 (t) Figure 1. (Color online) Schematic depiction of a general Λ-type system with a fourth unwanted level. The qubit is defined in the subspace of the two lower levels {|0 , |1 }. The coherent control is done through driving the transitions to the target level |t (with the detuning δ) which is separated from an unwanted level |u by the splitting ε. The off-resonant coupling to the unwanted level will cause low fidelities. We present resolutions for this problem by introducing a modification to the pulse and the detuning of the system.
First, we consider the case that the target and unwanted states are formed from a set of two basis states. We show that in this case the unwanted couplings lead to phase errors, and the problem can be solved exactly. In this case, quantum gates with minimal gate error can be implemented through a modification of the detuning. Second, we consider the case that no basis state structure is present. We show that in this case the errors are in the form of leakage and we develop a novel version of the Derivative Removal by Adiabatic Gates (DRAG) framework for this system that battles the leakage error. Former versions of the DRAG method have been widely established as a powerful tool in dealing with unwanted dynamics and leakage cancellations in superconducting qubits [32][33][34][35][36]. However, no version of this method applicable to Λ-type systems has yet been developed [16].
The paper is structured as follows: In Section II, we present an overview of the Λ-system with unwanted arXiv:2205.06945v1 [quant-ph] 14 May 2022 (Color online) Fidelity of an X-rotation by −π/2, RX (−π/2), in terms of the dimensionless parameter pulse bandwidth over splitting (σ/ε). Here, σ is the bandwidth of the sech pulse and ε is the splitting between the unwanted and target levels. The fidelities are shown for two cases of dependent (λ0 = − tan(π/4) = −1/λ1) and independent (λ0 = 1, λ1 = 1) couplings. Without any corrective measures and ignoring spontaneous emission, reasonable fidelities required for quantum information processing are only achievable by using extremely narrow bandwidth pulses (solid lines). However, upon inclusion of spontaneous emission (dashed lines), even narrow bandwidths pulses will not lead to perfect fidelities.
transitions and discuss the relevance of the composition of the target and unwanted levels in terms of a set of two basis states. In Sec. III we start with a hyperbolic secant uncorrected pulse and develop an exact analytical solution in the presence of such a basis state structure. In Sec. IV we present the case without such a basis state structure, discuss the DRAG methodology, and develop a new DRAG-based formalism that resolves the off-resonant coupling issue in this case. In Sec. V We analyze the effects of additional errors, namely the cross-talk between the transitions of the Λ-system and spontaneous emission. We conclude in Sec. VI. The appendices contain the technical details of the CPT scheme and our DRAG formalism.

II. Λ-SYSTEM WITH UNWANTED TRANSITIONS UNDER CPT
In an ideal three-level Λ-system under CPT, the transitions are driven using two external fields (e.g. E 0 (t) and E 1 (t) in Fig. 1). Under the two-photon resonance (equal detuning of the two fields from the respective transitions they drive), destructive interference of the transitions caused by the two drive fields leads to the formation of a dark state, i.e. a state that is completely decoupled from the dynamics of the other two levels. This allows us to define the CPT frame, in which the system is described in terms of two new states in the qubit subspace, the dark state and its orthogonal bright state (denoted by |D and |B , respectively), which are superpositions of the |0 and |1 states. Then, the three-level system reduces to a two-level system where transitions are driven between the target excited state and the bright state. When we implement CPT with pulses that have hyperbolic secant (sech) envelopes, a choice that yields an analytically solvable time-dependent Schrödinger equation in a twolevel system [37], we can design arbitrary single-qubit rotations for the Λ-system [3]. The parameters of the sech pulses that drive the bright and target state can be chosen such that the evolution is transitionless [2]: after the passage of a sech pulse, the population will always return to the bright state, with the bright state acquiring a non-trivial phase where σ is the bandwidth, and δ is the detuning. Considering that the bright state spans the full Bloch sphere through the driving field parameters, this leads to full SU(2) qubit control. The technical aspects of the CPT framework and the sech-pulse control are provided in Appendix A. In this work, we are concerned with a non-ideal version of this scheme: a Λ-system with couplings to an additional, unwanted excited level.
The system under consideration is depicted in Fig. 1; the four states are comprised of a pair of low-energy levels {|0 , |1 } that encode the qubit, an auxiliary level |t which is used to mediate the qubit rotation, and an unwanted level |u . In contrast to the ideal CPT scheme, the additional excited level |u (separated by an energy splitting ε from the |t state) introduces competing offresonant couplings to the qubit states. Our goal is to perform the control of the qubit states using the target level in a way that eliminates or reduces the detrimental effect of the unwanted transition. Throughout this work, for off-resonant drivings, we take the frequency of the control pulses to be smaller than the transition frequency of the target transitions (i.e., negative detuning δ); this choice minimizes the coupling to the unwanted level [38]. Figure 2 shows the fidelity of a −π/2 rotation about the x axis through direct application of the Λ-system CPT formalism where errors are caused by off-resonant coupling to the unwanted level. In the absence of corrective measures, high fidelity can be obtained only through the use of extremely narrow bandwidth pulses (thus using extremely long pulses in the time domain). In practice however, quantum gates should be implemented within the coherence time of the system and before it relaxes to its ground states through spontaneous emission. The relaxation time of relevant solid-state quantum emitters is usually short (e.g., ∼ one ns in quantum dots [39], 1.85 ns in silicon vacancy [40], 4.5 ns in tin-vacancy [41], and 10 ns in NV centers [42]), thus requiring broad bandwidth pulses in general. We develop a formalism that deals with the issue of off-resonant couplings without trading off the duration of the gates for selectivity, ensuring operations performed within the relaxation time. The effects of spontaneous emission are further discussed in Section V A.
As mentioned above, the two legs of the Λ-system are distinct and each transition is driven by a single drive field (E 0 (t) or E 1 (t)), as shown in Fig. 1. This can be satisfied by either polarization selection rules or large energy separation of the ground states, depending on the specifics of the platform considered. In the former case, the orthogonality of each transition dipole with one drive field ensures that each transition couples to a single drive. In the latter case, sufficient energy separation of the groundstate levels implies that the off-resonant couplings of the drive fields to the opposite Λ-transitions average out. We relax this assumption and discuss the effects of the crosstalk in Section V B.
We consider two separate cases for the composition of the target and unwanted levels that lead to different physics of the system. First, the case that the two excited states, |t and |u , are formed from superpositions of two basis states, |b 0 and |b 1 : |t = sin(η)|b 1 − cos(η)|b 0 and |u = cos(η)|b 1 + sin(η)|b 0 , where |0 couples only to |b 0 , and |1 couples only to |b 1 , with dipole moments d 0,b0 and d 1,b1 , respectively. In this case, the parameter η determines the weight of the composition of the two basis states. Since |0 couples only to |b 0 , we have the Rabi frequency Ω 0 (t) ≡ − cos(η)d 0,b0 E 0 . Similarly, since |1 couples only to |b 1 the Rabi frequency Ω 1 (t) ≡ sin(η)d 1,b1 E 1 . This choice translates into a set of couplings to the unwanted level given by λ 0 (for the |0 ↔ |u transition driven by E 0 ) and λ 1 (for the |1 ↔ |u transition driven by E 1 ), that are inversely related: λ 0 ≡ − tan(η) = −1/λ 1 .
Alternatively, we could consider the case that the target and unwanted level have bare couplings to the two qubit states. In this case we can set the couplings to the target state to unity (i.e., Ω 0 (t) ≡ d 0,t E 0 and Ω 1 (t) ≡ d 1,t E 1 ) and take the couplings to the unwanted state to be proportional to λ 0 and λ 1 as the bare coupling: λ 0 Ω 0 (t) (for the |0 ↔ |u transition) and λ 1 Ω 1 (t) (for the |1 ↔ |u transition). As such, the λ 0 and λ 1 parameters (coupling strengths) will be independent of each other.
In the following we discuss the errors caused by each of the cases above and develop methods that allow implementation of fast and high-fidelity gates. We leave the numerical simulations regarding the nature of the errors in each case (that is, phase error versus leakage), until Section VI where we have developed enough methodology. In Section III, we discuss that the errors from the case of dependent couplings (λ 0 ≡ − tan(η) = −1/λ 1 ) are mostly phase errors and we can achieve quantum gates with negligible error simply by modifying the drive frequency. We discuss the case of independent couplings in which the errors are due to leakage in Section IV where we develop a novel version of the DRAG method to achieve high fidelity gates in these systems.

"
|ti |ui |Bi |Di V t,B Figure 3. (Color online) Selection rules in the dressed basis of CPT in the case of equal couplings, i.e., η = π/4. In this case the problem reduces to two dissociated two-levels (bright and target, and dark and unwanted), each subject to a transitionless pulse.

III. DEPENDENT COUPLINGS: AN EXACT SOLUTION APPROACH
In this section, we first consider the Λ-system with an additional excited state formed from basis states that lead to dependent couplings. The additional excited state induces off-resonant couplings outside of the Λ-subspace, which in the CPT frame translates into transitions that link the bright and dark states to the unwanted level. In an ideal Λ-system, one would drive the target transitions with the fields E (t) = Ω ,o (t) cos(ω t) ( = 0, 1), and static detuning δ to implement the desired gate operation. Here o refers to the original drive fields, described by the unperturbed hyperbolic secant envelopes, and ω l is the drive frequency. The resonant frequencies for the transitions |0 ↔ |t and |1 ↔ |t are denoted as ω 0t and ω 1t respectively; the detuning δ is defined as δ = ω 0 − ω 0t = ω 1 − ω 1t . In the following, as is standard for CPT, we make use of the rotating wave approximation (RWA), i.e., we define E (t) = exp(−iω l t)Ω ,o (t) + c.c. .
Although our formalism is general and applicable to the design of arbitrary axis single-qubit gates, we choose to showcase our protocols in the particular context of Xgates (or equivalently Y -gates, up to a phase between the two drives). For this reason, we set the two Rabi frequencies of the target transitions (i.e. |0 ↔ |t and |1 ↔ |t ) to be equal, that is Ω o ≡ Ω 0,o = Ω 1,o . Under this condition and RWA, the Hamiltonian in the CPT frame rotating with the drive frequency is given by (see Appendix A for derivation): where Π m = |m m|. Notice that the CPT Hamilto-nian above is generic since the relation among the λ 0 and λ 1 have been kept implicit. However, in this section we consider the case of dependent couplings: λ 0 ≡ − tan(η) = −1/λ 1 . In this case, when η = π/4 (equivalently λ 0 = −1, λ 1 = 1), the four-level system in the CPT frame transforms into two independent two-level systems, each driven by a sech pulse (Fig. 3): the bright and target, and the dark and unwanted. For a single two-level system, the population returns to the ground state at the end of the evolution if the absolute value of the Rabi frequency is equal to the bandwidth of the pulse (see Appendix A). For η = π/4 this condition is satisfied for both two-level systems. Under this condition, we can therefore solve analytically the problem as it separates into two separate two-level problems. This enables us to find an exact solution, where the unwanted level is fully taken into account and its dynamics are incorporated into the gate design. The difference in the dynamics of the two two-level systems arises from the different detunings; each ground state acquires a different phase at the end of the evolution. Thus, the total phase (rotation angle in the bright-dark subspace) will be φ tot = φ B,t − φ D,u , where φ B,t = 2 arctan(σ/δ), and φ D,u = 2 arctan(σ/(δ − ε)). Therefore the error in implementation of the desired gate is caused by the deviation of φ tot from the intended phase in the bright-dark subspace. We can account for this phase difference by writing φ tot as and then obtain the detuning modification required for the exact implementation of any desired rotation angle. To implement a rotation by φ tot = −|φ|, we find By setting the detuning to either branch of the equation above, a −|φ| rotation with negligible gate error will be implemented. For instance, R X (−π/2) and R X (−π) are illustrated in Fig 4. These gate performances are quantified by averaging over all input states existing in the Hilbert space, which can be reduced to the following expression [43]: Here, the ρ j 's are the six cardinal states on the Bloch sphere, U i (ρ j ) is the evolution of the axial vectors under the actual evolution of the system, and i is either the original or the exact solution. The ideal evolution in the subspace of {|D , |B } is given by U ideal = diag(e −iφ , e iφ ), where φ is given by Eq. (1). We will use this quantification of the gate performances throughout the rest of this paper. It is noteworthy that, as Fig. 4 shows, without pulse shaping/chirping, or increased experimental overhead, we can realize unit fidelity operations by incorporating the unwanted level into our quantum control design. As we move further away from this particular setup of coupling strengths, the two two-level systems become coupled. Nevertheless, our analytic solution still gives rise to substantial gate improvement. We define the gate improvement as the ratio of the original gate error to the gate error using the improved solution. This is shown in Fig. 5. As seen from the figure, away from η = π/4 the gate performance improves, with of course the best improvement occurring near η ≈ π/4 where the two two-level system formation of CPT frame is valid.
The two two-level systems formation does not necessarily occur in the case of independent couplings. More generally, if the two couplings have the same sign, then the errors will not be in the form of phase errors anymore. For this reason, when the couplings are independent we resort to a different strategy based on perturbative expansion methods, as we describe in the following section.  4), in terms of the dimensionless parameter bandwidth over splitting (σ/ε). The gate times of the sech pulses are chosen to be 16/σ and the non-monotonic behavior of RX (−π) is due to negligible numerical instability rising from cutting off the sech pulses. The exact approach enables gate implementation with negligible gate error.

IV. INDEPENDENT COUPLINGS: DRAG FORMALISM WITH CPT
In this section, we consider the case of a Λ-system with an additional unwanted level, where the couplings to the auxiliary states are independent of each other but have the same sign. To gain insight into how the CPTtransformed Hamiltonian (Eq. (2)) in the system would look in a limiting case, we may consider when the couplings are equal to each other and have the same strength as the coupling to the target level: λ 0 = λ 1 = 1. It can be seen from Eq. (2) that in this case the system turns into a V-type system rather than two dissociated two-level systems. The errors, in this case, will not be a Figure 5. (Color online) Gate improvement of (a) RX (−π) and (b) RX (−π/2), in terms of the dimensionless parameter pulse bandwidth over splitting (σ/ε), and the parameter η which determines the weight distributions of the basis states in the target and unwanted levels |t = sin(η) |b1 − cos(η) |b0 and |u = cos(η) |b1 + sin(η) |b0 . The best improvements occur for when the two dissociated two-level approximation is valid (η ≈ π/4). However, even away from this value, the exact solution leads to substantial gate improvements, through a single modification of the detuning.
simple phase error that can be taken care of by modification of the detuning. In fact, as we will demonstrate in Section VI, the source of errors, in this case, is leakage to the unwanted level. Therefore, we develop a novel version of the DRAG method tailored to CPT that allows the implementation of high-fidelity gate operations.
To overcome the errors for the present case via the DRAG approach, we start by modulating the original pulses. We consider an additional corrective drive Ω l,c (t) ( = 0, 1) for each of the two fields, phase detuned from the original drive by π/2. We further set the frequency to be the same as that of the original drive, hence reducing the experimental overhead of an additional laser drive. From here on, we use the letter c to refer to any subsequent parameters of the corrective drive fields. We choose the total fields of the system to be Following the previous section, we apply our formalism to the context of X-rotations, so we set the two Rabi frequencies of the target transitions to be equal, that is Ω o ≡ Ω 0,o = Ω 1,o , with the additional condition of Ω c ≡ Ω 0,c = Ω 1,c . However, we note that these conditions can be lifted and the formalism we develop remains valid, with the difference that the rotation axis changes.
Under these conditions and after performing the RWA, the Hamiltonian in the CPT frame is given by (see Ap-pendix A) where to indicate the matrix elements for the generic transition |m ← → |n , we have defined σ o m,n ≡ |n m| + |m n| and σ c m,n ≡ i |n m|−i |m n|. Previous formulations of the DRAG method have focused on canceling out leakage errors in ladder-type systems (e.g., transmons) that occur between consecutive energy levels, making it inapplicable to Λ-systems. The complexity, in this case, arises as the qubit control is performed indirectly via the target level. Therefore, to mitigate the unwanted transitions of Λ-type structures we modify the DRAG method as we explain below.
Under the DRAG formalism, corrections to the driving fields are determined by utilizing a frame transformation [35]. In this new DRAG frame, the resulting Hamiltonian is constrained such that it implements the target evolution. The DRAG frame Hamiltonian for a generic Hamiltonian H(t), is generated by the transformation A(t) = e −iS(t) and is given by Ω(t) is Ω(t) = σ sech(σ(t − tg/2)), and its derivative corrective solution (in orange) modulated with the splitting and the couplings is . The bandwidth is taken to be σ = 0.02 meV, and gate time is tg = 16/σ. In our case, H(t) ≡H CPT (t). The operator S(t) can be an arbitrary Hermitian operator as long as it respects the boundary conditions of the transformation. That is, the frame transformation has to vanish at the beginning and end of the pulse (A(0) = A(t g ) = 1), such that the ideal gate we wish to design remains the same in both the CPT and DRAG frames. The generality of S(t) would in principle generate a wide range of pulse corrections. However, extracting such closed-form expressions is non-trivial for our four-level system (as the pulses are time-dependent), hence, we turn to a perturbative expansion of the transformation. To that end, we utilize the Schrieffer-Wolff (SW) transformation [44] and its perturbative expansion.
Additionally, we expand the CPT Hamiltonian of Eq. (6) into a power series. To do so, we first make the CPT Hamiltonian dimensionless by multiplying all quantities by the gate time t g [34]. The Hamiltonian is then expanded into a power series (Appendix B) of x = 1/εt g , where ε is the energy cost by which the unwanted subspace of |u is off-resonant with respect to the Λ-system. This way, one can collect the same orders of the expansion in the left-hand side and right-hand side of Eq. (7) and relate the elements of H DRAG to the corrective fields (involved inH CPT ) and to elements of S(t) order by order. Further details of this procedure are given in Appendix B.
Our goal is to constrain the DRAG-frame Hamiltonian such that it implements our target evolution. To this end, we define a target Hamiltonian, capable of performing arbitrary rotations within the qubit (dark-bright) subspace, where h i (t) and h z (t) are arbitrary control fields. Here we choose the h i (t) to be sech pulses, where h z (t) corresponds to the detuning. To ensure that the DRAG frame Hamiltonian implements the intended operation as dictated by H CPT target , we impose the following constraints: where i ∈ {o, c}. Notice that these constraints ensure that the zeroth order DRAG Hamiltonian remains the same as the ideal Hamiltonian of the CPT frame. Additionally, to enforce decoupling from the unwanted subspace in the new Hamiltonian, we impose the following constraints: where i ∈ {o, c}. The constraints can be solved consistently and lead to a set of corrective pulses for each order of the expansion. The details of the derivation for the corrective fields can be found in Appendix B. To the first order of expansion, for the two drives being subjected to an uncorrected pulse Ω(t), the simplest solution which respects the boundary conditions of the transformation is given by We depict the sech pulse envelopes for one specific set of parameters in Fig. 6. Due to the indirect control of the qubit in Λ-systems, the diagonal constraint of Eq. (9), unlike what happens in transmon qubits, does not lead to a global phase among all states for the first-order DRAG Hamiltonian. Instead, it resolves the off-resonant coupling at the cost of inducing a phase between the bright and dark states. Therefore, we need to investigate the form of the firstorder DRAG Hamiltonian to infer this phase between the dark and bright states. By restricting our attention to the Λ-system subspace, we find the zeroth-, and first-order DRAG frame Hamiltonian in the subspace of {|D , |B , |t } to be and Notice that for the limiting case of a V-system (i.e., λ 0 = λ 1 = 1), all elements of H (1) DRAG in Eq. (14) vanish, except for the diagonal entry of the bright state. While the zeroth-order DRAG Hamiltonian stays in the form enforced by Eqs. (9), we note that in the first-order DRAG Hamiltonian there is an induced phase between the dark state and the bright-target subspace. In this Hamiltonian (Eq. (14)), the bright and target states follow the same phase evolution, as dictated by the common diagonal element −1/(16ε)(λ 0 + λ 1 ) 2 Ω 2 (t). However, the qubit states are composed of the dark and bright states, which up to the first order (neglecting the off-diagonal elements of H (1) DRAG ) evolve with a different phase. Effectively, this implies that the DRAG corrections reduce the unintended off-resonant couplings to the unwanted level at the cost of inducing a relative phase between the qubit states in the CPT frame. This is an immediate consequence of the fact that we counteract the unwanted couplings indirectly; in the CPT frame, we design a target Hamiltonian that involves transitions between the bright and target levels.
Let us discuss the validity of the approximation above. As discussed, the unwanted subspace in Eq. (14) is neglected, since these elements constitute a higher-order error. This arises from the fact that in Eq. (7), the n-th order DRAG Hamiltonian includes contributions from S (n+1) (t). That is, the zeroth-order constraints ensure that H (0) DRAG ≡ H CPT target . The first-order target constraints give us the solutions for the pulse corrections (note that H (1) DRAG includes contributions from S (2) (t), which for the purpose of this analysis we set arbitrarily as S (2) (t) = 0). We note that the first-order target constraints do not mix with higher elements of S (2) (t), which is not the case for the first-order decoupling constraints. Terminating at n = 1 target constraints means that we leave secondorder errors in the unwanted subspace of H DRAG . If the goal is to find second-order pulse corrections, then one has to impose the second-order decoupling constraints, which will constrain certain elements of S (2) and disconnect the unwanted subspace of H (1) DRAG from the qubit subspace. Then, applying the second-order target constraints on H (2) DRAG gives rise to the second-order pulse corrections, by leaving 3rd order unwanted subspace errors in H (2) DRAG . This procedure can be repeated to obtain any subsequent order of pulse corrections as has been shown in Ref [34] for a transmon system. In the transmon case, the off-diagonal entries in the qubit subspace can be made to be zero in H DRAG . In the Λ-system case, due to the different selection rules, we find that there is no transformation that can make the off-diagonal qubit subspace elements to be zero and simultaneously respect the boundary conditions of the transformation. Hence, attempting to go to higher-order pulse corrections will leave non-zero off-diagonal entries in the qubit subspace and different diagonal phase evolutions, which are difficult to incorporate analytically into the gate design. Consequently, we terminate the search of pulse corrections at the first order, and we show later on that such pulses can still lead to enhanced gate performance, without demanding higher-order modifications. The gate improvement is possible by resolving only the different phase evolution of the dark and bright-target subspaces, which makes our approximation for deeming the off-diagonal qubit subspace entries to be higher-order errors valid.
We now highlight the procedure for finding the required detuning modification, on top of the pulse modulation, to correct for the phase error caused by the implementation of DRAG. We showcase this for the case of λ 0 = λ 1 = 1 because, as we discussed above, it only leads to an induced phase between the bright and dark states up to the first order. To this end, we denote the ideal evolution operator that corresponds to H which reduces to: where U is the evolution operator of H is the function that defines the relative phase shift between the dark and bright states. Solving the Schrödinger equation we find that the induced phase is given by θ = −(λ 2 0 + λ 2 1 − 6λ 0 λ 1 )σ/(4ε). Hence, in order to remove the θ-shift from the target evolution we modify the detuning: where φ is the target rotation angle. The detuning modification (17) together with Eqs. (11) and (12) complete our full set of solutions for the pulses. The final pulse shapes will depend on the specific details of the system such as the splitting of the two auxiliary states.  We demonstrate the performance of the DRAG solutions in Fig. 7. The figure shows the fidelity of rotations about the x axis, i.e. R X (−π) and R X (−π/2) gates, in terms of the dimensionless parameter given by the pulse bandwidth over the splitting between the target and unwanted levels. We consider the case in which both transitions are of equal strength larger than that with respect to the target level (λ 0 = λ 1 = 1.2) and the case of unequally-strong couplings with one larger and one smaller than the coupling to the target state (λ 0 = 0.8, λ 1 = 1.2). The top panels of Fig. 7 show the gate error of implemented operations and the bottom panels shows the deviation from unitarity of the operators, defined as |U U † − 1|. It is clear from these plots that our DRAG correction provides substantial improvement compared to the uncorrected gate, in some cases by several orders of magnitude. It can also be seen that the operators do not remain unitary for the whole range of bandwidths over splittings, which indicates that the populations are leaving the qubit subspace. This signals that the sources of errors in this case are leakage, rather than phase errors, unlike the case of dependent couplings in Section III. As such, our devised DRAG solution performs better at ranges of bandwidths over splittings where the gate implementation remains unitary, but it always provides an improvement throughout the full bandwidth range.
In Fig. 8, we show the improvement in the fidelity of the R X (−π/2) gate as a function of both the bandwidth over splitting and the ratio of the unwanted couplings for three different values of these couplings with respect to the coupling to the target state (which is set to unity in all cases). In all cases, the best DRAG improvements occur at points where the values of the couplings are similar to each other, i.e. λ 0 = λ 1 . The occurrence of minima near these values are related to the corresponding deviation from unitarity values: at these specific parameter values of the system, the transitionless condition is being satisfied and therefore the DRAG decoupling from the unwanted state becomes more efficient. Furthermore, notice that because our pulse corrections are functions of the couplings, improvements are better for the cases where the unwanted couplings are stronger compared to the coupling to the target state (i.e., λ 0(1) > 1). That is, for stronger couplings the corrections to the pulse shape will be more prominent and therefore DRAG becomes more effective. Consequently, in this formalism, we always pick the auxiliary state that has weaker dipole moments to be the target state.

V. EFFECTS OF ADDITIONAL ERRORS
In this section, we consider the effects of errors other than the off-resonant couplings to the unwanted level to examine how the performance of our solutions is affected. In particular, we focus on two important sources of error: spontaneous emission from the excited levels and crosstalk between the transitions of the Λ-system.

A. Spontaneous emission
A realistic optically active system is coupled to the environment and therefore relaxes through the spontaneous emission of photons. This relaxation process reduces the fidelity of our intended quantum gates, more severely impacting longer gates. We describe the dy-namics with a standard open quantum system approach to model the effects of spontaneous emission. We use the Liouville-von Neumann equation with Lindblad relaxation terms: The Lindblad operators account for the populations that leave the target and unwanted levels upon emission of a photon. We pick four Lindblad operators that correspond to emission from each excited state to each ground state: L ij = √ γ |i j|, where i ∈ {0, 1}, j ∈ {t, u}, and γ is the emission rate, which we take to be the same for all transitions.
We have considered the effects of spontaneous emission for both cases of dependent and independent couplings. Fig. 9 shows the effects of spontaneous emission on the corrected exact and DRAG solutions for different dimensionless values of γ/ε, for R X (−π/2). As evident from the figure, we see that two opposing effects are acting on the system. On one hand, the detrimental effect of the unwanted level favors short-bandwidth pulses. On the other hand, the spontaneous emission favors largebandwidth pulses (i.e. short pulses in time). This results in a minimum gate error for an intermediate value of the pulse bandwidth. (In the DRAG case, notice that this is independent of deviation from unitarity condition discussed in the previous section.) The insets in Fig. 9 show the maximum obtainable value of fidelities with their corresponding values of bandwidth over splitting. For both values of γ/ε, while the corrective solutions lead to a better maximum fidelity, they also occur at larger values of bandwidth over splitting. Therefore our formalism provides us with better gate performances, at more reasonable values of pulse bandwidths, suitable for systems with fast spontaneous emission rates.

B. Cross-talk
So far in our discussion we have assumed that either the polarization selection rules or the energy separation of the ground states allow to distinguish the two transitions of the Λ-system. Here we include the error from cross-talk of the two transitions. We model our interaction picture cross-talk Hamiltonian (for an X-rotation) in the frame rotating with the frequency of the pulses, as follows: Here ω s quantifies the ground state (GS) splitting (i.e., the splitting between the qubit states, |0 and |1 , in the lab frame). Following the assumption on the transition dipoles discussed in Section II, the cross-talk coup- With no corrective measures, best fidelities occur at narrow bandwidths. Our modified pulses lead to better fidelities at higher pulse bandwidths. Figure 10. (Color online) Gate improvement of RX (−π/2) for (a) the exact method and (b) using the DRAG method, including the cross-talks between the two transitions of the Λ-system. The vertical axis is the dimensionless parameter pulse bandwidth over the excited states splitting (σ/ε), and the horizontal axis is the dimensionless parameter ground states splitting over the excited states splitting (ωs/ε). The other parameter values are η = π/4 for the exact method, and λ0 = 1, λ1 = 1 (κ01 = 1 = 1/κ10) lings may also be implicitly dependent or independent of each other. When the E 0 (E 1 ) field drives the |1 → |t (|0 → |t ) transition, it leads to cross-talk coupling κ 01 = E 0 /E 1 (κ 10 = E 1 /E 0 ). Consequently, once the E 0 field drives the |1 → |u transition, the resulting crosstalk coupling will be λ 1 κ 01 , and similarly, the effect of E 1 driving the |0 → |u transition will be the coupling λ 0 κ 10 . Fig. 10 shows the effect of cross-talk on the gate improvement of R X (−π/2) for both the exact and the DRAG method. The values are shown for a range of dimensionless parameters of pulse bandwidths over ex- Figure 11. (Color online) RX (−π/2) gate improvement with cross-talks in the regime that ε ωs for (a) the exact method and (b) the DRAG method. In addition to the corrective measures, the Rabi frequencies are also multiplied by 1/2 due to the cross-talk modifications to the system Hamiltonian (see the text). The other parameter values are η = π/4 for the exact method, and λ0 = 1, λ1 = 1 (κ01 = 1 = 1/κ10). cited state (ES) splittings (σ/ε), and GS splitting over ES splitting (ω s /ε). In (a), for the exact method, we have set η = π/4 and in (b), for the DRAG method, we have set λ 0 = 1, λ 1 = 1. Furthermore, the ratio of the two electric fields is chosen to be unity E 1 /E 0 = 1 (i.e., κ 01 = 1 = 1/κ 10 ). As the figure indicates, the behavior of the improvement heavily relies on the specific parameters of the system. For example, the cross-talk Hamiltonian corresponding to Fig. 10(a) (i.e., λ 0 = −1, λ 1 = 1) in the CPT frame after RWA will be, This indicates that, firstly, the pulse profiles will be modulated with some time-dependent, periodic functions of ω s , therefore the pulse behavior will be non-trivial.
Additionally we notice that for ε ω s , the system reduces to that of two two-level systems. This means that in this regime our exact solution is applicable. Similarly, looking at the corresponding Hamiltonian to Fig. 10(b) (i.e., λ 0 = 1, λ 1 = 1) we find that in the same ε ω s limit, it reduces to that of the V-system for which we developed the DRAG formalism. Both of these cases require a slight correction of Rabi frequencies by a factor of 1/2 to account for additional factor of 2 arising from the terms 1 + cos(ω s t). We show the performance of these solutions in the presence of cross-talk in Fig. 11. As seen in both figures, improvements of several orders of magnitude are achievable through our approaches given that the condition ε ω s is satisfied. We also note that in the case of dependent couplings, because the errors are in the form of phase errors, in principle it is possible to track the phase change and correct it through detuning modification. Moreover, for a generic DRAG-based approach to avoid the cross-talks in Λ-systems we refer to the work in Ref. [45].

VI. DISCUSSION AND CONCLUSION
In this work we have developed a framework for highfidelity control of Λ-systems with unwanted transitions for two cases of dependent and independent couplings to the unwanted level. The case of the dependent couplings is motivated by the fact that off-resonant couplings are likely to arise due to the physical origin of the excited states as orthogonal superpositions of two basis states. For example, in quantum dot molecules where the qubit states are the spinor projection of a single hole, both the |t and |u states are created by the bonding and antibonding superpositions of two different charge and spin configurations, with the coupling mediated by tunneling.
In such an example, the degree of mixing can be controlled by applied electric fields [30]. We showed that in this case the Λ-system reduces to a combination of a two two-level-systems in the limit of equal basis weights (η = π/4). We further discussed that the errors are in the form of phase errors; as such we devised an exact solution that compensates for these errors in the form of detuning modification and allows for error-free gate implementations. On the other hand, in the absence of basis states, the system turn into a V-system in the limiting case of equal couplings with the same sign. In this system the source of errors has the form of leakage rather than phase errors. In Section IV we developed a novel version of the DRAG approach to battle the leakage in this case.
To see the different form of errors (i.e., phase error versus leakage), we briefly comment on the evolution of population of each level during the gate time. In Fig. 12 we have shown the populations of the system for a R X (−π/2) rotation, for the cases of both dependent and independent couplings. The transitionless nature of the pulses driving two independent two-level-systems is evident in Fig. 12(c): the final population of excited states are negligible for both uncorrected and corrected solutions. This indicates that the core of the errors in our system are the phase errors due to the off-resonant couplings to unwanted states, rather than leakage. Our formalism corrects these errors by populating the excited states appropriately and redistributing the populations on the qubit states such that the desired gates are implemented. On the other hand, as we discussed in Section IV, in the independent couplings case the evolution of the system is nonunitary. As such, by the end of the gate time some of the population leaves the qubit subspace, being stored in the unwanted level ( Fig. 12(d)). Our DRAG approach battles this leakage through modulation of the pulses and restores the full population to the qubit subspace. However, this comes at the price of inducing some phase errors which are compensated through the additional detuning modification.
In summary, the two methods we present above are capable of implementing high fidelity arbitrary gates in either case of a Λ-system with unwanted transitions. The exact method has the advantage that it only requires a simple modification of the detuning and therefore is less expensive from an experimental point of view. The DRAG pulse modification, on the other hand, is in the form of a corrective modification to the original pulse shape that drives each transition of the system.
It should be noted that the DRAG pulse modulation solution is not unique; we have made specific choices on the generative parameters of the DRAG frame in deriving our DRAG solutions (see Appendix B for details of the derivations). These generative elements are the free parameters of the system and can be set to arbitrary values as long as they satisfy the DRAG transformation condition S(t g ) = S(0) = 0. Therefore, alternative DRAG solutions based on different generative elements can also be devised. We have demonstrated such pulse shapes in Fig. 6. The corrective pulse is inversely proportional to the splitting ε and the coupling to the unwanted level for that transition (Eq. (12)). For smaller splittings, the correction pulse becomes more comparable to the original pulse. We note that even though the mathematical derivation of the pulse design is phrased in terms of an original (uncorrected) pulse and a correction, experimentally the full corrected pulse would be programmed and generated directly.
Finally, we discussed the implementation of X-, and Y -rotations. To achieve universal control of the qubit system, we also require implementations of Z-rotations. This can be done either through the same control scheme by driving only a single transition of the Λ-system with a sech pulse [2] or simply through 'virtual' Z gates [46]. The exact solution and the version of DRAG we have tailored to the case of a Λ-system with a fourth unwanted level are general and can be applied to a variety of optically active qubits, such as color centers (e.g., the NV center in diamond), trapped ions, rare-earth ions, selfassembled quantum dots, and quantum dot molecules. In this appendix, we present the mathematical details of the CPT framework. Two-level systems subject to a sech pulse with Rabi frequency Ω and bandwidth σ, can be solved analytically [37] and the solutions are in the form of hypergeometric functions. For the case of Ω/σ ∈ N, these pulses are transitionless [2], i.e., after the passage of the pulse the population will always return to the ground state with the ground state acquiring a non-trivial phase φ (Eq. (1)) through the process. For the specific case of Ω = σ, the Hamiltonian of a generic two-level system driven by a sech pulse in the rotating frame is given by where Ω(t) = σ sech(σ(t − t g /2)). For this system the unitary evolution operator at the end of the gate is U = diag(e −iφ , e iφ ). In the CPT scheme, the transitions of the system are excited using the drive field, E 0 f 0 (t)e iω0t + e iα E 1 f 1 (t)e iω1t + h.c. For identical temporal envelopes and detunings, with Rabi frequencies Ω 0 (t) and Ω 1 (t), the transformation of the original qubit states to the bright and dark states is given by: where sin θ 2 = Ω 0 /Ω eff , cos θ 2 = Ω 1 /Ω eff , and Ω 2 eff = Ω 2 0 + Ω 2 1 . In the CPT frame, the transition matrix elements between the target and the dark state vanishes and the bright and target state will have the matrix element defined by the effective Rabi frequency: V t,B = Ω eff f (t)e −iδt . For the case of both drives using a sech temporal envelope, i.e., f 0 (t) = f 1 (t) = sech(σt), the transitionless pulse with Ω eff = σ will induce the relative phase φ between the bright and dark states which translates to a rotation in the subspace of |D and |B . Therefore, by varying the drive parameters we can set the unitary transformation of Eq. (A2) to transform our original qubit states to the desired states in the CPT frame, effectively enabling the rotation about an arbitrary axis of rotationn = (sin θ cos α, sin θ sin α, cos θ): R n (φ) = e −iφn. σ .
The CPT framework can be applied to the case of Λsystem with unwanted level in a similar way. However, there will be additional transitions from the bright and dark state to the unwanted level. In the lab frame of a non-ideal system with an unwanted level, for the case of equal detunings δ, the Hamiltonian in the interaction frame after the RWA can be written as    From these constraints, we first find the zeroth-order solutions. According to the CPT framework, in order to implement X-rotations we set h   are set to zero, such that we obtain the desired evolution as dictated by the target Hamiltonian. Since we made our choice for the control fields h(t), by making use of H (0) extra = 0 we can now solve for the pulse envelopes, from which we find: This ensures that we have the right form in the brighttarget subspace of H DRAG to match the corresponding bright-target elements of the target Hamiltonian. (Notice that the additional factor of √ 2 is present due to the CPT frame transformation.) Next, we proceed with the zeroth-order decoupling constraints. These constraints together with the zeroth-order target constraints will fix certain elements of S (1) (t), such that we essentially satisfy H 2 (λ 0 + λ 1 )t g Ω(t). These two constraints ensure no transitions between darkunwanted and and bright-unwanted states, respectively.
The first order corrective controls are found by satisfying the first-order target constraints. In this first order, as we already mentioned we set all the first-order control fields to zero. Making use of Eqs. (B5), we find the following equations: To seek the simplest solution which satisfies the S (1) (0) = S (1) (t g ) = 0 (such that the implemented gate is the same in the DRAG and CPT frame), we pick the free parameters s (1) z,i 's and s (1) o,B,t all equal to zero. Substituting s (1) c,B,t into the first two equations, we find the pulse corrections given in Eqs. (11) and (12) (notice that we require an additional factor of √ 2 in the non-CPT frame to accommodate for the CPT transformation). Substituting these solutions and the choices we made above for the generative elements of S (1) (t) in the expansion Eq. (B1) will lead to the first-order form of the DRAG Hamiltonian H