Quantum transport in topological nodal-line semimetals

ABSTRACT Topological nodal-line semimetals offer an attractive research platform for exploring a variety of novel phenomena, which have attracted great research interest in the past decade. There are three unique features of the nodal-line semimetals: (i) band crossing along the closed loop that carries Berry phase, (ii) torus-shaped Fermi surface as the Fermi energy deviates from the nodal loop, and (iii) drumhead-like surface states induced by the band topology, which comprise the main scenarios for most of its electronic properties. Here, we review recent progress on the quantum transport theory of nodal-line semimetals. We show that the properties (i) and (ii) result in a fascinating interplay between the effective dimensionality of electron diffusion and the band topology, which may give rise to 3D weak localization or 2D weak anti-localization effect that can be probed by the magnetoconductivity induced by a weak magnetic field. Moreover, (i) and (ii) can also be revealed by the pattern of the Shubnikov–de Hass oscillation in a strong magnetic field. For the massive nodal-line semimetals with spin-orbit coupling, we show that the magnetic field can induce a topological Lifshitz transition of the Fermi torus from genus 1 to 0. Interestingly, the Lifshitz transition is in good coincidence with the sign reversal of the magnetoresistivity. Compared to the conventional Fermi sphere, the Fermi torus possesses multiple reflection channels, which may give rise to novel scattering processes such as anomalous Andreev reflection at the interface between the semimetal and a superconductor. Besides the bulk states, the drumhead surface states also possess interesting properties. It is shown that the spin polarization of the surface states can induce resonant spin-flipped reflection, which can be detected through both spin and charge transport measurements. We also briefly review the recent experimental progress on the transport properties of nodal-line semimetals and compare the results with the theoretical predictions. Graphical Abstract


