Microscopic formula for transport coefficients of causal hydrodynamics

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.

Microscopic derivation of coarse-grained, energy-conserving generalized Langevin dynamics

Properly simulating nonequilibrium phenomena such as thermal transport and shock wave propagation in complex condensed matter systems require the conservation of system's internal energy. This precludes the application of the coarse-grained (CG) generalized Langevin equation (GLE) dynamics due to the presence of dissipative interactions. Attempts to address this issue have been pursued both phenomenologically and from entropy-based first principles for dissipative particle dynamics (DPD, a Markovian variant of the CG GLE dynamics) by introducing an energy conserving extension of DPD (DPD-E). We present here a rigorous microscopic derivation of two energy conserving variants of the CG GLE dynamics by extending the CG equations of motion to include the GLE for certain internal energy observables of the microscopic system. We consider two choices of such observables: the total internal energy and a set of internal energies of the CG particles. The derivation is performed using the Mori-Zwanzig projection operator method in the Heisenberg picture for time evolution of thermodynamic expectations and the recently introduced interpretation of the Zwanzig projection operator [S. Izvekov, J. Chem. Phys. 146(12), 124109 (2017)] which allows an exact calculation of the memory and projected terms. We begin with equilibrium conditions and show that the GLE dynamics for the internal energy observables is purely dissipative. Our extension of the GLE dynamics to quasiequilibrium conditions (necessary to observe heat transport) is based on the generalized canonical ensemble approach and transport equation using the nonequilibrium statistical operator (NSO) method. We derive closed microscopic expressions for conductive heat transfer coefficients in the limit of neglecting dissipation in heat transfer and in the lowest order of deviation from equilibrium. After employing the Markov approximation, we compare the equations of motion to the published DPD-E equations. Our equations contain additional energy transfer terms not reported in the previous works. Additionally, we show that, despite neglecting dissipative processes in heat transport, the heat transfer coefficients and random force are related in a way reminiscent of the fluctuation-dissipation relation. The formalism presented here is sufficiently general for the rigorous formulation of the GLE dynamics for arbitrary microscopic phase space observables as well as sampling different microscopic ensembles in CG simulations.

Mori-Zwanzig theory for dissipative forces in coarse-grained dynamics in the Markov limit

We derive alternative Markov approximations for the projected (stochastic) force and memory function in the coarse-grained (CG) generalized Langevin equation, which describes the time evolution of the center-of-mass coordinates of clusters of particles in the microscopic ensemble. This is done with the aid of the Mori-Zwanzig projection operator method based on the recently introduced projection operator [S. Izvekov, J. Chem. Phys. 138, 134106 (2013)10.1063/1.4795091]. The derivation exploits the "generalized additive fluctuating force" representation to which the projected force reduces in the adopted projection operator formalism. For the projected force, we present a first-order time expansion which correctly extends the static fluctuating force ansatz with the terms necessary to maintain the required orthogonality of the projected dynamics in the Markov limit to the space of CG phase variables. The approximant of the memory function correctly accounts for the momentum dependence in the lowest (second) order and indicates that such a dependence may be important in the CG dynamics approaching the Markov limit. In the case of CG dynamics with a weak dependence of the memory effects on the particle momenta, the expression for the memory function presented in this work is applicable to non-Markov systems. The approximations are formulated in a propagator-free form allowing their efficient evaluation from the microscopic data sampled by standard molecular dynamics simulations. A numerical application is presented for a molecular liquid (nitromethane). With our formalism we do not observe the "plateau-value problem" if the friction tensors for dissipative particle dynamics (DPD) are computed using the Green-Kubo relation. Our formalism provides a consistent bottom-up route for hierarchical parametrization of DPD models from atomistic simulations.

Shear viscosity coefficient and relaxation time of causal dissipative hydrodynamics in QCD

The shear viscosity coefficient and the corresponding relaxation time for causal dissipative hydrodynamics are calculated based on the microscopic formula proposed in T. Koide and T. Kodama [Phys. Rev. E 78, 051107 (2008)10.1103/PhysRevE.78.051107]. Here, the exact formula is transformed into a more compact form and applied to evaluate these transport coefficients in the chiral perturbation theory and perturbative QCD. It is shown that in the leading order calculation, the causal shear viscosity coefficient eta reduces to that of the ordinary Green-Kubo-Nakano formula, and the relaxation time tau(pi) is related to eta and pressure P by a simple relationship, tau(pi)=eta/P.

Transport coefficients of non-Newtonian fluid and causal dissipative hydrodynamics

A formula to calculate the transport coefficients of the causal dissipative hydrodynamics is derived by using the projection operator method (Mori-Zwanzig formalism) [T. Koide, Phys. Rev. E 75, 060103(R) (2007)]. This is an extension of the Green-Kubo-Nakano (GKN) formula to the case of non-Newtonian fluids, which is the essential factor to preserve the relativistic causality in relativistic dissipative hydrodynamics. This formula is the generalization of the GKN formula in the sense that it can reproduce the GKN formula in a certain limit. In this work, we extend the previous work so as to apply to more general situations.