SAFT-γ force field for the simulation of molecular fluids: 4. A single-site coarse-grained model of water applicable over a wide temperature range

In this work, we develop coarse-grained (CG) force fields for water, where the effective CG intermolecular interactions between particles are estimated from an accurate description of the macroscopic experimental vapour-liquid equilibria data by means of a molecular-based equation of state. The statistical associating fluid theory for Mie (generalised Lennard-Jones) potentials of variable range (SAFT-VR Mie) is used to parameterise spherically symmetrical (isotropic) force fields for water. The resulting SAFT-γ CG models are based on the Mie (8-6) form with size and energy parameters that are temperature dependent; the latter dependence is a consequence of the angle averaging of the directional polar interactions present in water. At the simplest level of CG where a water molecule is represented as a single bead, it is well known that an isotropic potential cannot be used to accurately reproduce all of the thermodynamic properties of water simultaneously. In order to address this deficiency, we propose two CG potential models of water based on a faithful description of different target properties over a wide range of temperatures: our CGW1-vle model is parameterised to match the saturated-liquid density and vapour pressure; our other CGW1-ift model is parameterised to match the saturated-liquid density and vapour-liquid interfacial tension. A higher level of CG corresponding to two water molecules per CG bead is also considered: the corresponding CGW2-bio model is developed to reproduce the saturated-liquid density and vapour-liquid interfacial tension in the physiological temperature range, and is particularly suitable for the large-scale simulation of bio-molecular systems. A critical comparison of the phase equilibrium and transport properties of the proposed force fields is made with the more traditional atomistic models

raaSAFT: A framework enabling coarse-grained molecular dynamics simulations based on the SAFT-γ Mie force field

We describe here raaSAFT, a Python code that enables the setup and running of coarse-grained molecular dynamics simulations in a systematic and efficient manner. The code is built on top of the popular HOOMD-blue code, and as such harnesses the computational power of GPUs. The methodology makes use of the SAFT-γ Mie force field, so the resulting coarse grained pair potentials are both closely linked to and consistent with the macroscopic thermodynamic properties of the simulated fluid. In raaSAFT both homonuclear and heteronuclear models are implemented for a wide range of compounds spanning from linear alkanes, to more complicated fluids such as water and alcohols, all the way up to nonionic surfactants and models of asphaltenes and resins. Adding new compounds as well as new features is made straightforward by the modularity of the code. To demonstrate the ease-of-use of raaSAFT, we give a detailed walkthrough of how to simulate liquid–liquid equilibrium of a hydrocarbon with water. We describe in detail how both homonuclear and heteronuclear compounds are implemented. To demonstrate the performance and versatility of raaSAFT, we simulate a large polymer-solvent mixture with 300 polystyrene molecules dissolved in 42 700 molecules of heptane, reproducing the experimentally observed temperature-dependent solubility of polystyrene. For this case we obtain a speedup of more than three orders of magnitude as compared to atomistically-detailed simulations

Modelling the superspreading of surfactant-laden droplets with computer simulation

The surfactant-driven superspreading of droplets on hydrophobic substrates is considered. A key element of the superspreading mechanism is the adsorption of surfactant molecules from the liquid-vapour interface onto the substrate through the contact line, which must be coordinated with the replenishment of interfaces with surfactant from the interior of the droplet. We use molecular dynamics simulations with coarse-grained force fields to provide a detailed structural description of the droplet shape and surfactant dynamics during the superspreading process. We also provide a simple method for accurate estimation of the contact angle subtended by the droplets at the contact line

In Silico Determination of Gas Permeabilities by Non-Equilibrium Molecular Dynamics: CO2 and He through PIM-1.

