Date: Fri, 19 Jul 2024 | Total: 15

Nonequilibrium systems, in particular living organisms, are maintained by irreversible transformations of energy that drive diverse functions. Quantifying their irreversibility, as measured by energy dissipation, is essential for understanding the underlying mechanisms. However, existing techniques usually overlook experimental limitations, either by assuming full information or by employing a coarse-graining method that requires knowledge of the structure behind hidden degrees of freedom. Here, we study the inference of dissipation from finite-resolution measurements by employing a recently developed model-free estimator that considers both the sequence of coarse-grained transitions and the waiting time distributions: $\sigma_2=\sigma_2^\ell + \sigma_2^t$. The dominant term $\sigma_2^\ell$ originates from the sequence of observed transitions; we find that it scales with resolution following a power law. Comparing the scaling exponent with a previous estimator highlights the importance of accounting for flux correlations at lower resolutions. $\sigma_2^t$ comes from asymmetries in waiting time distributions, with its peak revealing characteristic scales of the underlying dissipative process. Alternatively, the characteristic scale can be detected in a crossover of the scaling of $\sigma_2^\ell$. This provides a novel perspective for extracting otherwise hidden characteristic dissipative scales directly from dissipation measurements. We illustrate these results in biochemical models as well as complex networks. Overall, this study highlights the significance of resolution considerations in nonequilibrium systems, providing insights into the interplay between experimental resolution, entropy production, and underlying complexity.

We report a new intermediate regime in the quench time, $\tau_{q}$, separating the usual validity of the Kibble-Zurek mechanism (KZM) and its breakdown for rapid quenches in open systems under finite quench protocols. It manifests in the power-law scaling of the transition time with $\tau_{q}$ as the system appears to enter the adiabatic regime, even though the ramp is already terminated and the final quench value is held constant. This intermediate regime, which we dub as the delayed KZM, emerges due to the dissipation preventing the system from freezing in the impulse regime. This results in a large delay between the actual time the system undergoes a phase transition and the time inferred from a threshold-based criterion for the order parameter, as done in most experiments. We demonstrate using the open Dicke model and its one-dimensional lattice version that this phenomenon is a generic feature of open systems that can be mapped onto an effective coupled oscillator model. We also show that the phenomenon becomes more prominent near criticality, and its effects on the transition time measurement can be further exacerbated by large threshold values for an order parameter. Due to this, we propose an alternative method for threshold-based criterion which uses the spatio-temporal information, such as the system's defect number, for identifying the transition time.

Odd-diffusive systems are characterized by transverse responses and exhibit unconventional behaviors in interacting systems. To address the dynamical interparticle rearrangements in a minimal system, we here exactly solve the problem of two hard disk-like interacting odd-diffusing particles. We calculate the probability density function (PDF) of the interacting particles in the Fourier-Laplace domain and find that oddness rotates all modes except the zeroth, resembling a ``mutual rolling'' of interacting odd particles. We show that only the first Fourier mode of the PDF, the polarization, enters the calculation of the force autocorrelation function (FACF) for generic systems with central-force interactions. An analysis of the polarization as a function of time reveals that the relative rotation angle between interacting particles overshoots before relaxation, thereby rationalizing the recently observed oscillating FACF in odd-diffusive systems.

