Renormalization group analysis of the Anderson model on random regular graphs.

We present a renormalization group (RG) analysis of the problem of Anderson localization on a random regular graph (RRG) which generalizes the RG of Abrahams, Anderson, Licciardello, and Ramakrishnan to infinite-dimensional graphs. The RG equations necessarily involve two parameters (one being the changing connectivity of subtrees), but we show that the one-parameter scaling hypothesis is recovered for sufficiently large system sizes for both eigenstates and spectrum observables. We also explain the nonmonotonic behavior of dynamical and spectral quantities as a function of the system size for values of disorder close to the transition, by identifying two terms in the beta function of the running fractal dimension of different signs and functional dependence. Our theory provides a simple and coherent explanation for the unusual scaling behavior observed in numerical data of the Anderson model on RRG and of many-body localization.

Localization and Melting of Interfaces in the Two-Dimensional Quantum Ising Model.

We study the nonequilibrium evolution of coexisting ferromagnetic domains in the two-dimensional quantum Ising model-a setup relevant in several contexts, from quantum nucleation dynamics and false-vacuum decay scenarios to recent experiments with Rydberg-atom arrays. We demonstrate that the quantum-fluctuating interface delimiting a large bubble can be studied as an effective one-dimensional system through a "holographic" mapping. For the considered model, the emergent interface excitations map to an integrable chain of fermionic particles. We discuss how this integrability is broken by geometric features of the bubbles and by corrections in inverse powers of the ferromagnetic coupling, and provide a lower bound to the timescale after which the bubble is ultimately expected to melt. Remarkably, we demonstrate that a symmetry-breaking longitudinal field gives rise to a robust ergodicity breaking in two dimensions, a phenomenon underpinned by Stark many-body localization of the emergent fermionic excitations of the interface.

Constraint-Induced Delocalization.

We study the impact of quenched disorder on the dynamics of locally constrained quantum spin chains, that describe 1D arrays of Rydberg atoms in both the frozen (Ising-type) and dressed (XY-type) regime. Performing large-scale numerical experiments, we observe no trace of many-body localization even at large disorder. Analyzing the role of quenched disorder terms in constrained systems we show that they act in two, distinct and competing ways: as an on-site disorder term for the basic excitations of the system, and as an interaction between excitations. The two contributions are of the same order, and as they compete (one towards localization, the other against it), one does never enter a truly strong disorder, weak interaction limit, where many-body localization occurs. Such a mechanism is further clarified in the case of XY-type constrained models: there, a term which would represent a bona fide local quenched disorder term acting on the excitations of the clean model must be written as a series of nonlocal terms in the unconstrained variables. Our observations provide a simple picture to interpret the role of quenched disorder that could be immediately extended to other constrained models or quenched gauge theories.

Impact of jamming criticality on low-temperature anomalies in structural glasses.

We present a mechanism for the anomalous behavior of the specific heat in low-temperature amorphous solids. The analytic solution of a mean-field model belonging to the same universality class as high-dimensional glasses, the spherical perceptron, suggests that there exists a cross-over temperature above which the specific heat scales linearly with temperature, while below it, a cubic scaling is displayed. This relies on two crucial features of the phase diagram: (i) the marginal stability of the free-energy landscape, which induces a gapless phase responsible for the emergence of a power-law scaling; and (ii) the vicinity of the classical jamming critical point, as the cross-over temperature gets lowered when approaching it. This scenario arises from a direct study of the thermodynamics of the system in the quantum regime, where we show that, contrary to crystals, the Debye approximation does not hold.

Many-Body Localization Dynamics from Gauge Invariance.

We show how lattice gauge theories can display many-body localization dynamics in the absence of disorder. Our starting point is the observation that, for some generic translationally invariant states, the Gauss law effectively induces a dynamics which can be described as a disorder average over gauge superselection sectors. We carry out extensive exact simulations on the real-time dynamics of a lattice Schwinger model, describing the coupling between U(1) gauge fields and staggered fermions. Our results show how memory effects and slow, double-logarithmic entanglement growth are present in a broad regime of parameters-in particular, for sufficiently large interactions. These findings are immediately relevant to cold atoms and trapped ion experiments realizing dynamical gauge fields and suggest a new and universal link between confinement and entanglement dynamics in the many-body localized phase of lattice models.

Diffusive and Subdiffusive Spin Transport in the Ergodic Phase of a Many-Body Localizable System.

We study high temperature spin transport in a disordered Heisenberg chain in the ergodic regime. By employing a density matrix renormalization group technique for the study of the stationary states of the boundary-driven Lindblad equation we are able to study extremely large systems (400 spins). We find both a diffusive and a subdiffusive phase depending on the strength of the disorder and on the anisotropy parameter of the Heisenberg chain. Studying finite-size effects, we show numerically and theoretically that a very large crossover length exists that controls the passage of a clean-system dominated dynamics to one observed in the thermodynamic limit. Such a large length scale, being larger than the sizes studied before, explains previous conflicting results. We also predict spatial profiles of magnetization in steady states of generic nondiffusive systems.

