Deconvolution of BIL-SFG and DL-SFG spectroscopic signals reveals order/disorder of water at the elusive aqueous silica interface

Phase of the second-order susceptibility in vibrational sum frequency generation spectroscopy: Origins, utility, and measurement techniques

Vibrational sum frequency generation can provide valuable structural information at surfaces and buried interfaces. Relating the measured spectra to the complex-valued second-order susceptibility χ(2) is at the heart of the technique and a requisite step in nearly all subsequent analyses. The magnitude and phase of χ(2) as a function of frequency reveal important information about molecules and materials in regions where centrosymmetry is broken. In this tutorial-style perspective, the origins of the χ(2) phase are first described, followed by the utility of phase determination. Finally, some practical methods of phase extraction are discussed.

Deprotonation of formic, acetic acids and bicarbonate ion in slit silica nanopores at infinite dilution and in the presence of electrolytes

Dielectric effects and the coupled electrostatics between the nanoconfined and the internal/external aqueous media contribute to the observed deviations of chemistry within the nanoconfined environment when compared with unconfined systems. A systematic understanding has remained elusive, especially with respect to background salt concentration and boundary condition effects like the nanopore surface chemistry and the reference state used to calculate free energies. We utilize molecular dynamics simulations along with thermodynamic integration to determine the free energy difference associated with acid-base chemistry in 2 nm and 4 nm slit pores open to a bulk-like reservoir. pKa increases are predicted when confining acetic acid, formic acid, and bicarbonate in the slits at infinite dilution conditions. We find that confinement weakens the acids, and the modulation of outer pore surface dipole magnitudes can tune the pKa shift values, suggesting that purely "intrinsic" electrostatic effect on confinement may not exist. At sufficiently high salt concentrations, the dielectric/electrostatic effects on pKa values diminish due to charge screening effects. These discoveries enable future modifications of nanopore chemistries to achieve desirable properties for industrial applications.

Lead (II) ions enable the ion-specific effects of monovalent anions on the molecular structure and interactions at silica/aqueous interfaces

HYPOTHESIS

The adsorption of heavy metal ions such as Pb(II) onto negatively charged minerals such as silica is expected to alter the structure and the interactions at the silica/aqueous interfaces. Besides the solution pH, the inner-sphere sorption of Pb(II) is expected to regulate the surface charge/potential, hypothesized to control the actions of monovalent anions in the aqueous environment. These complex pictures can be probed directly using surface-sensitive sum-frequency generation (SFG) spectroscopy.

EXPERIMENTS

The pH-dependent water structure within the double layer at silica/aqueous interfaces under the influence of different ions was examined using SFG. The recorded SFG spectra were deconvoluted into the Stern layer (SL) and diffuse layer (DL) using the maximum entropy method in conjunction with the electrical double-layer theory.

FINDINGS

Standalone monovalent sodium salts do not exhibit ion-specific effects on the silica/aqueous interfaces. However, the mixture of Pb(II) species and each of these salts display profound ion-specific effects on the structure of silica/aqueous interfaces, indicating the role of Pb(II) as an enabler of the ion-specificity of the investigated monovalent anions. The interesting effect arises from a complex interplay between the physical processes (i.e., electrostatic interactions, screening effects, etc.) and chemical processes such as the hydrolysis of Pb(II) ions, ion complexation, protonation and deprotonation of the surface silanol group.

Variable-Angle Surface Spectroscopy Reveals the Water Structure in the Stern Layer at Charged Aqueous Interfaces

At charged aqueous interfaces, the second-order nonlinear optical response originates from water molecules within the diffuse part of the electrical double layer, which are ordered by the surface field and from water that additionally experiences chemical and physical interactions with the surface in the Stern layer. These two environments can either reinforce or diminish the overall signal and can be disentangled by varying the coherence length of their interaction with external laser fields. Here, we demonstrate a method in which the angle of incidence is varied to afford a significant change in the coherence length. When this technique was applied to the silica-water interface, it was observed that water molecules in the Stern and diffuse layers direct their hydrogen atoms toward the mineral surface at a low ionic strength and neutral pH. A decrease in the signal with increasing ionic strength is attributed to hydrated cation adsorption that competes with free water for deprotonated silanol sites.

Water Film Structure and Wettability of Different Quartz Surfaces: Hydrogen Bonding Across Various Cutting Planes

