Mori-Zwanzig projection operator formalism for far-from-equilibrium systems with time-dependent Hamiltonians

Microscopic derivation of the thin film equation using the Mori-Zwanzig formalism

The hydrodynamics of thin films is typically described using macroscopic models whose connection to the microscopic particle dynamics is a subject of ongoing research. Existing methods based on density functional theory provide a good description of static thin films but are not sufficient for understanding nonequilibrium dynamics. In this work, we present a microscopic derivation of the thin film equation using the Mori-Zwanzig projection operator formalism. This method allows to directly obtain the correct gradient dynamics structure along with microscopic expressions for mobility and free energy. Our results are verified against molecular dynamics simulations for both simple fluids and polymers.

Derivation of the nonequilibrium generalized Langevin equation from a time-dependent many-body Hamiltonian

It has become standard practice to describe systems that remain far from equilibrium even in their steady state by Langevin equations with colored noise which is chosen independently from the friction contribution. Since these Langevin equations are typically not derived from first-principle Hamiltonian dynamics, it is not clear whether they correspond to physically realizable scenarios. By exact Mori projection in phase space we derive the nonequilibrium generalized Langevin equation (GLE) for an arbitrary phase-space dependent observable A from a generic many-body Hamiltonian with a time-dependent external force h(t) acting on the same observable A. This is the same Hamiltonian from which the standard fluctuation-dissipation theorem is derived, which reflects the generality of our approach. The observable A could, for example, be the position of an atom, of a molecule or of a macroscopic object, the distance between two such entities or a more complex phase-space function such as the reaction coordinate of a chemical reaction or of the folding of a protein. The Hamiltonian could, for example, describe a fluid, a solid, a viscoelastic medium, or even a turbulent inhomogeneous environment. The GLE, which is a closed-form equation of motion for the observable A, is obtained in explicit form to all orders in h(t) and without restrictions on the type of many-body Hamiltonian or the observable A. If the dynamics of the observable A corresponds to a Gaussian process, the resultant GLE has a similar form as the equilibrium Mori GLE, and in particular the friction memory kernel is given by the two-point autocorrelation function of the sum of the complementary and the external force h(t). This is a nontrivial and useful result, as many observables that characterize nonequilibrium systems display Gaussian statistics. For non-Gaussian nonequilibrium observables correction terms appear in the GLE and in the relation between the force autocorrelation and the friction memory kernel, which are explicitly given in terms of cubic correlation functions of A. Interpreting the external force h(t) as a stochastic process, we derive nonequilibrium corrections to the fluctuation-dissipation theorem and present methods to extract all GLE parameters from experimental or simulation time-series data, thus making our nonequilibrium GLE a practical tool to study and model general nonequilibrium systems.

Thermodynamics of an Empty Box

A gas in a box is perhaps the most important model system studied in thermodynamics and statistical mechanics. Usually, studies focus on the gas, whereas the box merely serves as an idealized confinement. The present article focuses on the box as the central object and develops a thermodynamic theory by treating the geometric degrees of freedom of the box as the degrees of freedom of a thermodynamic system. Applying standard mathematical methods to the thermodynamics of an empty box allows equations with the same structure as those of cosmology and classical and quantum mechanics to be derived. The simple model system of an empty box is shown to have interesting connections to classical mechanics, special relativity, and quantum field theory.

Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

Generalized Langevin dynamics simulation with non-stationary memory kernels: How to make noise

We present a numerical method to produce stochastic dynamics according to the generalized Langevin equation with a non-stationary memory kernel. This type of dynamics occurs when a microscopic system with an explicitly time-dependent Liouvillian is coarse-grained by means of a projection operator formalism. We show how to replace the deterministic fluctuating force in the generalized Langevin equation by a stochastic process, such that the distributions of the observables are reproduced up to moments of a given order. Thus, in combination with a method to extract the memory kernel from simulation data of the underlying microscopic model, the method introduced here allows us to construct and simulate a coarse-grained model for a driven process.

Modeling of biomolecular machines in non-equilibrium steady states

Numerical computations have become a pillar of all modern quantitative sciences. Any computation involves modeling-even if often this step is not made explicit-and any model has to neglect details while still being physically accurate. Equilibrium statistical mechanics guides both the development of models and numerical methods for dynamics obeying detailed balance. For systems driven away from thermal equilibrium, such a universal theoretical framework is missing. For a restricted class of driven systems governed by Markov dynamics and local detailed balance, stochastic thermodynamics has evolved to fill this gap and to provide fundamental constraints and guiding principles. The next step is to advance stochastic thermodynamics from simple model systems to complex systems with tens of thousands or even millions of degrees of freedom. Biomolecules operating in the presence of chemical gradients and mechanical forces are a prime example for this challenge. In this Perspective, we give an introduction to isothermal stochastic thermodynamics geared toward the systematic multiscale modeling of the conformational dynamics of biomolecular and synthetic machines, and we outline some of the open challenges.

Mori-Zwanzig Formalism for General Relativity: A New Approach to the Averaging Problem

Cosmology relies on a coarse-grained description of the universe, assumed to be valid on large length scales. However, the nonlinearity of general relativity makes coarse graining extremely difficult. We here address this problem by extending the Mori-Zwanzig projection operator formalism, a highly successful coarse-graining method from statistical mechanics, towards general relativity. Using the Buchert equations, we derive a new dynamic equation for the Hubble parameter which captures the effects of averaging through a memory function. This gives an empirical prediction for the cosmic jerk.

