Showing papers in "arXiv: Statistical Mechanics in 2020"

The nonperturbative functional renormalization group and its applications

Abstract: The renormalization group plays an essential role in many areas of physics, both conceptually and as a practical tool to determine the long-distance low-energy properties of many systems on the one hand and on the other hand search for viable ultraviolet completions in fundamental physics. It provides us with a natural framework to study theoretical models where degrees of freedom are correlated over long distances and that may exhibit very distinct behavior on different energy scales. The nonperturbative functional renormalization-group (FRG) approach is a modern implementation of Wilson's RG, which allows one to set up nonperturbative approximation schemes that go beyond the standard perturbative RG approaches. The FRG is based on an exact functional flow equation of a coarse-grained effective action (or Gibbs free energy in the language of statistical mechanics). We review the main approximation schemes that are commonly used to solve this flow equation and discuss applications in equilibrium and out-of-equilibrium statistical physics, quantum many-particle systems, high-energy physics and quantum gravity.

Measurement Protected Quantum Phases

Abstract: We introduce a class of hybrid quantum circuits, with random unitaries and projective measurements, which host long-range order in the area law entanglement phase of the steady state. Our primary example is circuits with unitaries respecting a global Ising symmetry and two competing types of measurements. The phase diagram has an area law phase with spin glass order, which undergoes a direct transition to a paramagnetic phase with volume law entanglement, as well as a critical regime. Using mutual information diagnostics, we find that such entanglement transitions preserving a global symmetry are in new universality classes. We analyze generalizations of such hybrid circuits to higher dimensions, which allow for coexistence of order and volume law entanglement, as well as topological order without any symmetry restrictions.

Exponentially faster cooling in a colloidal system

Abstract: As the temperature of a cooling object decreases as it relaxes to thermal equilibrium, it is intuitively assumed that a hot object should take longer to cool than a warm one. Yet, some 2,300 years ago, Aristotle observed that "to cool hot water quickly, begin by putting it in the sun". In the 1960s, this counterintuitive phenomenon was rediscovered as the statement that "hot water can freeze faster than cold water" and has become known as the Mpemba effect; it has since been the subject of much experimental investigation and some controversy. Although many specific mechanisms have been proposed, no general consensus exists as to the underlying cause. Here we demonstrate the Mpemba effect in a controlled setting - the thermal quench of a colloidal system immersed in water, which serves as a heat bath. Our results are reproducible and agree quantitatively with calculations based on a recently proposed theoretical framework. By carefully choosing parameters, we observe cooling that is exponentially faster than that observed using typical parameters, in accord with the recently predicted strong Mpemba effect. Our experiments outline the generic conditions needed to accelerate heat removal and relaxation to thermal equilibrium and support the idea that the Mpemba effect is not simply a scientific curiosity concerning how water freezes into ice - one of the many anomalous features of water - but rather the prototype for a wide range of anomalous relaxation phenomena of broad technological importance.

Topological Defects on the Lattice: Dualities and Degeneracies

Abstract: We construct topological defects in two-dimensional classical lattice models and quantum chains The defects satisfy local commutation relations guaranteeing that the partition function is independent of their path These relations and their solutions are extended to allow defect lines to fuse, branch and satisfy all the properties of a fusion category We show how the two-dimensional classical lattice models and their topological defects are naturally described by boundary conditions of a Turaev-Viro-Barrett-Westbury partition function These defects allow Kramers-Wannier duality to be generalized to a large class of models, explaining exact degeneracies between non-symmetry-related ground states as well as in the low-energy spectrum They give a precise and general notion of twisted boundary conditions and the universal behaviour under Dehn twists Gluing a topological defect to a boundary yields linear identities between partition functions with different boundary conditions, allowing ratios of the universal g-factor to be computed exactly on the lattice We develop this construction in detail in a variety of examples, including the Potts, parafermion and height models

Superuniversality of superdiffusion.

