First-principle protocol for calculating ionization energies and redox potentials of solvated molecules and ions: theory and application to aqueous phenol and phenolate

Surface accumulation and acid-base equilibrium of phenol at the liquid-vapor interface

We have investigated the surfactant properties of phenol in aqueous solution as a function of pH and bulk concentration using liquid-jet photoelectron spectroscopy (LJ-PES) and surface tension measurements. The emphasis of this work is on the determination of the Gibbs free energy of adsorption and surface excess of phenol and its conjugate base phenolate at the bulk pKa (9.99), which can be determined for each species using photoelectron spectroscopy. These values are compared to those obtained in measurements well below and well above the pKa, where pure phenol or phenolate, respectively, are the dominant species, and where the Gibbs free energy of adsorption determined from surface tension and LJ-PES data are in excellent agreement. At the bulk pKa the surface-sensitive LJ-PES measurements show a deviation of the expected phenol/phenolate ratio in favor of phenol, i.e., an apparent upward shift of the at the surface. In addition, the Gibbs free energies of adsorption determined by LJ-PES at the bulk pKa for phenol and phenolate deviate from those observed for the pure solutions. We discuss these observations in view of the different surface propensity of phenol and phenolate as well as potential cooperative interactions between them in the near-surface region.

Valence Electronic Structure of Interfacial Phenol in Water Droplets

Biochemistry and a large part of atmospheric chemistry occur in aqueous environments or at aqueous interfaces, where (photo)chemical reaction rates can be increased by up to several orders of magnitude. The key to understanding the chemistry and photoresponse of molecules in and "on" water lies in their valence electronic structure, with a sensitive probe being photoelectron spectroscopy. This work reports velocity-map photoelectron imaging of submicrometer-sized aqueous phenol droplets in the valence region after nonresonant (288 nm) and resonance-enhanced (274 nm) two-photon ionization with femtosecond ultraviolet light, complementing previous liquid microjet studies. For nonresonant photoionization, our concentration-dependent study reveals a systematic decrease in the vertical binding energy (VBE) of aqueous phenol from 8.0 ± 0.1 eV at low concentration (0.01 M) to 7.6 ± 0.1 eV at high concentration (0.8 M). We attribute this shift to a systematic lowering of the energy of the lowest cationic state with increasing concentration caused by the phenol dimer and aggregate formation at the droplet surface. Contrary to nonresonant photoionization, no significant concentration dependence of the VBE was observed for resonance-enhanced photoionization. We explain the concentration-independent VBE of ∼8.1 eV observed upon resonant ionization by ultrafast intermediate state relaxation and changes in the accessible Franck-Condon region as a consequence of the lowering of the intermediate state potential energy due to the formation of phenol excimers and excited phenol aggregates. Correcting for the influence of electron transport scattering in the droplets reduced the measured VBEs by 0.1-0.2 eV.

Machine Learning Approach to Vertical Energy Gap in Redox Processes

A straightforward approach to calculating the free energy change (ΔG) and reorganization energy of a redox process is linear response approximation (LRA). However, accurate prediction of redox properties is still challenging due to difficulties in conformational sampling and vertical energy-gap sampling. Expensive hybrid quantum mechanical/molecular mechanical (QM/MM) calculations are typically employed in sampling energy gaps using conformations from simulations. To alleviate the computational cost associated with the expensive QM method in the QM/MM calculation, we propose machine learning (ML) methods to predict the vertical energy gaps (VEGs). We tested several ML models to predict the VEGs and observed that simple models like linear regression show excellent performance (mean absolute error ∼0.1 eV) in predicting VEGs in all test systems, even when using features extracted from cheaper semiempirical methods. Our best ML model (extra trees regressor) shows a mean absolute error of around 0.1 eV while using features from the cheapest QM method. We anticipate our approach can be generalized to larger macromolecular systems with more complex redox centers.

Liquid-jet photoemission spectroscopy as a structural tool: site-specific acid-base chemistry of vitamin C

Liquid-jet photoemission spectroscopy (LJ-PES) directly probes the electronic structure of solutes and solvents. It also emerges as a novel tool to explore chemical structure in aqueous solutions, yet the scope of the approach has to be examined. Here, we present a pH-dependent liquid-jet photoelectron spectroscopic investigation of ascorbic acid (vitamin C). We combine core-level photoelectron spectroscopy and ab initio calculations, allowing us to site-specifically explore the acid-base chemistry of the biomolecule. For the first time, we demonstrate the capability of the method to simultaneously assign two deprotonation sites within the molecule. We show that a large change in chemical shift appears even for atoms distant several bonds from the chemically modified group. Furthermore, we present a highly efficient and accurate computational protocol based on a single structure using the maximum-overlap method for modeling core-level photoelectron spectra in aqueous environments. This work poses a broader question: to what extent can LJ-PES complement established structural techniques such as nuclear magnetic resonance? Answering this question is highly relevant in view of the large number of incorrect molecular structures published.

First and Second Reductions in an Aprotic Solvent: Comparing Computational and Experimental One-Electron Reduction Potentials for 345 Quinones

