Rare Region Effects Research Articles

Modelling on the very large-scale connectome

In this review, we discuss critical dynamics of simple nonequilibrium models on large connectomes, obtained by diffusion MRI, representing the white matter of the human brain. In the first chapter, we overview graph theoretical and topological analysis of these networks, pointing out that universality allows selecting a representative network, the KKI-18, which has been used for dynamical simulation. The critical and sub-critical behaviour of simple, two- or three-state threshold models is discussed with special emphasis on rare-region effects leading to robust Griffiths phases (GP). Numerical results of synchronization phenomena, studied by the Kuramoto model, are also shown, leading to a continuous analog of the GP, termed frustrated synchronization. The models presented here exhibit dynamical scaling behaviour with exponents in agreement with brain experimental data if local homeostasis is provided.

Phenomenological model of the third-harmonic magnetic response due to superconducting fluctuations: Application to Sr2RuO4

We employ the phenomenological Lawrence-Doniach model to compute the contributions of the superconducting fluctuations to the third-harmonic magnetic response, denoted here by $\overline{M_{3}}$, which can be measured in a precise way using ac magnetic fields and lock-in techniques. We show that, in an intermediate temperature regime, this quantity behaves as the third-order nonlinear susceptibility, which shows a power-law dependence with the reduced temperature $\epsilon=\frac{T-T_{c}}{T_{c}}$ as $\epsilon^{-5/2}$. Very close to $T_{c}$, however, $\overline{M_{3}}$ saturates due to the nonzero amplitude of the ac field. We compare our theoretical results with experimental data for three conventional superconductors -- lead, niobium, and vanadium -- and for the unconventional superconductor Sr$_{2}$RuO$_{4}$ (SRO). We find good agreement between theory and experiment for the elemental superconductors, although the theoretical values for the critical field systematically deviate from the experimental ones. In the case of SRO, however, the phenomenological model completely fails to describe the data, as the third-harmonic response remains sizable over a much wider reduced temperature range compared to Pb, Nb, and V. We show that an inhomogeneous distribution of $T_{c}$ can partially account for this discrepancy, since regions with a locally higher $T_{c}$ contribute to the fluctuation $\overline{M_{3}}$ significantly more than regions with the "nominal" $T_{c}$ of the clean system. However, the exponential temperature dependence of $\overline{M_{3}}$ first reported in Ref. [3] is not captured by the model with inhomogeneity. We conclude that, while inhomogeneity is an important ingredient to understand the superconducting fluctuations of SRO and other perovskite superconductors, additional effects may be at play, such as non-Gaussian fluctuations or rare-region effects.

Rare regions and avoided quantum criticality in disordered Weyl semimetals and superconductors

Disorder in Weyl semimetals and superconductors is surprisingly subtle, attracting attention and competing theories in recent years. In this brief review, we discuss the current theoretical understanding of the effects of short-ranged, quenched disorder on the low energy-properties of three-dimensional, topological Weyl semimetals and superconductors. We focus on the role of non-perturbative rare region effects on destabilizing the semimetal phase and rounding the expected semimetal-to-diffusive metal transition into a cross over. Furthermore, the consequences of disorder on the resulting nature of excitations, transport, and topology are reviewed. New results on a bipartite random hopping model are presented that confirm previous results in a p+ip Weyl superconductor, demonstrating that particle–hole symmetry is insufficient to help stabilize the Weyl semimetal phase in the presence of disorder. The nature of the avoided transition in a model for a single Weyl cone in the continuum is discussed. We close with a discussion of open questions and future directions.

Epidemic growth and Griffiths effects on an emergent network of excited atoms

Whether it be physical, biological or social processes, complex systems exhibit dynamics that are exceedingly difficult to understand or predict from underlying principles. Here we report a striking correspondence between the excitation dynamics of a laser driven gas of Rydberg atoms and the spreading of diseases, which in turn opens up a controllable platform for studying non-equilibrium dynamics on complex networks. The competition between facilitated excitation and spontaneous decay results in sub-exponential growth of the excitation number, which is empirically observed in real epidemics. Based on this we develop a quantitative microscopic susceptible-infected-susceptible model which links the growth and final excitation density to the dynamics of an emergent heterogeneous network and rare active region effects associated to an extended Griffiths phase. This provides physical insights into the nature of non-equilibrium criticality in driven many-body systems and the mechanisms leading to non-universal power-laws in the dynamics of complex systems.

