Group Contribution Coarse-Grained Molecular Simulations of Polystyrene Melts and Polystyrene Solutions in Alkanes Using the SAFT-γ Force Field

Simulation and Data-Driven Modeling of the Transport Properties of the Mie Fluid

This work reports the computation and modeling of the self-diffusivity (D*), shear viscosity (η*), and thermal conductivity (κ*) of the Mie fluid. The transport properties were computed using equilibrium molecular dynamics simulations for the Mie fluid with repulsive exponents (λr) ranging from 7 to 34 and at a fixed attractive exponent (λa) of 6 over the whole fluid density (ρ*) range and over a wide temperature (T*) range. The computed database consists of 17,212, 14,288, and 13,099 data points for self-diffusivity, shear viscosity, and thermal conductivity, respectively. The database is successfully validated against published simulation data. The above-mentioned transport properties are correlated using artificial neural networks (ANNs). Two modeling approaches were tested: a semiempirical formulation based on entropy scaling and an empirical formulation based on density and temperature as input variables. For the former, it was found that a unique formulation based on entropy scaling does not yield satisfactory results over the entire density range due to a divergent and incorrect scaling of the transport properties at low densities. For the latter empirical modeling approach, it was found that regularizing the data, e.g., modeling ρ*D* instead of D*, ln η* instead of η*, and ln κ* instead of κ*, as well as using the inverse of the temperature as an input feature, helps to ease the interpolation efforts of the artificial neural networks. The trained ANNs can model seen and unseen data over a wide range of density and temperature. Ultimately, the ANNs can be used alongside equations of state to regress effective force field parameters from volumetric and transport data.

Temperature-Transferable Coarse-Grained Models for Volumetric Properties of Poly(lactic Acid)

A new coarse-grained (CG) model, for which each monomer is mapped as one bead at its center of mass, was developed for simulating the volumetric properties of the polylactide (PLA) bulk. The three bonded CG potentials are first parametrized against the strain energies of the dimer, trimer, and tetramer, and the nonbonded CG potentials are then optimized to match the melt densities of the decamer. With the derived CG potentials, molecular dynamics (MD) simulations are found to reproduce thermal expansion and glass transition. By rescaling the dihedral and nonbonded potentials with temperature-independent factors, the glass transition temperature (Tg) is also satisfactorily restored with little modifications on the volumetric expansive coefficients at both rubbery and glassy states. Therefore, the finally optimized CG potentials exhibit excellent temperature transferability, as rationalized by the Simha-Boyer relation. Furthermore, it is confirmed that the dihedral torsions and nonbonded interactions play key roles in glass transition. Also, the simulated bulk moduli and conformational properties in a wide temperature range compare well with the referenced data. The proposed multiscale scheme has great potential in simulating thermo-mechanical properties of PLA and other polymers.

Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework

The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is utilized to study the effects of temperature on polymer brushes. The brushes are represented by a discrete energy spectrum, and energy degeneracies obtained through the replica-exchange Wang-Landau algorithm. The SEAQT equation of motion is applied to the density of states to establish a unique kinetic path from an initial thermodynamic state to a stable equilibrium state. The kinetic path describes the brush's evolution in state space, as it interacts with a thermal reservoir. The predicted occupation probabilities along the kinetic path are used to determine the expected thermodynamic and structural properties. The polymer density profile of a polystyrene brush in cyclohexane solvent is predicted using the equation of motion, and it agrees qualitatively with the experimental density profiles. The Flory-Huggins parameter chosen to describe brush-solvent interactions affects the solvent distribution in the brush but has a minimal impact on the polymer density profile. Three types of nonequilibrium kinetic paths with differing amounts of entropy production are considered: a heating path, a cooling path, and a heating-cooling path. Properties such as tortuosity, radius of gyration, brush density, solvent density, and brush chain conformations are calculated for each path.

Molecular modelling of the thermophysical properties of fluids: expectations, limitations, gaps and opportunities

This manuscript provides an overview of the current state of the art in terms of the molecular modelling of the thermophysical properties of fluids. It is intended to manage the expectations and serve as guidance to practising physical chemists, chemical physicists and engineers in terms of the scope and accuracy of the more commonly available intermolecular potentials along with the peculiarities of the software and methods employed in molecular simulations while providing insights on the gaps and opportunities available in this field. The discussion is focused around case studies which showcase both the precision and the limitations of frequently used workflows.

Continuum-scale modelling of polymer blends using the Cahn-Hilliard equation: transport and thermodynamics

The Cahn-Hilliard equation is commonly used to study multi-component soft systems such as polymer blends at continuum scales. We first systematically explore various features of the equation system, which give rise to a deep connection between transport and thermodynamics-specifically that the Gibbs free energy of mixing function is central to formulating a well-posed model. Accordingly, we explore how thermodynamic models from three broad classes of approach (lattice-based, activity-based and perturbation methods) can be incorporated within the Cahn-Hilliard equation and examine how they impact the numerical solution for two model polymer blends, noting that although the analysis presented here is focused on binary mixtures, it is readily extensible to multi-component mixtures. It is observed that, although the predicted liquid-liquid interfacial tension is quite strongly affected, the choice of thermodynamic model has little influence on the development of the morphology.

Oligomer/Polymer Blend Phase Diagram and Surface Concentration Profiles for Squalane/Polybutadiene: Experimental Measurements and Predictions from SAFT-γ Mie and Molecular Dynamics Simulations