Using reference reduction potentials of quinones recently measured relative to the saturated calomel electrode (SCE) in N,N-dimethylformamide (DMF), we benchmark absolute one-electron reduction potentials computed for 345 Q/Q•- and 265 Q•-/Q2- half-reactions using adiabatic electron affinities computed with density functional theory and solvation energies computed with four continuum solvation models: IEF-PCM, C-PCM, COSMO, and SM12. Regression analyses indicate a strong linear correlation between experimental and absolute computed Q/Q•- reduction potentials with Pearson's correlation coefficient (r) between 0.95 and 0.96 and the mean absolute error (MAE) relative to the linear fit between 83.29 and 89.51 mV for different solvation methods when the slope of the regression is constrained to 1. The same analysis for Q•-/Q2- gave a linear regression with r between 0.74 and 0.90 and MAE between 95.87 and 144.53 mV, respectively. The y-intercept values obtained from the linear regressions are in good agreement with the range of absolute reduction potentials reported in the literature for the SCE but reveal several sources of systematic error. The y-intercepts from Q•-/Q2- calculations are lower than those from Q/Q•- by around 320-410 mV for IEF-PCM, C-PCM, and SM12 compared to 210 mV for COSMO. Systematic errors also arise between molecules having different ring sizes (benzoquinones, naphthoquinones, and anthraquinones) and different substituents (titratable vs nontitratable). SCF convergence issues were found to be a source of random error that was slightly reduced by directly optimizing the solute structure in the continuum solvent reaction field. While SM12 MAEs were lower than those of the other solvation models for Q/Q•-, SM12 had larger MAEs for Q•-/Q2- pointing to a larger error when describing multiply charged anions in DMF. Altogether, the results highlight the advantages of, and further need for, testing computational methods using a large experimental data set that is not skewed (e.g., having more titratable than nontitratable substituents on different parent groups or vice versa) to help further distinguish between sources of random and systematic errors in the calculations.

Probing Antioxidant-Related Properties for Phenolic Compounds

In this work, properties related to antioxidant-potential mechanisms (such as the bond dissociation enthalpy, BDE, for the homolytic cleavage of the O-H bond and ionization energies, IEs) were determined for phenol, pyrocatechol, and gallic acid (GA). Both the protonated and deprotonated forms of GA were investigated. The Feller-Peterson-Dixon (FPD) composite method was employed with a variety of computational approaches, i.e., density functional theory, Möller-Plesset perturbation theory, and coupled-cluster-based methods, in combination with large correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. FPD results were compared to experimental and computational data available in the literature, presenting good agreement. For example, the FPD BDE (298 K) obtained for phenol, which was based on valence-correlated MP2/CBS calculations with contributions from correlating all electrons, was determined to be 87.56 kcal/mol, a value that is 0.42 kcal/mol lower than the result obtained in the most recent experiments, 87.98 ± 0.62. Calibration against coupled-cluster calculations was also carried out for phenol. We expect that the outcomes gathered here may help in establishing a general protocol for computational chemists that are interested in determining antioxidant-related properties for phenolic compounds with considerable accuracy as well as to motivate future IE measurements (particularly for GA) to be accomplished in the near future.

Voltage prediction of vanadium redox flow batteries from first principles

Global energy demand has been increasing for decades, which has created a necessity for large scale energy storage solutions for renewable energy sources. We studied the voltage of vanadium redox flow batteries (VRFBs) with density functional theory (DFT) and a newly developed technique usingab initiomolecular dynamics (AIMD). DFT was used to create cluster models to calculate the voltage of VRFBs. However, DFT is not suited for capturing the dynamics and interactions in a liquid electrolyte, leading to the need for AIMD, which is capable of accurately modeling such things. The molarities and densities of all systems were carefully considered to match experimental conditions. With the use of AIMD, we calculated a voltage of 1.23 V, which compares well with the experimental value of 1.26 V. The techniques developed using AIMD for voltage calculations will be useful for the investigation of potential future battery technologies or as a screening process for additives to make improvements to currently available batteries.

Ultrafast Vibrational Wave Packet Dynamics of the Aqueous Guanine Radical Anion Induced by Photodetachment

Studying the ultrafast dynamics of ionized aqueous biomolecules is important for gaining an understanding of the interaction of ionizing radiation with biological matter. Guanine plays an essential role in biological systems as one of the four nucleobases that form the building blocks of deoxyribonucleic acid (DNA). Guanine radicals can induce oxidative damage to DNA, particularly due to the lower ionization potential of guanine compared to the other nucleobases, sugars, and phosphate groups that are constituents of DNA. This study utilizes femtosecond optical pump-probe spectroscopy to observe the ultrafast vibrational wave packet dynamics of the guanine radical anion launched by photodetachment of the aqueous guanine dianion. The vibrational wave packet motion is resolved into 11 vibrational modes along which structural reorganization occurs upon photodetachment. These vibrational modes are assigned with the aid of density functional theory (DFT) calculations. Our work sheds light on the ultrafast vibrational dynamics following the ionization of nucleobases in an aqueous medium.

