Papers

We study Gibbs states projected onto the symmetric sector of contractible one-form symmetries. In a three-dimensional cluster-model interpolation, this projection stabilizes symmetry-protected topological and projected-paramagnetic phases, both sharply distinct from thermal disorder. A flux-twisted membrane invariant distinguishes these three phases by the values $-1,+1,0$, respectively. These values are exact on the endpoint, self-dual, zero-temperature, and infinite-temperature lines; elsewhere their quantization requires positive spatial-sheet, winding-line, and interface tensions. Quantum Monte Carlo supports the quantization through tension diagnostics and direct finite-size estimates.

The distribution of node depths in a network is crucial for analyzing network structure. Two measures, rank and layer, quantify how deep inside a network a node is. The rank is the node's distance to the closest leaf, a node of degree 1. The layer is the number of times all leaves must be stripped off for the node to become a leaf. We derive exact recursive expressions for the rank and layer distributions in random recursive trees. We show that the rank and layer distributions decay factorially and geometrically, respectively, and prove the self-averaging of the numbers of nodes with fixed rank or layer. Unlike previous studies, our approach does not depend on labels or a root node, providing a more versatile framework for analyzing rank and layer distributions in complex networks.

Data-driven discovery of governing equations from spatiotemporal data offers a promising route to obtaining coarse-grained descriptions of complex dynamical systems. Here, we investigate the performance of PDE-SINDy for discovering phase-ordering dynamics using the Allen–Cahn equation as a benchmark and the Ising model with Glauber spin-flip dynamics as a microscopic system. We systematically analyze the effects of data availability, size of the candidate library, and noise on the efficiency of the equation discovery. We find that stability-selection PDE-SINDy can robustly identify the relevant terms in the governing dynamics even under limited or noisy data, while the recovered coefficient values are substantially more sensitive to these factors. We further show that enlarging the candidate library can strongly affect both term identification and coefficient recovery. Incorporating library bagging with stability selection reduces this sensitivity and improves the efficiency of equation discovery. For the Glauber spin flip Ising model dynamics, the resulting coarse-grained equation reproduces the characteristic phase-separation and coarsening dynamics of the underlying microscopic system. Overall, our results demonstrate the potential of PDE-SINDy for phase-ordering systems while highlighting the importance of carefully assessing the factors that influence the efficiency of equation discovery.

We analyze the Gibbs variational principles associated with the probability distributions of (i) an isolated system and (ii) a system at constant temperature. We give an example of using the Gibbs inequality to obtain the free energy and analyze the phase diagram of an Ising model with competing interactions. We also review the Boltzmann irreversible theorem, and show how it is connected to the Gibbs variational principles. This connection is established by using the Kolmogorov equation for the evolution of the probability distribution, which predicts a monotonic increase of entropy for an isolated system, and a decrease of the free energy for a system in contact with a thermal reservoir.

Phase reduction provides a principled relationship between high-dimensional oscillator dynamics and low-dimensional phase models such as the Kuramoto model. However, connecting analytical results derived from the Kuramoto framework to experimental data remains an open challenge, as many physical platforms violate the weak-coupling assumptions underlying standard phase reduction. Here, we present an experimental platform based on electronic quadrature oscillators whose phase dynamics are, under certain assumptions, equivalent to the Kuramoto model. Crucially, this equivalence is achieved for both weak and strong coupling through a non-diffusive coupling circuit, and the design approach supports arbitrary (weighted and directed) network topologies while preserving scalability and oscillation regularity. Across different coupling schemes, we demonstrate that the platform reliably reproduces phase transitions predicted by the Kuramoto model for both global and cluster synchronization. This cost-effective platform enables systematic experimental validation of analytical and numerical results for synchronization in complex networks.

The subsystem-size dependence of ground-state entanglement entropy and its crossover to thermal entropy as a function of temperature are well understood for one-dimensional (1D) gapless systems described by conformal field theory (CFT), and for free fermions with a Fermi surface in any dimension. However, little is known about the entanglement entropy for gapless fermionic systems without quasi-particles, such as a strange metal. Here we study the entanglement entropy of fermions in a solvable large-$N$ 1D lattice model akin to the Yukawa-Sachdev-Ye-Kitaev (Yukawa SYK) model. In this model, two Fermi points are coupled to scalar bosons via spatially random Yukawa interactions, providing a solvable model of a strange metal when the bosons become critical at a quantum critical point. We exactly compute the second Rényi entropy of fermions in a spatial subregion in this model. Our results unravel crucial role of intra-subregion entanglement between fermionic and bosonic degrees of freedom along with the inter-subregion entanglement in understanding the ground states of such strongly coupled fermion-boson systems. We show that the crossover from thermal entropy to entanglement entropy, is captured by a single scaling ansatz, that collapses the second Rényi entropy of fermions for different subregion sizes and temperatures into a single universal curve. We further show that the universal scaling curve for the critical strange metal is well described by the standard CFT formula, albeit with an effective central charge substantially larger than the non-interacting value.

