Search

Dissecting the liquid water organization in contact with hydrophobic and hydrophilic surfaces is essential for understanding the chemical and physical properties of aqueous interfaces. Recently developed descriptors for microscopic hydrophobicity/hydrophilicity based either on molecular dynamics (MD) simulations or on surface-sensitive nonlinear optical techniques, such as sum frequency generation (SFG) spectroscopy, manage to capture and quantify the change in local molecular hydrophobicity at heterogeneous surfaces. However, the connections between the theoretical/structural descriptors and spectroscopic fingerprints have not been established yet. Here, we combine density functional theory-based MD simulations (DFT-MD) and both theoretical and experimental SFG spectroscopy to explore how the 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, which combines 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 pointing towards the surface 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.

We use second harmonic generation (SHG) spectroscopy, molecular dynamics simulation, and theoretical modeling to study the response of the neat liquid water-air interface to changes in the potential of an external electrode positioned near the liquid, but out of direct contact. We observe a parabolic dependence of second harmonic intensity on applied potential. Based on standard theory, we associate this dependence with the response of the diffuse layer water molecules to changes in interfacial potential profile. Taking the literature value for this response leads to the unexpected conclusion that the electric fields within the diffuse layer are opposite in sign from those originating from the electrodes. This conclusion implies that the traditional continuum-based models of interfacial screening lack the complexity necessary to properly describe the potential profile of the liquid water-vapor interface. Effects such as overscreening in the topmost interfacial layer and extended correlations in the interfacial hydrogen bonding network may play a role in governing the response of the water interface to external fields.

Glycerol is a major cryoprotective agent and is widely used to promote protein stabilization. Through a combined experimental and theoretical study, we show that global thermodynamic mixing properties of glycerol and water are dictated by local solvation motifs. We identify three water populations, i.e., bulk water, bound water H-bonded to hydrophilic groups of glycerol and wrap water hydrating hydrophobic moieties. Each population provides distinct spectroscopic fingerprints in the THz/FIR spectral range, which allow to quantify their respective abundance and their partial contributions to the mixing enthalpy. We uncover a 1:1 connection between the number of bound waters and the mixing enthalpy, as deduced from experiments as well as from simulations. The balance between local hydrophobic wrap and hydrophilic bound contributions at the molecular level dictates macroscopic thermodynamics of mixing. This offers opportunities to rationally design polyol water mixtures to optimize technological applications by tuning mixing enthalpy and entropy based on spectroscopic screening.

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.

Formation of biomolecular condensates through liquid-liquid phase separation (LLPS) has emerged as a pervasive principle in cell biology, allowing compartmentalization and spatiotemporal regulation of dynamic cellular processes. Proteins that form condensates under physiological conditions often contain intrinsically disordered regions with low-complexity domains. Among them, the RNA-binding proteins FUS and TDP-43 have been a focus of intense investigation because aberrant condensation and aggregation of these proteins is linked to neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal dementia. LLPS occurs when protein-rich condensates form surrounded by a dilute aqueous solution. LLPS is per se entropically unfavorable. Energetically favorable multivalent protein-protein interactions are one important aspect to offset entropic costs. Another proposed aspect is the release of entropically unfavorable preordered hydration water into the bulk. We used attenuated total reflection spectroscopy in the terahertz frequency range to characterize the changes in the hydrogen bonding network accompanying the FUS enrichment in liquid-liquid phase-separated droplets to provide experimental evidence for the key role of the solvent as a thermodynamic driving force. The FUS concentration inside LLPS droplets was determined to be increased to 2.0 mM independent of the initial protein concentration (5 or 10 μM solutions) by fluorescence measurements. With terahertz spectroscopy, we revealed a dewetting of hydrophobic side chains in phase-separated FUS. Thus, the release of entropically unfavorable water populations into the bulk goes hand in hand with enthalpically favorable protein-protein interaction. Both changes are energetically favorable, and our study shows that both contribute to the thermodynamic driving force in phase separation.