Core-ionization spectrum of liquid water

We present state-of-the-art calculations of the core-ionization spectrum of water. Despite significant progress in procedures developed to mitigate various experimental complications and uncertainties, the experimental determination of ionization energies of solvated species involves several non-trivial steps such as assessing the effect of the surface potential, electrolytes, and finite escape depths of photoelectrons. This provides a motivation to obtain robust theoretical values of the intrinsic bulk ionization energy and the corresponding solvent-induced shift. Here we develop theoretical protocols based on coupled-cluster theory and electrostatic embedding. Our value of the intrinsic solvent-induced shift of the 1sO ionization energy of water is -1.79 eV. The computed absolute position and the width of the 1sO peak in photoelectron spectrum of water are 538.47 eV and 1.44 eV, respectively, agreeing well with the best experimental values.

Ultrafast dynamics in polymeric carbon nitride thin films probed by time-resolved EUV photoemission and UV-Vis transient absorption spectroscopy

The ground- and excited-state electronic structures of four polymeric carbon nitride (PCN) materials have been investigated using a combination of photoemission and optical absorption spectroscopy. To establish the driving forces for photocatalytic water-splitting reactions, the ground-state data was used to produce a band diagram of the PCN materials and the triethanolamine electron scavenger, commonly implemented in water-splitting devices. The ultrafast charge-carrier dynamics of the same PCN materials were also investigated using two femtosecond-time-resolved pump-probe techniques: extreme-ultraviolet (EUV) photoemission and ultraviolet-visible (UV-Vis) transient absorption spectroscopy. The complementary combination of these surface- and bulk-sensitive methods facilitated photoinduced kinetic measurements spanning the sub-picosecond to few nanosecond time range. The results show that 400 nm (3.1 eV) excitation sequentially populates a pair of short-lived transient species, which subsequently produce two different long-lived excited states on a sub-picosecond time scale. Based on the spectro-temporal characteristics of the long-lived signals, they are assigned to singlet-exciton and charge-transfer states. The associated charge-separation efficiency was inferred to be between 65% and 78% for the different studied materials. A comparison of results from differently synthesized PCNs revealed that the early-time processes do not differ qualitatively between sample batches, but that materials of more voluminous character tend to have higher charge separation efficiencies, compared to exfoliated colloidal materials. This finding was corroborated via a series of experiments that revealed an absence of any pump-fluence dependence of the initial excited-state decay kinetics and characteristic carrier-concentration effects that emerge beyond few-picosecond timescales. The initial dynamics of the photoinduced charge carriers in the PCNs are correspondingly determined to be spatially localised in the immediate vicinity of the lattice-constituting motif, while the long-time behaviour is dominated by charge-transport and recombination processes. Suppressing the latter by confining excited species within nanoscale volumes should therefore affect the usability of PCN materials in photocatalytic devices.

Rapid and selective colorimetric determination of L-DOPA in human serum with silver nanoparticles

L-DOPA, or l-3,4-dihydroxyphenylalanine is an aromatic amino acid, which plays a significant role in human metabolism as a precursor of important neurotransmitters. We develop a fast and simple colorimetric method for the detection of L-DOPA in biological fluids. The method is based on the reduction of silver ions with L-DOPA and the subsequent formation of L-DOPA stabilized silver nanoparticles (Ag NPs). In this novel approach, L-DOPA works as both reducing and stabilizing agent, which provides selectivity and simplifies the procedure. HR-TEM images show very narrow Ag NPs distribution with an average size of 24 nm. Such sensor design is suggested for the first time. We also calculate vertical ionization potential, vertical electron affinity, and Gibbs free energy change of different ionic forms of L-DOPA and amino acids at the M06-2X/def2-TZVP level for the gas phase in comparison with that of silver. A model of silver ions reduction by aromatic amino acids is proposed: the ionic forms with charge -1 are suggested to reduce silver ions. High selectivity against aromatic amino acids, dopamine and serotonin is achieved by tuning pH and involving two L-DOPA forms with charged both hydroxyphenolate and carboxylate groups in the stabilization of uniform-sized Ag NPs. The method is applicable for the determination of L-DOPA in human serum with the 50 nM limit of detection and the linear range up to 5 μM. Ag NPs formation and coloring the solution proceeds in a few minutes. The suggested colorimetric method has potential application in clinical trials.

Accurate Vertical Electron Detachment Energies and Multiphoton Resonant Photoelectron Spectra of Biochromophore Anions in Aqueous Solution