Abstract: Anomalous finite-temperature transport has recently been observed in numerical studies of various integrable models in one dimension; these models share the feature of being invariant under a continuous non-abelian global symmetry. This work offers a comprehensive group-theoretic account of this elusive phenomenon. For an integrable quantum model invariant under a global non-abelian simple Lie group $G$, we find that finite-temperature transport of Noether charges associated with symmetry $G$ in thermal states that are invariant under $G$ is universally superdiffusive and characterized by dynamical exponent $z = 3/2$. This conclusion holds regardless of the Lie algebra symmetry, local degrees of freedom (on-site representations), Lorentz invariance, or particular realization of microscopic interactions: we accordingly dub it as superuniversal. The anomalous transport behavior is attributed to long-lived giant quasiparticles dressed by thermal fluctuations. We provide an algebraic viewpoint on the corresponding dressing transformation and elucidate formal connections to fusion identities amongst the quantum-group characters. We identify giant quasiparticles with nonlinear soliton modes of classical field theories that describe low-energy excitations above ferromagnetic vacua. Our analysis of these field theories also provides a complete classification of the low-energy (i.e., Goldstone-mode) spectra of quantum isotropic ferromagnetic chains.

Optimal mean first-passage time for a Brownian searcher subjected to resetting: experimental and theoretical results

Abstract: We study experimentally and theoretically the optimal mean time needed by a free diffusing Brownian particle to reach a target at a distance L from an initial position in the presence of resetting. Both the initial position and the resetting position are Gaussian distributed with width $\sigma$. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different values of the ratio b = L/$\sigma$ and find an interesting phase transtion at a critical value b = bc. For bc < b < $\infty$, there is a metastable optimum time which disappears for b < bc. The intrinsic diffculties in implementing these protocols in experiments are also discussed.

Symmetry resolved entanglement entropy of excited states in a CFT

Abstract: We report a throughout analysis of the entanglement entropies related to different symmetry sectors in the low-lying primary excited states of a conformal field theory (CFT) with an internal U(1) symmetry. Our findings extend recent results for the ground state. We derive a general expression for the charged moments, i.e. the generalised cumulant generating function, which can be written in terms of correlation functions of the operator that define the state through the CFT operator-state correspondence. We provide explicit analytic computations for the compact boson CFT (aka Luttinger liquid) for the vertex and derivative excitations. The Fourier transform of the charged moments gives the desired symmetry resolved entropies. At the leading order, they satisfy entanglement equipartition, as in the ground state, but we find, within CFT, subleading terms that break it. Our analytical findings are checked against free fermions calculations on a lattice, finding excellent agreement. As a byproduct, we have exact results for the full counting statistics of the U(1) charge in the considered excited states.

Local detailed balance

Abstract: We review the physical meaning and mathematical implementation of the condition of local detailed balance for a class of nonequilibrium mesoscopic processes. A central concept is that of fluctuating entropy flux for which the steady average gives the mean entropy production rate. We repeat how local detailed balance is essentially equivalent to the widely discussed fluctuation relations for that entropy flux and hence is at most "only half of the story."

Symmetry resolved entanglement in two-dimensional systems via dimensional reduction

Abstract: We report on the calculation of the symmetry resolved entanglement entropies in two-dimensional many-body systems of free bosons and fermions by \emph{dimensional reduction}. When the subsystem is translational invariant in a transverse direction, this strategy allows us to reduce the initial two-dimensional problem into decoupled one-dimensional ones in a mixed space-momentum representation. While the idea straightforwardly applies to any dimension $d$, here we focus on the case $d=2$ and derive explicit expressions for two lattice models possessing a $U(1)$ symmetry, i.e., free non-relativistic massless fermions and free complex (massive and massless) bosons. Although our focus is on symmetry resolved entropies, some results for the total entanglement are also new. Our derivation gives a transparent understanding of the well known different behaviours between massless bosons and fermions in $d\geq2$: massless fermions presents logarithmic violation of the area which instead strictly hold for bosons, even massless. This is true both for the total and the symmetry resolved entropies. Interestingly, we find that the equipartition of entanglement into different symmetry sectors holds also in two dimensions at leading order in subsystem size; we identify for both systems the first term breaking it. All our findings are quantitatively tested against exact numerical calculations in lattice models for both bosons and fermions.

Inertial self-propelled particles.