I. Introduction
In the past two decades, the topological phase of matter has become the frontier of condensed matter physics [1][2][3]. In stark contrast to the conventional Landau paradigm of phase transition in terms of symmetry breaking, different topological phases are characterized by the topological properties of the Bloch bands. According to their gapped and gapless nature, the topological phases can be roughly divided into insulating and semimetal phases, respectively [1][2][3][4][5]. A variety of topological semimetals have been theoretically predicted and experimentally confirmed, involving Weyl [6,7], Dirac [8], and nodal-line semimetals [9], etc., which correspond to different types of band crossing and the associated band topologies [5][6][7][8][9][10][11][12][13][14][15]. Specifically, the adjacent bands meet linearly at isolated points in the Brillouin zone for the Weyl and Dirac semimetals while intersect along closed loops for the nodal-line semimetals. Interestingly, the one-dimensional nodal loops can be arranged into various configurations, which yield different topological phases and rich physical effects [16][17][18]. Exotic surface states generally emerge at the open boundaries of topological semimetals due to the bulk-boundary correspondence. For example, the helicoid surface states in the Weyl and Dirac semimetals result in disconnected surface Fermi arcs, which connect the projections of the bulk nodal points in the surface Brillouin zone [7,8]. For the nodal-line semimetals, the so-called drumhead-like surface states (DSS) are induced by the π Berry phase carried by the nodal loop, which occupy the area enclosed by the nodal loop projected onto the surface Brillouin zone [19,20]. A common feature of topological semimetals is that both bulk and surface states can participate in various physical processes, making it an attractive research platform for exploring their sole and combined effects. This differs from the situation of topological insulators, in which only the topological surface states are usually of interest because the bulk bands stay far away from the Fermi energy. Various material candidates have been theoretically proposed or experimentally confirmed as nodal-line semimetals, such as HgCr 2 Se 4 [21], graphene networks [22], Cu 3 (Pd/Zn)N [23,24], SrIrO 3 [25,26], TiTaSe 2 [27], Ca 3 P 2 [19,28], CaTe [29], compressed black phosphorus [30], CaAg(P/As) [31], CaP 3 family [32], PdS monolayer [33], Zintl compounds [34], BaM X 3 (M = V,Nb, and Ta,X = S,Se) [35], rare earth monopnictides [36], alkaline-earth compounds [37][38][39], other carbon-based materials [40], and metallic rutile oxides XO 2 ,(X = Ir,Os,Rd) [41], PbTaSe 2 [42], ZrXY(X = Si,Ge,Y = S,Se,Te) [17,[43][44][45][46][47][48], and PtSn 4 [49]. The experimental approaches widely used for probing the electronic properties of topological materials involve the angle-resolved photoemission spectroscopy (ARPES), scanning tunneling microscopy (STM) and transport measurements. The ARPES directly probes the band structures of the bulk and surface states. The information of the Fermi surface can also be revealed by the STM through the quasiparticle interference process. Compared with ARPES and STM, transport measurements possess their own advantages. First, by elaborate design and fabrication of electronic devices and proper field manipulations, the message of the bulk and surface states can be selectively extracted; second, the transport techniques are especially pertinent to potential applications of topological materials in functional electronic devices. Great advances have been made in the transport studies on both Weyl/ Dirac and nodal-line semimetals. A special review for the former has been given in Ref [50], while that for the latter is still absent.
In this paper, we give a concise review of the recent progress in quantum transport theories of nodal-line semimetals. The rest of this paper is organized as follows. In Section II, we introduce the main properties of nodalline semimetals, which include the specific band topology, surface states and the peculiar torus-shaped Fermi surface. Several novel effects stemming from these properties are introduced in the following sections, which involve the weak (anti-)localization effect in Section III, quantum oscillation in Section IV, topological Lifshitz transition of the Fermi surface in Section V, anomalous Andreev reflection in Section VI, and the resonant spin-flipped reflection in Section VII. In Section VIII, we also review the recent experimental progress and compare it with theories. Finally, we draw our conclusion and give an outlook in Section IX.

II. Main properties of nodal-line semimetals
The terminology of nodal-line semimetal indicates a particular type of band crossing along one-dimensional lines or loops [9]. For most of the candidates, the Hamiltonian of the system can be linearized around the nodal-line in the reciprocal space. One can understand each point of the nodal-line as a single Dirac cone in the plane locally perpendicular to the nodal-line. Therefore, the nodal-line semimetal can be understood as a collection of two-dimensional massless Dirac Fermion along the line. As a result, the band topology can be interpreted by the Berry phase carried by the nodalline as shown in Figure 1), which is similar to that of the Dirac cone in graphene and surface states of the topological insulator. Besides the similarity to the two-dimensional Dirac fermion, the three-dimensional nodalline semimetal possesses its own features, such as the drumhead-like surface states and the torus-shaped Fermi surface when doped. For real materials, there may exist one or more nodal-lines. In most of cases, several isolated nodal-lines/loops can be analyzed separately so that the effective model with a single nodal-line is sufficient to capture the main physics therein.
To be specific, we adopt the following Hamiltonian to describe the nodalline semimetal as [9] HðkÞ where λ is the model parameter, v is the velocity along the z-direction, k 2 ¼ k 2 x þ k 2 y þ k 2 z is the total momentum squared and the Pauli matrices σ x;z operate in the spin or pseudo-spin space that will be specified if necessary. The Hamiltonian (1) yields the eigenvalues ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi v 2 k 2 z þ λ 2 ðk 2 0 À k 2 Þ 2 q and the corresponding eigenstates u � ðkÞ j i ¼ ð� hvk z = ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi The band degeneracy occurs as k 2 x þ k 2 y ¼ k 2 0 and k z ¼ 0, which defines a nodal loop in the momentum space, see Figure 1). As the Fermi energy E F deviates from zero that usually occurs for real materials, the nodal loop grows into a torus-shaped Fermi surface; see Figure 1). For E F � � hλk 2 0 , it is convenient to linearize and parameterize the Hamiltonian (1) into HðδkÞ ¼ � hv 0 κðcos φσ z þ sin φσ x Þ; (2) through the substitution k x ¼ ðk 0 þ κ cos φÞ cos θ; k y ¼ ðk 0 þ κ cos φÞ sin θ; k z ¼ κ=α, in which v 0 ¼ 2λk 0 and α ¼ v=v 0 is the ratio between the velocities along the z direction and in the x-y plane. The small radius κ, the toroidal angle θ and the poloidal angle φ are labeled in Figure 2). The small momentum δk ¼ κðcos φ; sin φ=αÞ is defined in a local poloidal plane in momentum space denoted by the toroidal angle θ.
To understand the band topology, let us consider the transverse wavevector k jj ¼ ðk x ; k y Þ as a parameter, and then the Hamiltonian (1) describes an effective 1D system in the z-direction. If k jj lies inside the Figure 2. (a) Torus-shaped Fermi surface, with major radius k 0 , minor radius κ, toroidal angle θ, and poloidal angle φ. (b) The impurity potentials results in 3D and 2D diffusion behaviors, respectively. (c) A coherent backscattering process from wave vector k to À k around the toroidal direction, which is a 3D diffusion behavior. (d) Backscattering from wave vector δk to À δk around poloidal direction, which is a 2D diffusion behavior [54]. nodal loop, i.e. jk jj j < k 0 , the effective 1D system is insulating with an energy gap Δðk jj Þ ¼ � hvk 0 0 opening around k z ¼ k 0 0 ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi k 2 0 À jk jj j 2 q .
One can verify that a nontrivial Berry phase can be obtained as [19,20]. Otherwise, the effective gap is opened at k z ¼ 0 and one has γ B ¼ 0 for jk jj j > k 0 . This result together with the fact that the Berry curvature vanishes anywhere away from the nodal loop indicates that the π-Berry phase is exactly carried by the nodal-line and the wave function will pick up the Berry phase along any closed path linked with the nodal-line. In the presence of chiral symmetry, such nontrivial topological winding results in zero-energy surface states encircled by the projection of the nodal loop onto the surface Brillouin zone [51,52], dubbed DSS [22]; see Figure 1). The above features of nodal-line semimetal can give rise to several unique transport properties that will be introduced in the following sections.

