Instantaneous liquid interfaces

Structure and Dynamics of Water in Polysaccharide (Alginate) Solutions and Gels Explained by the Core-Shell Model

In both biological and engineered systems, polysaccharides offer a means of establishing structural stiffness without altering the availability of water. Notable examples include the extracellular matrix of prokaryotes and eukaryotes, artificial skin grafts, drug delivery materials, and gels for water harvesting. Proper design and modeling of these systems require detailed understanding of the behavior of water confined in pores narrower than about 1 nm. We use molecular dynamics simulations to investigate the properties of water in solutions and gels of the polysaccharide alginate as a function of the water content and polymer cross-linking. We find that a detailed understanding of the nanoscale dynamics of water in alginate solutions and gels requires consideration of the discrete nature of water. However, we also find that the trends in tortuosity, permeability, dielectric constant, and shear viscosity can be adequately represented using the "core-shell" conceptual model that considers the confined fluid as a continuum.

How Thick is the Air-Water Interface?─A Direct Experimental Measurement of the Decay Length of the Interfacial Structural Anisotropy

The air-water interface is a highly prevalent phase boundary impacting many natural and artificial processes. The significance of this interface arises from the unique properties of water molecules within the interfacial region, with a crucial parameter being the thickness of its structural anisotropy, or "healing depth". This quantity has been extensively assessed by various simulations which have converged to a prediction of a remarkably short length of ∼6 Å. Despite the absence of any direct experimental measurement of this quantity, this predicted value has surprisingly become widely accepted as fact. Using an advancement in nonlinear vibrational spectroscopy, we provide the first measurement of this thickness and, indeed, find it to be ∼6-8 Å, finally confirming the prior predictions. Lastly, by combining the experimental results with depth-dependent second-order spectra calculated from ab initio parametrized molecular dynamics simulations, which are also in excellent agreement with this experimental result, we shed light on this surprisingly short correlation length of molecular orientations at the interface.

The wetting of H2O by CO2

Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO2-H2O interface, in particular, has been the subject of a number of studies on account of its importance for the carbon life cycle as well as carbon capture and sequestration schemes. Despite this attention, there remain a number of open questions on the nature of the CO2-H2O interface, particularly concerning the interfacial tension and phase behavior of CO2 at the interface. In this paper, we seek to address these ambiguities using ab initio-quality simulations. Harnessing the benefits of machine-learned potentials and enhanced statistical sampling methods, we present an ab initio-level description of the CO2-H2O interface. Interfacial tensions are predicted from 1 to 500 bars and found to be in close agreement with experiment at pressures for which experimental data are available. Structural analyses indicate the buildup of an adsorbed, saturated CO2 film forming at a low pressure (20 bars) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. The CO2 monolayer buildup coincides with a reduced structuring of water molecules close to the interface. This study highlights the predictive nature of machine-learned potentials for complex macroscopic properties of biphasic interfaces, and the mechanistic insight obtained into carbon dioxide aggregation at the water interface is of high relevance for geoscience, climate research, and materials science.

Effects of Surface Charge Distribution and Electrolyte Ions on the Nonlinear Spectra of Model Solid-Water Interfaces

Molecular dynamics simulations of model charged solid/water interfaces were carried out to provide insight about the relationship between the second-order nonlinear susceptibility χ(2) and the structure of the interfacial water layer. The results of the calculations reveal that the density fluctuations of water extend to about 12 Å from the surface regardless of the system, while the orientational ordering of molecules is long-ranged and is sensitive to the presence of electrolytes. The charge localization on the surface was found to affect only the high-frequency part of the Im[χ(2)] spectrum, and the addition of salt has very little effect on the spectrum of the first water layer. For solid/neat water interfaces, the spectroscopically active part of the liquid phase has a thickness largely exceeding the region of density fluctuations, and this long-ranged nonlinear activity is mediated by the electric field of the molecules. The electrolyte ions and their hydration shells act in a destructive way on the molecular field. This effect, combined with the screening of the surface charge by ions, drastically reduces the thickness of the spectroscopic diffuse layer. There is an electrolyte concentration at which the nonlinear response of the diffuse layer is suppressed and the χ(2) spectrum of the interface essentially coincides with that of the first water layer.

Water-air interface revisited by means of path-integral ab initio molecular dynamics

Although nuclear quantum effects (NQEs) have been considered on bulk liquid water, the impact of these latter on the air-water interface has not yet been reported. Herein, by performing and comparing ab initio molecular dynamics (AIMD) and path integral AIMD (PI-AIMD) simulations, we reveal the impact of NQEs on structural, dynamical and electronic properties as well as IR spectra of the air-water interface at room temperature. NQEs, being able to describe a more accurate proton delocalization in H-bonded system than AIMD, reveal a different structural arrangement and dynamical behaviour of both bulk and interfacial water molecules in comparison to AIMD results. A more de-structured and de-bound water arrangement and coordination are identified when the quantum nature of nuclei are considered for both bulk and interfacial water molecules. Structural properties, such as inter-/intra-molecular bond lengths, coordination numbers and H-bonding angles of bulk and interfacial water molecules here calculated, are affected by NQEs mitigating the overstructured description given by AIMD. Further evidences of an AIMD overstructured description of bulk water are in the computed IR spectra, where an increased absorption peak intensity and an increased strength of the hydrogen-bond network are alleviated by NQEs. In addition, NQEs show a valuable impact on the electronic structure of the air-water interface, reducing the total valence bandwidth and the electronic energy band-gap when passing from bulk to interfacial water. This work proves how NQEs significantly affect properties and features of the air-water interface, that are essential to accurately describe H-bonded interfacial systems.