We introduce a new methodology for calculating vertical electron detachment energies (VDEs) of biologically relevant chromophores in their deprotonated anionic forms in aqueous solution. It combines a large-scale mixed DFT/EFP/MD approach with the high-level multireference perturbation theory XMCQDPT2 and the Effective Fragment Potential (EFP) method. The methodology includes a multiscale flexible treatment of inner (∼1000 water molecules) and outer (∼18000 water molecules) water shells around a charged solute, capturing both the effects of specific solvation and the properties of bulk water. VDEs are calculated as a function of system size for getting a converged value at the DFT/EFP level of theory. The XMCQDPT2/EFP approach, adapted for calculating VDEs, supports the DFT/EFP results. When corrected for a solvent polarization contribution, the XMCQDPT2/EFP method yields the most accurate estimate to date of the first VDE for aqueous phenolate (7.3 ± 0.1 eV), which agrees well with liquid-jet X-ray photoelectron spectroscopy data (7.1 ± 0.1 eV). We show that the geometry of the water shell and its size are essential for accurate VDE calculations of aqueous phenolate and its biologically relevant derivatives. By simulating photoelectron spectra of aqueous phenolate upon two-photon excitation at wavelengths resonant with the S₀ → S₁ transition, we also provide interpretation of recent multiphoton UV liquid-microjet photoelectron spectroscopy experiments. We show that its first VDE is consistent with our estimate of 7.3 eV, when experimental two-photon binding energies are corrected for the resonant contribution.

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.

TDDFT versus GW/BSE Methods for Prediction of Light Absorption and Emission in a TADF Emitter

Design concepts for organic light emitting diode (OLED) emitters, which exhibit thermally activated delayed fluorescence (TADF) and thereby achieve quantum yields exceeding 25%, depend on singlet-triplet splitting energies of order kT to allow reverse intersystem crossing at ambient temperatures. Simulation methods for these systems must be able to treat relatively large organic molecules, as well as predict their excited state energies, transition energies, singlet-triplet splittings, and absorption and emission cross sections with reasonable accuracy, in order to prove useful in the design process. Here we compare predictions of TDDFT with M06-2X and ωB97X-D exchange-correlation functionals and a G_oW_o@HF/BSE method for these quantities in the well-studied DPTZ-DBTO2 TADF emitter molecule. Geometry optimization is performed for ground state (GS) and lowest donor-acceptor charge transfer (CT) state for each functional. Optical absorption and emission cross sections and energies are calculated at these geometries. Relaxation energies are on the order of 0.5 eV, and the importance of obtaining excited state equilibrium geometries in predicting delayed fluorescence is demonstrated. There are clear trends in predictions of G_oW_o@HF/BSE, and TDDFT/ωB97X-D and M06-2X methods in which the former method favors local exciton (LE) states while the latter favors DA CT states and ωB97X-D makes intermediate predictions. G_oW_o@HF/BSE suffers from triplet instability for LE states but not CT states relevant for TADF. Shifts in HOMO and LUMO levels on adding a conductor-like polarizable continuum model dielectric background are used to estimate changes in excitation energies on going from the gas phase to a solvated molecule.

UV Photoelectron Spectroscopy of Aqueous Solutions

Knowledge of the electronic structure of an aqueous solution is a prerequisite to understanding its chemical and biological reactivity and its response to light. One of the most direct ways of determining electronic structure is to use photoelectron spectroscopy to measure electron binding energies. Initially, photoelectron spectroscopy was restricted to the gas or solid phases due to the requirement for high vacuum to minimize inelastic scattering of the emitted electrons. The introduction of liquid-jets and their combination with intense X-ray sources at synchrotrons in the late 1990s expanded the scope of photoelectron spectroscopy to include liquids. Liquid-jet photoelectron spectroscopy is now an active research field involving a growing number of research groups. A limitation of X-ray photoelectron spectroscopy of aqueous solutions is the requirement to use solutes with reasonably high concentrations in order to obtain photoelectron spectra with adequate signal-to-noise after subtracting the spectrum of water. This has excluded most studies of organic molecules, which tend to be only weakly soluble. A solution to this problem is to use resonance-enhanced photoelectron spectroscopy with ultraviolet (UV) light pulses (hν ≲ 6 eV). However, the development of UV liquid-jet photoelectron spectroscopy has been hampered by a lack of quantitative understanding of inelastic scattering of low kinetic energy electrons (≲5 eV) and the impact on spectral lineshapes and positions.In this Account, we describe the key steps involved in the measurement of UV photoelectron spectra of aqueous solutions: photoionization/detachment, electron transport of low kinetic energy electrons through the conduction band, transmission through the water-vacuum interface, and transport through the spectrometer. We also explain the steps we take to record accurate UV photoelectron spectra of liquids with excellent signal-to-noise. We then describe how we have combined Monte Carlo simulations of electron scattering and spectral inversion with molecular dynamics simulations of depth profiles of organic solutes in aqueous solution to develop an efficient and widely applicable method for retrieving true UV photoelectron spectra of aqueous solutions. The huge potential of our experimental and spectral retrieval methods is illustrated using three examples. The first is a measurement of the vertical detachment energy of the green fluorescent protein chromophore, a sparingly soluble organic anion whose electronic structure underpins its fluorescence and photooxidation properties. The second is a measurement of the vertical ionization energy of liquid water, which has been the subject of discussion since the first X-ray photoelectron spectroscopy measurement in 1997. The third is a UV photoelectron spectroscopy study of the vertical ionization energy of aqueous phenol which demonstrates the possibility of retrieving true photoelectron spectra from measurements with contributions from components with different concentration profiles.

Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase

Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis. Theoretical methods can provide such quantification but are often hindered by the innate complexity of electron correlation and dynamic relaxation of solvent environments. We present and validate a protocol for calculating the redox potentials of the ferrocene/ferrocenium redox pair in acetonitrile. Equation-of-motion and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in the ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrates the ability of modern computational methods to predict absolute redox potentials quantitatively and to quantify the correlation of dynamic effects, which underlie their origin.