We study a $\mathscr{PT}$-symmetric non-Hermitian multichannel Kondo model consisting of a pair of spin-$\frac12$ impurities coupled to $n$ conduction-electron channels through complex-conjugate Kondo couplings. The impurity renormalization-group (RG) flow is characterized by the Kondo scale $T_K$ and a dimensionless non-Hermiticity parameter $\alpha$. As $\alpha$ increases, the exact Bethe Ansatz solution exhibits four impurity phases: overscreened Kondo, zero mode, Yu–Shiba–Rusinov (YSR), and local moment. The Kondo, zero-mode, and local-moment phases are $\mathscr{PT}$-unbroken, whereas the YSR phase spontaneously breaks $\mathscr{PT}$ symmetry. Using a generalized thermodynamic Bethe Ansatz, we determine the impurity free energy and Affleck–Ludwig $g$-function throughout the $\mathscr{PT}$-unbroken phases. In the Kondo phase, the defect RG flow connects the ultraviolet and infrared conformal fixed points, with the impurity entropy flowing from $2\ln2$ to $2\ln\left[2\cos\left(\frac{\pi}{n+2}\right)\right]$, in agreement with defect conformal field theory. In the zero-mode phase, zero-energy impurity strings reorganize the spectrum into multiple excitation towers, while in the local-moment phase, the RG flow becomes cyclic, returning to the unscreened local-moment fixed point. We conjecture that RG irreversibility, and hence a generalized Affleck–Ludwig $g$-theorem, survives throughout the Kondo phase. Our exact solution nevertheless shows that a real spectrum and defect entropies consistent with defect CFT do not guarantee RG irreversibility: the impurity entropy is non-monotonic in both the zero-mode and local-moment phases.

The minimum change principle provides an information-theoretic characterization of the Bayes reversal channel in classical probability theory and has recently been proposed as a framework for extending Bayes' rule to quantum information theory. Using quantum relative entropy, we investigate a minimum change principle for the setting of quantum statistical inference. Specifically, we consider a forward process based on a classical-to-quantum preparation channel and a reverse process based on a quantum-to-classical measurement channel. We establish a closed-form characterization of measurements that are optimal for this principle, and this optimal measurement can be found via a dual formulation involving a single unconstrained Hermitian variable. This perspective allows us to recover some notable measurements within the same framework, including pretty good measurements and Fermi-Dirac thermal measurements, and we use it to discover a novel family that we call softmin thermal measurements. We further show that softmin thermal measurements arise as optimal solutions to entropy-regularized semidefinite optimization problems, demonstrating that they play a role for measurements analogous to that of thermal states in statistical mechanics. Finally, we prove an additivity property for the relative-entropy minimum change principle and investigate the performance of Fermi-Dirac thermal measurements for quantum hypothesis testing.

Recently (Physica A 660 130370, 2025), emergence of superstatistical behavior in driven classical systems has been shown for systems with constant microcanonical heat capacity. As an application of this result, in this work we show that a system of $N$ particles with inverse gamma distribution of total kinetic energies must have kappa-distributed single-particle velocities in the limit $N \rightarrow \infty$. Our results provide insight into the nature of kappa distributions outside the theory of nonextensive statistical mechanics, while also bringing forward a practical method for the generation of kappa velocities via Monte Carlo Metropolis simulation in the inverse gamma ensemble.

A finite density of mutually avoiding directed polymers in a quenched random medium is a minimal model of glassy line matter. The dilute theory, solved by replica Bethe ansatz, predicts an interaction free energy proportional to $\rho^2$ and disorder cumulants with distinct power-law dependences on the density $\rho$, but direct numerical tests have been hindered by the combinatorially large many-polymer transfer matrix. We recast the problem as filling logarithmic eigenvalues of a single-polymer transfer-matrix product, obtaining the quenched free energy, its cumulants, and a disorder-induced linear spectral edge consistent with the replica prediction.

