Topological nodal chains in optical lattices

Bei-Bei Wang,1, 2 Jia-Hui Zhang,1, 2 Chuanjia Shan,3, ∗ Feng Mei,1, 2, † and Suotang Jia1, 2 1State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Laser Spectroscopy, , , Shanxi 030006, 2Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China 3College of Physics and Electronic , Hubei Normal University, Huangshi 435002, China (Dated: March 23, 2021) Topological nodal rings as the simplest topological nodal lines recently have been extensively studied in optical lattices. However, the realization of complex nodal line structures like nodal chains in this system remains a crucial challenge. Here we propose an experimental scheme to realize and detect topological nodal chains in optical Raman lattices. Specifically, we construct a three-dimensional optical Raman lattice which supports next nearest-neighbor spin-orbit couplings and hosts topological nodal chains in its energy spectra. Interestingly, the realized nodal chains are protected by mirror symmetry and could be tuned into a large variety of shapes, including the inner and outer nodal chains. We also demonstrate that the shapes of the nodal chains could be detected by measuring spin polarizations. Our study opens up the possibility of exploring topological nodal-chain semimetal phases in optical lattices.

I. INTRODUCTION of-the-art optical lattice technology developed in these experiments further paves the way for exploring high- dimensional gapless topological phases [44]. In particu- The study of topological nodal line semimetals have lar, a lot of optical lattice models have been constructed recently emerged as a new frontier in condensed mat- for realizing topological nodal rings [45–51]. Since optical ter physics [1,2]. Nodal line semimetals are systems lattices naturally have clean environments, the realized where the conduction and the valence bands touch along nodal rings are exactly located at the Fermi level and not closed nodal lines in the Brillouin zone [1,2]. Intrigu- mixed with trivial bulk bands. However, how to realize ingly, such gapless phase also possesses topological prop- more exotic nodal line structures like nodal chains in this erties, such as quantized Berry phase and topologically system is still unknown. protected drumhead surface states [3]. Nodal lines have numerous topological structures than nodal points, in- Very recently, two experiments reported the realiza- cluding nodal rings [4–11], nodal chains [12–19], nodal tion and detection of topological nodal lines [52] and links [20–24] and nodal knots [25, 26]. In particular, the Weyl points [53] with ultracold atoms trapped in opti- nodal chains are formed by the nodal rings located on cal Raman lattices. Motivated by these experiments, we mutually orthogonal high-symmetry planes and touching present a three-dimensional optical Raman lattice sys- each other at isolated points [12–14], which are the basics tem that could generate next nearest-neighbor spin-orbit of generating nodal links and knots [14, 19]. Thus, nodal couplings, which has not been reported before. More- chains have much richer topological properties waiting to over, the corresponding tight-binding Hamiltonian has be explored and observed, such as the unusual topological not been previously discovered in solid state material or electromagnetic responses [12] and distinct topological optical lattice systems. Distinct from previous nodal-ring surfaces states [14, 17, 18]. Experimental efforts recently studies in optical lattices [45–51], our system could host have been devoted to searching materials with different a lot of nodal rings which are located on mirror invari- nodal lines [27, 28]. Unfortunately, the nodal lines show- ant planes and protected by the mirror symmetry. By ing there are often not exactly located at the Fermi level tuning the laser intensities and frequencies to vary the as the ideal one, mixed with the trivial bulk bands, mak- hopping rates and effective Zeeman field, the nodal rings ing the experimental detection much challenge. on the mutually orthogonal mirror invariant planes are connected together to form various shapes of topological Meanwhile, ultracold atoms trapped in optical lat- nodal chains, including the inner and outer nodal chains. tices provides a clean and controllable platform for re- In addition, we also demonstrate that the shapes of the alizing, manipulating and probing topological phases of arXiv:2103.11291v1 [cond-mat.quant-gas] 21 Mar 2021 nodal chains could be detected through the measurement matter [29–31]. This has been well demonstrated by the of spin polarizations. recent experimental implementation of spin-orbit cou- plings [32–38] and topological insulator phases [39–43] The paper is organized as follows. SectionII presents in low-dimensional ultracold atomic systems. The state- an optical Raman lattice model with next nearest- neighbor spin-orbit couplings. SectionIII shows the en- ergy spectra hosts mirror symmetry protected topolog- ical nodal rings and chains. SectionIV discusses the ∗Electronic address: [email protected] experimental detection of the nodal chains. SectionV †Electronic address: [email protected] summarizes the main results of this paper. 2

