Hydrogen Bond Dynamics at the Water/Hydrocarbon Interface

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

Micro-Solvation of Propofol in Propylene Glycol-Water Binary Mixtures: Molecular Dynamics Simulation Studies

The water microstructure around propofol plays a crucial role in controlling their solubility in the binary mixture. The unusual nature of such a water microstructure can influence both translational and reorientational dynamics, as well as the water hydrogen bond network near propofol. We have carried out all-atom molecular dynamics simulations of five different compositions of the propylene glycol (PG)/water binary mixture containing propofol (PFL) molecules to investigate the differential behavior of water microsolvation shells around propofol, which is likely to control the propofol solubility. It is evident from the simulation snapshots for various compositions that the PG at high molecular ratio favors the water cluster and extended chainlike network that percolates within the PG matrix, where the propofol is in the dispersed state. We estimated that the radial distribution function indicates higher ordered water microstructure around propofol for high PG content, as compared to the lower PG content in the PG/water mixture. So, the hydrophilic PG regulates the stability of the water micronetwork around propofol and its solubility in the binary mixture. We observed that the translational and rotational mobility of water belonging to the propofol microsolvation shell is hindered for high PG content and relaxed toward the low PG molecular ratio in the PG/water mixture. It has been noticed that the structural relaxation of the hydrogen bond formed between the propofol and the water molecules present in the propofol microsolvation shell for all five compositions is found to be slower for high PG content and becomes faster on the way to low PG content in the mixture. Simultaneously, we calculated the intermittent residence time correlation function of the water molecules belonging to the microsolvation shell around the propofol for five different compositions and found a faster short time decay followed up with long time components. Again, the origin of such long time decay is primarily from the structural relaxation of the microsolvation shell around the propofol, where the high PG content shows the slower structural relaxation that turns faster as the PG content approaches to the other end of the compositions. So, our studies showed that the slower structural relaxation of the microsolvation shell around propofol for a high PG molecular ratio in the PG/water mixture correlate well with the extensive ordering of the water microstructure and restricted water mobility and facilitates the dissolution process of propofol in the binary mixture.

Resonance-Induced Reduction of Interfacial Tension of Water-Methane and Improvement of Methane Solubility in Water

Molecular dynamics simulations were performed to study the effect of the periodic oscillating electric field on the interface between water and methane. We propose a new strategy that utilizes oscillating electric fields to reduce the interfacial tension (IFT) between water and methane and increase the solubility of methane in water simultaneously. These are attributed to the hydrogen bond resonance induced by an electric field with a frequency close to the natural frequency of the hydrogen bond. The resonance breaks the hydrogen bond network among water molecules to the maximum, which destroys the hydration shell and reduces the cohesive action of water, thus resulting in the decrease of IFT and the increase of methane solubility. As the frequency of the electric field is close to the optimum resonant frequency of hydrogen bonds, IFT decreases from 56.43 to 5.66 mN/m; water and methane are miscible because the solubility parameter of water reduces from 47.63 to 2.85 MPa^1/2, which is close to that of methane (3.43 MPa^1/2). Our results provide a new idea for reducing the water-gas IFT and improving the solubility of insoluble gas in water and theoretical guidance in the fields of natural gas exploitation, hydrate generation, and nanobubble nucleation.

Rotational Dynamics of Water at the Phospholipid Bilayer Depending on the Head Groups Studied by Molecular Dynamics Simulations

Hydration states are a crucial factor that affect the self-assembly and properties of soft materials and biomolecules. Although previous experiments have revealed that the hydration state strongly depends on the chemical structure of lipid molecules, the mechanisms at the molecular level remain unknown. Classical and density-functional tight-binding (DFTB) molecular dynamics (MD) simulations are employed to determine the mechanisms underlying dissimilar water dynamics between lipid membranes with phosphatidylcholine (PC) and phosphatidylethanolamine (PE) head groups. The classical MD simulation shows that rotational relaxations of water are faster on the PE lipid than on the PC lipid, which is consistent with a previous experimental study using terahertz spectroscopy. Furthermore, DFTB-MD simulation of N(CH₃)₄⁺ and NH₄⁺ ions, which correspond to the different head groups in PC and PE, respectively, shows qualitative agreement with the classical MD simulation. Remarkably, the PE lipids and the NH₄⁺ ions break hydrogen bonds between water molecules in the secondary hydration shell. In contrast, the PC lipids and the N(CH₃)₄⁺ ions bind water molecules weakly in both the primary and secondary hydration shells and increase hydrogen bonds between water. Our atomistic simulations show that these changes in the hydrogen-bond network of water molecules cause the different rotational relaxation of water molecules between the two lipids.