Entropy production is a key property in stochastic thermodynamics. For partially observed and coarse-grained systems, its inference is challenging and typically rests on proven lower bounds. We derive two versions of such bounds based on the asymmetry of experimentally accessible two-time correlation functions of coarse-grained state observables. For non-equilibrium steady states, the bound is valid for arbitrary correlation lag. For time-dependent processes, it requires the limit of vanishing lag. These bounds hold true for any system that follows either a Markovian dynamics or a coupled set of overdamped Langevin equations on some underlying, unobservable level of description. We illustrate the bounds for both types of dynamics and discuss their optimization and potential tightness.

This paper investigates the numerical initialization of the one-dimensional Maxwell–Cattaneo–Vernotte (MCV) heat conduction model, addressing the critical mathematical challenge of defining the initial time derivative for non-equilibrium states. In modern high-frequency thermal engineering applications, traditional time-integration schemes and commercial finite-element solvers frequently introduce severe numerical artifacts when handling non-Fourier models. This study systematically evaluates three initialization strategies: a zero derivative, a spatially uniform non-zero derivative, and an exact space-dependent derivative. Using an explicit staggered finite-difference scheme, the transient responses to an exponentially distributed initial temperature field under adiabatic boundary conditions are compared against an analytical solution using the Galerkin method. The results demonstrate that assuming a zero or spatially uniform initial derivative introduces significant unphysical oscillatory deviations, leading to heat-flux prediction errors as the relaxation time increases. Conversely, mapping the exact space-dependent derivative onto the staggered grid preserves the local thermodynamic structure of the initial state, yielding robust transient responses that match the analytical benchmark without meaningful computational overhead. These findings establish a thermodynamically consistent initialization technique that can be extended to non-local models as well.

We introduce an extension of the Fibonacci number system, we call the fractional Fibonacci number system, which interpolates between the Fibonacci number system (used for natural numbers) and the irrational base-$\phi$ number system, which can be used to represent real numbers. The number system finds its use in interpreting the localization phase diagram of the many-body non-interacting Aubry-André model. For finite system sizes a localization criterion can be obtained if the particle density is written in the fractional Fibonacci number system. In the thermodynamic limit, the criterion remains, but in this case the particle density is expressed in base-$\phi$. The nature of the thermodynamic limit is also discussed.

The intermediate scattering function and its value at zero time, which is the static structure factor, are both characteristic functions. An exact equation of motion for the ISF carrying a local rate and a memory function is then obtained. To interpret this equation, the ISF is expressed in terms of relaxation modes. The local rate is fixed by a sum rule which reversibility turns into an equilibrium average. The de Gennes narrowing is establshed as a theorem. The memory function being non-negative, the closed form that results is a rigorous lower bound on the ISF at every momentum transfer and every time, and its accuracy is itself predicted. In the diffusive regime, the closed form keeps the single-exponential shape to which spin-echo data are routinely fitted, with both coefficients renormalized by the interaction: the amplitude becomes the SSF of the layer and the dephasing rate the dressed one-adsorbate rate divided by it, explicit in momentum transfer, time, coverage and temperature. When introducing the lateral interaction, a single correlation parameter fixes the equilibrium structure and the renormalized hop rate: exactly along a one-dimensional channel, and on the surface as the nearest-neighbour correlation of a pair approximation completed by a collective (random-phase) structure factor. The result contains no adjustable parameter and is exact at vanishing interaction for every coverage and momentum; the Darken relation comes with it, the correlation between successive moves separated off as a computable factor. Applications to Na/Cu(111) and H/Pt(111) are analysed and discussed.

By investigating the origin of the evolutionary discrepancy between the H-function and entropy, we elucidate that the H-function fails to serve as a valid arrow-of-time criterion in the rigorously derived Boltzmann equation for the hard-sphere gas model under the Boltzmann-Grad limit.

We develop a continuum theory for proliferating active matter starting from a microscopic stochastic model of self-propelled particles undergoing birth, death, and nonlocal competition. Beginning from the master equation, we derive mean-field evolution equations for the particle density and polarization fields and close the resulting hierarchy through a von Mises ansatz, providing a coupled hydrodynamic description applicable to a broad class of proliferating active systems. As an application, we study the recently introduced Active Brownian Bug model, in which the form of flocking emerges despite the absence of explicit alignment interactions. Linear stability analysis predicts both Turing and Hopf instabilities, whose analytical thresholds agree with numerical simulations. The continuum model reproduces the principal dynamical regimes of the underlying particle system, including homogeneous states, stationary periodic clusters, and coherently propagating flocking states. These results establish a general continuum framework for proliferating active matter and provide a physical interpretation of collective motion driven by the interplay between activity and population dynamics.