Abstract: We study a self-propelled particle moving in a solvent with the active Ornstein Uhlenbeck dynamics in the underdamped regime to evaluate the influence of the inertia. We focus on the properties of potential-free and harmonically confined underdamped active particles, studying how the single-particle trajectories modify for different values of the drag coefficient. In both cases, we solve the dynamics in terms of correlation matrices and steady-state probability distribution functions revealing the explicit correlations between velocity and active force. We also evaluate the influence of the inertia on the time-dependent properties of the system, discussing the mean square displacement and the time-correlations of particle positions and velocities. Particular attention is devoted to the study of the Virial active pressure unveiling the role of the inertia on this observable.

A reinforcement learning approach to rare trajectory sampling

Abstract: Very often when studying non-equilibrium systems one is interested in analysing dynamical behaviour that occurs with very low probability, so called rare events. In practice, since rare events are by definition atypical, they are often difficult to access in a statistically significant way. What are required are strategies to "make rare events typical" so that they can be generated on demand. Here we present such a general approach to adaptively construct a dynamics that efficiently samples atypical events. We do so by exploiting the methods of reinforcement learning (RL), which refers to the set of machine learning techniques aimed at finding the optimal behaviour to maximise a reward associated with the dynamics. We consider the general perspective of dynamical trajectory ensembles, whereby rare events are described in terms of ensemble reweighting. By minimising the distance between a reweighted ensemble and that of a suitably parametrised controlled dynamics we arrive at a set of methods similar to those of RL to numerically approximate the optimal dynamics that realises the rare behaviour of interest. As simple illustrations we consider in detail the problem of excursions of a random walker, for the case of rare events with a finite time horizon; and the problem of a studying current statistics of a particle hopping in a ring geometry, for the case of an infinite time horizon. We discuss natural extensions of the ideas presented here, including to continuous-time Markov systems, first passage time problems and non-Markovian dynamics.

Criticality and entanglement in non-unitary quantum circuits and tensor networks of non-interacting fermions

Abstract: Models for non-unitary quantum dynamics, such as quantum circuits that include projective measurements, have been shown to exhibit rich quantum critical behavior. There are many complementary perspectives on this behavior. For example, there is a known correspondence between d-dimensional local non-unitary quantum circuits and tensor networks on a D=(d+1)-dimensional lattice. Here, we show that in the case of systems of non-interacting fermions, there is furthermore a full correspondence between non-unitary circuits in d spatial dimensions and unitary non-interacting fermion problems with static Hermitian Hamiltonians in D=(d+1) spatial dimensions. This provides a powerful new perspective for understanding entanglement phases and critical behavior exhibited by non-interacting circuits. Classifying the symmetries of the corresponding non-interacting Hamiltonian, we show that a large class of random circuits, including the most generic circuits with randomness in space and time, are in correspondence with Hamiltonians with static spatial disorder in the ten Altland-Zirnbauer symmetry classes. We find the criticality that is known to occur in all of these classes to be the origin of the critical entanglement properties of the corresponding random non-unitary circuit. To exemplify this, we numerically study the quantum states at the boundary of Haar-random Gaussian fermionic tensor networks of dimension D=2 and D=3. We show that the most general such tensor network ensemble corresponds to a unitary problem of non-interacting fermions with static disorder in Altland-Zirnbauer symmetry class DIII, which for both D=2 and D=3 is known to exhibit a stable critical metallic phase. Tensor networks and corresponding random non-unitary circuits in the other nine Altland-Zirnbauer symmetry classes can be obtained from the DIII case by implementing Clifford algebra extensions for classifying spaces.

Hydrodynamics in lattice models with continuous non-Abelian symmetries

Abstract: We develop a systematic effective field theory of hydrodynamics for many-body systems on the lattice with global continuous non-Abelian symmetries. Models with continuous non-Abelian symmetries are ubiquitous in physics, arising in diverse settings ranging from hot nuclear matter to cold atomic gases and quantum spin chains. In every dimension and for every flavor symmetry group, the low energy theory is a set of coupled noisy diffusion equations. Independence of the physics on the choice of canonical or microcanonical ensemble is manifest in our hydrodynamic expansion, even though the ensemble choice causes an apparent shift in quasinormal mode spectra. We use our formalism to explain why flavor symmetry is qualitatively different from hydrodynamics with other non-Abelian conservation laws, including angular momentum and charge multipoles. As a significant application of our framework, we study spin and energy diffusion in classical one-dimensional SU(2)-invariant spin chains, including the Heisenberg model along with multiple generalizations. We argue based on both numerical simulations and our effective field theory framework that non-integrable spin chains on a lattice exhibit conventional spin diffusion, in contrast to some recent predictions that diffusion constants grow logarithmically at late times. We show that the apparent enhancement of diffusion is due to slow equilibration caused by (non-Abelian) hydrodynamic fluctuations.