Cysteine Oxidation in Proteins: Structure, Biophysics, and Simulation

Cysteine side chains can exist in distinct oxidation states depending on the pH and redox potential of the environment, and cysteine oxidation plays important yet complex regulatory roles. Compared with the effects of post-translational modifications such as phosphorylation, the effects of oxidation of cysteine to sulfenic, sulfinic, and sulfonic acid on protein structure and function remain relatively poorly characterized. We present an analysis of the role of cysteine reactivity as a regulatory factor in proteins, emphasizing the interplay between electrostatics and redox potential as key determinants of the resulting oxidation state. A review of current computational approaches suggests underdeveloped areas of research for studying cysteine reactivity through molecular simulations.

Accurate Vertical Ionization Energy of Water and Retrieval of True Ultraviolet Photoelectron Spectra of Aqueous Solutions

Ultraviolet (UV) photoelectron spectroscopy provides a direct way of measuring valence electronic structure; however, its application to aqueous solutions has been hampered by a lack of quantitative understanding of how inelastic scattering of low-energy (<5 eV) electrons in liquid water distorts the measured electron kinetic energy distributions. Here, we present an efficient and widely applicable method for retrieving true UV photoelectron spectra of aqueous solutions. Our method combines Monte Carlo simulations of electron scattering and spectral inversion, with molecular dynamics simulations of depth profiles of organic solutes in aqueous solution. Its application is demonstrated for both liquid water, and aqueous solutions of phenol and phenolate, which are ubiquitous biologically relevant structural motifs.

Ultrafast proton transfer of the aqueous phenol radical cation

Proton transfer (PT) reactions are fundamental to numerous chemical and biological processes. While sub-picosecond PT involving electronically excited states has been extensively studied, little is known about ultrafast PT triggered by photoionization. Here, we employ femtosecond optical pump-probe spectroscopy and quantum dynamics calculations to investigate the ultrafast proton transfer dynamics of the aqueous phenol radical cation (PhOH˙⁺). Analysis of the vibrational wave packet dynamics reveals unusually short dephasing times of 0.18 ± 0.02 ps and 0.16 ± 0.02 ps for the PhOH˙⁺ O-H wag and bend frequencies, respectively, suggestive of ultrafast PT occurring on the ∼0.1 ps timescale. The reduced potential energy surface obtained from ab initio calculations shows that PT is barrierless when it is coupled to the intermolecular hindered translation between PhOH˙⁺ and the proton-acceptor water molecule. Quantum dynamics calculations yield a lifetime of 193 fs for PhOH˙⁺, in good agreement with the experimental results and consistent with the PT reaction being mediated by the intermolecular O⋯O stretch. These results suggest that photoionization can be harnessed to produce photoacids that undergo ultrafast PT. In addition, they also show that PT can serve as an ultrafast deactivation channel for limiting the oxidative damage potential of radical cations.

DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding

We present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF4 ionic liquid were considered as ``simple' and ``complex' solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.

Why the Photochemical Reaction of Phenol Becomes Ultrafast at the Air-Water Interface: The Effect of Surface Hydration

Photochemical reactions at the air-water interface can show remarkably different rates from those in bulk water. The present study elucidates the reaction mechanism of phenol characteristic at the air-water interface by the combination of molecular dynamics simulation and quantum chemical calculations of the excited states. We found that incomplete hydrogen bonding to phenol at the air-water interface affects excited states associated with the conical intersection and significantly reduces the reaction barrier, resulting in the distinctively facilitated rate in comparison with the bulk phase. The present study indicates that the reaction dynamics can be substantially different at the interfaces in general, reflecting the difference in the stabilization energy of the electronic states in markedly different solvation at the interface.

Probing aqueous ions with non-local Auger relaxation

Non-local analogues of Auger decay are increasingly recognized as important relaxation processes in the condensed phase. Here, we explore non-local autoionization, specifically Intermolecular Coulombic Decay (ICD), of a series of aqueous-phase isoelectronic cations following 1s core-level ionization. In particular, we focus on Na+, Mg2+, and Al3+ ions. We unambiguously identify the ICD contribution to the K-edge Auger spectrum. The different strength of the ion-water interactions is manifested by varying intensities of the respective signals: the ICD signal intensity is greatest for the Al3+ case, weaker for Mg2+, and absent for weakly-solvent-bound Na+. With the assistance of ab initio calculations and molecular dynamics simulations, we provide a microscopic understanding of the non-local decay processes. We assign the ICD signals to decay processes ending in two-hole states, delocalized between the central ion and neighbouring water. Importantly, these processes are shown to be highly selective with respect to the promoted water solvent ionization channels. Furthermore, using a core-hole-clock analysis, the associated ICD timescales are estimated to be around 76 fs for Mg2+ and 34 fs for Al3+. Building on these results, we argue that Auger and ICD spectroscopy represents a unique tool for the exploration of intra- and inter-molecular structure in the liquid phase, simultaneously providing both structural and electronic information.

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments.