Interfacial structure in the liquid-liquid extraction of rare earth elements by phosphoric acid ligands: a molecular dynamics study

Solvent extraction (SX), wherein two immiscible liquids, one containing the extractant molecules and the other containing the solute to be extracted are brought in contact to effect the phase transfer of the solute, underpins metal extraction and recovery processes. The interfacial region is of utmost importance in the SX process, since besides thermodynamics, the physical and chemical heterogeneity at the interface governs the kinetics of the process. Yet, a fundamental understanding of this heterogeneity and its implications for the extraction mechanism are currently lacking. We use molecular dynamics (MD) simulations to study the liquid-liquid interface under conditions relevant to the SX of Rare Earth Elements (REEs) by a phosphoric acid ligand. Simulations revealed that the extractant molecules and varying amounts of acid and metal ions partitioned to the interface. The presence of these species had a significant effect on the interfacial thickness, hydrogen bond life times and orientations of the water molecules at the interface. Deprotonation of the ligands was essential for the adsorption of the metal ions at the interface, with these ions forming a number of different complexes at the interface involving one to three extractant molecules and four to eight water molecules. Although the interface itself was rough, no obvious 'finger-like' water protrusions penetrating the organic phase were seen in our simulations. While the results of our work help us gain fundamental insights into the sequence of events leading to the formation of a variety of interfacial complexes, they also emphasize the need to carry out a more detailed atomic level study to understand the full mechanism of extraction of REEs from the aqueous to organic phases by phosphoric acid ligands.

Tuning Solvent Miscibility: A Fundamental Assessment on the Example of Induced Methanol/ n-Dodecane Phase Separation

In this work, we assess the fundamental aspects of mutual miscibility of solvents by studying the mixing of two potential candidates, methanol and n-dodecane, for nonaqueous solvent extraction. To do so, ¹H NMR spectroscopy and molecular dynamics simulations are used jointly. The NMR spectra show that good phase separation can be obtained by adding LiCl and that the addition of a popular extractant (tri-n-butyl phosphate) yields the opposite effect. It is also demonstrated that in a specific case the poor phase separation is not due to the migration of n-dodecane into the more polar phase, but due to the transfer of the extractant into it, which is especially relevant when considering industrial applications of solvent extraction. With the aid of molecular dynamics simulations, explanations of this behavior are given. Specifically, an increase of all hydrogen-bond lifetimes is found to be consequent to the addition of LiCl which implies an indirect influence on the methanol liquid structure, by favoring a stronger hydrogen-bond network. Therefore, we found that better phase separation is not directly due to the presence of LiCl, but due to the “hardening” of the hydrogen-bond network.

A general method for determining molecular interfaces and layers

A general and direct computational scheme to locate the surface separating arbitrarily shaped domains made up of molecules (or any other particles) has been developed and is described and illustrated for several, both artificial and physical examples. The proposed scheme consists of two modules: (i) triangulation and (ii) assignment of simplices to domains. Three different triangulation methods are employed, viz., the Delaunay triangulation, regular triangulation, and quasi-triangulation. In the triangulated system, the assignment step is carried out in two different ways, one based on the characteristic metric of a particular triangulation procedure and the other on the concept of a touching sphere. Some of the combinations of the triangulation and assignment steps lead to methods already used by others to find interfacial or surface molecules, namely the alpha-shape-based method of Usabiaga nad Duque [Phys. Rev. E 79 (2009) 046709] and GITIM of Sega et al. [J. Chem. Phys. 138 (2013) 044110]. The resulting surface is defined not only as a discrete set of particles, but it is build up of facets of the triangulation forming a broken line in two dimensions or a polyhedral surface in three dimensions. Individual molecular layers are identified in a very straightforward manner, starting with the interfacial layer itself and proceeding into the interior of the phase. The proposed scheme is illustrated first by identifying border molecules of pre-sampled domains of several shapes in a plane and then applied to five physically meaningful examples: thin films, near critical water, liquid water slab in an electric field, liquid water at a solid wall, and water at condition of electric-field-induced jetting. Performance of the considered methods is critically assessed. Treatment of domains forming percolating clusters through periodic boundary conditions is also described along with the determination of their periodicity and dimensionality.