Quartz is ubiquitous in subsurface formations. The crystal faces have different atomic arrangements. Knowledge of the molecular structures on the surface of quartz is key in many processes. Molecular dynamics simulations are conducted to investigate the atomic arrangement effect on the water film structure, ion adsorption, and wettability at three different α-quartz surfaces. The interfacial structures depend on the quartz surface. Intrasurface hydrogen bonding between surface silanols differs in the α-quartz surface. At the (0001) surface, the OH density is 9.58 nm-2, consisting of Q2 units with two hydroxyl groups per silicone atom. At the (101̅0)-β surface, the OH density is 7.54 nm-2, consisting of Q3 units with one hydroxyl group per silicone atom; there is significant intrasurface hydrogen bonding. At the (101̅0)-α surface, the OH density is 7.54 nm-2, consisting of Q2 units; however, there is little intrasurface hydrogen bonding. The intrasurface hydrogen bonding results in the exposure of hydrogen-bond acceptors to the aqueous phase, causing water molecules to have an H-up (hydrogen toward surface) orientation. This orientation can be found at the (0001) and (101̅0)-β surfaces; it is related to the degree of ordering at the surface. The ordering at the (0001) and (101̅0)-β surfaces is higher than that at the (101̅0)-α surface. In aqueous systems with ions, cation adsorption is the most dominant at the (0001) surface due to the largest surface density of the intrasurface hydrogen bonding, providing interaction sites for cations to be adsorbed. We observe a pronounced decrease in water film thickness from the ions at the (0001) surface only, likely due to significant cation adsorption. In this work, we demonstrate that the hydrogen-bond network, which varies from the plane cut, affects the water film structure and ion adsorption. The contact is nearly zero despite the changes in the film thickness and molecular structure at the temperature of 318 K.

Microscopic properties of forces from ice solidification interface acting on silica surfaces based on molecular dynamics simulations

The origin of the forces acting on a silica surface from an ice solidification interface was investigated to understand the solidification phenomenon and its impact on nanometer-scale structures using molecular dynamics simulations. The microscopic forces were determined by appropriately averaging the forces acting on the silica wall from the water molecules in time and space; the time evolutions of these microscopic forces during the solidification processes were investigated for three types of silica surfaces. The results indicate that the microscopic forces fluctuate more after the solidification interface makes contact with the wall surface. To visualize the changes in the microscopic forces and hydrogen bonds due to solidification, their differences compared to the liquid state were calculated. When the solidification interface is near the wall, the changes in these microscopic forces and hydrogen bonds due to solidification are correlated. This tendency is more significant for an amorphous wall and a wall with a structure than for a crystalline wall. The changes in the microscopic force depend on the water molecules that behave as acceptors when forming the hydrogen bonds with the wall and on the configuration of the silanol groups on the silica surfaces.

Evolution of the electrical double layer with electrolyte concentration probed by second harmonic scattering

Investigating the electrical double layer (EDL) structure has been a long-standing challenge and has seen the emergence of several sophisticated techniques able to probe selectively the few molecular layers of a solid/water interface. While a qualitative estimation of the thickness of the EDL can be obtained using simple theoretical models, following experimentally its evolution is not straightforward and can be even more complicated in nano- or microscale systems, particularly when changing the ionic concentration by several orders of magnitude. Here, we bring insight into the structure of the EDL of SiO2 nanoparticle suspensions and its evolution with increasing ionic concentration using angle-resolved second harmonic scattering (AR-SHS). Below millimolar salt concentrations, we can successively characterize inner-sphere adsorption, diffuse layer formation, and outer-sphere adsorption. Moreover, we show for the first time that, by appropriately selecting the nanoparticle size, it is possible to retrieve information also in the millimolar range. There, we observe a decrease in the magnitude of the surface potential corresponding to a compression in the EDL thickness, which agrees with the results of several other electroanalytical and optical techniques. Molecular dynamics simulations suggest that the EDL compression mainly results from the diffuse layer compression rather than outer-sphere ions (Stern plane) moving closer to the surface.

Effects of Charged Surfactants on Interfacial Water Structure and Macroscopic Properties of the Air-Water Interface

Surfactants are used to control the macroscopic properties of the air-water interface. However, the link between the surfactant molecular structure and the macroscopic properties remains unclear. Using sum-frequency generation spectroscopy and molecular dynamics simulations, two ionic surfactants (dodecyl trimethylammonium bromide, DTAB, and sodium dodecyl sulphate, SDS) with the same carbon chain lengths and charge magnitude (but different signs) of head groups interact and reorient interfacial water molecules differently. DTAB forms a thicker but sparser interfacial layer than SDS. It is due to the deep penetration into the adsorption zone of Br- counterions compared to smaller Na+ ones, and also due to the flip-flop orientation of water molecules. SDS alters two distinctive interfacial water layers into a layer where H+ points to the air, forming strong hydrogen bonding with the sulphate headgroup. In contrast, only weaker dipole-dipole interactions with the DTAB headgroup are formed as they reorient water molecules with H+ point down to the aqueous phase. Hence, with more molecules adsorbed at the interface, SDS builds up a higher interfacial pressure than DTAB, producing lower surface tension and higher foam stability at a similar bulk concentration. Our findings offer improved knowledge for understanding various processes in the industry and nature.