Spectroscopic observation and ultrafast coherent vibrational dynamics of the aqueous phenylalanine radical

The phenylalanine radical (Phe˙) has been proposed to mediate biological electron transport (ET) and exhibit long-lived electronic coherences following attosecond photoionization. However, the coupling of ultrafast structural reorganization to the oxidation/ionization of biomolecules such as phenylalanine remains unexplored. Moreover, studies of ET involving Phe˙ are hindered by its hitherto unobserved electronic spectrum. Here, we report the spectroscopic observation and coherent vibrational dynamics of aqueous Phe˙, prepared by sub-6 fs photodetachment of phenylalaninate anions. Sub-picosecond transient absorption spectroscopy reveals the ultraviolet absorption signature of Phe˙. Ultrafast structural reorganization drives coherent vibrational motion involving nine fundamental frequencies and one overtone. DFT calculations rationalize the absence of the decarboxylation reaction, a photodegradation pathway previously identified for Phe˙. Our findings guide the interpretation of future attosecond experiments aimed at elucidating coherent electron motion in photoionized aqueous biomolecules and pave way for the spectroscopic identification of Phe˙ in studies of biological ET.

Bridging the Experiment-Calculation Divide: Machine Learning Corrections to Redox Potential Calculations in Implicit and Explicit Solvent Models

Prediction of redox potentials is essential for catalysis and energy storage. Although density functional theory (DFT) calculations have enabled rapid redox potential predictions for numerous compounds, prominent errors persist compared to experimental measurements. In this work, we develop machine learning (ML) models to reduce the errors of redox potential calculations in both implicit and explicit solvent models. Training and testing of the ML correction models are based on the diverse ROP313 data set with experimental redox potentials measured for organic and organometallic compounds in a variety of solvents. For the implicit solvent approach, our ML models can reduce both the systematic bias and the number of outliers. ML corrected redox potentials also demonstrate less sensitivity to DFT functional choice. For the explicit solvent approach, we significantly reduce the computational costs by embedding the microsolvated cluster in implicit bulk solvent, obtaining converged redox potential results with a smaller solvation shell. This combined implicit-explicit solvent model, together with GPU-accelerated quantum chemistry methods, enabled rapid generation of a large data set of explicit-solvent-calculated redox potentials for 165 organic compounds, allowing detailed investigation of the error sources in explicit solvent redox potential calculations.

Ultrafast vibrational wave packet dynamics of the aqueous tyrosyl radical anion induced by photodetachment

The ultrafast dynamics triggered by the photodetachment of the tyrosinate dianion in aqueous environment shed light on the elementary processes that accompany the interaction of ionizing radiation with biological matter. Photodetachment of the tryosinate dianion yields the tyrosyl radical anion, an important intermediate in biological redox reactions, although the study of its ultrafast dynamics is limited. Here, we utilize femtosecond optical pump-probe spectroscopy to investigate the ultrafast structural reorganization dynamics that follow the photodetachment of the tyrosinate dianion in aqueous solution. Photodetachment of the tyrosinate dianion leads to vibrational wave packet motion along seven vibrational modes that are coupled to the photodetachment process. The vibrational modes are assigned with the aid of density functional theory (DFT) calculations. Our results offer a glimpse of the elementary dynamics of ionized biomolecules and suggest the possibility of extending this approach to investigate the ionization-induced structural rearrangement of other aromatic amino acids and larger biomolecules.

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.

Following in Emil Fischer's Footsteps: A Site-Selective Probe of Glucose Acid-Base Chemistry

Liquid-jet photoelectron spectroscopy was applied to determine the first acid dissociation constant (pKa) of aqueous-phase glucose while simultaneously identifying the spectroscopic signature of the respective deprotonation site. Valence spectra from solutions at pH values below and above the first pKa reveal a change in glucose's lowest ionization energy upon the deprotonation of neutral glucose and the subsequent emergence of its anionic counterpart. Site-specific insights into the solution-pH-dependent molecular structure changes are also shown to be accessible via C 1s photoelectron spectroscopy. The spectra reveal a considerably lower C 1s binding energy of the carbon site associated with the deprotonated hydroxyl group. The occurrence of photoelectron spectral fingerprints of cyclic and linear glucose prior to and upon deprotonation are also discussed. The experimental data are interpreted with the aid of electronic structure calculations. Our findings highlight the potential of liquid-jet photoelectron spectroscopy to act as a site-selective probe of the molecular structures that underpin the acid-base chemistry of polyprotic systems with relevance to environmental chemistry and biochemistry.

The redox potential of a heme cofactor in Nitrosomonas europaea cytochrome c peroxidase: a polarizable QM/MM study

