Large spin-orbit coupling and helical spin textures in 2D heterostructure [Pb2BiS3][AuTe2].

Two-dimensional heterostructures with strong spin-orbit coupling have direct relevance to topological quantum materials and potential applications in spin-orbitronics. In this work, we report on novel quantum phenomena in [Pb2BiS3][AuTe2], a new 2D strong spin-orbit coupling heterostructure system. Transport measurements reveal the spin-related carrier scattering is at odds with the Abrikosov-Gorkov model due to strong spin-orbit coupling. This is consistent with our band structure calculations which reveal a large spin-orbit coupling gap of εso = 0.21 eV. The band structure is also characterized by helical-like spin textures which are mainly induced by strong spin-orbit coupling and the inversion symmetry breaking in the heterostructure system.

The motion of an electron through a lattice potential gives rise to a coupling between the electron’s momentum and its spin. This relativistic effect is known as spin-orbit coupling (SOC). SOC plays an important role in topologically nontrivial electronic structures123 as well as those which can realize Majorana fermions456. SOC also offers a unique route to tune the spin degree of freedom that has potential impact on spintronics devices789. Understanding SOC in confined systems such as interfacial boundaries and nanowires is particularly important in order to exert quantum control on electronic degrees of freedom in miniaturized devices. Extensive studies have been conducted in the 2D electron gas such as GaAs/Ga1-xAlxAs10 and HgTe/CdTe/HgTe quantum wells11. More emergent properties can be expected in similar systems with sufficiently strong SOC such as naturally formed crystalline heterostructures that adopt a structure similar to those of epitaxially grown heterostructure films. A strong advantage of naturally formed heterostructures is that a large number of organic and inorganic materials adopt this type of structure, for example the intergrowth chalcogenides12 and the inorganic-organic hybrid systems such as halides13. These crystal heterostructures provide a parallel platform to epitaxial films to uncover novel quantum phenomena derived from strong SOC and inversion symmetry breaking1011.

Recent work on the naturally formed 2D mineral [Pb2BiS3][AuTe2] shows heterostructure electronic states featuring a conductive [AuTe2] layer sandwiched between two insulating [Pb2BiS3] sheets14. The building blocks [Pb2BiS3]+1 and [AuTe2]−1 are alternately stacked along the crystal’s c-axis to form an intergrowth structure, shown in Fig. 1a 141516. A number of intriguing properties such as an extremely large electrical anisotropy, high mobility 2D carriers, and linear energy band dispersions were discovered14. The SOC strength in this material is expected to be strong, because many of the component materials are heavy elements. However, as we will discuss below, the SOC in [Pb2BiS3][AuTe2] is anomalous in that it lies beyond the expected scaling properties predicted by the Abrikosov-Gorkov model9. Furthermore, our calculations on atomically thin films unveil helical-like spin textures when the spin vectors traverse around the closed 2D Fermi surfaces. The helical-like spin textures could be an effect arising from the strong SOC and the unique heterostructure in [Pb2BiS3][AuTe2].

Figure 1

(a) The intergrowth structure of [Pb2BiS3][AuTe2] is constituted by two building blocks [Pb2BiS3]+ and [AuTe2]−, which stack alternatively along the c-axis. (b) Temperature dependent resistance and phase coherence length of a thin sheet crystal [Pb2BiS3][AuTe2]. (c) The field dependent conductance at elevated temperatures. The cusp shape curve at 1.6 K can be fitted with the 2D Hikami-Larkin-Nagaoka equation. (d) WAL of [Pb2BiS3][AuTe2] resembles that of Bi and TI Bi2Te3. The temperature is 2 K.

Results and Discussion

A direct measurement of a materials’ SOC strength is difficult because SOC is a relativistic effect. However, SOC manifests itself in quantum interference phenomena such as weak antilocalization (WAL)1718. WAL provides an effective approach to quantitatively investigate the SOC strength of a material using standard electrical characterization. In the diffusion regime and in the case of self-intersecting electron scattering paths, the interference between the electron’s time reversed paths is pronounced and tends to localize electrons. This phenomenon is known as weak localization (WL). WAL occurs when strong SOC suppresses WL effects.

Figure 1b shows the temperature dependent resistance of a cleaved single-crystal thin layer. The resistance upturn below 150 K demonstrates the poor metallic property of the samples. Below 40 K, the resistance increases linearly with decreasing temperature. This behavior could be related to electron-electron interaction that has been extensively discussed in disordered semiconductors and topological insulators1920. Our magnetoconductivity measurements focus on the low field regime, for example B < 0.6 T. In this regime, the Hall resistivity (ρxy) is far smaller than the longitudinal resistivity (ρxx) and the square conductance tensor takes the form . The magnetoconductance is defined as . Figure 1c is the field dependent conductivity at different temperatures. Both its cusp shape and its negative sign are hallmarks of WAL. Due to thermal fluctuations, the cusp shape gradually broadens with increasing temperature and diminishes at 40 K.