Spin hopping 푒 푧 Vsoσx, where Vso = Ωso cos(k0x) cos(k0y) cos(k0z) and (푎) (푏) 8 푒 −푡푠표 1 Ω = Ω¯ Ω¯ /∆. Then the total Hamiltonian of the sys- 푡푠표 so 1 2 푒ۧ 푦 tem reads| ∆ 푡푠표 푧 푒 푒6 Ω 푦 3 2 1 Ω2 p 푥 −푡푠표 X u  푡 Hs = + Vu + Vsoσx + mzσz, (1) 푒2 푠표 푥 2m 훿 u=x,y,z ۧ↑| −푡푠표 푒7 ۧ↓| 푡푠표 −푡푠표 where mz = δ/2 is the effective Zeeman field and δ = 푒 푒4 5 ω1 + ω0 − ω2 is the two-photon Raman detuning. Spin-orbit coupling In the second quantization, the continuum Hamilto- nian Hs is written into the following form FIG. 1: (a) Schematics of the laser configuration to realize the two-photon Raman푒Ԧ푧 coupling. Such coupling is generated Z 푒Ԧ 푒Ԧx † by two laser beams푦 with Rabi frequencies Ω1 and Ω2. ∆ is H = drψ (r)Hsψ(r) (2) the detuning from the auxiliary excited state and δ is the two-photon detuning in the Raman transition. (b) In the where ψ(r) = (ψ (r), ψ (r))T . Here we only study cubic optical lattice with the Raman coupling, the spin-orbit ↑ ↓ physics in the ground band. Then the field operator can couplings are generated only between the lattice sites j and 0 be expanded as j = j + ev (v = 1, 2, 3, 4, 5, 6, 7, 8). Moreover, the generated synthetic spin-orbit couplings are position-dependent, i.e., the X strengthes of the spin-orbit couplings along ev for v = 1, 2, 3, 4 ψσ(r) = W (r − j)Cjσ, (3) (v = 5, 6, 7, 8) are tso (−tso). j

where W (r − j) is the spin-independent Wannier func- tion for the ground band located at the lattice site II. OPTICAL RAMAN LATTICES WITH NEXT NEAREST-NEIGHBOR SPIN-ORBIT j = (jx, jy, jz). Throughout this work, we assume the COUPLINGS lattice spacing a = π/k0 = 1. Cjσ is the annihilation op- erator for an atom with spin σ =↑, ↓ at the lattice site j. With the above expansion, we can get the tight-binding We consider using the experimental setup reported for Hamiltonian observing ideal Weyl points [53]. This setup is made up of ultracold atoms trapped in optical Raman lattices. In X X † † H = − tu(Cj↑Cj+eu↑ + Cj↓Cj+eu↓ + h.c.) addition to conventional optical lattice trapping poten- j u=x,y,z tials, an effective Raman lattice potential is also applied X j,j0 † X † † in optical Raman lattice systems [54–56]. Suppose a spin + tso (Cj↑Cj0↓ + h.c) + mz (Cj↑Cj↑ − Cj↓Cj↓), is encoded by two atomic internal states. The effective j,j0 j Raman potential is generated through a two-photon Ra- (4) man transition which couples the spin up and the spin down [54–56]. Here we find that, by changing the Raman where the vectors ex = (1, 0, 0), ey = (0, 1, 0) and ez = potential applied in [53], the resulted optical Raman lat- (0, 0, 1), the nearest-neighbor hopping rate and the spin- tice could generate next nearest-neighbor spin-orbit cou- orbit coupling strength are expressed as plings and support various topological nodal rings and Z 2 ∗ pu chains. tu = − duWu (u)( + Vu)Wu(u − 1), (5) The experimental setup is shown in Fig.1(a), where 2m 0 0 87 j,j Y ju,ju ultracold Rb atoms are trapped in a three-dimensional tso = Ωso tso,u , (6) optical Raman lattice. The spin is encoded by two mag- u=x,y,z netic sublevels |↑i = |1, −1i and |↓i = |1, 0i [53], with an energy difference ω0. The conventional cubic optical lat- with tice is generated by applying three standing wave lasers 0 Z ju,ju ∗ 0 in the three real space directions, leading to a trapping tso,u = du Wu (u − ju) cos(k0u)Wu(u − ju) 2 potential (spin-independent) Vu=x,y,z = V¯u cos (k0u), 2 Z where V¯u ∝ |Eu| and Eu is the amplitude of the stand- ∗ 0 = du Wu (u) cos(k0(u + ju))Wu(u + ju − ju) ing wave laser in the u direction. Furthermore, the two-photon Raman transition between the spin up and Z = −(−1)ju du W ∗(u) sin(k u)W (u + j − j0 ). the spin down is generated by applying two standing u 0 u u u wave lasers with frequency ω1 in the x-y plane and one (7) standing wave laser with frequency ω2 in the z direc- tion. The corresponding Rabi frequencies are Ω1 = Here we use the identities W (r − j) = Wx(x − jx)Wy(y − 2 Ω¯ 1 cos(k0x) cos(k0y) and Ω2 = Ω¯ 2 cos(k0z) respectively, jy)Wz(z − jz) and cos (k0ju) = 0. The last identity leading to a Raman coupling potential (spin-dependent) is obtained because the bottom of the lattice trapping 3 potential is the lattice site. In optical lattices, we usually spin-orbit couplings (see Eq. (10)). Thanks to these par- 0 only consider ju − ju = 0, ±1, because the integral in ticular next nearest-neighbor spin-orbit couplings, as we 0 Eq. (7) for ju − ju > 1 is much smaller than the one for will demonstrate below, the nodal chains protected by 0 ju − ju = 0, ±1. mirror symmetry are generated in this system. Considering that the Wannier function Wu for the ground band has an even parity, basing on Eq. (7), we obtain III. MIRROR SYMMETRY PROTECTED NODAL RINGS AND CHAINS ju,ju tso,u = 0, (8) tju,ju+1 = −tju,ju−1. (9) Before seeing the nodal chains emerged in the so,u so,u energy spectra, we transfer the real space tight- Combining Eq. (6) and Eq. (8), we further get binding Hamiltonian into the momentum space to P † get H(k) = k c (k)h(k)c(k), where c(k) = T j,j [c↑(kx, ky, kz), c↓(kx, ky, kz)] , tso = 0,

