
Recent advances in hole-spin qubits 5
usual envelope-function approximation. Such pseudospin-electric coupling physically
originates from the interaction between the electric field and the dipole operator, which
takes the general form of Eq. (5) when expressed in the subspace of valence-band states.
A first-principles evaluation gives a relatively large value χ'0.2in GaAs [39].
For Si and Ge, the above Dresselhaus terms are absent in the bulk. However, it is
important to account for finite strain when modelling, e.g., high-mobility SiGe/Ge/SiGe
quantum wells [40, 41] as well as core-shell nanowires [42, 43]. In general, the effect of
strain is described by the well-known Bir-Pikus Hamiltonian Hε=−(av+5
4bv)(εxx +
εyy +εzz ) + bv(J2
xεxx + c.p.) + 2
√3dv({Jx, Jy}εxy + c.p.), where {εij }are the strain tensor
elements and av, bv, dvare material-specific deformation potentials [44]. Assuming that
only hydrostatic and uniaxial strain are present, i.e., εxy =εyz =εzx = 0 and εxx =εyy,
the effect of Hεis to induce different offsets for the heavy-hole (Jz=±3/2) and light-hole
(Jz=±1/2) states [40–42].
1.1. Spin-orbit interaction of confined holes
The description of confined hole systems can be based directly on the general four-band
formalism just described, eventually supplemented by the split-off band. However,
many of the hole-spin qubits of interest are fabricated from low-dimensional systems
such as quantum wells or nanowires, where additional gates realize the zero-dimensional
confinement. It is therefore useful to describe the effective Hamiltonians governing such
2D or 1D systems. Because such effective models are relatively simple, they are often
taken as starting points to characterize various important effects, such as spin relaxation
or electrical manipulation of spin qubits.
Quantum wells. For quasi-2D holes the potential in Eq. (2) may be taken as V(z)
if the growth direction is along a main crystal axis. Then, for B= 0 and zero in-
plane wavevector, kk= (kx, ky) = (0,0), Eq. (2) simplifies considerably, giving two
decoupled Schrödinger equations in 1D. The two sectors correspond to the Jz=±3/2
and Jz=±1/2subspaces, characterized by different effective masses m0/(γ1±2γ2).
The lowest-energy subband is formed by heavy-holes (Jz=±3/2) and the energy gap
to the lowest light-hole subband can be estimated as ∆LH '2γ2(π~/L)2/m0, assuming
an infinite potential well of width L. For more realistic scenarios, the value of ∆LH
depends of the detailed shape of V(z)and, as discussed already, is modified by uniaxial
strain. In general, the four-dimensional degeneracy of the Luttinger Hamiltonian has
been broken by the confining potential, generating naturally a two-level pseudospin. For
the low-energy states, we can identify the effective spin up/down eigenstates with the
Jz=±3/2eigenvalues, i.e., |⇑i =|+3/2iand |⇓i =|−3/2i.
Despite the simple decoupling described above (which, we stress, is exact only at
zero in-plane momentum), the light-hole/heavy-hole mixing terms still play a crucial role
in the effective 2D Hamiltonian. As seen from Eq. (1), finite values of kx,y will induce
matrix elements between the two subspaces, through the action of the Jx,y operators.
It is this coupling between light- and heavy-hole subbands what generates the desired