Reconsidering Dispersion Potentials: Reduced Cutoffs in Mesh-Based Ewald Solvers Can Be Faster Than Truncation

Impact of the Force Field on the Calculation of Density and Surface Tension of Epoxy-Resins

The molecular simulation of interfacial systems is a matter of debate because of the choice of many input parameters that can affect significantly the performance of the force field of reproducing the surface tension and the coexisting densities. After developing a robust methodology for the calculation of the surface tension on a Lennard-Jones fluid, we apply it with different force fields to calculate the density and surface tension of pure constituents of epoxy resins. By using the model that best reproduces the experimental density and surface tension, we investigate the impact of composition in mass fraction on uncured epoxy resins and the effects of degree of cross-linking on cured resins.

Crosslinking Rapidly Cured Epoxy Resin Thermosets: Experimental and Computational Modeling and Simulation Study

The power of computational modeling and simulation for establishing clear links between materials' intrinsic properties and their atomic structure has more and more increased the demand for reliable and reproducible protocols. Despite this increased demand, no one approach can provide reliable and reproducible outcomes to predict the properties of novel materials, particularly rapidly cured epoxy-resins with additives. This study introduces the first computational modeling and simulation protocol for crosslinking rapidly cured epoxy resin thermosets based on solvate ionic liquid (SIL). The protocol combines several modeling approaches, including quantum mechanics (QMs) and molecular dynamics (MDs). Furthermore, it insightfully provides a wide range of thermo-mechanical, chemical, and mechano-chemical properties, which agree with experimental data.

Interfacial Properties of Linear Alkane/Nitrogen Binary Mixtures: Molecular Dynamics Vapor-Liquid Equilibrium Simulations

Molecular dynamics simulations were used to investigate the vapor-liquid equilibria (VLE) and interfacial properties of binary mixtures of N2 with either ethane, propane, n-decane, or n-dodecane. Alkanes and N2 were modeled by using the TraPPE-UA and Rivera force fields, respectively. The typically used Lorentz-Berthelot combining rules resulted in liquid phases that are too N2-rich compared to experiment. To improve the accuracy of VLE predictions, the hydrocarbon-nitrogen interactions were fine-tuned, and these improved parameters were used to investigate interfacial properties. Scaling the interaction strength between nitrogen and -CH3 and -CH2- groups by factors of 0.95 and 0.85, respectively, relative to the Lorentz-Berthelot value, was found to minimize error in pressure-composition phase diagrams. These scaling parameters gave excellent agreement with experimental phase diagrams for mixtures of N2 with ethane, propane, or n-dodecane over a range of state points. For ethane/N2 and n-decane/N2 mixtures, trends in surface tension as a function of temperature and pressure are correctly reproduced, although the simulated values are slightly too high compared to experimental values. To assess how the accuracy of hydrocarbon-N2 interaction strength impacts interfacial property predictions, we have compared density profiles and surface tension using several different scaling factors. Using the Lorentz-Berthelot combining rules rather than optimized parameters gave the same qualitative trends, although some quantitative results, such as liquid-phase N2 mole fraction, were found to differ by a factor of ∼1.5. Using the optimized interaction parameters, interfacial behavior was examined by calculating density and free energy profiles. Nitrogen molecules preferentially adsorb at the interfacial region between the liquid and vapor phases. This interfacial adsorption becomes less energetically favorable as either the temperature, pressure, or length of the alkane chain increases.

First Hyperpolarizability of Water in Bulk Liquid Phase: Long-Range Electrostatic Effects Included via the Second Hyperpolarizability

The molecular first hyperpolarizability β contributes to second-order optical non-linear signals collected from molecular liquids. For the Second Harmonic Generation (SHG) response, the first hyperpolarizability β (2ω,ω,ω) often depends on...

Interfacial Properties of Hydrophobic Deep Eutectic Solvents with Water