Scaling of local persistence in the disordered contact process.

We study the time dependence of the local persistence probability during a nonstationary time evolution in the disordered contact process in d=1, 2, and 3 dimensions. We present a method for calculating the persistence with the strong-disorder renormalization group (SDRG) technique, which we then apply at the critical point analytically for d=1 and numerically for d=2,3. According to the results, the average persistence decays at late times as an inverse power of the logarithm of time, with a universal dimension-dependent generalized exponent. For d=1, the distribution of sample-dependent local persistence is shown to be characterized by a universal limit distribution of effective persistence exponents. Using a phenomenological approach of rare-region effects in the active phase, we obtain a nonuniversal algebraic decay of the average persistence for d=1 and enhanced power laws for d>1. As an exception, for randomly diluted lattices, the algebraic decay remains valid for d>1, which is explained by the contribution of dangling ends. Results on the time dependence of average persistence are confirmed by Monte Carlo simulations. We also prove the equivalence of the persistence with a return probability, a valuable tool for the argumentations.

Dipolar spin glass transition in three dimensions

Dilute dipolar Ising magnets remain a notoriously hard problem to tackle both analytically and numerically because of long-ranged interactions between spins as well as rare region effects. We study a new type of anisotropic dilute dipolar Ising system in three dimensions [Phys. Rev. Lett. {\bf 114}, 247207 (2015)] that arises as an effective description of randomly diluted classical spin ice, a prototypical spin liquid in the disorder-free limit, with a small fraction $x$ of non-magnetic impurities. Metropolis algorithm within a parallel thermal tempering scheme fails to achieve equilibration for this problem already for small system sizes. Motivated by previous work [Phys. Rev. X {\bf 4}, 041016 (2014)] on uniaxial random dipoles, we present an improved cluster Monte Carlo algorithm that is tailor-made for removing the equilibration bottlenecks created by clusters of {\it effectively frozen} spins. By performing large-scale simulations down to $x=1/128$ and using finite size scaling, we show the existence of a finite-temperature spin glass transition and give strong evidence that the universality of the critical point is independent of $x$ when it is small. In this $x \ll 1$ limit, we also provide a first estimate of both the thermal exponent, $\nu=1.27(8)$, and the anomalous exponent, $\eta=0.228(35)$.

Self-consistent Hartree-Fock approach to many-body localization

In this work, we develop a self-consistent Hartree-Fock approach to theoretically study the far-from-equilibrium quantum dynamics of interacting fermions, and apply this approach to explore the onset of many-body localization (MBL) in these systems. We investigate the dynamics of a state with a nonequilibrium density profile; we find that at weak disorder the density profile equilibrates rapidly, whereas for strong disorder it remains frozen on the accessible timescales. We analyze this behavior in terms of the Hartree-Fock self-energy. At weak disorder the self-energy fluctuates strongly and can be interpreted as a self-consistent noise process. By contrast, at strong disorder the self-energy evolves with a few coherent oscillations which cannot delocalize the system. Accordingly, the non-equilibrium site-resolved spectral function shows a broad spectrum at weak disorder and sharp spikes at strong disorder. Our Hartree-Fock theory incorporates spatial fluctuations and rare-region effects. As a consequence, we find subdiffusive relaxation in random systems; but, when the system is subjected to weak quasi-periodic potentials, the subdiffusive response ceases to exist, as rare region effects are absent in this case. This self-consistent Hartree-Fock approach can be regarded as a relatively simple theory that captures much of the MBL phenomenology.

Heterogeneity effects in power grid network models.