j,j+ex j,j+ey j,j+ez tso = tso = tso = 0, h(k) = hy(k)σy + hz(k)σz, (13) tj,j+ex+ey = tj,j+ex+ez = tj,j+ey +ez = 0. (10) so so so with From the above equations, we can see that the spin-orbit hy = 8tso sin(kx) sin(ky) sin(kz), couplings are mainly generated between j = (jx, jy, jz) 0 hz = mz − 2tx cos(kx) − 2ty cos(ky) − 2tz cos(kz). (14) and j = (jx ± 1, jy ± 1, jz ± 1). Specifically, as shown in Fig.1(b), the spin-orbit couplings act only along the The Bloch Hamiltonian h(k) respects both the time- vectors e = −e = (1, 1, 1), e = −e = (−1, 1, −1), 1 5 2 6 reversal symmetry and the inversion symmetry, i.e., e3 = −e7 = (−1, −1, 1) and e4 = −e8 = (1, −1, −1). Combining Eq. (6) and Eq. (9), we also can find that the T h(k)T −1 = h(−k) generated spin-orbit couplings are position-dependent −1 and have the following relationship Ph(k)P = h(−k) (15)

j,j+e j,j+e where T = K is the time-reversal symmetry opera- t s = −t t , (11) so so tor, with K being the complex conjugate operator, and where s = 1, 2, 3, 4 and t = 5, 6, 7, 8. By substituting Eqs. P = σz is the inversion symmetry operator. The time- (6,7,11) into Eq. (4) and performing a gauge transforma- reversal symmetry requires hx,z being even function of k jx+jy +jz and hy being odd function of k. In contrast, the inver- tion (−1) Cj↓ → Cj↓ [54, 55], the tight-binding Hamiltonian becomes sion symmetry requires hz being even function of k and hx,y being odd function of k. Then the consistence of X X † † these two symmetries in our system guarantees that hx H = − tu(Cj↑Cj+eu↑ − Cj↓Cj+eu↓ + h.c.) j u=x,y,z vanishes in the entire Brillouin zone. q 2 2 4 The energy spectra of h(k) is given by E = ± hy + hz X X † † + tso(Cj↑Cj+ev ↓ − Cj↑Cj−ev ↓ + h.c.) and consists of two energy bands. When hy(kx, ky, kz) = j v=1 hz(kx, ky, kz) = 0, the conduction and the valence bands X could cross each other to form a chain of connected nodal + m (C† C − C† C ), (12) z j↑ j↑ j↓ j↓ rings. For our system, the nodal rings are generated on j the mirror invariant planes and protected by the mirror symmetry. The reason is that the Bloch Hamiltonian where t = Ω t01 t01 t01 . The hopping rates, the so so so,x so,y so,z h(k) also respects the mirror symmetry spin-orbit coupling strengths and the effective Zeeman fields all can be flexibly tuned by changing laser intensi- −1 Mh(kx, ky, kz)M = h(−kx, ky, kz) ties and frequencies. = h(k , −k , k ) In contrast to previous experiments producing nearest- x y z neighbor spin-orbit couplings [42, 43, 52, 53], our pro- = h(kx, ky, −kz), (16) posed optical Raman lattice could generate the next nearest-neighbor spin-orbit couplings between the lattice where M = σz is the mirror symmetry operator. Note 0 that there are six mirror invariant planes kx,y,z = 0, π. sites j = (jx, jy, jz) and j = (jx ± 1, jy ± 1, jz ± 1), which are out of the x, y and z planes (see Fig.1(b)). The Now we briefly prove the existence of nodal rings on the mirror invariant planes kz = 0, π. The mirror symmetry reason is that, the Wannier function Wu=x,y,z for the ground band has an even parity and the Raman poten- in the z direction in Eq. (16) could lead to tial V centered at the lattice site has an odd parity (see so h (k , k , k ) = −h (k , k , −k ) Eq. (7)), which leads to the vanishing of the integrals in y x y z y x y z Eq. (6) for the on-site spin flips and the nearest-neighbor hz(kx, ky, kz) = hz(kx, ky, −kz), (17) 4

(푎) 휋 푘푧 푘z (푐) (푏) 휋 푘푦 (푑) Max 휋

푘푦 휋 0 휋 0 푬/풕 0 0 −휋 푘z 푘푦 푘 −휋 휋 푥 −휋 휋 휋 푘푦 0 푘z 0 휋 푘 푘푥 0 푘푥 푥 0 0 −휋 휋 0 휋 휋 −휋 0 −휋 −휋 −휋 0 −휋 −휋 −휋 (푒) (푓) (푔) (ℎ) 6 2 1 0.1 퐸/t 퐸/t 퐸/t 퐸/t −0.1 휋 −6 −2 −1 휋 휋 휋 휋 휋 휋 푘 푘 휋 푘 푦 푘 푦 푘 푘푦 푘푥 푦 푘 −휋 −휋 푥 −휋 −휋 푥 −휋 −휋 −휋 −휋 푥

(푖) 휋 (푗) 휋 (푘) 휋 (푙) 휋

휐 =-1 휐 =1 휐 =-1 휐 =1 휐 =-1 휐 =1 휐 =-1 휐 =1 푘 푘 푦 0 푦 0 푘푦 0 푘푦 0 휐 =1 휐 =-1 휐 =1 휐 =-1 휐 =1 휐 =-1 휐 =1 휐 =-1

−휋 −휋 −휋 −휋 −휋 0 휋 −휋 0 휋 −휋 0 휋 −휋 0 휋 푘푥 푘푥 푘푥 푘푥