How Molecular Orientation Affects the Static Permittivity Profile of the Polar and Nonpolar Liquid-Liquid Interface

The dielectric permittivity across the liquid-liquid interface presents an intrinsic response, with respect to the instantaneous interface reference. We hypothesize that dielectric responses across the nonpolar and polar liquid-liquid interfaces have different behaviors and underlying mechanisms. Molecular dynamics simulations were used to compare and contrast the dielectric response of a nonpolar (1,2-dichloroethane/water) and a polar (1-octanol/water) liquid-liquid interface system. We found that the enhanced dielectric permittivity at the nonpolar interface is attributed to the increased water dipole orientation and polarization density. In the case of the polar interface, strong association of the immiscible solvents inhibits the molecular dipole orientation, counteracting the effect from the enhanced surface water polarization density and resulting in a standard dielectric response. Detailed knowledge of the hydrogen bond networks and molecular dipole orientation with respect to the specific instantaneous interfacial and bulk regions reveals the effect of molecular proximity and the interaction with the opposing interfacial molecules on the mechanism of the dielectric permittivity response across the liquid-liquid interface phase boundary.

Influence of cholesterol on hydrogen-bond dynamics of water molecules in lipid-bilayer systems at varying temperatures

Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting a considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoylphosphatidylcholine (DPPC) and palmitoyl sphingomyelin; they are distinct with respect to the structures of the hydrogen-bond (H-bond) acceptors. Our investigation focuses on the dynamic properties and H-bonds of water molecules in the lipid-membrane systems, with a particular emphasis on the influence of Chol at varying temperatures. Notably, in the gel phase at 303 K, the presence of Chol extends the lifetimes of H-bonds of the oxygen atoms acting as H-bond acceptors within DPPC with water molecules by a factor of 1.5-2.5. In the liquid-crystalline phase at 323 K, on the other hand, H-bonding dynamics with lipid membranes remain largely unaffected by Chol. This observed shift in H-bonding states serves as a crucial key to unraveling the subtle control mechanisms governing water dynamics in lipid-membrane systems.

Spatial organization of the ions at the free surface of imidazolium-based ionic liquids

HYPOTHESIS

Experimental information on the molecular scale structure of ionic liquid interfaces is controversial, giving rise to two competing scenarios, namely the double layer-like and "chessboard"-like structures. This issue can be resolved by computer simulation methods, at least for the underlying molecular model. Systematically changing the anion type can elucidate the relative roles of electrostatic interactions, hydrophobic (or, strictly speaking, apolar) effects and steric restrictions on the interfacial properties.

SIMULATIONS

Molecular dynamics simulation is combined with intrinsic analysis methods both at the molecular and atomic levels, supplemented by Voronoi analysis of self-association.

FINDINGS

We see no evidence for the existence of a double-layer-type arrangement of the ions, or for their self-association at the surface of the liquid. Instead, our results show that cation chains associate into apolar domains that protrude into the vapour phase, while charged groups form domains that are embedded in this apolar environment at the surface. However, the apolar chains largely obscure the cation groups, to which they are bound, while the smaller and more mobile anions can more easily access the free surface, leading to a somewhat counterintuitive net excess of negative charge at the interface. Importantly, this excess charge could only be identified by applying intrinsic analysis.

Shape of AOT Reverse Micelles: The Mesoscopic Assembly Is More Than the Sum of the Parts

AOT reverse micelles are a common and convenient model system for studying the effects of nanoconfinement on aqueous solutions. The reverse micelle shape is important to understanding how the constituent components come together to form the coherent whole and the unique properties observed there. The shape of reverse micelles impacts the amount of interface present and the distance of the solute from the interface and is therefore vital to understanding interfacial properties and the behavior of solutes in the polar core. In this work, we use previously introduced measures of shape, the coordinate-pair eccentricity (CPE) and convexity, and apply them to a series of simulations of AOT reverse micelles. We simulate the most commonly used force field for AOT reverse micelles, the CHARMM force field, but we also adapt the OPLS force field for use with AOT, the first work to do so, in addition to using both 3- and 4-site water models. Altogether, these simulations are designed to examine the impact of the force field on the shape of the reverse micelles in detail. We also study the time autocorrelation of shape, the water rotational anisotropy decay, and how the CPE changes between the water pool and AOT tail groups. We find that although the force field changes the shape noticeably, AOT reverse micelles are always amorphous particles. The shape of the micelles changes on the order of 10 ns. The water rotational dynamics observed match the experiment and demonstrate slower dynamics relative to bulk water, suggesting a two-population model that fits a core/shell hypothesis. Taken together, our results indicate that it is likely not possible to create a perfect force field that can reproduce every aspect of the AOT reverse micelle accurately. However, the magnitude of the differences between simulations appears relatively small, suggesting that any reasonably derived force field should provide an acceptable model for most work on AOT reverse micelles.

Impact of interfacial curvature on molecular properties of aqueous interfaces

