Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2002.11722 · Physical review. B./Physical review. B

Avoided quantum criticality in exact numerical simulations of a single disordered Weyl cone

Justin H. Wilson, David A. Huse, S. Das Sarma, J. H. Pixley

Published Sep 8, 2020• 18 citations
arXivHTMLPDF

Abstract

Existing theoretical works differ on whether three-dimensional Dirac and Weyl semimetals are stable to a short-range-correlated random potential. Numerical evidence suggests the semimetal to be unstable, while some field-theoretic instanton calculations have found it to be stable. The differences go beyond method: the continuum field-theoretic works use a single, perfectly linear Weyl cone, while numerical works use tight-binding lattice models which inherently have band curvature and multiple Weyl cones. In this work, we bridge this gap by performing exact numerics on the same model used in analytic treatments, and we find that all phenomena associated with rare regions near the Weyl node energy found in lattice models persist in the continuum theory: The density of states is nonzero and exhibits an avoided transition. In addition to characterizing this transition, we find rare states and show that they have the expected behavior. The simulations utilize sparse matrix techniques with formally dense matrices; doing so allows us to reach Hilbert space sizes upwards of ${10}^{7}$ states, substantially larger than anything achieved before.

DOI
10.1103/physrevb.102.100201

Related papers

Uncovering the hidden quantum critical point in disordered massless Dirac and Weyl semimetals

J. H. Pixley, David A. Huse, S. Das Sarma

We study the properties of the avoided or hidden quantum critical point (AQCP) in three-dimensional Dirac and Weyl semimetals in the presence of short range potential disorder. By computing the averaged density of states (along with its second and fourth derivative at zero energy) with the kernel polynomial method (KPM) we systematically tune the effective length scale that eventually rounds out the transition and leads to an AQCP. We show how to determine the strength of the avoidance, establishing that it is not controlled by the long wavelength component of the disorder. Instead, the amount of avoidance can be adjusted via the tails of the probability distribution of the local random potentials. A binary distribution with no tails produces much less avoidance than a Gaussian distribution. We introduce a double Gaussian distribution to interpolate between these two limits. As a result we are able to make the length scale of the avoidance sufficiently large so that we can accurately study the properties of the underlying transition (that is eventually rounded out), unambiguously identify its location, and provide accurate estimates of the critical exponents $\ensuremath{\nu}=1.01\ifmmode\pm\else\textpm\fi{}0.06$ and $z=1.50\ifmmode\pm\else\textpm\fi{}0.04$. We also show that the KPM expansion order introduces an effective length scale that can also round out the transition in the scaling regime near the AQCP.

arXiv:1606.08860•Physical review. B./Physical review. B•Sep 8, 2016•50 citations

Stability of the topological quantum critical point between multi-Weyl semimetal and band insulator

Zhao-Kun Yang, Jing-Rong Wang, Guo-Zhu Liu

One could tune a topological double-Weyl semimetal or a topological triple-Weyl semimetal to become a topologically trivial insulator by opening a band gap. This kind of quantum phase transition is characterized by the change of certain topological invariant. A new gapless semimetallic state emerges at each topological quantum critical point. Here we perform a renormalization group analysis to investigate the stability of such critical points against perturbations induced by random scalar potential and random vector potential. We find that the quantum critical point between double-Weyl semimetal and band insulator is unstable and can be easily turned into a compressible diffusive metal by any type of weak disorder. The quantum critical point between triple-Weyl semimetal and band insulator flows to a stable strong-coupling fixed point if the system contains a random vector potential merely along the z-axis, but becomes a compressible diffusive metal when other types of disorders exist.

Chinese Physics B•Feb 10, 2023•0 citations

Upper critical dimension for dirty Weyl semimetal-to-metal quantum phase transitions

Yongtai Li, Bitan Roy

Weyl or Dirac fermions with the iconic linear energy-momentum relation and average density of states (ADOS) $ρ(E) \sim |E|^{d-1}$ at energy $E$ in $d$ spatial dimensions, constitute a unique setup to study the disorder-driven semimetal-to-metal quantum phase transition (QPT) between ballistic (realized for weak disorder) and diffusive (stabilized at stronger disorder) quasiparticles. Such a QPT takes place only for $d>2$ and falls beyond the realm of the Anderson metal-to-insulator transition. From numerically computed ADOS (using the kernel polynomial method) in dirty Weyl systems in $d=2$ to $6$, here we show that $d=2$ and $d=4$ are the lower ($d_\ell$) and upper ($d_u$) critical dimensions for such an unconventional QPT, respectively, as suggested from the solution of quasiparticle lifetime within the self-consistent Born approximation. Consequently, for $d \geq 4$ the associated correlation length exponent is found to be $ν\approx 0.5$ (within numerical accuracy). However, the dynamic scaling exponent at the quantum critical point is pinned close to $z \approx d/2$ (numerically) for any $d \geq 3$, which is shown to be an exact result from a field-theoretic renormalization group calculation. Therefore, Weyl semimetal-to-metal QPTs can be studied field theoretically around both $d_\ell$ and $d_u$.

arXiv:2609.30265•Sep 24, 2026•0 citations