FIG. 2: The nodal rings and chains in the energy spectra for (a) tx = ty = t, mz = 3t, (b) tx = ty = t, mz = 0, (c) tx = t, ty = 0, mz = 2t and (d) tx = ty = mz = 0. Specifically, the nodal rings (gap-closing points marked by black solid lines) are generated on the mirror invariant planes kx,y,z = 0, π. The topological features of nodal rings are characterized by the quantized Berry phases along a tiny loop enclosing the nodal lines. Our numerical results show that the quantized Berry phase is π (0) when the loop encloses one (two) nodal lines. By tuning the hopping rates and effective Zeeman field (see (a-d)), the nodal rings on the mutually orthogonal mirror invariant planes can be connected together to generate various shapes of nodal chains. With open boundary condition along the z direction and periodic boundary condition along the x and y directions, the corresponding energy spectra for the systematic parameters chosen in (a-d) are shown in (e-h) respectively, where the zero-energy drumhead surface states appear at the two surfaces of z-direction. The corresponding topological winding numbers defined in the Brillouin zone of z-direction are shown in (i-l) respectively. The other parameters are tso = tz = t and t is the energy unit.

which indicates that hy(kx, ky, kz) = 0 on the mir- X point is generated on the mirror invariant planes ror invariant planes kz = 0, π. Then the solution of kx,y,z = π; For −6t < mz < −2t, one nodal ring cen- hz(kx, ky, kz) = 0 determines the two bands crossing tered around the R point is generated on the mirror in- along rings located at the mirror invariant planes kz = variant planes kx,y,z = π, while there is no nodal ring on 0, π. Moreover, as we will show, these nodal rings are cen- the mirror invariant planes kx,y,z = 0. Here we assume tered around the high symmetry points. Since the mirror tx = ty = tz = tso = t. operator and the Hamiltonian have same eigenstates in As shown in Fig.2(a-d), via tuning the laser intensities the mirror invariant planes, the generated nodal rings and frequencies to vary tx, ty and mz, the nodal rings on are protected by the mirror symmetry [1]. Similarly, one the mutually orthogonal mirror invariant planes could be also can prove that nodal rings also emerge on the mirror connected together to form various nodal chains, includ- invariant planes kx,y = 0, π. ing the outer and inner nodal chains. The outer (inner) Specifically, for 2t < mz < 6t, one nodal ring cen- nodal chains are generated when two nodal rings are on tered around the Γ point is generated on the mirror in- opposite (same) sides of the touching point [14]. For the variant planes kx,y,z = 0, while there is no nodal ring outer case, there are three types of outer nodal chains, generated on the mirror invariant planes kx,y,z = π; For including nodal nets, nodal tubes and nodal cross-lines. −2t < mz < 2t, four nodal rings centered around the For example, Fig.2(a) shows that an inner nodal chain M point are generated on the mirror invariant planes is generated by three connected nodal rings on the mu- kx,y,z = 0, and one nodal ring centered around the tually orthogonal mirror invariant planes kx,y,z = 0. Fig. 5

휋 휋 1 휋

휎z 휎z 푘푦 푘푦 푘푦 푎 푘푧0 = 2휋/3 (푏) 푘푧0 = 휋/2 (푐) 푘푧0 = 0.28휋 휎z

푘푧 푘푦 푘푥 −휋 −휋 −1 c 휋 c c −휋 휋 1 휋

푘 푘 푘푦 푦 휎푦 푦 휎푦 휎푦

−휋 −휋 −1 −휋−휋 휋 −휋 휋 −휋 휋 푘푥 푘푥 푘푥 휋 휋 휋

휎z 휎z 푘푦 푘푦 푘푦 휎 푑 푘푧0 = −0.2휋 (푒) 푘푧0 = −0.6휋 푓 푘푧0 = −2휋/3 z

−휋 −휋 c −휋 c 휋 c 휋 휋

푘 푦 휎 푘푦 휎 푘푦 푦 푦 휎푦

−휋 −휋 −휋 휋 −휋 푘푥 −휋 휋 −휋 휋 푘푥 푘푥

FIG. 3: Measuring spin polarizations hσy,zi to detect an inner nodal chain. The left side in (a-d) shows the shape of an inner nodal chain (formed by three connected nodal rings on the mirror invariant planes kx,y,z = 0) and the corresponding nodal points (marked by black points) when kz is fixed to kz0. The right side in (a-d) shows the spin polarizations hσy,zi as a function of kx and ky for (a) kz0 = 2π/3, (b) kz0 = π/2, (c) kz0 = 0.28π, (d) kz0 = −0.2π, (e) kz0 = −0.6π and (f) kz0 = −2π/3, where hσyi = 0 (marked by orange straight solid lines) and hσzi = 0 (marked by orange rings) respectively give the locations of BIS1 and BIS2 defined in Eqs. (20, 21, 23). As we can see, the nodal points for kz = kz0 can be detected through the intersection of the BIS1 and BIS2. In this way, the shape of the inner nodal chain is mapped out by measuring the vanishing of spin polarizations hσy,zi varying with kz. The other parameters are tx = ty = tz = tso = t and mz = 3t.