III. Weak localization and anti-localization effects in nodal-line semimetals
An electronic system subject to weak disorder and at sufficiently low temperatures resides in the quantum diffusive regime, in which the mean free path is much less than the phase coherence length. The quantum interference effect results in additional quantum corrections to the conductivity [53]. It has been known that the quantum correction to the electronic conductivity exhibits a universal dependence on the dimensionality and band topology of the system. The nontrivial properties of nodal-line semimetal enrich the conventional results of quantum diffusion [54]. First, the torus-shaped Fermi surface brings additional backscattering channels compared with the conventional Fermi sphere, which strongly modifies the quantum diffusion behavior; see Figure 2. Specifically, the effective dimension of diffusion will be two and three for the long-range (LR) and shortrange (SR) disordered potentials, respectively. Second, whether the electron can feel the π-Berry phase carried by the nodal line is determined by its diffusion dimension, thus exhibiting a fascinating interplay between the diffusion dimensionality and the band topology. We start with the Hamiltonian (1) and (2) and make the substitutions σ z ! σ x ; σ x ! σ y for a convenient gauge choice, which yield H 0 ðkÞ ¼ � hvk z σ y þ � hλðk 2 À k 2 0 Þσ x ; H 0 ðδkÞ ¼ � hv 0 κðcos φσ x þ sin φσ y Þ: The Fermi energy E F is assumed to intersect the conduction band that is the only band concerned below. For H 0 , the dispersion of the conduction band is � k ¼ � hv 0 κ and the corresponding wave function is ψ k ðrÞ ¼ ½1; e iφ � T e ik�r = ffi ffi ffi ffi ffi ffi 2V p , with V the volume of the system. The density of states at the Fermi energy is ρ 0 ¼ KE F =ðαh 2 v 2 0 Þ, proportional to the circumference K ¼ 2πk 0 of the nodal loop.
The Hamiltonian H 0 ðkÞ and H 0 satisfy two types of anti-unitary symmetries as where T 1 ¼ K and T 2 ¼ iσ y K with K being the complex conjugation. Here, T 1 is a spinless time-reversal symmetry, while T 2 is a local antiunitary symmetry involving a spin/pseudospin inversion. The second equation of Equation (4) can be regarded as the spinful time-reversal symmetry for the 2D subsystem in the plane labeled by θ [see Figure 2)]. According to the symmetry classification, T 1 and T 2 correspond to the orthogonal and symplectic classes, respectively [55,56]. As a result, they will lead to weak localization (WL) and weak anti-localization (WAL), respectively, depending on which physical process in Figure 2 dominates the backscattering. There are two backscattering channels for the torus-shaped Fermi surface; see Figure 2. Which type of backscattering dominates the WL or WAL effects is determined by details of the disorder potential and size of the nodal loop. In this section, we consider two limiting cases of the disorder potentials with either SR or LR type relative to the size of the nodal loop. The SR limit is also referred to as the so-called white-noise disorder, which causes constant scattering strength in the reciprocal space so that the backscattering is dominated by the interference along the toroidal direction [see Figure  2)]. Such an interference loop does not encircle the nodal-line so that electrons can not feel the π Berry phase. The LR disorder may exist due to unconventional screening effects [57]. In this case, the backscattering is dominated by the interference along the poloidal direction locally in the reciprocal space. The interference loop encircles the nodal line [see Figure  2)] so that the π Berry phase plays an important role in the quantum diffusion.
The interference effect is manifested as the quantum correction to the conductivity, which can be obtained by calculating the Feynman diagrams in Figure 3). Three leading-order diagrams are involved in the calculation, namely, the bare Hikami box [56] and two dressed Hikami boxes [58,59]. The latter two jointly contribute −1/2 to the correction of the former for both the SR and LR limits, the same as that in graphene [58]. The overall quantum correction to the zero-temperature conductivity in the z-direction in both the SR and LR limits is given by where , e ¼ v 0 τ e is the mean free path and , ϕ is the phase coherence length. In the SR limit, a 3D scaling is obtained for the quantum correction of the conductivity and the overall minus sign means the WL effect [53]. In contrast, the LR correction is positive that exhibits a WAL behavior with a 2D scaling law. The prefactor K in σ L z ensures the unit of a 3D conductivity and also indicates an enhancement of the WAL effect contributed by a large number of 2D diffusion planes; see Figure 2). Because of the different dimensionality of diffusion, σ S z will saturate as , ϕ ! 1 while σ L z always increases as lnð, ϕ =, e Þ. The quantum correction to the conductivity can be suppressed by a magnetic field-induced dephasing effect [60]. It allows the probe of WL and WAL corrections using magnetoconductivity measurements. The magnetoconductivities of the nodal-line semimetal for the two types of disorder are defined by the difference between the conductivity with and without a magnetic field, i.e. δσ S;L z ðBÞ ; σ S;L z ðBÞ À σ S;L z ð0Þ, with The numerical results are shown in Figure 4. Specifically, in the SR regime, the WL effect leads to a positive magnetoconductivity [ Figure 4)]; In the LR regime, a negative magnetoconductivity corresponds to the WAL effect [ Figure 4)]. Moreover, the magnetoconductivities exhibit different Bdependence for the two regimes. Specifically, δσ S z ðBÞ / ffi ffi ffi B p in the SR regime while δσ L z ðBÞ / À ln B in the LR regime; see Figure 4. It is notable that δσ L z ðBÞ in the SR regime is three orders larger than δσ S z ðBÞ. This is because the system is composed of a large number of 2D planes of diffusion in the SR regime, as discussed above.
From the discussion, one can see that quantum diffusion in the nodal-line semimetal can lead to intriguing phenomena such as the 3D magnetoconductivity with a 2D scaling law and an extraordinarily large value. Relevant effects have been reported in a variety of nodal-line semimetals [61,62], which may be attributed to the scenario discussed in this section; see section VIII.