We analyze the dynamics of the Sisyphus random walk model, a discrete Markov chain in which the walkers may randomly return to their initial position $x_0$. In particular, we present a remarkably compact derivation of the time-dependent survival probability function $S(t;x_0)$ which characterizes the random walkers in the presence of an absorbing trap at the origin. The survival probabilities are expressed in a compact mathematical form in terms of the $x_0$-generalized Fibonacci-like numbers $G^{(x_0)}_t$. Interestingly, it is proved that, as opposed to the standard random walk model in which the survival probabilities depend linearly on the initial distance $x_0$ of the walkers from the trap and decay asymptotically as an inverse power of the time, in the Sisyphus random walk model the asymptotic survival probabilities decay exponentially in time and are characterized by a non-trivial (non-linear) dependence on the initial gap $x_0$ from the absorbing trap. We use the analytically derived results in order to analyze the underlying dynamics of the `survival-game', a highly risky investment strategy in which non-absorbed agents receive a reward $P(t)$ which gradually increases in time.

In this work, we show that the kinetically constrained quantum East model lies between a quantum scarred and a many-body localized system featuring an unconventional type of mobility edge in the spectrum. We name this scenario $\textit{state-dependent}$ mobility edge: while the system does not exhibit a sharp separation in energy between thermal and non-thermal eigenstates, the abundance of non-thermal eigenstates results in slow entanglement growth for $\textit{many}$ initial states, such as product states, below a finite energy density. We characterize the state-dependent mobility edge by looking at the complexity of classically simulating dynamics using tensor network for system sizes well beyond those accessible via exact diagonalization. Focusing on initial product states, we observe a qualitative change in the dynamics of the bond dimension needed as a function of their energy density. Specifically, the bond dimension typically grows $\textit{polynomially}$ in time up to a certain energy density, where we locate the state-dependent mobility edge, enabling simulations for long times. Above this energy density, the bond dimension typically grows $\textit{exponentially}$ making the simulation practically unfeasible beyond short times, as generally expected in interacting theories. We correlate the polynomial growth of the bond dimension to the presence of many non-thermal eigenstates around that energy density, a subset of which we compute via tensor network.

Landau's Fermi-liquid (FL) theory has been successful at the phenomenological description of the normal phase of many different Fermi systems. Using a dilute atomic Fermi fluid with tunable interactions, we investigate the microscopic basis of Landau's theory with a system describable from first principles. We study transport properties of an interacting Fermi gas by measuring its density response to a periodic external perturbation. In an ideal Fermi gas, we measure for the first time the celebrated Lindhard function. As the system is brought from the collisionless to the hydrodynamic regime, we observe the emergence of sound, and find that the experimental observations are quantitatively understood with a first-principle transport equation for the FL. When the system is more strongly interacting, we find deviations from such predictions. Finally, we observe the shape of the quasiparticle excitations directly from momentum-space tomography and see how it evolves from the collisionless to the collisional regime. Our study establishes this system as a clean platform for studying Landau's theory of the FL and paves the way for extending the theory to more exotic conditions, such as nonlinear dynamics and FLs with strong correlations in versatile settings.

The many-body mobility edge (MBME) in energy, which separates thermal states from many-body localization (MBL) states, is a critical yet controversial concept in many-body systems. Here we examine the quasiperiodic $t_1-t_2$ model that features a mobility edge. With the addition of nearest-neighbor interactions, we demonstrate the potential existence of a MBME. Then we investigate the impact of a type of bond dissipation on the many-body system by calculating the steady-state density matrix and analyzing the transport behavior, and demonstrate that dissipation can cause the system to predominantly occupy either the thermal region or the MBL region, irrespective of the initial state. Finally, we discuss the effects of increasing system size. Our results indicate that dissipation can induce transitions between thermal and MBL states, providing a new approach for experimentally determining the existence of the MBME.

Interactions between organisms are mediated by an intricate network of physico-chemical substances and other organisms. Understanding the dynamics of mediators and how they shape the population spatial distribution is key to predict ecological outcomes and how they would be transformed by changes in environmental constraints. However, due to the inherent complexity involved, this task is often unfeasible, from the empirical and theoretical perspectives. In this paper, we make progress in addressing this central issue, creating a bridge that provides a two-way connection between the features of the ensemble of underlying mediators and the wrinkles in the population density induced by a landscape defect (or spatial perturbation). The bridge is constructed by applying the Feynman-Vernon decomposition, which disentangles the influences among the focal population and the mediators in a compact way. This is achieved though an interaction kernel, which effectively incorporates the mediators' degrees of freedom, explaining the emergence of nonlocal influence between individuals, an ad hoc assumption in modeling population dynamics. Concrete examples are worked out and reveal the complexity behind a possible top-down inference procedure.

Diffusion is a dissipative transport phenomenon ubiquitously present in nature. Its details can now be analysed with modern effective field theory (EFT) techniques that use the closed-time-path (or Schwinger-Keldysh) formalism. We discuss the structure of the diffusive effective action appropriate for the analysis of stochastic or thermal loop effects, responsible for the so-called long-time tails, to all orders. We also elucidate and prove a number of properties of the EFT and use the theory to establish the analytic structure of the $n$-loop contributions to diffusive retarded two-point functions. Our analysis confirms a previously proposed result by Delacr\'{e}taz that used microscopic conformal field theory arguments. Then, we analyse a number of implications of these loop corrections to the dispersion relations of the diffusive mode and new, gapped modes that appear when the EFT is treated as exact. Finally, we discuss certain features of an all-loop model of diffusion that only retains a special subset of $n$-loop `banana' diagrams.

The strengths of glasses are intricately linked to their atomic-level heterogeneity. Atomistic simulations are frequently used to investigate the statistical physics of this relationship, compensating for the limited spatiotemporal resolution in experimental studies. However, theoretical insights are limited by the complexity of glass structures and the accuracy of the interatomic potentials used in simulations. Here, we investigate the strengths and fracture mechanisms of 2D silica, with all structural units accessible to direct experimental observation. We develop a neural network force field for fracture (NN-F$^{3}$) based on the deep potential-smooth edition (DeepPot-SE) framework. Representative atomic structures across crystals, nanocrystalline, paracrystalline, and continuous random network glasses are studied. We find that the virials or bond lengths control the initialization of bond-breaking events, creating nanoscale voids in the vitreous network. However, the voids do not necessarily lead to crack propagation due to a disorder-trapping effect, which is stronger than the lattice-trapping effect in a crystalline lattice and occurs over larger length and time scales. Fracture initiation proceeds with void growth and coalescence, and advances through a bridging mechanism. The fracture patterns are shaped by subsequent trapping and cleavage steps, often guided by voids forming ahead of the crack tip. These heterogeneous processes result in atomically smooth facets in crystalline regions and rough, amorphous edges in the glassy phase. These insights into 2D crystals and glasses, both sharing SiO$_{2}$ chemistry, highlight the pivotal role of atomic-level structures in determining fracture kinetics and path selection in materials.