2(b) shows that an outer nodal chain shaped like a net the surface states appear as flat bands at the zero energy. is generated on the mutually orthogonal mirror invariant Distinct from the drumhead surface states for nodal-ring planes kx,y,z = 0 and kx,y,z = π. Fig.2(c) shows that an phases, the surface states for nodal-chain phases have outer nodal chain shaped like a tube is generated on the abundant shapes. In particular, for the nodal cross-line mutually orthogonal mirror invariant planes kx,z = 0 and case shown in Fig.2(d), the regions where the surface ky = 0, ±π, where the nodal rings on the mirror invariant states exist cover the whole kx-ky plane. planes kx,z = 0 are reduced to straight lines. Fig.2(d) These surface states are protected by nontrivial topo- shows that an outer nodal chain shaped like cross-lines logical winding numbers defined in the Brillouin zone of is generated on the mutually orthogonal mirror invariant z-direction, i.e., planes kx,y = 0 and kx,y = ±π when kz = ±π/2, where Z π all the nodal rings are reduced to straight lines. 1 ν = (ny∂kz nz − nz∂kz ny)dkz (18) 2π −π The topological features of nodal chains are character- ized by the quantized Berry phases. It is well known that where (ny, nz) = (hy, hz)/|h|. Using the same param- the Berry phase around a single nodal line is π [3]. Conse- eters as shown in Figs.2(e-h), we calculate the values quently, the Berry phase around the chain point (enclos- of the topological winding numbers as a function of kx ing two nodal lines) becomes 0 (=π + π). According to and ky in Figs.2(i-l). The results show that the re- bulk-edge correspondence, such nontrivial topology im- gions where the topological winding number values are plies the existence of surface states under open bound- nontrivial (ν = ±1) coincide with the regions where the ary conditions [3]. Let us take open boundary condition surface states appear. It is worth pointing out that the along the z direction and periodic boundary condition sign of the nontrivial winding number values make no along the x and y directions as an example to show the difference to the presence of surface states. In these re- generated surface states. The corresponding energy spec- gions, the corresponding Hamiltonian H(kx, ky) supports tra is calculated in Fig.2(e-h), where we can find that one-dimensional inversion-symmetric topological insula- 6 tors protected by a Z2 topological invariant. kz-resolved spin-polarization measurements in the kx-ky plane. In experiment, the time-of-flight imaging measures the IV. EXPERIMENTAL DETECTION OF NODAL kz-integrated spin polarizations CHAINS 1 Z 2π hgσsi(kx, ky) = hσs(kx, ky, kz)i dkz. (24) The shapes of nodal chains can be detected by measur- 2π 0 ing spin polarizations. The detection principle is outlined as follows. First, on one hand, we know that the nodal Interestingly, the recent ultracold atoms experiments in chains are generated by the connected nodal rings which [52, 53] have demonstrated that, if the topological Hamil- tonian on the BISs with respect to a kz plane (for ex- are determined by c ample, kz = kz) respects a magnetic group symme- BIS BIS c try, the kz-resolved spin polarization hσz(k , k , k )i hy(kx, ky, kz) = hz(kx, ky, kz) = 0. (19) x y z BIS BIS equals the integral of hσz(kx , ky , kz) over kz, i.e., On the other hand, when kz is fixed to kz0, the nodal BIS BIS BIS BIS c hgσzi(kx , ky ) = hσz(kx , ky , kz)i. The experiments chain Bloch Hamiltonian h(kx, ky, kz) is reduced to a also developed a virtual slicing imaging technique to scan two-dimensional Hamiltonian h(kx, ky, kz0). In this case, c kz and measure the spin polarizations varying with kz. In there are two types of band inversion surfaces (BISs) in this way, the kz-resolved spin polarizations on the BISs the kx-ky plane, named as BIS1 and BIS2. The BISs are can be effectively detected by measuring the correspond- defined based on the vanishing of the vector fields in the ing kz-integrated spin polarizations. In the following, we Bloch Hamiltonian [57]. Specifically, the BIS1 and BIS2 will prove that the spin polarizations in our system also are defined by have such feature. Firstly, we prove that our system also respects a mag- BIS1 BIS1 hy(kx , ky , kz0) = 0, (20) c netic group symmetry with respect to the kz = π/2 BIS2 BIS2 hz(kx , ky , kz0) = 0. (21) plane, which ensures that the kz-resolved spin polariza- tions hσzi on the BIS2 can be effectively detected by Basing on Eqs. (14, 20, 21), we can obtain the locations measuring the corresponding kz-integrated spin polariza- of BIS1 and BIS2. For kz0 6= 0, π, the BIS1 is formed tions. For kz = π/2, the BIS2 is a ring determined by BIS1 BIS2 BIS2 by four straight lines kx,y = 0, π; While for kz0 = 0, π, mz = 2tx cos(kx )+2ty cos(ky ). For kz = π/2+δkz, the BIS1 covers the whole kx-ky plane. The BIS2 is a the nodal chain Hamiltonian h(kx, ky, kz) on this ring BIS2 ring determined by mz − 2tz cos(kz0) = 2tx cos(kx ) + (BIS2) is written as BIS2 2ty cos(ky ). From Eqs. (19, 20, 21), we can find that, BIS2 BIS2 the gap-closing points for kz = kz0 determined by Eq. h(δkz) =8tso sin(kx ) sin(ky ) cos(δkz)σy (19) satisfy both Eq. (20) and Eq. (21), which means + 2tz sin(δkz)σz, (25) that the nodal points for kz = kz0 are exactly located at the intersection of the BIS and BIS . As a consequence, BIS2 BIS2 1 2 where h(δkz) = h(kx , ky , π/2+δkz). As we can see, by scanning kz and measuring the corresponding BIS1 the nodal chain Hamiltonian respects a magnetic group and BIS and their intersection, one can precisely map c 2 symmetry with respect to the kz = π/2 plane, i.e., out the shape of the entire nodal chain. −1 The BISs are measured by detecting the vanishing of Gh(δkz)G = h(−δkz), (26) the spin polarizations in equilibrium. Specifically, to measure the BIS1 and BIS2, the spin polarizations re- where G = σxK is the magnetic group symmetry opera- quired to be detected are tor. On the other hand, we have