Multi-phonon resonance conditions underpin kinetic theories of phonon transport and lattice thermalization. We show that exact resonance matching, nonzero interaction coefficients, and network connectivity do not guarantee persistent energy diffusion. Symmetry-enforced balance relations drive exact-resonant collision currents to nonthermal zero-flux states, producing kinetic arrest from individual resonant sets to connected networks. Complete energy spreading is sustained by quasi-resonances. The thermodynamic and weak-nonlinearity limits do not commute: the leading kinetic behavior is recovered in the former, whereas at fixed finite size the thermalization time diverges through higher-order crossovers as the nonlinearity vanishes. Exact-resonance existence and connectivity are therefore kinematic, not sufficient dynamical, criteria for phonon energy diffusion.

In many-body physics, an explicit wavefunction often provides the most thorough understanding, yet for monitored random Clifford circuit it has remained missing. Here we develop a framework that grants direct access to the typical output states. As every stabilizer state is local-Clifford-equivalent to a graph state, the graph adjacency matrix, a classical bit matrix, provides a complete description of the quantum state. In the large-$N$ limit, we show that these graphs converge to the Erdős–Rényi random graph $G(N,1/2)$, which allows us to resolve the open problem of Greenberger–Horne–Zeilinger (GHZ) entanglement generated by deep random Clifford circuits. We obtain analytically the mean GHZ content $\langle g_3\rangle=1.204$ for even $N$ and $1.325$ for odd $N$. For monitored Clifford circuits with a 1D layout, we uncover an emergent Erdős–Rényi subgraph $G(N_{\mathrm{sub}},1/2)$ in the output states of the volume-law phase, where $N_{\mathrm{sub}}/N\approx \sqrt{1-p/p_c}$ with $p$ the measurement rate and $p_c$ the critical point of the measurement-induced phase transition (MIPT). The output state is thus equivalent to the output of an unmonitored random Clifford circuit on $N_{\mathrm{sub}}$ qubits, weakly perturbed by the remaining $N-N_{\mathrm{sub}}$ qubits carrying little entanglement. This result directly accounts for the quantum error-correcting capability of the volume-law phase, and implies the same GHZ statistics for the whole volume-law phase. We further identify a clustering effect for qubits in the dense subgraph, which we reproduce with an infection-recovery toy model that exhibits a measurement-induced absorbing-state phase transition. Finally, a mean-field argument on the graph locates the MIPT critical point at $p_c = 0.1608$, in excellent agreement with the numerical value $p_c\approx 0.16$.

We study simple inhomogeneous deformations of known two-dimensional classical lattice models described by conformal field theory (CFT) at large distances. The deformations are chosen to vary slowly at lattice scales, while preserving critical behavior. Globally, we find that such systems are described by a CFT in curved space, and identify the underlying space metric. Our two examples are the Ising model and the six vertex model with domain wall boundary conditions. In the latter boundaries are also inhomogeneous, which complicates the analysis. We nevertheless solve the free case using hydrodynamics. In the presence of interactions we determine the arctic curve which separates a critical fluctuating region from an ordered (frozen) phase.

The eigenstate thermalization hypothesis (ETH) provides the prevailing framework for understanding quantum thermalization and ergodicity in isolated many-body systems. Yet, no general theory describes the continuous onset of ETH violation between the conventional ETH and its complete breakdown. Here, we introduce a two-parameter ansatz for the ETH violation that unifies and distinguishes two mechanisms: fading ergodicity and trapped ergodicity. While fading ergodicity captures the established route to ergodicity breaking, trapped ergodicity describes a distinct scenario in which ETH is violated in finite systems but restored in the thermodynamic limit. We test this framework in the spin-1/2 $J_1$-$J_2$ chains with on-site disorder and linear potential. In both cases, we find that the observed ETH violation is consistent with trapped ergodicity.

We study the six-vertex model on the Sierpinski gasket, a four-coordinated hierarchical fractal with Hausdorff dimension $d_f=\log_23$. Given the importance of dimensionality for the long-wavelength behavior of such models, we specifically consider correlations and confinement as a function of the vertex weights, with the equal-weight point corresponding to the ice model. We calculate the partition function and correlators recursively to obtain a rich phase diagram hosting many different regimes. While the ice model shows entropic charge confinement, a particular four-vertex limit exhibits fractal deconfinement, with the string joining the deconfined charges itself a statistical fractal with fractal dimension $d_l=\log_2(5/2)\approx 1.3$. Finally, we propose a setup as an artificial spin ice to enable experimental study of the rich phenomenology of Sierpinski ice and its generalisations.