Multiscale Thermodynamics.

Abstract: Multiscale thermodynamics is a theory of relations among levels of investigation of complex systems. It includes the classical equilibrium thermodynamics as a special case but it is applicable to both static and time evolving processes in externally and internally driven macroscopic systems that are far from equilibrium and are investigated on microscopic, mesoscopic, and macroscopic levels. In this paper we formulate the multiscale thermodynamics, explain its origin, and illustrate it in mesoscopic dynamics that combines levels.

Diffusion with resetting in a logarithmic potential

Abstract: We study the effect of resetting on diffusion in a logarithmic potential. In this model, a particle diffusing in a potential $U(x) = U_0\log|x|$ is reset, i.e., taken back to its initial position, with a constant rate $r$. We show that this analytically tractable model system exhibits a series of phase transitions as a function of a single parameter, $\beta U_0$, the ratio of the strength of the potential to the thermal energy. For $\beta U_0 -1$ the potential is either weakly repulsive or attractive and the diffusing particle eventually reaches the origin. In this case, we provide a closed form expression for the subsequent first-passage time distribution and show that a resetting transition occurs at $\beta U_0=5$. Namely, we find that resetting can expedite arrival to the origin when $-1 5$. The results presented herein generalize results for simple diffusion with resetting -- a widely applicable model that is obtained from ours by setting $U_0=0$. Extending to general potential strengths, our work opens the door to theoretical and experimental investigation of a plethora of problems that bring together resetting and diffusion in logarithmic potential.

Large work extraction and the Landauer limit in a continuous Maxwell demon

Abstract: The relation between entropy and information dates back to the classical Maxwell demon (MD) paradox, a thought experiment proposed in 1867 by J. C. Maxwell to violate the second law of thermodynamics. A variant of the classical MD is the Szilard engine proposed by L. Szilard in 1926 in which the demon observes, at a given time, the compartment occupied by a single molecule in a vessel and extracts work by operating a pulley device. Here we introduce the Continuous Maxwell Demon (CMD), a device capable of extracting arbitrarily large amounts of work per cycle by repeated measurements of the state of a system, and experimentally test it in single DNA hairpin pulling experiments. In the CMD the demon monitors the state of the DNA hairpin (folded or unfolded) by observing it at equally spaced time intervals but extracts work only when the molecule changes state. We demonstrate that the average maximum work per cycle that can be extracted by the CMD is limited by the information-content of the stored sequences, in agreement with the second law. Work extraction efficiency is found to be maximal in the large information-content limit where work extraction is fuelled by rare events.

Interacting particle solutions of Fokker-Planck equations through gradient-log-density estimation

Abstract: Fokker-Planck equations are extensively employed in various scientific fields as they characterise the behaviour of stochastic systems at the level of probability density functions. Although broadly used, they allow for analytical treatment only in limited settings, and often is inevitable to resort to numerical solutions. Here, we develop a computational approach for simulating the time evolution of Fokker-Planck solutions in terms of a mean field limit of an interacting particle system. The interactions between particles are determined by the gradient of the logarithm of the particle density, approximated here by a novel statistical estimator. The performance of our method shows promising results, with more accurate and less fluctuating statistics compared to direct stochastic simulations of comparable particle number. Taken together, our framework allows for effortless and reliable particle-based simulations of Fokker-Planck equations in low and moderate dimensions. The proposed gradient-log-density estimator is also of independent interest, for example, in the context of optimal control.

Emergent memory and kinetic hysteresis in strongly driven networks

