Surface stratification determines the interfacial water structure of simple electrolyte solutions.

The distribution of ions at the air/water interface plays a decisive role in many natural processes. Several studies have reported that larger ions tend to be surface-active, implying ions are located on top of the water surface, thereby inducing electric fields that determine the interfacial water structure. Here we challenge this view by combining surface-specific heterodyne-detected vibrational sum-frequency generation with neural network-assisted ab initio molecular dynamics simulations. Our results show that ions in typical electrolyte solutions are, in fact, located in a subsurface region, leading to a stratification of such interfaces into two distinctive water layers. The outermost surface is ion-depleted, and the subsurface layer is ion-enriched. This surface stratification is a key element in explaining the ion-induced water reorganization at the outermost air/water interface.

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.

First-Principles Simulations of Tip Enhanced Raman Scattering Reveal Active Role of Substrate on High-Resolution Images.

Tip-enhanced Raman scattering (TERS) has emerged as a powerful tool to obtain subnanometer spatial resolution fingerprints of atomic motion. Theoretical calculations that can simulate the Raman scattering process and provide an unambiguous interpretation of TERS images often rely on crude approximations of the local electric field. In this work, we present a novel and first-principles-based method to compute TERS images by combining Time Dependent Density Functional Theory (TD-DFT) and Density Functional Perturbation Theory (DFPT) to calculate Raman cross sections with realistic local fields. We present TERS results on free-standing benzene and C60 molecules, and on the TCNE molecule adsorbed on Ag(100). We demonstrate that chemical effects on chemisorbed molecules, often ignored in TERS simulations of larger systems, dramatically change the TERS images. This observation calls for the inclusion of chemical effects for predictive theory-experiment comparisons and an understanding of molecular motion at the nanoscale.

Is Unified Understanding of Vibrational Coupling of Water Possible? Hyper-Raman Measurement and Machine Learning Spectra.

The impact of the vibrational coupling of the OH stretch mode on the spectra differs significantly between IR and Raman spectra of water. Unified understanding of the vibrational couplings is not yet achieved. By using a different class of vibrational spectroscopy, hyper-Raman (HR) spectroscopy, together with machine-learning-assisted HR spectra calculation, we examine the impact of the vibrational couplings of water through the comparison of isotopically diluted H2O and pure H2O. We found that the isotopic dilution reduces the HR bandwidths, but the impact of the vibrational coupling is smaller than in the IR and parallel-polarized Raman. Machine learning HR spectra indicate that the intermolecular coupling plays a major role in broadening the bandwidth, while the intramolecular coupling is negligibly small, which is consistent with the IR and Raman spectra. Our result clearly demonstrates a limited impact of the intramolecular vibration, independent of the selection rules of vibrational spectroscopies.

Charge Transfer-Mediated Dramatic Enhancement of Raman Scattering upon Molecular Point Contact Formation.

Charge-transfer enhancement of Raman scattering plays a crucial role in current-carrying molecular junctions. However, the microscopic mechanism of light scattering in such nonequilibrium systems is still imperfectly understood. Here, using low-temperature tip-enhanced Raman spectroscopy (TERS), we investigate how Raman scattering evolves as a function of the gap distance in the single C60-molecule junction consisting of an Ag tip and various metal surfaces. Precise gap-distance control allows the examination of two distinct transport regimes, namely tunneling regime and molecular point contact (MPC). Simultaneous measurement of TERS and the electric current in scanning tunneling microscopy shows that the MPC formation results in dramatic Raman enhancement that enables one to observe the vibrations undetectable in the tunneling regime. This enhancement is found to commonly occur not only for coinage but also transition metal substrates. We suggest that the characteristic enhancement upon the MPC formation is rationalized by charge-transfer excitation.

Multidimensional Hydrogen Tunneling in Supported Molecular Switches: The Role of Surface Interactions.

The nuclear tunneling crossover temperature (T_{c}) of hydrogen transfer reactions in supported molecular-switch architectures can lie close to room temperature. This calls for the inclusion of nuclear quantum effects (NQEs) in the calculation of reaction rates even at high temperatures. However, computations of NQEs relying on standard parametrized dimensionality-reduced models quickly become inadequate in these environments. In this Letter, we study the paradigmatic molecular switch based on porphycene molecules adsorbed on metallic surfaces with full-dimensional calculations that combine density-functional theory for the electrons with the semiclassical ring-polymer instanton approximation for the nuclei. We show that the double intramolecular hydrogen transfer (DHT) rate can be enhanced by orders of magnitude due to surface fluctuations in the deep-tunneling regime. We also explain the origin of an Arrhenius temperature dependence of the rate below T_{c} and why this dependence differs at different surfaces. We propose a simple model to rationalize the temperature dependence of DHT rates spanning diverse fcc [110] surfaces.

