Hydrogen Bond Dynamics at Water/Organic Liquid Interfaces

Design of deep eutectic solvents for multiple perfluoroalkyl substances removal: Energy-based screening and mechanism elucidation

The current landscape of perfluoroalkyl substances (PFAS) extraction methodologies presents significant challenges, particularly for multiple PFAS with different carbon chain lengths. This study introduced an energy-driven strategic approach for screening deep eutectic solvents (DESs) to effectively remove a diverse range of PFAS, including perfluoroalkylcarboxylic acids (PFCAs), perfluoroalkanesulfonic acids (PFSAs), and perfluoroalkyl amides (FAAs), from contaminated environments (total 13 target compounds). Utilizing energy-based screening, we identified DES candidates with high affinity for a spectrum of PFAS compounds from 1234 potential starting materials of eutectic systems. Key findings revealed the superior removal efficiency of tributylphosphineoxide/2-methylpiperazine system, exceeding 99 % for various PFAS with different carbon chain lengths in real environmental water samples. Additionally, we elucidated the molecular interactions between DESs and PFAS through ab initio molecular dynamics (AIMD) simulations, providing valuable insights into the mechanisms governing the removal process. The mechanism of extraction involves hydrogen bond network topology and structural organization, with DESs capable of extracting PFAS while maintaining a weakly aggregated state of target molecules and minimizing the impact on the intrinsic structures of DES. The proposed system forms a dynamic, complementary, and flexible non-covalent interaction network structure with PFAS. The study advances the understanding of DES as a designable, effective, and sustainable alternative to conventional solvents for PFAS remediation, offering a significant contribution to environmental chemistry and green technology.

Revisiting the Thickness of the Air-Water Interface from Two Extremes of Interface Hydrogen Bond Dynamics

The air-water interface plays a crucial role in many aspects of science because of its unique properties, such as a two-dimensional hydrogen bond (HB) network and completely different HB dynamics compared to bulk water. However, accurately determining the boundary of interfacial and bulk water, that is, the thickness of the air-water interface, still challenges experimentalists. Various simulation-based methods have been developed to estimate the thickness, converging on a range of approximately 3-10 (Å). In this study, we introduce a novel approach, grounded in density functional theory-based molecular dynamics and deep potential molecular dynamics simulations, to measure the air-water interface thickness, offering a different perspective based on prior research. To capture realistic HB dynamics in the air-water interface, two extreme scenarios of the interface HB dynamics are obtained: one underestimates the interface HB dynamics, while the other overestimates it. Surprisingly, our results suggest that the interface HB dynamics in both scenarios converges as the thickness of the air-water interface increases to 4 (Å). This convergence point, indicative of the realistic interface thickness, is also validated by our calculation of anisotropic decay of OH stretch and the free OH dynamics at the air-water interface.

Dynamic Community Detection Decouples Multiple Time Scale Behavior of Complex Chemical Systems

Although community or cluster identification is becoming a standard tool within the simulation community, traditional algorithms are challenging to adapt to time-dependent data. Here, we introduce temporal community identification using the Δ-screening algorithm, which has the flexibility to account for varying community compositions, merging and splitting behaviors within dynamically evolving chemical networks. When applied to a complex chemical system whose varying chemical environments cause multiple time scale behavior, Δ-screening is able to resolve the multiple time scales of temporal communities. This computationally efficient algorithm is easily adapted to a wide range of dynamic chemical systems; flexibility in implementation allows the user to increase or decrease the resolution of temporal features by controlling parameters associated with community composition and fluctuations therein.

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.

Ab Initio Molecular Dynamics Simulations of the Influence of Lithium Bromide on the Structure of the Aqueous Solution-Air Interface

We present the results of ab initio molecular dynamics simulations of the solution–air interface of aqueous lithium bromide (LiBr). We find that, in agreement with the experimental data and previous simulation results with empirical polarizable force field models, Br^– anions prefer to accumulate just below the first molecular water layer near the interface, whereas Li⁺ cations remain deeply buried several molecular layers from the interface, even at very high concentration. The separation of ions has a profound effect on the average orientation of water molecules in the vicinity of the interface. We also find that the hydration number of Li⁺ cations in the center of the slab N_{c,Li⁺–H2O} ≈ 4.7 ± 0.3, regardless of the salt concentration. This estimate is consistent with the recent experimental neutron scattering data, confirming that results from nonpolarizable empirical models, which consistently predict tetrahedral coordination of Li⁺ to four solvent molecules, are incorrect. Consequently, disruption of the hydrogen bond network caused by Li⁺ may be overestimated in nonpolarizable empirical models. Overall, our results suggest that empirical models, in particular nonpolarizable models, may not capture all of the properties of the solution–air interface necessary to fully understand the interfacial chemistry.