We apply the 2D Hikami-Larkin-Nagaoka (HLN) equation2122 to fit the data. Assuming that τ is much longer than both τso and τe, the HLN equation [see equation (2)] can be written

where α = 1/2, e is the electron charge and h is the Planck constant. The parameters τe, τso, and τ are the electronic elastic scatting time, spin-orbit scattering time and phase coherence time, respectively. The quantity Ψ(z) is the digamma function. We define a field strength B  = ħ/(4e) which is the characteristic field associated with the coherence length l, where and D is the diffusion constant. Figure 1c shows a fitting curve which virtually superposes on the experimental data. This fitting confirms that [Pb2BiS3][AuTe2] has strong SOC. At elevated temperatures (i. e. T = 35 K), l = 50 ± 10 nm, indicating robust quantum interference against thermal fluctuation in this system. Moreover, the agreement with a theoretical model of 2D transport confirms that the electronic degrees of freedom are confined to 2D in [Pb2BiS3][AuTe2]. Previous studies have shown large electrical anisotropy of Г~104 in this material14.

In Fig. 1d, we compare conductance data as a function of magnetic field of [Pb2BiS3][AuTe2] with that of Bi2Te323 and Bi24 at T = 2 K. For purposes of comparison, this data has been normalized to unity (see supporting information). The similarities among these various systems, especially between [Pb2BiS3][AuTe2] and Bi2Te3, is surprising. The topological insulator (TI) Bi2Te3 and semi-metal Bi are known for having unusually strong SOC among materials systems. This striking resemblance indicates that aspects of the SOC coupling in [Pb2BiS3][AuTe2] may be unusual. In order to gain a deeper understanding, we quantitatively investigated τso and τe of [Pb2BiS3][AuTe2] by analyzing the HLN equation that contains a form describing different scattering processes172122.

where B1 = Be + Bso + B, B2 = 4/3Bso + B + 2/3B, B3 = B + 2B, B4 = B2. The characteristic fields Bn,n = 1–4 are connected with the characteristic relaxation times τn by the relation Bn = ħ/4eDτn. The parameter Bis a characteristic field due to magnetic scattering and is zero in this case because [Pb2BiS3][AuTe2] is not a magnetic material with no magnetic impurities. It is worthy pointing out that Bergmann17 uses a special case of the formula derived by Hikami, Larkin, and Nagaoka21. Specifically, the corresponding decoherence times are the same for all three components of the spin. In addition, there is a typo in the HLN equation which has been corrected by Maekawa and Fukuyama22. Fitting data to equation (2) is challenging due to the large parameter space of Be, BSO and B. This problem can be addressed by using the value of B that was determined using equation (1), which assumes that τ ≪ τso, τe still holds at low temperatures. Our fitting yields Bso = 0.5 ± 0.2 T and Be = 1 ± 0.2 T. The reported Bso of Bi is 0.25 T24. The reported value of Bso of Bi2Te3 is greater than 8 T25. In contrast, the Bso of material systems with moderate SOC are generally far less than 0.1 T (see supporting information). The large Bso for these three systems qualitatively explains their similar WAL data. Our fitting parameters are reasonable: a value of Be = 1 ± 0.2 T corresponds to a mean free path of ~10 nm which is consistent with the sample’s poor metallic properties as revealed by the resistivity measurement presented in Fig. 1b.