Temperature dependence of the vibrational spectrum of porphycene: a qualitative failure of classical-nuclei molecular dynamics.

The temperature dependence of vibrational spectra can provide information about structural changes of a system and also serve as a probe to identify different vibrational mode couplings. Fully anharmonic temperature-dependent calculations of these quantities are challenging due to the cost associated with statistically converging trajectory-based methods, especially when accounting for nuclear quantum effects. Here, we train a high-dimensional neural network potential energy surface for the porphycene molecule based on data generated with DFT-B3LYP, including pairwise van der Waals interactions. In addition, we fit a kernel ridge regression model for the molecular dipole moment surface. The combination of this machinery with thermostatted path integral molecular dynamics (TRPMD) allows us to obtain well-converged, full-dimensional, fully-anharmonic vibrational spectra including nuclear quantum effects, without sacrificing the first-principles quality of the potential-energy surface or the dipole surface. Within this framework, we investigate the temperature and isotopologue dependence of the high-frequency vibrational fingerprints of porphycene. While classical-nuclei dynamics predicts a red shift of the vibrations encompassing the NH and CH stretches, TRPMD predicts a strong blue shift in the NH-stretch region and a smaller one in the CH-stretch region. We explain this behavior by analyzing the modulation of the effective potential with temperature, which arises from vibrational coupling between quasi-classical thermally activated modes and high-frequency quantized modes.

Elucidating the Nuclear Quantum Dynamics of Intramolecular Double Hydrogen Transfer in Porphycene.

We address the double hydrogen transfer (DHT) dynamics of the porphycene molecule, a complex paradigmatic system in which the making and breaking of H-bonds in a highly anharmonic potential energy surface require a quantum mechanical treatment not only of the electrons but also of the nuclei. We combine density functional theory calculations, employing hybrid functionals and van der Waals corrections, with recently proposed and optimized path-integral ring-polymer methods for the approximation of quantum vibrational spectra and reaction rates. Our full-dimensional ring-polymer instanton simulations show that below 100 K the concerted DHT tunneling pathway dominates but between 100 and 300 K there is a competition between concerted and stepwise pathways when nuclear quantum effects are included. We obtain ground-state reaction rates of 2.19 × 1011 s-1 at 150 K and 0.63 × 1011 s-1 at 100 K, in good agreement with experiment. We also reproduce the puzzling N-H stretching band of porphycene with very good accuracy from thermostated ring-polymer molecular dynamics simulations. The position and line shape of this peak, centered at around 2600 cm-1 and spanning 750 cm-1, stem from a combination of very strong H-bonds, the coupling to low-frequency modes, and the access to cis-like isomeric conformations, which cannot be appropriately captured with classical-nuclei dynamics. These results verify the appropriateness of our general theoretical approach and provide a framework for a deeper physical understanding of hydrogen transfer dynamics in complex systems.

Photophysics of Xanthene Dyes at High Concentrations in Solid Environments: Charge Transfer Assisted Triplet Formation.

The photophysical behavior of two xanthene dyes, Eosin Y and Phloxine B, included in microcrystalline cellulose particles is studied in a wide concentration range, with emphasis on the effect of dye concentration on fluorescence and triplet quantum yields. Absolute fluorescence quantum yields in the solid-state were determined by means of diffuse reflectance and steady-state fluorescence measurements, whereas absolute triplet quantum yields were obtained by laser-induced optoacoustic spectroscopy and their dependence on dye concentration was confirmed by diffuse reflectance laser flash photolysis and time-resolved phosphorescence measurements. When both quantum yields are corrected for reabsorption and reemission of radiation, ΦF values decrease strongly on increasing dye concentration, while a less pronounced decay is observed for ΦT . Fluorescence concentration quenching is attributed to the formation of dye aggregates or virtual traps resulting from molecular crowding. Dimeric traps are however able to generate triplet states. A mechanism based on the intermediacy of charge-transfer states is proposed and discussed. Calculation of parameters for photoinduced electron transfer between dye molecules within the traps evidences the feasibility of the proposed mechanism. Results demonstrate that photoactive energy traps, capable of yielding dye triplet states, can be formed even in highly-concentrated systems with random dye distributions.

Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature.

Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.

Quantum tunneling in real space: Tautomerization of single porphycene molecules on the (111) surface of Cu, Ag, and Au.

