Disorder effects on quantum transport and quantum phase transition in low-dimensional superconducting and topological systems

ABSTRACT Disorder effects inevitably exist in realistic samples, manifesting in various physical properties. In this paper, we review the recent progress in understanding the disorder effects on quantum transport and quantum phase transition properties in low-dimensional superconducting and topological systems. As a consequence of the pronounced quantum fluctuation in low-dimensional systems, rare events drastically change the physical characteristics and underlying microscopic transport process in these systems, which are beyond the traditional paradigms. Associating with recent experimental observations, we emphasize the microscopic mechanism for disordered Ising superconductivity, the quantum Griffiths singularity of superconductor metal transition and the discrete scale invariance in topological materials. GRAPHICAL ABSTRACT


Introduction
Recent years, the rapid developments of two-dimensional (2D) superconductor and topological materials have provided new platforms to investigate a variety of quantum phenomena [1][2][3][4][5][6][7]. Exploring these phenomena remains an important pursuit for both experimental and theoretical research. In the early days, the fundamental Hohenberg-Mermin-Wagner theorem points out that the long range ordered phase can not exist in finite temperature 2D system due to the thermal fluctuation [8,9]. However, the Berezinskii-Kosterlitz-Thouless (BKT) transition occurs in a realistic 2D superfluid and superconductor, which allows for the establishment of quasi-long-range correlation of the order parameter [10][11][12]. In recent years, with the advancements of the nanotechnology including molecular beam epitaxy, pulsed laser deposition, scanning tunneling microscopy/spectroscopy and so on, various types of atomic-scale superconducting materials have been fabricated and characterized in unprecedented progress [13][14][15]. Also, the prominent fluctuation and spin-orbit coupling have been proven indispensable in twodimensional (2D) superconductors. As a consequence, 2D superconductors exhibit many exotic properties different from conventional three-dimensional (3D) system, such as the Zeeman-protected Ising superconductivity in 2D superconductors protected by spin-orbit interaction [1,[16][17][18][19][20][21][22][23][24][25], the emergence phase of anomalous quantum metal [2,26], the quantum Griffiths singularity (QGS) [27] of superconductor-metal transitions (SMTs) etc. Explanations for these phenomena are beyond the scope of traditional paradigms.
Meanwhile, large emergence of topological materials have attracted considerable attention in the field of condensed matter physics [5][6][7]. Especially, certain topological semimetals have extremely low carrier concentrations, and are easy to reach the quantum limit under a strong magnetic field, showing anomalous transport properties [28][29][30][31][32]. In addition, topological Dirac or Weyl semimetals show linear dispersion around the Weyl points, and the supercritical Coulomb attraction around the charge impurity can lead to novel phenomena such as supercritical collapse behaviour [29,33]. In topological semimetals, the supercritical Coulomb attraction can lead to quasi-bound state spectra obeying discrete scale invariance, namely the binding energy is in geometric progression. Under magnetic field, novel type of quantum oscillations emerge with periodicity in ln(B), in sharp contrast to the well-known Shubnikov-de Haas oscillations. In general, the characteristics of quantum transport and quantum phase transition in topological semimetals and 2D superconductor are largely enriched by disorder effects, which can provide a platform to testing the new theoretical paradigms and to searching for possible original framework for understanding the anomalous transport features in recent experiments.

Disordered ising superconductivity
The research of superconductivity has always been a hot topic. In 1957, John Bardeen, Leon Cooper and Robert Schrieffer (BCS) proposed a microscopic theory of superconductivity that the attractive electron-electron interaction mediated by phonons gives rise to Cooper pairs and superconductivity can be described as a microscopic effect caused by the condensation of Cooper pairs [34]. Such a theory has many qualitative and a few quantitative points of agreement with related experiments, in particular with specific-heat data, energy gap measurements, and transition-temperature versus impurity curves [35]. To further prove the validity of BCS theory, P. W. Anderson studied the dirty superconductors with elastic scattering from impurities and demonstrated that other types of pairing except for the s-wave pairing are not compatible with the existence of dirty superconductors [36]. Meanwhile, A. A. Abrikosov and L. Gor'kov studied alloys containing paramagnetic impurities based on the quantum field theory [37], and had shown that the presence of paramagnetic impurities can diminish the s-wave pairing superconductivity. The dependence of T c on the paramagnetic impurity concentration can be derived via the Gor'kov Green function [38][39][40]. Moreover, apart from magnetic impurity, there exist other pair-breaking mechanisms of the conventional superconductor. Clogston and Chandrasekhar had illustrated under magnetic field the Zeeman energy exceeds the condensation energy corresponding to the so-called value of Pauli limit for the magnetic field [41][42][43]. Subsequently, Werthamer, Helfand and Hohenberg (WHH) generalized the Abrikosov and Gor'kov theory to include the orbital effect, the paramagnetic effect, and the spin-orbital scattering effect [44].
Nevertheless, the WHH theory mainly focus on the 3D superconductors, which is not suitable to describe low dimensional superconductors [16][17][18][19][20]. In layered superconductors, the reduction of dimensionality weakens the orbital effect when the magnetic field parallels the layered plane, thus the spin-orbit scattering can randomize the spin orientation and lead to enhancement of in-plane critical field B c 35 . Klemm-Luther-Beasley (KLB) had firstly proposed the microscopic theory applicable to layered superconductors in the dirty limit with strong spin-orbit scattering [45,46]. Afterward, with the advancements of film-growth techniques, various types of high-quality 2D crystalline superconducting films have been fabricated [13][14][15]. One striking discovery is the Ising pairing mechanism with pronounced enhancement of in-plane upper critical field B c beyond the Pauli limit, such as ion-gated MoS 2 , monolayer NbSe 2 . However, this enhancement of in-plane B c can not be explained by the KLB theory, because the value of the spin-orbit scattering time estimated from the KLB theory turns out to be much shorter than that of the total scattering time obtained from transport data [16]. In MoS 2 and NbSe 2 , the atomic spinorbital interaction (SOI) is very strong due to the heavy transition metal atoms. And inversion symmetry breaking and SOI lead to the spin polarizing perpendicular to the layer plane around the Dirac points, namely the type-I (inversion-asymmetric) Ising superconductors [1,[16][17][18][19][20][21][22][23]. However, in addition to the inversion-asymmetric Ising pairing, certain 2D inversionsymmetric materials with spin splitting around the Γ point due to intrinsic SOI are very recently found capable of superconductivity [47] [see Figure 1]. Systematic models are needed to investigate the pair-breaking mechanism and in-plane B c in both the inversion-symmetric and inversion-asymmetric Ising superconductors, taking into account of the combination effects of spin-orbit interaction and impurity scattering.

The microscopic theory of two-dimensional Ising superconductivity
To improve agreement between theory and experiment, a general microscopic theory for 2D Ising superconductors is proposed in Ref. [48], which generalize the Ising superconductivity theory containing the intrinsic SOI, the impurity scattering, and the Rashba SOI simultaneously. This theory stem from the Gor'kov equation of transition temperature T [44,49]: where the system has translation symmetry after the impurity average and the integral kernel function S ω ðpÞ; k B T Δ ð d�F ω ðpÞ is assumed to be isotropic. F ω ðpÞ is the anomalous Green's function in momentum space.

Type-I (Inversion-asymmetric) Ising superconductivity
We firstly discuss the inversion-asymmetric Ising superconductivity caused by the Zeeman-type intrinsic SOI. Taking 2D hexagonal lattice as an example [see Figure 2(a,b)], the intrinsic SOI serves as effective out-ofplane magnetic field which takes opposite value at valley K and K 0 , inducing the spin spiltting of energy bands, namely the Zeeman-type SOI. Around each valley, the Zeeman-type SOI can give rise to double Fermi surfaces (FSs) with almost the same radius and shape. Then the electrons at K þ k r and À K À k r (r ¼ 1; 2 labels the two FSs) with opposite out-of-plane spins form the inter-valley Cooper pairs obeying the time-reversal symmetry. When applying an in-plane magnetic field Bx, the Hamiltonian reads [50]: where m is the effective mass, σ and τ denote the real spin subspace and valley subspace, respectively. Then the relation between critical field B c and temperature T can be solved within the WHH framework. To be specific, if considering Zeeman-type SOI protected superconductivity with spin-independent scattering, the spin-independent scattering disorder gives a finite lifetime � h τ 0 ¼ 2πn i Nð0ÞV 2 with n i denoting the disorder concentration. Here, one can introduce effective Zeeman-type SOI parameter f β SO ¼ β SO 1þ� h=2πk B T C τ 0 , and the upper critical field B c ðTÞ for the Zeeman-type SOI protected superconductivity in the dirty limit reads [21,48]:  shows the results with a fixed effective Zeeman-type SOI and increasing effective Rashba-type SOI. The black curve presents the best fit of the experimental data with effective Zeeman-type SOI f β SO ¼ 3:16 meV and effective Rashba-type SOI g α R k F ¼ 0:22 meV. The olive line indicates a special case with only effective Rashba-type SOI ( g α R k F ¼ 30meV, f β SO ¼ 0 meV), which is far below the experimental results. Reproduced from Ref. [21] If taking both the spin-independent scattering and the spin-flipping scattering into account, the characteristic time satisfies τ À 1 ¼ τ À 1 0 þ τ À 1 SO , and τ, τ 0 , τ SO denote the mean free time for total scattering, spin-independent scattering and spin-flipping scattering, respectively [21]. Besides, Rashba-type SOI tends to polarize the spin to the in-plane direction, which weakens Zeeman-protection mechanism and changes the B c . The Hamiltonian including the Zeemantype SOI and Rashba-type SOI simultaneously at the K and K 0 points are [21]: Likewise, one can define the effective Rashba SOI g α R k F ; 1þ� h=2πk B T C τ 0 , the relationship between T=T c and B c can be given as well. This microscopic calculation is consistent with previous phenomenological formula [16] after the setting the disorder renormalized value of related parameters.

Type-II (inversion-symmetric) Ising superconductivity
Different from inversion-asymmetric Ising superconductors, like Figure 3(a), the energy valley locates at zero-momentum point Γ in inversion-symmetric Ising superconductors. The intrinsic SOI can split the 4-fold degenerate P þ x;y level into two doubly-degenerate levels, opening a gap at the Γ point [14,51]. The lower level crosses E F at two different Fermi wavevectors form two different-shaped FSs [14], each of which holds two states [see Figure 3(a, b)]. If the radius of the FS is small, the eigenstates of the lower energy level can be approximated by P þ xÀ iy;" , P þ xþiy;# . The SOI at valley Γ can be viewed as an out-of-plane magnetic field which takes opposite value À B effẑ and B effẑ on different orbits P þ xþiy and P þ xÀ iy , respectively. Thus, the system has TRS at B ¼ 0 and electrons with opposite momenta and spins on the same FS can form swave Cooper pairs. These pairs with out-of-plane spin-orbit locking give rise to the large in-plane B c . Based on previous study [52], the inversion-symmetric system can be represented by a four band model with basis (P þ xþiy;" , P þ xÀ iy;" , P þ xÀ iy;# , P þ xþiy;# ) with external in-plane magnetic field Bx, Here μ B is the effective Bohr magneton and A, M 0 , M 1 , v are fitting parameters. Considering the situation of lower band E À ðkÞ crossing E F at two different Fermi wave-vector and the spin-independent scattering disorder within each FS, the temperature dependence of the in-plane upper critical field read as [48]: 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 is the disorder renormalized SOC strength, k F , ψ, τ denote the Fermi momentum, the digamma function and the scattering time, respectively. Near T c , this result is consistent with the result from 2D Ginzburg-Landau (GL) theory [53] and the KLB theory [45,46]. . The arrows on the FSs denote spin directions along � z direction, k 1 in blue and k 2 in orange. Each FS is doubly degenerate, consisting of two different P þ x;y orbitals and out-of-plane spin directions. The pairing happens between P þ xÀ iy;" and P þ xþiy;# on the same FS if k r is small. (b) 4-fold degenerate P þ x;y level of stanene splits to two doubly-degenerate levels due to SOI. The lower level forms two FSs at k 1 (blue) and k 2 (orange). (c-e) show the temperature dependence of the in-plane upper critical fields in few-layer stanene samples. Data obtained from four samples with different stanene and PbTe substrate thicknesses. For example, 3 À Sn=6 À PbTe refers to a trilayer stanene grown on top of 6 layers of PbTe. The ratio of the in plane upper critical fields. i. e. the magnetic fields at which the sample resistance becomes 50% of the normal state resistance at a given temperature, to the Pauli limit field B ¼ 1:86T c;0 are plotted as circular symbols. Solid/dashed curves are theoretical fits by .Eq. (6). T c;0 is the zero-field transition temperature. β � SO is the intrinsic SOC strength renormalized by disorder. αk F denotes the renormalized Rashba SOC strength. Reproduced from Ref. [47]

Type-I Ising superconductors
It has been a significant discovery that ion gated MoS 2 and NbSe 2 can be superconducting [64,65]. Unlike the bulk transition metal dichalcogenides (TMDs) that preserve inversion symmetry, monolayer MoS 2 and NbSe 2 are inversion-asymmetric superconductors [66]. Moreover, the strong SOC in these systems give rise to Ising pairing mechanism. In addition to the aforementioned mechanism, the interface can also give rise to the type-I (inversion-asymmetric) Ising pairing. Ref [21]. has reported the observation of interface-induced Zeeman-protected superconductivity in ultrathin crystalline Pb films epitaxially grown on the striped incommensurate phase on a Si(111) substrate. The RðBÞ curves of a 6-ML sample are measured in a pulsed parallel magnetic field up to 47 T. The superconductivity survives under a high parallel magnetic field of 40 T at 1:7 K. The experimental data of B c =B p versus T=T c can be quantitatively fitted by Eq. (3). The fitting procedure gives effective Zeeman-type SOI f β SO ¼ 3:01 meV, as shown in Figure 2(c). The enhancement of in-plane critical field can be well explained by the type-I (inversion-asymmetric) Ising superconductivity mechanism. Moreover, in order to clarify the influence of Rashba-type SOI, Figure 2(d) displays a set of theoretical curves with a fixed effective Zeeman-type SOI f β SO ¼ 3:16 meV and increasing effective Rashba-type SOI g α R k F from 0.22 to 1.26 meV. The Rashba-type SOI bends the curves down at low temperatures and gives rise to a relatively small in-plane critical field, which deviates from the quasi-linear temperature dependence of measured data, indicating Rashba-type SOI alone cannot account for the observed large in-plane critical field.

Type-II Ising superconductor
Recent experiment has identified a pronounced enhancement of in-plane critical field B c at low temperature in epitaxial few layer α-Sn (111) [14,15,47]. Distinct from the previous system, few-layer stanene has no M z mirror symmetry and is centro-symmetric in the free standing case [51]. Stanene hosts bands around Γ-point instead of the K and K 0 points as in type-I Ising superconductor. Thus, it needs an alternative mechanism in stanene to produce the out-of-plane spin orientations.
According to the results shown in Figure 3, the experimental data matches well with Eq. (6), which captures the prominent up-turn feature in the low-T regime and helps to qualitatively understand the substrate and layer-thickness dependence of B c;== . Such pairing mechanisms behind Eq. (6), without the participation of inversion symmetry breaking, are termed type-II (inversionsymmetric) Ising superconductivity. In fact, the upturn of B c;== in low-T result from frozen spin orientation at T ! 0. When the temperature goes up, thermal activation results in a partially polarized system, suppressing the contribution of the spin-orbit induced spin split effect on B c;== . Figure 3 (c-e) show the temperature dependence of the in-plane upper critical fields in few-layer stanene samples. Comparing to 3Sn=12 À PbTe, the smoothening in 3 À Sn=6 À PbTe can be attributed to the variation of the spin locking strength as one moves away from the Γ point along the inverted Mexican hat band shape (inset to Figure 3 (c)). Spins of the j þ > and j À > orbitals are strongly locked out-of-plane at the Γ point. This Ising-like orientation becomes however less favorable at larger momenta. Lowering the Fermi level thus suppresses the spin polarization of the outer hole band, which can be simulated by an effective Rashba term in the Hamiltonian. Taking into account this effect, the modified formula nicely describes the upper critical fields of bilayer stanene (Figure 3d). Since stronger inversion symmetry breaking attribute to the missing up-turn feature, a pentalayer stanene should experience a weaker Rashba effect, giving rise to an apparent enhancement of B c;== at low T. This indicates a broader range of inversion symmetry materials host such pairing mechanisms. Recent experiment in ultra-thin crystalline PdTe2 films also demonstrates clear feature of type-II Ising pairing with in-plane critical field more than 7 times the Pauli limit [67].

Discussion
There may be multiple types of SOI, including the Zeeman type SOI and the Rashba type SOI both may affect the in-plane B c 16, [17], Zeeman-type SOI can cause an out-of-plane spin polarization of energy valleys; hence, it protects the superconductivity under a high in-plane magnetic field both in type-I and type-II Ising superconductor. In contrast, the Rashba-type SOI alone cannot account for the observed large in-plane critical field. Besides, the weak Rashba-type SOI tends to polarize the spin to the in-plane direction, making the Cooper pairs more susceptible to the in-plane magnetic, and destructs the upturn at very low temperature. Lastly, the combination effect large intrinsic SOI and quantum fluctuation may also dramatically influence the perpendicular critical field B c;2 , including the quantum Griffiths singularity along the anomalous phase boundary [68]. The influence of large intrinsic SOI on the quantum fluctuation needs to be considered with the second order corrections [68,69], which is beyond the WHH framework.

Disorder effect on superconductor metal transition: quantum griffiths singularity
The 2D superconductor not only help to explore new superconductivity characteristics, but also exhibit some properties largely deviated from the prediction of mean field theory, such as quantum Griffith singularity (QGS) associated with the superconductor-metal transition [1,27,[70][71][72][73][74][75][76][77][78][79][80]. Before introducing this novel property, we briefly introduce the basic theoretical concepts of continuous phase transition and critical phenomena. Phase transitions occur upon variation of an external control parameter. Classical phase transitions are characterized by temperature, and the macroscopic order is destroyed by thermal fluctuations at finite temperature. In contrast, quantum phase transitions take place at T ¼ 0 in which the control parameter is non-thermal, such as pressure, magnetic field, or chemical composition [81]. When approaching the critical point of continuous phase transition, the characteristic length scales become divergent [82]. Upon approaching the critical point, the correlation length � and the correlation time diverge as [83,84]: where r, B À B c ð Þ=B c j j is dimensionless parameter and B c is the critical parameter, v is the correlation length critical exponent and z is dynamical critical exponent. Power-law singularities in the length scales at the critical point generically lead to power-law singularities in observables [85]. The disorder effects may dramatically change the power-law singularity into other kind of universality class, i.e. the quantum Griffiths singularity, as shown in the following subsections. Before demonstrating this anomalous case, we firstly outline the general considerations of scaling theory for continuous phase transition.

Scaling theory for continuous phase transition
The first unification of phase transitions theories was achieved by L. D. Landau [86,87]. He proposed a mean-field order parameter m describing phase transition behavior at the transition point, m ¼ 0 represents a disordered state while m�0 represents an ordered state. The Landau theory successfully predicts the critical behavior of phase transitions in many systems and intuitively explains the universality of them [84]. Nevertheless, it often deviates from the experimental results, because Landau theory does not take into account the spatial fluctuations or the long-range correlations of order parameters of reality system. Hence, the Landau-Ginzburg-Wilson theory is proposed, which can be attacked successfully using modern renormalization group techniques [88].
The power law behavior at critical point implied the scaling theory. Widom firstly put forward as the scaling hypothesis [89]. The system is scale invariant close to criticality. If we rescale all lengths in the system by a common factor, the physical properties become unchanged. Therefore, the free energy density can be written as [82]: where r denotes the distance from the critical point, h is dimensionless measure of external magnetic field, d is the space dimensionality, v and y h are two critical exponents and b is an arbitrary positive number represents the scale factor. Since all thermodynamic observables can be determined from free energy density that only depends on two exponents, the various exponents are therefore connected by scaling relation [89]. Indeed, the critical exponents are determined only by the symmetries of the Hamiltonian and the spatial dimensionality of the system. The scaling theory indicates the universality of critical exponents [90], and the related concept of universality classes can be well explained within the framework of the modern renormalization group theory [91]. Scaling theory for classical phase transition can be generalized to quantum phase transition (QPT). In general, the kinetic and potential part of Hamiltonian does not commute in quantum statistical mechanics. Therefore, quantum Landau-Ginzburg-Wilson functional needs to be formulated in terms of space and time dependent fields. Then, the imaginary time variable τ acts as an additional coordinate [84]. Consider the fluctuations in the imaginary axis, � t / � z , if we scale length according to L ! bL, the imaginary time need scale like τ ! b z τ, and temperature regarded as an energy scales like an inverse time T ! b À z T. Consequently, the scaling form of the free energy density at a quantum critical point reads [83]: Here, T is the absolute temperature, b is the dimensionless parameter (i.e. the dimensionless magnetic field), and z is dynamics exponent. Comparing Eq. (9) to Eq. (8), one can directly find that a quantum phase transition in d dimensions is equivalent to a classical transition in d þ z spatial dimensions [92].

Quantum Griffiths singularity
Many realistic systems contain a certain amount of quenched disorder, viz., the impurities or defects that are spatially non-uniform but not varying over the relevant experimental time scales [93]. How does quenched disorder influence phase transitions and critical points? Harris [94] derived a criterion for the stability of a given critical behavior with respect to quenched disorder. For short range correlated disorder, if space dimensionality d and clean correlation length exponent v fulfill the Harris criterion dv > 2, weak disorder does not affect the critical behavior of the system, corresponding to the clean limit [74,[94][95][96]. On the contrary, if dv < 2, weak disorder becomes important and the critical exponent change, namely the dirty limit. Note that Harris criterion only deals with the average behavior of the disorder at large length scales. Furthermore, the spatial non-uniformity random disorder may lead to off-critical singularities at phase transitions [97][98][99][100], similarly to the features of random transverse field Ising model [101]. This kind of singularities are called the Griffiths singularities, because R. B. Griffiths firstly proved its existence in 1969 [102]. T. Vojta has provided the whole classification of disorder effect on the scaling behaviours, based on the analysis of the dimension of range regions [84]. Specifically, when Harris criterion is violated, besides the situation of phase transition destruction at strong disorder, there exist two possible critical points according to the behavior of the average disorder strength under coarse-grain [98].
(i) The disorder strength approaches a finite value under coarse-grain. It results in the divergence behavior of quantities near the critical point still obey power-law but with a different set of critical exponent (the new correlation length exponent v 0 fulfill the inequality dv 0 > 2), which differ from those of the clean system.
(ii) The disorder strength diverges under coarse-grain, the transition is characterized by an exotic infinite-randomness critical point. Surprisingly, the temporal correlation length � t and spatial correlation length satisfy a activated (exponential) scaling � t / e � ψ , distinctive from the relation � t / � z in power law scaling. Here ψ is the activated dynamical index.
The second special case is relevant to recent experiments in 2D superconductors. A typical model demonstrating the activated (exponential) scaling behaviour is the random transverse-field Ising model, as introduced in the following subsection.

Random transverse-field Ising model
Let us start by discussing the effects of rare strong spatial disorder fluctuations on phase transitions. Considering the example of a diluted classical ferromagnet, the dilution reduces the tendency towards magnetic longrange order and thus reduces T c from its clean bulk value T c;0 . In an infinite system, one can find arbitrarily large spatial regions that are devoid of impurities. These regions will show local magnetic order even though the bulk system is globally in the para-magnetic phase, called rare regions [84].
Long-range disorder correlations can immensely increase the rare region effects. This rare region effect was firstly well depicted in the McCoy-Wu model [103], a 2D classical Ising model with linear defects. In quantum phase transitions with quenched impurities, the systems are therefore always correlated at least in the temporal dimension. Thus, rare region effects in a quantum phase transition are expected to be stronger than in a classical one. One of the typical random quantum systems models displaying a quantum phase transition is the random transverse field Ising model, the Hamiltonian reads: where the J i;j are random positive interactions, the random transverse fields h i cause the quantum fluctuations. At temperature T, such a quantum system in d-dimensions is equivalent to a d þ 1 dimensional classical Ising model with the couplings independent of the "time" direction in which it has extent 1=T. The behavior of the random transverse field Ising model was not fully understood until D. S. Fisher [101,104,105] applied a strongdisorder renormalization group method to the problem. The main idea of the strong-disorder renormalization group is to take the strongest coupling ω ¼ maxðh i ; J i;j Þ in the system and coupling to the rest of the system is treated perturbatively [91]. After integrating the strongest coupling bond (or site), one can eliminate one bond (or site) and derive a new effective Hamiltonian with fewer degrees of freedom. The repeated application of this procedure establishes an exponential relationship between the length scale L of the clusters and effective bonds scale with the logarithm of the energy scale L ψ , ln ω 0 =ω ð Þ; with the exponent ψ, and ω 0 represents the basic energy scale. Thus, the random transverse field Ising model displays activated dynamical scaling rather than conventional power-law scaling. The critical fixed point in the random transverse field Ising model is also called infinite randomness fixed point. This peculiar activated scaling property originates from the property that the probability for finding a large spatial region of volume V decreases exponentially ,expðÀ pVÞ 84 . Moreover, if a system is pushed slightly off critical with a distance δ, the correlation lengths satisfies �,jδj À ν with ν. In the quantum Griffiths phase, the activated scaling law results in a typical energy scale ω, which obeys lnð ω ω 0 Þ,� ψ . If the conventional scaling analysis is applied to the finite temperature regime in the quantum Griffiths phase, one can obtain a diverging effective dynamical exponent z satisfying [105]: z,δ À ψν ; (12) which indicates the feature of Griffiths singularity. In other words, Griffiths phase is dominated by the rare regions in which the probability of correlation length � proportional to e À p� d , (p is parameter of system, d denotes dimensionality). The rare regions cause the non-analytic behavior of the thermodynamic potential function before reaching the critical point, thus, break the power law scaling law of the system [19,106]. Moreover, the activated scaling analysis can be applied directly to elucidate the quantum Griffiths phase [107][108][109].

Observations of quantum Griffiths singularity in SMT experiments
The superconductor insulator/metal transitions, representing typical examples of quantum phase transitions (QPTs), have been observed in several quasi-2D disordered superconducting systems [110][111][112]. Such QPTs are expected to be realized by tuning the film thickness or the applied magnetic field in these systems. The critical magnetic field B c of superconductor insulator/metal transition is considered equal to the crossing point of sheet resistance R S ðBÞ curves at different temperatures [74,113]. Recently, it is striking that two crossing points of R S ðBÞ curves were observed in superconductor metal transitions (SMTs) in 2D superconductor thin films [76,114]. Moreover, finite size scaling analysis elucidates two different values of critical exponents around these crossing points. This discordance limits the feasibility of universality class, and new theoretical paradigm is needed to clarified this inconsistency. In 2015, Xing, et. al. firstly discovered the QGS in 2D superconducting Ga films grown on a 3-mm GaN(0001) substrate in an ultrahigh-vacuum molecular beam epitaxy chamber [27]. In Figure 4, the magneto-resistance RðBÞ curve of the sample under different temperatures do not cross at one point [110,111,115], but forming a crossing line. Meanwhile, there is an upturn of B c at T ! 0 beyond the WHH theory. Utilizing the finite-size scaling (FSS) method for a bundle of R S ðBÞ curves around the multiple crossing points, It turns out the critical exponents zv exhibit a power-law divergence when B ! B c at T ! 0, which matches well with QGS feature [105]. In addition, Figure 4(c) sketch the B À T phase diagram of 2D SMT with quenched disorder in Ga film. The anomalous upturn of upper critical field as T ! 0 appears. In the conventional theory of homogeneous type II superconductors, the Abrikosov vortex lattice forms in the region of B 1 ðTÞ < B < B 2 ðTÞ, characterized by a continuous phase transition at the upper critical field B 2 ðTÞ [38,116]. In clean systems, thermal fluctuations give rise to phonon-like excitations of the vortex lattice. At T > T M regime, thermal fluctuations become sufficiently large, the system is expected to melt into the vortex liquid phase [117]. While at T < T M , the effect of quenched disorder Reproduced from Ref. [27] overtakes the thermal fluctuation effect and produces a vortex glass-like phase on length scales greater than characteristic length of disorder [118]. This vortex glass-like phase consists of spatially separated superconducting regions, where the local upper critical field exceeds the average mean-field value B 2 ðTÞ. In the ultralow temperature regime, these separate superconducting islands, namely, the rare regions can couple via long-range Josephson coupling and manifest global superconductivity in the system. These rare regions of superconductivity cause the peculiar upward turn of upper critical field as T ! 0. Shortly after, in 2016, Shen, et.al observed the SMT driven by a perpendicular magnetic field in superconducting 2D electron gas can form at the LaAlO 3 =SrTiO 3 ð110Þ interface and the dynamic critical exponent approaching the quantum critical point also exhibits QGS [119]. Subsequently, QGS has been reported in a wide range of crystalline 2D superconductors, such as NbSe 2 19 , ZrNCl and MoS 2 120 . Different from traditional superconductor Ga film with weak SOC, NbSe 2 monolayer film is Ising superconductor with strong SOC. The discovery of the QGS in the NbSe 2 thin film not only shows the universality of the QGS in the 2D superconducting system, but also shows that the QGS is not sensitive to the specific superconducting mechanism. The coexistence of Ising superconductivity and QGS at 2D limit could be an important topic in the future for further understanding 2D superconductivity [121][122][123].
Distinct from previous reports of QGS, in 2019, Liu, et. al. have observed the QGS associated with the anomalous phase boundary of SMT in ultrathin crystalline Pb films at ultra-low temperatures [68]. Usually, the SMT shows monotonic phase boundary separating the regions of dR=dT < 0 and dR=dT > 0. However, R s ðTÞ curves in 4-ML Pb films exhibit a remarkable reentrant behavior at lower temperatures in Figure 5(a). When applying the magnetic field slightly large than 3.3 T, the sheet resistance firstly decreases with decreasing temperature, reaches the minimum at T min and then rises at a lower temperature. This re-entrant behavior, including the fact that T min increases with increasing magnetic field, can be quantitatively reproduced after taking the superconducting fluctuation and the influence of large SOI into account [125,126]. Also, the R s ðBÞ curves cross each other in a relatively large transition region instead of a critical point. The measured SMT phase boundary in 4-ML Pb film is anomalous since it bends down at low temperatures > 2:44 K as shown in Figure 5(b), which can be ascribed to superconducting fluctuation as well. The SOI leads to the changes in the Cooperon propagator by introducing spin-triplet channel in addition to the spin-singlet channel [126]. When approaching the quantum critical point along the anomalous phase boundary, large superconducting islands emerge and dominate the dynamical property of the system, the divergent behavior of zv matches well with QGS behavior [97,99].
The QGS phase can be directly verified by the activated scaling analysis [108,109]. In practical analysis, certain irrelevant scaling variables influence the crossing point [127], and needs to be considered. The activated scaling with irrelevant corrections reads: Here ν and ψ are critical exponent and À y is related to the irrelevant parameter. Both the relevant and irrelevant parameters are associated with the activated scaling form with effective length scale ln T 0 T . The irrelevant parameter also changes the phase boundary of the system: The influence of irrelevant parameter is pronounced in related experiments, but the irrelevant parameter is ignored in the previous theoretical consideration. The activated scaling analysis with irrelevant critical exponents is utilized for analyzing recent SMT experiments in 2D superconducting thin films [108,109].

Discussion
In the previous subsection, we mainly mention the disorder effect on SMT, which is associated with the pair-breaking of cooper pairs and the corresponding critical resistance is much smaller than the resistance quantum On the other hand, there exist a large class of QPT with critical resistance near R Q , namely the SIT. The disorder and dissipation effects may dramatically change the quantum critical behaviour in SIT, and give rise to the anomalous metal phase in recent experiments. This very fascinating field is beyond the scope of this review article.

Discrete scale invariance in topological semimetals with charge impurity
Topological materials are hot topic comparable to superconductor in exploring novel quantum phenomena. The related studies can not only deepen our understanding of quantum behaviours in solids, but also provide possible applications for future information technology [7,128]. Recently, the discrete scale invariance (DSI) behavior in Dirac materials has attracted attention in related subfields of physics [28,[129][130][131][132]. Especially, a new magneto-resistance oscillation beyond the quantum limit [28,30,133] has been found in several high-quality topological materials. Different from the previously known periodic quantum oscillations (Shubnikov-de Haas oscillations in 1=B periodic [134,135], Aharonov-Bohm and Altshuler-Aronov-Spivak effects induced quantum oscillations in B periodic [136]), this novel oscillation is lnðBÞ periodic which manifests the appearance of DSI in solid-state system [28,30,133]. Such a peculiar DSI feature is considered to be universal in Dirac materials with Coulomb attraction [28,29], which may be closely related to the quasi-bound states formed by massless Dirac fermions and the longpursue atomic-collapse phenomenon [28,131,137]. Thus, it is desirable to explore the underlying mechanisms of this novel phenomenon.

Heuristic analysis for DSI in topological semimetal with charge impurity
Scale invariance is a prominent ingredient in the area of phase transitions. It means a field ϕðxÞ has the property that ϕðvxÞ ¼ v d ϕðxÞ, where v is a continuous parameter and d is the scaling dimension. When considering the quantization effect, the continuous scale invariance will be spontaneously broken down to the DSI, where v takes on discrete values [138]. Its applicability to bound-state problems in quantum mechanics was exploited by Vitaly Efimov, who found the 1=r 2 potential suited to Schrödinger equation and confirmed the DSI energy spectrum could occur in 3-body problem in the limit of resonant scattering [139,140] [see Figure 6]. Analogy to Efimov 3-body problem, topological semimetals hold Weyl equation that the energy and momentum satisfy E / k and DSI is also expected to achieve with 1=r potential [141].
Considering the Weyl fermions with Coulomb attraction near the Weyl point, the Hamiltonian read [29]: with potential V ¼ � hv F α r . Here α; Ze 2 2� hv F denotes fine structure constant, v F is the fermion velocity, Ze is the charge on the impurity and 2 is the dielectric constant. At B ¼ 0, then Eq. (15) can be solved analytically [142,143]. In the case of attractive potential α < κ (κ is a dimensionless quantum numbers of angular momentum), the energy levels of quasi-bound states exhibit DSI [129,130,144,145], At B�0, the magnetic field introduces a new characteristic length l B ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi � hc=eB p that breaks the DSI at large r > l B . The problem with the finite magnetic field is no longer exactly solvable but can be deal with based on Wentzel-Kramers Brillouin (WKB) approximation [29]. In this quasi-classical picture, the energy change occurs when l B ¼ � hv F γ jEj with γ ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi κ 2 À α 2 p . And since E satisfy Eq. (16), the characteristic values of B meet B nþ1 =B n ¼ e À 2 π γ . As shown in Figure 7, with the magnetic field increasing, the binding energy of resonances states cross the Fermi energy at discrete values of the magnetic field which constitute a geometric series. This indicates the lnðBÞ periodic transport property [29].

Influence of screening effect, finite mass and Zeeman effect on the DSI
Above studies neglect the influence of screening effect, finite mass and Zeeman effect on DSI which are exist in realistic system. We generally discuss the influence of these effects on the DSI. Considering screening effects, the impurity potential is V ¼ � hv F α r e À r=� within Thomas Fermi approximation, where the screening length � obeys the relation � À 2 ¼ 4πe 2 dn dμ , with n; μ and g denoting the carrier density, the chemical potential, and the degeneracy, respectively. The results point out the introduction of the length parameter � in V causes the breaking of DSI and when � ! 1, DSI is restored. Moreover, the finite mass Only the real part of five quasi-energies are shown. The broad quasi-bound states evolve into sharp resonances in the presence of the magnetic field, and approach the Fermi energy at certain sharp values of the magnetic field strength B n that form a geometric series. Reproduced from Ref. [29].
term also introduce a new length scale and it may break the DSI down. Solving the radial equation with mass term m 0 ¼ Mv F =� h, the spectrum still satisfies Eq. (16). The detailed results illustrate the DSI quasi-bound states are robust against a small finite mass gap at the Weyl point. Lastly, the Zeeman effect needs to be considered at large magnetic fields. Take the previous analysis in ZrTe 5 as an example [146], the Zeeman spin splitting can influence the lowest Landau level [29]. When magnetic field is applied perpendicularly on the layered ZrTe 5 sample. The Zeeman field gives rise to a spin-polarized conduction band and valence band. The influence of the Zeeman energy can be represented by the effective Fermi energy E F ¼ 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 . When magnetic field is applied parallel on the layered ZrTe 5 . The Zeeman field only shifts the wave vector parallel to the magnetic field and the Zeeman effect does not influence the Fermi energy in this case.

Ln(B) periodic oscillation
The lnðBÞ periodic transport property of system can be captured by the Tmatrix approximation for the Kubo formula [29]. According to the classical treatment of Bastin et al. [29,147], the formula for the transverse magnetoconductivity σ xx at zero temperature is [29]: where UðBÞ ¼ t 2 8π� hv, À 1 � ΓðBÞ . n s is the density of short-range scatterers, n C is the density of Coulomb scatterers, E n are the quasi-bound state energies as calculated and shown above, and ΓðBÞ is the width determined by overlap with the lowest Landau level. l � is the effective length along the magnetic field, t is the coupling strength between the bound states with the lowest Landau level states.
The above microscopic formula can be further simplified into an empirical form by parameterization which is more suitable for fitting the experimental data [29]: with fitting parameters η denoting an effective broadening factor and B 0 denoting a characteristic magnetic field. The longitudinal magneto-conductivity can be obtained by relation ρ xx � σ xx σ 2 xy . The conductivity in Eq. (18) contains two terms. The first term is short-range impurity scattering, giving linear-B magneto-resistance(MR), regarded as background. The second term is resonant scattering between mobile carrier and the quasi-bound state, which leads to lnðBÞ periodic MR oscillations. This new type of log-periodic MR oscillations was observed beyond the quantum limit in ultra-quantum topological materials ZrTe 5 . Figure 8(a) shows the extract oscillations in different ZrTe 5 samples (s6, s7, and s9) [28]. Using formula Eq. (18) can quantitatively fit these lnðBÞ periodic oscillations (black curves). To better presentation, the magnetic fields (B n ) corresponding to MR peaks/dips were extracted, and the dependence of lnðB n Þ on n are plotted in Figure 8(b). Here, the MR peaks and dips were indexed by integers ðnÞ and half-integers ðn À 0:5Þ. The well-defined linear behavior for the lnðB n Þ vs n curve indicates the clearest feature of DSI. This feature can be understood vividly like Figure 8(c). ZrTe 5 are high-quality single crystals with no signature of real-space fractal or strong disorder-induced multi fractal properties. Therefore, the observed DSI is very likely a notable aspect of the quasi-bound states with geometrical scaling. Owing to the ultra-low carrier density in ZrTe 5 samples, the absence of screening effect can give rise to the Coulomb attraction, and the small Fermi velocity further guarantees the supercritical collapse condition, which in combination, result in the two-body quasi-bound states with DSI.
Such peculiar work also inspires people to explore the log-periodic quantum oscillations and the DSI in other physical observables and topological materials. Ref [30]. reported the logB-periodic quantum oscillations in both the longitudinal MR ðρ xx Þ and the Hall resistance ðρ yx Þ of HfTe 5 crystals, which identify the lnðBÞ phenomenon in the total transport coefficients as a general quantum effect. Recently, in the small gap semiconductor tellurium with a relatively small carrier density and unique band topology [148] found similar logB-periodic longitudinal MR and the Hall resistance oscillations to those observed in HfTe 5 . The small gap does not change the DSI feature in this system, which is consistent with previous theoretical analysis [33]. Moreover, the MR quantum oscillations TaAs 133 is considered as a possible candidate for logB-periodic oscillation, revealing the universality of DSI property in Dirac/Weyl materials with supercritical Coulomb attraction.

Outlook
Dirac materials provide new platforms to explore the ultra-relativistic quantum phenomenon due to their massless spectrum. The effective speed of light in these Dirac materials are far less than the speed of light in nature, and thus the effective fine-structure constant can be greater than 1, providing an ideal platform to investigate the supercritical phenomena under tunable electric field and magnetic field [129,131]. Furthermore, the novel DSI of quasi bound states in massless Dirac particles manifests as a remarkable feature that is distinctive from the previous prediction of supercritical atomic collapse in massive Dirac particles. Both advantages may stimulate further investigations of relativistic quantum phenomena by transport measurement, spectral measurement etc. under extremely high magnetic field and ultra-low temperature in condensed matter systems.

Conclusions and outlook
This review mainly involves three peculiar features in low-dimensional superconducting and topological materials enriched by disorder effects, (i) 2D disordered Ising superconductivity mechanism in inversion-symmetric or asymmetric system; (ii) quantum Griffths singularity induced by quenched disorder in superconductor-metal transition; (iii) discrete scale invariance of quasi-bound states in topological materials. These anomalous features are beyond the traditional paradigms, provide clues for underlying microscopic processes, and reveal the profound characteristics of disordered quantum systems. Many important issues remain unsolved. For example, the mechanism for exotic anomalous metal phase associated with the superconductor-insulator transition is still under debate. Moreover, the systematic renormalization group treatment of supercritcal collapse for massless Dirac particles under magnetic field remains to be carried out. Future theoretical and experimental endeavors in these directions can largely enrich our knowledge of this exotic quantum garden, which may generate fruitful quantum technology.