Probing interfacial water structure induced by charge reversal and hydrophobicity of silica surface in the presence of divalent heavy metal ions using sum frequency generation spectroscopy

HYPOTHESIS

Adsorption of divalent heavy metal ions (DHMIs) at the mineral-water interfaces changes interfacial chemical species and charges, interfacial water structure, Stern (SL), and diffuse (DL) layers. These molecular changes can be detected by probing changing orientation and hydrogen-bond network of interfacial water molecules in response to changing local charges and hydrophobicity.

EXPERIMENTS

Sum-frequency generation (SFG) spectroscopy was used to probe changes in vibrational resonances of interfacial OH vs. DHMI concentration and pH. SFG spectra were deconvoluted using the measured surface potential and maximum entropy method in conjunction with the electrical double-layer theory for the SL and DL structures and correlated by hydrophobicity.

FINDINGS

Three surface charge reversals (CRs) were detected at low (CR1), medium (CR2), and high (CR3) pHs. Unlike CR1, SFG signals were minimized at CR2 and CR3 for DHMIs-silica systems highlighting considerable alterations in the structure of interfacial waters due to the inner-sphere sorption of metal hydroxo complexes. SFG results showed "hydrophobic-like" stretching modes at > 3600 cm^-1 for Pb-, Cu-, and Zn-treated silica. However, contact angle measurements revealed the hydrophobization of silica only in the presence of Pb(II), as confirmed by an in-depth SFG analysis of the hydrogen-bond network of the interfacial water molecules in the SL.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Elucidation of the pH-Dependent Electric Double Layer Structure at the Silica/Water Interface Using Heterodyne-Detected Vibrational Sum Frequency Generation Spectroscopy

The silica/water interface is one of the most abundant charged interfaces in natural environments, and the elucidation of the water structure at the silica/water interface is essential. In the present study, we measured the interface-selective vibrational (χ(2)) spectra in the OH stretch region of the silica/water interface in a wide pH range of pH 2.0-12.0 while changing the salt concentration by heterodyne-detected vibrational sum-frequency generation spectroscopy. With the help of singular value decomposition analysis, it is shown that the imaginary part of the χ(2) (Imχ(2)) spectra can be decomposed into the spectra of the diffuse Gouy-Chapman layer (DL) and the compact Stern layer (SL), which enables us to quantitatively analyze the spectra of DL and SL separately. The salt-concentration dependence of the DL spectra at different pH values is analyzed using the modified Gouy-Chapman theory, and the pH-dependent surface charge density and the pKa value (4.8 ± 0.2) of the silica/water interface are evaluated. Furthermore, it is found that the pH-dependent change of the SL spectra is quantitatively explained by three spectral components that represent the three characteristic water species appearing in different pH regions in SL. The quantitative understanding obtained from the analysis of each spectral component in the Imχ(2) spectra provides a clear molecular-level picture of the electric double layer at the silica/water interface.

On the Trail of Molecular Hydrophilicity and Hydrophobicity at Aqueous Interfaces

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

Nature of Cations Critically Affects Water at the Negatively Charged Silica Interface

Understanding the collective behavior of ions at charged surfaces is of paramount importance for geological and electrochemical processes. Ions screen the surface charge, and interfacial fields break the centro-symmetry near the surface, which can be probed using second-order nonlinear spectroscopies. The effect of electrolyte concentration on the nonlinear optical response has been semi-quantitatively explained by mean-field models based on the Poisson–Boltzmann equation. Yet, to explain previously reported ion-specific effects on the spectroscopic response, drastic ion-specific changes in the interfacial properties, including surface acidities and dielectric permittivities, or strong ion adsorption/desorption had to be invoked. Here, we use sum-frequency generation (SFG) spectroscopy to probe the symmetry-breaking of water molecules at a charged silica surface in contact with alkaline metal chloride solutions (LiCl, NaCl, KCl, and CsCl) at various concentrations. We find that the water response varies with the cation: the SFG response is markedly enhanced for LiCl compared to CsCl. We show that within mean-field models, neither specific ion–surface interactions nor a reduced dielectric constant of water near the interface can account for the variation of spectral intensities with cation nature. Molecular dynamics simulations confirm that the decay of the electrochemical potential only weakly depends on the salt type. Instead, the effect of different salts on the optical response is indirect, through the reorganization of the interfacial water: the salt-type-dependent alignment of water directly at the interface can explain the observations.