Mori-Zwanzig projection operator formalism: Particle-based coarse-grained dynamics of open classical systems far from equilibrium

We present a generalized Langevin equation (GLE) of motion that governs exactly the time evolution of phase-space observables in finite open systems described by classical Hamiltonians with explicitly time-dependent potentials. This formalism is based on the Mori-Zwanzig projection operator (PO) method with a time-independent Zwanzig PO within a Heisenberg (Lagrangian) picture and reduced description of Hamiltonian systems in terms of canonical relevant and irrelevant coordinates. We demonstrate that, similarly to closed systems, GLE dynamics in Hamiltonian systems in the presence of time-dependent potentials is determined by conservative, dissipative memory, and projected force fields, and that the memory functions relate to the projected force, which is a two-time process, in a way that is reminiscent of the equilibrium second fluctuation-dissipation relation. We further show that, in the most general case, the memory kernel depends on the relevant momentum gradients of the (Boltzmann) entropy of the irrelevant subsystem. Using two Zwanzig operators which are, respectively, functionals of the canonical and generalized canonical probability densities, we then derive what we call canonical and generalized canonical GLEs. Further, we can formulate the particle-based, coarse-grained (CG) GLE dynamics by transitioning to Jacobi coordinates which corresponds to a particle set partitioning of the Hamiltonian system. The obtained canonical CG GLE of motion for the relevant momenta is a generalization of the CG equation of motion known for closed systems. Also, using a Markovian approximation of the canonical CG GLE, we can extend the dissipative particle dynamics equation to open systems. A distinctive feature of our extension is a use of explicitly time-dependent frictions, which reflect the changes in the dissipation rate caused by time-dependent coupling to an external bath. Our GLE formalism and workflow constitute a general and viable framework that can be readily used as a starting point to rigorously formulate microscopically informed CG treatments for a variety of phenomena in externally forced systems far from equilibrium.

The five problems of irreversibility

Thermodynamics has a clear arrow of time, characterized by the irreversible approach to equilibrium. This stands in contrast to the laws of microscopic theories, which are invariant under time-reversal. Foundational discussions of this "problem of irreversibility" often focus on historical considerations, and do therefore not take results of modern physical research on this topic into account. In this article, I will close this gap by studying the implications of dynamical density functional theory (DDFT), a central method of modern nonequilibrium statistical mechanics not previously considered in philosophy of physics, for this debate. For this purpose, the philosophical discussion of irreversibility is structured into five problems, concerned with the source of irreversibility in thermodynamics, the definition of equilibrium and entropy, the justification of coarse-graining, the approach to equilibrium and the arrow of time. For each of these problems, it is shown that DDFT provides novel insights that are of importance for both physicists and philosophers of physics.

Comments on the validity of the non-stationary generalized Langevin equation as a coarse-grained evolution equation for microscopic stochastic dynamics

We recently showed that the dynamics of coarse-grained observables in systems out of thermal equilibrium are governed by the non-stationary generalized Langevin equation [H. Meyer, T. Voigtmann, and T. Schilling, J. Chem. Phys. 147, 214110 (2017); 150, 174118 (2019)]. The derivation we presented in these two articles was based on the assumption that the dynamics of the microscopic degrees of freedom were deterministic. Here, we extend the discussion to stochastic microscopic dynamics. The fact that the same form of the non-stationary generalized Langevin equation as derived for the deterministic case also holds for stochastic processes implies that methods designed to estimate the memory kernel, drift term, and fluctuating force term of this equation, as well as methods designed to propagate it numerically, can be applied to data obtained in molecular dynamics simulations that employ a stochastic thermostat or barostat.

Memory-dependent noise-induced resonance and diffusion in non-Markovian systems

We study random processes with nonlocal memory and obtain solutions of the Mori-Zwanzig equation describing non-Markovian systems. We analyze the system dynamics depending on the amplitudes ν and μ_{0} of the local and nonlocal memory and pay attention to the line in the (ν, μ_{0}) plane separating the regions with asymptotically stationary and nonstationary behavior. We obtain general equations for such boundaries and consider them for three examples of nonlocal memory functions. We show that there exist two types of boundaries with fundamentally different system dynamics. On the boundaries of the first type, diffusion with memory takes place, whereas on borderlines of the second type the phenomenon of noise-induced resonance can be observed. A distinctive feature of noise-induced resonance in the systems under consideration is that it occurs in the absence of an external regular periodic force. It takes place due to the presence of frequencies in the noise spectrum, which are close to the self-frequency of the system. We analyze also the variance of the process and compare its behavior for regions of asymptotic stationarity and nonstationarity, as well as for diffusive and noise-induced-resonance borderlines between them.

Effects of social distancing and isolation on epidemic spreading modeled via dynamical density functional theory

For preventing the spread of epidemics such as the coronavirus disease COVID-19, social distancing and the isolation of infected persons are crucial. However, existing reaction-diffusion equations for epidemic spreading are incapable of describing these effects. In this work, we present an extended model for disease spread based on combining a susceptible-infected-recovered model with a dynamical density functional theory where social distancing and isolation of infected persons are explicitly taken into account. We show that the model exhibits interesting transient phase separation associated with a reduction of the number of infections, and allows for new insights into the control of pandemics.