First-Principles Framework to Compute Sum-Frequency Generation Vibrational Spectra of Semiconductors and Insulators

Learning Electronic Polarizations in Aqueous Systems

The polarization of periodically repeating systems is a discontinuous function of the atomic positions, a fact which seems at first to stymie attempts at their statistical learning. Two approaches to build models for bulk polarizations are compared: one in which a simple point charge model is used to preprocess the raw polarization to give a learning target that is a smooth function of atomic positions and the total polarization is learned as a sum of atom-centered dipoles and one in which instead the average position of Wannier centers around atoms is predicted. For a range of bulk aqueous systems, both of these methods perform perform comparatively well, with the former being slightly better but often requiring an extra effort to find a suitable point charge model. As a challenging test, we also analyze the performance of the models at the air-water interface. In this case, while the Wannier center approach delivers accurate predictions without further modifications, the preprocessing method requires augmentation with information from isolated water molecules to reach similar accuracy. Finally, we present a simple protocol to preprocess the polarizations in a data-driven way using a small number of derivatives calculated at a much lower level of theory, thus overcoming the need to find point charge models without appreciably increasing the computation cost. We believe that the training strategies presented here help the construction of accurate polarization models required for the study of the dielectric properties of realistic complex bulk systems and interfaces with ab initio accuracy.

The Role of Sum-Frequency Generation Spectroscopy in Understanding On-Surface Reactions and Dynamics in Atmospheric Model-Systems

Surfaces, both water/air and solid/water, play an important role in mediating a multitude of processes central to atmospheric chemistry, particularly in the aerosol phase. However, the study of both static and dynamic properties of surfaces is highly challenging from an experimental standpoint, leading to a lack of molecular level information about the processes that take place at these systems and how they differ from bulk. One of the few techniques that has been able to capture ultrafast surface phenomena is time-resolved sum-frequency generation (SFG) spectroscopy. Since it is both surface-specific and chemically sensitive, the extension of this spectroscopic technique to the time domain makes it possible to study dynamic processes on the femtosecond time scale. In this Perspective, we will explore recent advances made in the field both in terms of studying energy dissipation as well as chemical reactions and the role the surface geometry plays in these processes.

Fully First-Principles Surface Spectroscopy with Machine Learning

Our current understanding of the structure and dynamics of aqueous interfaces at the molecular level has grown substantially due to the continuous development of surface-specific spectroscopies, such as vibrational sum-frequency generation (VSFG). As in other vibrational spectroscopies, we must turn to atomistic simulations to extract all of the information encoded in the VSFG spectra. The high computational cost associated with existing methods means that they have limitations in representing systems with complex electronic structure or in achieving statistical convergence. In this work, we combine high-dimensional neural network interatomic potentials and symmetry-adapted Gaussian process regression to overcome these constraints. We show that it is possible to model VSFG signals with fully ab initio accuracy using machine learning and illustrate the versatility of our approach on the water/air interface. Our strategy allows us to identify the main sources of theoretical inaccuracy and establish a clear pathway toward the modeling of surface-sensitive spectroscopy of complex interfaces.

Deciphering the Properties of Nanoconfined Aqueous Solutions by Vibrational Sum Frequency Generation Spectroscopy

When confined between walls at nanometer distances, water exhibits surprisingly different properties with reference to bare interfacial water. Based on computer simulations, we demonstrate how vibrational sum frequency generation (VSFG) spectroscopy can be used-even with very mild symmetry breaking-to discriminate multilayer water in wide slit pores from both bilayer and monolayer water confined within molecularly narrow pores. Applying the technique, the VSFG lineshapes of monolayer, bilayer, and multilayer water are found to differ in characteristic ways, which is explained by their distinct density stratifications giving rise to different H-bonding patterns in the respective solvation layers.

Interfacial Coupling and Modulation of van der Waals Heterostructures for Nanodevices

In recent years, van der Waals heterostructures (vdWHs) of two-dimensional (2D) materials have attracted extensive research interest. By stacking various 2D materials together to form vdWHs, it is interesting to see that new and fascinating properties are formed beyond single 2D materials; thus, 2D heterostructures-based nanodevices, especially for potential optoelectronic applications, were successfully constructed in the past few decades. With the dramatically increased demand for well-controlled heterostructures for nanodevices with desired performance in recent years, various interfacial modulation methods have been carried out to regulate the interfacial coupling of such heterostructures. Here, the research progress in the study of interfacial coupling of vdWHs (investigated by Photoluminescence, Raman, and Pump-probe spectroscopies as well as other techniques), the modulation of interfacial coupling by applying various external fields (including electrical, optical, mechanical fields), as well as the related applications for future electrics and optoelectronics, have been briefly reviewed. By summarizing the recent progress, discussing the recent advances, and looking forward to future trends and existing challenges, this review is aimed at providing an overall picture of the importance of interfacial modulation in vdWHs for possible strategies to optimize the device's performance.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Sum frequency generation, calculation of absolute intensities, comparison with experiments, and two-field relaxation-based derivation