Water Structure in the Electrical Double Layer and the Contributions to the Total Interfacial Potential at Different Surface Charge Densities

The electric double layer governs the processes of all charged surfaces in aqueous solutions; however, elucidating the structure of the water molecules is challenging for even the most advanced spectroscopic techniques. Here, we present the individual Stern layer and diffuse layer OH stretching spectra at the silica/water interface in the presence of NaCl over a wide pH range using a combination of vibrational sum frequency generation spectroscopy, heterodyned second harmonic generation, and streaming potential measurements. We find that the Stern layer water molecules and diffuse layer water molecules respond differently to pH changes: unlike the diffuse layer, whose water molecules remain net-oriented in one direction, water molecules in the Stern layer flip their net orientation as the solution pH is reduced from basic to acidic. We obtain an experimental estimate of the non-Gouy-Chapman (Stern) potential contribution to the total potential drop across the insulator/electrolyte interface and discuss it in the context of dipolar, quadrupolar, and higher order potential contributions that vary with the observed changes in the net orientation of water in the Stern layer. Our findings show that a purely Gouy-Chapman (Stern) view is insufficient to accurately describe the electrical double layer of aqueous interfaces.

First hyperpolarizability of water at the air-vapor interface: a QM/MM study questions standard experimental approximations

Surface Second-Harmonic Generation (S-SHG) experiments provide a unique approach to probe interfaces. One important issue for S-SHG is how to interpret the S-SHG intensities at the molecular level. Established frameworks commonly assume that each molecule emits light according to an average molecular hyperpolarizability tensor β(-2ω,ω,ω). However, for water molecules, this first hyperpolarizability is known to be extremely sensitive to their environment. We have investigated the molecular first hyperpolarizability of water molecules within the liquid-vapor interface, using a quantum description with explicit, inhomogeneous electrostatic embedding. The resulting average molecular first hyperpolarizability tensor depends on the distance relative to the interface, and it practically respects the Kleinman symmetry everywhere in the liquid. Within this numerical approach, based on the dipolar approximation, the water layer contributing to the Surface Second Harmonic Generation (S-SHG) intensity is less than a nanometer. The results reported here question standard interpretations based on a single, averaged hyperpolarizability for all molecules at the interface. Not only the molecular first hyperpolarizability tensor significantly depends on the distance relative to the interface, but it is also correlated to the molecular orientation. Such hyperpolarizability fluctuations may impact the S-SHG intensity emitted by an aqueous interface.

Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go

We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.

Gate Alignment of Liquid Water Molecules in Electric Double Layer

The behavior of liquid water molecules near an electrified interface is important to many disciplines of science and engineering. In this study, we applied an external gate potential to the silica/water interface via an electrolyte-insulator-semiconductor (EIS) junction to control the surface charging state. Without varying the ionic composition in water, the electrical gating allowed an efficient tuning of the interfacial charge density and field. Using the sum-frequency vibrational spectroscopy, we found a drastic enhancement of interfacial OH vibrational signals at high potential in weakly acidic water, which exceeded that from conventional bulk-silica/water interfaces even in strong basic solutions. Analysis of the spectra indicated that it was due to the alignment of liquid water molecules through the electric double layer, where the screening was weak because of the low ion density. Such a combination of strong field and weak screening demonstrates the unique tuning capability of the EIS scheme, and would allow us to investigate a wealth of phenomena at charged oxide/water interfaces.

Intermolecular Interactions at the Silica-Liquid Interface Modulate the Fermi Resonance Coupling in Surface Methanol

The buried solid/liquid interface between hydrophilic fused silica and binary solvent mixtures of acetonitrile (MeCN) and methanol (MeOH) was studied with vibrational sum-frequency generation (vSFG) spectroscopy. Our data showed that at high relative concentrations of methanol, the Fermi resonance peak in the vSFG spectrum is greatly suppressed, and it progressively gains intensity as methanol is diluted with perdeuterated acetonitrile. This phenomenon is quantified by the Fermi resonance coupling coefficient, W, extracted using a two-level model, as well as the experimental intensity ratio, R, of the methyl Fermi resonance band to that of the symmetric stretch. At a 1.0 MeOH mole fraction, W and R values were 10 ± 10 cm^-1 and 0.01 ± 0.02, respectively, whereas at a 0.1 mole fraction, W and R increased to 46 ± 4 cm^-1 and 0.43 ± 0.16, respectively. This indicates that solvation with acetonitrile effectively tunes the Fermi coupling of methanol vibrations at the silica/liquid interface.