On the disorder-driven quantum transition in three-dimensional relativistic metals

Thibaud Louvet, David Carpentier, Andrei A. Fedorenko

The Weyl semimetals are topologically protected from a gap opening against weak disorder in three dimensions. However, a strong disorder drives this relativistic semimetal through a quantum transition towards a diffusive metallic phase characterized by a finite density of states at the band crossing. This transition is usually described by a perturbative renormalization group in $d=2+\ensuremath{\varepsilon}$ of a $U(N)$ Gross-Neveu model in the limit $N\ensuremath{\rightarrow}0$. Unfortunately, this model is not multiplicatively renormalizable in $2+\ensuremath{\varepsilon}$ dimensions: An infinite number of relevant operators are required to describe the critical behavior. Hence its use in a quantitative description of the transition beyond one loop is at least questionable. We propose an alternative route, building on the correspondence between the Gross-Neveu and Gross-Neveu-Yukawa models developed in the context of high-energy physics. It results in a model of Weyl fermions with a random non-Gaussian imaginary potential which allows one to study the critical properties of the transition within a $d=4\ensuremath{-}\ensuremath{\varepsilon}$ expansion. We also discuss the characterization of the transition by the multifractal spectrum of wave functions.

arXiv:1605.02009•Physical review. B./Physical review. B•Dec 5, 2016•45 citations

Disorder-driven itinerant quantum criticality of three-dimensional massless Dirac fermions

J. H. Pixley, Pallab Goswami, S. Das Sarma

Progress in the understanding of quantum critical properties of itinerant electrons has been hindered by the lack of effective models which are amenable to controlled analytical and numerically exact calculations. Here we establish that the disorder-driven semimetal to metal quantum phase transition of three-dimensional massless Dirac fermions could serve as a paradigmatic toy model for studying itinerant quantum criticality, which is solved in this work by exact numerical and approximate field-theoretic calculations. As a result, we establish the robust existence of a non-Gaussian universality class, and also construct the relevant low-energy effective field theory that could guide the understanding of quantum critical scaling for many strange metals. Using the kernel polynomial method (KPM), we provide numerical results for the calculated dynamical exponent ($z$) and correlation length exponent ($\ensuremath{\nu}$) for the disorder-driven semimetal (SM) to diffusive metal (DM) quantum phase transition at the Dirac point for several types of disorder, establishing its universal nature and obtaining the numerical scaling functions in agreement with our field-theoretical analysis.

arXiv:1505.07938•Physical review. B./Physical review. B•Feb 1, 2016•61 citations

Direct topological insulator transitions in three dimensions are destabilized by nonperturbative effects of disorder

Yixing Fu, Justin H. Wilson, David A. Huse +1 more

We reconsider the phase diagram of a three-dimensional ${\mathbb{Z}}_{2}$ topological insulator in the presence of short-ranged potential disorder, with the insight that nonperturbative rare states destabilize the noninteracting Dirac semimetal critical point separating different topological phases. Based on our numerical data on the density of states, conductivity, and wave functions, we argue that the putative Dirac semimetal line is destabilized into a diffusive metal phase of finite extent due to nonperturbative effects of rare regions. We discuss the implications of these results for past and current experiments on doped topological insulators.