The curvature of soft interfaces plays a crucial role in determining their mechanical and thermodynamic properties, both at macroscopic and microscopic scales. In the case of air/water interfaces, particular attention has recently focused on water microdroplets, due to their distinctive chemical reactivity. However, the specific impact of curvature on the molecular properties of interfacial water and interfacial reactivity has so far remained elusive. Here, we use molecular dynamics simulations to determine the effect of curvature on a broad range of structural, dynamical, and thermodynamical properties of the interface. For a droplet, a flat interface, and a cavity, we successively examine the structure of the hydrogen-bond network and its relation to vibrational spectroscopy, the dynamics of water translation, rotation, and hydrogen-bond exchanges, and the thermodynamics of ion solvation and ion-pair dissociation. Our simulations show that curvature predominantly impacts the hydrogen-bond structure through the fraction of dangling OH groups and the dynamics of interfacial water molecules. In contrast, curvature has a limited effect on solvation and ion-pair dissociation thermodynamics. For water microdroplets, this suggests that the curvature alone cannot fully account for the distinctive reactivity measured in these systems, which are of great importance for catalysis and atmospheric chemistry.

Molecular Insights into Chemical Reactions at Aqueous Aerosol Interfaces

Atmospheric aerosols facilitate reactions between ambient gases and dissolved species. Here, we review our efforts to interrogate the uptake of these gases and the mechanisms of their reactions both theoretically and experimentally. We highlight the fascinating behavior of N2O5 in solutions ranging from pure water to complex mixtures, chosen because its aerosol-mediated reactions significantly impact global ozone, hydroxyl, and methane concentrations. As a hydrophobic, weakly soluble, and highly reactive species, N2O5 is a sensitive probe of the chemical and physical properties of aerosol interfaces. We employ contemporary theory to disentangle the fate of N2O5 as it approaches pure and salty water, starting with adsorption and ending with hydrolysis to HNO3, chlorination to ClNO2, or evaporation. Flow reactor and gas-liquid scattering experiments probe even greater complexity as added ions, organic molecules, and surfactants alter the interfacial composition and reaction rates. Together, we reveal a new perspective on multiphase chemistry in the atmosphere.

Computational insight into the effect of alkyl chain length in tetraalkylammonium-based deep eutectic solvents

The effect of the increase in the alkyl chain length of cation on the properties of deep eutectic solvents based on ethylene glycol has been investigated employing classical molecular dynamics simulations. The change in the structural and dynamic properties in both the bulk and liquid-vapor interface is explored through various analyses. The interaction between the anion and the ethylene glycol increases with an increase in the alkyl chain length of the cation, as observed in the increase of the lifetime of the hydrogen bond formed between the two. The terminal carbon atoms are found to be closer to each other when the cation changes from tetraethylammonium to tetrabutylammonium. The cations are located closer to the interface, and the association of the alkyl chains becomes more significant with increased alkyl chain length, decreasing the surface tension values.

Surface stratification determines the interfacial water structure of simple electrolyte solutions

The distribution of ions at the air/water interface plays a decisive role in many natural processes. Several studies have reported that larger ions tend to be surface-active, implying ions are located on top of the water surface, thereby inducing electric fields that determine the interfacial water structure. Here we challenge this view by combining surface-specific heterodyne-detected vibrational sum-frequency generation with neural network-assisted ab initio molecular dynamics simulations. Our results show that ions in typical electrolyte solutions are, in fact, located in a subsurface region, leading to a stratification of such interfaces into two distinctive water layers. The outermost surface is ion-depleted, and the subsurface layer is ion-enriched. This surface stratification is a key element in explaining the ion-induced water reorganization at the outermost air/water interface.

Identifying Local Interfaces

While interfacial regions often occupy a relatively small portion of a system, physical and chemical processes often proceed differently within them. It is therefore useful to identify interfacial regions to answer many questions in physical chemistry. Thermodynamic phases are often described by their density and local structure; therefore, interfacial regions can then be defined as regions with densities and structures that deviate from the properties of the neighboring phases. Using this perspective of local density and structure around an atom, we describe a "directed search cone" method that has proved useful in identifying atoms that sit at the interface between two regions of a system. We call the set of atoms found to be sitting on the surface "leading atoms", and we construct an interface from these atoms that we call the "leading layer interface". We demonstrate the leading layer interface on solid-vacuum, liquid-vacuum, and liquid-vapor systems. In addition to presenting our method and example calculations, we discuss some observations of local density fluctuations that may be useful for the analysis of heterogeneous systems. Depending on the circumstances, there are various perspectives of an interface that may be insightful, and our leading layer interface will be useful in situations where the correlation between interfacial dynamics and local molecular composition is investigated.

Predicting the Geometry of Core-Shell Structures: How a Shape Changes with Constant Added Thickness

The core-shell assembly motif is ubiquitous in chemistry. While the most obvious examples are core/shell-type nanoparticles, many other examples exist. The shape of the core/shell constructs is poorly understood, making it impossible to separate chemical effects from geometric effects. Here, we create a model for the core/shell construct and develop proof for how the eccentricity is expected to change as a function of the shell. We find that the addition of a constant thickness shell always creates a relatively more spherical shape for all shapes covered by our model unless the shape is already spherical or has some underlying radial symmetry. We apply this work to simulated AOT reverse micelles and demonstrate that it is remarkably successful at explaining the observed shapes of the chemical systems. We identify the three specific cases where the model breaks down and how this impacts eccentricity.

Spiers Memorial Lecture: Water at interfaces

In this article we discuss current issues in the context of the four chosen subtopics for the meeting: dynamics and nano-rheology of interfacial water, electrified/charged aqueous interfaces, ice interfaces, and soft matter/water interfaces. We emphasize current advances in both theory and experiment, as well as important practical manifestations and areas of unresolved controversy.

Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit

The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.

Relationships between Water's Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. However, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface-solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute-surface affinity; specifically, we find that zwitterionic chemistries produce solute-surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute-surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

Computational Investigations of the Water Structure at the α-Al2O3(0001)-Water Interface