© 2015 by the authors; licensee MDPI, Basel, Switzerland.We study the permeation dynamics of helium and carbon dioxide through an atomistically detailed model of a polymer of intrinsic microporosity, PIM-1, via non-equilibrium molecular dynamics (NEMD) simulations. This work presents the first explicit molecular modeling of gas permeation through a high free-volume polymer sample, and it demonstrates how permeability and solubility can be obtained coherently from a single simulation. Solubilities in particular can be obtained to a very high degree of confidence and within experimental inaccuracies. Furthermore, the simulations make it possible to obtain very specific information on the diffusion dynamics of penetrant molecules and yield detailed maps of gas occupancy, which are akin to a digital tomographic scan of the polymer network. In addition to determining permeability and solubility directly from NEMD simulations, the results shed light on the permeation mechanism of the penetrant gases, suggesting that the relative openness of the microporous topology promotes the anomalous diffusion of penetrant gases, which entails a deviation from the pore hopping mechanism usually observed in gas diffusion in polymers

Defect-Dependent Corrugation in Graphene

Graphene’s intrinsically corrugated and wrinkled topology fundamentally influences its electronic, mechanical, and chemical properties. Experimental techniques allow the manipulation of pristine graphene and the controlled production of defects which allows one to control the atomic out-of-plane fluctuations and thus tune graphene’s properties. Here, we perform large scale machine learning-driven molecular dynamics simulations to understand the impact of defects on the structure of graphene. We find that defects cause significantly higher corrugation leading to a strongly wrinkled surface. The magnitude of this structural transformation strongly depends on the defect concentration and specific type of defect. Analyzing the atomic neighborhood of the defects reveals that the extent of these morphological changes depends on the preferred geometrical orientation and the interactions between defects. While our work highlights that defects can strongly affect graphene’s morphology, it also emphasizes the differences between distinct types by linking the global structure to the local environment of the defects

Solubilities of pyrene in organic solvents: Comparison between chemical potential calculations using a cavity-based method and direct coexistence simulations

In this paper, we benchmark a cavity-based simulation method for calculating the relative solubility of large molecules in explicit solvents. The essence of the procedure is the accounting of the Gibbs energy change associated with an alchemical thermodynamic cycle where, in sequence, a cavity is created in a solvent, a solute is inserted in the cavity and the cavity is annihilated. The Gibbs energy change is equated to the excess chemical potential allowing the comparison of solubilities in different solvents. The results obtained using the cavity-based method are compared to direct large-scale molecular dynamics simulations performed using coarse-grained models for calculating the partition coefficient of pyrene between heptane and toluene. We demonstrate the applicability of this cavity-based technique under high pressure/temperature conditions.The authors gratefully acknowledge the generous funding and technical support for this work from BP Plc through the International Centre for Advanced Materials (ICAM) which made this research possible

Liquid crystal phase behaviour of attractive disc-like particles.

We employ a generalized van der Waals-Onsager perturbation theory to construct a free energy functional capable of describing the thermodynamic properties and orientational order of the isotropic and nematic phases of attractive disc particles. The model mesogen is a hard (purely repulsive) cylindrical disc particle decorated with an anisotropic square-well attractive potential placed at the centre of mass. Even for isotropic attractive interactions, the resulting overall inter-particle potential is anisotropic, due to the orientation-dependent excluded volume of the underlying hard core. An algebraic equation of state for attractive disc particles is developed by adopting the Onsager trial function to characterize the orientational order in the nematic phase. The theory is then used to represent the fluid-phase behaviour (vapour-liquid, isotropic-nematic, and nematic-nematic) of the oblate attractive particles for varying values of the molecular aspect ratio and parameters of the attractive potential. When compared to the phase diagram of their athermal analogues, it is seen that the addition of an attractive interaction facilitates the formation of orientationally-ordered phases. Most interestingly, for certain aspect ratios, a coexistence between two anisotropic nematic phases is exhibited by the attractive disc-like fluids

Spreading of aqueous droplets with common and superspreading surfactants. A molecular dynamics study