VI. Quantum oscillations in nodal-line semimetals
The shape of the Fermi surface as well as the geometry phase encoded in the wave functions can be probed through the quantum oscillation method [63,64], or specifically, Shubnikov-de Haas oscillation of resistivity that will be discussed here. The shape and size of the Fermi surface and the Berry phase information can be extracted through the periodicity and phase shift of the oscillation pattern, respectively [65]. The coexistence of the torus-shaped Fermi surface and the π Berry flux in nodal-line semimetal indicates that the fascinating interplay between the two may lead to interesting results in the quantum oscillation measurements [66][67][68]. A series of experiments have reported the results of quantum oscillations of the resistivity ρ  in nodal-line semimetals subjected to a strong magnetic field. Several theories have been developed on this problem and we will review the results in Ref [66]. here.
The main results of the Shubnikov-de Haas oscillation of resistivity are reflected by the Lifshitz-Kosevich formula [64] as where F is the oscillation frequency, ϕ is the phase shift and B is the magnetic field. The information of the Fermi surface and the geometric phase is revealed by F and ϕ. In experiments, both F and ϕ can be extracted from the Landau fan diagram. Specifically, the condition for the n th peak of the resistivity is n ¼ F=B þ ϕ. Therefore, F is the slope of the linear function and the intercept on the n axis corresponds to the phase shift ϕ.
According to the Onsager relation [92] F ¼ ð� h=2πeÞA, the oscillation frequency F is proportional to the extremal cross-section area A of the Fermi surface perpendicular to the magnetic field. By changing the directions of the external magnetic field, the frequency F varies and the shape and size of the Fermi surface can then be extracted. For nodal-line semimetal with light doping, the Fermi surface possesses a torus shape that is highly anisotropic, which is in stark contrast to the spherical Fermi surface for conventional electrons and Weyl fermions. For a magnetic field normal to the nodal-line plane, the extremal cross sections are determined by @A=@k z j E F ¼ 0, which yields two extremal cross sections labeled by α and β lying in the k z ¼ 0 plane; see Figure 5). The energy dispersions are E � ¼ � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi for electron (+) and hole ( À ), respectively. By combining the dispersion and Onsager relation, the high and low frequencies are found to be F α ¼ ð� hλk 2 0 þ E F Þ=2λe for the outer circle and F β ¼ ð� hλk 2 0 À E F Þ=2λe for the inner circle, respectively [66]. For the magnetic field in the plane of the nodal-line, the extremal cross sections labeled by γ and δ can be obtained similarly; see Figure 5.
The phase shift for each frequency component of the quantum oscillation satisfies the relation as where ϕ B is the Berry phase [93,94] and ϕ 3D ¼ �1=8 is a correction pertinent to three-dimensional systems. The rule for the " � " sign is as follows: for a maximum cross section, ϕ 3D ¼ À 1=8 for electrons and 1=8 for holes; for a minimum cross section, the opposite sign should be taken. For a conventional Fermi sphere, there is only a maximum cross section so that ϕ 3D ¼ À 1=8 for electrons and 1=8 for holes. The situation becomes more complex for the nodal-line semimetal with a torus-shaped Fermi surface [66]. The Fermi torus possesses both maximum and minimum cross sections; see Figure 5) for the magnetic field normal to the nodalline plane and Figure 5) for the in-plane magnetic field. Moreover, the Berry phase is determined by the fact whether the contours encircle the nodal line or not. For the out-of-plane magnetic field, the contours denoted by α and β do not encircle the nodal line and thus ϕ B ¼ 0 for both cross sections. For the in-plane magnetic field, it is easy to see ϕ B ¼ 0 for the contour γ while ϕ B ¼ π for δ. Combining the arguments above, one can obtain the overall phase shift for various cross sections as summarized in Table 1. For example, as the magnetic field is normal to the nodal-line plane, the cross sections are α and β and ϕ B ¼ 0 for both. The phase shift for the maximum cross section α is calculated as ϕ ¼ À 1=2 þ 0 À 1=8 ¼ À 5=8 for electrons and ϕ ¼ À 1=2 þ 0 þ 1=8 ¼ À 3=8 $ 5=8 for holes. The phase shift for the minimum cross section β is calculated as ϕ ¼ À 1=2 þ 0 þ 1=8 ¼ À 3=8 $ 5=8 for electrons and ϕ ¼ À 1=2 þ 0 À 1=8 ¼ À 5=8 for holes. The phase shift for the cross sections γ and δ can be obtained in a similar way.
Although the theory is based on the minimal model (1) of the nodal-line semimetal, the scenario for the phase shift should hold generally, which can be compared with experimental results of real materials with different configurations of nodal loops, such as Cu 3 PdN and ZrSiS; see the discussion in section VIII.