The α-Al2O3(0001)-water interface is investigated using ab initio molecular dynamics (AIMD) simulations. The spectral signatures of the vibrational sum frequency generation (vSFG) spectra of OH stretching mode for water molecules at the interface are related to the interfacial water orientation, hydrogen bond network, and water dissociation process at different water/alumina interfaces. Significant differences are found between alumina surfaces at different hydroxylation levels, namely, Al-terminated and O-terminated α-Al2O3(0001). By calculating the vibrational sum frequency generation spectrum and its imaginary component from AIMD results, the structure of interfacial waters as well as the termination of alumina slab are related to the spectral signatures of vSFG data.

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer icephilicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface-ice contact at the expense of surface-water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

Characterizing Interfaces by Voronoi Tessellation

The chemistry of interfaces differs markedly from that of the bulk. Calculation of interfacial properties depends strongly on the definition of the interface, which can lead to ambiguous results that vary between studies. There is a need for a method that can explicitly define the interfaces and boundaries in molecular systems. Voronoi tessellation offers an attractive solution to this problem through its ability to determine neighbors among specified groups of atoms. Here we discuss three cases where Voronoi tessellation combined with modeling of vibrational sum frequency generation (SFG) spectroscopy yields relevant insights: the breakdown of the air-water interface into clear and intuitive molecular layers, the study of the hydration shell in biological systems, and the acceleration of difficult spectral calculations where intermolecular vibrational couplings dominate. The utility of Voronoi tessellation has broad applications that extend beyond any single type of spectroscopy or system.

Microscopic insights on clathrate hydrate growth from non-equilibrium molecular dynamics simulations

Clathrate hydrates form and grow at interfaces. Understanding the relevant molecular processes is crucial for developing hydrate-based technologies. Many computational studies focus on hydrate growth within the aqueous phase using the 'direct coexistence method', which is limited in its ability to investigate hydrate film growth at hydrocarbon-water interfaces. To overcome this shortcoming, a new simulation setup is presented here, which allows us to study the growth of a methane hydrate nucleus in a system where oil-water, hydrate-water, and hydrate-oil interfaces are all simultaneously present, thereby mimicking experimental setups. Using this setup, hydrate growth is studied here under the influence of two additives, a polyvinylcaprolactam oligomer and sodium dodecyl sulfate, at varying concentrations. Our results confirm that hydrate films grow along the oil-water interface, in general agreement with visual experimental observations; growth, albeit slower, also occurs at the hydrate-water interface, the interface most often interrogated via simulations. The results obtained demonstrate that the additives present within curved interfaces control the solubility of methane in the aqueous phase, which correlates with hydrate growth rate. Building on our simulation insights, we suggest that by combining data for the potential of mean force profile for methane transport across the oil-water interface and for the average free energy required to perturb a flat interface, it is possible to predict the performance of additives used to control hydrate growth. These insights could be helpful to achieve optimal methane storage in hydrates, one of many applications which are attracting significant fundamental and applied interests.

Surface Affinity of Tetramethylammonium Iodide in Aqueous Solutions: A Combined Experimental and Computer Simulation Study

The surface affinity of tetramethylammonium iodide (TMAI) in aqueous solutions is investigated by surface tension measurements and molecular dynamics computer simulations. Experiments, performed in the entire composition range of solubility using the pendant drop method with two different setups, clearly reveal that TMAI is a weakly capillary active salt. Computer simulations performed with the AMBER force field reproduce the experimental data very well, while two other major force fields (i.e., CHARMM and OPLS) can still reproduce the experimental trend qualitatively; however, even qualitative reproduction of the experimental trend requires scaling down the ion charges according to the Leontyev-Stuchebrukhov correction. On the other hand, the GROMOS force field fails in reproducing the experimentally confirmed capillary activity of TMAI. Molecular dynamics simulation results show that, among the two ions, iodide has a clearly larger surface affinity than tetramethylammonium (TMA⁺). Further, the adsorption of the I^- anions is strictly limited to the first molecular layer beneath the liquid-vapor interface, which is followed by several layers of their depletion. On the other hand, the net negative charge of the surface layer, caused by the excess amount of I^- with respect to TMA⁺, is compensated by a diffuse layer of adsorbed TMA⁺ cations, extending to or beyond the fourth molecular layer beneath the liquid surface.

Local structure of liquid/vapour interfaces approaching the critical point

Investigating the structure of fluid interfaces at high temperatures is a particularly delicate task that requires effective ways of discriminating liquid from vapour and identifying the location of the liquid phase boundary, thereby allowing to distinguish intrinsic from capillary fluctuations. Several numerical approaches require introducing a coarse-graining length scale, often heuristically chosen to be the molecular size, to determine the location of the liquid phase boundary. Here, we propose an alternative rationale for choosing this coarse-graining length scale; we require the average position of the local liquid phase dividing surface to match its flat, macroscopic counterpart. We show that this approach provides additional insight into the structure of the liquid/vapour interface, suggesting the presence of another length scale beyond the bulk correlation one that plays an important role in determining the interface structure.

Molecular Dynamics Simulation Prediction of the Partitioning Constants (KH, Kiw, Kia) of 82 Legacy and Emerging Organic Contaminants at the Water-Air Interface