Surfactant-enhanced heterogeneity of the aqueous interface drives water extraction into organic solvents

Liquid/liquid extraction (LLE) is one of the most industrially relevant separations methods. Adsorbed surfactant is demonstrated to enhance interfacial heterogeneity and lead to water protrusions that form the basis for transport into the organic phase.

From classical to quantum and back: Hamiltonian adaptive resolution path integral, ring polymer, and centroid molecular dynamics

Path integral-based methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, including path-integral molecular dynamics, which allows for the calculation of quantum statistical properties, and ring-polymer and centroid molecular dynamics, which allow the calculation of approximate quantum dynamical properties. To this end, we derive a new integration algorithm that also makes use of multiple time-stepping. The scheme is validated via adaptive classical-path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.

Binary Solvent Organization at Silica/Liquid Interfaces: Preferential Ordering in Acetonitrile-Methanol Mixtures

Nonlinear vibrational spectroscopy experiments examined solvent organization at the silica/binary solvent interface where the binary solvent consisted of methanol and acetonitrile in varying mole fractions. Data were compared with surface vibrational spectra acquired from silica surfaces exposed to a vapor phase saturated with the same binary solvent mixtures. Changes in vibrational band intensities suggest that methanol ideally adsorbs to the silica/vapor interface but acetonitrile accumulates in excess relative to vapor-phase composition. At the silica/liquid interface, acetonitrile's signal increases until a solution phase mole fraction of ∼0.85. At higher acetonitrile concentrations, acetonitrile's signal decreases dramatically until only a weak signature persists with the neat solvent. This behavior is ascribed to dipole-paired acetonitrile forming a bilayer with the first sublayer associating with surface silanol groups and a second sublayer consisting of weakly associating, antiparallel partners. On the basis of recent simulations, we propose that the second sublayer accumulates in excess.

Thermo-molecular orientation effects in fluids of dipolar dumbbells

Plots of first-order (left) and novel second-order (right) thermomolecular orientation effects in fluids of dipolar dumbbells.

The unique molecular behavior of water at the chloroform-water interface

The molecular bonding and orientation of water at the chloroform-water interface has been examined in this study using vibrational sum-frequency spectroscopy (VSFS). The results provide a key puzzle piece towards our understanding of the systematic changes in the interfacial bonding and orientation of water that occur with variations in the polarity of the organic phase, especially when compared with previous studies of different liquid-liquid interfacial systems. In these VSFS studies the OH spectral responses of interfacial water molecules are used to characterize the interactions between water and the organic phase. The spectral analysis, aided by isotopic dilution studies, shows that the moderate polarity of the chloroform phase results in a mixed interfacial region with stronger organic-water bonding and fewer bonding interactions between adjacent water molecules than was previously found for studies of non-polar organic liquid-water interfaces. Even with the more mixed interfacial region and stronger organic-water interactions, interfacial water retains a significant amount of orientational ordering. These results are compared with recent predictions from molecular dynamics simulations about how molecules behave at the chloroform-water interface.

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.

Dynamical property of water droplets of different sizes adsorbed onto a poly(methyl methacrylate) surface

A molecular dynamics approach has been employed to study the dynamical behavior of a water droplet adsorbed on a poly(methyl methacrylate) (PMMA) surface. Several sizes of water droplets are considered in order to understand the size influence of the droplet on the dynamical properties of water molecules on the PMMA substrate. The local density profile of water molecules in the droplet upon impact with the PMMA surface is calculated, and the result shows an increase in water penetration with a decrease in the size of the droplet. By examining the velocity field, the regular motion of the water droplet is found both during the equilibrium process and after the droplet reaches the equilibrium state. The dynamical behavior of water molecule is studied by the velocity autocorrelation function (VACF) in different regions for different sizes of water droplets. The result shows that VACFs in different regions are significantly influenced for the droplet with 500 water molecules than for that with 2000 water molecules. Calculations in different regions are made for the vibrational spectrum of the oxygen atom, as well as for hydrogen bond dynamics, the lifetime, and the relaxation time of the hydrogen bond. The changes in the hydrogen bond dynamics are consistent with the change in the distribution of the hydrogen bond angle. We conclude that the dynamical properties of the water molecule are significantly affected by the region relative to the surface but only weakly influenced by the size of the droplet.