Structural and thermodynamic properties of water-membrane interphases: significance for peptide/membrane interactions

Water appears as a common intermediary in the mechanisms of interaction of proteins and polypeptides with membranes of different lipid composition. In this review, how water modulates the interaction of peptides and proteins with lipid membranes is discussed by correlating the thermodynamic response and the structural changes of water at the membrane interphases. The thermodynamic properties of the lipid-protein interaction are governed by changes in the water activity of monolayers of different lipid composition according to the lateral surface pressure. In this context, different water populations can be characterized below and above the phase transition temperature in relation to the CH₂ conformers' states in the acyl chains. According to water species present at the interphase, lipid membrane acts as a water state regulator, which determines the interfacial water domains in the surface. It is proposed that those domains are formed by the contact between lipids themselves and between lipids and the water phase, which are needed to trigger adsorption-insertion processes. The water domains are essential to maintain functional dynamical properties and are formed by water beyond the hydration shell of the lipid head groups. These confined water domains probably carries information in local units in relation to the lipid composition thus accounting for the link between lipidomics and aquaomics. The analysis of these results contributes to a new insight of the lipid bilayer as a non-autonomous, responsive (reactive) structure that correlates with the dynamical properties of a living system.

Molecular dynamics simulation of water/BHDC cationic reverse micelles. structural characterization, dynamical properties, and influence of solvent on intermicellar interactions

We report results obtained from molecular dynamics (MD) experiments of benzylhexadecyldimethylammonium chloride (BHDC) cationic reverse micelles (RMs). In particular we analyzed equilibrium and dynamical characteristics of water/BHDC RMs in pure benzene, at two different water/BHDC ratios (W0 = 5 and W0 = 10). The RMs appear as elliptical aggregates with eccentricities close to ∼0.9. Analysis of the different spatial correlations reveals three different spatial domains in the RMs: a water inner pool, the surfactant interface, and the external solvent. The calculated accessible surface areas for the aqueous inner cores suggest a strong penetration of solvent molecules within the micellar interface domains. Comparison between the density profiles of both RMs shows an increment of the broadness in the distributions of all species at the interface, along with an increasing overlap between the tail segments of the surfactant and benzene molecules as one considers larger micelles. For the dynamical side, the rotational characteristic time scale for the confined water was found to be 1 order of magnitude larger than that of the bulk water. A similar effect was also observed for hydrogen bond dynamics. Both retardation effects diminish with the size of the aggregate. To the estimate the influence of the external solvent on the intermicellar interactions, free energy profiles for the coalescence process between RMs of similar size in pure benzene and in a n-heptane/benzene mixture were also investigated. The results indicate that the association process is facilitated by the presence of n-heptane in the external nonpolar phase. Comparison with previous theoretical and experimental results is also carried out.

Temperature dependence of hydrophobic hydration dynamics: from retardation to acceleration

The perturbation induced by a hydrophobic solute on water dynamics is essential in many biochemical processes, but its mechanism and magnitude are still debated. A stringent test of the different proposed pictures is provided by recent NMR measurements by Qvist and Halle (J. Am. Chem. Soc. 2008, 130, 10345-10353) which showed that, unexpectedly, the perturbation changes in a non-monotonic fashion when the solution is cooled below room temperature. Here we perform and analyze molecular dynamics simulations of a small paradigm amphiphilic solute, trimethylamine N-oxide (TMAO), in dilute aqueous solutions over the 218-350 K temperature range. We first show that our simulations properly reproduce the non-monotonic temperature dependence. We then develop a model which combines our previously suggested entropic excluded-volume effect with a perturbation factor arising from the difference between local structural fluctuations in the shell and the bulk. Our model provides a detailed molecular understanding of the hydrophobic perturbation over the full temperature range investigated. It shows that the excluded-volume factor brings a dominant temperature-independent contribution to the perturbation at all temperatures, and provides a very good approximation at room temperature. The non-monotonic temperature dependence of the perturbation is shown to arise from the structural factor and mostly from relative shifts between the shell and bulk distributions of local structures, whose amplitude remains very small compared to the widths of those distributions.

Connected and isolated CH2 populations in acyl chains and its relation to pockets of confined water in lipid membranes as observed by FTIR spectrometry

