Research Article

Realization of two-dimensional spin-orbit coupling for Bose-Einstein condensates

See allHide authors and affiliations

Science  07 Oct 2016:
Vol. 354, Issue 6308, pp. 83-88
DOI: 10.1126/science.aaf6689

Spin-orbit coupling in an optical lattice

Studying topological matter in cold-atom systems may bring fresh insights, thanks to the intrinsic purity and controllability of this experimental setting. However, the necessary spin-orbit coupling can be tricky to engineer. Wu et al. conceived and experimentally demonstrated a simple scheme that involves only a single laser source and can be continuously tuned between one- and two-dimensional spin-orbit coupling (see the Perspective by Aidelsburger). Although this experiment used bosonic atoms, it is expected that the setup would also work for fermions.

Science, this issue p. 83; see also p. 35


Cold atoms with laser-induced spin-orbit (SO) interactions provide a platform to explore quantum physics beyond natural conditions of solids. Here we propose and experimentally realize two-dimensional (2D) SO coupling and topological bands for a rubidium-87 degenerate gas through an optical Raman lattice, without phase-locking or fine-tuning of optical potentials. A controllable crossover between 2D and 1D SO couplings is studied, and the SO effects and nontrivial band topology are observed by measuring the atomic cloud distribution and spin texture in momentum space. Our realization of 2D SO coupling with advantages of small heating and topological stability opens a broad avenue in cold atoms to study exotic quantum phases, including topological superfluids.

The spin-orbit (SO) interaction of an electron is a relativistic quantum mechanics effect that characterizes the coupling between the motion and spin of the electron when moving in an electric field. In the rest frame, the electron experiences a magnetic field proportional to the electron velocity and couples to its spin by the magnetic dipole interaction, rendering the SO coupling. The SO interaction plays an essential role in topological insulators, which have been predicted and experimentally discovered in two-dimensional (2D) and 3D materials (1, 2), and topological superconductors (3, 4), which host exotic zero-energy states called Majorana fermions (5) and still necessitate rigorous experimental verification. For topological insulators, the strong SO interaction leads to band inversion, which drives topological phase transitions in such systems. In superconductors, triplet p-wave pairing may occur when SO coupling is present and results in topologically nontrivial superconductivity under proper conditions (6).

Recently, there has been considerable interest in emulating SO effects and topological phases with cold atoms, driven by the fact that cold atoms can offer extremely clean platforms with full controllability to explore such exotic physics. In cold atoms, the synthetic SO interaction can be generated by Raman coupling schemes that flip atom spins and transfer momentum simultaneously (7, 8); 1D SO interaction (9) has been successfully demonstrated in experiment for both cold boson (10, 11) and fermion degenerate gases (12, 13). With the 1D SO coupling, which corresponds to an Abelian gauge potential (14), one can study effects such as the magnetized or stripe ground states for bosons (1518), spin dynamics (19, 20), and 1D insulating topological states for fermions (21). Realizing higher dimensional SO couplings, which correspond to non-Abelian gauge potentials (7, 8), can enable the study of a broader range of nontrivial quantum states such as topological insulators driven by 2D and 3D SO interactions (1, 2). Furthermore, a 2D SO interaction is the minimal requirement to reach a gapped topological superfluid phase through a conventional s-wave superfluid state (22, 23).

Several schemes have been proposed for generating 2D and 3D SO couplings (7, 8, 2426). Notable progress was recently achieved when 2D SO couplings were demonstrated for pseudospins defined by two dark states in an empty tripod system (27). However, realizing 2D SO coupling for quantum degenerate atom gases remains challenging. Very recently, it was proposed that 2D SO coupling can be realized by a simple optical Raman lattice scheme that applies two pairs of light beams to create the lattice and Raman potentials simultaneously (28). However, this scheme requires the challenging realization of two Raman transitions with a locked relative phase. Here we propose a minimal scheme that overcomes these challenges and realizes 2D SO coupling with 87Rb Bose-Einstein condensates (BECs).

Theoretical proposal

We aim to realize 2D SO coupling and topological band for ultracold atoms on a square lattice with the HamiltonianEmbedded Image(1)where ħ is Planck’s constant h divided by 2π, k is the wave vector that represents the momentum of the atoms, 1 is the 2-by-2 unit matrix, Embedded Image are Pauli matrices acting on the spins, m is the mass of an atom, Embedded Image denotes the lattice potential in the xz plane, Embedded Image are periodic Raman coupling potentials, and Embedded Image is a tunable Zeeman constant. The lattice potential Embedded Image is spin-independent and can induce nearest-neighbor hopping that conserves the atom spin, whereas Embedded Image induce hopping that flips atom spin. This is different from the laser-assisted tunneling scheme without spin-flip (2931). The overall effect of hopping along Embedded Image and Embedded Image directions results in 2D SO coupling, which can lead to nontrivial topological bands for the square lattice.