arXiv:2309.09857•Physical review. B./Physical review. B•May 1, 2024•3 citations
Existing theoretical works differ on whether three-dimensional Dirac and Weyl semimetals are stable to a short-range-correlated random potential. Numerical evidence suggests the semimetal to be unstable, while some field-theoretic instanton calculations have found it to be stable. The differences go beyond method: the continuum field-theoretic works use a single, perfectly linear Weyl cone, while numerical works use tight-binding lattice models which inherently have band curvature and multiple Weyl cones. In this work, we bridge this gap by performing exact numerics on the same model used in analytic treatments, and we find that all phenomena associated with rare regions near the Weyl node energy found in lattice models persist in the continuum theory: The density of states is nonzero and exhibits an avoided transition. In addition to characterizing this transition, we find rare states and show that they have the expected behavior. The simulations utilize sparse matrix techniques with formally dense matrices; doing so allows us to reach Hilbert space sizes upwards of 107{10}^{7}107 states, substantially larger than anything achieved before.
We study the properties of the avoided or hidden quantum critical point (AQCP) in three-dimensional Dirac and Weyl semimetals in the presence of short range potential disorder. By computing the averaged density of states (along with its second and fourth derivative at zero energy) with the kernel polynomial method (KPM) we systematically tune the effective length scale that eventually rounds out the transition and leads to an AQCP. We show how to determine the strength of the avoidance, establishing that it is not controlled by the long wavelength component of the disorder. Instead, the amount of avoidance can be adjusted via the tails of the probability distribution of the local random potentials. A binary distribution with no tails produces much less avoidance than a Gaussian distribution. We introduce a double Gaussian distribution to interpolate between these two limits. As a result we are able to make the length scale of the avoidance sufficiently large so that we can accurately study the properties of the underlying transition (that is eventually rounded out), unambiguously identify its location, and provide accurate estimates of the critical exponents $\ensuremath{\nu}=1.01\ifmmode\pm\else\textpm\fi{}0.06$ and $z=1.50\ifmmode\pm\else\textpm\fi{}0.04$. We also show that the KPM expansion order introduces an effective length scale that can also round out the transition in the scaling regime near the AQCP.
One could tune a topological double-Weyl semimetal or a topological triple-Weyl semimetal to become a topologically trivial insulator by opening a band gap. This kind of quantum phase transition is characterized by the change of certain topological invariant. A new gapless semimetallic state emerges at each topological quantum critical point. Here we perform a renormalization group analysis to investigate the stability of such critical points against perturbations induced by random scalar potential and random vector potential. We find that the quantum critical point between double-Weyl semimetal and band insulator is unstable and can be easily turned into a compressible diffusive metal by any type of weak disorder. The quantum critical point between triple-Weyl semimetal and band insulator flows to a stable strong-coupling fixed point if the system contains a random vector potential merely along the z-axis, but becomes a compressible diffusive metal when other types of disorders exist.
Weyl or Dirac fermions with the iconic linear energy-momentum relation and average density of states (ADOS) ρ(E)∼∣E∣d−1ρ(E) \sim |E|^{d-1}ρ(E)∼∣E∣d−1 at energy EEE in ddd spatial dimensions, constitute a unique setup to study the disorder-driven semimetal-to-metal quantum phase transition (QPT) between ballistic (realized for weak disorder) and diffusive (stabilized at stronger disorder) quasiparticles. Such a QPT takes place only for d>2d>2d>2 and falls beyond the realm of the Anderson metal-to-insulator transition. From numerically computed ADOS (using the kernel polynomial method) in dirty Weyl systems in d=2d=2d=2 to 666, here we show that d=2d=2d=2 and d=4d=4d=4 are the lower (dℓd_\elldℓ​) and upper (dud_udu​) critical dimensions for such an unconventional QPT, respectively, as suggested from the solution of quasiparticle lifetime within the self-consistent Born approximation. Consequently, for d≥4d \geq 4d≥4 the associated correlation length exponent is found to be ν≈0.5ν\approx 0.5ν≈0.5 (within numerical accuracy). However, the dynamic scaling exponent at the quantum critical point is pinned close to z≈d/2z \approx d/2z≈d/2 (numerically) for any d≥3d \geq 3d≥3, which is shown to be an exact result from a field-theoretic renormalization group calculation. Therefore, Weyl semimetal-to-metal QPTs can be studied field theoretically around both dℓd_\elldℓ​ and dud_udu​.
The Weyl semimetals are topologically protected from a gap opening against weak disorder in three dimensions. However, a strong disorder drives this relativistic semimetal through a quantum transition towards a diffusive metallic phase characterized by a finite density of states at the band crossing. This transition is usually described by a perturbative renormalization group in $d=2+\ensuremath{\varepsilon}$ of a U(N)U(N)U(N) Gross-Neveu model in the limit $N\ensuremath{\rightarrow}0$. Unfortunately, this model is not multiplicatively renormalizable in $2+\ensuremath{\varepsilon}$ dimensions: An infinite number of relevant operators are required to describe the critical behavior. Hence its use in a quantitative description of the transition beyond one loop is at least questionable. We propose an alternative route, building on the correspondence between the Gross-Neveu and Gross-Neveu-Yukawa models developed in the context of high-energy physics. It results in a model of Weyl fermions with a random non-Gaussian imaginary potential which allows one to study the critical properties of the transition within a $d=4\ensuremath{-}\ensuremath{\varepsilon}$ expansion. We also discuss the characterization of the transition by the multifractal spectrum of wave functions.
Progress in the understanding of quantum critical properties of itinerant electrons has been hindered by the lack of effective models which are amenable to controlled analytical and numerically exact calculations. Here we establish that the disorder-driven semimetal to metal quantum phase transition of three-dimensional massless Dirac fermions could serve as a paradigmatic toy model for studying itinerant quantum criticality, which is solved in this work by exact numerical and approximate field-theoretic calculations. As a result, we establish the robust existence of a non-Gaussian universality class, and also construct the relevant low-energy effective field theory that could guide the understanding of quantum critical scaling for many strange metals. Using the kernel polynomial method (KPM), we provide numerical results for the calculated dynamical exponent (zzz) and correlation length exponent ($\ensuremath{\nu}$) for the disorder-driven semimetal (SM) to diffusive metal (DM) quantum phase transition at the Dirac point for several types of disorder, establishing its universal nature and obtaining the numerical scaling functions in agreement with our field-theoretical analysis.
We reconsider the phase diagram of a three-dimensional Z2{\mathbb{Z}}_{2}Z2​ topological insulator in the presence of short-ranged potential disorder, with the insight that nonperturbative rare states destabilize the noninteracting Dirac semimetal critical point separating different topological phases. Based on our numerical data on the density of states, conductivity, and wave functions, we argue that the putative Dirac semimetal line is destabilized into a diffusive metal phase of finite extent due to nonperturbative effects of rare regions. We discuss the implications of these results for past and current experiments on doped topological insulators.