2025-05-02 | | Total: 16
The general epidemic process is a paradigmatic model in non-equilibrium statistical physics displaying a continuous phase transition between active and absorbing states.The dynamic isotropic percolation universality class captures its universal properties, which we aim to quantitatively study by means of the field-theoretic formulation of the model augmented with a perturbative renormalization group analysis. The main purpose of this work consists in determining the critical dynamic exponent z to the three-loop approximation. This allows us to finalize the quantitative description of the dynamic isotropic percolation class to this order of perturbation theory. The calculations are performed within the dimensional regularization with the minimal subtraction scheme and actual perturbative expansions are carried out in a formally small parameter ϵ, where ϵ=6−d is a deviation from the upper critical dimension dc=6.
We apply the operator approach to a stochastic system belonging to a class of death-birth processes, which we introduce utilizing the master equation approach. By employing Doi- Peliti formalism we recast the master equation in the form of a Schrödinger-like equation. Therein appearing pseudo-Hamiltonian is conveniently expressed in a suitable Fock space, constructed using bosonic-like creation and annihilation operators. The kernel of the associated time evolution operator is rewritten using a functional integral, for which we propose an approximate method that allows its analytical treatment. The method is based on the expansion in eigenfunctions of the Hamiltonian generating given functional integral. In this manner, we obtain approximate values for the probabilities of the system being in the first and second states for the case of the pure birth process.
We consider the effect of non-reciprocity in a binary mixture of self-propelled particles with anti-aligning interactions, where a particle of type A reacts differently to a particle of type B than vice versa. Starting from a well-known microscopic Langevin-model for the particles, setting up the corresponding exact N-particle Fokker-Planck equation and making Boltzmann's assumptions of low density and one-sided molecular chaos, the non-linear active Boltzmann equation with the exact collision operator is derived. In this derivation, the effect of phase-space compression and the build-up of pair-correlations during binary interactions is explicitly taken into account, leading to a theoretical description beyond mean-field. This extends previous results for reciprocal interactions, where it was found that orientational order can emerge in a system with purely anti-aligning interactions. Although the equations of motion are more complex than in the reciprocal system, the theory still leads to analytical expressions and predictions. Comparisons with agent-based simulations show excellent quantitative agreement of the dynamic and static behavior in the low density and/or small coupling limit.
In open systems with strong coupling, the interaction energy between the system and the environment is significant, so thermodynamic quantities cannot be reliably obtained by traditional statistical mechanics methods. The Hamiltonian of mean force H∗β offers an in principle accurate theoretical basis by explicitly accounting for the interaction energy. However, calculating the Hamiltonian of mean force is challenging both theoretically and computationally. We demonstrate that when the condition VarE0(e−βVSE)=0 is met, the dependence of thermodynamic variables can be shifted from {Pβ(xS),H∗β(xS)} to {Pβ(xS),P(VSE)}. This change simplifies thermodynamic measurements. As a central result, we derive a general equality that holds for arbitrary coupling strengths and from which an inequality follows - aligned with Jensen's inequality applied to the Gibbs-Bogoliubov-Feynman bound. This equality, analogous in importance to the Jarzynski equality, offers deeper insight into free energy differences in strongly coupled systems. Finally, by combining our result with said Jarzynski equality, we derive additional relations that further clarify thermodynamic behavior in strongly coupled open systems.
Bose Einstein condensate of optically inactive (dark) intervalley excitons must be inherently multi-component because of the point group and the time-reversal symmetries of a crystal. Accordingly, a number of the condensate components Nv is determined by the symmetry. Since the valleys hosting such excitons are separated by large quasi-momenta, the minimal inter-component Josephson-type coupling can only be established between pairs of excitons from the time-reversed valleys. As a result, a paired condensate can emerge at finite temperature, that is, the phase where individual valleys are not characterized by the condensate order, while the order exists for the pairs from the time-reversed valleys. This prediction follows from the elementary mean field analysis regardless of the dimensionality. However, as Monte Carlo simulations show, no such a phase exists in 3D crystals. Instead, the Nv-component condensation proceeds as the Ist order transition from the normal state. The paired phase does exist in 2D for Nv≥6. It forms by Berezinskii-Kosterlitz-Thouless transition from the high temperature (normal) phase. Upon further lowering temperature, a second transition transforms the paired phase into the Nv-component condensate.
We investigate the out-of-equilibrium dynamics of quantum information in one-dimensional systems undergoing a quantum quench using a local perspective based on the information lattice. This framework provides a scale- and space-resolved decomposition of quantum correlations, enabling a hydrodynamic description of the information flow through well-defined local densities -- termed local information -- and currents. We apply this framework to three local quenches in noninteracting fermionic chains: (i) the release of a single particle into an empty tight-binding chain, (ii) the connection of two critical chains via the removal of a central barrier, and (iii) the coupling of a topological Kitaev chain to a critical chain. In each case, the information lattice reveals the local structure of correlation buildup and information interface effects, going beyond global measures such as the von Neumann entropy. In particular, through the information lattice we uncover the signatures in the local information flow associated with topological edge modes and analytically explain the fractional von Neumann entropy values observed in Majorana quench protocols. Our approach is general and applicable to interacting, disordered, and open systems, providing a powerful tool for characterizing quantum information dynamics.
In this paper, we consider the quantum Langevin equation for the Caldeira-Leggett model with an arbitrary time-dependent coupling constant. We solve this equation exactly by employing a train of Dirac-delta switchings. This method also enables us to visualize the memory effect in the environment. Furthermore, we compute the two-time correlation functions of the system's quadratures and show that the discrete-time Fourier transform is well-suited for defining spectral densities, as the Dirac-delta switchings turn continuous functions into discretized samples.
The design of superionic conductors for all-solid-state batteries often faces a fundamental trade-off between stability and ionic conductivity. Random Substitutional Crystals (RSCs), where atomic species are randomly distributed throughout a crystal lattice, present a promising route to overcome this competitive relation. Although extensive studies have focused on local ionic hopping, the role of mesoscale structural organization in determining macroscopic conductivity remains poorly understood, limiting the rational design of optimal compositions. Here, we systematically investigate the ionic conductivity of NaCl-type RSCs as a function of composition using molecular dynamics simulations. We find that ionic conductivity increases sharply once the carrier ion concentration exceeds a critical threshold, without disrupting the underlying crystal structure. Strikingly, this threshold aligns with the site percolation threshold predicted by percolation theory. Our findings establish ion percolation as a universal design principle that reconciles the trade-off between conductivity and stability, offering a simple and broadly applicable strategy for the development of robust, high-performance solid electrolytes.
We introduce a unified statistical framework for quantifying system-environment coupling by treating the interaction energy VSE as a stochastic variable. Using a reference-particle decomposition, we derive exact, closed-form expressions for the mean and variance of VSE in terms of the single-particle density and up to four-body correlation functions. When VSE is approximately Gaussian, these two moments suffice to compute the free energy shift of the strongly coupled system. To validate our framework, we ran explicit Monte Carlo simulations of the full system-environment configurations across a range of system sizes, generating reference distributions of the interaction energy VSE. We then applied our derived analytical formulas to predict these distributions and found excellent agreement in both the weak- and strong-coupling regimes.
For disordered systems within the random first-order transition (RFOT) universality class, such as structural glasses and certain spin glasses, the role played by activated relaxation processes is rich to the point of perplexity. Over the last decades, various efforts have attempted to formalize and systematize such processes in terms of instantons similar to the nucleation droplets of first- order phase transitions. In particular, Kirkpatrick, Thirumalai, and Wolynes proposed in the late '80s an influential nucleation theory of relaxation in structural glasses. Already within this picture, however, the resulting structures are far from the compact objects expected from the classical droplet description. In addition, an altogether different type of single-particle hopping-like instantons has recently been isolated in molecular simulations. Landscape studies of mean-field spin glass models have further revealed that simple saddle crossing does not capture relaxation in these systems. We present here a landscape-agnostic study of rare dynamical events, which delineates the richness of instantons in these systems. Our work not only captures the structure of metastable states, but also identifies the point of irreversibility, beyond which activated relaxation processes become a fait accompli. An interpretation of the associated landscape features is articulated, thus charting a path toward a complete understanding of RFOT instantons.
It is known that Green's functions can be expressed as continued fractions; the content at the n-th level of the fraction is encoded in a coefficient bn, which can be recursively obtained using the Lanczos algorithm. We present a theory concerning errors in approximating Green's functions using continued fractions when only the first N coefficients are known exactly. Our focus lies on the stitching approximation (also known as the recursion method), wherein truncated continued fractions are completed with a sequence of coefficients for which exact solutions are available. We assume a now standard conjecture about the growth of the Lanczos coefficients in chaotic many-body systems, and that the stitching approximation converges to the correct answer. Given these assumptions, we show that the rate of convergence of the stitching approximation to a Green's function depends strongly on the decay of staggered subleading terms in the Lanczos cofficients. Typically, the decay of the error term ranges from 1/poly(N) in the best case to 1/poly(logN) in the worst case, depending on the differentiability of the spectral function at the origin. We present different variants of this error estimate for different asymptotic behaviours of the bn, and we also conjecture a relationship between the asymptotic behavior of the bn's and the smoothness of the Green's function. Lastly, with the above assumptions, we prove a formula linking the spectral function's value at the origin to a product of continued fraction coefficients, which we then apply to estimate the diffusion constant in the mixed field Ising model.
We show a surprising relation between quantum learning theory and algorithmic hardness. We demonstrate that finding near-ground states of certain sparse disordered quantum systems is average-case hard for "Lipschitz" quantum algorithms if there exists an efficient, local learning algorithm -- such as the classical shadows algorithm -- for estimating the energy of a state of the system. A corollary of our result is that many standard quantum algorithms fail to find near-ground states of these systems, including short-time Lindbladian dynamics, short-time quantum annealing, phase estimation, and shallow-depth variational quantum algorithms. To achieve this, we introduce a topological property of quantum systems that we call the quantum overlap gap property (QOGP). This property is only satisfied by systems with an efficient local learning algorithm for the energy. We prove that systems which exhibit this topological property in their low-energy space are intractable for quantum algorithms whose outputs are stable under perturbations to their inputs. We then prove that the QOGP is satisfied for a sparsified variant of the quantum p-spin model, giving the first known algorithmic hardness-of-approximation result for quantum algorithms in finding the ground state of a non-stoquastic, noncommuting quantum system. Our resulting lower bound for quantum algorithms optimizing this model using Lindbladian evolution matches the best-known time lower bound for classical Langevin dynamics optimizing classical p-spin models. For this reason we suspect that finding ground states of typical quantum p-spin models using quantum algorithms is, in practice, as intractable as the classical p-spin model is for classical algorithms. Inversely, we show that the Sachdev--Ye--Kitaev (SYK) model does not exhibit the QOGP, consistent with previous evidence that the model is rapidly mixing at low temperatures.
We numerically study ZN lattice gauge theories in 4D as prototypical models of systems with ZN 1-form symmetry. For N≥3, we provide evidence that such systems exhibit not only the expected phases with spontaneously broken/restored symmetry but also a third photon phase. When present, the 1-form symmetry provides a precise notion of confinement, and it is commonly believed that confinement ensues due to the proliferation of extended, string-like objects known as center vortices, which carry a ZN flux. However, this picture is challenged by the three-phase scenario investigated here. We show that both the confined and the photon phases are associated with the proliferation of center vortices and that the key difference between them lies in whether or not vortex-junctions, the monopoles, proliferate.
Matrix product states are useful representations for a large variety of naturally occurring quantum states. Studying their typical properties is important for understanding universal behavior, including quantum chaos and thermalization, as well as the limits of classical simulations of quantum devices. We show that the usual ensemble of sequentially generated random matrix product states (RMPS) using local Haar random unitaries is not uniform when viewed as a restriction of the full Hilbert space. As a result, the entanglement across the chain exhibits an anomalous asymmetry under spatial inversion. We show how to construct an unbiased measure starting from the left-canonical form and design a Metropolis algorithm for sampling random states. Some properties of this new ensemble are investigated both analytically and numerically, such as the resulting resolution of identity over matrix product states and the typical entanglement spectrum, which is found to differ from the sequentially generated case.
We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with N qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of N non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one w0=1, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths p≳, and all decoherence strengths. We also show that a generic initial condition w_0 >1 leads to a metastable state that eventually collapses to the w_0=1 case after a lifetime \sim \log(N/w_0). Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.
Metals at the brink of electronic quantum phase transitions display high-temperature superconductivity, competing orders, and unconventional charge transport, revealing strong departures from conventional Fermi liquid behavior. Investigation of these fascinating intertwined phenomena has been at the center of research across a variety of correlated materials over the past many decades. A ubiquitous experimental observation is the emergence of a universal timescale that governs electrical transport and momentum relaxation. In this work, we analyze an equally important theoretical question of how the energy contained in the electronic degrees of freedom near a quantum phase transition relaxes to the environment via their coupling to acoustic phonons. Assuming that the bottleneck for energy dissipation is controlled by the coupling between electronic degrees of freedom and acoustic phonons, we present a universal theory of the temperature dependence of the energy relaxation rate in a marginal Fermi liquid. We find that the energy relaxation rate exhibits a complex set of temperature-dependent crossovers controlled by emergent energy scales in the problem. We place these results in the context of recent measurements of the energy relaxation rate via non-linear optical spectroscopy in the normal state of hole-doped cuprates.