Here we propose to realize the Hamiltonian (Eq. 1) through a minimal scheme, which is generic and applicable to both boson and fermion atoms. Figure 1 illustrates the realization in a 87Rb Bose gas, with Embedded Image and Embedded Image; the hyperfine state Embedded Image can be removed by a sufficiently large two-photon detuning. The minimal ingredients of the realization include a blue-detuned square lattice created with two light components denoted by the blue lines, and the periodic Raman potentials generated together with additional light components denoted by the red lines (Fig. 1A). Both ingredients can be achieved with a single in-plane (xz) linearly polarized laser source. The initial phases of the light beams have no effect on this optical Raman lattice scheme (32), and we neglect them in the following discussion. The optical lattice is generated by Embedded Image and Embedded Image (blue lines), which are incident from horizontal (x) and vertical (z) directions, respectively (Fig. 1A), and can be created from a single light of frequency Embedded Image by a beam splitter. The beams are reflected by two mirrors (M1 and M2) and form standing waves in the intersecting area described by Embedded Image and Embedded Image, whereEmbedded Image are amplitudes and the phase Embedded Image is acquired through the optical path L from the intersecting point to mirror M1, then to M2, and back to the intersecting point, with Embedded Image. For alkali atoms, we can show that the optical potential generated by linearly polarized lights is spin-independent when the detuning Δ is much larger than the hyperfine structure splittings (32). The lattice potential then takes the form Embedded Image(2)where Embedded Image. The Rabi frequency amplitudes of the standing waves Embedded Image, where Embedded Image, and the effective dipole matrix Embedded Image takes into account the transitions from a ground state Embedded Image to all relevant excited states in D1 and D2 lines (Fig. 1B). The lattice potential induces spin-conserved hopping, as illustrated in Fig. 1D.

Fig. 1 Proposal of the optical Raman lattice scheme.

(A) Sketch of the setup for realization. The light components Embedded Image (blue lines) form a spin-independent square optical lattice in the intersecting area and generate two periodic Raman potentials, together with the light components Embedded Image (red lines). (B) Optical transitions to generate lattice potentials by Embedded Image for the states Embedded Image Embedded Image, with Embedded Image and Embedded Image, including all relevant D1 and D2 transitions. Here, F is the quantum number of hyperfine states, mF is the quantum number of magnetic substates, and Δs denotes the fine-structure splitting. (C) Two periodic Raman potentials are generated through a double-Λ type configuration by Embedded Image and Embedded Image, respectively. (D to F) Profiles of lattice potential for Embedded Image (D) and Raman potentials Embedded Image (E) and Embedded Image (F) for Embedded Image. The color intensity characterizes the relative height of the potentials. The Raman potentials Embedded Image and Embedded Image are antisymmetric (or symmetric) with respect to the lattice site along the x (or z) and z (or x) directions, respectively.

The Raman couplings are induced when another beam Embedded Image of frequency Embedded Image is incident from the z direction. The light components Embedded Image and Embedded Image can be generated from a single laser source via an acoustic-optic modulator (AOM), which controls their frequency difference Embedded Image and amplitude ratio Embedded Image. The light Embedded Image generates plane-wave fields Embedded Image and Embedded Image, where the irrelevant initial phase is neglected (32). The relative phase Embedded Image, acquired by Embedded Image, is a crucial parameter, which can be precisely manipulated by changing the optical path L or Embedded Image, and it controls the dimensionality of the realized SO coupling. The standing-wave and plane-wave beams form a double-Λ type configuration (Fig. 1C), with Embedded Image and Embedded Image generating one Raman potential via Embedded Image in the form Embedded Image, and Embedded Image and Embedded Image producing another one via Embedded Image as Embedded Image. Note that lattice sites of a blue-detuned lattice are located around zeros of optical fields. It follows that terms such as Embedded Image, which is antisymmetric with respect to each site in both x and z directions, have small contribution to low-band physics. Neglecting such terms yieldsEmbedded Image(3)Embedded Image(4)Here Embedded Image, with Embedded Image (32). Together with an effective Zeeman term Embedded Image, which is controlled by tuning the two-photon detuning δ (Fig. 1C), we reach the effective Hamiltonian (Eq. 1). Note that Embedded Image Embedded Image is antisymmetric with respect to each lattice site along the Embedded Image Embedded Image direction (Fig. 1, E and F). This feature has an important consequence that Embedded Image Embedded Image leads to spin-flipped hopping only along x (z) direction. Moreover, the phase difference Embedded Image governs the relative strength of the Embedded Image and Embedded Image terms and thus determines the dimensionality of the SO coupling. For example, setting δω = 50 MHz yields Embedded Image for L = 1.5 m, resulting in optimal 2D SO coupling. Further increasing the optical path to L = 3.0 m gives Embedded Image, and the SO coupling becomes 1D. This enables a fully controllable study of the crossover between 2D and 1D SO couplings by tuning Embedded Image and provides a comparison measurement to confirm the realization of 2D SO interaction.