The tendency of organic contaminants (OCs) to partition between different phases is a key set of properties that underlie their human and ecological health impacts and the success of remediation efforts. A significant challenge associated with these efforts is the need for accurate partitioning data for an ever-expanding list of OCs and breakdown products. All-atom molecular dynamics (MD) simulations have the potential to help generate these data, but existing studies have applied these techniques only to a limited variety of OCs. Here, we use established MD simulation approaches to examine the partitioning of 82 OCs, including many compounds of critical concern, at the water-air interface. Our predictions of the Henry's law constant (K_H) and interfacial adsorption coefficients (K_iw, K_ia) correlate strongly with experimental results, indicating that MD simulations can be used to predict K_H, K_iw, and K_ia values with mean absolute deviations of 1.1, 0.3, and 0.3 logarithmic units after correcting for systematic bias, respectively. A library of MD simulation input files for the examined OCs is provided to facilitate future investigations of the partitioning of these compounds in the presence of other phases.

Ångstrom-Depth Resolution with Chemical Specificity at the Liquid-Vapor Interface

The determination of depth profiles across interfaces is of primary importance in many scientific and technological areas. Photoemission spectroscopy is in principle well suited for this purpose, yet a quantitative implementation for investigations of liquid-vapor interfaces is hindered by the lack of understanding of electron-scattering processes in liquids. Previous studies have shown, however, that core-level photoelectron angular distributions (PADs) are altered by depth-dependent elastic electron scattering and can, thus, reveal information on the depth distribution of species across the interface. Here, we explore this concept further and show that the experimental anisotropy parameter characterizing the PAD scales linearly with the average distance of atoms along the surface normal obtained by molecular dynamics simulations. This behavior can be accounted for in the low-collision-number regime. We also show that results for different atomic species can be compared on the same length scale. We demonstrate that atoms separated by about 1 Å along the surface normal can be clearly distinguished with this method, achieving excellent depth resolution.

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface─related to its hydrophobicity or tendency to repel water─can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating the relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon-66, and poly(vinyl alcohol). We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation─even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. These results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

Adsorption Studies at the Graphene Oxide-Liquid Interface: A Molecular Dynamics Study

The adsorption of organic aromatic molecules, namely aniline, onto graphene oxide is investigated using molecular simulations. The effect of the oxidation level of the graphene oxide sheet as well as the presence of two different halide salts, sodium chloride and sodium iodide, were examined. The aniline molecule in the more-reduced graphene oxide case, in the absence of added salt, showed a slightly greater affinity for the graphene oxide-water interface as compared to the oxidized form. The presence of the iodide ion increased the affinity of the aniline molecule in the reduced case but had the opposite effect for the more-oxidized form. The effect of oxidation and added salt on the interfacial water layer was also examined.

Do Machine-Learning Atomic Descriptors and Order Parameters Tell the Same Story? The Case of Liquid Water

Machine-learning (ML) has become a key workhorse in molecular simulations. Building an ML model in this context involves encoding the information on chemical environments using local atomic descriptors. In this work, we focus on the Smooth Overlap of Atomic Positions (SOAP) and their application in studying the properties of liquid water both in the bulk and at the hydrophobic air-water interface. By using a statistical test aimed at assessing the relative information content of different distance measures defined on the same data space, we investigate if these descriptors provide the same information as some of the common order parameters that are used to characterize local water structure such as hydrogen bonding, density, or tetrahedrality to name a few. Our analysis suggests that the ML description and the standard order parameters of the local water structure are not equivalent. In particular, a combination of these order parameters probing local water environments can predict SOAP similarity only approximately, and vice versa, the environments that are similar according to SOAP are not necessarily similar according to the standard order parameters. We also elucidate the role of some of the metaparameters in the SOAP definition in encoding chemical information.

Quantifying the Molecular Polarization Response of Liquid Water Interfaces at Heterogeneously Charged Surfaces

The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water's hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water's collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein.

Machine learning assisted coarse-grained molecular dynamics modeling of meso-scale interfacial fluids

A hallmark of meso-scale interfacial fluids is the multi-faceted, scale-dependent interfacial energy, which often manifests different characteristics across the molecular and continuum scale. The multi-scale nature imposes a challenge to construct reliable coarse-grained (CG) models, where the CG potential function needs to faithfully encode the many-body interactions arising from the unresolved atomistic interactions and account for the heterogeneous density distributions across the interface. We construct the CG models of both single- and two-component polymeric fluid systems based on the recently developed deep coarse-grained potential [Zhang et al., J. Chem. Phys. 149, 034101 (2018)] scheme, where each polymer molecule is modeled as a CG particle. By only using the training samples of the instantaneous force under the thermal equilibrium state, the constructed CG models can accurately reproduce both the probability density function of the void formation in bulk and the spectrum of the capillary wave across the fluid interface. More importantly, the CG models accurately predict the volume-to-area scaling transition for the apolar solvation energy, illustrating the effectiveness to probe the meso-scale collective behaviors encoded with molecular-level fidelity.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces

Uncovering microscopic hydrophilicity and hydrophobicity at heterogeneous aqueous interfaces is essential as it dictates physico/chemical properties such as wetting, the electrical double layer, and reactivity. Several molecular and spectroscopic descriptors were proposed, but a major limitation is the lack of connections between them. Here, we combine density functional theory-based MD simulations (DFT-MD) and SFG spectroscopy to explore how interfacial water responds in contact with self-assembled monolayers (SAM) of tunable hydrophilicity. We introduce a microscopic metric to track the transition from hydrophobic to hydrophilic interfaces. This metric combines the H/V descriptor, a structural descriptor based on the preferential orientation within the water network in the topmost binding interfacial layer (BIL) and spectroscopic fingerprints of H-bonded and dangling OH groups of water carried by BIL-resolved SFG spectra. This metric builds a bridge between molecular descriptors of hydrophilicity/hydrophobicity and spectroscopically measured quantities and provides a recipe to quantitatively or qualitatively interpret experimental SFG signals.