Materials with large SOC strength may have potential applications as platforms for spin-orbitronics7. Comparison of the SOC related carrier scattering in various systems is thus worthwhile. Here we use the model of Abrikosov and Gorkov which predicts τe/τso ≈ (αZ)49 where α ≈ 1/137 is the fine structure constant and Z represents the atomic number. This model reveals that the electrons’ spin-orbit scattering is not only determined by the atomic number, but also is influenced by the elastic scattering of electrons by impurities. This is particularly important for spin-orbitronics applications since it predicts that τso can be tuned by the material’s purity. Here, we use the weighted arithmetic mean to calculate the effective atomic number of materials composed of multiple elements. For [Pb2BiS3][AuTe2], only the elements Au and Te are accounted for in the effective Z since the density of states at the Fermi surface are mainly contributed by Au 3d electrons and Te p electrons14. Figure 2 summarizes the relations between τe/τso and the effective atomic mass of a variety of systems whose SOC strength has been investigated. Consistent with the Abrikosov-Gorkov prediction, a majority of the materials follow the trend that τe/τso is proportional to Z4, which is depicted by the red solid line in Fig. 2. The topological insulator Bi2Te3 is an exception due to its extremely large SOC26, the topology-related Berry’s phase and the backscattering immunity27 which always make τso extremely short and τe exceptionally long. In topological insulators, the elastic scattering rates can be weaker than the spin orbit scattering rates825, which results in τe/τso being greater than unity. The topological insulators are thus considered to be excellent spin-generator materials7. We categorize the TI material Bi2Te3 into the unconventional class of systems and the rest such as Au28 into the conventional class. Our [Pb2BiS3][AuTe2] sample falls between these two classes, suggesting the strong SOC, heterostructured systems may have potential for spin generation. The violation of the Abrikosov-Gorkov scattering model in our system is intriguing, and must be related to strong SOC.

Figure 2

τe/τso versus effective atomic number (Z) for different systems.

Weighted arithmetic mean is used to calculate the effective atomic number of materials composed of multiple elements. For [Pb2BiS3][AuTe2], only the elements Au and Te are accounted for the Z because the density of states at the Fermi surface are mainly contributed by Au 3d electrons and Te p electrons14. The red solid line is a simulation using the relation Z4. TI Bi2Te3 violates the Abrikosov and Gorkov’s prediction τe/τso~ Z4 due to the topology-related Berry’s phase and the extremely large SOC. [Pb2BiS3][AuTe2] does not follow the Z4 relation. The τso and τe of Al, InP/InGaAs, In2O3, Hg0.76Cd0.24Te, Au and Ag are from references2835363738.

The large SOC can be understood by our ab initio calculations. The electronic structure near the Fermi level of [Pb2BiS3][AuTe2] purely originates from the [AuTe2]− layer. Figure 3a is a contour plot of the local density of state (LDOS) of [Pb2BiS3][AuTe2] near the Fermi level, which clearly shows that the electronic states are tightly confined to the [AuTe2]− layer. Figure 3b exhibits the energy band dispersions of the [AuTe2]− layer and the single-unit-thick films. The vertical arrow marks the gap opening due to strong SOC. Figure 3c shows the detailed information of the gap opening in 3D reciprocal space. Without SOC in the [AuTe2]− layer, a Dirac-like gapless state appears at the hole-band between the Г and X points. SOC opens the gapless Dirac-like state to form a spin-orbit gap of εSO = 0.21 eV. This value is comparable to the SOC-induced band inversion energy in TIs and the spin-orbit gap of the topology-nontrivial Sb-bilayer system2629. The SOC gap in the [AuTe2]− layer also occurs at the Г point near the Fermi Energy (see the top panel of Figure S6 in supporting information). The large spin-orbit gap elucidates the large SOC and the violation of the Abrikosov-Gorkov scattering model. Since [Pb2BiS3][AuTe2] is not a TI, our studies thus suggest that materials can violate the Abrikosov-Gorkov scattering model as long as their SOC strengths are sufficiently large.

Figure 3

(a) Contour plot of the local density of state (LDOS) of [Pb2BiS3][AuTe2] near the Fermi level. The [AuTe2]− layer dominates the electronic structure. (b) Energy bands dispersion of one-unit-cell [Pb2BiS3][AuTe2] film (light blue), and [AuTe2]− single layer (light red). A hole pocket lies at Г and an electron pocket is formed between Г-Y. (c) A Dirac-like gapless state appears at the hole-band without SOC (upper panel). SOC opens the Dirac-like gapless state with a spin-orbit gap of 0.21 eV (lower panel). The gap opening is denoted in (b) by the black arrow.

