Dissipative particle dynamics with reactions: Application to RDX decomposition

A coarse-grain reactive model of RDX: Molecular resolution at the μm scale

Predictive models for the thermal, chemical, and mechanical response of high explosives at extreme conditions are important for investigating their performance and safety. We introduce a particle-based, reactive model of 1,3,5-trinitro-1,3,5-triazinane (RDX) with molecular resolution utilizing generalized energy-conserving dissipative particle dynamics with reactions. The model is parameterized with respect to the data from atomistic molecular dynamics simulations as well as from quantum mechanical calculations, thus bridging atomic processes to the mesoscales, including microstructures and defects. It accurately captures the response of RDX under a range of thermal loading conditions compared to atomistic simulations. In addition, the Hugoniot response of the CG model in the overdriven regime reasonably matches atomistic simulations and experiments. Exploiting the model's high computational efficiency, we investigate mesoscale systems involving millions of molecules and characterize size-dependent criticality of hotspots in RDX. The combination of accuracy and computational efficiency of our reactive model provides a tool for investigation of mesoscale phenomena, such as the role of microstructures and defects in the shock-to-deflagration transition, through particle-based simulation.

Generalized Energy-Conserving Dissipative Particle Dynamics with Mass Transfer. Part 1: Theoretical Foundation and Algorithm

An extension of the generalized energy-conserving dissipative particle dynamics method (GenDPDE) that allows mass transfer between mesoparticles via a diffusion process is presented. By considering the concept of the mesoparticles as property carriers, the complexity and flexibility of the GenDPDE framework were enhanced to allow for interparticle mass transfer under isoenergetic conditions, notated here as GenDPDE-M. In the formulation, diffusion is described via the theory of mesoscale irreversible processes based on linear relationships between the fluxes and thermodynamic forces, where their fluctuations are described by Langevin-like equations. The mass exchange between mesoparticles is such that the mass of the mesoparticle remains unchanged after the transfer process and requires additional considerations regarding the coupling with other system properties such as the particle internal energy. The proof-of-concept work presented in this article is the first part of a two-part article series. In Part 1, the development of the GenDPDE-M theoretical framework and the derivation of the algorithm are presented in detail. Part 2 of this article series is targeted for practitioners, where applications, demonstrations, and practical considerations for implementing the GenDPDE-M method are presented and discussed.

Generalized Energy-Conserving Dissipative Particle Dynamics with Mass Transfer. Part 2: Applications and Demonstrations

We present the second part of a two-part paper series intended to address a gap in computational capability for coarse-grain particle modeling and simulation, namely, the simulation of phenomena in which diffusion via mass transfer is a contributing mechanism. In part 1, we presented a formulation of a dissipative particle dynamics method to simulate interparticle mass transfer, termed generalized energy-conserving dissipative particle dynamics with mass transfer (GenDPDE-M). In the GenDPDE-M method, the mass of each mesoparticle remains constant following the interparticle mass exchange. In part 2 of this series, further verification and demonstrations of the GenDPDE-M method are presented for mesoparticles with embedded binary mixtures using the ideal gas (IG) and van der Waals (vdW) equation-of-state (EoS). The targeted readership of part 2 is toward practitioners, where applications and practical considerations for implementing the GenDPDE-M method are presented and discussed, including a numerical discretisztion algorithm for the equations-of-motion. The GenDPDE-M method is verified by reproducing the particle distributions predicted by Monte Carlo simulations for the IG and vdW fluids, along with several demonstrations under both equilibrium and non-equilibrium conditions. GenDPDE-M can be generally applied to multi-component mixtures and to other fundamental EoS, such as the Lennard-Jones or Exponential-6 models, as well as to more advanced EoS models such as Statistical Associating Fluid Theory.

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Generalized Energy-Conserving Dissipative Particle Dynamics with Reactions