2D SO coupling and topological band

The Hamiltonian (Eq. 1) has an inversion symmetry defined by Embedded Image, where the 2D spatial operator Embedded Image transforms the Bravais lattice vector Embedded Image. For the s-band, the Bloch Hamiltonian is given by Embedded Image, which, around the Γ point, takes the form Embedded Image, with the SO coefficient Embedded Image being tunable by varying Raman coupling strength (32), unlike in the previous schemes (1013, 27). Here, t0 and tso denote the spin-conserved and spin-flip hopping coefficients, respectively, and a is the lattice constant. This is a quantum anomalous Hall model driven by SO coupling (28), which cannot be exactly realized in solid-state materials. It was shown (33) that the topology of inversion symmetric Chern bands can be determined by the product of the spin-polarizations Embedded Image at four highly symmetric momenta Embedded Image, with the momenta Embedded Image. The topological (or trivial) phase corresponds to Embedded Image (or +1). Two typical examples are shown in Fig. 2 by exactly diagonalizing H, with Embedded Image, Embedded Image, Embedded Image, and Embedded Image (Fig. 2, A, B, E, and F) or Embedded Image (Fig. 2, C and D), where the recoil energy is Embedded Image. For the chosen parameters, the lowest two subbands are gapped (Fig. 2, A to D). When Embedded Image, the spin polarizations at the Γ and M points are opposite (Fig. 2, B, E, and F), implying that the band is topologically nontrivial. In contrast, the polarizations are the same for Embedded Image (Fig. 2D), and the band is trivial.

Fig. 2 Band structure and spin texture with 2D SO interaction.

This figure shows an example of gapped band structure with nontrivial band topology (A), spin texture along the loop Γ-X-M-Γ (B), and spin polarization distributions Embedded Image of the lowest band (E) and the second band (F) for Embedded Image. (C and D) Example of a trivial band with gapped band structure (C) and spin texture along the loop Γ-X-M-Γ (D) for Embedded Image. For all panels, we take Embedded Image, Embedded Image, and Embedded Image.

The present scheme displays several essential advantages: (i) Fluctuations, such as those caused by the mirror oscillations, have a tiny effect on L and, thus, do not affect Embedded Image appreciably. The initial phases of light beams globally shift the optical Raman lattice but cannot affect the relative configuration between Embedded Image and the Raman potentials (32). Thus, the present scheme is intrinsically immune to any phase fluctuations in the setting, which avoids the need for phase-locking, a challenging task in practical realizations. (ii) As long as Embedded Image and Embedded Image, which are easily accomplished, the system becomes uniform in the x and z directions: Embedded Image and Embedded Image. No fine-tuning of optical potentials is required. (iii) All of the coupling beams can be created from only a single laser source, simplifying the experimental layout. (iv) Compared with the generation of 1D SO coupling, the present double-Λ Raman configuration (Fig. 1C) does not suffer additional instability or heating in the realization. These advantages render the present scheme immediately feasible in ultracold atom experiments with the current technology.

Experimental setup