Analysis of the band corresponding to the frequency of vibrational symmetric stretching mode of methylene groups in the lipid acyl chains and the bands of water below and above the phase transition of different lipids by Fourier transform infrared spectroscopy gives strong support to the formation of confined water pockets in between the lipid acyl chains. Our measures and analysis consolidate the mechanism early proposed by Traüble, in the sense that water is present in kinks formed by trans-gauche isomers along the hydrocarbon tails. The formation of these regions depends on the acyl lipid composition, which determines the presence of different populations of water species, characterized by its degree of H bond coordination in fluid saturated or unsaturated lipids. The free energy excess due to the reinforcement of the water structure along few water molecules in the adjacencies of exposed membrane residues near the phase transition is a reasonable base to explain the insertion and translocation of polar peptides and amino acid residues through the biomembrane on thermodynamic and structural grounds.

Water jump reorientation: from theoretical prediction to experimental observation

Liquid water is remarkably labile in reorganizing its hydrogen-bond (HB) network through the breaking and forming of HBs. This rapid restructuring, which occurs on the picosecond time scale, is critical not only for many of the pure liquid's special features but also for a range of aqueous media phenomena, including chemical reactions and protein activity. An essential part of the HB network reorganization is water molecule reorientation, which has long been described as Debye rotational diffusion characterized by very small angular displacements. Recent theoretical work, however, has presented a starkly contrasting picture: a sudden, large-amplitude jump mechanism, in which the reorienting water molecule rapidly exchanges HB partners in what amounts to an activated chemical reaction. In this Account, we first briefly review the jump mechanism and then discuss how it is supported by a series of experiments. These studies range from indirect indications to direct characterization of the jumps through pioneering two-dimensional infrared spectroscopy (2D-IR), the power of which accords it a special focus here. The scenarios in which experimental signatures of the jump mechanism are sought increase in complexity throughout the Account, beginning with pure water. Here 2D-IR in combination with theory can give a glimpse of the jumps, but the tell-tale markers are not pronounced. A more fruitful arena is provided by aqueous ionic solutions. The difference between water-water and water-anion HB strengths provides the experimental handle of differing OH stretch frequencies; in favorable cases, the kinetic exchange of a water between these two sites can be monitored. Sole observation of this exchange, however, is insufficient to establish the jump mechanism. Fortunately, 2D-IR with polarized pulses has demonstrated that HB exchange is accompanied by significant angular displacement, with an estimated jump angle similar to theoretical estimates. The Janus-like character of amphiphilic solutes, with their hydrophobic and hydrophilic faces, presents a special challenge for theory and experiment. Here a consensus on the 2D-IR interpretation has not yet been achieved; this lack of accord impedes the understanding of, for example, biochemical solutes and interfaces. We argue that the influence of hydrophobic groups on water jumps is only modest and well accounted for by an excluded volume effect in the HB exchange process. Conversely, hydrophilic groups have an important influence when their HB strength with water differs significantly from that of the water-water HB. The power of 2D-IR is argued to be accompanied by subtleties that can lead to just the opposite and, in our view, erroneous conclusion. We close with a prediction that a hydrophobic surface offers an arena in which the dynamics of "dangling" water OHs, bereft of a HB, could provide a 2D-IR confirmation of water jumps.

Glassy nature of water in an ultraconfining disordered material: the case of calcium-silicate-hydrate

We present the structural and dynamic nature of water ultraconfined in the quasi-two-dimensional nanopores of the highly disordered calcium-silicate-hydrate (C-S-H), the major binding phase in cement. Our approach is based on classical molecular simulations. We demonstrate that the C-S-H nanopore space is hydrophilic, particularly because of the nonbridging oxygen atoms on the disordered silicate chains which serve as hydrogen-bond acceptor sites, directionally orienting the hydrogen atoms of the interfacial water molecules toward the calcium-silicate layers. The water in this interlayer space adopts a unique multirange structure: a distorted tetrahedral coordination at short range up to 2.7 Å, a disordered structure similar to that of dense fluids and supercooled phases at intermediate range up to 4.2 Å, and persisting spatial correlations through dipole-dipole interactions up to 10 Å. A three-stage dynamics governs the mean square displacement (MSD) of water molecules, with a clear cage stage characteristic of the dynamics in supercooled liquids and glasses, consistent with its intermediate-range structure identified here. At the intermediate time scales corresponding to the β-relaxation of glassy materials, coincident with the cage stage in MSD, the non-Gaussian parameter indicates a significant heterogeneity in the translational dynamics. This dynamic heterogeneity is induced primarily because of the heterogeneity in the distribution of hydrogen bond strengths. The strongly attractive interactions of water molecules with the calcium silicate walls serve to constrain their motion. Our findings have important implications on describing the cohesion and mechanical behavior of cement from its setting to its aging.

