SAFT-γ force field for the simulation of molecular fluids: 8. Hetero-segmented coarse-grained models of perfluoroalkylalkanes assessed with new vapour–liquid interfacial tension data

Origin of the central pit in hemimicelles of semifluorinated alkanes: How molecular dipoles and substrate deformation can determine supra-molecular morphology

HYPOTHESIS

Semifluorinated alkanes amphiphiles spontaneously form highly monodispersed hemimicelles at the surface of water. The origin of the formation and complex structure of these surprising supramolecular aggregates were only recently clarified using molecular dynamics simulations (MD). The existence of a pit at the center of these aggregates made up of almost 3000 molecules was indeed reproduced by the MD simulations, but not understood.

METHOD

A careful strategy of atomistic MD simulations comparing non-electrostatic force fields with force fields that include electrostatic forces, thus bearing an implicit or explicit dipole, allowed demonstrating the roles of dipolar interactions and interactions with the liquid subphase on the morphology of the aggregates.

FINDINGS

The simulation results clearly show that within the hemimicelles the strong molecular dipoles located at the CH2-CF2 junctions tend to align, leading to a collective shift of the PFAA molecules relatively to each other. This shift is responsible for the curvature of the hemimicelles and originates the central pit, provided the possibility of deforming the surface of the water sub-phase. Comparisons with non-electrostatic force field results further contribute to understand the origin of the self-assembling process. The results directly connect for the first time a molecular property with a mesoscopic structural feature. Given the molecular simplicity of these "primitive" amphiphiles compared to the common hydrophilic/hydrophobic surfactants, the results contribute to understand self-assembly in general.

Multiscale simulation of fluids: coupling molecular and continuum

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.

Molecular dynamics simulation of the coalescence of surfactant-laden droplets

We investigate the coalescence of surfactant-laden water droplets by using several different surfactant types and a wide range of concentrations by means of a coarse-grained model obtained by the statistical associating fluid theory. Our results demonstrate in detail a universal mass transport mechanism of surfactant across many concentrations and several surfactant types during the process. Coalescence initiation is seen to occur via a single pinch due to aggregation of surface surfactant, and its remnants tend to become engulfed in part inside the forming bridge. Across the board we confirm the existence of an initial thermal regime with constant bridge width followed by a later inertial regime with bridge width scaling roughly as the square root of time, but see no evidence of an intermediate viscous regime. Coalescence becomes slower as surfactant concentration grows, and we see evidence of the appearance of a further slowdown of a different nature for several times the critical concentration. We anticipate that our results provide further insights in the mechanisms of coalescence of surfactant-laden droplets.

Breaking the Structure of Liquid Hydrogenated Alcohols Using Perfluorinated tert-Butanol: A Multitechnique Approach (Infrared, Raman, and X-ray Scattering) Analyzed by DFT and Molecular Dynamics Calculations

The state of aggregation at room temperature of tert-butanol (TBH) and perfluoro tert-butanol (TBF) liquid mixtures is assessed by vibrational spectroscopy (Raman and infrared) and X-ray diffraction and analyzed using density functional theory (DFT) and molecular dynamics (MD) simulations. It is shown that larger clusters (mostly tetramers) of TBH are destroyed upon dilution with TBF. Small oligomers, monomers, and mainly heterodimers are present at the equimolar concentration. At variance with slightly interacting solvents, the signature of hetero-oligomers is shown by the appearance of a new broad band detected in the infrared region. The same spectral observation is detected for mixtures of other hydrogenated alcohols (methanol and 1-butanol). The new infrared feature is unaffected by dilution in a polar solvent (CDCl₃) in a high-concentration domain, allowing us to assign it to the signature of small hetero-oligomers. MD simulations are used to assess the nature of the species present in the mixture (monomers and small hetero-oligomers) and to follow the evolution of their population upon the dilution. Combining MD simulations with DFT calculations, the infrared spectral profile is successfully analyzed in equimolecular mixtures. This study shows that TBF is a structure breaker of hydrogen-bonded alcohol networks and that the TBF (donor)-TBH (acceptor) heterodimer is the dominant species in an extended range of concentration, centered in the vicinity of the equimolar fraction.

Use of Boundary-Driven Nonequilibrium Molecular Dynamics for Determining Transport Diffusivities of Multicomponent Mixtures in Nanoporous Materials

The boundary-driven molecular modeling strategy to evaluate mass transport coefficients of fluids in nanoconfined media is revisited and expanded to multicomponent mixtures. The method requires setting up a simulation with bulk fluid reservoirs upstream and downstream of a porous media. A fluid flow is induced by applying an external force at the periodic boundary between the upstream and downstream reservoirs. The relationship between the resulting flow and the density gradient of the adsorbed fluid at the entrance/exit of the porous media provides for a direct path for the calculation of the transport diffusivities. It is shown how the transport diffusivities found this way relate to the collective, Onsager, and self-diffusion coefficients, typically used in other contexts to describe fluid transport in porous media. Examples are provided by calculating the diffusion coefficients of a Lennard-Jones (LJ) fluid and mixtures of differently sized LJ particles in slit pores, a realistic model of methane in carbon-based slit pores, and binary mixtures of methane with hypothetical counterparts having different attractions to the solid. The method is seen to be robust and particularly suited for the study of study of transport of dense fluids and liquids in nanoconfined media.

SAFT-γ-Mie Cross-Interaction Parameters from Density Functional Theory-Predicted Multipoles of Molecular Fragments for Carbon Dioxide, Benzene, Alkanes, and Water