Abstract: Stochastic network-dynamics are typically assumed to be memory-less. Involving prolonged dwells interrupted by instantaneous transitions between nodes such Markov networks stand as a coarse-graining paradigm for chemical reactions, gene expression, molecular machines, spreading of diseases, protein dynamics, diffusion in energy landscapes, epigenetics and many others. However, as soon as transitions cease to be negligibly short, as often observed in experiments, the dynamics develops a memory. That is, state-changes depend not only on the present state but also on the past. Here, we establish the first thermodynamically consistent -- dissipation-preserving -- mapping of continuous dynamics onto a network, which reveals ingrained dynamical symmetries and an unforeseen kinetic hysteresis. These symmetries impose three independent sources of fluctuations in state-to state kinetics that determine the `flavor of memory'. The hysteresis between the forward/backward in time coarse-graining of continuous trajectories implies a new paradigm for the thermodynamics of active molecular processes in the presence of memory, that is, beyond the assumption of local detailed balance. Our results provide a new understanding of fluctuations in the operation of molecular machines as well as catch-bonds involved in cellular adhesion.

Spectral statistics in constrained many-body quantum chaotic systems

Abstract: We study the spectral statistics of spatially-extended many-body quantum systems with on-site Abelian symmetries or local constraints, focusing primarily on those with conserved dipole and higher moments. In the limit of large local Hilbert space dimension, we find that the spectral form factor $K(t)$ of Floquet random circuits can be mapped exactly to a classical Markov circuit, and, at late times, is related to the partition function of a frustration-free Rokhsar-Kivelson (RK) type Hamiltonian. Through this mapping, we show that the inverse of the spectral gap of the RK-Hamiltonian lower bounds the Thouless time $t_{\mathrm{Th}}$ of the underlying circuit. For systems with conserved higher moments, we derive a field theory for the corresponding RK-Hamiltonian by proposing a generalized height field representation for the Hilbert space of the effective spin chain. Using the field theory formulation, we obtain the dispersion of the low-lying excitations of the RK-Hamiltonian in the continuum limit, which allows us to extract $t_{\mathrm{Th}}$. In particular, we analytically argue that in a system of length $L$ that conserves the $m^{th}$ multipole moment, $t_{\mathrm{Th}}$ scales subdiffusively as $L^{2(m+1)}$. We also show that our formalism directly generalizes to higher dimensional circuits, and that in systems that conserve any component of the $m^{th}$ multipole moment, $t_{\mathrm{Th}}$ has the same scaling with the linear size of the system. Our work therefore provides a general approach for studying spectral statistics in constrained many-body chaotic systems.

Experimental evidence of symmetry breaking of transition-path times

Abstract: While Kramers' rates have been studied for almost a century, the transition path time between states has only recently received attention. Transition paths between different energy levels are expected to be indistinguishable in shape and have equal uphill and downhill lengths. This fundamental symmetry often prevents directional sensing in experiments. Here, we report experimental evidence for transition path time symmetry and its breakdown on the single-molecule and mesoscale. In automated Brownian dynamics experiments, we establish first-passage time symmetries of colloids driven by femtoNewton-range forces in holographically-created optical landscapes confined in microchannels. Conversely, we show that transitions which couple in a path-dependent manner to fluctuating forces exhibit asymmetry. We reproduce this asymmetry in folding transitions of DNA-hairpins driven out of equilibrium and suggest a topological mechanism for the symmetry breakdown. Our results can reveal directionality in molecular transitions or translocations through membrane channels and nanopores.

Topological order and criticality in (2+1)D monitored random quantum circuits

Abstract: It has recently been discovered that random quantum circuits provide an avenue to realize rich entanglement phase diagrams, which are hidden to standard expectation values of operators. Here we study (2+1)D random circuits with random Clifford unitary gates and measurements designed to stabilize trivial area law and topologically ordered phases. With competing single qubit Pauli-Z and toric code stabilizer measurements, in addition to random Clifford unitaries, we find a phase diagram involving a tricritical point that maps to (2+1)D percolation, a possibly stable critical phase, topologically ordered, trivial, and volume law phases, and lines of critical points separating them. With Pauli-Y single qubit measurements instead, we find an anisotropic self-dual tricritical point, with dynamical exponent $z \approx 1.46$, exhibiting logarithmic violation of the area law and an anomalous exponent for the topological entanglement entropy, which thus appears distinct from any known percolation fixed point. The phase diagram also hosts a measurement-induced volume law entangled phase in the absence of unitary dynamics. We also propose a definition of equivalent area law entangled phases of quantum random circuits which is motivated by our numerical results.

Boundary effects on symmetry resolved entanglement

Abstract: We study the symmetry resolved entanglement entropies in one-dimensional systems with boundaries. We provide some general results for conformal invariant theories and then move to a semi-infinite chain of free fermions. We consider both an interval starting from the boundary and away from it. We derive exact formulas for the charged and symmetry resolved entropies based on theorems and conjectures about the spectra of Toeplitz+Hankel matrices. En route to characterise the interval away from the boundary, we prove a general relation between the eigenvalues of Toeplitz+Hankel matrices and block Toeplitz ones. An important aspect is that the saddle-point approximation from charged to symmetry resolved entropies introduces algebraic corrections to the scaling that are much more severe than in systems without boundaries.

From tunnels to towers: quantum scars from Lie Algebras and q-deformed Lie Algebras

Abstract: We present a general symmetry-based framework for obtaining many-body Hamiltonians with scarred eigenstates that do not obey the eigenstate thermalization hypothesis. Our models are derived from parent Hamiltonians with a non-Abelian (or q-deformed) symmetry, whose eigenspectra are organized as degenerate multiplets that transform as irreducible representations of the symmetry (`tunnels'). We show that large classes of perturbations break the symmetry, but in a manner that preserves a particular low-entanglement multiplet of states -- thereby giving generic, thermal spectra with a `shadow' of the broken symmetry in the form of scars. The generators of the Lie algebra furnish operators with `spectrum generating algebras' that can be used to lift the degeneracy of the scar states and promote them to equally spaced `towers'. Our framework applies to several known models with scars, but we also introduce new models with scars that transform as irreducible representations of symmetries such as SU(3) and $q$-deformed SU(2), significantly generalizing the types of systems known to harbor this phenomenon. Additionally, we present new examples of generalized AKLT models with scar states that do not transform in an irreducible representation of the relevant symmetry. These are derived from parent Hamiltonians with enhanced symmetries, and bring AKLT-like models into our framework.