Redox reactions are crucial to biological processes that protect organisms against oxidative stress. Metalloenzymes, such as peroxidases which reduce excess reactive oxygen species into water, play a key role in detoxification mechanisms. Here we present the results of a polarizable QM/MM study of the reduction potential of the electron transfer heme in the cytochrome c peroxidase of Nitrosomonas europaea. We have found that environment polarization does not substantially affect the computed value of the redox potential. Particular attention has been given to analyzing the role of electrostatic interactions within the protein environment and the solvent on tuning the redox potential of the heme co-factor. We have found that the electrostatic interactions predominantly explain the fluctuations of the vertical ionization/attachment energies of the heme for the sampled configurations, and that the long range electrostatic interactions (up to 40 Å) contribute substantially to the absolute values of the vertical energy gaps.

Accurate vertical ionization energy and work function determinations of liquid water and aqueous solutions

The absolute-scale electronic energetics of liquid water and aqueous solutions, both in the bulk and at associated interfaces, are the central determiners of water-based chemistry. However, such information is generally experimentally inaccessible. Here we demonstrate that a refined implementation of the liquid microjet photoelectron spectroscopy (PES) technique can be adopted to address this. Implementing concepts from condensed matter physics, we establish novel all-liquid-phase vacuum and equilibrated solution–metal-electrode Fermi level referencing procedures. This enables the precise and accurate determination of previously elusive water solvent and solute vertical ionization energies, VIEs. Notably, this includes quantification of solute-induced perturbations of water's electronic energetics and VIE definition on an absolute and universal chemical potential scale. Defining and applying these procedures over a broad range of ionization energies, we accurately and respectively determine the VIE and oxidative stability of liquid water as 11.33 ± 0.03 eV and 6.60 ± 0.08 eV with respect to its liquid-vacuum-interface potential and Fermi level. Combining our referencing schemes, we accurately determine the work function of liquid water as 4.73 ± 0.09 eV. Further, applying our novel approach to a pair of exemplary aqueous solutions, we extract absolute VIEs of aqueous iodide anions, reaffirm the robustness of liquid water's electronic structure to high bulk salt concentrations (2 M sodium iodide), and quantify reference-level dependent reductions of water's VIE and a 0.48 ± 0.13 eV contraction of the solution's work function upon partial hydration of a known surfactant (25 mM tetrabutylammonium iodide). Our combined experimental accomplishments mark a major advance in our ability to quantify electronic–structure interactions and chemical reactivity in liquid water, which now explicitly extends to the measurement of absolute-scale bulk and interfacial solution energetics, including those of relevance to aqueous electrochemical processes.

A generalised liquid-phase photoelectron spectroscopy approach is reported, allowing accurate, absolute energy scale ionisation energies of liquid water and aqueous solutions, as well as liquid water's work function to be reported.

Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory

Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.

Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface

Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.

The electronic structure of the aqueous permanganate ion: aqueous-phase energetics and molecular bonding studied using liquid jet photoelectron spectroscopy

Permanganate aqueous solutions, MnO4-(aq.), were studied using liquid-micro-jet-based soft X-ray non-resonant and resonant photoelectron spectroscopy to determine valence and core-level binding energies. To identify possible differences in the energetics between the aqueous bulk and the solution-gas interface, non-resonant spectra were recorded at two different probing depths. Similar experiments were performed with different counter ions, Na+ and K+, with the two solutions yielding indistinguishable anion electron binding energies. Our resonant photoelectron spectroscopy measurements, performed near the Mn LII,III- and O K-edges, selectively probed valence charge distributions between the Mn metal center, O ligands, and first solvation shell in the aqueous bulk. Associated resonantly-enhanced solute ionisation signals revealed hybridisation of the solute constituents' atomic orbitals, including the inner valence Mn 3p and O 2s. We identified intermolecular coulombic decay relaxation processes following resonant X-ray excitation of the solute that highlight valence MnO4-(aq.)-H2O(l) electronic couplings. Furthermore, our results allowed us to infer oxidative reorganisation energies of MnO4˙(aq.) and adiabatic valence ionisation energies of MnO4-(aq.), revealing the Gibbs free energy of oxidation and permitting estimation of the vertical electron affinity of MnO4˙(aq.). Finally, the Gibbs free energy of hydration of isolated MnO4- was determined. Our results and analysis allowed a near-complete binding-energy-scaled MnO4-(aq.) molecular orbital and a valence energy level diagram to be produced for the MnO4-(aq.)/MnO4˙(aq.) system. Cumulatively, our mapping of the aqueous-phase electronic structure of MnO4- is expected to contribute to a deeper understanding of the exceptional redox properties of this widely applied aqueous transition-metal complex ion.

An experimental and computational study of the effect of aqueous solution on the multiphoton ionisation photoelectron spectrum of phenol

We revisit the photoelectron spectroscopy of aqueous phenol in an effort to improve our understanding of the impact of inhomogeneous broadening and inelastic scattering on solution-phase photoelectron spectra.