We have compared the phase synchronization transition of the second-order Kuramoto model on two-dimensional (2D) lattices and on large, synthetic power grid networks, generated from real data. The latter are weighted, hierarchical modular networks. Due to the inertia the synchronization transitions are of first-order type, characterized by fast relaxation and hysteresis by varying the global coupling parameter K. Finite-size scaling analysis shows that there is no real phase transition in the thermodynamic limit, unlike in the mean-field model. The order parameter and its fluctuations depend on the network size without any real singular behavior. In case of power grids the phase synchronization breaks down at lower global couplings, than in case of 2D lattices of the same sizes, but the hysteresis is much narrower or negligible due to the low connectivity of the graphs. The temporal behavior of desynchronization avalanches after a sudden quench to low K values has been followed and duration distributions with power-law tails have been detected. This suggests rare region effects, caused by frozen disorder, resulting in heavy-tailed distributions, even without a self-organization mechanism as a consequence of a catastrophic drop event in the couplings.

Smeared nematic quantum phase transitions due to rare-region effects in inhomogeneous systems

The concept of a vestigial nematic order emerging from a "mother" spin or charge density-wave state has been applied to describe the phase diagrams of several systems, including unconventional superconductors. In a perfectly clean system, the two orders appear simultaneously via a first-order quantum phase transition, implying the absence of quantum criticality. Here, we investigate how this behavior is affected by impurity-free droplets that are naturally present in inhomogeneous systems. Due to their quantum dynamics, finite-size droplets sustain long-range nematic order but not long-range density-wave order. Interestingly, rare droplets with moderately large sizes undergo a second-order nematic transition even before the first-order quantum transition of the clean system. This gives rise to an extended regime of inhomogeneous nematic order, which is followed by a density-wave quantum Griffiths phase. As a result, a smeared quantum nematic transition, separated from the density-wave quantum transition, emerges in moderately disordered systems.

Theory of avoided criticality in quantum motion in a random potential in high dimensions

The density of states of a three dimensional Dirac equation with a random potential as well as in other problems of quantum motion in a random potential placed in sufficiently high spatial dimensionality appears to be singular at a certain critical disorder strength. This was seen numerically in a variety of studies as well as supported by detailed renormalization group calculations. At the same time it was suggested by a number of arguments accompanied by detailed numerical simulations that this singularity is rounded off by the rare region fluctuations of random potential, and that tuning the disorder past its critical value is not a genuine phase transition but rather a crossover. Here we develop an analytic theory which explains how rare region effects indeed lead to rounding off of the singularity and to the crossover replacing the transition.

Rare‐region effects and dynamics near the many‐body localization transition

The low‐frequency response of systems near the many‐body localization phase transition, on either side of the transition, is dominated by contributions from rare regions that are locally “in the other phase”, i.e., rare localized regions in a system that is typically thermal, or rare thermal regions in a system that is typically localized. Rare localized regions affect the properties of the thermal phase, especially in one dimension, by acting as bottlenecks for transport and the growth of entanglement, whereas rare thermal regions in the localized phase act as local “baths” and dominate the low‐frequency response of the MBL phase. We review recent progress in understanding these rare‐region effects, and discuss some of the open questions associated with them: in particular, whether and in what circumstances a single rare thermal region can destabilize the many‐body localized phase.image

Rare-Region-Induced Avoided Quantum Criticality in Disordered Three-Dimensional Dirac and Weyl Semimetals

We numerically study the effect of short ranged potential disorder on massless noninteracting three-dimensional Dirac and Weyl fermions, with a focus on the question of the proposed quantum critical point separating the semimetal and diffusive metal phases. We determine the properties of the eigenstates of the disordered Dirac Hamiltonian ($H$) and exactly calculate the density of states (DOS) near zero energy, using a combination of Lanczos on $H^2$ and the kernel polynomial method on $H$. We establish the existence of two distinct types of low energy eigenstates contributing to the disordered density of states in the weak disorder semimetal regime. These are (i) typical eigenstates that are well described by linearly dispersing perturbatively dressed Dirac states, and (ii) nonperturbative rare eigenstates that are weakly-dispersive and quasi-localized in the real space regions with the largest (and rarest) local random potential. Using twisted boundary conditions, we are able to systematically find and study these two types of eigenstates. We find that the Dirac states contribute low energy peaks in the finite-size DOS that arise from the clean eigenstates which shift and broaden in the presence of disorder. On the other hand, we establish that the rare quasi-localized eigenstates contribute a nonzero background DOS which is only weakly energy-dependent near zero energy and is exponentially small at weak disorder. We find that the expected semimetal to diffusive metal quantum critical point is converted to an {\it avoided} quantum criticality that is "rounded out" by nonperturbative effects, with no signs of any singular behavior in the DOS near the Dirac energy. We discuss the implications of our results for disordered Dirac and Weyl semimetals, and reconcile the large body of existing numerical work showing quantum criticality with the existence of the rare region effects.