Random walks and community detection in hypergraphs

Abstract: We propose a one parameter family of random walk processes on hypergraphs, where a parameter biases the dynamics of the walker towards hyperedges of low or high cardinality. We show that for each value of the parameter the resulting process defines its own hypergraph projection on a weighted network. We then explore the differences between them by considering the community structure associated to each random walk process. To do so, we generalise the Markov stability framework to hypergraphs and test it on artificial and real-world hypergraphs.

How to "measure" a structural relaxation time that is too long to be measured?

Abstract: It has recently become possible to prepare ultrastable glassy materials characterised by structural relaxation times which vastly exceed the duration of any feasible experiment. Similarly, new algorithms have led to the production of ultrastable computer glasses. Is it possible to obtain a reliable estimate of a structural relaxation time that is too long to be measured? We review, organise, and critically discuss various methods to estimate very long relaxation times. We also perform computer simulations of three dimensional ultrastable hard spheres glasses to test and quantitatively compare some of these methods for a single model system. The various estimation methods disagree significantly and it is not yet clear how to accurately estimate extremely long relaxation times.

Emergence of molecular friction in liquids: bridging between the atomistic and hydrodynamic pictures

Abstract: Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with rigorous theoretical arguments, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids; as a consequence, its origin is non-local in time. Concomitantly, the molecules experience Brownian forces that display persistent correlations and long-lasting memory. A critical test of the generalised Stokes-Einstein relation, mapping the friction of single molecules to the viscoelastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield novel insights into vitrification mechanisms and the intriguing mechanical properties of soft materials.

Nonlinear stochastic modeling with Langevin regression

