Abstract
Strong spin-orbit interactions make hole quantum dots central to the quest for electrical spin qubit manipulation enabling fast, low-power, scalable quantum computation. Yet it is important to establish to what extent spin-orbit coupling exposes qubits to electrical noise, facilitating decoherence. Here, taking Ge as an example, we show that group IV gate-defined hole spin qubits generically exhibit optimal operation points, defined by the top gate electric field, at which they are both fast and long-lived: the dephasing rate vanishes to first order in the electric field noise along with all directions in space, the electron dipole spin resonance strength is maximized, while relaxation is drastically reduced at small magnetic fields. The existence of optimal operation points is traced to group IV crystal symmetry and properties of the Rashba spin-orbit interaction unique to spin-3/2 systems. Our results overturn the conventional wisdom that fast operation implies reduced lifetimes and suggest group IV hole spin qubits as ideal platforms for ultra-fast, highly coherent scalable quantum computing.
Similar content being viewed by others
Introduction
Quantum computing architectures require reliable qubit initialization, robust single-qubit operations, long coherence times, and a clear pathway towards scaling up. Solid-state platforms are supported by the well-developed solid-state device industry, with mature microfabrication and miniaturization technologies. Among solid-state platforms, semiconductor quantum dot (QD) spin qubits have been actively pursued1, with an energetic recent focus on hole spins in diamond and zincblende nano-structures2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26.
The primary motivation for this focus is the strong hole spin-orbit interaction, which enables qubit control via electron dipole spin resonance (EDSR), making quantum computing platforms faster, more power-efficient, and easier to operate27,28,29,30,31,32,33,34. Electric fields are much easier to apply and localize than magnetic fields used in electron spin resonance. Only a global static magnetic field is required to split the qubit levels. In addition, the p-symmetry of the hole wave function causes the contact hyperfine interaction to vanish, and complications involving valley degrees of freedom are absent35,36,37,38,39. Initial studies indicate that hole spins may possess sufficiently long coherence times for quantum computing40,41,42,43,44. Meanwhile, much progress has been made in the initialization and readout of hole spin qubits8,11,14,17,45,46.
The existential question that will determine the future of hole QD spin qubits is: Does the strong spin-orbit interaction that allows fast qubit operation also enhance undesired couplings to stray fields such as phonons and charge noise leading to intractable relaxation and dephasing? In this paper, we demonstrate theoretically that this is emphatically not the case for hole spin qubits in group IV materials taking Ge as the most prominent example.
In fact, the unique properties of the hole Rashba interaction overturn the conventional understanding of qubit coherence in spin-orbit coupled systems, which states that, as long as the qubit is described by an effective 2 × 2 Hamiltonian, holes behave in the same way as electrons. That is, given that the EDSR rate is linear in the spin-orbit strength, while the relaxation and dephasing rates increase as the square of the spin-orbit strength, the Rabi ratio can be enhanced by operating the qubit at points at which the spin-orbit interaction is weaker27,28. In contrast group IV hole qubits achieve the best coherence when the electrical driving rate is at its maximum. In all other systems, one has to choose between creating long-lived qubits by isolating them from the environment and accepting slower gate times, or designing fast qubits strongly coupled to the environment, but which decohere rapidly.
The key realization is that holes in group IV materials are qualitatively different from group III–V materials. They have tremendous potential for qubit coherence, with Ge and Si possessing isotopes with no hyperfine interaction, as well as a near-inversion symmetry that eliminates piezo-electric phonons. This near-inversion symmetry will eliminate the Dresselhaus interaction, the interface inversion asymmetry terms are expected to be negligible in the system47,48. Which enables spin resonance in group III–V materials27. On the other hand, we show that strong cubic-symmetry terms enable a fully-tunable Rashba coupling. Unlike the Dresselhaus interaction, and unlike electron systems, the hole Rashba interaction evolves nonlinearly as a function of the gate electric field, a fact traced to the hole spin-3/2. The qualitative difference between the Rashba and Dresselhaus interactions for holes is vital for qubit coherence. Thanks to this nonlinearity, dephasing due to electric field fluctuations in all spatial directions can be essentially eliminated at specific optimal operation points defined by the gate electric field15,16,38,49,50,51,52,53,54. At these points electrical qubit rotations are at their most efficient, with the spin resonance Rabi gate time at a minimum. The relaxation rate due to phonons can be made as small as desired at small magnetic fields of the order of 0.1 T, which allows 106 − 107 operations in one relaxation time for an in-plane alternating field EAC ~ 103 Vm−1.
Our focus in this paper is on Ge, which has witnessed enormous recent progress6,7,9,10,55. Holes in planar Ge quantum wells have a very large out of plane Landé g-factor, g ≈ 20, enabling operation at very small magnetic fields, which would not impede coupling to a superconducting resonator. The low resistivity of Ge when contacting with metals makes couplings between other devices such as superconductors easier9,56,57. In the past decade some spectacular results have been reported, for example, EDSR detection techniques18,55, structures of quantum confinement systems7,58,59,60, the anisotropy of g-tensors10,58, spin-orbit couplings, and transport phenomena in two-dimensional hole systems4,6,10,61,62. We focus a on single dot throughout this work. A prototype device, including a neighboring dot, is shown in Fig. 1. The Hamiltonian describing a single hole quantum dot has the general form H = HLK + HBP + HZ + Hph + Hconf, where HLK represents the Luttinger-Kohn Hamiltonian, HZ is the Zeeman interaction between the hole and an external magnetic field, and Hph the hole-phonon interaction. Hconf is the confinement potential including the vertical and lateral confinement. The vertical confinement is achieved by applying a gate electric field Fz in the growth direction, leading to a term eFzz in the Hamiltonian; the lateral confinement is modeled as an in-plane parabolic potential well. The Bir-Pikus Hamiltonian HBP represents strain18, which appears naturally as part of the quantum well growth process. A typical configuration of holes in Ge is achieved by growing a thin strained Ge layer (usually about 10 nm to 20 nm) between SiGe layers such that, if the barrier between the two layers is high enough, a quantum well can be formed. We consider SixGe1−x, where x = 0.154,6,7,58, with other values of x discussed in Supplementary Table 1.
We start from the bulk band structure of holes as derived by Luttinger and Kohn63. The spinor basis is formed by the eigenstates of Jz, \(\left\{\left|+\frac{3}{2}\right\rangle ,\left|-\frac{3}{2}\right\rangle ,\left|+\frac{1}{2}\right\rangle ,\left|-\frac{1}{2}\right\rangle \right\}\). For a 2D hole gas grown along \(\hat{z}\parallel (001)\), we write the Luttinger-Kohn Hamiltonian as:
where \(P=\frac{{\hslash }^{2}{\gamma }_{1}}{2{m}_{0}}\left({k}^{2}+{k}_{\,\text{z}\,}^{2}\right)\), \(Q=-\frac{{\hslash }^{2}{\gamma }_{2}}{2{m}_{0}}\left(2{k}_{\,\text{z}\,}^{2}-{k}^{2}\right)\), \(L=-\sqrt{3}\frac{{\hslash }^{2}{\gamma }_{3}}{2{m}_{0}}{k}_{-}{k}_{\text{z}}\), \(M=-\frac{\sqrt{3}{\hslash }^{2}}{2{m}_{0}}\left(\overline{\gamma }{k}_{-}^{2}+\delta {k}_{+}^{2}\right)\), and m0 is the free electron mass, γ1, γ2, γ3 are Luttinger parameters which are determined by the band structure. The in-plane wave vector will be \({k}^{2}={k}_{x}^{2}+{k}_{y}^{2}\), k± = kx ± iky. In this manuscript, we use the symmetric gauge A = (B/2)(−y, x, 0) so that k → − i∇ + eA/\(\hbar\) . We have also used \(\bar{\gamma }=({\gamma }_{2}+{\gamma }_{3})/2\) and δ = (γ3 − γ2)/2 to simplify the algebra. In Ge \(\delta /\bar{\gamma }\,<\,0.15\), hence δ can be treated perturbatively, while bulk Dresselhaus terms are absent. Although interface inversion asymmetry terms with the same functional form may exist47, at the strong gate fields considered here will be overwhelmed by the Rashba interaction and are not discussed in detail. The diagonal terms of HBP in the HH manifold are \({P}_{\varepsilon }+{Q}_{\varepsilon }=-{a}_{\text{v}}\left({\varepsilon }_{xx}+{\varepsilon }_{yy}+{\varepsilon }_{zz}\right)\), while in the LH manifold they are \({P}_{\varepsilon }-{Q}_{\varepsilon }=-({b}_{\text{v}}/2)\left({\varepsilon }_{\text{xx}}+{\varepsilon }_{\text{yy}}-2{\varepsilon }_{\text{zz}}\right)\), where av = −2 eV and bv = −2.3 eV are deformation potential constants18. In our chosen configuration εxx = εyy = −0.006, the minus sign indicates that the germanium is compressed in xy-plane. In the \(\hat{\,\text{z}\,}\)-direction, the Ge layer will be stretched, and εzz = (−2C12/C11)εxx = 0.0042, with C12 = 44 GPa, C11 = 126 GPa for Ge. The diagonal terms of the strain-relaxed barrier configuration will change the HH-LH energy splitting by a constant, which is approximately 50 meV.
The growth direction provides the spin quantization axis, with the heavy hole states (HHs) representing the ±3/2 angular momentum projection onto this axis, while the light hole states (LHs) represent ±1/2. In 2D hole systems, the HHs are the ground state7,27,51,64,65. To define a quantum dot a series of gates are added on top of the 2D hole gas confinement, as in Fig. 1, and we ultimately seek an effective Hamiltonian describing the two lowest-lying HH states in a quantum dot. Since we expect the HH-LH splitting to be much larger than the quantum dot confinement energy, we proceed with the standard assumptions of k ⋅ p theory, retaining at first only terms containing kz, with kx and ky initially set to zero. This determines the approximate eigenstates ψH,L(z) corresponding to the growth-direction. These are described by two variational Bastard wave functions ψH and ψL48,66,
where the dimensionless variational parameters βH,L are sensitive to the gate electric field due to the term eFzz, and d is the width of the quantum well in the growth direction, which is an input parameter. The orthogonality of the HH and LH states is ensured by the spinors. This wave function is suitable for inversion layers, as well as accumulation layers, although our focus will be primarily on the latter. For inversion layers, the Bastard wave function will also be appropriate, because in experiments the electric field can be made large enough such that the hole gas sticks around the top of the quantum well.
In the xy-plane, we model the confinement potential using a harmonic oscillator
where mp = m0/(γ1 + γ2) is the in-plane effective mass of the heavy holes, the subscript ∥ refers to the xy-plane, ω0 is the oscillator frequency, a0 the QD radius which satisfies, \(a^{2}_{0}=\hbar/(m_{\mathrm{p}}\omega_{1})\) i.e., a magnetic field will narrow the QD radius. The solutions are the well-known Fock-Darwin wave functions \(\left|{\phi }_{{n}_{1},{n}_{2}}\right\rangle\)18 with eigen-energies \({\varepsilon }_{{n}_{1},{n}_{2}}=\hslash ({n}_{1}+{n}_{2}+1){\omega }_{l}+\frac{1}{2}\hslash ({n}_{2}-{n}_{1}){\omega }_{\text{c}}\), where \({\omega }_{\text{l}}=\sqrt{{\omega }_{0}^{2}+{\omega }_{\text{c}\,}^{2}/4}\), ωc = eB/mp is the cyclotron frequency. The Bastard wave functions account for the perpendicular confinement, in the \(\hat{{\bf{z}}}\)-direction. The Fock-Darwin wave functions \(\left|{\phi }_{{n}_{1},{n}_{2}}\right\rangle\) account for the quantum dot confinement in the xy-plane. The Bastard and Fock-Darwin wave functions would be formally the same for electrons, while in a hole gas is that separate Bastard wave functions are required for the heavy and light holes.
Finally, the hole-phonon interaction is67,68,69:
where q is the phonon wave vector, Vc is the unit cell volume, NVc is the crystal volume, \({\hat{e}}_{\text{s}}\) is the polarization direction vector. The density of the material is denoted by ρ, Dα,β represents the deformation potential matrix, and \({\hat{a}}^{\dagger }\) and \(\hat{a}\) are the phonon creation and the annihilation operators. The details are presented in the Supplementary Information Eq. (1)36,70,71,72,73,74.
Our approach is semi-analytical. To incorporate the contributions from the LHs, and the Rashba spin-orbit couplings, we start from the 4 × 4 Luttinger Hamiltonian and project it onto the product states of the out-of-plane sub-bands HH1, LH1, each with two spin projections, and the first four orbital levels of the in-plane confinement, such that our Hamiltonian matrix is 40 × 40. This refers to states \(\left\{\left|\frac{3}{2},S\right\rangle \otimes \left|{\phi }_{n,m}\right\rangle \right\}\), where S = 3/2, −3/2, 1/2, −1/2 represent the HHs or LHs, n = 0, 1, 2, 3 and m = 0, ±1, 2, ±3 denote the in-plane Fock-Darwin states. We have checked that addition of the HH2, LH2 sub-bands does not modify the results, this is attributed to the significant energy gaps separating them from HH1, LH1. Given the high computational cost of adding these sub-bands, we have not taken them into account in the results presented here. To obtain the matrix elements required for the dephasing and relaxation times, as well as for the EDSR Rabi frequency, we perform a 3rd order Schrieffer-Wolff transformation on the 40 × 40 Hamiltonian. This transformation takes into account all the spin-orbit terms that are not separable in the spatial coordinates, which are precisely the terms leading to the Rashba interaction. To ensure the accuracy of the Schrieffer-Wolff method we compare the results for the Larmor frequency with a full numerical diagonalization of the 40 × 40 matrix.
Results and discussion
Qubit Zeeman splitting
The qubit Larmor frequency has been plotted in Fig. 2 as a function of the gate electric field. The Schrieffer-Wolff method agrees well with the numerical diagonalization: the location of the optimal operation point differs by only 2% in the two approaches. We note the non-monotonic behavior as a function of the gate field, leading to an optimal operation point in the range of 30–50 MV/m. Electric fields of such magnitude are used routinely in quantum computing experiments7,46,75,76. The non-monotonic behavior is directly related to the behavior of Rashba spin-orbit coupling discussed below. We note that the spatial dimensionality of the qubit is determined by the relevant energy scales, namely the \(\hat{{\bf{z}}}\)-sub-band energy spacing compared to the energy splitting of the lateral wave functions. The heavy hole– light hole splitting, given by the perpendicular confinement, is many times larger than the in-plane qubit confinement energy, determined by the in-plane confinement, so the system is in the quasi-2D limit. Nevertheless, our findings, such as trends with the top gate field, can be interpreted qualitatively by analogy with the Rashba interaction in the asymptotic 2D limit d → 077, which we also obtain from Eq. (1). For a system with cubic symmetry this contains two terms with different rotational properties:
where \({\hat{\sigma }}_{\pm }\equiv ({\hat{\sigma }}_{\text{x}}\pm i{\hat{\sigma }}_{\text{y}})/2\). The coefficients are evaluated as:
where EH and EL are the energies of the lowest-lying HH and LH states, respectively, and are strong functions of the gate electric field. These formulas explain three main features.
Comparison of the qubit Zeeman splitting between Schrieffer-Wolff transformation (to the third order) and exact numerical diagonalization for four different configurations. When the gate electric field is turned off, the qubit Zeeman splitting g0μBB ≈ 110 μeV. In all these figures, the out-of-plane magnetic field is B = 0.1 T. We can notice that the sweet spot does not change much as a function of the quantum dot radius, but the size of the qubit Zeeman splitting will be smaller for a larger quantum dot size. In all of these plots, we have \(\hbar\)ωl ≫ g0μBB. Numerical diagonalization is the red curve, Schrieffer-Wolff method is the blue curve. a d = 11 nm, a0 = 50 nm. b d = 11 nm, a0 = 60 nm. c d = 15 nm, a0 = 50 nm. d d = 15 nm, a0 = 60 nm.
Firstly, the optimal operation point reflects the interplay of the quadrupole degree of freedom with the gate electric field unique to spin-3/2 systems. The behavior of the qubit Zeeman splitting and Rashba coefficients is understood by recalling that the Rashba effect for the HH sub-bands is primarily driven by the off-diagonal matrix element L in Eq.(1) connecting the HH and LH sub-bands. This term, which is ∝kzk+, increases with the top gate field. At small gate fields, the Rashba spin-orbit constants increase monotonically due to the increase in the kz overlap integral. This continues until a critical top gate field is reached at which the HH–LH splittings, determined by the matrix element Q, begin to increase faster than the off-diagonal matrix element L. The heavy hole-light hole splitting induced by the confinement potential and the gate electric field is traced to the different effective masses for heavy and light holes. This physics has been shown previously by Winkler and collaborators78,79,80. Beyond this critical field, the Rashba terms decrease, resulting in a relatively broad optimal operation region at which the qubit is insensitive to background electric field fluctuations in the \(\hat{z}\)-direction and the dephasing rate vanishes to first order in the \(\hat{z}\)-electric field. As we show below, electric field fluctuations in the \(\hat{z}\)-direction are by far the most damaging to the qubit and are the key source of decoherence to be avoided. The breadth and smoothness of the extreme make the tuning of the electric field to reach the optimal operation point easier.
Secondly, the sweet spot shifts slightly with the dot radius64, and this is fully captured by our Schrieffer-Wolff results. The reason for this is that, in a 2D system, while the angular form of the Rashba interaction is dictated by rotational symmetry, the Rashba coupling constants α2 and α3 are functions of the magnitude of the wave vector. This implies they are functions of the density, and this is vital at large densities. Hence a quantum dot can be envisaged as having Rashba parameters that are functions of the in-plane radius a0, and their dependence is less pronounced at larger a0, since this corresponds to smaller densities.
Thirdly, each of the two spin-orbit coupling terms in Eq. (5) can be envisaged as the interaction of the hole spin with an effective spin-orbit field that depends on the momentum. In the absence of a magnetic field, the α2-Rashba spin-orbit field winds around the Fermi surface three times, whereas the α3-Rashba spin-orbit field winds only once. In the strict 2D limit it is the α3-term that enables EDSR. Although the quantum dot is not in the exact 2D limit, it still holds that EDSR is enabled by the cubic symmetry terms ∝δ. Setting δ = 0 in our calculations causes the EDSR frequency to vanish.
where u are a long expression that can be found in the Supplementary Information Eq. (17).
Dephasing time
The main dephasing mechanisms are fluctuating electrical fields such as charge noise. We focus on random telegraph noise (RTN) due to charge defects, noting that a similar discussion can be presented for 1/f noise, which is typically caused by an incoherent superposition of RTN sources. For this reason, we expect the trends for the two types of noise to be similar, while reliable numbers for 1/f noise must await the experimental determination of the noise spectral density S(ω) for hole qubits. To begin with, we estimate the dephasing time \({T}_{2}^{* }\), which is expected to be primarily determined by fluctuations in the Larmor frequency of the qubit induced by charge noise. The electric potential induced at the qubit by a defect located at rD, which may give rise to RTN, can be modeled as a quasi-2D screened Coulomb potential:
where ϵ0 is the vacuum permeability, ϵr is the relative permeability for Ge, qTF is the Thomas-Fermi wave vector, and kF is the Fermi wave vector49. In a dilution refrigerator, the high energy modes of the Coulomb potential are negligible, therefore the q > 2kF part is ignored. Another source of dephasing is dipole defects due to the asymmetry in bond polarities.
where RD is the distance between the dot and the unscreened charge dipole. p is the dipole moment of the charge p = el, the size of the dipole is about 1Å.
As a worst-case estimate of the dephasing time, we use the motional narrowing result49,50, the dephasing time \({T}_{2}^{* -1}={(\delta \omega )}^{2}\tau /2\), where δω is the change in qubit Larmor frequency due to the fluctuator, and we consider τ = 103tRabi, where tRabi is the single-qubit operation time (the inverse of the EDSR frequency), which can be found from Fig. 3a. Because of the weak coupling between the spin degree of freedom and external reservoirs, slower fluctuators can be eliminated via pulse sequences and the spin-echo techniques81. We consider two sample defects separately. One is a single-charge defect located 100 nm away from the quantum dot in the plane of the dot as a worst-case scenario for a charge trap. We use rD = 100 nm since regions inside this range will be depleted by the top gate, and charge traps will not be active. We also consider a dipole defect immediately under the gate and above the dot, with RD = 20 nm in the \(\hat{z}\) direction. This is because within the depleted region the most relevant defects are charge dipoles, whose orientation fluctuates. To estimate the pure dephasing time at the optimal operation point due to such a defect, we first note that the in-plane electric field will not contribute to dephasing. An in-plane electric field enters the QD Hamiltonian as E∥ ⋅ r∥. This in-plane electric field term does not couple states with different spin orientations. When we consider the qubit Zeeman splittings, the corrections to the effective quantum dot levels due to the in-plane electric field will read the same for H1,1 and H2,2 up to the second-order, therefore, fluctuations in qubit Zeeman splitting H1,1 − H2,2 will not depend on the in-plane electric field. However, higher-order terms in the expansion of the electrostatic potential of the defects will lead to dephasing, and these are responsible for dephasing at the optimal operation point itself. To determine their effect, we write the ground state energy as ELK + E0 + Ez + v0 where E0 is the lateral confinement energy, Ez is the Zeeman energy, and v0 is the energy correction due to the defect.
In all plots, the quantum well width is d = 11 nm and dot radius a0 = 50 nm. a Dephasing time in the motional narrowing regime. b The allowable number of single-qubit operations in one dephasing time in motional narrowing regime. c Dephasing time in the quasi-static limit. d The allowable number of single-qubit operations in one dephasing time in the quasi-static limit.
We estimate the approximate qubit window of operation around the optimal operation point. Away from the optimal operation point, due to the fluctuating electric potential of the defect, the energy levels of the quantum dot will gain a correction, i.e., \(\left\langle {\phi }_{{n}_{1},{m}_{1}}\right|{U}_{\text{sc}}\left|{\phi }_{{n}_{2},{m}_{2}}\right\rangle\). With these assumptions, the dephasing time is plotted as a function of the gate electric field in Fig. 3. At the optimal operation point, the dephasing time due to the out-of-plane fluctuations is calculated to the second-order, since the first-order fluctuation vanishes, the in-plane fluctuations will dominate the dephasing. Away from the optimal operation point, the motional narrowing result is much smaller than the quasi-static limit result. This is because the first-order variation of the qubit Zeeman splitting will weaken the correlation time, while the quasi-static limit does not consider any correlations. However, as the gate electric field approaches the optimal operation point, the variation of qubit Zeeman splitting decreases; at the optimal operation point, compared with the quasi-static limit result, a longer correlation time will lead to a larger dephasing time. We also determine the pure dephasing time in the quasi-static limit, where the switching time is the longest time scale in the system. This is essentially given by T2 = 2π/(δω), and is plotted in Fig. 3b.
Relaxation time and EDSR
We briefly discuss electrically driven spin resonance. An in-plane oscillating electric field represented in the Hamiltonian by eEAC(t)x drives spin-conserving transitions between the QD states. For a multiple occupied hole dot, the excited state structure may be more complex but the argument above remains valid because the α2 and α3 Rashba terms couple the ground state to different excited states. The spin resonance Rabi time is the time taken to accomplish an operation. The Rabi frequency can be tuned by changing the gate electric field and with it the Rashba spin-orbit coupling constant. However, note that because the two Rashba terms directly determine the correction to the g-factor, the Rashba interaction and the g-factor cannot be tuned independently at present.
Given that the spin resonance frequency is a maximum at the optimal operation point, it follows that the qubit can be tuned to have maximum coherence and maximum electrical driving simultaneously. The nonlinearity in the hole Rashba interaction as a function of the gate field that enables this feature has no counterpart in electron systems. In GaAs hole systems27 this nonlinearity does not lead to optimal operation points. This is because, firstly, the spin resonance Rabi frequency in GaAs hole qubits is driven by the Dresselhaus interaction, which is not tunable via the gate electric field, while α3 is negligible in GaAs. Secondly, GaAs qubits are exposed to decoherence through the hyperfine interaction, piezoelectric phonons, and the Dresselhaus interaction, none of which can be mitigated.
Since the Rabi frequency is maximized at the optimal operation point, the relaxation time T1 is minimized there. For the qubit to be operated efficiently it is vital to determine the ratio of the EDSR and relaxation rates. Hyperfine interactions and phonon-hole interactions are two major factors affecting the relaxation time, hence the quality of the qubit. However, the p-type symmetry of the valence band excludes the contact hyperfine interaction. There is no bulk inversion asymmetry in group IV elements; this leads to no Dresselhaus spin-orbit coupling. However, there is still the Rashba spin-orbit coupling due to the structure inversion asymmetry, which couples the heavy-hole states to the light-hole states. Neither the spin nor the orbital angular momentum will be a good quantum number, as the admixture of the spin-down and the spin-up states will modify the wave functions. We emphasize that, whereas EDSR comes only from the α3-Rashba term, the qubit relaxation is caused by both the α2- and the α3-Rashba terms.
The relaxation time evaluated using Fermi’s golden rule is shown in Fig. 4. For completeness, we also consider two-phonon relaxation processes, which include virtual emission and absorption of a phonon between two heavy hole states, since in the first-order relaxation calculation there is no direct matrix element between the two heavy-hole states. However, the two-phonon process calculation returns a negligible relaxation rate, which will not contribute significantly to the relaxation time. The relaxation rate will depend on the external magnetic field as (1/T1) ∝ B7 for the α3-Rashba term and (1/T1) ∝ B9 for the α2-Rashba term. This is shown in Fig. 4a.
In all plots d = 11 nm, a0 = 50 nm, the external magnetic field is B = 0.1 T. The density of Ge ρ = 5.33 × 103 kg/m3. The phonon propagation speed along the transverse direction is vt = 3.57 × 103 m/s, along the longitudinal direction it is vl = 4.85 × 103 m/s. a Relaxation time and EDSR Rabi time as a function of the gate electric field. b The allowable number of single-qubit operations in one relaxation time.
We also plot the ratio between the relaxation time and the EDSR time, demonstrating that the system allows for a large number of operations. The allowable number of single-qubit operations is calculated by evaluating the ratio of the relaxation time and the EDSR time, i.e., the Rabi ratio. The in-plane electric field we used is EAC = 103 V/m. In Fig. 4, we plot the relaxation time, EDSR Rabi time, comparison of the magnitude of the relaxation time and EDSR time and an estimation of red the allowable number of single-qubit operations as the function of the gate electric field at a magnetic field B = 0.1 T which is parallel to the growth direction. The relaxation time calculations mainly consider the hole-phonon interactions and the details can be found in the Supplementary Information Eq. (22). Both the relaxation time and the EDSR time will depend on the spin-orbit coupling coefficients, therefore, their extrema coincide. The relaxation time increases for smaller dot sizes because that corresponds to larger confinement energy, while both the phonon and the spin-orbit coupling terms connect the orbital ground state to higher excited states. From Fig. 4a, we can see that the Ge hole quantum dot has a long relaxation time and large Rabi ratio at dilution refrigerator temperatures. It is also useful to study the relaxation time at slightly higher temperatures, e.g., 4 K, at which both phonon absorption and emission must be taken into account. The phonon occupation number is given by the Bose-Einstein distribution \(N={({e}^{\hslash \omega /({k}_{\text{B}}T)}-1)}^{-1}\), where N is the occupation number, ω = qv, q is the phonon wave vector and v is the phonon propagation velocity, T is the temperature, kB is the Boltzmann constant. More details can be found in the Supplementary Information Eq. (22), where a plot of the temperature dependence of the relaxation rate is presented as well. For T = 4 K, the relaxation time is 17 ms, suggesting that the qubit can easily be operated at this temperature.
Applicability and implementation
Although we have used a simple parabolic model for the in-plane QD confinement, our conclusions are very general. Firstly, the dephasing optimal operation point will be present for potentials of arbitrary complexity (for example hut wire geometries)7,55,56, since it is due to the fundamental interplay between the HH and LH that gives rise to the Rashba spin-orbit coupling in the HH manifold. Secondly, we have examined the possibility that the insensitivity of the g-factor to in-plane electric fields is an artifact of the model. We have tested three deviations from parabolicity and found that none of them exposes the qubit to dephasing by fluctuating in-plane electric fields. This implies (i) that the dot does not have to be perfectly parabolic allowing for some flexibility in the gate structure; (ii) that in-plane electric field fluctuations generally have a negligible effect on the g-factor, while out-of-plane electric field fluctuations cause fluctuations in the Rashba spin-orbit coupling and affect the g-factor, therefore it is most important to avoid the effect of the out-of-plane field; and (iii) that dephasing at the optimal operation point itself comes about primarily from higher-order terms in the electrical potential, i.e., electrical quadrupole and higher. Our results hold qualitatively in Si as well, where the spin-orbit interaction is weaker than in Ge, while δ is larger. However, the large δ and frequent failure of the Schrieffer-Wolff approximation in Si calls for fully numerical treatment48.
Experimentally, the configuration we describe requires a double-gated device with separate plunger gates and barrier gates allowing the number density and the gate electric-field (and spin-orbit coupling) to be controlled independently58. The numerical estimates above suggest that, in general, a smooth and broad optimal operation point will enable the Ge hole qubit to work insensitively to the charge noise inside a large range of gate electric fields accessible to experiment. Exchange-based two-qubit gates should be possible for hole QDs, and their speed depends on the values of exchange obtained, which are expected to be tunable by gates. Moreover, it is likely to simplify the coupling between the two qubits since the valley degree of freedom is absent in hole systems. However, a two-qubit gate in the setup discussed here is not optimized for long-distance coupling, which leads to the two-qubit gate time is of the order of microseconds for dipole-dipole interactions and hundreds of microseconds for circuit QED, limited by the Ge Luttinger parameters. They can be sped up by enhancing the spin-orbit interaction, but we defer the discussion to a future publication.
A smaller g-factor will lead to a smaller Rabi frequency, a smaller change in the qubit Zeeman splitting due to the spin-orbit interaction, and a shorter dephasing time but a longer relaxation time and an improved Rabi ratio. The optimal operation point will not change its location, which is determined only by the effective mass and the width of the quantum well. A larger quantum dot radius would make the confinement energy smaller, increasing the effect of the spin-orbit interaction and resulting in a faster Rabi frequency, but also shorter T1 and T2*. Nevertheless, the Rabi ratio decreases with increasing dot radius. Moreover, since the confinement energy decreases as the square of the radius, it is preferable to work at smaller radii to ensure the thermal broadening is overcome. Increasing x increases the heavy hole - light hole splitting, leading to a reduced Rashba spin-orbit coupling and a smaller change in the qubit Zeeman splitting. The change in the Zeeman splitting will be large for smaller x e.g., 0.05 − 0.10, while at x = 0.3 it is essentially not noticeable.
We have demonstrated that electrostatically defined hole quantum dot spin qubits naturally exhibit an optimal operation point at which sensitivity to charge noise is minimized while the speed of electrical operation is maximized. The location of the optimal operation point can be determined from the width of the quantum well and the strain tensors applied. Relaxation times are long even at 4 K, while dephasing is determined by higher-order terms in the expansion of the electrostatic potential due to charge defects, but are expected to allow for a large window of operation around the optimal operation point. Our results provide a theoretical guideline for achieving fast, highly coherent, low-power electrically operated spin qubits experimentally. Future studies must consider in-plane magnetic fields, which interact much more weakly with HH spins and are more complicated to treat theoretically.
Methods
Numerical diagonalization
Most of the results are obtained by theoretical analysis, we also present a numerical simulation by diagonalizing the full Hamiltonian. The full Hamiltonian is projected onto all the heavy-hole states and light-hole states and all the Fock-Darwin state.
Data availability
Data sharing not applicable to this article as no data sets were generated or analyzed during this study.
References
Loss, D. & DiVincenzo, D. P. Quantum computation with quantum dots. Phys. Rev. A 57, 120–126 (1998).
Fischer, J., Trif, M., Coish, W. & Loss, D. Spin interactions, relaxation and decoherence in quantum dots. Solid State Commun. 149, 1443 – 1450 (2009).
Chesi, S., Wang, X. J. & Coish, W. A. Controlling hole spins in quantum dots and wells. Eur. Phys. J. Plus 129, 86 (2014).
Lodari, M. et al. Light effective hole mass in undoped Ge/SiGe quantum wells. Phys. Rev. B 100, 041304 (2019).
Liles, S. D. et al. Spin and orbital structure of the first six holes in a silicon metal-oxide-semiconductor quantum dot. Nat. Commun. 9, 3255 (2018).
Watzinger, H. et al. A germanium hole spin qubit. Nat. Commun. 9, 3902 (2018).
Watzinger, H. et al. Heavy-hole states in germanium hut wires. Nano Lett. 16, 6879–6885 (2016).
Vukušić, L. et al. Single-shot readout of hole spins in Ge. Nano Lett. 18, 7141–7145 (2018).
Ares, N. et al. Sige quantum dots for fast hole spin Rabi oscillations. Appl. Phys. Lett. 103, 263113 (2013).
Mizokuchi, R., Maurand, R., Vigneau, F., Myronov, M. & De Franceschi, S. Ballistic one-dimensional holes with strong g-factor anisotropy in germanium. Nano Lett. 18, 4861–4865 (2018).
Bohuslavskyi, H. et al. Pauli blockade in a few-hole PMOS double quantum dot limited by spin-orbit interaction. Appl. Phys. Lett. 109, 193101 (2016).
Zwanenburg, F. A., van Rijmenam, C. E. W. M., Fang, Y., Lieber, C. M. & Kouwenhoven, L. P. Spin states of the first four holes in a silicon nanowire quantum dot. Nano Lett. 9, 1071–1079 (2009).
Spruijtenburg, P. C. et al. Single-hole tunneling through a two-dimensional hole gas in intrinsic silicon. Appl. Phys. Lett. 102, 192105 (2013).
Brauns, M. et al. Anisotropic Pauli spin blockade in hole quantum dots. Phys. Rev. B 94, 041411 (2016).
Salfi, J., Mol, J. A., Culcer, D. & Rogge, S. Charge-insensitive single-atom spin-orbit qubit in silicon. Phys. Rev. Lett. 116, 246801 (2016).
Salfi, J., Tong, M., Rogge, S. & Culcer, D. Quantum computing with acceptor spins in silicon. Nanotechnology 27, 244001 (2016).
Hung, J.-T., Marcellina, E., Wang, B., Hamilton, A. R. & Culcer, D. Spin blockade in hole quantum dots: Tuning exchange electrically and probing Zeeman interactions. Phys. Rev. B 95, 195316 (2017).
Terrazos, L. A. et al. Qubits Based on Hole Quantum Dots in Strained Ge. Phys. Rev. B Preprint at http://arxiv.org/abs/1803.10320. (2018).
van der Heijden, J. et al. Readout and control of the spin-orbit states of two coupled acceptor atoms in a silicon transistor. Sci. Adv. 4, eaat9199 (2018).
Chekhovich, E. A. et al. Nuclear spin effects in semiconductor quantum dots. Nat. Mater. 12, 494 EP – (2013).
Kyrychenko, F. V. & Kossut, J. Diluted magnetic semiconductor quantum dots: an extreme sensitivity of the hole Zeeman splitting on the aspect ratio of the confining potential. Phys. Rev. B 70, 205317 (2004).
Nadj-Perge, S., Frolov, S. M., Bakkers, E. P. A. M. & Kouwenhoven, L. P. Spin-orbit qubit in a semiconductor nanowire. Nature 468, 1084 EP (2010).
Golovach, V. N., Khaetskii, A. & Loss, D. Spin relaxation at the singlet-triplet crossing in a quantum dot. Phys. Rev. B 77, 045328 (2008).
Rossi, A. et al. Detection of variable tunneling rates in silicon quantum dots. Appl. Phys. Lett. 98, 133506 (2011).
Yoneda, J. et al. A quantum-dot spin qubit with coherence limited by charge noise and fidelity higher than 99.9%. Nat. Nanotechnol. 13, 102–106 (2018).
Stano, P. et al. Orbital effects of a strong in-plane magnetic field on a gate-defined quantum dot. Phys. Rev. B 99, 085308 (2019).
Bulaev, D. V. & Loss, D. Electric dipole spin resonance for heavy holes in quantum dots. Phys. Rev. Lett. 98, 097202 (2007).
Golovach, V. N., Borhani, M. & Loss, D. Electric-dipole-induced spin resonance in quantum dots. Phys. Rev. B 74, 165319 (2006).
Coish, W. A., Golovach, V. N., Egues, J. C. & Loss, D. Measurement, control, and decay of quantum-dot spins. Phys. Status Solid. 243, 3658–3672 (2006).
Szumniak, P., Bednarek, S., Pawłowski, J. & Partoens, B. All-electrical control of quantum gates for single heavy-hole spin qubits. Phys. Rev. B 87, 195307 (2013).
Budich, J. C., Rothe, D. G., Hankiewicz, E. M. & Trauzettel, B. All-electric qubit control in heavy hole quantum dots via non-abelian geometric phases. Phys. Rev. B 85, 205425 (2012).
Salis, G. et al. Electrical control of spin coherence in semiconductor nanostructures. Nature 414, 619–622 (2001).
Pribiag, V. S. et al. Electrical control of single hole spins in nanowire quantum dots. Nat. Nanotechnol. 8, 170 EP – (2013).
Khaetskii, A. V. Spin relaxation in semiconductor mesoscopic systems. Phys. E10, 27 – 31 (2001).
Gvozdić, D. M. & Ekenberg, U. Superiority of p-type spin transistors. Phys. Scripta T126, 21–26 (2006).
Woods, L. M., Reinecke, T. L. & Kotlyar, R. Hole spin relaxation in quantum dots. Phys. Rev. B 69, 125330 (2004).
Hsieh, C.-Y., Cheriton, R., Korkusinski, M. & Hawrylak, P. Valence holes as Luttinger spinor based qubits in quantum dots. Phys. Rev. B 80, 235320 (2009).
Kloeffel, C., Trif, M. & Loss, D. Strong spin-orbit interaction and helical hole states in Ge/Si nanowires. Phys. Rev. B 84, 195314 (2011).
Machnikowski, P., Gawarecki, K. & Cywiński, L. Hyperfine interaction for holes in quantum dots: k ⋅ p model. Phys. Rev. B 100, 085305 (2019).
Szumniak, P., Bednarek, S., Partoens, B. & Peeters, F. M. Spin-orbit-mediated manipulation of heavy-hole spin qubits in gated semiconductor nanodevices. Phys. Rev. Lett. 109, 107201 (2012).
Khaetskii, A. V. & Nazarov, Y. V. Spin relaxation in semiconductor quantum dots. Phys. Rev. B 61, 12639–12642 (2000).
Khaetskii, A. V., Loss, D. & Glazman, L. Electron spin decoherence in quantum dots due to interaction with nuclei. Phys. Rev. Lett. 88, 186802 (2002).
Fischer, J. & Loss, D. Hybridization and spin decoherence in heavy-hole quantum dots. Phys. Rev. Lett. 105, 266603 (2010).
Kobayashi, T. et al. Engineering long spin coherence times of spin–orbit qubits in silicon. Nat. Mater. 20, 38–42, https://doi.org/10.1038/s41563-020-0743-3 (2021).
Yamahata, G. et al. Magnetic field dependence of Pauli spin blockade: a window into the sources of spin relaxation in silicon quantum dots. Phys. Rev. B 86, 115322 (2012).
Hendrickx, N. W., Franke, D. P., Sammak, A., Scappucci, G. & Veldhorst, M. Fast two-qubit logic with holes in germanium. Nature 577, 487 (2020).
Durnev, M. V., Glazov, M. M. & Ivchenko, E. L. Spin-orbit splitting of valence subbands in semiconductor nanostructures. Phys. Rev. B 89, 075430 (2014).
Marcellina, E., Hamilton, A. R., Winkler, R. & Culcer, D. Spin-orbit interactions in inversion-asymmetric two-dimensional hole systems: A variational analysis. Phys. Rev. B 95, 075305 (2017).
Bermeister, A., Keith, D. & Culcer, D. Charge noise, spin-orbit coupling, and dephasing of single-spin qubits. Appl. Phys. Lett. 105, 192102 (2014).
Culcer, D., Hu, X. & Das Sarma, S. Dephasing of si spin qubits due to charge noise. Appl. Phys. Lett. 95, 073102 (2009).
Fischer, J., Coish, W. A., Bulaev, D. V. & Loss, D. Spin decoherence of a heavy hole coupled to nuclear spins in a quantum dot. Phys. Rev. B 78, 155329 (2008).
Csontos, D., Brusheim, P., Zülicke, U. & Xu, H. Q. Spin-\(\frac{3}{2}\) physics of semiconductor hole nanowires: valence-band mixing and tunable interplay between bulk-material and orbital bound-state spin splittings. Phys. Rev. B 79, 155323 (2009).
Meier, L. et al. Measurement of Rashba and Dresselhaus spin–orbit magnetic fields. Nat. Phys. 3, 650 EP – (2007).
Governale, M. & Zülicke, U. Spin accumulation in quantum wires with strong Rashba spin-orbit coupling. Phys. Rev. B 66, 073311 (2002).
Li, S.-X. et al. Measuring hole spin states of single quantum dot in germanium hut wire. Appl. Phys. Lett. 110, 133105 (2017).
Li, Y. et al. Coupling a germanium hut wire hole quantum dot to a superconducting microwave resonator. Nano Lett. 18, 2091–2097 (2018).
Kloeffel, C., Trif, M., Stano, P. & Loss, D. Circuit QED with hole-spin qubits in Ge/Si nanowire quantum dots. Phys. Rev. B 88, 241405 (2013).
Sammak, A. et al. Shallow and undoped germanium quantum wells: a playground for spin and hybrid quantum technology. Adv. Funct. Mater. 29, 1807613 (2019).
Hardy, W. J. et al. Single and double hole quantum dots in strained Ge/SiGe quantum wells. Nanotechnology 30, 215202 (2019).
Takeda, K. et al. Characterization and suppression of low-frequency noise in Si/SiGe quantum point contacts and quantum dots. Appl. Phys. Lett. 102, 123113 (2013).
Katsaros, G. et al. Observation of spin-selective tunneling in SiGe nanocrystals. Phys. Rev. Lett. 107, 246601 (2011).
Rol, M. et al. Fast, high-fidelity conditional-phase gate exploiting leakage interference in weakly anharmonic superconducting qubits. Phys. Rev. Lett. 123, https://doi.org/10.1103/PhysRevLett.123.120502. (2019).
Luttinger, J. M. & Kohn, W. Motion of electrons and holes in perturbed periodic fields. Phys. Rev. 97, 869–883 (1955).
Marcellina, E. et al. Electrical control of the Zeeman spin splitting in two-dimensional hole systems. Phys. Rev. Lett. 121, 077701 (2018).
Liu, H., Marcellina, E., Hamilton, A. R. & Culcer, D. Strong spin-orbit contribution to the hall coefficient of two-dimensional hole systems. Phys. Rev. Lett. 121, 087701 (2018).
Bastard, G., Mendez, E. E., Chang, L. L. & Esaki, L. Variational calculations on a quantum well in an electric field. Phys. Rev. B 28, 3241–3245 (1983).
Golovach, V. N., Khaetskii, A. & Loss, D. Phonon-induced decay of the electron spin in quantum dots. Phys. Rev. Lett. 93, 016601 (2004).
Trif, M., Simon, P. & Loss, D. Relaxation of hole spins in quantum dots via two-phonon processes. Phys. Rev. Lett. 103, 106601 (2009).
Bulaev, D. V. & Loss, D. Spin relaxation and decoherence of holes in quantum dots. Phys. Rev. Lett. 95, 076805 (2005).
Maier, F., Kloeffel, C. & Loss, D. Tunable g factor and phonon-mediated hole spin relaxation in Ge/Si nanowire quantum dots. Phys. Rev. B 87, 161305 (2013).
Woods, L. M., Reinecke, T. L. & Lyanda-Geller, Y. Spin relaxation in quantum dots. Phys. Rev. B 66, 161318 (2002).
Climente, J. I., Segarra, C. & Planelles, J. Spin–orbit-induced hole spin relaxation in InAs and GaAs quantum dots. New Journal of Physics 15, 093009 (2013).
Kornich, V., Kloeffel, C. & Loss, D. Phonon-mediated decay of singlet-triplet qubits in double quantum dots. Phys. Rev. B 89, 085410 (2014).
Ares, N. et al. Nature of tunable hole g factors in quantum dots. Phys. Rev. Lett. 110, 046602 (2013).
Laucht, A. et al. Electrically controlling single-spin qubits in a continuous microwave field. Sci. Adv. 1, e1500022 (2015).
Li, R., Hudson, F. E., Dzurak, A. S. & Hamilton, A. R. Pauli spin blockade of heavy holes in a silicon double quantum dot. Nano Lett. 15, 7314–7318 (2015).
Bravyi, S., DiVincenzo, D. P. & Loss, D. Schrieffer-Wolff transformation for quantum many-body systems. Ann. Phys. 326, 2793 – 2826 (2011).
Winkler, R. Rashba spin splitting in two-dimensional electron and hole systems. Phys. Rev. B 62, 4245–4248 (2000).
Winkler, R., Papadakis, S. J., De Poortere, E. P. & Shayegan, M. Highly anisotropic g-factor of two-dimensional hole systems. Phys. Rev. Lett. 85, 4574–4577 (2000).
Habib, B. et al. Negative differential Rashba effect in two-dimensional hole systems. Appl. Phys. Lett. 85, 3151–3153 (2004).
Press, D. et al. Ultrafast optical spin echo in a single quantum dot. Nat. Photon. 4, 367 (2010).
Acknowledgements
We thank Mark Friesen and András Pályi for a multitude of engaging and educational discussions. We are also grateful to Thaddeus Ladd, Sue Coppersmith, and Andrè Saraiva for stimulating feedback. This research is supported by the Australian Research Council Centre of Excellence in Future Low-Energy Electronics Technologies (project CE170100039) and funded by the Australian Government.
Author information
Authors and Affiliations
Contributions
Z.W. performed all the calculations presented in this work in collaboration with E.M. J.H.C. assisted with the theoretical calculations. J.S., A.H., and S.R. contributed to the discussion of experimental implementation. D.C. and J.S. devised the theoretical model, supervised the calculations, and wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, Z., Marcellina, E., Hamilton, A.R. et al. Optimal operation points for ultrafast, highly coherent Ge hole spin-orbit qubits. npj Quantum Inf 7, 54 (2021). https://doi.org/10.1038/s41534-021-00386-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41534-021-00386-2