Electrostatic interactions are central to the structure and function of nucleic acids, including their folding, condensation, and interaction with proteins and other charged molecules. These interactions are profoundly affected by ions surrounding nucleic acids, the constituents of the so-called ion atmosphere. Here, we report precise Fourier Transform-Terahertz/Far-Infrared (FT-THz/FIR) measurements in the frequency range 30–500 cm(−1) for a 24-bp DNA solvated in a series of alkali halide (NaCl, NaF, KCl, CsCl, and CsF) electrolyte solutions which are sensitive to changes in the ion atmosphere. Cation excess in the ion atmosphere is detected experimentally by observation of cation modes of Na(+), K(+), and Cs(+) in the frequency range between 70–90 cm(−1). Based on MD simulations, we propose that the magnitude of cation excess (which is salt specific) depends on the ability of the electrolyte to perturb the water network at the DNA interface: In the NaF atmosphere, the ions reduce the strength of interactions between water and the DNA more than in case of a NaCl electrolyte. Here, we explicitly take into account the solvent contribution to the chemical potential in the ion atmosphere: A decrease in the number of bound water molecules in the hydration layer of DNA is correlated with enhanced density fluctuations, which decrease the free energy cost of ion-hydration, thus promoting further ion accumulation within the DNA atmosphere. We propose that taking into account the local solvation is crucial for understanding the ion atmosphere.

International audience; We show that phase-sensitive vibrational sum frequency generation (SFG) spectra of solid/water and air/water interfaces, neutral and charged, can be successfully predicted using classical molecular dynamics (CMD) simulations in combination with simple nonpolarizable force fields (FFs). This can be achieved when employing velocity–velocity autocorrelation functions weighted by parameterized Raman and atomic polar tensors for the computation of the SFG. This procedure avoids computing polarizability tensors and dipole moments using either costly ab initio molecular dynamics (AIMD) simulations or CMD simulations with more complex and computationally demanding FFs. Such a methodology paves the way to a broad usage and computationally low-cost theoretical SFG spectroscopy, as even flexible nonpolarizable water models and common FFs for inorganic surfaces can provide good predictions of the SFG spectra, in rather good qualitative agreement with AIMD and/or experiments. The strongly reduced computational cost in our approach opens the possibility to study larger systems for long periods of time, for example, allowing a detailed characterization of the electric double-layer formation at interfaces with “environmentally relevant” ionic concentrations (mM), extracting fingerprints by theoretical CMD–SFG spectroscopy.

International audience; 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”.

We use second harmonic generation (SHG) spectroscopy, molecular dynamics simulation, and theoretical modeling to study the response of the neat liquid water-air interface to changes in the potential of an external electrode positioned above the liquid out of contact. We observe a parabolic dependence of second harmonic intensity on applied potential. This dependence is reminiscent of bulk-phase electric field induced second harmonic (EFISH) but more complicated because it combines the second-order response of the topmost water layer and the potential dependent response of the interfacial electrical double-layer. Based on the literature values for these contributions, we derive a physical interpretation of our measurements that reveals new insight into the response of the neat water interface to external electric fields. Specifically, we find that the net dipolar orientation of water molecules within the double-layer is primarily responsive to the internal fields generated by the excess surface concentrations of OH- and H3O+ that arise to screen the external potential. Notably, this interpretation implies that the orientational response of water dipoles at the interface can actually oppose the direction of the external field, a subtle effect that is not captured by traditional models.

Hydration free energies are dictated by a subtle balance of hydrophobic and hydrophilic interactions. which is crucial for many biological processes and technological applications, such as protein folding and molecular recognition. Whereas so far the overall entropy and enthalpy are experimentally determined based on equilibrium measurements using a calorimeter, we present here a pure spectroscopic access to these important observables, which give direct access to the underlying molecular mechanism that determines these driving forces. Using THz calorimetry the contributions due to cavity formation and hydrophilic interactions can be traced back to changes in the intermolecular hydrogen bond stretching region around 150-200 cm−1 and spectroscopic changes due to strong solute-water interactions in the frequency range of the librational modes, i.e. between 540 and 600 cm−1. Thus, we are able to link the thermodynamic model of the Lum-Chandler-Weeks theory, which was a pure ”Gedankenexperiment”, directly to experimental observables. We show that alcohol hydration can be described by a sum of a free energy cost of forming and wrapping a cavity around the solute (which is entropic for small alcohols) and an enthalpic gain due to the hydrogen bonds formed between the alcohol OH group and bound water molecules around it. In the future, our approach will allow to quantify entropic cost and enthalpic gain not only in equilibrium but also in non-equilibrium processes.