Molecular Fingerprints of Hydrophobicity at Aqueous Interfaces from Theory and Vibrational Spectroscopies

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.

Size dependence of hydrophobic hydration at electrified gold/water interfaces

Hydrophobic hydration at metal/water interfaces actively contributes to the energetics of electrochemical reactions, e.g. [Formula: see text] and [Formula: see text] reduction, where small hydrophobic molecules are involved. In this work, constant applied potential molecular dynamics is employed to study hydrophobic hydration at a gold/water interface. We propose an adaptation of the Lum-Chandler-Weeks (LCW) theory to describe the free energy of hydrophobic hydration at the interface as a function of solute size and applied voltage. Based on this model we are able to predict the free energy cost of cavity formation at the interface directly from the free energy cost in the bulk plus an interface-dependent correction term. The interfacial water network contributes significantly to the free energy, yielding a preference for outer-sphere adsorption at the gold surface for ideal hydrophobes. We predict an accumulation of small hydrophobic solutes of sizes comparable to CO or [Formula: see text], while the free energy cost to hydrate larger hydrophobes, above 2.5-Å radius, is shown to be greater at the interface than in the bulk. Interestingly, the transition from the volume dominated to the surface dominated regimes predicted by the LCW theory in the bulk is also found to take place for hydrophobes at the Au/water interface but occurs at smaller cavity radii. By applying the adapted LCW theory to a simple model addition reaction, we illustrate some implications of our findings for electrochemical reactions.

Structure and sum-frequency generation spectra of water on neutral hydroxylated silica surfaces

Structural organization and vibrational sum-frequency generation (VSFG) spectra of water on crystalline and amorphous neutral silica surfaces were investigated by classical molecular dynamics simulations. The liquid phase represented with neat water and 1 M NaCl solution was analysed in terms of bonded interfacial layer (BIL), diffuse layer (DL) and bulk region. The simulations show that the structure of BIL depends on the surface morphology and density of surface OH groups. The water-silanol H-bond network and BIL structure are mainly insensitive to the presence of ions in the liquid phase. Molecules in DL of SiO₂/neat water interfaces preferentially orient their OH bonds towards the surfaces. This effect is directly related to an effective negative charge of formally neutral surfaces. Ions of the electrolyte solution affect the intermolecular structure in DL by screening the surface electric field and by the chaotropic effect. Calculated phase-sensitive VSFG (Im[χ⁽²⁾]) spectrum of BIL features low-frequency negative and high-frequency positive bands. Characteristics of the positive band reflect the strength of water-surface interactions and surface crystallinity, while the position and shape of the negative band are common to all interfaces. The Im[χ⁽²⁾] spectrum of DL is dominated by a contribution from the third-order χ⁽³⁾ susceptibility with the sign of the contribution directly related to the sign of electrostatic potential in the interfacial region. The DL spectrum is strongly affected by the presence of solvated ions. The computed intensity and Im[χ⁽²⁾] spectra of the amorphous silica/NaCl solution interface are in a good agreement with the conventional and phase-sensitive experimental VSFG spectra of fused SiO₂/water system at low pH, in contrast to the spectra of the amorphous silica/neat water interface. Origins of the discrepancy are discussed.

Evaluation of Molecular Polarizability and of Intensity Carrying Modes Contributions in Circular Dichroism Spectroscopies

We re-examine the theory of electronic and vibrational circular dichroism spectroscopy in terms of the formalism of frequency-dependent molecular polarizabilities. We show the link between Fermi’s gold rule in circular dichroism and the trace of the complex electric dipole–magnetic dipole polarizability. We introduce the C++ code polar to compute the molecular polarizability complex tensors from quantum chemistry outputs, thus simulating straightforwardly UV-visible absorption (UV-Vis)/electronic circular dichroism (ECD) spectra, and infrared (IR)/vibrational circular dichroism (VCD) spectra. We validate the theory and the code by referring to literature data of a large group of chiral molecules, showing the remarkable accuracy of density functional theory (DFT) methods. We anticipate the application of this methodology to the interpretation of vibrational spectra in various measurement conditions, even in presence of metal surfaces with plasmonic properties. Our theoretical developments aim, in the long run, at embedding the quantum-mechanical details of the chiroptical spectroscopic response of a molecule into the simulation of the electromagnetic field distribution at the surface of plasmonic devices. Such simulations are also instrumental to the interpretation of the experimental spectra measured from devices designed to enhance chiroptical interactions by the surface plasmon resonance of metal nanostructures.