We present an extension of the generalized energy-conserving dissipative particle dynamics method (J. Bonet Avalos, et al., Phys Chem Chem Phys, 2019, 21, 24891-24911) to include chemical reactivity, denoted GenDPDE-RX. GenDPDE-RX provides a means of simulating chemical reactivity at the micro- and mesoscales, while exploiting the attributes of density- and temperature-dependent many-body force fields, which include improved transferability and scalability compared to two-body pairwise models. The GenDPDE-RX formulation considers intra-particle reactivity via a coarse-grain reactor construct. Extent-of-reaction variables assigned to each coarse-grain particle monitor the temporal evolution of the prescribed reaction mechanisms and kinetics assumed to occur within the particle. Descriptions of the algorithm, equations of motion, and numerical discretization are presented, followed by verification of the GenDPDE-RX method through comparison with reaction kinetics theoretical model predictions. Demonstrations of the GenDPDE-RX method are performed using constant-volume adiabatic heating simulations of three different reaction models, including both reversible and irreversible reactions, as well as multistep reaction mechanisms. The selection of the demonstrations is intended to illustrate the flexibility and generality of the method but is inspired by real material systems that span from fluids to solids. Many-body force fields using analytical forms of the ideal gas, Lennard-Jones, and exponential-6 equations of state are used for demonstration, although application to other forms of equation of states is possible. Finally, the flexibility of the GenDPDE-RX framework is addressed with a brief discussion of other possible adaptations and extensions of the method.

Interactions of Crosslinked Polyacrylic Acid Polyelectrolyte Gels with Nonionic and Ionic Surfactants

The morphology and stability of surfactant-loaded polyelectrolyte gels are of great interest for a variety of personal care, cosmetic, and pharmaceutical products. However, the mechanisms of surfactant interactions with gel-forming polymers are poorly understood and experimentally challenging. The aim of this work is to explore in silico the specifics of surfactant absorption within polyelectrolyte gels drawing on the examples of typical non-ionic octaethylene glycol monooctyl ether (C₈E₈) and anionic sodium dodecyl sulfate (SDS) surfactants and polyacrylic acid modified with hydrophobic sidechains mimicking the practically important Carbopol polymer. Using the systematically parameterized coarse-grained dissipative particle dynamics models, we generate and characterize the equilibrium conformations and swelling of the polymer films in aqueous solutions with the surfactant concentrations varied up to the critical micelle concentration (cmc). We discover the striking difference in interactions of Carbopol-like polymers with nonionic and ionic surfactants under mildly acidic conditions. The sorption of C₈E₈ within the polymer film is found substantial. As the surfactant concentration increases, the polymer film swells and, close to cmc, becomes unstable due to the formation and growth of water pockets filled with surfactant micelles. Sorption of SDS at the same bulk concentrations is found much lower, with only about 1% of surfactant mass fraction achieved at cmc. As the SDS concentration increases further, a lamellae structure is formed within the film, which remains stable. Reduced swelling and higher stability indicate better prospects of using SDS-type surfactants with Carbopol-based gels in formulations for detergents and personal care products.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Generalized energy-conserving dissipative particle dynamics revisited: Insight from the thermodynamics of the mesoparticle leading to an alternative heat flow model

Recently we introduced the generalized energy-conserving dissipative particle dynamics method (GenDPDE) [J. Bonet Avalos, M. Lísal, J. P. Larentzos, A. D. Mackie, and J. K. Brennan, Phys. Chem. Chem. Phys. 21, 24891 (2019)]PPCPFQ1463-907610.1039/C9CP04404C, which has been formulated for an emerging class of density- and temperature-dependent coarse-grain models. In the original work, GenDPDE was formulated to ensure a fundamental link is maintained with the underlying physical system at the higher resolution scale. In this paper, we revisit the formulation of the GenDPDE method, and rederive the particle thermodynamics to ensure consistency at the opposing scale extreme, i.e., between the local thermodynamics in the mesoscopic systems and the corresponding macroscopic properties. We demonstrate this consistency by introducing unambiguous, physically meaningful definitions of the heat and work, which lead to the formulation of an alternative heat flow model that is analogous to Fourier's law of heat conduction. We present further analysis of the internal, unresolved degrees-of-freedom of the mesoparticles by considering the thermodynamics of an individual mesoparticle within the GenDPDE framework. Several key outcomes of the analysis include: (i) demonstration that the choice of the independent variables alters the particle thermodynamic description; (ii) demonstration that the mesoscopic thermodynamic transformations introduce additional terms of the order of the size of the local fluctuations, which prevent an unambiguous definition of both the heat and work; (iii) an emphasis on the importance of the choice of the proper estimators of the thermodynamic properties that are embedded in the chosen thermodynamic description; and (iv) a clearly defined path for determining any thermodynamic quantity dressed by the fluctuations. The further insight provided by this deeper analysis is useful for both readers interested in the GenDPDE theoretical framework, as well as readers interested in the practical ramifications of the analysis, namely, the alternative heat flow model.

Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties

Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.