V. Magnetic field induced Lifshitz transition
Transport measurements are a direct manifestation of the properties of the Fermi surface. Therefore, modulation of the Fermi surface by pressure, doping or external fields can be well reflected by the transport signatures. Interestingly, apart from continuous deformations, the Fermi surface can also undergo abrupt changes in its topology, the so-called topological Lifshitz transition [95]. Such a transition may lead to the discontinuity of physical quantities. The torus-shaped Fermi surface of the nodal-line semimetal possesses genus one and thus provides an ideal platform to investigate its Lifshitz transition. Moreover, the interplay between the topology of the Fermi surface and that of the wave functions may yield unique physical effects.
It is convenient to change the Fermi surface by an external field. For this purpose, we consider a massive nodal-line semimetal with finite spin-orbit coupling which exists in many material candidates. We revise the linearized Hamiltonian H 0 ðδkÞ in Equation (3) by involving a mass term induced by the spin-orbit coupling, which yields [96] H mass ¼ H 0 þ Δτ z σ z (9)  with the sign of the mass term depending on the spin τ z while σ z is the pseudo-spin. The energy spectra become ε � k κ ð Þ ¼ � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi � hv 0 κ ð Þ 2 þ Δ 2 q in which a gap 2Δ is opened along the nodal loop. The gap induces a finite spreading of the Berry curvature in the reciprocal space that introduces the self-rotation and thus a finite orbit magnetic moment (OMM) of the wave packet. Then the coupling between the OMM and a magnetic field can be implemented [94,98]. We focus on the regime that the Zeeman coupling between the magnetic field and the electronic spin is weak so that the main results do not rely on the spin or equivalently, the sign of the mass. We set τ z ¼ 1 henceforth.
Without loss of generality, we consider a positive Fermi energy, such that only the conduction band contributes to the transport and the wave function is found to be ju k i ¼ ðκe À iφ = ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ; γ= ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi κ 2 þ γ 2 p Þ with γ ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi m 2 þ κ 2 p À m and m ¼ Δ=ð� hv 0 Þ. The Berry curvature defined by Ω k ¼ ihÑ k u k j � jÑ k u k i has only a nonzero value along the toroidal direction in the parametric coordinates as Ω k ¼ À m=2ðκ 2 þ m 2 Þ 3 2ê θ which reduces to a Berry flux density as m ! 0. In the present twoband system, the OMM is related to the Berry curvature by , which is also a vector field along the toroidal direction. One can infer that when an external magnetic field is imposed, only its components in the x-y plane couple to the OMM. We consider a magnetic field B ¼ Bx in the x-direction, which gives rise to the coupling energy as δP k ðκ; θÞ ¼ À M k � B ¼ À ðmev 0 BÞ sin θ=½2ðκ 2 þ m 2 Þ�: The whole energy of the conduction band becomes � k ðκ; θÞ ¼ � þ k ðκÞ þ δ k ðκ; θÞ, which now depends on the toroidal angle θ. A direct consequence of this anisotropy is a deformation of the Fermi torus, as shown in Figure 6. We will demonstrate in the following that such a deformation plays a pivotal role in the magnetotransport.
In the weak field limit B � � 2 F =ð� hev 2 0 Þ the impact of Landau level can be neglected. By using the semiclassical Boltzmann formalism, the longitudinal conductivity can be obtained as [97,99] Analytical results can be obtained in the κ F � m regime as σ x where κ F is the Fermi vector, σ 0 is the Drude conductivity and δσ is the additional magnetic conductivity induced by the OMM. Different from the situation of Weyl semimetal [100] or topological insulator [99], the magnetoconductivity in the massive nodal-line semimetal cannot be induced by the pure Berry curvature effect without taking into account the OMM effect [97]. That is to say, the magnetoconductivity is exactly caused by the deformation of the Fermi torus [see Figure 6] due to the coupling between the OMM and the magnetic field.
As the coupling between the OMM and the magnetic field becomes stronger, which can be achieved by increasing either B or the ration m=κ, the deformation of the Fermi torus becomes more severe. Surprisingly, the  magnetoresistivity (MR) does not obey a monotonic dependence on B; see Figure 7, which first goes up to its maximum and then drops rapidly as B increases. As B increases further, the MR can even change its sign. The sign reversal of MR induced by the Berry curvature has not been found in other systems [99,100], which indicates a possible new scenario in the massive nodal-line semimetal. The unique response of MR with respect to B in nodal-line semimetals is deeply rooted in the fact that the torus-shaped Fermi surface can undergo a topological Lifshitz transition from genus one to zero; see Figure 6, which has no counterpart in a conventional Fermi sphere. The topological Lifshitz transition provides us with a possible explanation of the sign reversal of MR in massive nodal-line semimetal, which is proved to be correct.
In the limit κ F � m, the Lifshitz transition can be induced by a small magnetic field. The numerical results show that the sign reversal points of MR are in good coincidence with those of the Lifshitz transition (open circles); See Figure 7. Such correspondence can also be analytical proved by expanding the energy as � F ¼ � hv 0 ðm þ ηÞ given that η � m. The critical field of Lifshitz transition can be obtained as B � ¼ 2� hmη=e. By expanding Equation (11) to the linear order of η, we obtain σðB � Þ ¼ σð0Þ, which means that the correspondence holds strictly in this regime [97].