The compatibility and surface behavior of squalane-polybutadiene mixtures are studied by experimental cloud point and neutron reflectivity measurements, statistical associating fluid theory (SAFT), and molecular dynamics (MD) simulations. A SAFT-γ Mie model is shown to be successful in capturing the cloud point curves of squalane-polybutadiene and squalane-cis-polybutadiene binary mixtures, and the same SAFT-γ Mie model is used to develop a thermodynamically consistent top-down coarse-grained force field to describe squalane-polybutadiene. Coarse-grained molecular dynamics simulations are performed to study surface behavior for different concentrations of squalane, with the system exhibiting surface enrichment and a wetting transition. Simulated surface profiles are compared with those obtained by fitting to neutron reflectivity data obtained from thin films composed of deuterated squalane (d-sq)-polybutadiene. The presented top-down parametrization methodology is a fast and thermodynamically reliable approach for predicting properties of oligomer-polymer mixtures, which can be challenging for either theory or MD simulations alone.

Molecular Dynamics Simulation of the Superspreading of Surfactant-Laden Droplets. A Review

Superspreading is the rapid and complete spreading of surfactant-laden droplets on hydrophobic substrates. This phenomenon has been studied for many decades by experiment, theory, and simulation, but it has been only recently that molecular-level simulation has provided significant insights into the underlying mechanisms of superspreading thanks to the development of accurate force-fields and the increase of computational capabilities. Here, we review the main advances in this area that have surfaced from Molecular Dynamics simulation of all-atom and coarse-grained models highlighting and contrasting the main results and discussing various elements of the proposed mechanisms for superspreading. We anticipate that this review will stimulate further research on the interpretation of experimental results and the design of surfactants for applications requiring efficient spreading, such as coating technology.

Aggregation Behavior of Model Asphaltenes Revealed from Large-Scale Coarse-Grained Molecular Simulations

Fully atomistic simulations of models of asphaltenes in simple solvents have allowed the study of trends in aggregation phenomena to understand the underlying role played by molecular structure. The detail included at this scale of molecular modeling is, however, at odds with the required spatial and temporal resolution needed to fully understand asphaltene aggregation. The computational cost required to explore the relevant scales can be reduced by employing coarse-grained (CG) models, which consist of lumping a few atoms into a single segment that is characterized by effective interactions. In this work, CG force fields developed via the statistical associating fluid theory (SAFT-γ) [ Müller , E. A. ; Jackson , G. Annu. Rev. Chem. Biomol. Eng. 5 , 2014 , 405 - 427 ] equation of state (EoS) provide a reliable pathway to link the molecular description with macroscopic thermophysical data. A recent modification of the SAFT-VR EoS [ Müller , E. A. ; Mejía , A. Langmuir 33 , 2017 , 11518 - 11529 ], which allows for the parameterization of homonuclear rings, is selected as the starting point to develop CG models for polycyclic aromatic hydrocarbons. The new aromatic-core models, along with others published for simpler organic molecules, are adopted for the construction of asphaltene models by combining different chemical moieties in a group-contribution fashion. We apply the procedure to two previously reported asphaltene models and perform molecular dynamics simulations to validate the coarse-grained representation against benchmark systems of 27 asphaltenes in a pure solvent (toluene or heptane) described in a fully atomistic fashion. An excellent match between both levels of description is observed for the cluster size, radii of gyration, and relative-shape-anisotropy-factor distributions. We exploit the advantages of the CG representation by simulating systems containing up to 2000 asphaltene molecules in an explicit solvent investigating the effect of asphaltene concentration, solvent composition, and temperature on aggregation. By studying large systems facilitated by the use of CG models, we observe stable continuous distributions of molecular aggregates at conditions away from the two-phase precipitation point. As a further example application, a widely accepted interpretation of cluster-size distributions in asphaltenic systems is challenged by performing system-size tests, reversibility checks, and a time-dependence analysis. The proposed coarse-graining procedure is seen to be general and predictive and, hence, can be applied to other asphaltenic molecular structures.

Fluid-solid phase transition of n-alkane mixtures: Coarse-grained molecular dynamics simulations and diffusion-ordered spectroscopy nuclear magnetic resonance

Wax appearance temperature (WAT), defined as the temperature at which the first solid paraffin crystal appears in a crude oil, is one of the key flow assurance indicators in the oil industry. Although there are several commonly-used experimental techniques to determine WAT, none provides unambiguous molecular-level information to characterize the phase transition between the homogeneous fluid and the underlying solid phase. Molecular Dynamics (MD) simulations employing the statistical associating fluid theory (SAFT) force field are used to interrogate the incipient solidification states of models for long-chain alkanes cooled from a melt to an arrested state. We monitor the phase change of pure long chain n-alkanes: tetracosane (C₂₄H₅₀) and triacontane (C₃₀H₆₂), and an 8-component surrogate n-alkane mixture (C₁₂-C₃₃) built upon the compositional information of a waxy crude. Comparison to Diffusion Ordered Spectroscopy Nuclear Magnetic Resonance (DOSY NMR) results allows the assessment of the limitations of the coarse-grained models proposed. We show that upon approach to freezing, the heavier components restrict their motion first while the lighter ones retain their mobility and help fluidize the mixture. We further demonstrate that upon sub-cooling of long n-alkane fluids and mixtures, a discontinuity arises in the slope of the self-diffusion coefficient with decreasing temperature, which can be employed as a marker for the appearance of an arrested state commensurate with conventional WAT measurements.