arXiv:1909.00775 [physics.flu-dyn] preprint [v1] Mon, 2 Sep 2019 15:55:46 UTC (8,139 KB) submitted to Colloids and Surfaces A: Physicochemical and Engineering Aspects 581, 123810 (2019) DOI: https://doi.org/10.1016/j.colsurfa.2019.12381The surfactant-driven spreading of droplets is an essential process in many applications ranging from coating flow technology to enhanced oil recovery. Despite the significant advancement in describing spreading processes in surfactant-laden droplets, including the exciting phenomena of superspreading, many features of the underlying mechanisms require further understanding. Here, we have carried out molecular dynamics simulations of a coarse-grained model with a force-field obtained from the statistical associating fluid theory to study droplets laden with common and superspreading surfactants. We have confirmed the important elements of the superspreading mechanism, i.e. the adsorption of surfactant at the contact line and the fast replenishment of surfactant from the bulk. Through a detailed comparison of a range of droplets with different surfactants, our analysis has indicated that the ability of surfactant to adsorb at the interfaces is the key feature of the superspreading mechanism. To this end, surfactants that tend to form aggregates and have a strong hydrophobic attraction in the aggregated cores prevent the fast replenishment of the interfaces, resulting in reduced spreading ability. We also show that various surfactant properties, such as diffusion and architecture, play a secondary role in the spreading process. Moreover, we discuss various drop properties such as the height, contact angle, and surfactant surface concentration, highlighting differences between superspreading and common surfactants. We anticipate that our study will provide further insight for applications requiring the control of wetting.This project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 778104. E.A.M. acknowledges partial support from the Engineering and Physical Sciences Research Council (EPSRC) of the U.K. through grants EP/I018212, EP/J010502, and EP/R013152. This research was supported in part by PLGrid Infrastructure

A Langevin model for fluctuating contact angle behaviour parametrised using molecular dynamics

Molecular dynamics simulations are employed to develop a theoretical model to predict the fluid-solid contact angle as a function of wall-sliding speed incorporating thermal fluctuations. A liquid bridge between counter-sliding walls is studied, with liquid-vapour interface-tracking, to explore the impact of wall-sliding speed on contact angle. The behaviour of the macroscopic contact angle varies linearly over a range of capillary numbers beyond which the liquid bridge pinches off, a behaviour supported by experimental results. Nonetheless, the liquid bridge provides an ideal test case to study molecular scale thermal fluctuations, which are shown to be well described by Gaussian distributions. A Langevin model for contact angle is parametrised to incorporate the mean, fluctuation and auto-correlations over a range of sliding speeds and temperatures. The resulting equations can be used as a proxy for the fully-detailed molecular dynamics simulation allowing them to be integrated within a continuum-scale solver

What we observe is biased by what other people tell us: beliefs about the reliability of gaze behavior modulate attentional orienting to gaze cues

For effective social interactions with other people, information about the physical environment must be integrated with information about the interaction partner. In order to achieve this, processing of social information is guided by two components: a bottom-up mechanism reflexively triggered by stimulus-related information in the social scene and a top-down mechanism activated by task-related context information. In the present study, we investigated whether these components interact during attentional orienting to gaze direction. In particular, we examined whether the spatial specificity of gaze cueing is modulated by expectations about the reliability of gaze behavior. Expectations were either induced by instruction or could be derived from experience with displayed gaze behavior. Spatially specific cueing effects were observed with highly predictive gaze cues, but also when participants merely believed that actually non-predictive cues were highly predictive. Conversely, cueing effects for the whole gazed-at hemifield were observed with non-predictive gaze cues, and spatially specific cueing effects were attenuated when actually predictive gaze cues were believed to be non-predictive. This pattern indicates that (i) information about cue predictivity gained from sampling gaze behavior across social episodes can be incorporated in the attentional orienting to social cues, and that (ii) beliefs about gaze behavior modulate attentional orienting to gaze direction even when they contradict information available from social episodes