Water dynamics around T0 vs R4 of hemoglobin from local hydrophobicity analysis

The local hydration around tetrameric hemoglobin (Hb) in its T₀ and R₄ conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α₁β₂ and α₂β₁ interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod-Wyman-Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas "buried surface" reports on the available space for solvent. For simulations with Hb frozen in its T₀ and R₄ states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41α, Tyr42α, Tyr140α, Trp37β, Glu101β (for T₀) and Thr38α, Tyr42α, Tyr140α (for R₄). The number of water molecules at the interface is found to increase by ∼25% for T₀ → R₄, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.

Highly Altered State of Proton Transport in Acid Pools in Charged Reverse Micelles

Transport mechanisms of solvated protons of 1 M HCl acid pools, confined within reverse micelles (RMs) containing the negatively charged surfactant sodium bis(2-ethylhexyl) sulfosuccinate (NaAOT) or the positively charged cetyltrimethylammonium bromide (CTABr), are analyzed with reactive force field simulations to interpret dynamical signatures from TeraHertz absorption and dielectric relaxation spectroscopy. We find that the forward proton hopping events for NaAOT are further suppressed compared to a nonionic RM, while the Grotthuss mechanism ceases altogether for CTABr. We attribute the sluggish proton dynamics for both charged RMs as due to headgroup and counterion charges that expel hydronium and chloride ions from the interface and into the bulk interior, thereby increasing the pH of the acid pools relative to the nonionic RM. For charged NaAOT and CTABr RMs, the localization of hydronium near a counterion or conjugate base reduces the Eigen and Zundel configurations that enable forward hopping. Thus, localized oscillatory hopping dominates, an effect that is most extreme for CTABr in which the proton residence time increases dramatically such that even oscillatory hopping is slow.

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.

Vibrational Sum Frequency Generation Spectra of Water-Vapor Interfaces Covered by Alcohols: Effects of Surface Coverage and Coupling between Oscillators

The present study deals with the effects of varying coverage of water surface by alcohols on the vibrational sum frequency generation (VSFG) spectrum of interfacial water. We have considered two different alcohols: Tertiary butyl alcohol (TBA) whose alkyl part is fully branched and stearyl alcohol (STA) which has a long linear alkyl chain with larger hydrophobic surface area than that of TBA. With increase of the alcohol concentration, the hydrogen bonded OH stretch region of the VSFG spectrum is found to change following a regular trend for the STA-water system, whereas non-monotonic variation of the VSFG spectrum is observed for the TBA-water system which can be correlated with the presence of very different interactions of TBA molecules at different concentrations. On increasing the concentration of TBA, the hydrophobic groups get more tilted towards the water phase and significant hydrophobic interactions are introduced at higher concentrations. Whereas, for STA, there is a gradual increase in the hydrophilic interaction. Because of stacking interactions between the long chain alkyl groups, the hydrophobic parts stay outward from the water phase at higher concentrations and a regular change in the VSFG spectrum is observed. We have also presented a computationally efficient scheme to calculate the VSFG spectrum of interfacial systems for coupled oscillators which is expected to be beneficial for the treatment of coupling where the interfacial system size is inherently large.

Spontaneous Formation of Hydrogen Peroxide in Water Microdroplets

There is accumulating evidence that many chemical reactions are accelerated by several orders of magnitude in micrometer-sized aqueous or organic liquid droplets compared to their corresponding bulk liquid phase. However, the molecular origin of the enhanced rates remains unclear as in the case of spontaneous appearance of 1 μM hydrogen peroxide in water microdroplets. In this Letter, we consider the range of ionization energies and whether interfacial electric fields of a microdroplet can feasibly overcome the high energy step from hydroxide ions (OH^-) to hydroxyl radicals (OH^•) in a primary H₂O₂ mechanism. We find that the vertical ionization energies (VIEs) of partially solvated OH^- ions are greatly lowered relative to the average VIE in the bulk liquid, unlike the case of the Cl^- anion which shows no reduction in the VIEs regardless of solvation environment. Overall reduced hydrogen-bonding and undercoordination of OH^- are structural features that are more readily present at the air-water interface, where the energy scale for ionization can be matched by statistically probable electric field values.

Analysis of Gas Nanoclusters in Water Using All-Atom Molecular Dynamics

The Young-Laplace equation suggests that nanosized gas clusters would dissolve under the effects of perturbation. The fact that nanobubbles are observed raises questions as to the mechanism underlying their stability. In the current study, we used all-atom molecular dynamics simulations to investigate the gas-water interfacial properties of gas clusters. We employed the instantaneous coarse-graining method to define the fluctuating boundaries and analyze the deformation of gas clusters. Fourier transform analysis of the cluster morphology revealed that the radius and morphology deformation variations exhibit power law relationships with the vibrational frequency, indicating that the surface energy dissipated through morphology variations. Increasing pressure in the liquid region was found to alter the network of water molecules at the interface, whereas increasing pressure in the gas region did not exhibit this effect. The overall gas concentration was oversaturated and proportional to the gas density inside the clusters. However, the result of comparison with Henry's law reveals that the gas pressure at the interface reduced by the interfacial effects is much lower than that inside the gas region, thus reducing the demanding degree of oversaturation. Originating from the interfacial charge allocation, the magnitude of the electrostatic stress is greater than that of the gas pressure inside the cluster. However, the magnitude of the reversed tension induced by electrostatic stress is far below the value of interfacial tension. The potential of mean force (PMF) profiles revealed that a barrier potential at the interface hindered gas particles from escaping the cluster. Several effects contribute to stabilizing the gas clusters in water, including high-frequency morphological deformation, electrostatic stress, reduced interfacial tension, and gas oversaturation conditions. Our results suggest that gas clusters can exist in water under gas oversaturation conditions in the absence of hydrophobic contaminants or pinning charges at interfaces.