VI. Anomalous Andreev reflection
One important feature of the torus-shaped Fermi surface is the existence of one additional backscattering channel compared with the conventional Fermi sphere, which has nontrivial effects on the quantum diffusion as discussed in Section III. In the ballistic regime, the torus-shaped Fermi surface can be revealed by the reflection process. Given that the current due to normal reflection compensates with the injected current, it is beneficial to search for other types of reflection processes. Andreev reflection (AR) is referred to as the electron-to-hole reflection at the interface between a metal and a superconductor, by which a pair of electrons penetrate into the superconductor and form a Cooper pair. Therefore, the AR process gives rise to a measurable current.
There exist two types of AR previously. The conventional retro-AR (RAR) indicates that all velocity components of the reflected hole are inverted so that the reflected hole retraces the path of the incident electron; see Figure 8). The specular AR (SAR) occurs if the incident electron and the reflected hole are from the conduction and valence bands, respectively [102]. During SAR only the velocity perpendicular to the interface is inverted; see Figure 8, which resembles the normal reflection apart from the electron-hole conversion. The most promising candidates for the implementation of SAR are the massless Dirac Fermion in graphene and topological semimetals [102][103][104], in which no energy gap needs to be overcome during the SAR progress.
According to the sign inversion of the velocities during the AR processes, it can be inferred that the RAR and SAR constitute a complete set of AR in two dimensions, in which only two velocity components exist. Interestingly, in three dimensions, the types of AR can be further enriched. A novel type of AR dubbed 'anomalous-trajectory Andreev reflection' (AAR) was proposed in Ref [101]., in which only one of the two lateral velocities is inverted; see Figure 8). Three types of AR processes are compared by their sign inversion during the AR in Table 2. Note that the type of AR is determined by the band structure of the normal metal or semimetal rather than that of the superconductor, because AR takes place on the normal side of the junction. As a result, it is sufficient to obtain the AAR by adopting a conventional superconductor.   [101].
Consider a nodal-line semimetal-superconductor junction oriented along x-direction with the toroidal circle of the Fermi torus perpendicular to the interface of the junction as shown in Figure 8. A spindegenerate nodal-line semimetal lies in the x < 0 region and a conventional s-wave superconductor occupies the x > 0 region. To study the problem of electron-hole scattering, it is convenient to write the Hamiltonian in the Nambu space. The Hamiltonian of the whole junction can be expressed as where ΘðxÞ is the unit step function and τ x;z are the Pauli matrices in the Nambu space, H NL and H S correspond to the nodal-line semimetal and the conventional s-wave superconductor, respectively, H 0 ðkÞ captures a spin-degenerate nodal-line semimetal and μ is its chemical potential, and UδðxÞτ z is a δ-function interface barrier at x ¼ 0. The normal state energy of the superconductor is ε k ¼ � h 2 k 2 =ð2mÞ À μ S (measured from the chemical potential μ S ) with Δ the pair potential that is diagonal in the pseudo-spin space and m the effective mass. For a nodal-line semimetal with slight doping such that μ � � hλk 2 0 , the Fermi surface possesses a torus shape. Unlike previous sections where only the conduction band is considered, in this section, both conduction and valence bands will be involved in the discussion. The excitation energies for the electron and hole in the nodal-line semimetal are E e ¼ �� hv 0 κ À μ and E h ¼ �� hv 0 κ þ μ, with " � "(" � ") corresponding to the conduction and valence band, respectively.
The nontrivial effect of the torus-shaped Fermi surface on AR is the doubling of reflection channels. By involving both reflection waves and solving the scattering problem at the interface, one can obtain the probabilities for various processes [101]. The types of the AR depend on the magnitude of the chemical potential of the nodal-line semimetal. Specifically, in the limit E; Δ � μ, the electron and AR hole belong to the same band, and two types of AR are the RAR and AAR (with inversion of v x and v z ); for E; Δ � μ, the electron and AR hole reside in different bands, and the AR types become the SAR and AAR (with inversion of v x and v y ). One can see that AAR occurs in both regimes.
For a transparent interface with U ¼ 0, the AR amplitudes in the limit of E; Δ � μ can be obtained as r R ¼ 2χ½2χ cos β cos 2 φþi sin βðcos 2 φþχ 2 Þ� ½i sin βðcos 2 φþχ 2 Þþ2χ cos β� 2 À 4χ 2 sin 2 φ r xz A ¼ À 2χe iφ sin β sin φðcos 2 φÀ χ 2 Þ ½i sin βðcos 2 φþχ 2 Þþ2χ cos β� 2 À 4χ 2 sin 2 φ (13) where r R and r xy A correspond to the RAR and AAR, respectively, with v e x as the x-direction velocity of electron in the nodal-line semimetal and v F x ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi 2mμ S À � h 2 ðk 2 y þ k 2 z Þ q =m as that in the superconductor. In this limit, the RAR process dominates the electron transport for E < Δ. For E ¼ Δ and thus β ¼ 0, r R ¼ 1 and r xy A ¼ 0, meaning that perfect RAR takes place for all incident angles ðθ; φÞ.
For the opposite limit E; Δ � μ, the AR amplitudes reduce to where r S and r xy A correspond to the SAR and AAR processes, respectively. In this case, one can see that the AAR dominates the electron transport with the inversion of v x and v y .
In order to probe the AAR, one should focus on the inter-band AR regime, such that the AAR dominates the transport. The feature of the AAR is its unique reflection trajectory as shown in Figure 8). It was then proposed to use the nonlocal scanning tunneling spectroscopy (STM) to extract such a trajectory of the reflected hole. The AAR can be manifested as a spatially resolved ridge structure of the nonlocal conductance as shown in Figure 9, which is in stark contrast to the peak structure of the RAR [101].

