We present a corresponding states correlation based on the description of fluid phase properties by means of an Mie intermolecular potential applied to tangentially bonded spheres. The macroscopic properties of this model fluid are very accurately represented by a recently proposed version of the Statistical Associating Fluid Theory (the SAFT-γ version). The Mie potential can be expressed in a conformal manner in terms of three parameters that relate to a length scale, σ, an energy scale, ε, and the range or functional form of the potential, λ, while the nonsphericity or elongation of a molecule can be appropriately described by the chain length, m. For a given chain length, correlations are given to scale the SAFT equation of state in terms of three experimental parameters: the acentric factor, the critical temperature, and the saturated liquid density at a reduced temperature of 0.7. The molecular nature of the equation of state is exploited to make a direct link between the macroscopic thermodynamic parameters used to characterize the equation of state and the parameters of the underlying Mie potential. This direct link between macroscopic properties and molecular parameters is the basis to propose a top-down method to parametrize force fields that can be used in the coarse grained molecular modeling (Monte Carlo or molecular dynamics) of fluids. The resulting correlation is of quantitative accuracy and examples of the prediction of simulations of vapor−liquid equilibria and interfacial tensions are given. In essence, we present a recipe that allows one to obtain intermolecular potentials for use in the molecular simulation of simple and chain fluids from macroscopic experimentally determined constants, namely the acentric factor, the critical temperature, and a value of the liquid density at a reduced temperature of 0.7.

A c c e p t e d M a n u s c r i p t Highlights -A framework to obtain accurate coarse-grained force fields is described -Force field parameters are presented for compounds of interest -A methodology to simulate bubble points of complex mixtures is proposed. -The force fields are shown to be robust and transferrable Highlights (for review) AbstractThe atomistically-detailed molecular modelling of petroleum fluids is challenging, amongst other aspects, due to the very diverse multicomponent and asymmetric nature of the mixtures in question. Complicating matters further, the time scales for many important processes can be much larger that the current and foreseeable capacity of modern computers running fully-atomistic models. To overcome these limitations, a coarse grained (CG) model is proposed where some of the less-important degrees of freedom are safely integrated out, leaving as key parameters the average energy levels, the molecular conformations and the range of the Mie intermolecular potentials employed as the basis of the model. The parametrization is performed by using an analytical equation of state of the statistical associating fluid theory (SAFT) family to link the potential parameters to macroscopically observed thermophysical properties. The parameters found through this top-down approach are used directly in molecular dynamics simulations of multi-component multi-phase systems. The procedure is exemplified by calculating the phase envelope of the methane-decane binary and of two synthetic light condensate mixtures. A procedure based on a discrete expansion of a mixture is used to determine the bubble points of these latter mixtures, with an excellent agreement to experimental data. The model presented is entirely predictive and an abridged table of parameters for some fluids of interest is provided.

A molecular model for simulating the aggregation of asphaltenes and resins in crude oils on a mesoscale is proposed. The asphaltene molecules are treated as discotic seven-center Lennard-Jonesium molecules, the resins are modeled as single spheres, and the surrounding crude oil is modeled as a continuum, characterized by a screening factor, and defined using a combination of its Hamaker and dielectric constants. The parameters for the model are obtained by coarse-graining the potential energy surface obtained from model atomistic simulations of pairs of asphaltenes and resins. Canonical Monte Carlo simulations are performed with this model, and effects of temperature, asphaltene, and resin concentration are studied parametrically. The results agree with experimentally observed tendencies. The asphaltene is seen not to conform to a linear aggregation model, but exhibits a more complex multimodal aggregation pattern. The screening constant of the crude oil, which ultimately controls the aggregation, can itself be related to other measurable quantities such as the refractive index.

Abstract:A methodology for the determination of the solid-fluid contact angle, to be employed within molecular dynamics (MD) simulations, is developed and systematically applied. The calculation of the contact angle of a fluid drop on a given surface, averaged over an equilibrated MD trajectory, is divided in three main steps: (i) the determination of the fluid molecules that constitute the interface, (ii) the treatment of the interfacial molecules as a point cloud data set to define a geometric surface, using surface meshing techniques to compute the surface normals from the mesh, (iii) the collection and averaging of the interface normals collected from the post-processing of the MD trajectory. The average vector thus found is used to calculate the Cassie contact angle (i.e., the arccosine of the averaged normal z-component). As an example we explore the effect of the size of a drop of water on the observed solid-fluid contact angle. A single coarsegrained bead representing two water molecules and parameterized using the SAFT-γ Mie equation of state (EoS) is employed, meanwhile the solid surfaces are mimicked using integrated potentials. The contact angle is seen to be a strong function of the system size for small nano-droplets. The thermodynamic limit, corresponding to the infinite size (macroscopic) drop is only truly recovered when using an excess of half a million water coarse-grained beads and/or a drop radius of over 26 nm.