Determining unlike-pair interaction parameters, whether for group contribution equation of state or molecular simulations, is a challenge for the prediction of thermodynamic properties. As the number of components and their respective complexity increase, it becomes impractical to fit all the unlike interactions. Lorentz-Berthelot combining rules work well for systems, where the main interactions are dispersion forces, but they do not account for electrostatics. In this work, we derive predictive combining rules within the SAFT-γ-Mie framework. In the resulting model, the unlike-pair interactions account for the effect of ionization energies, partial charges, dipole moments, and quadrupole moments. We then estimate these properties for molecular fragments using density functional theory calculations and demonstrate their use to obtain realistic cross-interaction energies without the need for experimental data. An open-source python package, Multipole Approach to Predictively Scale Cross-Interactions, is included to facilitate use of the methods presented in this work. A good qualitative agreement was obtained for all phase equilibria calculations of binary mixtures containing carbon dioxide with propane, hexane, benzene, and water, as well as mixtures of hexane and benzene. Finally, we discuss future improvements to our methodology, including the use of physical insights when fitting self-interaction parameters.

Development of accurate coarse-grained force fields for weakly polar groups by an indirect parameterization strategy

Coarse-grained (CG) molecular dynamics simulations are widely used to predict morphological structures and interpret mechanisms of mesoscopic behavior between the scope of traditional experiments and all-atom simulations. However, most current CG force fields (FFs) are not precise enough, especially for polar molecules or functional groups. A main obstacle in developing accurate CG FFs for polar molecules is the freezing problem met at room temperature. In this work, we introduce an indirect parametrization strategy for weakly polar groups by considering their short-chain homologs to avoid freezing. Here, a polar group containing three to four heavy atoms is mapped into one CG bead that is connected to one alkyl bead composed of three or four carbons. The CG beads interact via 4-parameter nonbonded Morse potentials and harmonic bonded potentials. An efficient meta-multilinear interpolation parameterization algorithm, as recently developed by us, is used to rigorously optimize the force parameters. Satisfactory accuracy is witnessed in terms of the density, heat of vaporization, surface tension, and solvation free energy of the homologs of twelve polar molecules, all deviating from the experiment by less than 5%. The transferability of the current FF is indicated by the predicted density, heat of vaporization, and end-to-end distance distributions of fatty acid methyl esters composed of multiple functional groups parameterized in this work.

Effect of Poly(vinyl butyral) Comonomer Sequence on Adhesion to Amorphous Silica: A Coarse-Grained Molecular Dynamics Study

Modulating a comonomer sequence, in addition to the overall chemical composition, is the key to unlocking the true potential of many existing commercial copolymers. We employ coarse-grained molecular dynamics (MD) simulations to study the behavior of random-blocky poly(vinyl butyral-co-vinyl alcohol) (PVB) melts in contact with an amorphous silica surface, representing the interface found in laminated safety glass. Our two-pronged coarse-graining approach utilizes both macroscopic thermophysical data and all-atom MD simulation data. Polymer-polymer nonbonded interactions are described by the fused-sphere SAFT-γ Mie equation of state, while bonded interactions are derived using Boltzmann inversion to match the bond and angle distributions from all-atom PVB chains. Spatially dependent polymer-surface interactions are mapped from a hydroxylated all-atom amorphous silica slab model and all-atom monomers to an external potential acting on the coarse-grained sites. We discovered an unexpected complex relationship between the blockiness parameter and the adhesion energy. The adhesion strength between PVB copolymers with intermediate VA content and silica was found to be maximal for random-blocky copolymers with a moderately high degree of blockiness rather than for diblock copolymers. We attribute this to two main factors: (1) changes in morphology, which dramatically alter the number of VA beads interacting with the surface and (2) a non-negligible contribution of vinyl butyral (VB) monomers to adhesion energy because of their preference to adsorb to zones with low hydroxyl density on the silica surface.

From nano-seggregation to mesophases: probing the liquid structure of perfluoroalkylalkanes with 129Xe NMR spectroscopy

In this work we demonstrate that pure perfluoroalkylalkane diblock molecules are not isotropic liquids and self organize forming domains at the nanometric scale. 129Xe NMR spectra were obtained as a function of temperature for seven liquid perfluoroalkylalkanes, covering a range of relative lengths of the hydrogenated and fluorinated segments. The results support the presence of domains richer in the hydrogenated groups, in which xenon is preferentially dissolved. The average local concentration within the xenon coordination sphere is estimated to be 0.05 mole fraction higher in hydrogenated groups than the stoichiometric proportion. Atomistic molecular dynamics simulations support this analysis and allow a detailed analysis of the liquid structure. Furthermore, 129Xe NMR spectra in perfluorohexylhexane (F6H6) and perfluorohexyloctane (F6H8) obtained as a function of temperature, clearly detect the existence of two distinct environments in the fluid, one richer in hydrogenated groups and another richer in fluorinated groups, consistent with the formation of mesophases. It is important to stress that nano-segregation is this case observed in liquids interacting exclusively through dispersion forces, unlike most common examples of segregation which are determined by hydrogen bonding and polarity. Given the simple molecular structure and interactions of the studied PFAA, we believe that the present results can have a general impact in understanding the early mechanisms of segregation, phase separation and self-assembly.

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.

Spontaneous self-assembly and structure of perfluoroalkylalkane surfactant hemimicelles by molecular dynamics simulations

Fully atomistic molecular-dynamics (MD) simulations of perfluoroalkylalkane molecules at the surface of water show the spontaneous formation of aggregates whose size and topography closely resemble the experimentally observed hemimicelles for this system. Furthermore, the grazing incidence X-ray diffraction (GIXD) pattern calculated from the simulation trajectories reproduces the experimental GIXD spectra previously obtained, fully validating the MD simulation results. The detailed analysis of the internal structure of the aggregates obtained by the MD simulations supports a definite rational explanation for the spontaneous formation, stability, size, and shape of perfluoroalkylalkane hemimicelles at the surface of water.