hσsi = hψG(kx, ky, kz0)|σs|ψG(kx, ky, kz0)i, (22) h(±δkz)|ψG(±δkz)i = E−|ψG(±δkz)i, (27)

where |ψ (±δk )i = |ψ (kBIS2 , kBIS2 , π/2±δk )i are the where s = y, z, |ψG(kx, ky, kz0)i denotes the ground state G z G x y z ground states of h(±δk ), and E = E(k = π/2+δk ) = of h(kx, ky, kz0). The locations of BIS1 and BIS2 are re- z − z z spectively determined by the vanishing of the spin polar- E(kz = π/2 − δkz). Combining Eq. (26) and Eq. (27), izations [57] we can obtain

BIS1 BIS1 |ψG(δkz)i = G|ψG(−δkz)i. (28) hσy(kx , ky , kz0)i = 0,

BIS2 BIS2 hσz(kx , ky , kz0)i = 0. (23) By substituting Eq. (28) into Eq. (22), we can get

The reason is that, on the BIS1 (BIS2), |ψGi is re- hσz(δkz)i = hψG(δkz)|σz|ψG(δkz)i duced to the ground state of σz (σy), which exactly gives −1 = hψG(−δkz)|G σzG|ψG(−δkz)i hσyi = 0 (hσzi = 0). Therefore, for getting the BISs and = −hσ (−δk )i, (29) their intersection varying with kz, we need to perform z z 7

BIS2 BIS2 where hσz(±δkz)i = hσz(kx , ky , π/2 ± δkz)i. Then be applied to detect the shapes of the outer nodal chains we have illustrated in Fig.2(b-d).

BIS2 BIS2 BIS2 BIS2 c hgσzi(kx , ky ) = hσz(kx , ky , kz = π/2)i. (30) Next, we prove that the mirror symmetry with respect to the mirror invariant plane kz = 0 guarantees that the V. SUMMARY kz-resolved spin polarizations hσyi on the entire kx-ky plane (not just on the BIS ) have the similar feature. 1 In summary, we have presented an experimental For k = δk , according to Eq. (16), our system respects z z scheme to realize topological nodal chain semimetal a mirror symmetry with respect to the kc = 0 plane, i.e., z phases with ultracold atoms trapped in optical Raman 0 −1 0 lattices. In our scheme, we have constructed a three- Mh (δkz)M = h (−δkz), (31) dimensional optical Raman lattice that could produce 0 where h (δkz) = h(kx, ky, δkz). Similarly, we can prove next nearest-neighbor spin-orbit couplings and host var- that ious shapes of mirror symmetry protected nodal chains in its energy spectra. Finally, we have shown that the 0 0 |ψG(δkz)i = M|ψG(−δkz)i, shapes of the nodal chains could be detected by mea- suring spin polarizations. Thus, our study has demon- hσy(δkz)i = −hσy(−δkz)i, (32) strated the potential of optical Raman lattice systems as where hσy(±δkz)i = hσy(kx, ky, ±δkz)i. Then we get a versatile platform for exploring topological nodal-chain semimetal phases. c Moreover, our study sets an example for realizing com- hgσyi(kx, ky) = hσy(kx, ky, k = 0)i. (33) z plex spin-orbit couplings in optical Raman lattices by Therefore, basing on Eqs. (30, 33) and the virtual slicing applying suitable Raman potentials. This strategy is ex- imaging technique [52, 53], the kz-resolved spin polariza- pected to overcome the challenges faced in generating the tions on the BISs could be effectively measured. From complex spin-orbit couplings required for realizing high- this measurement, we can acquire the BISs and their in- dimensional topological phases in optical lattice systems. tersection and map out the shapes of nodal chains. A straightforward example is that the topological nodal As a concrete example, Fig.3 presents how the shape links and knots [19] could be realized by adding pertur- of the inner nodal chain (see Fig.2(a)) is detected bations into the Raman potentials applied in our nodal by measuring kz-resolved spin polarizations. The chain platform. In the future, it would be quite inter- corresponding spin polarizations hσy,zi versus kx and esting to apply this strategy in optical Raman lattices ky for different kz0 are numerically calculated in Figs. to construct the spin-orbit couplings for realizing the 3(a-f). For kz0 = 2π/3, the results in Fig.3(a) show four-band topological phases, including the Z2 topologi- that hσyi = 0 when kx,y = 0, π and hσzi = 0 when cal insulator phases [58, 59], topological Dirac semimetal kx = ky = 0. Then we can find that the BIS1 locates at phases [60, 61] and higher-order topological phases [62]. BIS1 BIS2 BIS2 kx,y = 0, π and the BIS2 at kx = ky = 0. Through the intersection of the BIS1 and BIS2, we can detect the nodal point as (kx, ky, kz)=(0, 0, 2π/3). For kz0 = π/2, VI. ACKNOWLEDGMENT the results in Fig.3(b) show that the BIS 1 still locates at BIS1 kx,y = 0, π, while the BIS2 locates on a ring determined B.B.W. and J.H.Z. contributed equally to this work. BIS2 BIS2 by cos(kx ) + cos(ky ) = 3/2. The intersection of This work was supported by the National Key Research the BIS1 and BIS2 allows us to detect the four nodal and Development Program of China (2017YFA0304203), points as (kx, ky, kz)=(0, ±π/3, π/2),(±π/3, 0, π/2). The National Natural Science Foundation of China (NSFC) results in Figs.3(c-f) also demonstrate this point, i.e., (12034012, 12074234), Research Program of the Ed- the nodal points associated with different kz0 can be ucation Department of Hubei province (D20181903), directly mapped out by measuring the kz-resolved spin Changjiang Scholars and Innovative Research Team in polarizations and extracting the intersection of the BIS1 University of Ministry of of China (PC- and BIS2. In this way, the shape of the inner nodal SIRT)(IRT 17R70), Fund for Shanxi 1331 Project Key chain is detected. This method is general and also can Subjects Construction, and 111 Project (D18001).

