Spectroscopy and dynamics of the hydrated electron at the water/air interface.

The hydrated electron, e-(aq), has attracted much attention as a central species in radiation chemistry. However, much less is known about e-(aq) at the water/air surface, despite its fundamental role in electron transfer processes at interfaces. Using time-resolved electronic sum-frequency generation spectroscopy, the electronic spectrum of e-(aq) at the water/air interface and its dynamics are measured here, following photo-oxidation of the phenoxide anion. The spectral maximum agrees with that for bulk e-(aq) and shows that the orbital density resides predominantly within the aqueous phase, in agreement with supporting calculations. In contrast, the chemistry of the interfacial hydrated electron differs from that in bulk water, with e-(aq) diffusing into the bulk and leaving the phenoxyl radical at the surface. Our work resolves long-standing questions about e-(aq) at the water/air interface and highlights its potential role in chemistry at the ubiquitous aqueous interface.

Detection and Correction of Delocalization Errors for Electron and Hole Polarons Using Density-Corrected DFT.

Modeling polaron defects is an important aspect of computational materials science, but the description of unpaired spins in density functional theory (DFT) often suffers from delocalization error. To diagnose and correct the overdelocalization of spin defects, we report an implementation of density-corrected (DC-)DFT and its analytic energy gradient. In DC-DFT, an exchange-correlation functional is evaluated using a Hartree-Fock density, thus incorporating electron correlation while avoiding self-interaction error. Results for an electron polaron in models of titania and a hole polaron in Al-doped silica demonstrate that geometry optimization with semilocal functionals drives significant structural distortion, including the elongation of several bonds, such that subsequent single-point calculations with hybrid functionals fail to afford a localized defect even in cases where geometry optimization with the hybrid functional does localize the polaron. This has significant implications for traditional workflows in computational materials science, where semilocal functionals are often used for structure relaxation. DC-DFT calculations provide a mechanism to detect situations where delocalization error is likely to affect the results.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package.

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments.

The method of increments and frozen natural orbital (MI-FNO) framework is introduced to help expedite the application of noisy, intermediate-scale quantum (NISQ) devices for quantum chemistry simulations. The MI-FNO framework provides a systematic reduction of the occupied and virtual orbital spaces for quantum chemistry simulations. The correlation energies of the resulting increments from the MI-FNO reduction can then be solved by various algorithms, including quantum algorithms such as the phase estimation algorithm and the variational quantum eigensolver (VQE). The unitary coupled-cluster singles and doubles VQE framework is used to obtain correlation energies for the case of small molecules (i.e., BeH2, CH4, NH3, H2O, and HF) using the cc-pVDZ basis set. The quantum resource requirements are estimated for a constrained geometry complex catalyst that is utilized in industrial settings for the polymerization of α-olefins. We show that the MI-FNO approach provides a significant reduction in the quantum bit (qubit) requirements relative to the full system simulations. We propose that the MI-FNO framework can create scalable examples of quantum chemistry problems that are appropriate for assessing the progress of NISQ devices.

Quantum chemistry in arbitrary dielectric environments: Theory and implementation of nonequilibrium Poisson boundary conditions and application to compute vertical ionization energies at the air/water interface.

Widely used continuum solvation models for electronic structure calculations, including popular polarizable continuum models (PCMs), usually assume that the continuum environment is isotropic and characterized by a scalar dielectric constant, ε. This assumption is invalid at a liquid/vapor interface or any other anisotropic solvation environment. To address such scenarios, we introduce a more general formalism based on solution of Poisson's equation for a spatially varying dielectric function, ε(r). Inspired by nonequilibrium versions of PCMs, we develop a similar formalism within the context of Poisson's equation that includes the out-of-equilibrium dielectric response that accompanies a sudden change in the electron density of the solute, such as that which occurs in a vertical ionization process. A multigrid solver for Poisson's equation is developed to accommodate the large spatial grids necessary to discretize the three-dimensional electron density. We apply this methodology to compute vertical ionization energies (VIEs) of various solutes at the air/water interface and compare them to VIEs computed in bulk water, finding only very small differences between the two environments. VIEs computed using approximately two solvation shells of explicit water molecules are in excellent agreement with experiment for F-(aq), Cl-(aq), neat liquid water, and the hydrated electron, although errors for Li+(aq) and Na+(aq) are somewhat larger. Nonequilibrium corrections modify VIEs by up to 1.2 eV, relative to models based only on the static dielectric constant, and are therefore essential to obtain agreement with experiment. Given that the experiments (liquid microjet photoelectron spectroscopy) may be more sensitive to solutes situated at the air/water interface as compared to those in bulk water, our calculations provide some confidence that these experiments can indeed be interpreted as measurements of VIEs in bulk water.

The Hydrated Electron.

Existence of a hydrated electron as a byproduct of water radiolysis was established more than 50 years ago, yet this species continues to attract significant attention due to its role in radiation chemistry, including DNA damage, and because questions persist regarding its detailed structure. This work provides an overview of what is known in regards to the structure and spectroscopy of the hydrated electron, both in liquid water and in clusters [Formula: see text], the latter of which provide model systems for how water networks accommodate an excess electron. In clusters, the existence of both surface-bound and internally bound states of the excess electron has elicited much debate, whereas in bulk water there are questions regarding how best to understand the structure of the excess electron's spin density. The energetics of the equilibrium species e-(aq) and its excited states, in bulk water and at the air/water interface, are also addressed.

The Hydrated Electron at the Surface of Neat Liquid Water Appears To Be Indistinguishable from the Bulk Species.

Experiments have suggested that the aqueous electron, e(-)(aq), may play a significant role in the radiation chemistry of DNA. A recent measurement of the energy (below vacuum level) of the putative "interfacial" hydrated electron at the water/vacuum interface, performed using liquid microjet photoelectron spectroscopy, has been interpreted to suggest that aqueous electrons at the water/biomolecule interface may possess the appropriate energetics to induce DNA strand breaks, whereas e(-)(aq) in bulk water lies too far below the vacuum level to induce such reactions. Other such experiments, however, find no evidence of a long-lived feature at low binding energy. We employ a variety of computational strategies to demonstrate that the energetics of the hydrated electron at the surface of neat liquid water are not significantly different from those of e(-)(aq) in bulk water and as such are incompatible with dissociative electron attachment reactions in DNA. We furthermore suggest that no stable interfacial species may exist at all, consistent with the interpretation of certain surface-sensitive spectroscopy measurements, and that even if a short-lived, metastable species does exist at the vacuum/water interface, it would be extremely difficult to distinguish, experimentally, from e(-)(aq) in bulk water, using either optical absorption or photoelectron spectroscopy.

Optical spectroscopy of the bulk and interfacial hydrated electron from ab initio calculations.

The optical spectrum of the hydrated (aqueous) electron, e(aq)(­), is the primary observable by means of which this species is detected, monitored, and studied. In theoretical calculations, this spectrum has most often been simulated using one-electron models. Here, we present ab initio simulations of that spectrum in both bulk water and, for the first time, at the water/vapor interface, using density functional theory and its time-dependent variant. Our results indicate that this approach provides a reliable description, and quantitative agreement with the experimental spectrum for the bulk species is obtained using a "tuned" long-range corrected functional. The spectrum of the interfacial electron is found to be very similar to the bulk spectrum.