Spin-orbit interaction induced anisotropic property in interacting quantum wires

We investigate theoretically the ground state and transport property of electrons in interacting quantum wires (QWs) oriented along different crystallographic directions in (001) and (110) planes in the presence of the Rashba spin-orbit interaction (RSOI) and Dresselhaus SOI (DSOI). The electron ground state can cross over different phases, e.g., spin density wave, charge density wave, singlet superconductivity, and metamagnetism, by changing the strengths of the SOIs and the crystallographic orientation of the QW. The interplay between the SOIs and Coulomb interaction leads to the anisotropic dc transport property of QW which provides us a possible way to detect the strengths of the RSOI and DSOI. PACS numbers: 73.63.Nm, 71.10.Pm, 73.23.-b, 71.70.Ej


Introduction
All-electrical manipulation of spin degree of freedom is one of the central issues and the ultimate goal of spintronics field. The spin-orbit interaction (SOI) is a manifestation of special relativity. An electric field in the laboratory frame can transform into an effective magnetic field in the moving frame of electron and consequently leads to electron spin splitting. Therefore, the SOI provides us an efficient way to control electron spin electrically and has attracted tremendous interest because of its potential application in all-electrical spintronic devices [1,2]. The spin degeneracy can be lifted by applying magnetic field to break the time-reversal symmetry and/or by applying electric fields to break the spatial inversion symmetry. However, the latter could be more easily realized in spintronics devices. In semiconductors the spatial inversion symmetry can be broken by the structural inversion asymmetry and bulk crystal inversion asymmetry, named, Rashba SOI (RSOI) and Dresselhaus SOI (DSOI), respectively [3,4]. Usually, the RSOI in semiconductor quantum well is much stronger than that of DSOI, and therefore, most of the previous theoretical and experimental studies focused on the RSOI and its consequence on the spin transport properties in twodimensional electron gas (2DEG) [5][6][7]. In thin quantum wells, the strength of the DSOI is comparable to that of the RSOI since the strength of the DSOI depends significantly on the thickness of quantum wells. The interplay between the RSOI and DSOI leads to interesting phenomena in 2DEG, e.g., the anisotropic photogalvanic effect [8], and the persistent spin helix [9]. However, the interplay between the RSOI and DSOI in quantum wires (QWs) remains relatively unexplored.
Very recently, the anisotropic behavior of transport property in semiconductor QWs was proposed to detect the relative strength between the RSOI and DSOI in a quasi-one-dimensional (Q1D) semiconductor QW system [10,11]. However, the effect of the Coulomb interaction on the transport property is not addressed. Since the Coulomb interaction becomes very strong in Q1 D electron systems where electrons are strongly correlated, and therefore the conventional Fermi liquid theory breaks down. There are no fermionic quasi-particle in Q1 D electron gas, and the elementary excitations are bosonic collective charge and spin fluctuations with different propagating velocities. The Luttinger liquid (LL) theory [12] is of fundamental importance because it is one of a very few strongly correlated non-Fermi liquid systems that can be solved analytically. The LL displays very unique properties, e.g., the spin and charge separation and the power-law behavior of the correlation functions. The unique behavior was observed experimentally in many Q1 D systems, for instance, narrow QW formed in semiconductor heterostructures [13], carbon nanotube [14], graphene nanoribbon [15], as well as the edge states of the fractional Quantum Hall liquid [16]. Recent studies have found that the RSOI would lead to the mixing between the spin and charge excitations in QW with the RSOI alone [17][18][19][20]. It is interesting to study the interplay between the RSOI and DSOI on the ground state and transport property of QWs in the presence of the Coulomb interaction.
In this study, based on the LL theory, we study the effect of the interplay between the Coulomb interaction and SOIs on the electron ground state and transport property of an interacting QW oriented along different crystallographic directions in different planes. The electron ground state can display the transitions among the different phases, e.g., the spin density wave (SDW), charge density wave (CDW), singlet superconductivity (SS), and metamagnetism (MM), by tuning the crystallographic plane and orientation of the QW, the strengths of SOIs and the Coulomb interaction. The anisotropy of the dc conductivity of interacting QW is induced by the interplay between the Coulomb interaction and SOIs, which could be used for detecting the strengths of the RSOI and DSOI.