In our experiment, a BEC of about 1.5 × 105 87Rb atoms in the state Embedded Image is prepared in a crossed optical dipole trap with trapping frequencies of Embedded ImageHz, which can suppress the antitrapping effect of blue detuned lights. A bias magnetic field of 49.6 G is applied along the Embedded Image direction to generate the Zeeman splitting and determine the quantization axis. As shown in Fig. 3A, three laser beams (wavelength = 767 nm) in the xz plane illuminate the atoms for the generation of the Hamiltonian (Eq. 1). Among these laser beams, a pair of counterpropagating lasers with the same frequency Embedded Image (the blue lines in Fig. 3A labeled as “lattice lasers”) produces the two-dimensional optical lattice. These two lasers are incident along the Embedded Image and Embedded Image directions, respectively, and are reflected by two mirrors (M1 and M2) to form the standing waves in both directions. The polarizations are set in the xz plane so that the interference between the Embedded Image and Embedded Image directions is automatically avoided. The third laser with frequency Embedded Image (the red line in Fig. 3A labeled as “Raman laser”) is a running wave, which is incoming along the Embedded Image direction with the same polarization as the lattice lasers. All three laser beams are generated from the same Ti:sapphire laser, and the frequencies and amplitudes of these beams are controlled by two phase-locked AOMs. Thus, the phase coherence is automatically kept, and no additional phase-locking is needed. The Raman and lattice lasers are also coupled into the same optical fiber and then lead to the science chamber, which helps to avoid the phase noise due to the imperfect overlap in propagation. The frequency difference Embedded Image is set to Embedded Image MHz to match the Zeeman splitting between the Embedded Image and Embedded Image states. The Embedded Image state is effectively suppressed because of a large quadratic Zeeman splitting, and the system can be treated as a two-level system. The detuning mz can be adjusted by tuning the bias magnetic field. By controlling the intensities of the lattice and Raman lights, we set the lattice depth as Embedded Image and the Raman coupling strength as Embedded Image.

Fig. 3 Experimental realization of 2D SO interaction and 1D-2D crossover.

(A) Experimental setup. B is the biased magnetic field, which generates the Zeeman splitting and gives the quantum axis of the atoms. (B) Spin-resolved TOF images of BEC atoms for Embedded Image, Embedded Image, Embedded Image, and Embedded Image. The other parameters are measured as Embedded Image, Embedded Image, and Embedded Image. (C) Measured imbalance W between the Raman coupling–induced atoms in the two diagonal directions as a function of the relative phase Embedded Image, compared to a cosine curve Embedded Image. The results are averaged over ~30 TOF images.

In the experiment, the BEC is first prepared in the dipole trap with the bias magnetic field being switched on. Then, the intensities of the lattice and the Raman beams are simultaneously ramped up to the setting value in 40 ms. As a consequence, the BEC atoms are adiabatically loaded in the local minimum of the lowest band at the Γ point. The phase difference Embedded Image in the Hamiltonian can be achieved by setting the propagating length between the two mirrors M1 and M2. The detection is performed in the same way as described in (11, 18, 34). The spin-resolved time-of-flight (TOF) imaging is taken after all laser beams and the bias magnetic field are suddenly turned off and the gas has expanded freely for 24 ms within a gradient magnetic field to resolve both the momentum and spin.

Experimental results

To demonstrate the realization of 2D SO coupling, we study the crossover effect in the BEC regime by tuning Embedded Image. At mz = 0 and by preparing the atoms in the spin-up state, we can adiabatically load the 87Rb condensate into the Γ point. Then, we perform the spin-resolved TOF expansion, which projects Bloch states onto free momentum states with fixed spin polarizations. Figure 3B shows the TOF images for various values of Embedded Image. For the spin-up Embedded Image state, five atom clouds are observed: The major portion of the BEC cloud remains at momentum Embedded Image, whereas four small fractions are transferred to momenta Embedded Image and Embedded Image by the first-order transitions due to the lattice potential Embedded Image. Depending on Embedded Image, two or four small BEC clouds are formed in the Embedded Image state at the diagonal corners with momenta Embedded Image. This is a consequence of SO coupling; the atom clouds are generated by the Raman transitions, which flip spin and transfer momenta of magnitude Embedded Image along the diagonal directions. As given in Eqs. 3 and 4, the Raman terms Embedded Image and Embedded Image depend on Embedded Image. For Embedded Image, four small clouds in the Embedded Image state with TOF momentum Embedded Image are observed (Fig. 3B), reflecting the 2D SO coupling. On the other hand, by tuning the relative phase to Embedded Image, the population of atom clouds in the two diagonal directions becomes imbalanced. Furthermore, the system reduces to 1D SO couplings when Embedded Image and Embedded Image, with Embedded Image and Embedded Image. In this case, the Raman pumping generates only a single diagonal pair of BEC clouds, as shown in Fig. 3B for Embedded Image. This is similar to the 1D SO coupling in the free space in (11), where the Raman coupling flips the atom spin and generates a pair of atom clouds with opposite momenta. Figure 3B also shows that there is a difference of distribution between the lower left and upper right BEC clouds at Embedded Image, which is due to non–tight-binding correction. A simple analysis reveals that although the fully antisymmetric Raman terms Embedded Image have negligible effects in the tight-binding limit of the lattice, they give finite contributions in the moderate lattice regime and are responsible for such difference of distribution (32). To quantify the crossover effect, we define Embedded Image to characterize the imbalance of the Raman coupling–induced atom clouds, with Embedded Image denoting the atom number of the two BEC clouds along the diagonal Embedded Image direction. W can be fitted by a simple cosine curve Embedded Image (Fig. 3C), reflecting the crossover between the 2D and 1D SO couplings realized in the present BEC regime.