The experimental sum frequency generation (SFG) spectrum is the response to an infrared pulse and a visible pulse and is a highly surface-sensitive technique. We treat the surface dangling OH bonds at the air/water interface and focus on the absolute SFG intensities for the resonant terms, a focus that permits insight into the consequences of some approximations. For the polarization combinations, the calculated linewidths for the water interface dangling OH SFG band at 3,700 [Formula: see text] are, as usual, too large, because of the customary neglect of motional narrowing. The integrated spectrum is used to circumvent this problem and justified here using a Kubo-like formalism and theoretical integrated band intensities rather than peak intensities. Only relative SFG intensities are usually reported. The absolute integrated SFG intensities for three polarization combinations for sum frequency, visible, and infrared beams are computed. We use molecular dynamics and the dipole and the polarizability matrix elements obtained from infrared and Raman studies of [Formula: see text]O vapor. The theoretical expressions for two of the absolute susceptibilities contain only a single term and agree with experiment to about a factor of 1.3, with no adjustable parameters. The Fresnel factors are included in that comparison. One of the susceptibilities contains instead four positive and negative terms and agrees less well. The expression for the SFG correlation function is normally derived from a statistical mechanical formulation using a time-evolving density matrix. We show how a derivation based on a two-field relaxation leads to the same final result.

Hydrogen order at the surface of ice Ih revealed by vibrational spectroscopy

A combination of heterodyne-detected sum frequency generation spectroscopy and theoretical modeling elucidates that the surface of ice I_h at 100 K has hydrogen order with the OH group pointing upward to the air (“H-up” orientation).

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

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.

Coupling of Surface Chemistry and Electric Double Layer at TiO2 Electrochemical Interfaces

Surfaces of metal oxides at working conditions are usually electrified because of the acid-base chemistry. The charged interface compensated with counterions forms the so-called electric double layer. The coupling of surface chemistry and the electric double layer is considered to be crucial but is poorly understood because of the lack of information at the atomistic scale. Here, we used the latest development in density functional theory-based finite-field molecular dynamics simulation to investigate the pH dependence of the Helmholtz capacitance at electrified rutile TiO₂(110)-NaCl electrolyte interfaces. It is found that, because of competing forces from surface adsorption and from the electric double layer, water molecules have a stronger structural fluctuation at high pH, and this leads to a much larger capacitance. It is also seen that interfacial proton transfers at low pH increase significantly the capacitance value. These findings elucidate the microscopic origin of the same trend observed in titration experiments.

First-principles study on sum-frequency generation spectroscopy of methanol adsorbed on TiO2(110) surface: Effects of substrate and molecular coverages

In this work, starting from the general theory of sum-frequency generation (SFG), we proposed a computational strategy utilizing density functional theory with periodic boundary conditions to simulate the vibrational SFG of molecules/solid surface adsorption system. The method has been applied to the CH₃OH/TiO₂(110) system successfully. Compared with the isolated molecule model, our theoretical calculations showed that the TiO₂ substrate can significantly alter the second-order susceptibilities of a methanol molecule which is directly related to the SFG intensity. In addition, the SFG spectra have obvious changes while the methanol coverage increases, especially for the OH vibration peaks. Our theoretical spectra agree reasonably well with experimental measurements at 1 ML coverage, and an interesting peak which is absent in the theoretical spectra is tentatively assigned to some CH₃ stretch vibration of methanol adsorbed on the oxygen vacancy of TiO₂.

Ultrafast direct electron transfer at organic semiconductor and metal interfaces

The ability to control direct electron transfer can facilitate the development of new molecular electronics, light-harvesting materials, and photocatalysis. However, control of direct electron transfer has been rarely reported, and the molecular conformation-electron dynamics relationships remain unclear. We describe direct electron transfer at buried interfaces between an organic polymer semiconductor film and a gold substrate by observing the first dynamical electric field-induced vibrational sum frequency generation (VSFG). In transient electric field-induced VSFG measurements on this system, we observe dynamical responses (<150 fs) that depend on photon energy and polarization, demonstrating that electrons are directly transferred from the Fermi level of gold to the lowest unoccupied molecular orbital of organic semiconductor. Transient spectra further reveal that, although the interfaces are prepared without deliberate alignment control, a subensemble of surface molecules can adopt conformations for direct electron transfer. Density functional theory calculations support the experimental results and ascribe the observed electron transfer to a flat-lying polymer configuration in which electronic orbitals are found to be delocalized across the interface. The present observation of direct electron transfer at complex interfaces and the insights gained into the relationship between molecular conformations and electron dynamics will have implications for implementing novel direct electron transfer in energy materials.

Proton Order toward the Surface of Ice Ih Revealed by Heterodyne-Detected Sum Frequency Generation Spectroscopy