VIII. Surface states induced resonant spin-flipped reflection
The manifestation of the bulk-boundary correspondence principle in nodalline semimetals is the existence of drumhead-like surface states (DSS) [22], as we described in section II. Although the spectra of the DSS can be probed by the ARPES measurements, it is impossible to extract its spatial distribution in this way. On the other hand, it is hopeful to reveal the spatial distribution of the DSS through transport measurements by fabricating electrodes directly coupled to the surface. No matter whether the full energy spectra are gapped or gapless, the surface states are always composed of bound states in the corresponding transverse channels labeled by the lateral momenta k jj . For topological semimetal phases, only the transverse channels that coincide with the nodal points or nodal-lines are gapless, while the other channels are still gapped. Specifically, the DSS in nodal-line semimetals are a set of bound states shown in Figure 1) and Figure 10    nodal-line semimetals with spin-orbit coupling, the spin-polarized DSS correspond to the resonant spin-flipped reflection at the interface between a normal metal and a nodal-line semimetal [105].
Consider the heterojunction in Figure 10) in the z-direction that is composed of a normal metal with spin degeneracy and a nodal-line semimetal with spin-orbit coupling. In the z < 0 region, the normal metal is described by the Hamiltonian H NM ¼ Cjkj 2 À μ 0 , where C is a mass-dependent parameter and μ 0 is the chemical potential that relates to the Fermi wavevector through jk F j ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi μ 0 =C p . In the z > 0 region, the nodal-line semimetal is described by Hamiltonian (1), with the Pauli matrices being the real spin here. The δfunction barrier UδðzÞ is imposed at the interface. Solving the scattering problem at the interface for a given k jj inside the nodal loop yields the spinflipped reflection amplitude as where η ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi v NM =v SM p is the square root of the ratio of perpendicular velocities in the normal metal and nodal-line semimetal, with v NM ¼ 2Cjk F j cos θ and v SM ¼ 2� hλk 0 0 , respectively, θ is the electron's incident angle (relative to the z-axis), Z 0 ¼ 2U= ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi v NM v SM p is a dimensionless interface barriar height, and ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi k 0 0 2 À ð� h 2 v 2 À ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi Þ=ð2� h 2 λ 2 Þ q . In the tunneling limit Z 0 � 1, the probability of the spin-flipped reflection R f ¼ jr f j 2 exhibits a sharp resonant peak around zero energy; see Figure  10). This can be understood by expanding Equation (15) to a series of Feynman paths constructed by multiple reflection processes [105]. It can be shown that the condition of resonant spin-flipped reflection is equivalent to the Bohr-Sommerfeld quantization condition for the surface bound states, meaning that such a resonance is exactly induced by the DSS and thus provides a direct evidence for its detection.
Two transport schemes were proposed to detect such an effect [105]; see Figure 11) and (d). For the spin transport scheme in Figure 11), spinpolarized electrons are injected from a ferromagnetic lead into the normal metal and then get reflected along with spin-flipping at the interface with the nodal-line semimetal. It generates a nearly pure spin current (zero charge current) around the resonance energy as shown in Figure 11, which can be measured by the spin Hall voltage V SH [106] taking advantage of the inverse spin Hall effect [107,108]. For the charge transport scheme in Figure 11), the charge current flowing between two anti-parallel magnetic terminals directly reveals the resonant spin-flipped reflection, which can be manifested as the resonant peak in the differential conductance; see Figure 11).
The spin-flipped reflection is induced by all the bound states that comprise the DSS so that any dispersion of the DSS will result in a broadening of the resonant peak; see Figure 11) and (f). Moreover, the interface disorder has similar effect on the resonant peak [105]. Therefore, in order to observe a sharp resonant peak, nodal-line semimetals with flat DSS are favorable. Nevertheless, even the DSS possess finite dispersion, the spin current will still be strongly enhanced by the resonant spin-flipped reflection, which may lead to interesting application in spintronics.
Finally, the resonant spin-flipped reflection does not rely on the exact form of the Hamiltonian (1), which differs from that of most materials. The key point for such a scenario is the spin polarization of the DSS. It was shown that the effective model of HgCr 2 Se 4 that gives rise to spiral spin texture of the DSS will yield the same results as well [105]. For the nodal-line semimetals with negligible spin-orbit coupling, the DSS may undergo a Stoner ferromagnetic instability [109,110]. As a result, the DSS split into two bands with opposite spin polarization. Whether such spin-polarized DSS will induce resonant spin-flipped reflection remains to be explored.