Hydrophobic deep eutectic solvents (DESs) have recently gained much attention as water-immiscible solvents for a wide range of applications. However, very few studies exist in which the hydrophobicity of these DESs is quantified. In this work, the interfacial properties of hydrophobic DESs with water were computed at various temperatures using molecular dynamics simulations. The considered DESs were tetrabutylammonium chloride–decanoic acid (TBAC–dec) with a molar ratio of 1:2, thymol–decanoic acid (Thy–dec) with a molar ratio of 1:2, and dl-menthol–decanoic acid (Men–dec) with a molar ratio of 2:1. The following properties were investigated in detail: interfacial tensions, water-in-DES solubilities (and salt-in-water solubilities for TBAC–dec/water), density profiles, and the number densities of hydrogen bonds. Different ionic charge scaling factors were used for TBAC–dec. Thy–dec and Men–dec showed a high level of hydrophobicity with negligible computed water-in-DES solubilities. For charge scaling factors of 0.7 and 1 for the thymol and decanoic acid components of Thy–dec, the computed interfacial tensions of the DESs are in the following order: TBAC–dec (ca. 4 mN m^–1) < Thy–dec (20 mN m^–1) < Men–dec (26 mN m^–1). The two sets of charge scaling factors for Thy–dec did not lead to different density profiles but resulted in considerable differences in the DES/water interfacial tensions due to different numbers of decanoic acid–water hydrogen bonds at the interfaces. Large peaks were observed for the density profiles of (the hydroxyl oxygen of) decanoic acid at the interfaces of all DES/water mixtures, indicating a preferential alignment of the oxygen atoms of decanoic acid toward the aqueous phase.

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor β(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

We use molecular dynamics to calculate the rotational and vibrational energy relaxation of C₂H₆ in Ar, Kr, and Xe bath gases over a pressure range of 10-400 atm and at temperatures of 300 and 800 K. The C₂H₆ is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a nonlinear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures does not yet show any nonlinearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure- and temperature-dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the nonlinear pressure dependence of rotational relaxation.

Evaluation of nine condensed-phase force fields of the GROMOS, CHARMM, OPLS, AMBER, and OpenFF families against experimental cross-solvation free energies

Experimental solvation free energies are nowadays commonly included as target properties in the validation of condensed-phase force fields, sometimes even in their calibration. In a previous article [Kashefolgheta et al., J. Chem. Theory. Comput., 2020, 16, 7556-7580], we showed how the involved comparison between experimental and simulation results could be made more systematic by considering a full matrix of cross-solvation free energies . For a set of N molecules that are all in the liquid state under ambient conditions, such a matrix encompasses N×N entries for considering each of the N molecules either as solute (A) or as solvent (B). In the quoted study, a cross-solvation matrix of 25 × 25 experimental value was introduced, considering 25 small molecules representative for alkanes, chloroalkanes, ethers, ketones, esters, alcohols, amines, and amides. This experimental data was used to compare the relative accuracies of four popular condensed-phase force fields, namely GROMOS-2016H66, OPLS-AA, AMBER-GAFF, and CHARMM-CGenFF. In the present work, the comparison is extended to five additional force fields, namely GROMOS-54A7, GROMOS-ATB, OPLS-LBCC, AMBER-GAFF2, and OpenFF. Considering these nine force fields, the correlation coefficients between experimental values and simulation results range from 0.76 to 0.88, the root-mean-square errors (RMSEs) from 2.9 to 4.8 kJ mol-1, and average errors (AVEEs) from -1.5 to +1.0 kJ mol-1. In terms of RMSEs, GROMOS-2016H66 and OPLS-AA present the best accuracy (2.9 kJ mol-1), followed by OPLS-LBCC, AMBER-GAFF2, AMBER-GAFF, and OpenFF (3.3 to 3.6 kJ mol-1), and then by GROMOS-54A7, CHARM-CGenFF, and GROMOS-ATB (4.0 to 4.8 kJ mol-1). These differences are statistically significant but not very pronounced, and are distributed rather heterogeneously over the set of compounds within the different force fields.

Molecular modeling of aqueous electrolytes at interfaces: Effects of long-range dispersion forces and of ionic charge rescaling

Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions-described by a truncated Lennard-Jones (LJ) potential-and electrostatic interactions-described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges [electronic continuum correction (ECC)], possibly complemented with rescaling of LJ parameters [ECC rescaled (ECCR)], have shown promising results in bulk, but their performance at interfaces has been less explored. Here, we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk accurately describes image charge repulsion at the liquid-vapor interface, and the rescaling of LJ parameters in ECCR models-aimed at capturing local ion-ion and ion-water interactions in bulk- describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.

Insight into Cellulose Dissolution with the Tetrabutylphosphonium Chloride-Water Mixture using Molecular Dynamics Simulations

All-atom molecular dynamics simulations are utilized to determine the properties and mechanisms of cellulose dissolution using the ionic liquid tetrabutylphosphonium chloride (TBPCl)–water mixture, from 63.1 to 100 mol % water. The hydrogen bonding between small and large cellulose bundles with 18 and 88 strands, respectively, is compared for all concentrations. The Cl, TBP, and water enable cellulose dissolution by working together to form a cooperative mechanism capable of separating the cellulose strands from the bundle. The chloride anions initiate the cellulose breakup, and water assists in delaying the cellulose strand reformation; the TBP cation then more permanently separates the cellulose strands from the bundle. The chloride anion provides a net negative pairwise energy, offsetting the net positive pairwise energy of the peeling cellulose strand. The TBP–peeling cellulose strand has a uniquely favorable and potentially net negative pairwise energy contribution in the TBPCl–water solution, which may partially explain why it is capable of dissolving cellulose at moderate temperatures and high water concentrations. The cellulose dissolution declines rapidly with increasing water concentration as hydrogen bond lifetimes of the chloride–cellulose hydroxyl hydrogens fall below the cellulose’s largest intra-strand hydrogen bonding lifetime.

Thermodynamic and Transport Properties of Tetrabutylphosphonium Hydroxide and Tetrabutylphosphonium Chloride-Water Mixtures via Molecular Dynamics Simulation

Thermodynamic, structural, and transport properties of tetrabutylphosphonium hydroxide (TBPH) and tetrabutylphosphonium chloride (TBPCl)–water mixtures have been investigated using all-atom molecular dynamics simulations in response to recent experimental work showing the TBPH–water mixtures capability as a cellulose solvent. Multiple transitional states exist for the water—ionic liquid (IL) mixture between 70 and 100 mol% water, which corresponds to a significant increase in water hydrogen bonds. The key transitional region, from 85 to 92.5 mol% water, which coincides with the mixture’s maximum cellulose solubility, reveals small and distinct water veins with cage structures formed by the TBP⁺ ions, while the hydroxide and chloride ions have moved away from the P atom of TBP⁺ and are strongly hydrogen bonded to the water. The maximum cellulose solubility of the TBPH–water solution at approximately 91.1 mol% water, appears correlated with the destruction of the TBP’s interlocking structure in the simulations, allowing the formation of water veins and channeling structures throughout the system, as well as changing from a subdiffusive to a near-normal diffusive regime, increasing the probability of the IL’s interaction with the cellulose polymer. A comparison is made between the solution properties of TBPH and TBPCl with those of alkylimidazolium-based ILs, for which water appears to act as anti-solvent rather than a co-solvent.

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.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.

Sparsification of long range force networks for molecular dynamics simulations

Atomic interactions in solid materials are described using network theory. The tools of network theory focus on understanding the properties of a system based upon the underlying interactions which govern their dynamics. While the full atomistic network is dense, we apply a spectral sparsification technique to construct a sparse interaction network model that reduces the computational complexity while preserving macroscopic conservation properties. This sparse network is compared to a reduced network created using a cut-off radius (threshold method) that is commonly used to speed-up computations while approximating interatomic forces. The approximations used to estimate the total forces on each atom are quantified to assess how local interatomic force errors propagate errors at the global or continuum scale by comparing spectral sparsification to thresholding. In particular, we quantify the performance of the spectral sparsification algorithm for the short-range Lennard-Jones potential and the long-range Coulomb potential. Spectral sparsification of the Lennard–Jones potential yields comparable results to thresholding while spectral sparsification yields improvements when considering a long-range Coulomb potential. The present network-theoretic formulation is implemented on two sample problems: relaxation of atoms near a surface and a tensile test of a solid with a circular hole.

Importance of the tail corrections on surface tension of curved liquid-vapor interfaces

We report molecular simulations of the liquid-vapor cylindrical interface of methane. We apply the truncated Lennard-Jones potential and specific long-range corrections for the surface tension developed especially for cylindrical interfaces. We investigate the impact of the cutoff on the radial density profile, the intrinsic and long-range correction parts to the surface tension, and Tolman length. We also study the curvature dependence of the surface tension as a function of the cutoff used. In this work we shed light that both density and Tolman length are cutoff-dependent whereas the total surface tension is slightly curvature and cutoff dependent.

Atomistic and Coarse-Grained Modeling of the Adsorption of Graphene Nanoflakes at the Oil-Water Interface

The high interfacial tension between two immiscible liquids can provide the necessary driving force for the self-assembly of nanoparticles at the interface. Particularly, the interface between water and oily liquids (hydrocarbon chains) has been exploited to prepare networks of highly interconnected graphene sheets of only a few layers thickness, which are well suited for industrial applications. Studying such complex systems through particle-based simulations could greatly enhance the understanding of the various driving forces in action and could possibly give more control over the self-assembly process. However, the interaction potentials used in particle-based simulations are typically derived by reproducing bulk properties and are therefore not suitable for describing systems dominated by interfaces. To address this issue, we introduce a methodology to derive solid-liquid interaction potentials that yield an accurate representation of the balance between interfacial interactions at atomistic and coarse-grained resolutions. Our approach is validated through its ability to lead to the adsorption of graphene nanoflakes at the interface between water and n-hexane. The development of accurate coarse-grained potentials that our approach enables will allow us to perform large-scale simulations to study the assembly of graphene nanoparticles at the interface between immiscible liquids. Our methodology is illustrated through a simulation of many graphene nanoflakes adsorbing at the interface.

Comparison of Additive and Polarizable Models with Explicit Treatment of Long-Range Lennard-Jones Interactions Using Alkane Simulations

Long-range Lennard-Jones (LJ) interactions have a significant impact on the structural and thermodynamic properties of nonpolar systems. While several methods have been introduced for the treatment of long-range LJ interactions in molecular dynamics (MD) simulations, increased accuracy and extended applicability is required for anisotropic systems such as lipid bilayers. The recently refined Lennard-Jones particle-mesh Ewald (LJ-PME) method extends the particle-mesh Ewald (PME) method to long-range LJ interactions and is suitable for use with anisotropic systems. Implementation of LJ-PME with the CHARMM36 (C36) additive and CHARMM Drude polarizable force fields improves agreement with experiment for density, isothermal compressibility, surface tension, viscosity, translational diffusion, and ¹³C T₁ relaxation times of pure alkanes. Trends in the temperature dependence of the density and isothermal compressibility of hexadecane are also improved. While the C36 additive force field with LJ-PME remains a useful model for liquid alkanes, the Drude polarizable force field with LJ-PME is more accurate for nearly all quantities considered. LJ-PME is also preferable to the isotropic long-range correction for hexadecane because the molecular order extends to nearly 20 Å, well beyond the usual 10-12 Å cutoffs used in most simulations.

Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method

The smooth particle mesh Ewald (SPME) method is the standard method for computing the electrostatic interactions in the molecular simulations. In this work, we develop the multiple staggered mesh Ewald (MSME) method, which averages the SPME forces computed on, for example, M, staggered meshes. We prove, from a theoretical perspective, that the MSME is as accurate as the SPME, but uses M² times less mesh points in a certain parameter range. In the complementary parameter range, the MSME is as accurate as the SPME with twice the interpolation order. The theoretical conclusions are numerically validated both by a uniform and uncorrelated charge system, and by a three-point-charge water system that is widely used as solvent for the biomacromolecules.

Classification of precursors in nanoscale droplets

Molecular precursors, ultrathin films that precede spreading droplets, are still far from being understood, despite intensive study. The inherent microscopic length scales make small-scale experimental techniques and molecular simulation ideal methods to study this phenomenon. Previous work on molecular precursors using nanoscale droplets, however, consistently suffers from incorrect measurement of the dimensions of the precursor film. An alternative method to accurately characterize the precursor film is presented here. In contrast to previous measures, this method (i) allows for easy detection and characterization of precursors and (ii) yields wetting dynamics that agree with experimental observations. Finally, we briefly comment on previous studies whose conclusions may merit reconsideration in light of the present work.

Requirements for the Formation and Shape of Microscopic Precursors in Droplet Spreading

While the mass transport mechanisms and dynamics of molecular precursors, ultrathin films that precede spreading droplets, are well understood, the requirements for their formation and the reasons for the occurrence of different precursor shapes remain unclear. In this work, we study these requirements using molecular dynamics simulations of spreading droplets and extensive free energy computations. For a simple simulation model, we demonstrate that with growing droplet-substrate attraction, spreading passes succesively through regimes with no precursor, a monolayer precursor, and a continuously growing precursor. We show that the onset of layer formation and the changes in the precursor shape correlate with the free energy of layer formation. On the basis of these findings, we show that a positive spreading coefficient is sufficient but not necessary for precursor formation.

Computer modelling of the surface tension of the gas–liquid and liquid–liquid interface

This review presents the state of the art in molecular simulations of interfacial systems and of the calculation of the surface tension from the underlying intermolecular potential.

Thermophysical properties of imidazolium tricyanomethanide ionic liquids: experiments and molecular simulation

The low-viscous tricyanomethanide ([TCM]⁻)-based ionic liquids (ILs) are gaining increasing interest as attractive fluids for a variety of industrial applications.

Direct-Space Corrections Enable Fast and Accurate Lorentz-Berthelot Combination Rule Lennard-Jones Lattice Summation

Long-range lattice summation techniques such as the particle-mesh Ewald (PME) algorithm for electrostatics have been revolutionary to the precision and accuracy of molecular simulations in general. Despite the performance penalty associated with lattice summation electrostatics, few biomolecular simulations today are performed without it. There are increasingly strong arguments for moving in the same direction for Lennard-Jones (LJ) interactions, and by using geometric approximations of the combination rules in reciprocal space, we have been able to make a very high-performance implementation available in GROMACS. Here, we present a new way to correct for these approximations to achieve exact treatment of Lorentz-Berthelot combination rules within the cutoff, and only a very small approximation error remains outside the cutoff (a part that would be completely ignored without LJ-PME). This not only improves accuracy by almost an order of magnitude but also achieves absolute biomolecular simulation performance that is an order of magnitude faster than any other available lattice summation technique for LJ interactions. The implementation includes both CPU and GPU acceleration, and its combination with improved scaling LJ-PME simulations now provides performance close to the truncated potential methods in GROMACS but with much higher accuracy.

Smoothing of contact lines in spreading droplets by trisiloxane surfactants and its relevance for superspreading

Superspreading, the greatly enhanced spreading of aqueous solutions of trisiloxane surfactants on hydrophobic substrates, is of great interest in fundamental physics and technical applications. Despite numerous studies in the last 20 years, the superspreading mechanism is still not well understood, largely because the molecular scale cannot be resolved appropriately either experimentally or using continuum simulations. The absence of molecular-scale knowledge has led to a series of conflicting hypotheses based on different assumptions of surfactant behavior. We report a series of large-scale molecular dynamics simulations of aqueous solutions of superspreading and non-superspreading surfactants on different substrates. We find that the transition from the liquid-vapor to the solid-liquid interface is smooth for superspreading conditions, allowing direct adsorption through the contact line. This finding complements a study [Karapetsas et al., J. Fluid Mech., 2011, 670, 5-37], which predicts that superspreading can occur if this adsorption path is possible. Based on the observed mechanism, we provide plausible explanations for the influence of the substrate hydrophobicity, the surfactant chain length, and the surfactant concentration on the superspreading phenomenon. We also briefly address that the observed droplet shape is a mechanism to overcome the Huh-Scriven paradox of infinite viscous dissipation at the contact line.

Atomistic potentials for trisiloxane, alkyl ethoxylate, and perfluoroalkane-based surfactants with TIP4P/2005 and application to simulations at the air-water interface

The mechanism of superspreading, the greatly enhanced spreading of water droplets facilitated by trisiloxane surfactants, is still under debate, largely because the role and behavior of the surfactants cannot be sufficiently resolved by experiments or continuum simulations. Previous molecular dynamics studies have been performed with simple model molecules or inaccurate models, strongly limiting their explanatory power. Here we present a force field dedicated to superspreading, extending existing quantum-chemistry-based models for the surfactant and the TIP4P/2005 water model ( Abascal et al. J. Chem. Phys. , 2005 , 123 , 234505 ). We apply the model to superspreading trisiloxane surfactants and nonsuperspreading alkyl ethoxylate and perfluoroalkane surfactants at various concentrations at the air-water interface. We show that the developed model accurately predicts surface tensions, which are typically assumed important for superspreading. Significant differences between superspreading and traditional surfactants are presented and their possible relation to superspreading discussed. Although the force field has been developed for superspreading problems, it should also perform well for other simulations involving polymers or copolymers with water.