Theory
The Q1 D QW system is shown schematically in Figure 1a where electrons are confined laterally and move freely along the x-axis. The RSOI can be generated in the region below the top gate, and the DSOI always exists in the conventional zincblende semiconductors lacking the spatial inversion symmetry, e.g., GaAs, InAs, and InSb. The Hamiltonian of noninteracting electrons in a Q1 D QW oriented along different crystallographic directions in (001) plane is [11] where θ is the angle between the orientation of the QW and the [100] axis, m* is the electron effective mass, s i (i = x, y, z) are the Pauli matrices, and a and b are the strengths of the RSOI and DSOI, respectively. The linearized non-interacting electron Hamiltonian of the QW with both RSOI and DSOI is [17,20] where the operators ψ gs (g = -1(L),1(R); s = -1(↓),1(↑)) annihilate spin-down (↓) or spin-up (↑) electrons near the left (L) and right (R) Fermi points, respectively. v s γ = γ v F − sδv are the four different Fermi velocities, where v F is the bare Fermi velocity of right-and left-moving non-interacting electrons, δv = α 2 + β 2 + 2αβ sin 2θ /h , usually δv ≪ v F . Note that the RSOI and DSOI split the spin-up and spin-down subbands and make the electron Fermi velocities become orientation dependent. The total Hamiltonian of the system including the Coulomb interac- The Umklapp scattering process is neglected because the Fermi energy in QWs formed in semiconductor heterostructure is far from the half-filled case, and the electron-electron backscattering can be negligible for a sufficiently long interacting region [17]. Using the bosonization technique [21], the Hamiltonian becomes Where ϑ r and ϑ s are the phase fields for the charge and spin degrees of freedom, respectively, and Π r and Π σ are the corresponding conjugate momenta. v r and v s are the propagation velocities of the decoupled charge and spin-collective modes in the absence of the SOIs (δv = 0), respectively.
For a QW embedded in (110) plane, the Hamiltonian of the noninteracting electrons reads [11] where the parameters a, b, and θ are the same as in Equation (1).
induced by the DSOI is perpendicular to (110) plane, which could generate an out-of-plane spin-polarized current; and (3) the anisotropy is mainly determined by the DSOI. These differences would lead to the distinct anisotropic behavior of QWs embedded in (110) plane. For (111) plane, the DSOI shows the same formulism as the RSOI; the Hamiltonian is H 0 =h does not contain any θ-dependent term, and therefore the anisotropic behavior disappear. This isotropic character can be easily seen from the constant energy surface (see Figure 1b). We will not discuss the (111) plane any more.
The correlation functions at zero temperature in an interacting QW in the presence of the SOIs behave as R i (x) ∼ |x| −2+α i , and the exponents a i 's are determined by [18] α CDW = α SDW, z = 2 − K ρ ν ρ − K σ ν σ , Here Where δv 2 are the propagation velocities of coupled collective modes that depend on the crystallographic orientation θ. From the above equations, one can see that the interplay between the RSOI and DSOI will lead to anisotropy, i.e., the dependence on the crystallographic direction. This is the dominant difference between this work and the previous study [18]. Interestingly, for a QW in (001) plane, the spin and charge excitations can be decoupled again when θ = 3π/4 and a = b, θ = π/4 and a = -b, since the coupling between the spin and charge excitation disappear when δv = 0 (see the third term in Equation (4)).
Considering a point-like density-density interaction for the electron-electron interaction [22,23], we have v r,s = v F /K r,s and the parameter K r/s being defined as is the electron-electron interaction potential. We consider an infinitely long interacting QW, i.e., a homogenous interacting QW, driven by a time-dependent electric field E(x, t), which could be realized by applying a microwave radiation.
describes the effect of the ac electric field on the charge excitation. The total Hamiltonian becomes H = H 0 + H int + H ac . By minimizing the action functional of such a Q1 D system, we get the equation of motion and consequently obtain the non-local charge conductivity using the linear response theory: Numerical results and discussions First, we study how the interplay among the Coulomb interaction, the RSOI, and the DSOI affects the ground state of Q1 D electron gas. It is worth to note that if one does not assume any specific form for electronelectron interactions, then all four parameters v r (s) and K r (s) (or equivalently g 2|| , g 2⊥ , g 4|| and g 4⊥ ) are independent [18]. Figure 2 describes the phase diagram of Q1 D QW embedded in (001) plane as function of the strength of the Coulomb interaction K p and crystallographic orientations θ for the different strengths of SOIs a and b when K s = 0.7. We should stress that the different phases are determined by which correlation function decays most slowly for |x| ∞ when the other correlation functions become negligible. We define the phase by the dominant correlation function, e.g., CDW or SDW. Here the CDW and SDW are actually not the pure charge and spin fluctuations, but a mixed state of them induced by the SOIs. Interestingly, one can see that the ground state of the system displays a crossover between the CDW and SDW by tuning the crystallographic orientation for the fixed parameters, e.g., K r , a, and b. Note that the ground state of Q1 D electron system in the presence of RSOI alone cannot be affected by tuning the crystallographic orientation [18]. We find that there are crossovers of the ground state of Q1 D electron gas between CDW and SDW when we tune the strengths of the Coulomb interaction and the SOIs. For a fixed strength of the DSOI, the CDW phase region will be squeezed and die away gradually, with increase in the strengths of the RSOI. The ground state of the system is always SDW at strong RSOI for the arbitrary orientation of the QW. On the contrary, as the strength of the RSOI decreases for a given strength of the DSOI, the CDW phase region expands, and the ground state of the system becomes the CDW eventually. When δv >δv s , i.e., in the region below the dotted line, the ground state is the so-called MM phase, which was earlier observed in the Q1 D systems, e.g., Ba 3 Cu 2 O 4 Cl 2 [24], and attributed to the next-nearest-neighbor coupling in the XXZ model [25]. By tuning the parameters properly, the ground state of the system exhibits the crossovers among the SDW, CDW, SS, and MM.
The phase diagram of the ground state for the QW embedded in (110) plane is very different from that of the QW in (001) plane for the same parameters v r , v s and K s . Tuning the strength of the RSOI at a fixed strength of the DSOI, the ground state can also transit from CDW to SDW, but in very narrow region in the K r -θ space (see Figure 3). This is because the anisotropy is weaker compared to that of the QW in (001) plane.
Next, we turn to study the transport property of electron in an interacting QW. In the dc case, Equation (9) shows that the charge conductivity of an infinitely long interacting QW for the different parameters g and δv.
For the noninteracting QW structure, there is only single transverse mode, and thus no anisotropy of the dc conductivity can be found. This result can also be obtained from Equation (9) by taking g = 0, i.e., noninteracting case, u 2 1,2 = (v F ± δv) 2 ; we easily obtain s r = 2e 2 /h. When the Coulomb interaction is included but without the SOIs, i.e., a = b = 0, we can obtain an isotropic conductivity s r = 2K r e 2 /h from Equation (9). The spin and charge excitations propagate independently at the velocities v r and v s in this case. In the presence of both the Coulomb interaction and the SOIs, the SOIs mix the spin and charge excitations, leading to anisotropic velocities of the collective excitations u 1 and u 2 (see Equation 8), or equivalently the anisotropic interaction parameter K ρ , consequently resulting in the anisotropic conductivity σ ρ = 2K ρ e 2 /h. Therefore, the dc conductivity of the infinitely long-interacting QW depends sensitively on the crystallographic direction θ of the QW, i.e., the anisotropic transport behavior (see Figure 4). Interestingly, the dc conductivity oscillates with varying the angle θ with a periodicity π. With increasing the strengths of the DSOI and/or RSOI (see Figure 4a) or the strengths of the Coulomb interaction (see Figure  4b), the oscillation of the conductivity for the infinitely long-interacting QW becomes stronger. This is because the anisotropy of the interaction parameter K ρ depends on the anisotropic velocities of the mixed spin and charge excitations, which increases when the strengths of the SOIs or the Coulomb interaction increase.
For an interacting QWs embedded in (110) plane, the anisotropy of the dc conductivity becomes different (see Figure 4c,d): (1) The RSOI becomes less important for the anisotropy of the dc conductivity comparing with that in (001) plane, the increase of the strength of RSOI only affects the anisotropy slightly (see the dashed green line and dotted blue line in Figure 4c), since the anisotropy is mainly determined by the DSOI. This can be understood from the formulism of the Fermi velocities v s γ = γ v F − sδv , where δv = α 2 + β 2 sin 2 θ /4/h . It is noted that for (001) plane, the anisotropy of dc conductivity will disappear with the RSOI or DSOI alone, while for (110) plane, the anisotropy will disappear only when the DSOI is absent, which can be understood from the expression of δv. (2) The effect of the Coulomb interaction on the anisotropy of dc conductivity is weakened significantly when compared to that in (001) plane (see Figure 4b,d). Finally, we discuss how to detect the relative strength of the SOIs utilizing the anisotropic dc conductivity. We consider a QW embedded in (001) plane. One can tune the strength of RSOI by adjusting the gate voltage. For a QW oriented along the [110] axis, i.e., θ = 3π/4 When a = b the dc conductivity becomes s r = 2K r e 2 /h, where the Coulomb interaction parameter K r can be deduced from the experiments which ranges from 0.4 to 0.7 in semiconductor QWs [26,27]. Therefore, one determines the relative strength of the RSOI and DSOI from the dc conductivity. In our calculation, we take the electron effective mass m = 0.067m e , the strengths of SOIs a and b are about 1 × 10 -11 eVm which is the typical strength of SOI in conventional semiconductors [28]. Electrons only occupy the lowest subband in GaAs QW assuming the Fermi wavevector is k F = 0.01 nm -1 .

Conclusions
In conclusion, we investigate theoretically the effect of the interplay between the Coulomb interaction and the SOIs on the electron ground state and charge transport property of interacting QWs oriented along different crystallographic directions in different planes. We find that the ground state of electrons in the QWs can transit among the different phases, e.g., the SDW, CDW, SS, and MM, by tuning the plane and orientation of the QW, the strengths of SOIs and the Coulomb interaction. The anisotropy of the dc conductivity in an interacting QW is induced by the interplay between the Coulomb interaction and SOIs. This anisotropy enables us to detect the strengths of RSOI and DSOI, which are very important for the comprehensive understanding of spin decoherence and constructing all-electrical spintronic device.
Author details