Next, we focus on the 2D isotropic SO coupling with Embedded Image, measure the spin distribution in the first Brillouin zone, and detect the topology of the bands by varying the bias magnetic field to tune δ, which governs mz. For this purpose, we need a cloud of atoms with a temperature such that the lowest band is occupied by a sufficient number of atoms, whereas the population of atoms in the higher bands is small. A similar procedure used in the above BEC measurement is followed, except that the atoms are cooled to relatively higher temperatures, which are measured a posteriori using the momentum distribution of hot atoms. After a TOF expansion, we obtain the atom distributions of both spin-up and spin-down states in the momentum space and then map them back to the Bloch momentum space according to the plane-wave expansion of eigenfunctions. We define the spin polarization Embedded Image, with Embedded Image being the density of atoms of the corresponding spin state in the first Brillouin zone. Figure 4, A and B, shows the numerical results and experimentally measured spin polarizations at different temperatures, respectively, for Embedded Image, Embedded Image, and Embedded Image. In performing numerical simulation, the finite temperature effect is taken into account based on the Bose-Einstein statistics Embedded Image, with kB the Boltzmann constant and Eq given by band energies of the Hamiltonian (Eq. 1), plus the kinetic energy Embedded Image due to the motion in the out-of-lattice plane (y) direction. The average atom density is taken as n = 3 × 1019 m–3, which determines the chemical potential μ. There is good agreement between the theoretical and experimental results, which demonstrates the feasibility and reliability of the spin polarization measurement. Furthermore, the results in Fig. 4 suggest that a temperature around T = 100 nK is optimal to extract the spin-texture information of the lowest band. In comparison, if the temperature is too high, atoms are distributed over several bands and the visibility of the spin polarization will be greatly reduced, whereas too low a temperature can also reduce the experimental resolution because the atoms will be mostly condensed at the band bottom.

Fig. 4 Spin texture at different temperatures with Embedded Image.

(A) Numerical calculations for temperatures T = 250, 175, 125, and 75 nK. P is the polarization of the atomic spin, where +1 means spin up and –1 means spin down. (B) Experimental measurements of spin polarization at different measured temperatures. Parameters for the numerical calculation are the same as the experimental parameters: Embedded Image, Embedded Image, and Embedded Image.

We then measure the spin polarization as a function of detuning mz to reveal the topology of the lowest-energy band, with Embedded Image and Embedded Image. The numerical calculations and TOF-measured images of Embedded Image are in agreement (Fig. 5A and Fig. 5, B to D, respectively). In Fig. 5E, we plot the values of polarization Embedded Image for the four highly symmetric momenta Γ, X1, M, and X2. Embedded Image and Embedded Image always have the same sign, whereas the signs of Embedded Image and Embedded Image are opposite for small Embedded Image and the same for large Embedded Image, with a transition occurring at the critical value of Embedded Image that is a bit larger than Embedded Image. At transition points, the spin polarization Embedded Image or Embedded Image vanishes due to the gap closing and thermal equilibrium. From the measured spin polarizations, the product Embedded Image and the corresponding Chern number, given by Embedded Image (32, 33), can be read off (Fig. 5F). A numerical calculation using exact diagonalization (32) in the present experimental parameter regime can show two transitions between the topologically trivial and nontrivial bands near Embedded Image, according to the theory in (33), which agrees with the experimental observation. Note that around mz = 0, the spin polarizations at X1,2 change sign through zero, implying the gap closing at X1,2 and a change of Chern number by 2. This confirms that for the 2D SO-coupled system realized in the present experiment, the energy band is topologically nontrivial when Embedded Image, whereas it is trivial for Embedded Image.