Hydration free energies are dictated by a subtle balance of hydrophobic and hydrophilic interactions. We present here a spectroscopic approach, which gives direct access to the two main contributions: Using THz-spectroscopy to probe the frequency range of the intermolecular stretch (150-200 cm

Significance For centuries the double layer at the solid/electrolyte interface has been a central concept in electrochemistry. Today, it is still crucial for virtually all renewable energy storage and conversion technologies. Here, the double-layer formation is probed by THz spectroscopy with ultrabright synchrotron light as a source. Our results capture the molecular details of double-layer formation at positively/negatively charged Au electrodes for an NaCl electrolyte. We reveal a contrasting response applying positive versus negative bias, which is dictated by the interfacial water network and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations. While Na+ is directly attracted toward the negatively charged electrode, stripping of the Cl− hydration shell is observed only at larger potential values., The double layer at the solid/electrolyte interface is a key concept in electrochemistry. Here, we present an experimental study combined with simulations, which provides a molecular picture of the double-layer formation under applied voltage. By THz spectroscopy we are able to follow the stripping away of the cation/anion hydration shells for an NaCl electrolyte at the Au surface when decreasing/increasing the bias potential. While Na+ is attracted toward the electrode at the smallest applied negative potentials, stripping of the Cl− hydration shell is observed only at higher potential values. These phenomena are directly measured by THz spectroscopy with ultrabright synchrotron light as a source and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations.

The electrical double-layer plays a key role in important interfacial electrochemical processes from catalysis to energy stor-age and corrosion. Therefore, understanding its structure is crucial for the progress of sustainable technologies. We extract new physico-chemical information on the capacitance and structure of the electrical double-layer of platinum and gold nanoparticles at the molecular level, employing single nanoparticle electrochemistry. We reveal that the charge storage ability of the solid/liquid interface is larger by one order-of-magnitude than predicted by the traditional mean-field models of the double-layer such as the Gouy-Chapman-Stern-model. Performing Molecular Dynamics simulations, we investigate the possible relationship between the measured high capacitance and adsorption strength of the water adlayer formed at the metal surface. These insights may launch the active tuning of solid-solvent and solvent-solvent interactions as innovative design strategy to transform energy technologies towards superior performance and sustainability.

The double layer at the solid/electrolyte interface is a key concept in electrochemistry. Here, we present an experimental study combined with simulations, which provides a molecular picture of the double-layer formation in operando processes. By THz spectroscopy we are able to follow the stripping off of the cation/anion hydration shells for a NaCl electrolyte at the Au surface when decreasing/increasing the bias potential. While Na+ is attracted toward the electrode already at the smallest applied negative potentials, stripping-off of the Cl- hydration shell is observed only at higher potential values. These phenomena are directly measured by in operando THz spectroscopy with ultra-bright synchrotron light as a source and rationalized by accompanying molecular-dynamics simulations and electronic-structure calculations.

Significance Using a supramolecular assembly that catalyzes substrates in water, we show that the nanocage encapsulates a water cluster that is structurally and dynamically distinct from any known phase of water. It plays an important role in the driving force for guest encapsulation: The release of the highly unusual encapsulated water droplet creates a strong thermodynamic drive for the high-affinity binding of guests in aqueous solution for the supramolecular construct., Water under nanoconfinement at ambient conditions has exhibited low-dimensional ice formation and liquid–solid phase transitions, but with structural and dynamical signatures that map onto known regions of water’s phase diagram. Using terahertz (THz) absorption spectroscopy and ab initio molecular dynamics, we have investigated the ambient water confined in a supramolecular tetrahedral assembly, and determined that a dynamically distinct network of 9 ± 1 water molecules is present within the nanocavity of the host. The low-frequency absorption spectrum and theoretical analysis of the water in the Ga4L612− host demonstrate that the structure and dynamics of the encapsulated droplet is distinct from any known phase of water. A further inference is that the release of the highly unusual encapsulated water droplet creates a strong thermodynamic driver for the high-affinity binding of guests in aqueous solution for the Ga4L612− supramolecular construct.

The structure and ultrafast dynamics of the electric double layer (EDL) are central to chemical reactivity and physical properties at solid/aqueous interfaces. While the Gouy-Chapman-Stern model is widely used to describe EDLs, it is solely based on the macroscopic electrostatic attraction of electrolytes for the charged surfaces. Structure and dynamics in the Stern layer are, however, more complex because of competing effects due to the localized surface charge distribution, surface-solvent-ion correlations, and the interfacial hydrogen bonding environment. Here, we report combined time-resolved vibrational sum frequency generation (TR-vSFG) spectroscopy with ab initio DFT-based molecular dynamics simulations (AIMD/DFT-MD) to get direct access to the molecular-level understanding of how ions change the structure and dynamics of the EDL. We show that innersphere adsorbed ions tune the hydrophobicity of the silica-aqueous interface by shifting the structural makeup in the Stern layer from dominant water-surface interactions to water-water interactions. This drives an initially inhomogeneous interfacial water coordination landscape observed at the neat interface toward a homogeneous, highly interconnected in-plane 2D hydrogen bonding (2D-HB) network at the ionic interface, reminiscent of the canonical, hydrophobic air-water interface. This ion-induced transformation results in a characteristic decrease of the vibrational lifetime (T1) of excited interfacial O-H stretching modes from T1 ∼ 600 fs to T1 ∼ 250 fs. Hence, we propose that the T1 determined by TR-vSFG in combination with DFT-MD simulations can be widely used for a quantitative spectroscopic probe of the ion kosmotropic/chaotropic effect at aqueous interfaces as well as of the ion-induced surface hydrophobicity.

International audience; 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 χ(2)(ω) non-linear SFG (sum frequency generation) response into χ(2)BIL(ω) and χ(2)DL(ω) contributions, thus providing a detailed molecular interpretation of these signals and of the measured total SFG. We furthermore show that the χ(2)DL(ω) spectrum arises from the χ(3)(ω) non-linear third order contribution of bulk liquid water, here calculated for several charged interfaces and shown to be universal. The χ(2)DL(ω) contribution therefore has the same origin in terms of molecular normal modes at any charged interface. The molecular interpretation of χ(2)BIL(ω) is hence at the heart of the unambiguous molecular comprehension and interpretation of the measured total SFG signal at any charged interface.

Metal/water interfaces catalyze a large variety of chemical reactions, which often involve small hydrophobic molecules. In the present theoretical study we show that hydrophobic hydration at the Au(100)/water interface actively contributes to the reaction free energy by up to several hundreds of meV. This occurs either in adsorption/desorption reaction steps, where the vertical distance from the surface changes in going from reactants to products, or in addition and elimination reaction steps, where two small reactants merge into a larger product and viceversa. We find that size and position effects cannot be captured by treating them as independent variables. Instead, their simultaneous evaluation allows to map the important contributions, and we provide examples of their combinations for which interfacial reactions can be either favoured or disfavoured. By taking a N2 and a CO2 reduction pathway as test cases, we show that explicitly considering hydrophobic effects is important for the selectivity and rate of these relevant interfacial processes.

DFT-based molecular dynamics simulations of the electrified air-liquid water interface are presented, where a homogeneous field is applied parallel to the surface plane. We unveil the field intensity for the onset of proton transfer and molecular dissociation; the protonic current/proton conductivity is measured as a function of the field intensity/voltage. The air-water interface is shown to exhibit a proton conductivity twice the one in the liquid water for field intensities below 0.40 V angstrom(-1). We show that this difference arises from the very specific organization of water in the binding interfacial layer (BIL, i.e. the air-water interface region) into a 2D-HBond-network that is maintained and enforced at the electrified interface. Beyond fields of 0.40 V angstrom(-1), water in the BIL and in the bulk liquid are aligned in the same way by the rather intense fields, hence leading to the same proton conductivity in both BIL and bulk water.

Significance The optimization of “green” electrochemical processes is one of the most important challenges in the transition toward renewable energy technologies. In many of these processes, including, e.g., C O 2 and N 2 reduction, small hydrophobic molecules are formed and react at the interface, and their hydration free energy modulates the associated thermodynamics. Here, we use molecular dynamics simulations to elucidate the mechanisms and energetics of hydrophobic hydration at an electrified gold/water interface. We propose an adaptation of the Lum–Chandler–Weeks theory that maps the changes in hydration free energies at the interface as a function of solute size and applied potential.

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. Following our recent work [Phys. Chem. Chem. Phys., 2018, 20, 5190-5199], we show that once the interfacial water is decomposed into BIL (Binding Interfacial Layer) and DL (Diffuse Layer) interfacial regions, the SFG signals can be deconvolved and unambiguously interpreted, and a global microscopic understanding on silica/water interfaces can be obtained. By comparing crystalline quartz/water and amorphous (fused) silica/water interfaces, the dependence of interfacial structural and spectroscopic properties on the degree of surface crystallinity is established, while by adding KCl electrolytes at the quartz/water interface, the chaotropic effect of ions on the interfacial molecular arrangement is unveiled. The evolution of structure and SFG spectra of silica/water interfaces with respect to increasing surface deprotonation, i.e., with respect to pH conditions, is also evaluated. Spectroscopic BIL-SFG markers that experimentally allow one detect the water order/disorder in the BIL as a function of surface hydroxylation and ion concentration are revealed, while the pH-induced modulations in the experimentally recorded SFG spectra are rationalized in terms of changes in both BIL and DL SFG signatures.

Within the general context of the electrochemical oxygen evolution reaction of the water oxidation/electrolysis, we focus on one essential aspect of electrochemical interfaces, i.e., the comprehension of the interaction and organisation of liquid water at the (semiconductor) (110)-Co

Significance Silica, 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., Interfaces between water and silicates are ubiquitous and relevant for, among others, geochemistry, atmospheric chemistry, and chromatography. The molecular-level details of water organization at silica surfaces are important for a fundamental understanding of this interface. While silica is hydrophilic, weakly hydrogen-bonded OH groups have been identified at the surface of silica, characterized by a high O-H stretch vibrational frequency. Here, through a combination of experimental and theoretical surface-selective vibrational spectroscopy, we demonstrate that these OH groups originate from very weakly hydrogen-bonded water molecules at the nominally hydrophilic silica interface. The properties of these OH groups are very similar to those typically observed at hydrophobic surfaces. Molecular dynamics simulations illustrate that these weakly hydrogen-bonded water OH groups are pointing with their hydrogen atom toward local hydrophobic sites consisting of oxygen bridges of the silica. An increased density of these molecular hydrophobic sites, evident from an increase in weakly hydrogen-bonded water OH groups, correlates with an increased macroscopic contact angle.

International audience; Within the general context of the electrochemical oxygen evolution reaction of the water oxidation/electrolysis, we focus on one essential aspect of electrochemical interfaces, i.e., the comprehension of the interaction and organisation of liquid water at the (semiconductor) (110)-Co3O4 surface using density functional theory-molecular dynamics simulations. A detailed characterization of the chemical and physical properties of the aqueous interface is provided in terms of structure, dynamics, electric field, work function, and spectroscopy, as a preliminary step into the modelling of the (110)-Co3O4 aqueous surface in more relevant electrochemical conditions. The water at the aqueous B-termination is, in particular, shown more dynamical than that at the A-termination and more “undisciplined”: the water is indeed mostly an HB-acceptor with the solid, with an orientation of their dipole moments found opposite the field generated by the negative surface charge. At both aqueous interfaces, the work function is twice lower than that at the bare (non-hydroxylated) surfaces. The SFG (Sum Frequency Generation) spectroscopy is shown dominated by the water in the diffuse layer, while the SFG signal from the binding interfacial layer reflects the single orientation of water at the aqueous A-termination and the two orientations of water at the aqueous B-termination.

International audience; Graph theory algorithms have been proposed in order to identify, follow in time, and statistically analyze the changes in conformations that occur along molecular dynamics (MD) simulations. The atomistic granularity level of the MD simulations is maintained within the graph theoric algorithms proposed here, isomorphism is a key component together with keeping the chemical nature of the atoms. Isomorphism is used to recognize conformations and construct the graphs of transitions, and the reduction in complexity of the isomorphism has been achieved by the introduction of "orbits" and "reference snapshots." The proposed algorithms are applied to MD trajectories of gas phase molecules and clusters as well as condensed matter. The changes in conformations followed over time are hydrogen bond(s), proton transfer(s), coordination number(s), covalent bond(s), multiple fragmentation(s), and H-bonded membered rings. The algorithms provide an automatic analysis of multiple trajectories in parallel, and can be applied to ab initio and classical MD trajectories alike, and to more coarse grain representations.

International audience; An extensive structural characterization of the 2D-HB-Network formed by water molecules at the air-water interface has been carried out by means of DFT-based and classical molecular dynamics simulations. SPC/E and SPC/Fw force fields commonly used for modelling liquid water are shown to correctly reproduce the 2D-HB-Network structure obtained from the reference ab-initio simulation. For both DFT-MD and classical MD representations, identical results have been obtained by increasing the size and time-scale of the simulations, starting from a simulation box of 256 water molecules in the liquid phase, simulated for 20 ps up to 25 ns. One pivotal result is that ∼90% of the water molecules in the interfacial layer are connected by a collective and extended net of HBs oriented parallel to the surface, and this 2D-HB-Network is built upon adjacent water rings mostly formed by 4, 5 or 6 water molecules. This finding in particular rationalizes previous non linear SFG spectroscopic results suggesting that liquid water at the interface with hydrophobic media could arrange in rings.

Following our recent work [Phys. Chem. Chem. Phys. 20:5190&ndash, 99 (2018)] that provided the means to unambigously define and extract the three water regions at any charged interface (solid&ndash, liquid and air&ndash, 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) &chi, 2(&omega, ) spectroscopic contributions from Density Functional Theory (DFT)-based ab initio molecular dynamics simulations (DFT-MD/AIMD), we show here that the &chi, DL2(&omega, ) 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 &chi, ) 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) &alpha, &ndash, quartz&ndash, water interface and on the air&ndash, 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 &chi, ) alone with the protocols given in this work, or by making use of the interplay between experiments and simulations, as described in this work.

International audience; Following our previous work where the existence of a special 2-Dimensional H-Bond (2D-HB)-Network was revealed at the air-water interface [S. Pezzotti et al., J. Phys. Chem. Lett. 8, 3133 (2017)], we provide here a full structural and dynamical characterization of this specific arrangement by means of both Density Functional Theory based and Force Field based molecular dynamics simulations. We show in particular that water at the interface with air reconstructs to maximize H-Bonds formed between interfacial molecules, which leads to the formation of an extended and non-interrupted 2-Dimensional H-Bond structure involving on average ∼90% of water molecules at the interface. We also show that the existence of such an extended structure, composed of H-Bonds all oriented parallel to the surface, constrains the reorientional dynamics of water that is hence slower at the interface than in the bulk. The structure and dynamics of the 2D-HB-Network provide new elements to possibly rationalize several specific properties of the air-water interface, such as water surface tension, anisotropic reorientation of interfacial water under an external field, and proton hopping.

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 χ

International audience; 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.