Anomalous properties in the potential energy landscape of a monatomic liquid across the liquid-gas and liquid-liquid phase transitions

As a liquid approaches the gas state, the properties of the potential energy landscape (PEL) sampled by the system become anomalous. Specifically, (i) the mechanically stable local minima of the PEL [inherent structures (IS)] can exhibit cavitation above the so-called Sastry volume, vS, before the liquid enters the gas phase. In addition, (ii) the pressure of the liquid at the sampled IS [i.e., the PEL equation of state, PIS(v)] develops a spinodal-like minimum at vS. We perform molecular dynamics simulations of a monatomic water-like liquid and verify that points (i) and (ii) hold at high temperatures. However, at low temperatures, cavitation in the liquid and the corresponding IS occurs simultaneously and a Sastry volume cannot be defined. Remarkably, at intermediate/high temperatures, the IS of the liquid can exhibit crystallization, i.e., the liquid regularly visits the regions of the PEL that belong to the crystal state. The model liquid considered also exhibits a liquid-liquid phase transition (LLPT) between a low-density and a high-density liquid (LDL and HDL). By studying the behavior of PIS(v) during the LLPT, we identify a Sastry volume for both LDL and HDL. The HDL Sastry volume marks the onset above which IS are heterogeneous (composed of LDL and HDL particles), analogous to points (i) and (ii) above. However, the relationship between the LDL Sastry volume and the onset of heterogeneous IS is less evident. We conclude by presenting a thermodynamic argument that can explain the behavior of the PEL equation of state PIS(v) across both the liquid-gas phase transition and LLPT.

Connecting Non-Gaussian Water Density Fluctuations to the Lengthscale Dependent Crossover in Hydrophobic Hydration

We connect density fluctuations in liquid water to lengthscale dependent crossover in hydrophobic hydration. Specifically, we employ indirect umbrella sampling (INDUS) simulations to characterize density fluctuations in observation volumes of various sizes and shapes in water and as a function of temperature and salt concentration. Consistent with previous observations, density fluctuations are Gaussian in small molecular scale volumes, but they display non-Gaussian "low-density fat tails" in larger volumes. These non-Gaussian tails are indicative of the proximity of water to its liquid to vapor phase transition and have implications on biomolecular interactions and function. We show that the onset of non-Gaussian fluctuations in large volumes is accompanied by the formation of a cavity in the observation volume. We develop a model that uses the physics of cavity-water interface formation as a key ingredient and show that it captures the nature of non-Gaussian density fluctuations over a broad region in water and in salt solutions. We discuss the limitations of this model in the very low density region of the distribution. Our calculations provide new insights into the origins of non-Gaussian density fluctuations in water and their connections to lengthscale dependent crossover in hydrophobic hydration.

Individual adsorption of low volatility pheromones: Amphiphilic molecules on a clean water-air interface

Environmental conditions can alter olfactory scent and chemical communication among biological species. In particular, odorant molecules interact with aerosols. Thermodynamics variables governing the adsorption from air to water surface of bombykol, the most studied pheromone, and of three derivative molecules, bombykal, bombykoic acid, and bombykyle acetate, are computed by steered and un-biased molecular dynamics in order to compare the role of their polar head group on adsorption on aqueous aerosols. When adsorbed, the molecule center of mass stands at about 1.2 Å from the interface and oscillates on the same length scale, trapped in an energy well. Gibbs energy of adsorption and desorption time of bombykol are found to be 9.2 k_BT and 59 µs, respectively. The following ordering between the molecules is observed, reading from the more to the least adsorbed: bombykoic acid > bombykol > bombykoic acetate > bombykal. It originates from a complex interplay of entropy and enthalpy. The entropy and enthalpy of adsorption are discussed in the light of structural arrangement, H-bonding, and hydrophilic tail positioning of the molecules at the interface. Our results show that, when dispersed in the air, pheromones adsorb on aqueous aerosols. However, the individual residence time is quite short on pure water surfaces. Aerosols can, therefore, only have a decisive influence on chemical communication through collective effects or through their chemical composition that is generally more complex than that of a pure water surface.

The dielectric function profile across the water interface through surface-specific vibrational spectroscopy and simulations

The dielectric properties of interfacial water on subnanometer length scales govern chemical reactions, carrier transfer, and ion transport at interfaces. Yet, the nature of the interfacial dielectric function has remained under debate as it is challenging to access the interfacial dielectric with subnanometer resolution. Here we use the vibrational response of interfacial water molecules probed using surface-specific sum-frequency generation (SFG) spectra to obtain exquisite depth resolution. Different responses originate from water molecules at different depths and report back on the local interfacial dielectric environment via their spectral amplitudes. From experimental and simulated SFG spectra at the air/water interface, we find that the interfacial dielectric constant changes drastically across an ∼1 Å thin interfacial water region. The strong gradient of the interfacial dielectric constant leads, at charged planar interfaces, to the formation of an electric triple layer that goes beyond the standard double-layer model.

Structure, Thermodynamics, and Dynamics of Thiocyanate Ion Adsorption and Transfer across the Water/Toluene Interface