We study post-measurement ensembles of ground states of tricritical and critical 1D quantum Ising Hamiltonians subjected, respectively, to weak energy and spin measurements without post-selection. These measurements act as relevant perturbations about the unmeasured critical ground states. Using finite-size renormalization group (RG) crossover analyses, we characterize their universal properties through the entanglement effective central charge, effective Affleck-Ludwig boundary entropy, and signatures of multifractality from moments of measurement-averaged correlation functions. In both cases, we find evidence for "measurement-dominated" or "measurement-altered" fixed points governed by the underlying Born-rule randomness. For critical Ising, we find a direct RG flow to a projective-measurement fixed point with area-law entanglement, whereas for the tricritical Ising model, we find evidence for a weak-measurement fixed point with logarithmic entanglement. These results clarify the RG-flow structure of weakly measured multicritical Ising ground states and show how intrinsic measurement-induced randomness can generate complex and rich universal long-distance scaling behavior in the post-measurement ensembles, accessible to controlled analytical RG and numerical finite-size RG crossover analyses.

Pseudomode approaches allow for an exact and unapproximated description of a quantum system interacting arbitrarily strongly with a bath. In general, a system coupled to pseudomodes will not thermalize to the system's Gibbs state: This is to be expected when the system-bath coupling is non-perturbative, but conflicts with common thermodynamic intuition when the system-bath coupling is asymptotically weak. We explore under which circumstances pseudomode models satisfy detailed balance (and subsequently thermalize to the system Gibbs state) and how specific choices of parameters can force "weak" detailed balance that is restricted to a limited frequency range. A combination of Hermitian and non-Hermitian pseudomodes that yields a flat effective-temperature profile is also considered. The results and criteria established here are relevant for the construction of pseudomode models in contexts where thermodynamic consistency is required.

We investigate the spread complexity of quantum quenches in generalized and long-range Aubry-Andre-Harper (AAH) models, encompassing regimes with and without mobility edges. In particular, in the generalized AAH models supporting energy-dependent mobility edges, we demonstrate that the long-time averaged spread complexity exhibits nonanalytic behavior when the post-quench quasiperiodic potential crosses the mobility edge associated with the energy of the initial eigenstate, thereby accurately identifying the mobility-edge transition. Such a behavior is supported by the spreading of local density of states. We further derive analytical expressions for the moments and the corresponding Lanczos coefficients for quenches between the limits of vanishing and strong quasiperiodic potentials. The Lanczos coefficients display qualitatively distinct behavior depending on the presence of mobility edges - they exhibit an initial plateau followed by a decay with the Krylov basis index, in contrast to the nearly constant behavior of the conventional AAH model without mobility edges. For LR hopping, the coefficients decay with the Krylov basis index for quenches from the localized to the extended phase, while they coincide with the AAH results for quenches in the opposite direction.

A common assumption in theoretical models of Bayesian inference is that the signal has i.i.d. components. To study the effect of correlations in the signal prior, we consider a minimal model: the planted spin glass on random regular graphs, where the signal is sampled from an Ising model with coupling $\kappa$. Depending on the phase of the prior, we find that adding structure in the signal can either help or hinder inference. In the paramagnetic regime, correlations in the signal lower the reconstruction threshold, so that weaker signal strength is sufficient for recovery. In the ferromagnetic regime, the prior alone already enables partial recovery, and we identify the threshold above which the observations provide additional information. When the prior itself is in a replica symmetry breaking (RSB) phase, we detect a static RSB transition in the posterior under Nishimori conditions. This provides an example where a non-separable, correlated prior leads to static RSB in a Bayes-optimal inference problem. We discuss the consequences of this glassy phase for algorithmic performance, in particular for Belief Propagation.