VIII. Transport experiments
ARPES measurements provide an effective way to obtain the band structures of new materials. This method combined with the first principle's calculation has become a powerful tool to detect nodal-line semimetals [42][43][44][45][46][47][48][111][112][113][114][115][116][117][118][119][120][121][122][123]. Different from topological insulators in which the surface- state spectra are well separated from the bulk ones [1,2], the bulk states in topological semimetals always coexist with the surface states, making the latter difficult to identify. Apart from the band structures, other information such as the geometric phase encoded in the wave functions cannot be extracted by ARPES either.
Transport measurements provide another powerful tool for the investigation of the physical properties of topological materials. Although the current flowing through the sample cannot tell the detailed band structures, it can provide other important information of electronic properties, such as the spatial distribution of the surface states, the geometric phase of the wave function, and the topology of the Fermi surface as discussed in previous sections. The transport measurements have two advantages: first, various field manipulation schemes can be applied in transport experiments; second, some transport signatures such as the scaling of MR stemming from WL or WAL effects and the geometric phase involved in the pattern of quantum oscillation are universal and thus do not rely on the details of the sample.
The MR in a weak magnetic field is dominated by the WL/WAL effects or the Berry curvature effect. Several experiments have reported the MR results that are consistent with the WAL theory of nodal-line semimetal [54] as introduced in section III. Zhou et al [61] reported that the MR of Mg 3 Bi 2 films can be well fitted by the MR formula defined by Equation (6b), which is dedicated for the WAL in nodal-line semimetals in the presence of LR disorder potential [54]. The magnetotransport experiment in the nodal-line semimetal YbCdGe shows that the magnetoconductivity due to WAL effect is five orders larger than that for 2D systems [124], which is also consistent with the theoretical predictions. Similar results with large negative MR was also reported in nodal-line semimetal SrAs 3 [125]. The Berry curvature effect on the low-field magnetotransport in the nodal-line semimetals remains to be explored. To have a notable effect, the Fermi energy and the strength of the spin-orbit coupling should be comparable. This stems from the fact that only the electrons at the Fermi surface contribute to the longitudinal conductance. For many nodal-line semimetals, the spin-orbit coupling is small compared with the Fermi energy and so the Berry curvature effect on the MR is negligible.
As the magnetic field increases, the system enters the quantum regime in which Landau levels dominate the physics. The MR in the stong-field limit exhibits quantum oscillation, known as Shubnikov-de Hass oscillation. By analyzing the oscillation pattern one can extract the shape of the Fermi surface as well as the Berry phase therein, as discussed in section IV. A series of experiments of quantum oscillation on nodal-line semimetals have been performed [69,70,[84][85][86][87][88][89][90][91]. Several works have reported the large unsaturated magnetoresistance [84,85,[87][88][89] as well as an unusual butterfly-shaped titanic angular magnetoresistance [85,87,88]; see for example, Figure 12, which reveals the high anisotropy of the nodal-line semimetal. Moreover, Ali et al. [85] showed that the oscillation may or may not exhibit a nontrivial π Berry phase as the magnetic field is set to different directions. Such directional dependence of the geometric phase is consistent with the theory in Ref [66], as discussed in section IV. Similar results were also reported in Ref [87]. Other experiments on HfSiS [70], SrAr 3 [125], ZrP 2 [126] all observed nontrivial Berry phase extracted from the oscillation patterns. The nodal lines can have rich configurations, which explains the high anisotropy of the physical results. Meanwhile, the Berry phase is a robust effect that is not sensitive to the modulation of the nodal lines, which is the reason why it can be observed in various materials.

IX. Summary and outlook
We have reviewed a variety of interesting transport properties of the nodalline semimetal. These transport phenomena are induced by the unique features of the nodal-line semimetal, which are the π Berry phase carried by the nodal line, the torus-shaped Fermi surface for slightly doped samples and the existence of the DSS. Transport experiments on the WL/WAL effects and quantum oscillation have shown consistency with the theoretical predictions. Among all the transport properties, the results of the WL/WAL effect and quantum oscillations are quite universal, which do not rely on high-quality samples. Actually, the WL/WAL effect resides in the quantum diffusive regime and the scaling of the resultant MR can reveal both the effective dimension of electron diffusion and the Berry phase of the wave function. Quantum oscillation can be achieved as long as the magnetic length is much smaller than the mean free path. Meanwhile, other effects such as the AAR and the spin-flipped reflection do require samples and devices of high quality. These two effects arise only when the system is in the ballistic regime. Otherwise, the conservation of the transverse momentum breaks down and the AAR turns into diffuse reflection and the specific trajectory disappears. Similarly, the signal of the resonant spin-flipped reflection will be strongly suppressed by the disorder in the bulk.
Apart from the quality of samples, the transport signatures rely on the intrinsic properties of the nodal-line semimetals as well. For theoretical studies, a minimal model as Equation (1) is favorable. However, only a few real materials can be captured by such an effective model [9]. In reality, the nodal line may have various configurations, such as closed loops with deformation and open nodal lines [127][128][129][130]. There may also exist finite dispersion along the nodal lines, which makes the Fermi surface deviate from the standard torus shape [131]. The natural question is whether the theoretical results obtained from the minimal model are applicable or not. We give some remarks on this point below.
Complex configurations of the nodal line and finite dispersions along the nodal line can easily break the T 1 symmetry in Equation (4). As a result, the WL effect cannot be achieved even when the disorder potential is of the SR type. At the same time, the local T 2 symmetry, as well as the WAL effect may persist apart from some modifications of the MR formula (6b) by taking into account the different Fermi surfaces for various 2D subsystems. Different configurations of the nodal lines will also change the quantitative results of the quantum oscillations given that the latter is a direct manifestation of the Fermi surface. In principle, the Berry phase can still be extracted from the oscillation patterns by noting that the variation of the Fermi surface only changes the periodicity of the oscillation, i.e. the F factor in Equation (7). However, the complex Fermi surface will make the analysis of the oscillation patterns quite challenging, especially for the Fermi surfaces with nodes, in which the phase factor ϕ 3D in Equation (8) may change alternately along the nodal lines from electron to hole. The nodes of Fermi surfaces may also induce magnetic breakdown for certain magnetic orbits that intersect with the nodes, which needs further exploration. The Berry curvature effect on the magnetotransport of the massive nodal-line semimetals should still exist. However, the coincidence between the Lifshitz transition and the sign reversal of the MR will probably fail.
The change of the nodal lines and accordingly, the torus-shaped Fermi surface will directly affect the AAR effect. The unique trajectory of the AAR in Figure 8) is a direct result of the Fermi torus. As it deforms, the group velocities of the quasiparticles will change dramatically as well. Therefore, in order to implement the AAR, nodal-line semimetals with torus-shaped Fermi surface such as CaTX (T = Ag,Cd; X = As,Ge) [132] and CaAgAs [133,134] are favorable. On the contrary, the resonant spin-flipped reflection induced by the DSS does not rely on the configurations of the nodal lines. Such an effect is directly related to the DSS itself. As long as the spin polarization is established for the DSS, the resonant spin-flipped reflection should be achieved.
Apart from the conventional nodal-line semimetals, the so-called type-II nodal-line semimetals have been proposed recently [135][136][137][138][139]. Similar to the type-II Weyl semimetal [140], the bands of the type-II nodal-line semimetal are tilted as well. The change of the density of states at the nodal line indicates that the transport properties of the type-II nodal-line semimetal should be quite different from the conventional nodal-line semimetals. Very recently, another type of nodal-line semimetal with linked nodal loops [17,18] has been reported [141,142]. The results and scenarios reviewed here may not be applicable to this case because one cannot analyze different nodal loops separately anymore. The nodal links encircle with each other which is topologically distinct from the isolated nodal loops. The transport properties of the nodal-link semimetals such as WL/WAL effect and quantum oscillation are interesting open questions for further study.