2025-12-29 | | Total: 17
In the approaches based on matrix-product states (MPSs), such as the density-matrix renormalization group (DMRG) method, the ordering of the sites crucially affects the computational accuracy. We investigate the performance of an algorithm that searches for the optimal site order by iterative local site rearrangement. We improve the algorithm by expanding the range of site rearrangement and apply it to a one-dimensional quantum Heisenberg model with random site permutation. The results indicate that increasing the range of the site rearrangement significantly improves the computational accuracy of the DMRG method. In particular, increasing the rearrangement range from two to three sites reduces the average relative error in the ground-state energy by 65% to 94% in the cases we tested. We also discuss the computational cost of the algorithm and its application as a preprocessing for MPS-based calculations.
Diffusion in rugged free-energy landscapes is central to diverse problems in chemical physics, biomolecular dynamics, polymer transport and numerous disordered systems. Zwanzig's well-known classic mean-field theory predicts that roughness reduces the diffusion coefficient by an exponential factor determined solely by the variance of the disorder. However, the numerical studies of Banerjee, Seki, and Bagchi (BSB) showed that this result fails for uncorrelated Gaussian-distributed site energies because rare but deep three-site traps dominate long-time transport. BSB introduced Gaussian \emph{spatial} correlations - originally developed in astrophysics to model turbulent density fluctuations - and demonstrated that even modest correlations suppress these pathological traps and restore Zwanzig's exponential scaling. Here we present here a unified theoretical framework clarifying (i) why Zwanzig's local averaging, may be viewed as a Gaussian cumulant expansion, could breakdown. In particular, how its validity is destroyed by uncorrelated disorder, and (ii) how Gaussian spatial correlations reshape roughness increments, eliminate asymmetric multi-site traps, and thereby recover mean-field diffusion, and (iii) a derivation showing exactly how Gaussian spatial correlations modify roughness increments, trap statistics, and ultimately the diffusion constant. We also provide explicit numerical triplet examples illustrating the dramatic reduction of escape times produced by spatial correlations.
Complex behaviors often sit at a critical threshold between order and disorder. But not all disorder is created equal. Disorder can be trivial or constrained, and correlated disorder can even be topological. Crucially, constrained disorder can harbor memory, leading to non-trivial, sequence-dependent responses to external manipulations. And yet the fascinating subject of "memory of disorder" remains poorly explored, as memory is often associated to the retention of metastable order. In recent years artificial frustrated materials -- in particular arrays of frustrated nanomagnets known as artificial spin ices -- have been employed to study complex disorders and its wealth of exotic behaviors, yet their memory properties have received much less attention. Here, we investigate both analytically and numerically the sequence-dependent responses of two somehow opposite yet related artificial spin ices: the Landau-ordered square spin ice and the disordered but topologically-ordered Shakti spin ice. We find that Shakti exhibits a pronounced sequence-dependent response, whereas in the square lattice, such path dependence is absent. Within Shakti, even the minimal periodic supercell demonstrates both deterministic and stochastic forms of sequence memory, depending on the interaction strength. Extending our study to cyclic driving, we find that retracing the same input path leads to enhanced memory retention. These results open new perspectives on how topological constraints and correlated disorder generate robust memory effects in frustrated artificial materials, hitherto examined mainly in terms of their ground-state kinetics and thermodynamics.
In this paper, we have studied the critical temperature $T_c$ of continuous spin $2d$ square-lattice Ising model using Monte-Carlo simulation. We have considered spins $s$ in a bounded interval, where $s \in [-1,+1]$ in square-lattice configuration with periodic boundary condition. We have observed that the critical temperature $T_c$ is approximately $0.925$, showing a clear second order phase transition. Considering finite size scaling, we have also obtained the critical exponents associated with susceptibility, specific heat, magnetization and we find that these values are in good agreement with the corresponding values obtained for the standard $2d$ Ising universality class.
This study is presented in a series of two papers. The first paper is an analytical study. This is the second paper, and here we numerically study the thermal energy harvesting capability of two electronic circuits. The first circuit consists of a diode and capacitor in series. We solve the time-dependent Fokker-Planck equation and show the capacitor initially charges and then discharges to zero. The peak charge on the capacitor increases with temperature, capacitance, and diode quality. The second circuit has two current loops with one small capacitor, two storage capacitors, and two diodes wired in opposition. When the diodes are held at different temperatures we observe a non-zero steady-state charge is accumulated on both storage capacitors. The magnitude of the stored charges are nearly equal but the signs are opposite. When resistors are used in place of diodes there is no transient and no steady-state charge buildup. Numerical studies for the time-independent Fokker-Planck equation are presented and confirm the steady state charges.
Nonlinear elements in a rectifying circuit can be used to harvest energy from thermal fluctuations either steadily or transitorily. We study an energy harvesting system comprising a small variable capacitor (e.g., free standing graphene) wired to two diodes and two storage capacitors that may be kept at different temperatures (or at a single one) and use two current loops. The system reaches very rapidly a quasi stationary state with constant overall charge while the difference of the charges at the storage capacitors evolves much more slowly to its stationary value. In this paper, we extract an exponentially small factor out of the solution of the Fokker-Planck equation and use a Chapman-Enskog procedure to describe the long evolution of the marginal probability density for the charge difference, from the quasi stationary state to the final stationary state (thermal equilibrium for equal temperatures). The second paper of this series shows that the results of the perturbation procedure compare well with direct numerical simulations. For a specific form of the diodes' nonlinear mobilities, we can approximate the quasi stationary state by Gaussian functions and further study the evolution of the marginal probability density. The latter adopts the shape of a slowly expanding pulse (comprising left and right moving wave fronts whose fore edges become sharper as time elapses) in the space of charge differences that leaves the final stationary state behind it.
Quantum critical phases are extended regions of phase space characterized by a diverging correlation length. By analogy, we define an information critical phase as an extended region of a mixed state phase diagram where the Markov length, the characteristic length scale governing the decay of the conditional mutual information (CMI), diverges. We demonstrate that such a phase arises in decohered $\mathbb{Z}_{N}$ Toric codes by assessing both the CMI and the coherent information, the latter quantifying the robustness of the encoded logical qudits. For $N>4$, we find that the system hosts an information critical phase intervening between the decodable and non-decodable phases where the coherent information saturates to a fractional value in the thermodynamic limit, indicating that a finite fraction of logical information is still preserved. We show that the density matrix in this phase can be decomposed into a convex sum of Coulombic pure states, where gapped anyons reorganize into gapless photons. We further consider the ungauged $\mathbb{Z}_{N}$ Toric code and interpret its mixed state phase diagram in the language of strong-to-weak spontaneous symmetry breaking. We argue that in the dual model, the information critical phase arises because the spontaneously broken off-diagonal $\mathbb{Z}_{N}$ symmetry gets enhanced to a U(1) symmetry, resulting in a novel superfluid phase whose gapless modes involve coherent excitations of both the system and the environment. Finally, we propose an optimal decoding protocol for the corrupted $\mathbb{Z}_{N}$ Toric code and evaluate its effectiveness in recovering the fractional logical information preserved in the information critical phase. Our findings identify a gapless analog for mixed-state phases that still acts as a fractional topological quantum memory, thereby extending the conventional paradigm of quantum memory phases.
We investigate the statistical physics of quantum error correction in ${\rm U}(1)$ symmetry-enriched topological quantum memories. Starting from a phenomenological error model of charge-conserving noise, we study the optimal decoder assuming the local charges of each anyon can be measured. The error threshold of the optimal decoder corresponds to a continuous phase transition in a disordered two-dimensional integer loop model on the Nishimori line. Using an effective replica field theory analysis and Monte Carlo numerics, we show that the optimal decoding transition exhibits Berezinskii-Kosterlitz-Thouless universality with a modified universal jump in winding number variance. We further generalize the model beyond the Nishimori line, which defines a large class of suboptimal decoders. At low nonzero temperatures and strong disorder, we find numerical evidence of a disorder-dominated loop-glass phase which corresponds to a "confidently incorrect" decoder. The zero-temperature limit defines the minimum-cost flow decoder, which serves as the ${\rm U}(1)$ analog of minimum-weight perfect matching in $\mathbb{Z}_2$ topological codes. Both the optimal and minimum-cost flow decoders are shown to dramatically outperform the charge-agnostic optimal decoder in symmetry-enriched topological codes.
We present an interacting branching model of neural network dynamics, incorporating key biological features such as inhibition with several types of inhibitory interactions. We establish a hierarchy of analytical mean-field approximations to the model, which characterizes nonequilibrium phase transitions between disorder and ordered phases, and perform a stability analysis. Generically, inhibitory neurons increase the stability of the model dynamics. The model is consistent with the quasi-criticality hypothesis in that it displays regions of maximal dynamical susceptibility and maximal mutual information predicated on the strength of the external stimuli. Directed percolation emerges as the universality class of the critical transition of the model, consistent with some previous experimental data and models. In the unstable phase, chaotic dynamics emerge, which may be linked to the occurrence of epileptic seizures.
A powerful method for analyzing quantum error-correcting codes is to map them onto classical statistical mechanics models. Such mappings have thus far mostly focused on static codes, possibly subject to repeated syndrome measurements. Recent progress in quantum error correction, however, has prompted new paradigms where codes emerge from stabilizer circuits in spacetime -- a unifying perspective encompassing syndrome extraction circuits of static codes, dynamically generated codes, and logical operations. We show how to construct statistical mechanical models for stabilizer circuits subject to independent Pauli errors, by mapping logical equivalence class probabilities of errors to partition functions using the spacetime subsystem code formalism. We also introduce a modular language of spin diagrams for constructing the spin Hamiltonians, which we describe in detail focusing on independent circuit-level X-Z error channels. With the repetition and toric codes as examples, we use our approach to analytically rank logical error rates and thresholds between code implementations with standard and dynamic syndrome extraction circuits, describe the effect of transversal logical Clifford gates on logical error rates, and perform Monte Carlo simulations to estimate maximum likelihood thresholds. Our framework offers a universal prescription to analyze, simulate, and compare the decoding properties of any stabilizer circuit, while revealing the innate connections between dynamical quantum systems and noise-resilient phases of matter.
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 (ED), we show that the model with spin $S\ge2$ 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 $\mathcal{O}(2^L)$ number of SSB ground states originating from the $\mathcal{O}(2^L)$ number of exponential $U(1)$ charge sectors, leading to a finite entropy density $\ln 2$. This enforces a first order SSB phase transition into this phase, as observed in ED 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 non-decaying 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 consider both unitary and nonunitary A-D-E minimal models on the cylinder with topological defects along the non-contractible cycle of the cylinder. We define the coset graph $A \otimes G/\mathbb{Z}_2$ and argue that it encodes not only the (i) coset graph fusion algebra, but also (ii) the Affleck-Ludwig boundary g-factors; (iii) the defect g-factors (quantum dimensions) and (iv) the relative symmetry resolved entanglement entropy. By studying A-D-E restricted solid-on-solid models, we find that these boundary conformal field theory structures are also present on the lattice: defects (seams) are implemented by face weights with special values of the spectral parameter. Integrability allows the study of lattice transfer matrix T- and Y-system functional equations to reproduce the fusion algebra of defect lines. The effective central charges and conformal weights are expressed in terms of dilogarithms of the braid and bulk asymptotics of the Y-system expressed in terms of the quantum dimensions.
We formulate a time-dependent Fluctuating Local Field (TD-FLF) method for correlated fermion dynamics, extending the stationary FLF approach. The wavefunction is approximated as an ensemble of non-interacting states subject to a classical fluctuating field, with dynamics encoded in the field's time-dependent distribution. This reduces the time-dependent Schrödinger equation to a generalized eigenvalue problem in a significantly reduced basis. Applied to half-filled 2D Hubbard lattices, TD-FLF yields highly accurate results, outperforming mean-field theory and capturing oscillation frequencies and amplitudes in good agreement with exact diagonalization. Its low computational cost and flexibility make TD-FLF a promising tool for simulating driven correlated systems.
We consider a one-dimensional mixture of bosons and spinless fermions with contact interactions. In this system, the excitations at low energies are described by four linearly dispersing modes characterized by two excitation velocities. Here we study the velocities in a system with equal interaction strengths and equal masses of bosons and fermions. The resulting model is integrable and admits an exact Bethe-ansatz solution. We analyze it and analytically derive various exact results, which include the Drude weight matrix. We show that the excitation velocities can be calculated from the knowledge of the matrices of compressibility and the Drude weights, as their squares are the eigenvalues of the product of the two matrices. The elements of the Drude weight matrix obey certain sum rules as a consequence of Galilean invariance. Our results are consistent with the presence of a momentum-momentum coupling term between the two subsystems of bosons and fermions in the effective low-energy Hamiltonian. The analytical method developed in the present study can be extended to other models that possess a nested Bethe-ansatz structure.
We demonstrate the existence of a pair of almost dissipationless oscillating modes at low temperatures in both the shear and sound channels of a hybrid quantum system, comprised of a weakly self-interacting perturbative sector coupled to strongly self-interacting holographic degrees of freedom described by a black hole geometry. We argue that these modes realize prethermal time-crystal behavior in semiholographic systems without fine-tuning and can be observed by measuring operators that probe either the hard (perturbative) or the soft (holographic) sector. We also find novel Gregory Laflamme type instabilities that lead to the formation of inhomogeneities at higher temperatures. These results provide evidence that black holes with planar horizons and dynamical boundary conditions can develop both inhomogeneous and metastable time crystal phases over a wide range of temperatures set by an intermediate scale. Furthermore, they suggest that such phases can be realized without external driving in nonAbelian plasmas of asymptotically free gauge theories in the large-$N$ limit.
The notion of quantum symmetry has recently been extended to include reduced-dimensional transformations and algebraic structures beyond groups. Such generalized symmetries lead to exotic phases of matter and excitations that defy Landau's original paradigm. Here, we develop an algebraic framework for systematically deriving the fusion rules of topological defects in higher-dimensional lattice systems with non-invertible generalized symmetries, and focus on a (2+1)-dimensional quantum Ising plaquette model as a concrete illustration. We show that bond-algebraic automorphisms, when combined with the so-called half-gauging procedure, reveal the structure of the non-invertible duality symmetry operators, which can be explicitly represented as a sequential quantum circuit. The resulting duality defects are constrained by the model's rigid higher symmetries (lower-dimensional subsystem symmetries), leading to restricted mobility. We establish the fusion algebra of these defects. Finally, in constructing the non-invertible duality transformation, we explicitly verify that it acts as a partial isometry on the physical Hilbert space, thereby satisfying a recent generalization of Wigner's theorem applicable to non-invertible symmetries.
While quantum simulators promise to explore quantum many-body physics beyond classical computation, their capabilities are limited by the available native interactions in the hardware. On many platforms, accessible Hamiltonians are largely restricted to one- and two-body interactions, limiting access to multi-body Hamiltonians and to systems governed by simultaneous, non-commuting interaction terms that are central to condensed matter, quantum chemistry, and high-energy physics. We introduce and experimentally demonstrate a hybrid digital-analog protocol that overcomes these limitations by embedding analog evolution between shallow entangling-gate layers. This method produces effective Hamiltonians with simultaneous non-commuting three- and four-body interactions that are generated non-perturbatively and without Trotter error -- capabilities not practically attainable on near-term hardware using purely digital or purely analog schemes. We implement our scheme on a trapped-ion quantum processor and use it to realize a topological spin chain exhibiting prethermal strong zero modes persisting at high temperature, as well as models featuring three- and four-body interactions. Our hardware-agnostic and scalable method opens new routes to realizing complex many-body physics across quantum platforms.