Using heterodyne-detected sum frequency generation (HD-SFG) spectroscopy, we investigated surface proton order at the basal, primary prism, and secondary prism faces of single-crystalline ice I_h at ca. 130 K. The complex phase of the obtained spectra clearly indicates that second-order nonlinear polarization from which the HD-SFG signal arises is generated exclusively at the surfaces. This suggests surface proton ordering along the normal, whereas the bulk remains proton-disordered, as is well known for ice I_h. A strong positive peak observed in the HD-SFG spectra enables us to determine the "direction" of the surface proton order as "H-up", that is, the hydrogen atom of the OH group pointing away from the bulk, irrespective of the ice faces. Reliable HD-SFG measurements carried out in the present study have greatly advanced our understanding of surface structure of ice I_h.

Excess Hydrogen Bond at the Ice-Vapor Interface around 200 K

Phase-resolved sum-frequency generation measurements combined with molecular dynamics simulations are employed to study the effect of temperature on the molecular arrangement of water on the basal face of ice. The topmost monolayer, interrogated through its nonhydrogen-bonded, free O-H stretch peak, exhibits a maximum in surface H-bond density around 200 K. This maximum results from two competing effects: above 200 K, thermal fluctuations cause the breaking of H bonds; below 200 K, the formation of bulklike crystalline interfacial structures leads to H-bond breaking. Knowledge of the surface structure of ice is critical for understanding reactions occurring on ice surfaces and ice nucleation.

Observation and Identification of a New OH Stretch Vibrational Band at the Surface of Ice

We study the signatures of the OH stretch vibrations at the basal surface of ice using heterodyne-detected sum-frequency generation and molecular dynamics simulations. At 150 K, we observe seven distinct modes in the sum-frequency response, five of which have an analogue in the bulk, and two pure surface-specific modes at higher frequencies (∼3530 and ∼3700 cm^-1). The band at ∼3530 cm^-1 has not been reported previously. Using molecular dynamics simulations, we find that the ∼3530 cm^-1 band contains contributions from OH stretch vibrations of both fully coordinated interfacial water molecules and water molecules with two donor and one acceptor hydrogen bond.

Modelling heterogeneous interfaces for solar water splitting

The generation of hydrogen from water and sunlight offers a promising approach for producing scalable and sustainable carbon-free energy. The key of a successful solar-to-fuel technology is the design of efficient, long-lasting and low-cost photoelectrochemical cells, which are responsible for absorbing sunlight and driving water splitting reactions. To this end, a detailed understanding and control of heterogeneous interfaces between photoabsorbers, electrolytes and catalysts present in photoelectrochemical cells is essential. Here we review recent progress and open challenges in predicting physicochemical properties of heterogeneous interfaces for solar water splitting applications using first-principles-based approaches, and highlights the key role of these calculations in interpreting increasingly complex experiments.

Sum Frequency Generation of Acetonitrile on a Rutile (110) Surface from Density Functional Theory-Based Molecular Dynamics

We present the calculation of vibrational sum frequency generation for molecules adsorbed on a semiconductor surface via density functional theory-based molecular dynamics. Using an efficient approach relying on the Gaussian and plane waves method and density functional perturbation theory, acetonitrile adsorbed on a rutile (110) surface has been studied, the vapor-solid interface of which has recently been investigated experimentally. Further analysis of the orientation of the acetonitrile molecules directly adsorbed on the rutile (110) surface agrees well with parameters derived from experiment. This opens a promising way for detailed study of semiconductor interfaces, which is of particular importance with respect to numerous applications such as, for instance, in materials design.

Dissecting the Molecular Structure of the Air/Water Interface from Quantum Simulations of the Sum-Frequency Generation Spectrum

The molecular characterization of the air/water interface is a key step in understanding fundamental multiphase phenomena ranging from heterogeneous chemical processes in the atmosphere to the hydration of biomolecules. The apparent simplicity of the air/water interface, however, masks an underlying complexity associated with the dynamic nature of the water hydrogen-bond network that has so far hindered an unambiguous characterization of its microscopic properties. Here, we demonstrate that the application of quantum many-body molecular dynamics, which enables spectroscopically accurate simulations of water from the gas to the condensed phase, leads to a definitive molecular-level picture of the interface region. For the first time, excellent agreement is obtained between the simulated vibrational sum-frequency generation spectrum and the most recent state-of-the-art measurements, without requiring any empirical frequency shift or ad hoc scaling of the spectral intensity. A systematic dissection of the spectral features demonstrates that a rigorous representation of nuclear quantum effects as well as of many-body energy and electrostatic contributions is necessary for a quantitative reproduction of the experimental data. The unprecedented accuracy of the simulations presented here indicates that quantum many-body molecular dynamics can enable predictive studies of aqueous interfaces, which by complementing analogous experimental measurements will provide unique molecular insights into multiphase and heterogeneous processes of relevance in chemistry, biology, materials science, and environmental research.