Tautomerization in single porphycene molecules is investigated on Cu(111), Ag(111), and Au(111) surfaces by a combination of low-temperature scanning tunneling microscopy (STM) experiments and density functional theory (DFT) calculations. It is revealed that the trans configuration is the thermodynamically stable form of porphycene on Cu(111) and Ag(111), whereas the cis configuration occurs as a meta-stable form. The trans → cis or cis → trans conversion on Cu(111) can be induced in an unidirectional fashion by injecting tunneling electrons from the STM tip or heating the surface, respectively. We find that the cis ↔ cis tautomerization on Cu(111) occurs spontaneously via tunneling, verified by the negligible temperature dependence of the tautomerization rate below ∼23 K. Van der Waals corrected DFT calculations are used to characterize the adsorption structures of porphycene and to map the potential energy surface of the tautomerization on Cu(111). The calculated barriers are too high to be thermally overcome at cryogenic temperatures used in the experiment and zero-point energy corrections do not change this picture, leaving tunneling as the most likely mechanism. On Ag(111), the reversible trans ↔ cis conversion occurs spontaneously at 5 K and the cis ↔ cis tautomerization rate is much higher than on Cu(111), indicating a significantly smaller tautomerization barrier on Ag(111) due to the weaker interaction between porphycene and the surface compared to Cu(111). Additionally, the STM experiments and DFT calculations reveal that tautomerization on Cu(111) and Ag(111) occurs with migration of porphycene along the surface; thus, the translational motion couples with the tautomerization coordinate. On the other hand, the trans and cis configurations are not discernible in the STM image and no tautomerization is observed for porphycene on Au(111). The weak interaction of porphycene with Au(111) is closest to the gas-phase limit and therefore the absence of trans and cis configurations in the STM images is explained either by the rapid tautomerization rate or the similar character of the molecular frontier orbitals of the trans and cis configurations.

Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures.

We present molecular dynamics simulation results describing proton/deuteron exchange equilibria along hydrogen bonds at the vicinity of HX acids (X = F, I) in aqueous clusters at low temperatures. To allow for an adequate description of proton transfer processes, our simulation scheme resorted on the implementation of a multistate empirical valence bond hamiltonian coupled to a path integral scheme to account for effects derived from nuclear quantum fluctuations. We focused attention on clusters comprising a number of water molecules close to the threshold values necessary to stabilize contact-ion-pairs. For X = F, our results reveal a clear propensity of the heavy isotope to lie at the bond bridging the halide to the nearest water molecule. Contrasting, for X = I, the thermodynamic stability is reversed and the former connectivity is preferentially articulated via the light isotope. These trends remain valid for undissociated and ionic descriptions of the stable valence bond states. The preferences are rationalized in terms of differences in the quantum kinetic energies of the isotopes which, in turn, reflect the extent of the local spatial confinements prevailing along the different hydrogen bonds in the clusters. In most cases, these features are also clearly reflected in the characteristics of the corresponding stretching bands of the simulated infrared spectra. This opens interesting possibilities to gauge the extent of the isotopic thermodynamic stabilizations and the strengths of the different hydrogen bonds by following the magnitudes and shifts of the spectral signals in temperature-controlled experiments, performed on mixed clusters combining H2O and HOD.

Tuning the concentration of dye loaded polymer films for maximum photosensitization efficiency: phloxine B in poly(2-hydroxyethyl methacrylate).

Fluorescence and singlet molecular oxygen ((1)O2) quantum yields for phloxine B loaded poly(2-hydroxyethyl methacrylate) thin films are determined at dye concentrations from 0.015 to 22 wt%. Fluorescence self-quenching and the fall off of the (1)O2 quantum yield observed above 0.1 wt% are attributed to very weakly interacting close-lying dye molecules acting as energy traps arising from molecular confinement. The maximum singlet oxygen generation efficiency (quantum yield × absorption factor) lies at concentrations around 2 wt%, where fluorescence self quenching amounts to more than 80%. Data are fitted quantitatively by using a quenching radius model involving energy migration and trapping with rQ = 1.2 nm. The present results constitute a proof of concept for the rational design of heterogeneous photosensitizers in general and, particularly, for applications in which the antimicrobial activity of singlet oxygen is central.

Effect of concentration on the formation of rose bengal triplet state on microcrystalline cellulose: a combined laser-induced optoacoustic spectroscopy, diffuse reflectance flash photolysis, and luminescence study.

Laser-induced optoacoustic spectroscopy (LIOAS), diffuse reflectance laser flash photolysis (DRLFP), and laser-induced luminescence (LIL) have been applied in conjunction to the determination of triplet state quantum yields of Rose Bengal (RB) supported on microcrystalline cellulose, a strongly light-scattering solid. Among the three used methods, the only one capable of providing absolute triplet quantum yields is LIOAS, but DRLFP and LIL aid in demonstrating that the LIOAS signal arises in fact from the triplet state and confirm the trend found with RB concentration. The coherence found for the three techniques demonstrates the usefulness of the approach. Observed triplet quantum yields are nearly constant within a limited concentration range, after which they decay strongly due to the generation of inactive dye aggregates or energy trapping centers. When quantum yields are divided by the fraction of absorbed light exciting the dye, the quotient falls off steadily with concentration, following the same trend as the observed fluorescence quantum yield. The conditions that maximize triplet formation are determined as a compromise between the rising light absorption and the decrease of quantum yield with RB concentration.