Griffiths effects and slow dynamics in nearly many-body localized systems

The low-frequency response of systems near a many-body localization transition can be dominated by rare regions that are locally critical or "in the other phase". It is known that, in one dimension, these rare regions can cause the d.c. conductivity and diffusion constant to vanish even inside the delocalized thermal phase. Here, we present a general analysis of such Griffiths effects in the thermal phase near the many-body localization transition: we consider both one-dimensional and higher-dimensional systems, subject to quenched randomness, and discuss both linear response (including the frequency- and wavevector-dependent conductivity) and more general dynamics. In all the regimes we consider, we identify observables that are dominated by rare-region effects. In some cases (one-dimensional systems and Floquet systems with no extensive conserved quantities), essentially all long-time local observables are dominated by rare-region effects; in others, generic observables are instead dominated by hydrodynamic long-time tails throughout the thermal phase, and one must look at specific probes, such as spin echo, to see Griffiths behavior.

Enhanced rare-region effects in the contact process with long-range correlated disorder.

We investigate the nonequilibrium phase transition in the disordered contact process in the presence of long-range spatial disorder correlations. These correlations greatly increase the probability for finding rare regions that are locally in the active phase while the bulk system is still in the inactive phase. Specifically, if the correlations decay as a power of the distance, the rare-region probability is a stretched exponential of the rare-region size rather than a simple exponential as is the case for uncorrelated disorder. As a result, the Griffiths singularities are enhanced and take a non-power-law form. The critical point itself is of infinite-randomness type but with critical exponent values that differ from the uncorrelated case. We report large-scale Monte Carlo simulations that verify and illustrate our theory. We also discuss generalizations to higher dimensions and applications to other systems such as the random transverse-field Ising model, itinerant magnets, and the superconductor-metal transition.

Localization transition, Lifschitz tails, and rare-region effects in network models.

Effects of heterogeneity in the suspected-infected-susceptible model on networks are investigated using quenched mean-field theory. The emergence of localization is described by the distributions of the inverse participation ratio and compared with the rare-region effects appearing in simulations and in the Lifschitz tails. The latter, in the linear approximation, is related to the spectral density of the Laplacian matrix and to the time dependent order parameter. I show that these approximations indicate correctly Griffiths phases both on regular one-dimensional lattices and on small-world networks exhibiting purely topological disorder. I discuss the localization transition that occurs on scale-free networks at γ=3 degree exponent.

Rare region effects dominate weakly disordered three-dimensional Dirac points

We study three-dimensional Dirac fermions with weak finite-range scalar potential disorder. We show that even though disorder is perturbatively irrelevant at 3D Dirac points, nonperturbative effects from rare regions give rise to a nonzero density of states and a finite mean free path, with the transport at the Dirac point being dominated by hopping between rare regions. As one moves in chemical potential away from the Dirac point, there are interesting intermediate-energy regimes where the rare regions produce scattering resonances that determine the DC conductivity. We also discuss the interplay of disorder with interactions at the Dirac point. Attractive interactions drive a transition into a granular superconductor, with a critical temperature that depends strongly on the disorder distribution. In the presence of Coulomb repulsion and weak retarded attraction, the system can be a Bose glass. Our results apply to all 3D systems with Dirac points, including Weyl semimetals, and overturn a thirty year old consensus regarding the irrelevance of weak disorder at 3D Dirac points.