Satisfiability-unsatisfiability transition in the adversarial satisfiability problem.

Adversarial satisfiability (AdSAT) is a generalization of the satisfiability (SAT) problem in which two players try to make a Boolean formula true (resp. false) by controlling their respective sets of variables. AdSAT belongs to a higher complexity class in the polynomial hierarchy than SAT, and therefore the nature of the critical region and the transition are not easily parallel to those of SAT and worthy of independent study. AdSAT also provides an upper bound for the transition threshold of the quantum satisfiability problem (QSAT). We present a complete algorithm for AdSAT, show that 2-AdSAT is in P, and then study two stochastic algorithms (simulated annealing and its improved variant) and compare their performances in detail for 3-AdSAT. Varying the density of clauses α we claim that there is a sharp SAT-UNSAT transition at a critical value whose upper bound is αc≲1.5, suggesting a much stricter upper bound for the QSAT transition than those previously found.

How many eigenvalues of a Gaussian random matrix are positive?

We study the probability distribution of the index N(+), i.e., the number of positive eigenvalues of an N×N Gaussian random matrix. We show analytically that, for large N and large N(+) with the fraction 0≤c=N(+)/N≤1 of positive eigenvalues fixed, the index distribution P(N(+)=cN,N)~exp[-ßN(2)Φ(c)] where ß is the Dyson index characterizing the Gaussian ensemble. The associated large deviation rate function Φ(c) is computed explicitly for all 0≤c≤1. It is independent of ß and displays a quadratic form modulated by a logarithmic singularity around c=1/2. As a consequence, the distribution of the index has a Gaussian form near the peak, but with a variance Δ(N) of index fluctuations growing as Δ(N)~lnN/ßπ(2) for large N. For ß=2, this result is independently confirmed against an exact finite-N formula, yielding Δ(N)=lnN/2π(2)+C+O(N(-1)) for large N, where the constant C for even N has the nontrivial value C=(γ+1+3ln2)/2π(2)≃0.185 248… and γ=0.5772… is the Euler constant. We also determine for large N the probability that the interval [ζ(1),ζ(2)] is free of eigenvalues. Some of these results have been announced in a recent letter [Phys. Rev. Lett. 103, 220603 (2009)].

Statistical properties of determinantal point processes in high-dimensional Euclidean spaces.

The goal of this paper is to quantitatively describe some statistical properties of higher-dimensional determinantal point processes with a primary focus on the nearest-neighbor distribution functions. Toward this end, we express these functions as determinants of NxN matrices and then extrapolate to N-->infinity . This formulation allows for a quick and accurate numerical evaluation of these quantities for point processes in Euclidean spaces of dimension d . We also implement an algorithm due to Hough for generating configurations of determinantal point processes in arbitrary Euclidean spaces, and we utilize this algorithm in conjunction with the aforementioned numerical results to characterize the statistical properties of what we call the Fermi-sphere point process for d=1-4 . This homogeneous, isotropic determinantal point process, discussed also in a companion paper [S. Torquato, A. Scardicchio, and C. E. Zachary, J. Stat. Mech.: Theory Exp. (2008) P11019.], is the high-dimensional generalization of the distribution of eigenvalues on the unit circle of a random matrix from the circular unitary ensemble. In addition to the nearest-neighbor probability distribution, we are able to calculate Voronoi cells and nearest-neighbor extrema statistics for the Fermi-sphere point process, and we discuss these properties as the dimension d is varied. The results in this paper accompany and complement analytical properties of higher-dimensional determinantal point processes developed in a prior paper.

Index distribution of gaussian random matrices.

We compute analytically, for large N, the probability distribution of the number of positive eigenvalues (the index N+) of a random N x N matrix belonging to Gaussian orthogonal (beta=1), unitary (beta=2) or symplectic (beta=4) ensembles. The distribution of the fraction of positive eigenvalues c=N+/N scales, for large N, as P(c,N) approximately = exp[-betaN(2)Phi(c)] where the rate function Phi(c), symmetric around c=1/2 and universal (independent of beta), is calculated exactly. The distribution has non-Gaussian tails, but even near its peak at c=1/2 it is not strictly Gaussian due to an unusual logarithmic singularity in the rate function.

Parity doubling and SU(2)LxSU(2)R restoration in the hadron spectrum.

We construct the most general nonlinear representation of chiral SU(2)LxSU(2)R broken down spontaneously to the isospin SU(2), on a pair of hadrons of same spin and isospin and opposite parity. We show that any such representation is equivalent, through a hadron field transformation, to two irreducible representations on two hadrons of opposite parity with different masses and axial-vector couplings. This implies that chiral symmetry realized in the Nambu-Goldstone mode does not predict the existence of degenerate multiplets of hadrons of opposite parity nor any relations between their couplings or masses.