Hydration interactions beyond the first solvation shell in aqueous phenolate solution

We investigate the orientational dynamics of water molecules solvating phenolate ions using ultrafast vibrational spectroscopy and density functional theory-based molecular dynamics simulations.

Design and characterization of a recirculating liquid-microjet photoelectron spectrometer for multiphoton ultraviolet photoelectron spectroscopy

A new recirculating liquid-microjet photoelectron spectrometer for multiphoton ultraviolet photoelectron spectroscopy is described. A recirculating system is essential for studying samples that are only available in relatively small quantities. The reduction in background pressure when using the recirculating system compared to a liquid-nitrogen cold-trap results in a significant improvement in the quality of the photoelectron spectra. Moreover, the recirculating system results in a negligible streaming potential. The instrument design, operation, and characterization are described in detail, and its performance is illustrated by comparing a photoelectron spectrum of aqueous phenol recorded using the recirculating system with one recorded using a liquid nitrogen cold-trap.

Role of solvent-anion charge transfer in oxidative degradation of battery electrolytes

Electrochemical stability windows of electrolytes largely determine the limitations of operating regimes of lithium-ion batteries, but the degradation mechanisms are difficult to characterize and poorly understood. Using computational quantum chemistry to investigate the oxidative decomposition that govern voltage stability of multi-component organic electrolytes, we find that electrolyte decomposition is a process involving the solvent and the salt anion and requires explicit treatment of their coupling. We find that the ionization potential of the solvent-anion system is often lower than that of the isolated solvent or the anion. This mutual weakening effect is explained by the formation of the anion-solvent charge-transfer complex, which we study for 16 anion-solvent combinations. This understanding of the oxidation mechanism allows the formulation of a simple predictive model that explains experimentally observed trends in the onset voltages of degradation of electrolytes near the cathode. This model opens opportunities for rapid rational design of stable electrolytes for high-energy batteries.

Ultrafast structural rearrangement dynamics induced by the photodetachment of phenoxide in aqueous solution

The elementary processes that accompany the interaction of ionizing radiation with biologically relevant molecules are of fundamental importance. However, the ultrafast structural rearrangement dynamics induced by the ionization of biomolecules in aqueous solution remain hitherto unknown. Here, we employ femtosecond optical pump-probe spectroscopy to elucidate the vibrational wave packet dynamics that follow the photodetachment of phenoxide, a structural mimic of tyrosine, in aqueous solution. Photodetachment of phenoxide leads to wave packet dynamics of the phenoxyl radical along 12 different vibrational modes. Eight of the modes are totally symmetric and support structural rearrangement upon electron ejection. Comparison to a previous photodetachment study of phenoxide in the gas phase reveals the important role played by the solvent environment in driving ultrafast structural reorganization induced by ionizing radiation. This work provides insight into the ultrafast molecular dynamics that follow the interaction of ionizing radiation with molecules in aqueous solution.

The interaction of biomolecules with ionizing radiation induces structural changes which are still largely unknown. The authors use femtosecond wave packet spectroscopy to observe ultrafast structural dynamics that follow the photodetachment of phenoxide in aqueous solution.

Discrete Solvent Reaction Field Calculations for One- and Two-Photon Absorptions of Solution-Phase Dimethylaminonitrostilbene Molecule

Based on the configurations generated by molecular dynamics (MD) simulations using the on-the-fly density-functional tight-bonding (DFTB) force field, we investigated performance of the discrete solvent reaction field (DRF) model coupled to time-dependent density functional theory (TD-DFT) for solvatochromic effect of one- and two-photon absorption phenomena. Dimethylaminonitrostilbene (DANS) molecule solvated in chloroform, dichloromethane, and dimethyl sulfoxide solvents was selected as a model system for our research purpose. For every selected MD/DFTB configuration, within the context of the DRF, solute molecule is represented by TD-DFT and solvent molecules are described by atomic charges and polarizabilities. The calculated one-photon absorption energies reproduce well the positive solvatochromic behavior of solvated DANS and are in good agreement with available experimental data. For the two-photon absorption cross section, even though our approach overshot the experimental data by about 20% in absolute magnitude, experimentally observed solvatochromic change was captured qualitatively in this work. At last, we examined the contributions of atomic charges and polarizabilities of solvent molecules to the solvatochromic shifts of properties of interest.

Polarizable embedding for simulating redox potentials of biomolecules

Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana. The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (-0.15 V).

To Be or Not To Be a Molecular Ion: The Role of the Solvent in Photoionization of Arginine

Application of photoionization mass spectroscopy, a technique capable of assessing protonation states in complex molecules in the gas phase, is challenging for arginine due to its fragility. We report photoionization efficiencies in the valence region of aqueous aerosol particles produced from arginine solutions under various pH and vaporization conditions. By using ab initio calculations, we investigate the stability of different conformers. Our results show that neutral arginine fragments upon ionization in the gas phase but solvation stabilizes the molecular ion, resulting in different photoionization dynamics. We also report the valence-band photoelectron spectra of the aerosol solutions obtained at different pH values.