Transport coefficients from Einstein-Helfand relations using standard and energy-conserving dissipative particle dynamics methods.

Abstract

In this article we demonstrate that contrary to general belief, the standard Einstein-Helfand (EH) formulas are valid for the evaluation of transport coefficients of systems containing dissipative and random forces provided that for these mesoscopic systems: (i) the corresponding conservation laws are satisfied, and (ii) the transition probabilities satisfy detailed balance. Dissipative particle dynamics (DPD) and energy-conserving DPD methods (DPDE), for instance, are archetypical of such mesoscopic approaches satisfying these properties. To verify this statement, we have derived a mesoscopic heat flux form for the DPDE method, suitable for the calculation of the thermal conductivity from an EH expression. We have compared EH measurements against non-equilibrium simulation values for different scenarios, including many-body potentials, and have found excellent agreement in all cases. The expressions are valid notably for systems with density- and temperature-dependent potentials, such as the recently developed generalised DPDE method (GenDPDE) [Avalos et al., Phys. Chem. Chem. Phys., 2019, 21, 24891]. We thus demonstrate that traditional EH formulas in equilibrium simulations can be widely used to obtain transport coefficients, provided that the appropriate expression for the associated flux is used.