The large SOC, combined with the unique heterostructure in [Pb2BiS3][AuTe2], may lead to helical spin textures that have been observed in topological insulators and semiconductor heterostructure films2730. We carried out theoretical investigations on the single-unit-cell film [Pb2BiS3][AuTe2] films where the twofold spin-degeneracy is lifted due to inversion symmetry breaking. Figure 4a presents the Fermi surfaces of the electron and hole bands and their projections at the Fermi level. The dual projection sheets are induced by broken inversion symmetry. Both electron bands and hole bands are found to have spin structures whose direction changes as a function of crystal momentum k. The spin textures of the electron and hole pockets exhibit different topological properties. For the electron pockets, the spin does not undergo a full winding as one traverses around the Fermi surface while for the hole pockets the spin winds around 3 times. The helical-like spin texture in the hole pocket can be attributed to the strong SOC and the inversion-asymmetry induced spin splitting in a heterostructure system. For example the Ga1-xInxAs/InP quantum well has shown helical spin texture in the presence of SOC30. The spin textures in the hole pocket exhibit spin-flip between the hole carriers with opposite momentum k and –k, as a consequence of time-reversal symmetry. The shape of the Fermi surface reflects the orthorhombic crystalline symmetry of [Pb2BiS3][AuTe2]. In addition, the spin-flip in the hole pocket may indicate a non-vanishing Berry’s phase. Direct detection of the Berry’s phase is beyond the scope of this study. However, the WAL in our sample is consistent with a nontrivial Berry’s phase which gives opposite signs to time-reversed electron paths. This type of destructive quantum interference can also result in WAL182731. Our modeled helical spin textures in single-unit-cell films could be experimentally pursued as we demonstrate that single-unit-cell films naturally occur on the cleaved surfaces of a bulk crystal (Figure S1 in supporting information). Further experiments such as spin-resolved photoemission spectroscopy are needed to confirm our theoretical predictions of the helical spin textures.

Figure 4

(a) A 3D plot of the electronic band structure of the one-unit-cell [Pb2BiS3][AuTe2] film near the Fermi level. The 2D projection of electron bands and hole bands are depicted on the top and at the bottom, respectively. Black solid lines indicate the Fermi level. The dual projected sheets are due to inversion symmetry breaking on the surfaces of the one-unit-cell film. (b,c) shows spin direction along the closed 2D Fermi surfaces for electron and hole bands, respectively. Left and right panels of (b,c) correspond to two non-degenerated bands due to inversion symmetry breaking in the one-unit-cell film. The area of (b,c) in momentum space are denoted by purple-colored rectangles shown in (a).

Conclusions

To summarize, large SOC strength was discovered in the 2D heterostructure [Pb2BiS3][AuTe2]. This large SOC is induced by a large spin-orbit gap εSO = 0.21 eV. More broadly, our work suggests that materials can violate the Abrikosov-Gorkov scattering model for sufficiently strong spin-orbit coupling. Ab initio calculations reveal helical-like spin textures and spin-flips at the Fermi surfaces. These predictions can be attributed to the effect of strong SOC and inversion symmetry breaking in the heterostructure as well as the constraint of crystalline symmetry in the AuTe2 layer. More generally, our work points out that naturally forming heterostructures made of heavy atoms provide a new direction for exploring novel quantum phenomena at the atomic scale. Given the large number of naturally formed organic and inorganic heterostructures and their hybrids1213, more discoveries can be expected in this direction.

Methods

Bulk crystal synthesis follows a self-flux method by melting stoichiometric compositions of [Pb2BiS3][AuTe2]. Details of the crystal synthesis were reported in our previous paper14. Crystals were exfoliated on SiO2/Si wafers using the so-called scotch tape method and were post-annealed at 350 °C in argon gas to remove chemical residuals on the cleaved surface. The crystals’ thickness is determined by AFM (PSIA Corporation). Morphology studies on cleaved surfaces were conducted by a MultiMode scanning probe microscope (Veeco) operated in a peak-force tapping mode that has vertical spatial resolution up to 50 pm. Cleaved crystals with thickness 40 nm were selected for transport characterization. Standard photolithography (LaserWirter, MICROTECH) and magnetron sputtering were employed to pattern contacts on cleaved crystals. In-plane electrical characterization was conducted using the standard four probe method in a liquid 4He variable temperature cryostat equipped with a triple-axis vector magnet system (AMI).

Electronic structures of [Pb2BiS3][AuTe2] bulk, one-unit-cell thick [Pb2BiS3][AuTe2] film and the [AuTe2]− single layer were calculated within density functional theory using Projector Augmented Wave method32 implemented in the VASP code33. Perdew-Burke-Ernzerhof type34 generalized gradient approximation was utilized for the exchange correlation functional. Spin-orbit coupling was included in the non-collinear form. Experimentally resolved crystal structure was employed for [Pb2BiS3][AuTe2], and crystal structure of [AuTe2]− single layer was prepared from the bulk crystal structure of [Pb2BiS3][AuTe2] by removing [Pb2BiS3]+ layers. The crystal structure of the monolayer [Pb2BiS3][AuTe2] is calculated assuming a vacuum layer of 3 nm.

Additional Information

How to cite this article: Fang, L. et al. Large spin-orbit coupling and helical spin textures in 2D heterostructure [Pb2BiS3][AuTe2]. Sci. Rep. 6, 35313; doi: 10.1038/srep35313 (2016).