This work provides unambiguous definitions from theoretical simulations of the two interfacial regions named the BIL (binding interfacial layer) and DL (diffuse layer) at charged solid/water and air/water interfaces. The BIL and DL nomenclature follows the pioneering work of Wen et al. [Phys. Rev. Lett. 2016, 116, 016101]. Our definitions are based on the intrinsic structural properties of water only. Knowing the BIL and DL interfacial regions, one is then able to deconvolve the χ(ω) non-linear SFG (sum frequency generation) response into χ(ω) and χ(ω) contributions, thus providing a detailed molecular interpretation of these signals and of the measured total SFG. We furthermore show that the χ(ω) spectrum arises from the χ(ω) non-linear third order contribution of bulk liquid water, here calculated for several charged interfaces and shown to be universal. The χ(ω) contribution therefore has the same origin in terms of molecular normal modes at any charged interface. The molecular interpretation of χ(ω) is hence at the heart of the unambiguous molecular comprehension and interpretation of the measured total SFG signal at any charged interface.

SignificanceSilica, the most abundant mineral on Earth, is exploited in many technologies and naturally occurring geological and atmospheric processes. The physical and chemical interactions between silica and water are the fundamental driving forces for water purification systems, oil extraction, and coatings. Characterizing the silica/water interface is therefore important to improve existing technologies, in particular for silica coatings, which rely on wettability and thermal-resistant properties to remain effective. We investigated the silica/water interface using a mixture of macroscopic and microscopic techniques, including experimental and theoretical surface-specific sum frequency generation spectroscopy and contact angle measurements. Strikingly, we observed the presence of water molecules non–hydrogen bonded to the nominally hydrophilic silica surface.

We provide a detailed description of the structure of water at the interface with the air (liquid-vapor LV interface) from state-of-the-art DFT-based molecular dynamics simulations. For the first time, a two-dimensional (2D) H-bond extended network has been identified and fully characterized, demonstrating that interfacial water is organized into a 2D sheet with H-bonds oriented parallel to the instantaneous surface and following its spatial and temporal oscillations. By analyzing the nonlinear vSFG (vibrational sum frequency generation) spectrum of the LV interface in terms of layer-by-layer signal, we demonstrate that the 2D water sheet is solely responsible for the spectral signatures, hence providing the interfacial 3.5 Å thickness effectively probed in nonlinear interfacial spectroscopy. The 2D H-bond network unraveled here is the essential key to rationalize macroscopic properties of water-air interfaces, as demonstrated here for spectroscopy and the surface potential.

Water, being the universal solvent, acts as a competing agent in fundamental processes, such as folding, aggregation or biomolecular recognition. A molecular understanding of hydrophobic hydration is of central importance to understanding the subtle free energy differences, which dictate function. Ab initio and classical molecular dynamics simulations yield two distinct hydration water populations in the hydration shell of solvated tert -butanol noted as “HB-wrap” and “HB-hydration2bulk”. The experimentally observed hydration water spectrum can be dissected into two modes, centered at 164 and 195 cm –1 . By comparison to the simulations, these two bands are attributed to the “HB-wrap” and “HB-hydration2bulk” populations, respectively. We derive a quantitative correlation between the population in each of these two local water coordination motifs and the temperature dependence of the solvation entropy. The crossover from entropy to enthalpy dominated solvation at elevated temperatures, as predicted by theory and observed experimentally, can be rationalized in terms of the distinct temperature stability and thermodynamic signatures of “HB-wrap” and “HB-hydration2bulk”.

Through the prism of the rather controversial and elusive silica/water interface, ab initio DFT-based molecular dynamics simulations of the structure and non-linear SFG spectroscopy of the interface are analysed.

Hydrophobicity/hydrophilicity of aqueous interfaces at the molecular level results from a subtle balance in the water–water and water–surface interactions. This is characterized here via density functional theory–molecular dynamics (DFT-MD) coupled with vibrational sum frequency generation (SFG) and THz-IR absorption spectroscopies. We show that water at the interface with a series of weakly interacting materials is organized into a two-dimensional hydrogen-bonded network (2D-HB-network), which is also found above some macroscopically hydrophilic silica and alumina surfaces. These results are rationalized through a descriptor that measures the number of “vertical” and “horizontal” hydrogen bonds formed by interfacial water, quantifying the competition between water–surface and water–water interactions. The 2D-HB-network is directly revealed by THz-IR absorption spectroscopy, while the competition of water–water and water–surface interactions is quantified from SFG markers. The combination of SFG and THz-IR spectroscopies is thus found to be a compelling tool to characterize the finest details of molecular hydrophobicity at aqueous interfaces.

Following our recent work [Phys. Chem. Chem. Phys. 20:5190–99 (2018)] that provided the means to unambigously define and extract the three water regions at any charged interface (solid–liquid and air–liquid alike), denoted the BIL (Binding Interfacial Layer), DL (Diffuse Layer) and Bulk, and how to calculate their associated non-linear Sum Frequency Generation Spectroscopy (SFG) χ2(ω) spectroscopic contributions from Density Functional Theory (DFT)-based ab initio molecular dynamics simulations (DFT-MD/AIMD), we show here that the χDL2(ω) signal arising from the DL water region carries a wealth of essential information on the microscopic and macroscopic properties of interfaces. We show that the χDL2(ω) signal carries information on the surface potential and surface charge, the isoelectric point, EDL (Electric Double Layer) formation, and the relationship between a nominal electrolyte solution pH and surface hydroxylation state. This work is based on DFT-MD/AIMD simulations on a (0001) α–quartz–water interface and on the air–water interface, with various surface quartz hydroxylation states and various electrolyte concentrations. The conclusions drawn make use of the interplay between experiments and simulations. Most of the properties listed above can now be extracted from experimental χDL2(ω) alone with the protocols given in this work, or by making use of the interplay between experiments and simulations, as described in this work.

Four models for the calculation of the IR spectrum of gas phase molecules and clusters from molecular dynamics simulations are presented with the aim to reduce the computational cost of the usual Fourier transform (FT) of the time correlation function of the dipole moment. These models are based on the VDOS, FT of time correlation function of velocities, and atomic polar tensors (APT). The models differ from each other by the number of APTs inserted into the velocities correlation function. Excellent accuracy is achieved by the model adopting a weighted linear combination of a few selected APTs adapted for the rotation of the molecule (model D). The achieved accuracy relates to band positions, band shapes, and band intensities. Depending on the degree of actual dynamics of the molecule, rotational motion, conformational isomerization, and large amplitude motions that can be seen during the finite temperature trajectory, one could also apply one of the other models (models A, B, or C), but with caution. Model D is therefore found simple and accurate, with appealing computational cost and should be systematically applied. Its generalization to condensed phase systems should be straightforward.