Fig. 5 Spin texture and band topology with Embedded Image.

(A and B) Spin texture at different Embedded Image by tuning the two-photon detuning. Experimental measurements (B) are compared to numerical calculations at T = 100 nK (A). (C and D) Measured spin texture in topologically trivial bands at Embedded Image (C) and Embedded Image (D). (E and F) Measured spin polarization Embedded Image at the four symmetric momenta Embedded Image as a function of Embedded Image (E), and the product Embedded Image (F), which determines the Chern number Ch1 and characterizes the topology of the band. In all the cases, we set Embedded Image and Embedded Image.

Estimation of heating

The heating rate of the dipole trap is measured to be 18 nK/s, mainly owing to the photon scattering and the intensity noise. This results in a BEC lifetime of about 10 s. The heating rate of the lattice and the Raman lights caused by photon scattering is about four times that of the dipole trap, in the regime for Embedded Image and Embedded Image (32). Nevertheless, in the current experiment, residual heating is induced by the fluctuation of the bias magnetic field, which drives additional spin-flip dynamics in the presence of resonant Raman couplings. This contribution to the heating is about one order higher that of the dipole trap, reducing the lifetime of the SO-coupled BEC to just above 300 ms. This lifetime is sufficient to explore both single-particle and interacting physics for the 87Rb BEC system. Moreover, stabilizing the bias magnetic field may result in an even longer lifetime of seconds in appropriate parameter regimes.

Discussion and outlook

The 2D SO coupling we realized here is for real spins (hyperfine eigenstates) of atoms, which can be precisely measured and engineered experimentally. In comparison, a 2D SO coupling via a tripod system (7, 24, 27) or ring-coupling scheme (25) corresponds to pseudospins defined by superpositions of multiple hyperfine levels with superposition coefficients being spatially dependent. This conceptual difference manifests the advantages of the present realization for future broad studies of SO effects and interacting physics. Furthermore, owing to the realization in the optical lattice, the present 2D SO coupling can bring about much richer physics than a pure 2D Rashba correspondence. In the s-band regime, the present Bloch Hamiltonian describes a quantum anomalous Hall model driven by SO coupling, which cannot be exactly realized in solid-state materials. Thus, even in the single-particle regime, our realization leads to nontrivial topological bands, whereas a single-particle 2D Rashba system is topologically trivial. Moreover, even richer physics can be obtained if considering the higher-band (e.g., p-band) regimes.

Many experimental studies—including the measurement of topological Hall effects, Berry phase mechanism, and k-space monopole—can be performed on the basis of the present realization. On the other hand, with the high controllability of the present realization, the SO interaction can be readily switched on and off and can be adjusted between 1D and 2D limits. This may lead to rich quench spin dynamics in the optical lattice with nontrivial band topology. Moreover, with the present SO coupling in the optical lattice, one may explore states of matter [such as SO-coupled Mott insulators with interacting bosons (35, 36)] that have no analog in solids.

Furthermore, the present optical Raman lattice scheme is generic and can be immediately applied to fermion systems (e.g., 40K), in which case, the quantum anomalous Hall effect in the single-particle regime and topological superfluid (28) or novel magnetic phases (37) in the interacting regimes will be especially noteworthy. In particular, the topological superfluid phase is highly sought after because it hosts Majorana quasi-particles, which obey non-Abelian statistics (38) and have attracted attention in both condensed matter and cold atom physics (5). Finally, although the present study is focused on a 2D lattice system, generalizing our scheme to 3D optical lattices may lead to the realization of topological phases in 3D systems, including the Weyl topological semimetals (39, 40).


Supplementary Text

Figs. S1 to S6

References (41, 42)


  1. See supplementary materials for more details about model realization, SO coupling, nontrivial topological bands, band topological measurement, and estimation of heating.
  2. Acknowledgments: We thank J. Ho, T.-F. J. Poon, and G.-B. Jo for helpful discussions. This work has been supported by the Ministry of Science and Technology of China (under grants 2016YFA0301601 and 2016YFA0301604), National Natural Science Foundation of China, the CAS, the National Fundamental Research Program (under grant 2013CB922001), Fundamental Research Funds for the Central Universities (under grants 2030020028 and 2340000034), and Peking University Initiative Scientific Research Program. X.-J.L. acknowledges support from the National Natural Science Foundation of China (grant 11574008). X.-J.L. and L.Z. are also supported by the Thousand-Young-Talent Program of China.
View Abstract

Stay Connected to Science

Navigate This Article