A molecular dynamics study of nanoindentation on a methyl methacrylate ultrathin film on a Au(111) substrate: interface and thickness effects

The molecular dynamics simulation model of nanoindentation is proposed in order to study the mechanical and structural deformation properties of an ultrathin MMA (methyl methacrylate) film on a Au(111) surface. First, the significant differences in the structural arrangement of MMA thin films with different thicknesses are observed. Two layers are apparent in the thinnest MMA thin film next to the Au(111) surface, while three layer structures are apparent in the thicker film. Second, this study examines the indentation tip that penetrates the MMA thin film into the Au(111) substrate in order to understand the influence of the interface on the properties and deformation behavior in both the thin film and substrate. The result shows that the indentation force is influenced both by the layer structure and by the thickness of the MMA film. The thinnest case exhibits different deformation behavior from that of the thicker cases. In addition, the deformation of MMA molecules becomes significant at the interface between the MMA film and the Au(111) surface with the increase of film thickness, and detailed deformation behavior of the Au surface for different thicknesses of MMA film is reported in this paper. Finally, both the rigid and the active models for the indentation tip are utilized in the simulation to examine the interaction differences between the tip and the film and the deformation mechanism.

Water/hydrocarbon interfaces: effect of hydrocarbon branching on single-molecule relaxation

Water/hydrocarbon interfaces are studied using molecular dynamics simulations in order to understand the effect of hydrocarbon branching on the dynamics of the system at and away from the interface. A recently proposed procedure for studying the intrinsic structure of the interface in such systems is utilized, and dynamics are probed in the usual laboratory frame as well as the intrinsic frame. The use of these two frames of reference leads to insight into the effect of capillary waves at the interface on dynamics. The systems were partitioned into zones with a width of 5 A, and a number of quantities of dynamical relevance, namely, the residence times, mean squared displacements, the velocity auto correlation functions, and orientational time correlations for molecules of both phases, were calculated in the laboratory and intrinsic frames at and away from the interface. For the aqueous phase, translational motion is found to be (a) diffusive at long times and not anomalous as in proteins or micelles, (b) faster at the interface than in the bulk, and (c) faster upon reduction of the effect of capillary waves. The rotational motion of water is (a) more anisotropic at the interface than in the bulk and (b) dependent on the orientation of the covalent O-H bond with respect to the plane of the interface. The effect of hydrocarbon branching on aqueous dynamics was found to be small, a result similar to the effect on the interfacial water structure. The hydrocarbon phase shows a larger variation for all dynamical probes, a trend consistent with their interfacial structure.

Hydrogen-Bond Structure and Dynamics at the Interface between Water and Carboxylic Acid-Functionalized Self-Assembled Monolayers

Molecular dynamics computer simulations are used to study hydrogen-bond structure and dynamics at the interface between water and carboxylic acid-functionalized self-assembled monolayers (CAFSAMs). Water-water, water-CAFSAM, and internal CAFSAM hydrogen bonds are examined. Roughly half of all adjacent carboxylic acid-terminated hydrocarbon chains are hydrogen-bonded to one another. This is consistent with experimental results reflecting two pKa values for CAFSAMs. Hydrogen-bond dynamics are expressed in terms of hydrogen-bond population autocorrelation functions and are found to be nonexponential. The water-water hydrogen-bond dynamics are slower at the interface than in the bulk, which is similar to what was found at the interface between water and weakly polar liquids such as nitrobenzene. The water-CAFSAM hydrogen bonds are found to be long-lived, on the order of tens of picoseconds. Internal CAFSAM chain-chain hydrogen bonds show almost no relaxation on the simulation time scale.

Depth Profiling of Water Molecules at the Liquid−Liquid Interface Using a Combined Surface Vibrational Spectroscopy and Molecular Dynamics Approach