Disordered systems theory provides powerful tools to analyze the generic behaviors of highdimensional systems, such as species-rich ecological communities or neural networks. By assuming randomness in their interactions, universality ensures that many microscopic details are irrelevant to system-wide dynamics; but the choice of a random ensemble still limits the generality of results. We show here, in the context of ecological dynamics, that these analytical tools do not require a specific choice of ensemble, and that solutions can be found based only on a fundamental rotational symmetry in the interactions, encoding the idea that traits can be recombined into new species without altering global features. Dynamical outcomes then depend on the spectrum of the interaction matrix as a free parameter, allowing us to bridge between results found in different models of interactions, and extend beyond them to previously unidentified behaviors. The distinctive feature of ecological models is the possibility of species extinctions, which leads to an increased universality of dynamics as the fraction of extinct species increases. We expect that these findings can inform new developments in theoretical ecology as well as for other families of complex systems.

The theory of spin models intersects with complex systems, graph theory, combinatorial optimization, computational complexity and neural networks. Many ensuing applications rely on the fact that complicated spin models can be transformed to simpler ones. What is the structure of such transformations? Here, we provide a framework to study and construct emulations between spin models. A spin model is a set of spin systems, and emulations are efficiently computable simulations with arbitrary energy cut-off, where a source spin system simulates a target system if, below the cut-off, the target Hamiltonian is encoded in the source Hamiltonian. We prove that emulations preserve important properties, as they induce reductions between computational problems such as computing ground states, approximating partition functions and approximate sampling from Boltzmann distributions. Emulations are modular (they can be added, scaled and composed), and allow for universality, i.e. certain spin models have maximal reach. We prove that a spin model is universal if and only if it is scalable, closed and functional complete. Because the characterization is constructive, it provides a step-by-step guide to construct emulations. We prove that the 2d Ising model with fields is universal, for which we also provide two new crossing gadgets. Finally, we show that simulations can be computed by linear programs. While some ideas of this work are contained in [1], we provide new definitions and theorems. This framework provides a toolbox for applications involving emulations of spin models.

Dissipative many-body quantum dynamics can feature strong symmetries which give rise to conserved quantities. We discuss here how a strong symmetry in conjunction with a nonequilibrium phase transition allows to devise a protocol for measuring collective many-body observables. To demonstrate this idea we consider a collective spin system whose constituents are governed by a dissipative dynamics that conserves the total angular momentum. We show that by continuously monitoring the system output the value of the total angular momentum can be inferred directly from the time-integrated emission signal, without the need of repeated projective measurements or reinitializations of the spins. This may offer a route towards the measurement of collective properties in qubit ensembles, with applications in quantum tomography, quantum computation and quantum metrology.

It is well established that the brain spontaneously traverses through a very large number of states. Nevertheless, despite its relevance to understanding brain function, a formal description of this phenomenon is still lacking. To this end, we introduce a machine learning based method allowing for the determination of the probabilities of all possible states at a given coarse-graining, from which all the thermodynamics can be derived. This is a challenge not unique to the brain, since similar problems are at the heart of the statistical mechanics of complex systems. This paper uncovers a linear scaling of the entropies and energies of the brain states, a behaviour first conjectured by Hagedorn to be typical at the limiting temperature in which ordinary matter disintegrates into quark matter. Equivalently, this establishes the existence of a Zipf law scaling underlying the appearance of a wide range of brain states. Based on our estimation of the density of states for large scale functional magnetic resonance imaging (fMRI) human brain recordings, we observe that the brain operates asymptotically at the Hagedorn temperature. The presented approach is not only relevant to brain function but should be applicable for a wide variety of complex systems.

We present a self-consistent theory for sound propagation in a simple model of a disordered solid. The solid is modeled as a collection of randomly distributed particles connected by harmonic springs with strengths that depend on the interparticle distances, i.e the Euclidean random matrix model of Mezard et al. [Nucl. Phys. 559B, 689 (1999)]. The derivation of the theory combines two exact projection operator steps and a factorization approximation. Within our approach the square of the speed of sound is non-negative. The unjamming transition manifests itself through vanishing of the speed of sound.