We introduce a quantum stochastic resetting protocol with uniform memory, in which each resetting event returns the system to a state visited at a time chosen uniformly from its entire history. The resulting dynamics is nonunitary, non-Markovian and a direct quantum generalization of the classical preferential relocation model. Working in the energy eigenbasis, we derive the exact evolution of every density-matrix element for an arbitrary time-independent Hamiltonian and show that the Hamiltonian enters the dynamics only through the corresponding Bohr frequencies. This leads to a natural distinction between two classes of quantum systems: gapped and gapless. In gapped systems (systems with a discrete energy spectrum), while the diagonal elements remain unchanged, the off-diagonal elements of the density matrix in the energy eigenbasis decay algebraically with a continuously varying exponent and with an amplitude that oscillates periodically in $\log t$. The system therefore approaches a stationary state that is independent of the resetting rate and retains a strong memory of the initial state. In gapless systems (systems with a continuous energy spectrum), arbitrarily small Bohr frequencies prevent stationarity. Instead, the position distribution spreads on the universal (ultra-slow) scale $\log(rt)/r$, independently of the initial state and of the details of the Hamiltonian. We illustrate these results with a two-level system, a harmonic oscillator, and a free quantum particle, and contrast them with their classical counterparts.

The Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equation provides a paradigmatic example of front propagation and asymptotic-state selection. Using a unified renormalization-group (RG) framework, we show that its traveling-wave solutions form a continuous family of RG fixed points and that all fronts with v >= 2 are locally RG-stable. The resulting RG stability criterion reproduces the classical linear-stability analysis of Canosa, according to which all fronts with v >= 2 are stable. Nevertheless, the asymptotically selected front is v=2, determined by the leading-edge structure of the initial condition rather than by RG stability. The FKPP equation thus provides an explicit example in which RG fixed-point stability does not determine asymptotic-state selection, demonstrating that stability and selection constitute distinct aspects of RG behavior.

Decentralized populations of Large Language Model (LLM) agents can spontaneously reach consensus on shared conventions, yet the microscopic mechanisms by which their internal stochasticity shapes macroscopic ordering remain unexplored. We study a minimal LLM Naming Game in which the listener's decision is a single-token LLM call at decoding temperature $T$, replacing the inventory check of the deterministic Naming Game. Each interaction decomposes into an in-inventory and an out-inventory channel with conditional rates $\pi(T)\!\equiv\!P(\text{YES}\mid w\in P_j)$ and $\phi(T)\!\equiv\!P(\text{YES}\mid w\notin P_j)$, whose balance controls an ordering-disordering drift. A mean-field theory of the two-rate dynamics yields an analytical ordering condition that generalizes the consensus threshold of the stochastic Naming Game to a critical line in the $(\pi,\phi)$ plane. Across three open-weight architectures, consensus is always reached, but through three distinct listener regimes: permissive (repaint-noise dominated), near-deterministic, and conservative (missed-collapse dominated). The effective finite-size exponent $\beta(T)$ in $t_{\rm conv}\!\sim\!N^{\beta}$ shifts with temperature, and the temperature-sensitivity $\alpha$ in $t_c\!\sim\!e^{\alpha T}$ ranges from ${\approx}\,0.67$ to ${\approx}\,0$ across architectures. Decoding temperature thus emerges as an architecture-dependent control parameter for decentralized LLM populations, quantitatively characterized by the statistical-physics toolkit.

We investigate the two-dimensional inhomogeneous Ising model (2DIM) on the kagom'e lattice by mapping it onto a particular non-symmetric eight-vertex model and constructing the corresponding $R$-matrix. Using a fermionic representation, we evaluate the partition function and derive explicit expressions for the main thermodynamic quantities. In the thermodynamic limit, we obtain an exact equation for the critical surface determining the phase transition of the model. We also calculate the free energy, specific heat, and spontaneous magnetization in the ferromagnetic case. Furthermore, we show that when one or two coupling constants vanish, the model reduces, respectively, to the square-lattice and one-dimensional Ising models. In both limits, our results reproduce the corresponding exact critical couplings and free energies.

We prove that the entanglement of formation of every two-mode Gaussian state coincides with its Gaussian restriction, solving a longstanding open problem in continuous variable quantum information theory. The key result is a sharp affine relation between entanglement and a generalized Einstein–Podolsky–Rosen observable, valid for arbitrary pure two-mode states, including non-Gaussian ones. We further prove that the entanglement of formation is additive on arbitrary tensor products of two-mode Gaussian states, and hence coincides with the entanglement cost, settling the latter for all two-mode Gaussian states. As a consequence, we obtain in closed form the entanglement cost of thermal-loss, thermal-amplifier and additive-noise bosonic channels. Our results, which also extend to bisymmetric multimode states, establish that no two-mode state is less entangled than the Gaussian state with the same covariance matrix, yielding measurable lower bounds on both the entanglement of formation and the entanglement cost of arbitrary non-Gaussian states.