The studies presented here combine experimental and computational approaches to provide new insights into how water structures and penetrates into the organic phase at two different liquid-liquid systems: the interfaces of carbon tetrachloride-water (CCl4-H2O) and 1,2-dichloroethane-water (DCE-H2O). In particular, molecular dynamics simulations are performed to generate computational spectral intensities of the CCl4-H2O and DCE-H2O interfaces that are directly comparable with experimental measurements. These simulations are then applied toward the generation of spectral profiles, responses that vary as functions of both frequency and interfacial depth. These studies emphasize the similarities and differences in the structure, orientation, and bonding of interfacial water as a function of interfacial depth for these two liquid-liquid systems and demonstrate the differing behavior of water monomers that penetrate into the organic phase.

Dynamics of Dilute Water in Carbon Tetrachloride

A dilute solution of water in a hydrophobic solvent, such as carbon tetrachloride (CCl4), presents an opportunity to study the rotational properties of water without the complicating effects of hydrogen bonds. We report here the results of theoretical, experimental, and semiempirical studies of a 0.03 mole percent solution of water in CCl4. It is shown that for this solution there are negligible water-water interactions or water-CCl4 interactions; theoretical and experimental values for proton NMR chemical shifts (deltaH) are used to confirm the minimal interactions between water and the CCl4. Calculated ab initio values and semiempirical values for oxygen-17 and deuterium quadrupole coupling constants (chi) of water/CCl4 clusters are reported. Experimental values for the 17O, 2H, and 1H NMR spin-lattice relaxation times, T1, of 0.03 mole percent water in dilute CCl4 solution at 291 K are 94+/-3 ms, 7.0+/-0.2 s, and 12.6+/-0.4 s, respectively. These T1 values for bulk water are also referenced. "Experimental" values for the quadrupole coupling constants and relaxation times are used to obtain accurate, experimental values for the rotational correlation times for two orthogonal vectors in the water molecule. The average correlation time, tauc, for the position vector of 17O (orthogonal to the plane of the molecule) in monomer water, H2(17)O, is 91 fs. The average value for the deuterium correlation time for the deuterium vector in 2H2O is 104 fs; this vector is along the OD bond. These values indicate that the motion of monomer water in CCl4 is anisotropic. At 291 K, the oxygen rotational correlation time in bulk 2H2(17)O is 2.4 ps, the deuterium rotational correlation time in the same molecule is 3.25 ps. (Ropp, J.; Lawrence, C.; Farrar, T. C.; Skinner, J. L. J. Am. Chem. Soc. 2001, 123, 8047.) These values are a factor of about 20 longer than the tauc value for dilute monomer water in CCl4.

Importance of Interfacial Adsorption in the Biphasic Hydroformylation of Higher Olefins Promoted by Cyclodextrins: A Molecular Dynamics Study at the Decene/Water Interface

We report herein a molecular dynamics study of the main species involved in the hydroformylation of higher olefins promoted by cyclodextrins in 1-decene/water biphasic systems at a temperature of 350 K. The two liquids form a well-defined sharp interface of approximately 7 A width in the absence of solute; the decene molecules are generally oriented "parallel" to the interface where they display transient contacts with water. We first focused on rhodium complexes bearing water-soluble TPPTS(3-) ligands (where TPPTS(3-) represents tris(m-sulfonatophenyl)phosphine) involved in the early steps of the reaction. The most important finding concerned the surface activity of the "active" form of the catalyst [RhH(CO)(TPPTS)(2)](6-), the [RhH(CO)(2)(TPPTS)(2)](6-) complex, and the key reaction intermediate [RhH(CO)(TPPTS)(2)(decene)](6-) (with the olefin pi-coordinated to the metal center) which are adsorbed at the water side of the interface in spite of their -6 charge. The free TPPTS(3-) ligands themselves are also surface-active, whereas the -9 charged catalyst precursor [RhH(CO)(TPPTS)(3)](9-) prefers to be solubilized in water. The role of cyclodextrins was then investigated by performing simulations on 2,6-dimethyl-beta-cyclodextrin ("CD") and its inclusion complexes with the reactant (1-decene), a reaction product (undecanal), and the corresponding key reaction intermediate [RhH(CO)(TPPTS)(2)(decene)](6-) as guests; they were all shown to be surface-active and prefer the interface over the bulk aqueous phase. These results suggest that the biphasic hydroformylation of higher olefins takes place "right" at the interface and that the CDs promote the "meeting" of the olefin and the catalyst in this peculiar region of the solution by forming inclusion complexes "preorganized" for the reaction. Our results thus point to the importance of adsorption at the liquid/liquid interface in this important phase-transfer-catalyzed reaction.

Hydrogen Bond Lifetime Dynamics at the Interface of a Surfactant Monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.