Abstract: Many physical systems characterized by nonlinear multiscale interactions can be effectively modeled by treating unresolved degrees of freedom as random fluctuations. However, even when the microscopic governing equations and qualitative macroscopic behavior are known, it is often difficult to derive a stochastic model that is consistent with observations. This is especially true for systems such as turbulence where the perturbations do not behave like Gaussian white noise, introducing non-Markovian behavior to the dynamics. We address these challenges with a framework for identifying interpretable stochastic nonlinear dynamics from experimental data, using both forward and adjoint Fokker-Planck equations to enforce statistical consistency. If the form of the Langevin equation is unknown, a simple sparsifying procedure can provide an appropriate functional form. We demonstrate that this method can effectively learn stochastic models in two artificial examples: recovering a nonlinear Langevin equation forced by colored noise and approximating the second-order dynamics of a particle in a double-well potential with the corresponding first-order bifurcation normal form. Finally, we apply the proposed method to experimental measurements of a turbulent bluff body wake and show that the statistical behavior of the center of pressure can be described by the dynamics of the corresponding laminar flow driven by nonlinear state-dependent noise.

Quantum-mechanical exploration of the phase diagram of water

Abstract: The phase diagram of water harbours many mysteries: some of the phase boundaries are fuzzy, and the set of known stable phases may not be complete. Starting from liquid water and a comprehensive set of 50 ice structures, we compute the phase diagram at three hybrid density-functional-theory levels of approximation, accounting for thermal and nuclear fluctuations as well as proton disorder. Such calculations are only made tractable because we combine machine-learning methods and advanced free-energy techniques. The computed phase diagram is in qualitative agreement with experiment, particularly at pressures $\lesssim$8000 bar, and the discrepancy in chemical potential is comparable with the subtle uncertainties introduced by proton disorder and the spread between the three hybrid functionals. None of the hypothetical ice phases considered is thermodynamically stable in our calculations, suggesting the completeness of the experimental water phase diagram in the region considered. Our work demonstrates the feasibility of predicting the phase diagram of a polymorphic system from first principles and provides a thermodynamic way of testing the limits of quantum-mechanical calculations.

General description for nonequilibrium steady states in periodically driven dissipative quantum systems

Abstract: Laser technology has developed and accelerated photo-induced nonequilibrium physics from both scientific and engineering viewpoints. The Floquet engineering, i.e., controlling material properties and functionalities by time-periodic drives, is a forefront of quantum physics of light-matter interaction, but limited to ideal dissipationless systems. For the Floquet engineering extended to a variety of materials, it is vital to understand the quantum states emerging in a balance of the periodic drive and energy dissipation. Here we derive the general description for nonequilibrium steady states (NESS) in periodically driven dissipative systems by focusing on the systems under high-frequency drive and time-independent Lindblad-type dissipation with the detailed balance condition. Our formula correctly describes the time-average, fluctuation, and symmetry property of the NESS, and can be computed efficiently in numerical calculations. Our approach will play fundamental roles in Floquet engineering in a broad class of dissipative quantum systems such as atoms and molecules, mesoscopic systems, and condensed matter.

Active Brownian Heat Engines

Abstract: When do non-equilibrium forms of disordered energy qualify as heat? \textcolor{blue}{We address this question in the context of cyclically operating heat engines in contact with a non-equilibrium energy reservoir that defies the zeroth law of thermodynamics. To consistently address the latter as a heat bath requires the existence of a precise mapping to an equivalent cycle with an equilibrium bath at a time-dependent effective temperature. We identify the most general setup for which this can generically be ascertained and thoroughly discuss an analytically tractable, experimentally relevant scenario}: a Brownian particle confined in a \textcolor{blue}{periodically} modulated harmonic potential and coupled to some non-equilibrium bath of variable activity. We deduce formal limitations for its thermodynamic performance, including maximum efficiency, efficiency at maximum power, and maximum efficiency at fixed power. They can guide the design of new micro-machines and clarify how much these can outperform passive-bath designs, which has been a debated issue for recent experimental realizations. To illustrate the general principles for practical quasi-static and finite-rate protocols, we further analyze a specific realization of such an active heat engine based on the paradigmatic Active Brownian Particle (ABP) model. This reveals some non-intuitive features of the explicitly computed dynamical effective temperature, illustrates various conceptual and practical limitations of the effective-equilibrium mapping, and clarifies the operational relevance of various coarse-grained measures of dissipation.