Slow, bursty dynamics as a consequence of quenched network topologies.

Bursty dynamics of agents is shown to appear at criticality or in extended Griffiths phases, even in case of Poisson processes. I provide numerical evidence for a power-law type of intercommunication time distributions by simulating the contact process and the susceptible-infected-susceptible model. This observation suggests that in the case of nonstationary bursty systems, the observed non-Poissonian behavior can emerge as a consequence of an underlying hidden Poissonian network process, which is either critical or exhibits strong rare-region effects. On the contrary, in time-varying networks, rare-region effects do not cause deviation from the mean-field behavior, and heterogeneity-induced burstyness is absent.

Spectral analysis and slow spreading dynamics on complex networks

The susceptible-infected-susceptible (SIS) model is one of the simplest memoryless systems for describing information or epidemic spreading phenomena with competing creation and spontaneous annihilation reactions. The effect of quenched disorder on the dynamical behavior has recently been compared to quenched mean-field (QMF) approximations in scale-free networks. QMF can take into account topological heterogeneity and clustering effects of the activity in the steady state by spectral decomposition analysis of the adjacency matrix. Therefore, it can provide predictions on possible rare-region effects, thus on the occurrence of slow dynamics. I compare QMF results of SIS with simulations on various large dimensional graphs. In particular, I show that for Erdős-Rényi graphs this method predicts correctly the occurrence of rare-region effects. It also provides a good estimate for the epidemic threshold in case of percolating graphs. Griffiths Phases emerge if the graph is fragmented or if we apply a strong, exponentially suppressing weighting scheme on the edges. The latter model describes the connection time distributions in the face-to-face experiments. In case of a generalized Barabási-Albert type of network with aging connections, strong rare-region effects and numerical evidence for Griffiths Phase dynamics are shown. The dynamical simulation results agree well with the predictions of the spectral analysis applied for the weighted adjacency matrices.

Slow epidemic extinction in populations with heterogeneous infection rates

We explore how heterogeneity in the intensity of interactions between people affects epidemic spreading. For that, we study the susceptible-infected-susceptible model on a complex network, where a link connecting individuals i and j is endowed with an infection rate β(ij)=λw(ij) proportional to the intensity of their contact w(ij), with a distribution P(w(ij)) taken from face-to-face experiments analyzed in Cattuto et al. [PLoS ONE 5, e11596 (2010)]. We find an extremely slow decay of the fraction of infected individuals, for a wide range of the control parameter λ. Using a distribution of width a we identify two large regions in the a-λ space with anomalous behaviors, which are reminiscent of rare region effects (Griffiths phases) found in models with quenched disorder. We show that the slow approach to extinction is caused by isolated small groups of highly interacting individuals, which keep epidemics alive for very long times. A mean-field approximation and a percolation approach capture with very good accuracy the absorbing-active transition line for weak (small a) and strong (large a) disorder, respectively.

Rare regions of the susceptible-infected-susceptible model on Barabási-Albert networks

I extend a previous work to susceptible-infected-susceptible (SIS) models on weighted Barabási-Albert scale-free networks. Numerical evidence is provided that phases with slow, power-law dynamics emerge as the consequence of quenched disorder and tree topologies studied previously with the contact process. I compare simulation results with spectral analysis of the networks and show that the quenched mean-field (QMF) approximation provides a reliable, relatively fast method to explore activity clustering. This suggests that QMF can be used for describing rare-region effects due to network inhomogeneities. Finite-size study of the QMF shows the expected disappearance of the epidemic threshold λ(c) in the thermodynamic limit and an inverse participation ratio ~0.25, meaning localization in case of disassortative weight scheme. Contrarily, for the multiplicative weights and the unweighted trees, this value vanishes in the thermodynamic limit, suggesting only weak rare-region effects in agreement with the dynamical simulations. Strong corrections to the mean-field behavior in case of disassortative weights explains the concave shape of the order parameter ρ(λ) at the transition point. Application of this method to other models may reveal interesting rare-region effects, Griffiths phases as the consequence of quenched topological heterogeneities.