Molecular dynamics simulations are used to examine in detail the structure, thermodynamics, and dynamics involve in the adsorption and transfer of the thiocyanate ion (SCN^-) across the water/toluene interface. Free energy, hydration structure, and several dynamical properties as a function of the ion location along the interface normal are calculated and contrasted with recent experiments. The free energy profile exhibits a local minimum near the interface corresponding to adsorption free energy relative to bulk water of -6 kJ/mol, in reasonable agreement with experiments. The simulations provide insight into the water surface fluctuations that are coupled to the ion transfer, demonstrating formation of water finger-like structures assisting the transfer process.

Local Ice-like Structure at the Liquid Water Surface

Experiments and computer simulations have established that liquid water's surfaces can deviate in important ways from familiar bulk behavior. Even in the simplest case of an air-water interface, distinctive layering, orientational biases, and hydrogen bond arrangements have been reported, but an overarching picture of their origins and relationships has been incomplete. Here we show that a broad set of such observations can be understood through an analogy with the basal face of crystalline ice. Using simulations, we demonstrate a number of structural similarities between water and ice surfaces, suggesting the presence of domains at the air-water interface with ice-like features that persist over 2-3 molecular diameters. Most prominent is a shared characteristic layering of molecular density and orientation perpendicular to the interface. Lateral correlations of hydrogen bond network geometry point to structural similarities in the parallel direction as well. Our results bolster and significantly extend previous conceptions of ice-like structure at the liquid's boundary and suggest that the much-discussed quasi-liquid layer on ice evolves subtly above the melting point into a quasi-ice layer at the surface of liquid water.

Polarization-Dependent Sum-Frequency Generation Spectroscopy for Ångstrom-Scale Depth Profiling of Molecules at Interfaces

The three-dimensional spatial distribution of molecules at soft matter interfaces is crucial for processes ranging from membrane biophysics to atmospheric chemistry. While several techniques can access surface composition, obtaining information on the depth distribution is challenging. We develop a noninvasive, polarization-resolved, surface-specific sum-frequency generation spectroscopy providing quantitative depth information. We demonstrate the technique on formic acid molecules at the air-water interface. With increasing molar fraction from 2.5% to 10%, the formic acid molecules shift, on average, ∼0.9  Å into the bulk. The consistency with the simulation data manifests that the technique allows for probing the Ångstrom-scale depth profile.

Hydration of Linear Alkanes is Governed by the Small Length-Scale Hydrophobic Effect

Length-scale dependence of the hydrophobic effect is well understood for apolar spherical solutes: for small solutes (diameter, d ≲ 0.8 nm), the hydration free energy is entropically driven, while for larger solutes (d ≳ 2 nm), it is enthalpically driven. The nature of the hydrophobic effect in the case of anisotropic molecules such as linear alkanes is not understood yet. In this work, we have calculated the hydration free energy of linear alkanes going from methane to octadecane and of a spherical decane droplet of d ≈ 3 nm using molecular simulations. We show that the hydration free energies of alkanes, irrespective of their size, are governed by the small length-scale hydrophobic effect. That is, unlike the case of large spherical solutes, the hydration free energies of linear alkanes are entropically driven.

The Effects of Flexibility on dsDNA–dsDNA Interactions

A detailed understanding of the physical mechanism of ion-mediated dsDNA interactions is important in biological functions such as DNA packaging and homologous pairing. We report the potential of mean force (PMF) or the effective solvent mediated interactions between two parallel identical dsDNAs as a function of interhelical separation in 0.15 M NaCl solution. Here, we study the influence of flexibility of dsDNAs on the effective interactions by comparing PMFs between rigid models and flexible ones. The role of flexibility of dsDNA pairs in their association is elucidated by studying the energetic properties of Na⁺ ions as well as the fluctuations of ions around dsDNAs. The introduction of flexibility of dsDNAs softens the vdW contact wall and induces more counterion fluctuations around dsDNAs. In addition, flexibility facilitates the Na⁺ ions dynamics affecting their distribution. The results quantify the extent of attraction influenced by dsDNA flexibility and further emphasize the importance of non-continuum solvation approaches.

The Intrinsic Fragility of the Liquid-Vapor Interface: A Stress Network Perspective

The evolution of the liquid-vapor interface of a Lennard-Jones fluid is examined with molecular dynamics simulations using the intrinsic sampling method. Results suggest clear damping of the intrinsic profiles with increasing temperature. Investigating the surface stress distribution, we have identified a linear variation of the space-filling nature (fractal dimension) of the stress clusters at the intrinsic surface with increasing surface tension or, equivalently, with decreasing temperature. A percolation analysis of these stress networks indicates that the stress field is more disjointed at higher temperatures. This leads to more fragile (or poorly connected) interfaces which result in a reduction in surface tension.

How do interfaces alter the dynamics of supercooled water?

The structure of liquid water in the proximity of an interface can deviate significantly from that of bulk water, with surface-induced structural perturbations typically converging to bulk values at about ∼1 nm from the interface. While these structural changes are well established it is, in contrast, less clear how an interface perturbs the dynamics of water molecules within the liquid. Here, through an extensive set of molecular dynamics simulations of supercooled bulk and interfacial water films and nano-droplets, we observe the formation of persistent, spatially extended dynamical domains in which the average mobility varies as a function of the distance from the interface. This is in stark contrast with the dynamical heterogeneity observed in bulk water, where these domains average out spatially over time. We also find that the dynamical response of water to an interface depends critically on the nature of the interface and on the choice of interface definition. Overall these results reveal a richness in the dynamics of interfacial water that opens up the prospect of tuning the dynamical response of water through specific modifications of the interface structure or confining material.