[1] C. Fang, H. Weng, X. Dai, and Z. Fang, Chin. Phys. B [3] A. A. Burkov, M. D. Hook, and L. Balents, Phys. Rev. 25, 117106 (2016). B 84, 235126 (2011). [2] G. E. Volovik, The Universe in a Helium Droplet (Oxford [4] H. Weng, Y. Liang, Q. Xu, R. Yu, Z. Fang, X. Dai, and University Press, 2003). Y. Kawazoe, Phys. Rev. B 92, 045108 (2015). 8

[5] Y. Kim, B. J. Wieder, C. L. Kane, and A. M. Rappe, [33] P. Wang, Z.-Q. Yu, Z. Fu, J. Miao, L. Huang, S. Chai, H. Phys. Rev. Lett. 115, 036806 (2015). Zhai, and J. Zhang, Phys. Rev. Lett. 109, 095301 (2012). [6] R. Yu, H. Weng, Z. Fang, X. Dai, and X. Hu, Phys. Rev. [34] V. Galitski and I. B. Spielman, Nature (London) 494, 49 Lett. 115, 036807 (2015). (2013). [7] T. T. Heikkila and G. E. Volovik, New J. Phys. 17, [35] J. Zhang, H. Hu, X.-J. Liu, H. Pu, Annual Re- 093019 (2015). view of Cold Atoms and Molecules, pp. 81-143 (2014); [8] C. Fang, Y. Chen, H.-Y. Kee, and L. Fu, Phys. Rev. B arXiv:1411.3043v1. 92, 081201(R) (2015). [36] N. Goldman, G. Juzeli¯unas,P. Ohberg,¨ and I. B. Spiel- [9] G.E. Volovik, JETP Lett. 107, 324–326 (2018). man, Rep. Prog. Phys. 77, 126401 (2014). [10] V.B. Eltsov, T. Kamppinen, J. Rysti, G.E. Volovik, [37] H. Zhai, Reports on Progress in Physics 78, 026001 arXiv:1908.01645. (2015). [11] J. Nissinen, Tero T. Heikkila, G.E. Volovik, [38] L. Zhang and X.-J. Liu, Synthetic Spin-Orbit Coupling arXiv:2008.02158 in Cold Atoms, pp. 1-87 (2018) (World Scientific, Singa- [12] T. Bzduˇsek,Q.S. Wu, A. R¨uegg,M. Sigrist, and A. A. pore); arXiv:1806.05628v1. Soluyanov, Nature (London) 538, 75 (2016). [39] M. Atala, M. Aidelsburger, J. T. Barreiro, D. Abanin, [13] R. Yu, Q. Wu, Z. Fang, and H. Weng, Phys. Rev. Lett. T. Kitagawa, E. Demler and I. Bloch, Nat. Phys. 9, 795 119, 036401 (2017). (2013). [14] G. Chang, S.-Y. Xu, X. Zhou, S.-M. Huang, B. Singh, B. [40] G. Jotzu, M. Messer, R. Desbuquois, M. Lebrat, T. Wang, I. Belopolski, J. Yin, S. Zhang, A. Bansil, H. Lin, Uehlinger, D. Greif and T. Esslinger, Nature 515, 237 and M. Z. Hasan, Phys. Rev. Lett. 119, 156401 (2017). (2014). [15] S. S. Wang, Y. Liu, Z. M. Y. Yu, X. L. Sheng, S. A. Yang, [41] M. Aidelsburger, M. Lohse, C. Schweizer, M. Atala, J. Nat. Commun. 8, 1844 (2017). T. Barreiro, S. Nascimb`ene,N. R. Cooper, I. Bloch and [16] Q. Yan, R. Liu, Z. Yan, B. Liu, H. Chen, Z. Wang, and N. Goldman, Nat. Phys. 11, 162 (2014). L. Lu, Nat. Phys. 14, 461 (2018). [42] Z. Wu, L. Zhang, W. Sun, X.-T. Xu, B.-Z. Wang, S.-C. [17] C. J. Yi, B. Q. Lv, Q. S. Wu, B. B. Fu, X. Gao, M. Yang, Ji, Y. Deng, S. Chen, X.-J. Liu, J.-W. Pan, Science 354, X. L. Peng, M. Li, Y. B. Huang, P. Richard, M. Shi, G. 83 (2016). Li, O. V. Yazyev, Y. G. Shi, T. Qian, H. Ding, Phys. [43] B. Song, L. Zhang, C. He, T. F. J. Poon, E. Hajiyev, Rev. B 97, 201107 (2018). S. Zhang, X.-J. Liu and G.-B. Jo, Science advances 4, [18] M. Xiao, X. Q Sun, and S. Fan, Phys. Rev. B 99, 094206 eaao4748 (2018). (2019). [44] Y. Xu, Front. Phys. 14, 43402 (2019). [19] Z. Yang, C.-K. Chiu, C. Fang, and J. Hu, Phys. Rev. [45] Y. Xu and C. Zhang, Phys. Rev. A 93, 063606 (2016). Lett. 124, 186402 (2020). [46] D.-W. Zhang, Y. X. Zhao, R.-B. Liu, Z.-Y. Xue, S.-L. [20] W. Chen, H.-Z. Lu, and J.-M. Hou, Phys. Rev. B 96, Zhu, and Z. D. Wang, Phys. Rev. A 93, 043617 (2016). 041102(R) (2017). [47] C. Yin, L. Li, and S. Chen, Phys. Rev. A 97, 013604 [21] Z. Yan, R. Bi, H. Shen, L. Lu, S.-C. Zhang, and Z. Wang, (2018). Phys. Rev. B 96, 041103(R) (2017). [48] B. Liu, P. Zhang, R. Zhang, H. Gao, and F. Li, Phys. [22] P.-Y. Chang and C.-H. Yee, Phys. Rev. B 96, 081114(R) Rev. A 98, 063610 (2018). (2017). [49] W.-Y. He, D.-H. Xu, B. T. Zhou, Q. Zhou, and K. T. [23] Y. Zhou, F. Xiong, X. Wan, and J. An, Phys. Rev. B 97, Law, Phys. Rev. A 97, 043618 (2018). 155140 (2018). [50] T. Zhou, Z. Yu, Z. Li, X. Chen, X. Zhou, Phys. Rev. A [24] L. Li, C. H. Lee, and J. Gong, Phys. Rev. Lett. 121, 102, 023328 (2020). 036401 (2018). [51] Y.-H. Lu, B.-Z. Wang, X.-J. Liu, Science Bull. 65, 2080 [25] R. Bi, Z. Yan, L. Lu, and Z. Wang, Phys. Rev. B 96, (2020). 201305(R) (2017). [52] B. Song, C. He, S. Niu, L. Zhang, Z. Ren, X.-J. Liu, [26] M. Ezawa, Phys. Rev. B 96, 041202 (2017). G.-B. Jo, Nature Phys. 15, 911 (2019). [27] G. Bian, T.-R. Chang, R. Sankar, S.-Y. Xu, H. Zheng, [53] Z.-Y. Wang, X.-C. Cheng, B.-Z. Wang, J.-Y. Zhang, Y.- T. Neupert, C.-K. Chiu, S.-M. Huang, G. Chang, I. Be- H. Lu, C.-R. Yi, S. Niu, Y. Deng, X.-J. Liu, S. Chen, lopolski, D. S. Sanchez, M. Neupane, N. Alidoust, C. Liu, J.-W. Pan, arXiv:2004.02413v1. B. K. Wang, C.-C. Lee, H.-T. Jeng, C. Zhang, Z. Yuan, [54] X.-J. Liu, Z.-X. Liu, and M. Cheng, Phys. Rev. Lett. S. Jia, A. Bansil, F. Chou, H. Lin and M. Z. Hasan, Nat. 110, 076401 (2013). Commun. 7, 10556 (2016). [55] X.-J. Liu, K. T. Law, and T. K. Ng, Phys. Rev. Lett. [28] M. Neupane, I. Belopolski, M. M. Hosen, D. S. Sanchez, 112, 086401 (2014). R. Sankar, M. Szlawska, S.-Y. Xu, K. Dimitri, N. Dhakal, [56] B.-Z. Wang,Y.-H. Lu, W. Sun, S. Chen, Y. Deng and P. Maldonado, P. M. Oppeneer, D. Kaczorowski, F. X.-J. Liu, Phys. Rev. A 97, 011605 (R) (2018). Chou, M. Z. Hasan and T. Durakiewicz, Phys. Rev. B [57] L. Zhang, L. Zhang, S. Niu, and X.-J. Liu, Science Bull. 93, 201104(R) (2016). 63, 1835 (2018). [29] N. Goldman, J. C. Budich, and P. Zoller, Nat. Phys. 12, [58] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 639 (2016). 106803 (2007). [30] D.-W. Zhang, Y.-Q. Zhu, Y. X. Zhao, H. Yan, and S.-L. [59] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). Zhu, Adv. Phys. 67, 253 (2019). [60] S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. [31] N. R. Cooper, J. Dalibard, and I. B. Spielman, Rev. Mod. Mele, and A. M. Rappe, Phys. Rev. Lett. 108, 140405 Phys. 91, 015005 (2019). (2012). [32] Y.-J. Lin, K. Jim´enez-Garc´ıa,and I. B. Spielman, Nature [61] S. M. Young and C. L. Kane, Phys. Rev. Lett. 115, (London) 471, 83 (2011). 126803 (2015). 9

[62] F. Schindler, A. M. Cook, M. G. Vergniory, Z. Wang, 4, eaat0346 (2018). S.S.P. Parkin, B. A. Bernevig, and T. Neupert, Sci. Adv.