We propose a quantum model of fermions simulating the electrical breakdown process of dielectrics. The model consists of
M sites with
N fermion modes per site and has a conserved charge
Q. It has an on-site chemical potential $\ensuremath{\mu}$ with disorder
W and an interaction of strength
J restricting each fermion to excite two more fermions when moving forward by one site. We show that the
N=3 model with disorder
W=0 shows a Hilbert space fragmentation and is exactly solvable except for very few Krylov subspaces. The analytical solution shows that the
N=3 model exhibits many-body localization (MBL) as $M\ensuremath{\rightarrow}\ensuremath{\infty}$, which is stable against
W>0 as our exact diagonalization (ED) shows. At
N>3, our ED suggests an MBL to quantum chaos crossover at small
W as
M/N decreases across 1, and persistent MBL at large
W. At
W=0, an exactly solvable many-body scar flat band exists in many charge
Q sectors, which has a nonzero measure in the thermodynamic limit. We further calculate the time evolution of a fermion added to the particle vacuum, which shows a breakdown (dielectric) phase when $\ensuremath{\mu}/J<1/2$ ($\ensuremath{\mu}/J>1/2$) if $W\ensuremath{\ll}J$, and no breakdown if $W\ensuremath{\gg}J$.
We study the role of particle-hole symmetry on the universality class of various quantum phase transitions corresponding to the onset of superfluidity at zero temperature of bosons in a quenched random medium. To obtain a model with an exact particle-hole symmetry it is necessary to use the Josephson junction array, or quantum rotor, Hamiltonian, which may include disorder in both the site energies and the Josephson couplings between wave function phase operators at different sites. The functional integral formulation of this problem in
d spatial dimensions yields a
(d+1)-dimensional classical
XY model with extended disorder, constant along the extra imaginary time dimension---the so-called random rod problem. Particle-hole symmetry may then be broken by adding nonzero site energies, which may be uniform or site dependent. We may distinguish three cases: (i) exact particle-hole symmetry, in which the site energies all vanish; (ii) statistical particle-hole symmetry, in which the site energy distribution is symmetric about zero, vanishing on average; and (iii) complete absence of particle-hole symmetry in which the distribution is generic. We explore in each case the nature of the excitations in the nonsuperfluid Mott insulating and Bose glass phases. We show, in particular, that, since the boundary of the Mott phase can be derived exactly in terms of that for the pure, nondisordered system, there can be no direct Mott-superfluid transition. Recent Monte Carlo data to the contrary can be explained in terms of rare region effects that are inaccessible to finite systems. We find also that the Bose glass compressibility, which has the interpretation of a temporal spin stiffness or superfluid density, is positive in cases (ii) and (iii), but that it vanishes with an essential singularity as full particle-hole symmetry is restored. We then focus on the critical point and discuss the relevance of type (ii) particle-hole symmetry-breaking perturbations to the random rod critical behavior, identifying a nontrivial crossover exponent. This exponent cannot be calculated exactly but is argued to be positive and the perturbation therefore relevant. We argue next that a perturbation of type (iii) is irrelevant to the resulting type (ii) critical behavior: The statistical symmetry is restored on large scales close to the critical point, and case (ii) therefore describes the dirty boson fixed point. Using various duality transformations we verify all of these ideas in one dimension. To study higher dimensions, we attempt, with partial success, to generalize the Dorogovtsev--Cardy--Boyanovsky double-epsilon expansion technique to this problem. We find that when the dimension of time ${ϵ}_{\ensuremath{\tau}}<{ϵ}_{\ensuremath{\tau}}^{c}\ensuremath{\simeq}\frac{8}{29}$ is sufficiently small a type (ii) symmetry-breaking perturbation is irrelevant, but that for sufficiently large ${ϵ}_{\ensuremath{\tau}}>{ϵ}_{\ensuremath{\tau}}^{c}$ particle-hole asymmetry is a relevant perturbation and a new stable fixed point appears. Furthermore, for ${ϵ}_{\ensuremath{\tau}}>{ϵ}_{\ensuremath{\tau}}^{c2}\ensuremath{\approx}\frac{2}{3}$, this fixed point is stable also to perturbations of type (iii): at $ϵ={ϵ}_{\ensuremath{\tau}}^{c2}$ the generic type (iii) fixed point merges with the new fixed point. We speculate, therefore, that this new fixed point becomes the dirty boson fixed point when ${ϵ}_{\ensuremath{\tau}}=1$. We point out, however, that ${ϵ}_{\ensuremath{\tau}}=1$ may be quite special. Thus, although the qualitative renormalization group flow picture the double-epsilon expansion technique provides is quite compelling, one should remain wary of applying it quantitatively to the dirty boson problem.
NOTE: Text or symbols not renderable in plain ASCII are indicated by [...]. Abstract is included in .pdf document. We study the nature of various quantum phase transitions corresponding to the onset of superfluidity, at zero temperature, of bosons in a quenched medium. Particle-hole symmetry plays an essential role in determining the universality class of the transitions. To obtain a model with an exact particle-hole symmetry it is necessary to use the Josephson junction array Hamiltonian, which may include disorder in the Josephson couplings between phases at different sates. The functional integral formulation of this problem in d spatial dimensions yields a (d + 1)-dimensional classical XY-model with extended disorder, constant along the extra imaginary time dimension -- the so-called random rod problem. Particle-hole symmetry may then be broken by adding nonzero chemical potentials or site energies, which may also be site dependent and random. We may then distinguish three cases: (i) exact particle-hole symmetry, in which the site energies all vanish, (ii) statistical particle-hole symmetry in which the site energy distribution is symmetric about zero and hence vanishes on average, and (iii) complete absence of particle-hole symmetry in which the distribution is generic. We explore in each case the nature of the excitations in the nonsuperfluid Pose glass phase. We find, for example, that the compressibility, which has the interpretation of a temporal spin stiffness or superfluid density, is positive in cases (ii) and (iii), but that it vanishes with an essential singularity as full particle-hole symmetry is restored. We then focus on the critical point and discuss the validity of various scaling arguments. In particular, we argue that the dynamical exponent z could be different from d, and the arguments leading to their equality are incorrect. We then discuss the relevance of a type (ii) particle-hole symmetry breaking perturbation to the random rod critical behavior, identifying a nontrivial crossover exponent. This exponent cannot be calculated exactly but is argued to be positive and the symmetry breaking perturbation therefore relevant. We argue next that a perturbation of type (iii) is irrelevant to the resulting type (ii) critical behavior: the statistical symmetry is restored on large scales close to the critical point, and case (ii) therefore describes the dirty boson fixed point. Using various duality transformations we verify all of these ideas in one dimension. To study higher dimensions we attempt, with partial success, to generalize the Dorogovtsev-Cardy-Boyonovsky double epsilon expansion technique to this problem. We find that when the dimension of time [...] is sufficiently small the symmetry breaking perturbation of type (ii) is irrelevant, but that for sufficiently large [...] this is a relevant perturbation and a new stable commensurate fixed point appears. We speculate that this new fixed point becomes the dirty boson fixed point when [...] = 1. We also show that for [...], there exists a particle-hole asymmetric fixed point of type (iii), but we provide evidence that it merges with the commensurate fixed point for some finite [...]. This tends to confirm symmetry restoration at the physical [...] = 1.
We study the phase diagram of a one-dimensional spin quantum breakdown model, which has an exponential U(1) symmetry with charge unit decaying as 2^{-j} with site position j. By exact diagonalization and density matrix renormalization group, we show that the model with spin S≥2 exhibits an exponential U(1) spontaneous symmetry-breaking (SSB) phase dubbed a quantum breakdown condensate. It exhibits a bulk gap violating the Goldstone theorem and an edge mode only on the left edge if in open boundary condition. In a length L lattice, the condensate has O(2^{L}) number of SSB ground states originating from the O(2^{L}) number of exponential U(1) charge sectors, leading to a finite entropy density ln2. This enforces a first-order SSB phase transition into this phase, as observed numerically and verified in the large S limit on an exactly solvable Rokhsar-Kivelson line. The condensate has an SSB order parameter being the local in-plane spin, which points in angles related by the chaotic Bernoulli (dyadic) map and thus is effectively random. Moreover, we show the condensate exhibits nondecaying local autocorrelations and does not have an off-diagonal long-range order. The quantum breakdown condensate thus behaves as a disorder-free quantum glass and is beyond the existing classifications of phases of matter.
We introduce and study a disorder-free version of the quantum breakdown model with all-to-all interactions. The Hamiltonian factorizes into the product of the zero-momentum-mode occupation number and a quadratic Hamiltonian including only pairing terms. This structure makes the model exactly solvable and produces an extensive zero-energy degeneracy. We obtain the spectrum and partition function in closed form and analyze the thermodynamics, spectral form factor, two-point functions, and a regulated out-of-time-ordered correlator (OTOC). The OTOC is time independent at leading order in the large-
N limit, while its time-dependent contributions first appear at order
1/N. At infinite temperature, the time-dependent correction reduces to a two-point function arising from the zero-momentum component of the fermion operators. The model therefore provides a controlled setting for isolating the consequences of the breakdown interaction in the absence of disorder, spatial structure, and environmental coupling.