Influence of heterogeneity of confined water on photophysical behavior of acridine with amines: a time-resolved fluorescence and laser flash photolysis study

The photophysical behavior of acridine (Acr) shows facilitated water-assisted protonation equilibrium between its deprotonted (Acr* ∼ 10 ns) and protonated forms (AcrH(+*) ∼ 28 ns) within confined region of ordered water molecules inside AOT/H(2)O/n-heptane reverse micelles (RMs). The time-resolved-area-normalized-emission spectra confirm both Acr* and AcrH(+*), while time-resolved-emission spectra depict time evolution between them. Quenching of AcrH(+*) with N,N-dimethylaniline (DMA) is a purely diffusion-controlled bimolecular quenching with linear Stern-Volmer (S-V) plot, while nonlinearity arises with triethylamine (TEA) that forms ground state complex with AcrH(+) (AcrH(+)··H(2)O··TEA) indicating both static and dynamic quenching. Transient intermediates, DMA(•+) and AcrH(•) infer photoinduced electron transfer from DMA to Acr, while those from AcrH(+)··H(2)O··TEA complex suggest water mediated excited-state proton transfer (ESPT) between AcrH(+) and TEA. The ESPT becomes faster in larger RMs due to enhanced mobility of hydronium ions in AcrH(+)··H(2)O··TEA, which reduces in smaller RMs as water becomes much more constrained owing to stronger complexation by excess confinement.

Hydrogen-bond structure at the interfaces between water/poly(methyl methacrylate), water/poly(methacrylic acid), and water/poly(2-aminoethylmethacrylamide)

The molecular dynamics approach was employed to study the structural characteristics at the interface of water/poly(methyl methacrylate) (PMMA), water/poly(methacrylic acid) (PMAA), and poly(2-aminoethylmethacrylamide) (PAEMA). It is found that the water on the PMAA surface shows a significant increase in the density at the interface, with a greater number of water molecules permeating into the bulk of the substrate region. The structure of hydrogen bonds of water and the radial distribution function for given polar atoms in the polymer substrate are calculated. We found that a network structure of hydrogen bonding between water and the polar atom of the polymer forms at the interface. PMAA exhibits a more hydrophilic property than PMMA and PAEMA because it generates a shell-like structure of water molecules around its functional group. Finally, the hydrogen bond numbers of PMMA, PMAA, and PAEMA are also analyzed. The results detail the hydrogen bond structure of each specific atom and find that, in all three cases, the carboxyl oxygen attracts the greatest number of water molecules compared with other atoms.

Role of electrostatics in modulating hydrophobic interactions and barriers to hydrophobic assembly

Hydrophobic effects continue to be an active area of research due to implications for a wide range of physicochemical phenomena. Molecular dynamics simulations have been used extensively in the study of such effects using various water potential models, with few studies addressing the differences between models. In particular, studies considering the explicit treatment of water polarizability are underrepresented in the literature. We present results from molecular dynamics simulations that systematically compare the dependence of large-scale hydrophobic effects on the water model. We consider three common nonpolarizable models (SPC/E, TIP3P, and TIP4P) and two common polarizable models (TIP4P-FQ and SWM4-NDP). Results highlight the similarities and differences of the different water models in the vicinity of two large hydrophobic plates. In particular, profiles of average density, density fluctuations, orientation, and hydrogen bonding show only minor differences among the water models studied. However, the potential of mean force for the hydrophobe dimerization is significantly reduced in the polarizable water systems. TIP4P-FQ shows the deepest minimum of approximately -54(+/-3) kcal/mol compared to -40(+/-3), -40(+/-2), -42(+/-3), and -45(+/-5) kcal/mol for TIP4P, TIP3P, SPC/E, and SWM4-NDP (all relative to the dissociated state). We discuss the relationship between hydrophobic association and the strength of water-water interactions in the liquid phase. Results suggest that models treating polarizability (both implicitly and explicitly) influence a stronger driving force toward hydrophobic assembly. Implications of these results, as well as prospectives on future work, are discussed.