Local density dependent potential for compressible mesoparticles

Ultra-coarse-graining of homopolymers in inhomogeneous systems

We develop coarse-grained (CG) models for simulating homopolymers in inhomogeneous systems, focusing on polymer films and droplets. If the CG polymers interact solely through two-body potentials, then the films and droplets either dissolve or collapse into small aggregates, depending on whether the effective polymer-polymer interactions have been determined from reference simulations in the bulk or at infinite dilution. To address this shortcoming, we include higher order interactions either through an additional three-body potential or a local density-dependent potential (LDP). We parameterize the two- and three-body potentials via force matching, and the LDP through relative entropy minimization. While the CG models with three-body interactions fail at reproducing stable polymer films and droplets, CG simulations with an LDP are able to do so. Minor quantitative differences between the reference and the CG simulations, namely a slight broadening of interfaces accompanied by a smaller surface tension in the CG simulations, can be attributed to the deformation of polymers near the interfaces, which cannot be resolved in the CG representation, where the polymers are mapped to spherical beads.

Coarse-graining of polyisoprene melts using inverse Monte Carlo and local density potentials

Bottom-up coarse-graining of polymers is commonly performed by matching structural order parameters such as distribution of bond lengths, bending and dihedral angles, and pair distribution functions. In this study, we introduce the distribution of nearest-neighbors as an additional order parameter in the concept of local density potentials. We describe how the inverse-Monte Carlo method provides a framework for forcefield development that is capable of overcoming challenges associated with the parameterization of interaction terms in polymer systems. The technique is applied on polyisoprene melts as a prototype system. We demonstrate that while different forcefields can be developed that perform equally in terms of matching target distributions, the inclusion of nearest-neighbors provides a straightforward route to match both thermodynamic and conformational properties. We find that several temperature state points can also be addressed, provided that the forcefield is refined accordingly. Finally, we examine both the single-particle and the collective dynamics of the coarse-grain models, demonstrating that all forcefields present a similar acceleration relative to the atomistic systems.

Generalised dissipative particle dynamics with energy conservation: density- and temperature-dependent potentials

We present a generalised, energy-conserving dissipative particle dynamics (DPDE) method appropriate for the non-isothermal simulation of particle interaction force fields that are both density- and temperature-dependent. A detailed derivation is formulated in a bottom-up manner by considering the thermodynamics of small systems with the appropriate consideration of the fluctuations. Connected to the local volume is a local density and corresponding local pressure, which is determined from an equation-of-state based force field that depends also on a particle temperature. Compared to the original DPDE method, the formulation of the generalised DPDE method requires a change in the independent variable from the particle internal energy to the particle entropy. As part of the re-formulation, the terms dressed particle entropy and the corresponding dressed particle temperature are introduced, which depict the many-body contributions in the local volume. The generalised DPDE method has similarities to the energy form of the smoothed dissipative particle dynamics method, yet fundamental differences exist, which are described in the manuscript. The basic dynamic equations are presented along with practical considerations for implementing the generalised DPDE method, including a numerical integration scheme based on the Shardlow-like splitting algorithm. Demonstrations and validation tests are performed using analytical equation-of-states for the van der Waals and Lennard-Jones fluids. Particle probability distributions are analysed, where excellent agreement with theoretical estimates is demonstrated. As further validation of the generalised DPDE method, both equilibrium and non-equilibrium simulation scenarios are considered, including adiabatic flash heating response and vapour-liquid phase separation.