2026-10-20 | | Total: 52
We show that quantum systems on non-amenable lattices, whose boundaries constitute an extensive fraction of the system, host novel forms of bulk ordering with no Euclidean counterpart. Combining boundary-retaining iDMRG with an exact symmetry-fixing duality, we investigate $p<5$ state clock models on the coordination-3 infinite Cayley tree. We find that between the ordered and disordered phases lies an emergent intermediate phase with bulk order induced by an arbitrarily weak boundary field rather than by conventional spontaneous symmetry-breaking. This establishes boundary-sensitive order as a distinct quantum ordering mechanism defying Landau's spontaneous symmetry-breaking paradigm. The intermediate phase realizes a bulk-boundary correspondence in which the unpolarized boundary selects a specific ordered bulk state.
Near equilibrium, the Boltzmann distribution is corrected by the transient irreversible entropy flux, a result known as the McLennan formula. By expanding around equiaccessible steady states instead, we derive an analogue that can hold arbitrarily far from equilibrium. At zeroth order, steady-state occupation is determined solely by local escape rates, yielding a strong form of the active-matter principle of "low rattling." The first-order correction measures the impact of traffic imbalances on the accessibility of each state, marking the onset of nonlocal determinants of relative occupation.
We investigate quantum-state texture during the thermal relaxation of an anisotropic spin-$1/2$ Heisenberg dimer subject to a longitudinal magnetic field and a Dzyaloshinskii-Moriya (DM) interaction. Starting from local transverse and longitudinal couplings to a bosonic bath, we construct a global weak-coupling Lindblad equation, identify its transition channels, and discuss the applicability and limitations of the secular approximation. The Gibbs state is obtained explicitly in the local spin basis, revealing that the DM coupling controls both the modulus and phase of the single-excitation coherence by modifying the energy splitting and rotating the chiral phase. For heating and cooling protocols, we compare rugosity, fidelity-based texture, and trace-distance texture with coherence, vector spin chirality, and purity. The relaxation dynamics separates population thermalization from the decay and precession of deliberately prepared energy-basis coherence. Their combined contribution to local coherence produces DM-dependent texture overshoots and, during cooling, damped oscillations and sign changes of the coherence quadratures. These transients differ from the evolution of coherence magnitude and purity, while trace-distance texture weights them differently from the overlap-based quantifiers. Calculations at reduced dissipation preserve the texture maxima within the secular model. The results identify a phase-sensitive diagnostic of relaxation.
We study the phase diagram of a three-dimensional Heisenberg model with Dzyaloshinski--Moriya interaction under the presence of Potts-like anisotropy. The phase transition between the winding phase and the ferromagnetically ordered phase exhibits a critical phase transition with anomalous features such as critical exponents consistent with that of a first-order transition while maintaining continuity of physical quantities. Furthermore, at the critical point, we see a growing number of free energy local minima, supporting the picture of infinitely many discontinuous transitions forming criticality. Our numerical result unify de Gennes' theory of nucleation-type transitions and Dzyaloshinskii's calculation of chiral magnets by demonstrating that such anomalous criticality can be realized in microscopic models beyond mean-field.
We introduce statistical decision theory in stochastic thermodynamic systems. Statistical decision theory integrates several statistical tasks, such as estimation, hypothesis testing, and classification. In this study, we consider statistical decision theory in stochastic thermodynamics, which unifies sensing, classification, and uncertainty relations. In particular, we focus on the minimax risk, the smallest achievable expected loss under the worst-case parameter choice, even with an optimal estimator. We derive lower bounds on the minimax risk using Le Cam's two-point method. The derived relationships can be used across a broader range of thermodynamic scenarios. To illustrate these relations, we study classification and distribution estimation problems in a Markov process. We demonstrate that the minimax risk provides lower bounds given by the dynamical activity in these problems. Furthermore, by examining particular dynamics, we establish that the resulting bounds simplify to uncertainty relations involving dynamical activity or entropy production.
Models of colliding self-propelled rods are abstractions of many real systems such as gliding bacteria and biofilaments displaced by molecular motors. They have been studied mostly in two space dimensions and shown to display a wealth of different collective regimes with many parameters involved in the ensuing phase diagrams. Continuum-level theories faithful to this phenomenology are largely lacking since coarse-graining models of explicit rods is notoriously difficult. Here we study collections of particles in which the rods geometry is taken into account effectively, by implementing the nematic alignment of their velocities and two types of repulsion. These Vicsek-style models are simple enough to be starting points for systematic coarse-graining. We thus derive hydrodynamic equations that we show to have solutions qualitatively faithful to their particle-level origin.
Gain, loss, and nonlinearity shape the phase-space landscape of driven-dissipative systems, such as Kerr parametric oscillator (KPO) networks realized in superconducting, nanomechanical, and photonic arrays. Networks of KPOs exploit the stationary states of this landscape for Ising machines, bosonic error-correcting codes, and sensing. Yet how transients reach these states from a given initial condition remains largely unexplored. Here we show that KPO networks with nonlinear couplings can develop avalanches in their transient dynamics: rapid cascades of amplification between modes reminiscent of the relaxation avalanches of glassy systems. The mechanism is that parametric gain first activates the modes whose detuning it can overcome, and their nonlinear frequency shifts then pull dormant modes into the gain window one at a time. Dormant modes keep decaying while they wait for their turn, so every activation starts from a smaller seed and the intervals between activations grow geometrically. The network thus ages: the longer it has evolved, the slower it changes, and activations accumulate only logarithmically in time. In a 20-mode network, each dwell lasts about twice as long as the previous one. While aging delays the approach to the stationary state, additive noise sets a floor on the seeds and unsticks the network, shortening the longest dwells first while the activation sequence survives to noise several orders of magnitude stronger.
Motivated by gaps in our current knowledge of the Vicsek model, we present a major update of its phase diagram in two space dimensions, summarizing comprehensive numerical simulations at length- and timescales never explored before. Our main conclusions are that multiple phase separation scenarios are at play, not just the simple microphase separation put forward in [Phys. Rev. Lett. 114, 068101 (2015)], and that the coexistence phase is subdivided in two. We also map out the region of metastability of the homogeneous Toner-Tu liquid phase.
The Vicsek model continues to be central in statistical physics studies of active matter, with major recent surprises such as the partial metastability of the famous Toner-Tu flocking phase. Here we report on the diversity of band patterns observed in its coexistence phase --some of them never observed before--, but we also argue that most disappear in the infinite-size limit, leaving nevertheless a different picture than anticipated previously. We finally show that these rather regular, globally ordered band configurations can be destroyed by a localized perturbation: the coexistence phase is also at least partially metastable, and thus must be, asymptotically, a disordered phase made of a variety of chaotically interacting traveling bands.
Higher-order Ehrenfest phase transitions are exceptionally rare and poorly understood. We present an exactly solvable one-dimensional model exhibiting a sequence of Ehrenfest transitions of arbitrary $n\geq2$ order. These transitions occur without symmetry breaking and even the second-order transition lacks conventional critical behavior, although the density derivative of the correlation length diverges. Our results reveal a distinct class of noncritical phase transitions and help clarify the scope of the Ehrenfest classification.
Fluctuating interfaces are typically characterized by their global, spectral, and local roughness exponents. Faceted anomalous scaling, characterized by $α_s>1$, $α_s\neqα$, and $α_{\mathrm{loc}}=1$, owes its name to the tendency of such interfaces to develop long, nearly linear segments called facets. Its statistical origin, and whether facets are required, have remained unclear. We show that this exponent structure is determined by the spatial statistics of the slope field. Slope increments split into typical fluctuations and rare, larger jumps. This anomalous scaling arises either when the typical slope increment displays genuine system-size-dependent anomalous scaling or when the rare jumps dominate the second moment. The two contributions can coexist, and their competition leads to bifractal scaling of the slopes. We establish this framework analytically in the generalized anharmonic Larkin model and the columnar Kardar--Parisi--Zhang equation, and we derive exactly the global, spectral, and local roughness exponents in both cases. The Larkin model realizes both mechanisms and produces anomalous interfaces with or without facets, whereas columnar KPZ belongs to the jump-dominated case and yields $α_s=3/2$.
Mean First Passage Time (MFPT) is a crucial quantity that characterizes stochastic processes in various stochastic dynamical systems across different disciplines. The inverse of the MFPT, known as the escape rate, is helpful in the context of Brownian motion and chemical processes. The MFPT dynamics of a passive particle in the presence of active noises presents a fascinating avenue. In this work, we study the dynamics of two different systems: first, a single Brownian particle and second, a dimer (two Brownian particles connected by a spring). In both cases, the systems are kept in contact with a thermal bath and an active bath, and are trapped in a bistable potential. The active bath produces kicks that follow the Poisson distribution. The MFPT of the particle and the dimer is investigated in various scenarios. Analytical results are presented for the small kick variance approximation only. In this approximation, we formulate the MFPT in two scaling regimes: a weak-noise regime characterized by classical Arrhenius scaling for escape rate and a strong-noise regime in which the MFPT scales as a power law in the effective noise strength. For the general case, the problem is solved numerically. Next, under the same approximation, the MFPT of a dimer in the small-bond-strength regime is studied using a sequential-hopping mechanism, and in the large-bond-strength regime, the MFPT is formulated using the adiabatic elimination approach. Motivated by recent studies on speeding up the escape mechanism, the MFPT of the particle and the dimer in the presence of an intermediate barrier is also examined. We deduced from the numerical results that the MFPT can be enhanced and further optimized by selecting a suitable intermediate potential profile.
We investigate the non-equilibrium stationary states of a one-dimensional active gas composed of run-and-tumble particles (RTPs) undergoing motility-induced phase separation (MIPS) and subjected to local stochastic resetting at the origin. Unlike global resetting, local resetting maintains a mass-injection flux at the origin while preserving large-scale spatial self-organization. We demonstrate that above a critical resetting rate $r_c$, the active system undergoes a transition toward a stationary two-phase state featuring a dense liquid droplet centered at the origin coexisting with a surrounding gas phase. Using kinetic evolution equations, we derive the exact steady-state density profile expressed in terms of the Lambert function. The size of the droplet has a non-monotonic dependence on the resetting rate $r$, and attains its maximum. Furthermore, we analyze the dynamical stability of the supersaturated gas phase and show that a nucleated droplet relaxes exponentially toward its steady-state size with a relaxation timescale of $1/r$. We extend our framework to model a finite resetting velocity, proving that the stationary droplet size is invariant with respect to the resetting velocity.
Classical density-functional theory is a powerful theoretical framework for studying the equilibrium density distribution of complex many-body systems. The core of the framework is to bridge microscopic and macroscopic descriptions of fluids through a variational principle that minimises the system's free energy. For dynamical settings, dynamic DFT (DDFT) extends the equilibrium framework to capture the time evolution of an order parameter, typically density. However, DDFT is often insufficient for mesoscopic systems far from equilibrium where fluctuations become significant, necessitating a stochastic extension: fluctuating DDFT. Here we derive a generalised FDDFT from first principles starting from the microscopic Hamiltonian dynamics using coarse-graining techniques, including Zwanzig's projection operator method. The resulting FDDFT applies to arbitrary observable fields, where the dynamics are driven by a thermodynamic force expressed in terms of an entropy functional. The dynamics consist of reversible, dissipative, and stochastic contributions, each directly linked to the underlying microscopic dynamics via the Liouville operator, Green-Kubo relations, and the fluctuation-dissipation theorem. To highlight the difference in the level of description between DDFT and FDDFT, we demonstrate how further coarse-graining via averaging leads to a generalised DDFT, valid at the macroscale. To showcase the effectiveness and utility of the generalised FDDFT, we apply it to two prototypical fluid systems: colloidal and simple viscous fluids. In the former case, we utilise the new formalism to re-derive the recently found colloidal FDDFT. In the latter case, we rigorously recover the non-isothermal extension of the Landau-Lifshitz fluctuating Navier-Stokes equations, providing fundamental derivation of this well-known phenomenological model, enabling to clarify its underlying principles.
Using the public database of Google Scholar article citations of scientists and universities we apply the algorithm for constructing Lorenz curves to the resulting citation distributions. Using Lorenz and Pareto curves, broadly used for economic analyses of wealth inequality across countries, we show that this approach gives a useful insight of citation distributions of world-leading scientists characterized by the Gini coefficient $G$ obtained from the Lorenz curves. For articles of 15 Nobel laureates (including A. Einstein and R. Feynman) we find the average $G=0.81$ with 50% of articles capturing only 2.7% of all citations and the top 10% (1%) of articles capturing 71% (30%) of all citations. Furthmore, we determine for world-leading universities the Gini coefficient via the citations of affiliated scientists. We show that about 65 percent of all citations belonging to a university (e.g. Princeton, Cambridge) are generated only by 10 percent of their university scientists while 50% of scientists produce only 3% of all university citations. The citation distributions are well described by the concept of Rayleigh-Jeans thermalization and condensation with similarities traced to the world wealth inequality.
Motivated by the conflicting findings in prior studies regarding the nature of the phase transition for the three-dimensional antiferromagnetic $J_1$-$J_2$-$J_3$ Ising model on the cubic lattice with $J_1 = J_2 = J_3 = 1.0$, in this study we conduct a detailed investigation of this model using both Monte Carlo simulations and neural networks (NN). Here $J_1$ and $J_2$ are the intralayer nearest-neighbor and next-nearest-neighbor couplings, respectively, and $J_3$ is the nearest-neighbor interlayer interaction. By investigating the data obtained from much larger lattices than those used in previous studies, we find that the phase transition for $J_1 = J_2 = J_3 = 1.0$ is second order without any ambiguity. Hence we bring a resolution to this debate. Apart from $J_1 = J_2 = J_3 = 1.0$, we study the critical phenomenon associated with $J_2/J_1 = 4/7$ and $J_3 = J_1 = 1.0$ as well. In addition, by employing a dedicated training and testing strategy, our NN computation not only is much more efficient than those utilizing the standard procedures commonly used in the literature, but also requires only a per-mille fraction of the storage capacity needed by a conventional NN calculation when a large linear system size $L$ is considered. The method employed here provides an extremely resource-efficient route for applying NN techniques to study the critical behaviors of higher-dimensional systems.
The \emph{kappa} distribution is widely used to model heavy-tailed, non-equilibrium velocity distributions in space and laboratory plasmas. Its correct use in particle-in-cell (PIC) simulations requires more than sampling the single-particle marginal: the $N$-particle \emph{kappa} distribution predicts non-trivial correlations between the kinetic energies of particles. Building on the superstatistical structure, we implement a Gibbs sampler whose stationary distribution is the joint distribution of the velocities and $β$, so that the sampled velocities follow the $N$-particle \emph{kappa} distribution and reproduce these correlations by construction.
The frustrated square-lattice J1-J2 Ising model exhibits rich magnetic phases, yet its phase transitions remain difficult to characterize due to critical slowing down in conventional Monte Carlo and truncation errors inherent to tensor renormalization group methods.We develop a tensor-network Monte Carlo (TNMC) algorithm for this frustrated model, which supports both open and periodic boundary conditions. For periodic lattices, we use a sweeping block-update strategy that updates spin blocks while fixing surrounding spins.To boost efficiency, we couple our sampler with histogram reweighting, and test each technique standalone to validate these independent optimizations.We perform extensive finite-size scaling analyses at g=0.55 and g=1.0, extracting apparent critical temperatures and effective thermal exponents from Binder-cumulant data.Our method offers an efficient, nearly bias-free path for large-scale simulations of J1-J2-type frustrated magnets under periodic boundaries and generalizes readily to other frustrated spin lattices.
We study the dynamical scaling properties of the two-dimensional fully frustrated Ising model by simulated annealing. The model has a highly degenerate, critical ground state with well known equilibrium properties. To study its out-of-equilibrium behavior, we develop an annealing protocol that is convenient for extending the Kibble-Zurek mechanism to the case of an exponentially divergent correlation length in the limit of the temperature $T \to T_c = 0$. We find excellent agreement with the scaling forms derived, including a dependence on the initial temperature that is not present in the conventional Kibble-Zurek mechanism. The dynamic exponent $z$ is consistent among all quantities studied; fits of its dependence on the initial temperature give $z = 1.99(2)$ (susceptibility), $z = 2.072(8)$ (excess energy), and $z = 2.09(3)$ (Edwards-Anderson order parameter). These values are close to $z=2$, which is consistent with previous studies of relaxation after instantaneous quenches and of time dependent correlations in equilibrium, including the estimate $z = 2.004(9)$ from equilibrium autocorrelation functions.
We investigate the distributions of off-diagonal scattering matrix elements and cross-sections in the Ericson regime for systems that are time-reversal and rotationally invariant. In previous works, we developed techniques to analyze the Ericson transition by means of asymptotic expansions in the cases of unitarily and symplectically invariant Hamiltonians. Here, we prove that Gaussian and exponential behavior is found as the leading-order term of an asymptotic expansion of those distributions also in the orthogonal case, which is often the relevant one for analyzing experiments. Moreover, we derive the subleading-order term and compare with numerical simulations. Finally, we discuss similarities and differences of the Ericson transition for all three Wigner-Dyson classes.
We consider a Brownian particle with position $x(t)$ whose diffusivity changes at Poissonian switching times, with each new value drawn independently from a bounded distribution $W(D)$. Assuming $W(D)\sim C(D_{\max}-D)^ν$ near its upper edge, we study the large-deviation regime $x(t)\simeq yt$ with $y = O(1)$ at large times, far beyond the typical diffusive scale. In a recent work Guéneau et al., Phys. Rev. Lett. 135, 067102 (2025), the rate function of the displacement was obtained through a renewal approach and it was shown to display several singularities. In this paper, we carry out a detailed study of the corresponding dynamical phase transitions and trajectory condensation. For moderate $y$, the imposed displacement is shared among $O(t)$ displacement increments accumulated during successive intervals between switches. Beyond a first threshold $y_-$, one interval produces an increment of order $t$ carrying a finite fraction of the total displacement, signaling dynamical condensation. For $ν>1$, a second threshold $y_+>y_-$ is reached when this macroscopic increment carries asymptotically the entire imposed displacement, whereas for $0<ν\leq1$ the two thresholds merge, $y_-=y_+$. In turn, the moment-generating function of $x(t)$ displays a complementary transition from collective contributions of many diffusivities to selection of states near $D_{\max}$. We show that this transition is governed by the same edge mechanism as the Baik-Ben Arous-Péché (BBP) transition, indicating that diffusivity selection and trajectory condensation are two manifestations of the same phenomenon. We validate our predictions for the dynamical condensation against numerical results from rare event sampling.
A method for integrating the Langevin equations of motion for run-and-tumble particles is introduced to compute the joint distribution of their position and active velocity, with and without a harmonic external force field. The key step is the construction of the moments of the distribution and the correlation functions, conditional on the particles being in a specific state at the final time. For one-dimensional run-and-tumble particles in a harmonic external potential, the exact steady-state distributions for the tumble state, as well as for the run state conditional on a specified velocity direction, are obtained. This enables us to demonstrate their equivalence to the exact solutions of the corresponding Fokker-Planck equations.
The long distance behavior of two-point functions of critical lattice statistical mechanics models depend on detailed information of the underlying quantum field theory. Yet these universal properties can be difficult to extract because generic lattice operators couple to the infinite hierarchy of continuum scaling operators with unknown coefficients. Recently in Physical Review E 114, 024146 (2026), one approach to disentangle some of this information was proposed using the eigen-decomposition of the operator covariance matrix constructed from cross correlations of a finite set of $N_\ell$ lattice operators. In this work we analyze this method in detail. We provide an explicit proof that under the assumption of scale invariance the leading decay of the eigenvalues are controlled by the smallest $N_\ell$ scaling dimensions, even after taking irrelevant perturbations into account. Based on this understanding, we propose a new method by which the eigenvectors of the covariance matrix can be used to construct the amplitudes of the continuum two-point correlation functions. This allows for a complete determination of the structure of the continuum two-point functions from the lattice covariance matrix. We demonstrate our ideas both on synthetic data sets (drawn from the scaling hypothesis) and numerically exact data for the two dimensional Ising model.
Frustration can enhance fluctuations and create new and interesting phase behavior. When combined with the ability to host topological defects, frustration can increase the complexity of a system's phase diagram through topological and frustrated topological transitions. To understand the interplay of frustration and topology on phase behavior, we introduce and study solute-frustrated $q$-clock (SF$q$C) models, in which solute sites are added to the $q$-clock model, a lattice model where each spin can take one of $q$ orientations. The $q$-clock model can exhibit spontaneous symmetry breaking and topological transitions, and the interaction between the spins and a solute frustrates the system, in essence by creating a fixed lower bound on the concentration of positive topological defects or vortices. We perform Monte Carlo simulations and develop theoretical descriptions of SF$q$C models with $q=4$ and $q=6$ to understand their phase behavior. We find that these systems exhibit rich phase diagrams that depend on the choice of parameters. In particular, SF$q$C models have a wide range of interesting low-temperature phases, including the presence of strongly frustrated domains and the absence of long-range order. We also find topological effects even at low temperatures; for instance, we find that topological defects can drive crystallization. We anticipate these results to guide the development of materials that take advantage of frustration and topology to tune phase behavior and fluctuations.
We study spontaneous particle-hole symmetry breaking in a disordered two-site quantum breakdown model. The model is solvable in the limit of a large number of fermionic modes on each site. For a small ratio of the numbers of modes on the two sites, we derive the free energy analytically and find a finite-temperature symmetry-breaking transition. This transition changes from continuous to first order as the number of fermions produced in each conversion increases. Numerical results show that both types of transition persist beyond this regime, with an order-parameter exponent close to $1/2$ for continuous transitions.