![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
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
![](/cms/asset/e3de4d86-b483-4095-aee5-1a724961398b/tapx_a_1884133_uf0001_oc.jpg)
1. 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 [Citation1–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 [Citation8,Citation9]. 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 [Citation10–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 [Citation13–15]. Also, the prominent fluctuation and spin-orbit coupling have been proven indispensable in two-dimensional (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 [Citation1,Citation16–25], the emergence phase of anomalous quantum metal [Citation2,Citation26], the quantum Griffiths singularity (QGS) [Citation27] 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 [Citation5–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 [Citation28–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 [Citation29,Citation33]. 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.
2. 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 [Citation34]. 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 [Citation35]. 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 [Citation36]. Meanwhile, A. A. Abrikosov and L. Gor’kov studied alloys containing paramagnetic impurities based on the quantum field theory [Citation37], and had shown that the presence of paramagnetic impurities can diminish the s-wave pairing superconductivity. The dependence of on the paramagnetic impurity concentration can be derived via the Gor’kov Green function [Citation38–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 [Citation41–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 [Citation44].
Nevertheless, the WHH theory mainly focus on the 3D superconductors, which is not suitable to describe low dimensional superconductors [Citation16–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 Citation35. Klemm-Luther-Beasley (KLB) had firstly proposed the microscopic theory applicable to layered superconductors in the dirty limit with strong spin-orbit scattering [Citation45,Citation46].
Afterward, with the advancements of film-growth techniques, various types of high-quality 2D crystalline superconducting films have been fabricated [Citation13–15]. One striking discovery is the Ising pairing mechanism with pronounced enhancement of in-plane upper critical field beyond the Pauli limit, such as ion-gated
, monolayer
. However, this enhancement of in-plane
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 [Citation16]. In
and
, the atomic spin-orbital 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 [Citation1,Citation16–23]. However, in addition to the inversion-asymmetric Ising pairing, certain 2D inversion-symmetric materials with spin splitting around the
point due to intrinsic SOI are very recently found capable of superconductivity [Citation47] [see ]. Systematic models are needed to investigate the pair-breaking mechanism and in-plane
in both the inversion-symmetric and inversion-asymmetric Ising superconductors, taking into account of the combination effects of spin-orbit interaction and impurity scattering.
Figure 1. Mechanisms for an enhanced in-plane upper critical field. (a) Type-I Ising superconductivity: pairing of electrons on opposite spin split valleys. Here only one pair of electron pockets centered on and
points are highlighted. (b) Type-II Ising superconductivity: pairing of charge carriers on orbits around Γ point with their spins aligned in the out-of-plane orientation. Hole bands are illustrated as an example. The red and blue circles indicate two energetically degenerate bands with opposite spin orientations, each of which have spin split counterparts below the Fermi level (indicated by the dashed circles). Reproduced from Ref.[Citation47]
![Figure 1. Mechanisms for an enhanced in-plane upper critical field. (a) Type-I Ising superconductivity: pairing of electrons on opposite spin split valleys. Here only one pair of electron pockets centered on K and K ′ points are highlighted. (b) Type-II Ising superconductivity: pairing of charge carriers on orbits around Γ point with their spins aligned in the out-of-plane orientation. Hole bands are illustrated as an example. The red and blue circles indicate two energetically degenerate bands with opposite spin orientations, each of which have spin split counterparts below the Fermi level (indicated by the dashed circles). Reproduced from Ref.[Citation47]](/cms/asset/ef9d2276-22f2-48c2-a463-0aa992bab9b1/tapx_a_1884133_f0001_oc.jpg)
2.1. 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.[Citation48], 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 [Citation44,Citation49]:
where the system has translation symmetry after the impurity average and the integral kernel function is assumed to be isotropic.
is the anomalous Green’s function in momentum space.
2.1.1. 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 )], the intrinsic SOI serves as effective out-of-plane magnetic field which takes opposite value at valley and
, 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
and
(
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
, the Hamiltonian reads [Citation50]:
Figure 2. (Colour online) (a) Brillouin zone of ultrathin Pb film (inversion-asymmetric Ising pairing), with double non-degenerate Fermi surfaces (FSs) at each valley (green dashed circles). Inter-valley Cooper pairing form between electrons with the same color and opposite momentum. (b) Valley structure of Pb film in the vicinity of . Each valley has two electron pockets, and Zeeman-type SOI polarizes electron spin oppositely around
and
. (c-d)Theoretical fitting of
as a function of
for 6-ML Pb film. (c) shows good agreement with the microscopic theory of Zeeman-protected superconductivity [black solid line, .EquationEq. (3
(3)
(3) )] and yielding an effective Zeeman-type SOI
meV. (d) 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
meV and effective Rashba-type SOI
meV. The olive line indicates a special case with only effective Rashba-type SOI (
meV,
meV), which is far below the experimental results. Reproduced from Ref. [Citation21]
![Figure 2. (Colour online) (a) Brillouin zone of ultrathin Pb film (inversion-asymmetric Ising pairing), with double non-degenerate Fermi surfaces (FSs) at each valley (green dashed circles). Inter-valley Cooper pairing form between electrons with the same color and opposite momentum. (b) Valley structure of Pb film in the vicinity of EF. Each valley has two electron pockets, and Zeeman-type SOI polarizes electron spin oppositely around K and K ′. (c-d)Theoretical fitting of Bc/Bp as a function of T/Tc for 6-ML Pb film. (c) shows good agreement with the microscopic theory of Zeeman-protected superconductivity [black solid line, .EquationEq. (3(3) lnTcT=μB2B2βSO˜2+μB2B2Reψ12+iβSO˜2+μB2B22πkBT−ψ12.(3) )] and yielding an effective Zeeman-type SOI βSO˜=3.01 meV. (d) 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 βSO˜=3.16 meV and effective Rashba-type SOI αRkF˜=0.22 meV. The olive line indicates a special case with only effective Rashba-type SOI (αRkF˜=30meV, βSO˜=0 meV), which is far below the experimental results. Reproduced from Ref. [Citation21]](/cms/asset/2e873f36-c56d-428c-ac00-66c0f39e4b57/tapx_a_1884133_f0002_oc.jpg)
where is the effective mass,
and
denote the real spin subspace and valley subspace, respectively. Then the relation between critical field
and temperature
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 with
denoting the disorder concentration. Here, one can introduce effective Zeeman-type SOI parameter
, and the upper critical field
for the Zeeman-type SOI protected superconductivity in the dirty limit reads [Citation21,Citation48]:
If taking both the spin-independent scattering and the spin-flipping scattering into account, the characteristic time satisfies , and
,
,
denote the mean free time for total scattering, spin-independent scattering and spin-flipping scattering, respectively [Citation21]. Besides, Rashba-type SOI tends to polarize the spin to the in-plane direction, which weakens Zeeman-protection mechanism and changes the
. The Hamiltonian including the Zeeman-type SOI and Rashba-type SOI simultaneously at the
and
points are [Citation21]:
Likewise, one can define the effective Rashba SOI and effective Zeeman-type SOI parameter
, the relationship between
and
can be given as well. This microscopic calculation is consistent with previous phenomenological formula [Citation16] after the setting the disorder renormalized value of related parameters.
2.1.2. Type-II (inversion-symmetric) Ising superconductivity
Different from inversion-asymmetric Ising superconductors, like ), the energy valley locates at zero-momentum point in inversion-symmetric Ising superconductors. The intrinsic SOI can split the 4-fold degenerate
level into two doubly-degenerate levels, opening a gap at the
point [Citation14,Citation51]. The lower level crosses
at two different Fermi wavevectors form two different-shaped FSs [Citation14], each of which holds two states [see )]. If the radius of the FS is small, the eigenstates of the lower energy level can be approximated by
,
. The SOI at valley
can be viewed as an out-of-plane magnetic field which takes opposite value
and
on different orbits
and
, respectively. Thus, the system has TRS at
and electrons with opposite momenta and spins on the same FS can form s-wave Cooper pairs. These pairs with out-of-plane spin-orbit locking give rise to the large in-plane
. Based on previous study [Citation52], the inversion-symmetric system can be represented by a four band model with basis (
,
,
,
) with external in-plane magnetic field
,
Figure 3. Mechanisms for an enhanced in-plane upper critical field and superconducting properties of trilayer stanene. (a) Two doubly-degenerate FSs (green dashed circle) around the point in the Brillouin zone of stanene (inversion-symmetric Ising). The arrows on the FSs denote spin directions along
direction,
in blue and
in orange. Each FS is doubly degenerate, consisting of two different
orbitals and out-of-plane spin directions. The pairing happens between
and
on the same FS if
is small. (b) 4-fold degenerate
level of stanene splits to two doubly-degenerate levels due to SOI. The lower level forms two FSs at
(blue) and
(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
substrate thicknesses. For example,
refers to a trilayer stanene grown on top of 6 layers of
. The ratio of the in plane upper critical fields. i.e. the magnetic fields at which the sample resistance becomes
of the normal state resistance at a given temperature, to the Pauli limit field
are plotted as circular symbols. Solid/dashed curves are theoretical fits by .EquationEq. (6)
(6)
(6) .
is the zero-field transition temperature.
is the intrinsic SOC strength renormalized by disorder.
denotes the renormalized Rashba SOC strength. Reproduced from Ref. [Citation47]
![Figure 3. Mechanisms for an enhanced in-plane upper critical field and superconducting properties of trilayer stanene. (a) Two doubly-degenerate FSs (green dashed circle) around the Γ point in the Brillouin zone of stanene (inversion-symmetric Ising). The arrows on the FSs denote spin directions along ±z direction, k1 in blue and k2 in orange. Each FS is doubly degenerate, consisting of two different Px,y+ orbitals and out-of-plane spin directions. The pairing happens between Px−iy,↑+ and Px+iy,↓+ on the same FS if kr is small. (b) 4-fold degenerate Px,y+ level of stanene splits to two doubly-degenerate levels due to SOI. The lower level forms two FSs at k1 (blue) and k2 (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.86Tc,0 are plotted as circular symbols. Solid/dashed curves are theoretical fits by .EquationEq. (6)(6) lnTcTc,0=−μB2B2βso∗2+μB2B2Reψ12+iβso∗2+μB2B22πkBTc−ψ12,(6) . Tc,0 is the zero-field transition temperature. βSO∗ is the intrinsic SOC strength renormalized by disorder. αkF denotes the renormalized Rashba SOC strength. Reproduced from Ref. [Citation47]](/cms/asset/f64e0b30-08e4-431c-a417-4cc58652b862/tapx_a_1884133_f0003_oc.jpg)
Here is the effective Bohr magneton and
,
,
,
are fitting parameters. Considering the situation of lower band
crossing
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 [Citation48]:
where is the disorder renormalized SOC strength,
,
,
denote the Fermi momentum, the digamma function and the scattering time, respectively. Near
, this result is consistent with the result from 2D Ginzburg-Landau (GL) theory [Citation53] and the KLB theory [Citation45,Citation46].
2.2. Recent experimental observations of Ising superconductivity
With the advancements of the nanotechnology including molecular beam epitaxy, pulsed laser deposition etc, various types of 2D superconducting materials have been fabricated [Citation13–15]. The high-quality 2D superconductors give rise to many novel physical phenomena due to the prominent influence of intrinsic spin-orbit coupling and quantum fluctuations. Therefore, these systems provide an insightful platform to explore the essence of superconductivity [Citation4,Citation27,Citation54–63].
2.2.1. Type-I Ising superconductors
It has been a significant discovery that ion gated and
can be superconducting [Citation64,Citation65]. Unlike the bulk transition metal dichalcogenides (TMDs) that preserve inversion symmetry, monolayer
and
are inversion-asymmetric superconductors [Citation66]. 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[Citation21]. 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
curves of a 6-ML sample are measured in a pulsed parallel magnetic field up to
T. The superconductivity survives under a high parallel magnetic field of
T at
K. The experimental data of
versus
can be quantitatively fitted by EquationEq. (3)
(3)
(3) . The fitting procedure gives effective Zeeman-type SOI
meV, as shown in ). 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, ) displays a set of theoretical curves with a fixed effective Zeeman-type SOI
meV and increasing effective Rashba-type SOI
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.
2.2.2. Type-II Ising superconductor
Recent experiment has identified a pronounced enhancement of in-plane critical field at low temperature in epitaxial few layer
-Sn(111) [Citation14,Citation15,Citation47]. Distinct from the previous system, few-layer stanene has no
mirror symmetry and is centro-symmetric in the free standing case [Citation51]. Stanene hosts bands around
-point instead of the
and
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 , the experimental data matches well with EquationEq. (6(6)
(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
. Such pairing mechanisms behind EquationEq. (6
(6)
(6) ), without the participation of inversion symmetry breaking, are termed type-II (inversion-symmetric) Ising superconductivity. In fact, the upturn of
in low-T result from frozen spin orientation at
. 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
. ) show the temperature dependence of the in-plane upper critical fields in few-layer stanene samples. Comparing to
, the smoothening in
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 )). Spins of the
and
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 (). Since stronger inversion symmetry breaking attribute to the missing up-turn feature, a penta-layer stanene should experience a weaker Rashba effect, giving rise to an apparent enhancement of
at low
. 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 [Citation67].
2.3. Discussion
There may be multiple types of SOI, including the Zeeman type SOI and the Rashba type SOI both may affect the in-plane Citation16,[Citation17], 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
, including the quantum Griffiths singularity along the anomalous phase boundary [Citation68]. The influence of large intrinsic SOI on the quantum fluctuation needs to be considered with the second order corrections [Citation68,Citation69], which is beyond the WHH framework.
3. 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 [Citation1,Citation27,Citation70–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 in which the control parameter is non-thermal, such as pressure, magnetic field, or chemical composition [Citation81]. When approaching the critical point of continuous phase transition, the characteristic length scales become divergent [Citation82]. Upon approaching the critical point, the correlation length
and the correlation time diverge as [Citation83,Citation84]:
where is dimensionless parameter and
is the critical parameter,
is the correlation length critical exponent and
is dynamical critical exponent. Power-law singularities in the length scales at the critical point generically lead to power-law singularities in observables [Citation85]. 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.
3.1. Scaling theory for continuous phase transition
The first unification of phase transitions theories was achieved by L. D. Landau [Citation86,Citation87]. He proposed a mean-field order parameter describing phase transition behavior at the transition point,
represents a disordered state while
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 [Citation84]. 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 [Citation88].
The power law behavior at critical point implied the scaling theory. Widom firstly put forward as the scaling hypothesis [Citation89]. 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 [Citation82]:
where denotes the distance from the critical point,
is dimensionless measure of external magnetic field,
is the space dimensionality,
and
are two critical exponents and
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 [Citation89]. 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 [Citation90], and the related concept of universality classes can be well explained within the framework of the modern renormalization group theory [Citation91].
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 [Citation84]. Consider the fluctuations in the imaginary axis,
, if we scale length according to
, the imaginary time need scale like
, and temperature regarded as an energy scales like an inverse time
. Consequently, the scaling form of the free energy density at a quantum critical point reads [Citation83]:
Here, is the absolute temperature,
is the dimensionless parameter (i.e. the dimensionless magnetic field), and
is dynamics exponent. Comparing EquationEq. (9)
(9)
(9) to EquationEq. (8
(8)
(8) ), one can directly find that a quantum phase transition in
dimensions is equivalent to a classical transition in
spatial dimensions [Citation92].
3.2. 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 [Citation93]. How does quenched disorder influence phase transitions and critical points? Harris [Citation94] derived a criterion for the stability of a given critical behavior with respect to quenched disorder. For short range correlated disorder, if space dimensionality and clean correlation length exponent
fulfill the Harris criterion
, weak disorder does not affect the critical behavior of the system, corresponding to the clean limit [Citation74,Citation94–96]. On the contrary, if
, 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 [Citation97–100], similarly to the features of random transverse field Ising model [Citation101]. This kind of singularities are called the Griffiths singularities, because R. B. Griffiths firstly proved its existence in 1969 [Citation102]. T. Vojta has provided the whole classification of disorder effect on the scaling behaviours, based on the analysis of the dimension of range regions [Citation84]. 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 [Citation98].
(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 fulfill the inequality
), 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 and spatial correlation length satisfy a activated (exponential) scaling
, distinctive from the relation
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.
3.2.1. 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 long-range order and thus reduces from its clean bulk value
. 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 [Citation84].
Long-range disorder correlations can immensely increase the rare region effects. This rare region effect was firstly well depicted in the McCoy-Wu model [Citation103], 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 are random positive interactions, the random transverse fields
cause the quantum fluctuations. At temperature T, such a quantum system in d-dimensions is equivalent to a
dimensional classical Ising model with the couplings independent of the ”time” direction in which it has extent
. The behavior of the random transverse field Ising model was not fully understood until D. S. Fisher [Citation101,Citation104,Citation105] applied a strong-disorder renormalization group method to the problem.
The main idea of the strong-disorder renormalization group is to take the strongest coupling in the system and coupling to the rest of the system is treated perturbatively [Citation91]. 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
of the clusters and effective bonds scale with the logarithm of the energy scale
with the exponent , and
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
decreases exponentially
Citation84. Moreover, if a system is pushed slightly off critical with a distance
, the correlation lengths satisfies
with
. In the quantum Griffiths phase, the activated scaling law results in a typical energy scale
, which obeys
. 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
satisfying [Citation105]:
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
, (
is parameter of system,
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 [Citation19,Citation106]. Moreover, the activated scaling analysis can be applied directly to elucidate the quantum Griffiths phase [Citation107–109].
3.3. 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 [Citation110–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 of superconductor insulator/metal transition is considered equal to the crossing point of sheet resistance
curves at different temperatures [Citation74,Citation113]. Recently, it is striking that two crossing points of
curves were observed in superconductor metal transitions (SMTs) in 2D superconductor thin films [Citation76,Citation114]. 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 [Citation27]. In , the magneto-resistance
curve of the sample under different temperatures do not cross at one point [Citation110,Citation111,Citation115], but forming a crossing line. Meanwhile, there is an upturn of
at
beyond the WHH theory. Utilizing the finite-size scaling (FSS) method for a bundle of
curves around the multiple crossing points,
Figure 4. Quantum Griffiths singularity in Ga film. (a) Magneto-resistance isotherms. The magnetic field dependence of at various temperatures ranging from 0.025 K to 2.80 K, with the same sweep direction and the same sweep rate, exhibits an SMT region. Inset: the critical magnetic fields
(T). Crossing points of
(B) curves at every two adjacent temperatures are denoted as black dots on the transition boundary; the red stars come from the temperature plateaus on
(T) curves at given magnetic fields. The blue line is the fitting curve using WHH theory and the dashed blue line (below 1.00 K) are guides to the eye. (b) The activated quantum scaling behavior: exponent
as a function of magnetic field
. When approaching the zero-temperature limit,
rapidly increases with no sign of saturation. The magenta line shows a fitting based on the activated scaling law. (c) Sketch of the
phase diagram of 2D SMT with quenched disorder.
is the upper critical field and would be the QCP for a clean system. When temperature decreases to
, quenched disorder overtakes thermal fluctuation and gives rise to a vortex glass–like phase. Rare regions of superconducting islands emerge and result in activated scaling behavior around such an infinite-randomness QCP
. Reproduced from Ref.[Citation27]
![Figure 4. Quantum Griffiths singularity in Ga film. (a) Magneto-resistance isotherms. The magnetic field dependence of RS at various temperatures ranging from 0.025 K to 2.80 K, with the same sweep direction and the same sweep rate, exhibits an SMT region. Inset: the critical magnetic fields Bc(T). Crossing points of RS (B) curves at every two adjacent temperatures are denoted as black dots on the transition boundary; the red stars come from the temperature plateaus on RS (T) curves at given magnetic fields. The blue line is the fitting curve using WHH theory and the dashed blue line (below 1.00 K) are guides to the eye. (b) The activated quantum scaling behavior: exponent zv as a function of magnetic field B. When approaching the zero-temperature limit, zv rapidly increases with no sign of saturation. The magenta line shows a fitting based on the activated scaling law. (c) Sketch of the B−T phase diagram of 2D SMT with quenched disorder. B2 is the upper critical field and would be the QCP for a clean system. When temperature decreases to T<TM, quenched disorder overtakes thermal fluctuation and gives rise to a vortex glass–like phase. Rare regions of superconducting islands emerge and result in activated scaling behavior around such an infinite-randomness QCP Bc∗. Reproduced from Ref.[Citation27]](/cms/asset/508582aa-fd5d-45ba-a258-97d34ae9d83d/tapx_a_1884133_f0004_oc.jpg)
It turns out the critical exponents exhibit a power-law divergence when
at
,
which matches well with QGS feature [Citation105]. In addition, ) sketch the phase diagram of 2D SMT with quenched disorder in Ga film. The anomalous upturn of upper critical field as
appears. In the conventional theory of homogeneous type II superconductors, the Abrikosov vortex lattice forms in the region of
, characterized by a continuous phase transition at the upper critical field
[Citation38,Citation116]. In clean systems, thermal fluctuations give rise to phonon-like excitations of the vortex lattice. At
regime, thermal fluctuations become sufficiently large, the system is expected to melt into the vortex liquid phase [Citation117]. While at
, the effect of quenched disorder overtakes the thermal fluctuation effect and produces a vortex glass–like phase on length scales greater than characteristic length of disorder [Citation118]. This vortex glass–like phase consists of spatially separated superconducting regions, where the local upper critical field exceeds the average mean-field value
. 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
.
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 interface and the dynamic critical exponent approaching the quantum critical point also exhibits QGS [Citation119]. Subsequently, QGS has been reported in a wide range of crystalline 2D superconductors, such as
Citation19,
and
Citation120. Different from traditional superconductor Ga film with weak SOC,
monolayer film is Ising superconductor with strong SOC. The discovery of the QGS in the
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 [Citation121–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 ultra-thin crystalline Pb films at ultra-low temperatures [Citation68]. Usually, the SMT shows monotonic phase boundary separating the regions of and
. However,
curves in 4-ML Pb films exhibit a remarkable re-entrant behavior at lower temperatures in ). When applying the magnetic field slightly large than 3.3 T, the sheet resistance firstly decreases with decreasing temperature, reaches the minimum at
and then rises at a lower temperature. This re-entrant behavior, including the fact that
increases with increasing magnetic field, can be quantitatively reproduced after taking the superconducting fluctuation and the influence of large SOI into account [Citation125,Citation126]. Also, the
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
K as shown in ), 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 [Citation126]. 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
matches well with QGS behavior [Citation97,Citation99].
Figure 5. Anomalous quantum Griffiths singularity of SMT in Pb film. (a) curves at various magnetic fields show SMT and re-entrant behavior. (b) The crossing points from neighboring
curves and the onset
from
curves determine the phase boundary of SMT. The phase boundary abnormally bends down at lower temperatures. (c) Scaling analysis indicates anomalous QGS behavior in 3.5-ML Pb film. The critical exponent diverges as
with temperature approaching zero and magnetic field tending to
. (d) Schematic of the
phase diagram of 2D SMT with pronounced fluctuation effect.
is the value of upper critical field with mean field theory. Under the influence of superconducting fluctuations with strong SOI, the mean field phase boundary (blue dashed line) buckles outward to the solid red line. The anomalous phase boundary exists in the regime
. When approaching the real infinite-randomness quantum critical point
along this anomalous phase boundary, the system exhibits anomalous QGS behavior. Reproduced from Ref.[Citation124]
![Figure 5. Anomalous quantum Griffiths singularity of SMT in Pb film. (a)Rs(T) curves at various magnetic fields show SMT and re-entrant behavior. (b) The crossing points from neighboring Rs(B) curves and the onset Tc from Rs(T) curves determine the phase boundary of SMT. The phase boundary abnormally bends down at lower temperatures. (c) Scaling analysis indicates anomalous QGS behavior in 3.5-ML Pb film. The critical exponent diverges as zν∝B−Bc∗−0.6 with temperature approaching zero and magnetic field tending to Bc∗. (d) Schematic of the B−T phase diagram of 2D SMT with pronounced fluctuation effect. B2 is the value of upper critical field with mean field theory. Under the influence of superconducting fluctuations with strong SOI, the mean field phase boundary (blue dashed line) buckles outward to the solid red line. The anomalous phase boundary exists in the regime Bc∗,Bc′. When approaching the real infinite-randomness quantum critical point Bc∗ along this anomalous phase boundary, the system exhibits anomalous QGS behavior. Reproduced from Ref.[Citation124]](/cms/asset/39a689e4-a42d-42ab-b921-338c71b5cc27/tapx_a_1884133_f0005_oc.jpg)
The QGS phase can be directly verified by the activated scaling analysis [Citation108,Citation109]. In practical analysis, certain irrelevant scaling variables influence the crossing point [Citation127], and needs to be considered. The activated scaling with irrelevant corrections reads: . Here
and
are critical exponent and
is related to the irrelevant parameter. Both the relevant and irrelevant parameters are associated with the activated scaling form with effective length scale
. 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 [Citation108,Citation109].
3.4. 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
, 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.
4. 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 [Citation7,Citation128]. Recently, the discrete scale invariance (DSI) behavior in Dirac materials has attracted attention in related subfields of physics [Citation28,Citation129–132]. Especially, a new magneto-resistance oscillation beyond the quantum limit [Citation28,Citation30,Citation133] has been found in several high-quality topological materials. Different from the previously known periodic quantum oscillations (Shubnikov–de Haas oscillations in periodic [Citation134,Citation135], Aharonov-Bohm and Altshuler-Aronov-Spivak effects induced quantum oscillations in B periodic [Citation136]), this novel oscillation is
periodic which manifests the appearance of DSI in solid-state system [Citation28,Citation30,Citation133]. Such a peculiar DSI feature is considered to be universal in Dirac materials with Coulomb attraction [Citation28,Citation29], which may be closely related to the quasi-bound states formed by massless Dirac fermions and the long-pursue atomic-collapse phenomenon [Citation28,Citation131,Citation137]. Thus, it is desirable to explore the underlying mechanisms of this novel phenomenon.
4.1. 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 has the property that
, where
is a continuous parameter and
is the scaling dimension. When considering the quantization effect, the continuous scale invariance will be spontaneously broken down to the DSI, where
takes on discrete values [Citation138]. Its applicability to bound-state problems in quantum mechanics was exploited by Vitaly Efimov, who found the
potential suited to Schrödinger equation and confirmed the DSI energy spectrum could occur in 3-body problem in the limit of resonant scattering [Citation139,Citation140] [see ]. Analogy to Efimov 3-body problem, topological semimetals hold Weyl equation that the energy and momentum satisfy
and DSI is also expected to achieve with
potential [Citation141].
Figure 6. System with discrete-scale invariance (a) M. C. Escher’s painting: Circle LimitIII represents a famous fractal structure: Poincarë disk model. (b) Mandelbrot set (c) Efimov 3-body states
![Figure 6. System with discrete-scale invariance (a) M. C. Escher’s painting: Circle LimitIII represents a famous fractal structure: Poincarë disk model. (b) Mandelbrot set (c) Efimov 3-body states](/cms/asset/958268ab-8331-4841-8dfa-f43b8985987f/tapx_a_1884133_f0006_oc.jpg)
Considering the Weyl fermions with Coulomb attraction near the Weyl point, the Hamiltonian read [Citation29]:
with potential . Here
denotes fine structure constant,
is the fermion velocity,
is the charge on the impurity and
is the dielectric constant. At
, then EquationEq. (15)
(15)
(15) can be solved analytically [Citation142,Citation143]. In the case of attractive potential
(
is a dimensionless quantum numbers of angular momentum), the energy levels of quasi-bound states exhibit DSI [Citation129,Citation130,Citation144,Citation145],
At , the magnetic field introduces a new characteristic length
that breaks the DSI at large
. The problem with the finite magnetic field is no longer exactly solvable but can be deal with based on Wentzel-Kramers Brillouin (WKB) approximation [Citation29]. In this quasi-classical picture, the energy change occurs when
with
. And since
satisfy EquationEq. (16
(16)
(16) ), the characteristic values of
meet
. As shown in , 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
periodic transport property [Citation29].
Figure 7. The DSI properties of quasi-bound states. (a) Without magnetic field, quasi-bound states form around the charge center. The Klein tunneling with electron-hole excitation gives rise to finite widths to the quasi-bound state spectra. Due to the small FS, screening is relatively weak. (b) The local density of states (LDOS) at certain location as a function of the dimensionless energy showing the broad quasi-bound states. The exact analytical solutions (solid line) and the WKB-type solutions (dashed line) are shown. (c) With a magnetic field, the potential changes in the large region, and the quasi-bound states move to the Fermi level at certain field values
. (d) The evolution of the quasi-energy spectrum versus the dimensionless inverse magnetic length
. 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
that form a geometric series. Reproduced from Ref.[Citation29]
![Figure 7. The DSI properties of quasi-bound states. (a) Without magnetic field, quasi-bound states form around the charge center. The Klein tunneling with electron-hole excitation gives rise to finite widths to the quasi-bound state spectra. Due to the small FS, screening is relatively weak. (b) The local density of states (LDOS) at certain location as a function of the dimensionless energy showing the broad quasi-bound states. The exact analytical solutions (solid line) and the WKB-type solutions (dashed line) are shown. (c) With a magnetic field, the potential changes in the large r region, and the quasi-bound states move to the Fermi level at certain field values Bn. (d) The evolution of the quasi-energy spectrum versus the dimensionless inverse magnetic length lB. 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 Bn that form a geometric series. Reproduced from Ref.[Citation29]](/cms/asset/5b201709-6b66-492b-8865-62dd423e0fd6/tapx_a_1884133_f0007_oc.jpg)
4.2. 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 within Thomas Fermi approximation, where the screening length
obeys the relation
, with
and
denoting the carrier density, the chemical potential, and the degeneracy, respectively. The results point out the introduction of the length parameter
in
causes the breaking of DSI and when
, DSI is restored. Moreover, the finite mass term also introduce a new length scale and it may break the DSI down. Solving the radial equation with mass term
, the spectrum still satisfies EquationEq. (16)
(16)
(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
as an example [Citation146], the Zeeman spin splitting can influence the lowest Landau level [Citation29]. When magnetic field is applied perpendicularly on the layered
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
. When magnetic field is applied parallel on the layered
. 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.
4.3. Ln(B) periodic oscillation
The periodic transport property of system can be captured by the T-matrix approximation for the Kubo formula [Citation29]. According to the classical treatment of Bastin et al. [Citation29,Citation147], the formula for the transverse magneto-conductivity
at zero temperature is [Citation29]:
where .
is the density of short-range scatterers,
is the density of Coulomb scatterers,
are the quasi-bound state energies as calculated and shown above, and
is the width determined by overlap with the lowest Landau level.
is the effective length along the magnetic field,
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 [Citation29]:
with fitting parameters denoting an effective broadening factor and
denoting a characteristic magnetic field. The longitudinal magneto-conductivity can be obtained by relation
. The conductivity in EquationEq. (18)
(18)
(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
periodic MR oscillations.
This new type of log-periodic MR oscillations was observed beyond the quantum limit in ultra-quantum topological materials . ) shows the extract oscillations in different
samples (s6, s7, and s9) [Citation28]. Using formula EquationEq. (18)
(18)
(18) can quantitatively fit these
periodic oscillations (black curves). To better presentation, the magnetic fields (
) corresponding to MR peaks/dips were extracted, and the dependence of
on
are plotted in ). Here, the MR peaks and dips were indexed by integers
and half-integers
. The well-defined linear behavior for the
vs
curve indicates the clearest feature of DSI. This feature can be understood vividly like ).
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
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.
Figure 8. DSI in ultra-quantum . (a) Index plot for the log-periodic oscillations. (b) Quantitative fitting (black curves) of the log-periodic oscillations in
. arb.u, arbitrary units. (c) Schematic of the two-body quasi-bound states composed of a Dirac-type massless hole and charged center via Coulomb attraction. Reproduced from Ref.[Citation28]
![Figure 8. DSI in ultra-quantum ZrTe5. (a) Index plot for the log-periodic oscillations. (b) Quantitative fitting (black curves) of the log-periodic oscillations in ZrTe5. arb.u, arbitrary units. (c) Schematic of the two-body quasi-bound states composed of a Dirac-type massless hole and charged center via Coulomb attraction. Reproduced from Ref.[Citation28]](/cms/asset/b5e283e0-d16f-402d-8a21-d768b4185808/tapx_a_1884133_f0008_oc.jpg)
Such peculiar work also inspires people to explore the log-periodic quantum oscillations and the DSI in other physical observables and topological materials. Ref[Citation30]. reported the logB-periodic quantum oscillations in both the longitudinal MR and the Hall resistance
of
crystals, which identify the
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 [Citation148] found similar logB-periodic longitudinal MR and the Hall resistance oscillations to those observed in
. The small gap does not change the DSI feature in this system, which is consistent with previous theoretical analysis [Citation33]. Moreover, the MR quantum oscillations
Citation133 is considered as a possible candidate for logB-periodic oscillation, revealing the universality of DSI property in Dirac/Weyl materials with supercritical Coulomb attraction.
4.4. 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 [Citation129,Citation131]. 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.
5. 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.
Acknowledgements
We thank Jian Wang, Qi-Kun Xue, Robert Joynt, Ding Zhang, Jie Xiong, Faxian Xiu, Xi Lin, Ziqiang Wang, Xu-cun Ma, Yi Liu, Huichao Wang, Hongchao Liu, Chao Yang and Ce Huang for continuous collaboration of related works in this review. This work is financially supported by the National Basic Research Program of China (Grants No. 2017YFA0303301, No. 2017YFA0304600, and No. 2019YFA0308403) and the National Natural Science Foundation of China (Grants No. 11534001, No. 11674028, No. 11822407, and No. 11921005).
* Haiwen Liu: [email protected]
† X. C. Xie: [email protected]
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Saito Y, Kasahara Y, Jianting Y, et al. Metallic ground state in an ion-gated two-dimensional superconductor. Science. 2015;350:409–32.
- Kapitulnik A, Kivelson SA, Spivak B. Colloquium: anomalous metals: failed superconductors. Rev Mod Phys. 2019 Jan;91:011002.
- Reyren N, Thiel S, Caviglia AD, et al. Superconducting interfaces between insulating oxides. Science. 2007;317:1196–1199.
- Zhang T, Cheng P, Wenjuan L, et al. Superconductivity in one-atomic-layer metal films grown on si(111). Nat Phys. 2010;6:104–108.
- Weng H, Xi D, Zhong F. Topological semimetals predicted from first-principles calculations. J Phys. 2016;28:303001.
- Wang S, Lin B-C, Wang A-Q, et al. Quantum transport in dirac and weyl semimetals: a review. Adv Phys X. 2017;2:518–544.
- Ren Y, Qiao Z, Niu Q. Topological phases in two-dimensional materials: a review. Rep Prog Phys. 2016;79:066501.
- Hohenberg PC. Existence of long-range order in one and two dimensions. Phys Rev. 1967 Jun;158:383–386.
- Mermin ND, Wagner H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic heisenberg models. Phys Rev Lett. 1966;17:1133–1136.
- Kosterlitz JM, Thouless DJ. Ordering, metastability and phase transitions in two-dimensional systems. J Phys C. 1973;6:1181–1203.
- Berezinsky VL. Destruction of long range order in one-dimensional and two-dimensional systems having a continuous symmetry group. i. classical systems. J Exp Theor Phys. 1970;32:493–500.
- Goldman A. Percolation, localization, and superconductivity. Vol. 109. Springer Science & Business Media, Berlin; 2013.
- Uchihashi T. Two-dimensional superconductors with atomic-scale thickness. Superconductor SciTechnol. 2017;30:013002.
- Liao M, Zang Y, Guan Z, et al. Superconductivity in few-layer stanene. Nat Phys. 2018;14:344–348.
- Zang Y, Jiang T, Gong Y, et al. Realizing an epitaxial decorated stanene with an insulating bandgap. Adv Funct Mater. 2018;28:1802723.
- Jianming L, Zheliuk O, Leermakers I, et al. Evidence for two-dimensional ising superconductivity in gated mos2. Science. 2015;350:1353–1357.
- Bahramy NY, M. S. Kohama Y, Kasahara YJ, et al. Superconductivity protected by spin–valley locking in ion-gated mos2. Nat Phys. 2016;12:144–149.
- Xiaoxiang X, Wang Z, Zhao W, et al. Ising pairing in superconducting nbse2 atomic layers. Nat Phys. 2016;12:139–143.
- Xing Y, Zhao K, Shan P, et al. Ising superconductivity and quantum phase transition in macro-size monolayer nbse2. Nano Lett. 2017;17:6802–6807.
- Saito Y, Nojima T, Iwasa Y. Highly crystalline 2d superconductors. Nat Rev Mater. 2017;2:1–18.
- Liu Y, Wang Z, Zhang X, et al. Interface-induced zeeman-protected superconductivity in ultrathin crystalline lead films. Phys Rev X. 2018;8:021002.
- Jianming L, Zheliuk O, Chen Q, et al. Full superconducting dome of strong ising protection in gated monolayer ws2. Proc Natl Acad Sci U S A. 2018;115:3551–3556.
- De La Barrera SC, Sinko MR, Gopalan DP, et al. Tuning Ising superconductivity with layer and spin–orbit coupling in two-dimensional transition-metal dichalcogenides. Nat Commun. 2018;9:1427.
- Sergio C, Sinko MR, Gopalan DP, et al. Tuning ising superconductivity with layer and spin–orbit coupling in two-dimensional transition-metal dichalcogenides. Nat Commun. 2018;9:1–8.
- Cui J, Peiling L, Zhou J, et al. Transport evidence of asymmetric spin–orbit coupling in few-layer superconducting 1td-mote2. Nat Commun. 2019;10:2044.
- Yang C, Liu Y, Yang W, et al. Intermediate bosonic metallic state in the superconductor-insulator transition. Science. 2019;366:1505–1509.
- Xing Y, Zhang H, Hailong F, et al. Quantum griffiths singularity of superconductor-metal transition in ga thin films. Science. 2015;350:542–545.
- Wang H, Liu H, Yanan L, et al. Discovery of log-periodic oscillations in ultraquantum topological materials. Sci Adv. 2018;4:eaau5096.
- Liu H, Jiang H, Wang Z, et al. Discrete scale invariance in topological semimetals. arXiv preprint arXIv:1807.02459. 2018.
- Wang H, Liu Y, Liu Y, et al. Log-periodic quantum magneto-oscillations and discrete-scale invariance in topological material hfte5. Natl Sci Rev. 2019;6:914–920.
- Wang CM, Sun H-P, Hai-Zhou L, et al. 3d quantum hall effect of fermi arcs in topological semimetals. Phys Rev Lett. 2017 Sep;119:136806.
- Hailong L, Liu H, Jiang H, et al. 3d quantum hall effect and a global picture of edge states in weyl semimetals. Phys Rev Lett. 2020 Jul;125:036602.
- Liu H, Liu H, Joynt R, et al. Effect of coulomb screening on the discrete scale invariance of quasibound states in three-dimensional topological semimetals. Phys Rev B. 2019;100:195140.
- Bardeen J, Cooper LN, Schrieffer JR. Theory of superconductivity. Phys Rev. 1957;108:1175–1204.
- Cooper DEFLN. BCS: 50 years. 1 edition ed. World Scientific, Singapore; 2011.
- Anderson PW. Theory of dirty superconductors. J Phys Chem Solids. 1959;11:26–30.
- Abrikosov AA, Gorkov LP. On the theory of superconducting alloys. 1. the electrodynamics of alloys at absolute zero. Sov Phys JETP. 1959;8:1090–1098.
- Abrikosov AA, Gorkov LP. Contribution to the theory of superconducting alloys with paramagnetic impurities. Sov Phys JETP. 1960;12:1243.
- Maki K. The magnetic properties of superconducting alloys. Physics Physique Fizika. 1964;1:21.
- De Gennes PG. Behavior of dirty superconductors in high magnetic fields. Physik der kondensierten Materie. 1964;3:79–90.
- Clogston AM. Upper limit for the critical field in hard superconductors. Phys Rev Lett. 1962;9:266–267.
- Chandrasekhar BS. A note on the maximum critical field of high‐field superconductors. Appl Phys Lett. 1962;1:7–8.
- Maki K. Pauli paramagnetism and superconducting state. Prog Theoretical Phys. 1964;32:29–36.
- Werthamer NR, Helfand E, Hohenberg PC. Temperature and purity dependence of the superconducting critical field, h c 2 . iii. electron spin and spin-orbit effects. Phys Rev. 1966;147:295–302.
- Klemm RA, Luther A, Beasley MR. Theory of the upper critical field in layered superconductors. Phys Rev B. 1975;12:877–891.
- Tedrow PM, Meservey R. Critical magnetic field of very thin superconducting aluminum films. Phys Rev B. 1982;25:171–178.
- Falson J, Yong X, Liao M, et al. Type-ii ising pairing in few-layer stanene. Science. 2020;367:1454–1457.
- Liu H, Liu H, Zhang D, et al. Microscopic theory of in-plane critical field in two-dimensional ising superconducting systems. Phys Rev B. 2020;102:174510.
- Gorkov L, Abrikosov P, A A, et al. Methods of quantum field theory in statistical physics. Courier, Corporation, Massachusets; 2012.
- Yuan NFQ, Mak KF, Law KT. Possible topological superconducting phases of mos 2. Phys Rev Lett. 2014;113:097001.
- Yong X, Yan B, Zhang H, et al. Large-gap quantum spin hall insulators in tin films.. Phys Rev Lett. 2013;111:136804––1368045.
- Wang C, Lian B, Guo X, et al. Type-ii ising superconductivity in two-dimensional materials with spin-orbit coupling. Phys Rev Lett. 2019;123:126402.
- Tinkham M. Effect of fluxoid quantization on transitions of superconducting films. Phys Rev. 1963;129:2413–2422.
- Guo Y, Zhang Y, Bao X, et al. Superconductivity modulated by quantum size effects. Science. 2004;306:1915–1917.
- Ozer MM, Thompson JR, Weitering HH. Hard superconductivity of a soft metal in the quantum regime. Nat Phys. 2006;2:173–176.
- Eom D, Qin S, Chou MY, et al. Persistent superconductivity in ultrathin pb films: a scanning tunneling spectroscopy study. Phys Rev Lett. 2006;96:027005.
- Qin S, Kim J, Niu Q, et al. Superconductivity at the two-dimensional limit. Science. 2009;324:1314–1317.
- Uchihashi T, Mishra P, Aono M, et al. Macroscopic superconducting current through a silicon surface reconstruction with indium adatoms: si(111)-(√7 × √3)-in. Phys Rev Lett. 2011;107:207001.
- Yamada M, Hirahara T, Hasegawa S. Magnetoresistance measurements of a superconducting surface state of in-induced and pb-induced structures on si(111). Phys Rev Lett. 2013;110:237001.
- Brun C, Cren T, Cherkez V, et al. Remarkable effects of disorder on superconductivity of single atomic layers of lead on silicon. Nat Phys. 2014;10:444–450.
- Zhang H, Sun Y, Wei L, et al. Detection of a superconducting phase in a two-atom layer of hexagonal ga film grown on semiconducting gan(0001). Phys Rev Lett. 2015;114:107003.
- Wang Q, Zhi L, Zhang W, et al. Interface-induced high-temperature superconductivity in single unit-cell fese films on srtio3. Chin Phys Lett. 2012;29:037402.
- Zhang W, Sun Y, Zhang J, et al. Direct observation of high-temperature superconductivity in one-unit-cell fese films. Chin Phys Lett. 2014;31:017401.
- Akashi R, Bahramy MS, Arita R, et al. Superconducting dome in a gate tuned band insulator. Science. 2012;338:11193.
- Shimotani H, Taniguchi K, Matsumoto A, et al. Electric-field-induced superconductivity at 9.4 k in a layered transition metal disulphide mos2.. Appl Phys Lett. 2012;101:042603.
- Yuan NFQ, Zhou BT, He WY, et al. Ising superconductivity in transition metal dichalcogenides. 2016.
- Liu Y, Yong X, Sun J, et al. Type-ii ising superconductivity and anomalous metallic state in macro-size ambient-stable ultrathin crystalline films. Nano Lett. 2020 June;20:5728–5734.
- Liu Y, Wang Z, Shan P, et al. Anomalous quantum griffiths singularity in ultrathin crystalline lead films. Nat Commun. 2019;10:1–6.
- Varlamov AA, Galda A, Glatz A. Fluctuation spectroscopy: from rayleigh-jeans waves to abrikosov vortex clusters. Rev Mod Phys. 2018 Mar;90:015009.
- Mason N, Kapitulnik A. Dissipation effects on the superconductor-insulator transition in 2d superconductors. Phys Rev Lett. 1999;82:5341–5344.
- Sachdev S. Quantum Phase Transitions. Cambridge University Press, Cambridge; 2011.
- Christiansen C, Hernandez LM, Goldman AM. Evidence of collective charge behavior in the insulating state of ultrathin films of superconducting metals. Phys Rev Lett. 2002;88:037004.
- Del Maestro A, Rosenow B, Sachdev S. Theory of the pairbreaking superconductor-metal transition in nanowires. Ann Phys. 2009;324:523–583.
- Goldman AM. Superconductor-insulator transitions. Int J Modern Phys B. 2010;24:4081–4101.
- Eley S, Gopalakrishnan S, Goldbart PM, et al. Approaching zero-temperature metallic states in mesoscopic superconductor-normal-superconductor arrays. Nat Phys. 2012;8:59–62.
- Biscaras J, Bergeal N, Hurand S, et al. Multiple quantum criticality in a two-dimensional superconductor. Nat Mater. 2013;12:542–548.
- Shi X, Lin PV, Sasagawa T, et al. Two-stage magnetic-field-tuned superconductor–insulator transition in underdoped la 2−x sr x cuo 4. Nat Phys. 2014;10:437–443.
- Breznay N, Kapitulnik A. Particle-hole symmetry reveals failed superconductivity in the metallic phase of two-dimensional superconducting films. Sci Adv. 2017;3:9.
- Bottcher CGL, Nichele F, Kjaergaard M, et al. Superconducting, insulating and anomalous metallic regimes in a gated two-dimensional semiconductor–superconductor array. Nat Phys. 2018;14:1138–1144.
- Chen Z, Swartz A, Yoon H, et al. Carrier density and disorder tuned superconductor-metal transition in a two-dimensional electron system.. Nat Commun. 2018;9:4008.
- Sondhi SL, Girvin SM, Carini JP, et al. Continuous quantum phase transitions. Rev Mod Phys. 1997;69:315–333.
- Ma SK. Modern theory of critical phenomena. Routledge, Oxfordshire; 2018.
- Herbut I. A Modern Approach to Critical Phenomena. Cambridge University Press; 2007.
- Vojta T. Rare region effects at classical, quantum and nonequilibrium phase transitions. J Phys A. 2006;39:R143-R205.
- Belitz D, Kirkpatrick TR, Vojta T. How generic scale invariance influences quantum and classical phase transitions. Rev Mod Phys. 2005;77:579–632.
- Landau LD. On the theory of phase transitions.. Ukr J Phys. 1937;(11):19–32.
- Landau LD. Broken symmetry and phase transitions.. Phys Z Sowjetunion. 1937;26:542.
- Wilson KG, Kogut JB. The renormalization group and the epsilon expansion. Phys Rep. 1973;12:75–200.
- Widom B. Surface tension and molecular correlations near the critical point. J Chem Phys. 1965;43:3892–3897.
- Gitterman M. Igor herbut: a modern approach to critical phenomena. J Stat Mech. 2008;131:173.
- Vojta T. Phases and phase transitions in disordered quantum systems. In AIP Conference Proceedings,1550, 188–247. American Institute of Physics,Salerno; 2013.
- Vojta M. Quantum phase transitions. Rep Prog Phys. 2003;66:2069–2110.
- Fisher MPA. Quantum phase transitions in disordered two-dimensional superconductors. Phys Rev Lett. 1990;65:923.
- Harris AB. Effect of random defects on the critical behaviour of ising models. J Phys C. 1974;7:1671–1692.
- Chayes J, Chayes L, Fisher DS, et al. Finite-size scaling and correlation lengths for disordered systems. Phys Rev Lett. 1986;57:2999–3002.
- Vojta T, Hoyos JA. Criticality and quenched disorder: harris criterion versus rare regions. Phys Rev Lett. 2014;112:075702.
- Vojta T, Farquhar A, Mast J. Infinite-randomness critical point in the two-dimensional disordered contact process. Phys Rev E. 2009;79:011111.
- Motrunich OI, Mau S, Huse DA, et al. Infinite-randomness quantum ising critical fixed points. Phys Rev B. 2000;61:1160–1172.
- Kovacs I, Igloi F. Renormalization group study of the two-dimensional random transverse-field ising model. Phys Rev B. 2010;82:054437.
- Pich C, Young AP, Rieger H, et al. Critical behavior and griffiths-mccoy singularities in the two-dimensional random quantum ising ferromagnet. Phys Rev Lett. 1998;81:5916–5919.
- Fisher DS. Critical behavior of random transverse-field ising spin chains. Phys Rev B. 1995;51:6411–6461.
- Griffiths RB. Nonanalytic behavior above the critical point in a random ising ferromagnet. Phys Rev Lett. 1969;23:17–19.
- Mccoy BM, Wu TT. Theory of a two-dimensional ising model with random impurities. i. thermodynamics.. Phys Rev. 1968;176:631–643.
- Fisher DS. Random transverse field ising spin chains. Phys Rev Lett. 1992;69:534–537.
- Fisher DS. Phase transitions and singularities in random quantum systems. Physica A-statistical Mechanics and Its Applications. 1999;263:222–233.
- Del Maestro A, Rosenow B, Muller M, et al. Infinite randomness fixed point of the superconductor-metal quantum phase transition. Phys Rev Lett. 2008;101:035701.
- Lewellyn NA, Percher IM, Nelson JJ, et al. Infinite-randomness fixed point of the quantum superconductor-metal transitions in amorphous thin films. Phys Rev B. 2019 Feb;99:054515.
- Huang C, Zhang E, Zhang Y, et al. The discovery of tunable universality class in superconducting beta-w thin films. arXiv:2010.12775.
- Liu Y, Shichao Q, Sun CL, et al. The observation of in-plane quantum griffiths singularity in two-dimensional crystalline superconductors. arXiv:2011.00268.
- Haviland DB, Liu Y, Goldman AM. Onset of superconductivity in the two-dimensional limit. Phys Rev Lett. 1989;62:2180–2183.
- Hebard AF, Paalanen MA. Magnetic-field-tuned superconductor-insulator transition in two-dimensional films. Phys Rev Lett. 1990;65:927–930.
- Yazdani A, Kapitulnik A. Superconducting-insulating transition in two-dimensional a-moge thin films. Phys Rev Lett. 1995;74:3037–3040.
- Markovic N, Christiansen C, Goldman AM. Thickness–magnetic field phase diagram at the superconductor-insulator transition in 2d. Phys Rev Lett. 1998;81:5217–5220.
- Shi X, Lin PV, Sasagawa T, et al. Two-stage magnetic-field-tuned superconductor–insulator transition in underdoped la 2- x sr x cuo 4. Nat Phys. 2014;10:437–443.
- Ephron D, Yazdani A, Kapitulnik A, et al. Observation of quantum dissipation in the vortex state of a highly disordered superconducting thin film. Phys Rev Lett. 1996;76:1529–1532.
- Maki K. Critical fluctuation of the order parameter in a superconductor. i. Prog Theoretical Phys. 1968;40:193–200.
- Fisher DS, Fisher MPA, Huse DA. Thermal fluctuations, quenched disorder, phase transitions, and transport in type-ii superconductors. Phys Rev B. 1991;43:130.
- Galitski V, Larkin AI. Disorder and quantum fluctuations in superconducting films in strong magnetic fields. Phys Rev Lett. 2001;87:087001.
- Shen S, Xing Y, Wang P, et al. Observation of quantum griffiths singularity and ferromagnetism at the superconducting laal o 3/srti o 3 (110) interface.. Phys Rev B. 2016;94:144517.
- Saito Y, Nojima T, Iwasa Y. Quantum phase transitions in highly crystalline two-dimensional superconductors. Nat Commun. 2018;9:778.
- Hsu Y, Vaezi A, Fischer MH, et al. Topological superconductivity in monolayer transition metal dichalcogenides. Nat Commun. 2017;8:14985.
- Wenyu H, Zhou BT, He JJ, et al. Magnetic field driven nodal topological superconductivity in monolayer transition metal dichalcogenides. Commun Phys. 2018;1:40.
- Zhou BT, Yuan NFQ, Jiang H, et al. Ising superconductivity and majorana fermions in transition metal dichalcogenides. Phys Rev B. 2016;93:180501.
- Liu Y, Wang Z, Shan P, et al. Anomalous quantum griffiths singularity in ultrathin crystalline lead films. Nat Commun.2019;10:1-6.
- Aslamazov LG, Larkin AI. Effect of fluctuations on properties of a superconductor above critical temperature.. Fiz. Tverd. Tela (Leningr.). 1968;10.
- Galitski V, Larkin AI. Superconducting fluctuations at low temperature. Phys Rev B. 2001;63:174506.
- Slevin K, Ohtsuki T. Corrections to scaling at the anderson transition. Phys Rev Lett. 1999 Jan;82:382–385.
- Ando Y, Liang F. Topological crystalline insulators and topological superconductors: from concepts to materials. Ann Rev Condens Matter Phys. 2015;6:361–381.
- Shytov AV, Katsnelson MI, Levitov LS. Atomic collapse and quasi-rydberg states in graphene. Phys Rev Lett. 2007;99:246802.
- Nishida Y. Vacuum polarization of graphene with a supercritical coulomb impurity: low-energy universality and discrete scale invariance. Phys Rev B. 2014;90:165414.
- Ovdat O, Mao J, Jiang Y, et al. Observing a scale anomaly and a universal quantum phase transition in graphene. Nat Commun. 2017;8:1–6.
- Zhang P, Zhai H. Efimov effect in dirac semi-metals. Front Phys. 2018;13:137204.
- Zhang C, Tong B, Yuan Z, et al. Signature of chiral fermion instability in the weyl semimetal taas above the quantum limit. Phys Rev B. 2016;94:205120.
- Shoenberg D. Magnetic oscillations in metals. Cambridge university press, Cambridge; 2009.
- Zhao Y, Liu H, Zhang C, et al. Anisotropic fermi surface and quantum limit transport in high mobility three-dimensional dirac semimetal cd 3 as 2. Phys Rev X. 2015;5:031037.
- Webb RA, Washburn S, Umbach CP, et al. Observation of h e aharonov-bohm oscillations in normal-metal rings. Phys Rev Lett. 1985;54:2696.
- Wang Y, Wong D, Shytov AV, et al. Observing atomic collapse resonances in artificial nuclei on graphene. Science. 2013;340:734–737.
- Braaten E, Hammer HW. Universality in few-body systems with large scattering length. Phys Rep. 2006;428:259–390.
- Efimov V. Energy levels arising from resonant two-body forces in a three-body system. Phys Lett B. 1970;33:563–564.
- Efimov V. Energy levels of three resonantly interacting particles.. Nucl Phys. 1973;210:157–188.
- Wan X, Turner AM, Vishwanath A, et al. Topological semimetal and fermi-arc surface states in the electronic structure of pyrochlore iridates. Phys Rev B. 2011;83:205101.
- Zeldovich YB, Popov VS. Electronic structure of superheavy atoms. Physics-Uspekhi. 1972;14:673–694.
- Greiner W Quantum electrodynamics of strong fields. In: Hadrons and heavy ions. Springer, Berlin; 1985. p. 95–226.
- Shytov AV, Katsnelson MI, Levitov LS. Vacuum polarization and screening of supercritical impurities in graphene. Phys Rev Lett. 2007;99:236801.
- Pereira VM, Nilsson J, Castro Neto AH. Coulomb impurity problem in graphene. Phys Rev Lett. 2007;99:166802.
- Chen RY, Chen ZG, Song X, et al. Magnetoinfrared spectroscopy of landau levels and zeeman splitting of three-dimensional massless dirac fermions in zrte(5). Phys Rev Lett. 2015;115:176404.
- Ruder H, Wunner G, Herold H, et al. Atoms in strong magnetic fields. Springer-Verlag, Berlin; 1994.
- Zhang N, Zhao G, Lin L, et al. Magnetotransport signatures of weyl physics and discrete scale invariance in the elemental semiconductor tellurium. Proc Nat Acad Sci. 2020 May;117:202002913.