Modeling Liquid Photoemission Spectra: Path-Integral Molecular Dynamics Combined with Tuned Range-Separated Hybrid Functionals

Deconvolution of the X-ray absorption spectrum of trans-1,3-butadiene with resonant Auger spectroscopy

High-resolution carbon K-edge X-ray photoelectron, X-ray absorption, non-resonant and resonant Auger spectra are presented of gas phase trans-1,3-butadiene alongside a detailed theoretical analysis utilising nuclear ensemble approaches and vibronic models to simulate the spectroscopic observables. The resonant Auger spectra recorded across the first pre-edge band reveal a complex evolution of different electronic states which remain relatively well-localised on the edge or central carbon sites. The results demonstrate the sensitivity of the resonant Auger observables to the weighted contributions from multiple electronic states. The gradually evolving spectral features can be accurately and feasibly simulated within nuclear ensemble methods and interpreted with the population analysis.

Jahn-Teller effects in initial and final states: high-resolution X-ray absorption, photoelectron and Auger spectroscopy of allene

Carbon K-edge resonant Auger spectra of gas-phase allene following excitation of the pre-edge 1s → π* transitions are presented and analysed with the support of EOM-CCSD/cc-pVTZ calculations. X-Ray absorption (XAS), X-ray photoelectron (XPS), valence band and non-resonant Auger spectra are also reanalysed with a series of computational approaches. The results presented demonstrate the importance of including nuclear ensemble effects for simulating X-ray observables and as an effective strategy for capturing Jahn-Teller effects in spectra.

Ultrafast Charge and Proton Transfer in Doubly Ionized Ammonia Dimers

We investigate the ultrafast energy and charge transfer processes between ammonia molecules following ionization reactions initiated by electron impact. Exploring ionization-induced processes in molecular clusters provides us with a detailed insight into the dynamics using experiments in the energy domain. We ionize the ammonia dimer with 200 eV electrons and apply the fragment ions coincident momentum spectroscopy and nonadiabatic molecular dynamics simulations. We identify two mechanisms leading to the doubly charged ammonia dimer. In the first one, a single molecule is ionized. This initiates an ultrafast proton transfer process, leading to the formation of the NH₂⁺ + NH₄⁺ pair. Alternatively, a dimer with a delocalized charge is formed dominantly via the intermolecular Coulombic decay, forming the NH₃⁺·NH₃⁺ dication. This dication further dissociates into two NH₃⁺ cations. The ab initio calculations have reproduced the measured kinetic energy release of the ion pairs and revealed the dynamical processes following the double ionization.

Different timescales during ultrafast stilbene isomerization in the gas and liquid phases revealed using time-resolved photoelectron spectroscopy

Directly contrasting ultrafast excited-state dynamics in the gas and liquid phases is crucial to understanding the influence of complex environments. Previous studies have often relied on different spectroscopic observables, rendering direct comparisons challenging. Here, we apply extreme-ultraviolet time-resolved photoelectron spectroscopy to both gaseous and liquid cis-stilbene, revealing the coupled electronic and nuclear dynamics that underlie its isomerization. Our measurements track the excited-state wave packets from excitation along the complete reaction path to the final products. We observe coherent excited-state vibrational dynamics in both phases of matter that persist to the final products, enabling the characterization of the branching space of the S₁-S₀ conical intersection. We observe a systematic lengthening of the relaxation timescales in the liquid phase and a red shift of the measured excited-state frequencies that is most pronounced for the complex reaction coordinate. These results characterize in detail the influence of the liquid environment on both electronic and structural dynamics during a complete photochemical transformation.

Combining Localized Orbital Scaling Correction and Bethe-Salpeter Equation for Accurate Excitation Energies

We applied localized orbital scaling correction (LOSC) in Bethe-Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar or better than the GW Green's function approach and with much less computational cost. The QP energies from LOSC instead of commonly used G0 W0 and ev GW are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/ G0W0 approach for predicting excitations with different characters. For the calculations for Truhlar-Gagliardi test set containing valence, charge transfer (CT) and Rydberg excitations, BSE/LOSC with the Tamm-Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/ev GW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/ G0W0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is thus a promising and efficient approach to calculate excitation energies for molecular systems.

Calculating Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds

Characterizing the photochemical reactivity of transient volatile organic compounds (VOCs) in our atmosphere begins with a proper understanding of their abilities to absorb sunlight. Unfortunately, the photoabsorption cross-sections for a large number of transient VOCs remain unavailable experimentally due to their short lifetime or high reactivity. While structure-activity relationships (SARs) have been successfully employed to estimate the unknown photoabsorption cross-sections of VOCs, computational photochemistry offers another promising strategy to predict not only the vertical electronic transitions of a given molecule but also the width and shape of the bands forming its absorption spectrum. In this work, we focus on the use of the nuclear ensemble approach (NEA) to determine the photoabsorption cross-section of four exemplary VOCs, namely, acrolein, methylhydroperoxide, 2-hydroperoxy-propanal, and (microsolvated) pyruvic acid. More specifically, we analyze the influence that different strategies for sampling the ground-state nuclear density-Wigner sampling and ab initio molecular dynamics with a quantum thermostat-can have on the simulated absorption spectra. We highlight the potential shortcomings of using uncoupled harmonic modes within Wigner sampling of nuclear density to describe flexible or microsolvated VOCs and some limitations of SARs for multichromophoric VOCs. Our results suggest that the NEA could constitute a powerful tool for the atmospheric community to predict the photoabsorption cross-section for transient VOCs.

On-the-Fly Ring-Polymer Molecular Dynamics Calculations of the Dissociative Photodetachment Process of the Oxalate Anion

The dissociative photodetachment dynamics of the oxalate anion, C₂O₄H^- + hν → CO₂ + HOCO + e^-, were theoretically studied using the on-the-fly path-integral and ring-polymer molecular dynamics methods, which can account for nuclear quantum effects at the density-functional theory level in order to compare with the recent experimental study using photoelectron-photofragment coincidence spectroscopy. To reduce computational time, the force acting on each bead of ring-polymer was approximately calculated from the first and second derivatives of the potential energy at the centroid position of the nuclei beads. We find that the calculated photoelectron spectrum qualitatively reproduces the experimental spectrum and that nuclear quantum effects are playing a role in determining spectral widths. The calculated coincidence spectrum is found to reasonably reproduce the experimental spectrum, indicating that a relatively large energy is partitioned into the relative kinetic energy between the CO₂ and HOCO fragments. This is because photodetachment of the parent anion leads to Franck-Condon transition to the repulsive region of the neutral potential energy surface. We also find that the dissociation dynamics are slightly different between the two isomers of the C₂O₄H^- anion with closed- and open-form structures.

Determination of the valence band edge of Fe oxide nanoparticles dispersed in aqueous solution through resonant photoelectron spectroscopy from a liquid microjet

We use X-ray photoemission and a near ambient pressure with a liquid microjet setup to investigate the electronic structure of FeOOH nanoparticles dispersed in aqueous solution. In particular, we show that by using X-ray resonant photoemission in dilute solutions, we can overcome the limits of conventional photoemission such as low nanoparticle-to-solvent signal ratio, and local nanoparticle charging and measure the valence band structure of FeOOH nanoparticles in aqueous solution with chemical specificity. The resonant photoemission signal across the Fe 2p_3/2 absorption edge is measured for 2 wt% aqueous solutions of FeOOH nanoparticles (NPs) and the valence band maximum (VBM) of the hydrated FeOOH nanoparticles is determined. We compare the obtained VBM value in aqueous solution to that of FeOOH NPs in the dry phase. We show that the valence band edge position of NPs in the liquid phase can be accurately predicted from the values obtained in the dry phase provided that a simple potential shift due to solution chemistry is applied. Our results demonstrate the suitability of resonant photoemission in measuring the electronic structure of strongly diluted nanosystems where the conventional non-resonant photoemission technique fails.

Electronic Structure of Water from Koopmans-Compliant Functionals

Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW + f_xc). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that different MD approaches lead to up to 1 eV change in the average band gap; thus, we focused on the band gap dependence with the geometrical properties of a system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance and HOH angle, as well as the intermolecular hydrogen bond O···O distance and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O···O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.

Quantum calculations of the photoelectron spectra of the OH-·NH3 anion: implications for OH + NH3→ H2O + NH2 reaction dynamics

We present the results of quantum dynamics calculations for analyzing the experimentally measured photoelectron spectra of the OH^-·NH₃ anion complex. Detachment of an excess electron of OH^-·NH₃ initially produces a molecular arrangement, which is close to the transition-state structure of the neutral OH + NH₃→ H₂O + NH₂ hydrogen abstraction reaction due to the Franck-Condon principle, and thus finally leads to the OH + NH₃ or H₂O + NH₂ asymptotic channel. We used both the path integral method and the reduced-dimensionality quantum wave packet method to simulate the photoelectron spectra of the OH^-·NH₃ anion. The calculated spectra were found to be in qualitative agreement with the experimental spectra. It was found that the photodetached complex mainly dissociates into the OH + NH₃ channel; however, we found that the hydrogen exchange process also contributes to the photodetachment spectra.

Vibrationally Mediated Stabilization of Electrons in Nonpolar Matter

We explore solvation of electrons in nonpolar matter, here represented by butadiene clusters. Isolated butadiene supports only the existence of transient anions (resonances). Two-dimensional electron energy loss spectroscopy shows that the resonances lead to an efficient vibrational excitation of butadiene, which can result into the almost complete loss of energy of the interacting electron. Cluster-beam experiments show that molecular clusters of butadiene form stable anions, however only at sizes of more than 9 molecular units. We have calculated the distribution of electron affinities of clusters using classical and path integral molecular dynamics simulations. There is almost a continuous transition from the resonant to the bound anions with an increase in cluster size. The comparison of the classical and quantum dynamics reveals that the electron binding is strongly supported by molecular vibrations, brought about by nuclear zero-point motion and thermal agitation. We also inspected the structure of the solvated electron, finding it well localized.

Photoemission from non-polar aromatic molecules in the gas and liquid phase

The photoelectron spectra of both liquid and gas phase aromatic molecules are reported. The spectra were obtained using a 34.1 eV source produced by high harmonic generation and analysed with the help of high-level ab initio simulations using the reflection principle combined with path integral molecular dynamics simulations accounting for nuclear quantum effects for the gas phase. We demonstrate the suitability of three trimethylbenzenes (1,3,5-trimethylbenzene, 1,2,3-trimethylbenzene and 1,2,4-trimethylbenzene) as a solvent for liquid photoelectron spectroscopy of solute species. We also discuss the electrokinetic charging of a non-polar liquid jet.

Beyond Koopmans' theorem: electron binding energies in disordered materials

The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.

Do water's electrons care about electrolytes?

Ions have a profound effect on the geometrical structure of liquid water and an aqueous environment is known to change the electronic structure of ions. Here we combine photoelectron spectroscopy measurements from liquid microjets with molecular dynamical and quantum chemical calculations to address the reverse question, to what extent do ions affect the electronic structure of liquid water? We study aqueous solutions of sodium iodide (NaI) over a wide concentration range, from nearly pure water to 8 M solutions, recording spectra in the 5 to 60 eV binding energy range to include all water valence and the solute Na+ 2p, I- 4d, and I- 5p orbital ionization peaks. We observe that the electron binding energies of the solute ions change only slightly as a function of electrolyte concentration, less than 150 ± 60 meV over an ∼8 M range. Furthermore, the photoelectron spectrum of liquid water is surprisingly mildly affected as we transform the sample from a dilute aqueous salt solution to a viscous, crystalline-like phase. The most noticeable spectral changes are a negative binding energy shift of the water 1b2 ionizing transition (up to -370 ± 60 meV) and a narrowing of the flat-top shape water 3a1 ionization feature (up to 450 ± 90 meV). A novel computationally efficient technique is introduced to calculate liquid-state photoemission spectra using small clusters from molecular dynamics (MD) simulations embedded in dielectric continuum. This theoretical treatment captured the characteristic positions and structures of the aqueous photoemission peaks, reproducing the experimentally observed narrowing of the water 3a1 feature and weak sensitivity of the water binding energies to electrolyte concentration. The calculations allowed us to attribute the small binding energy shifts to ion-induced disruptions of intermolecular electronic interactions. Furthermore, they demonstrate the importance of considering concentration-dependent screening lengths for a correct description of the electronic structure of solvated systems. Accounting for electronic screening, the calculations highlight the minimal effect of electrolyte concentration on the 1b1 binding energy reference, in accord with the experiments. This leads us to a key finding that the isolated, lowest-binding-energy, 1b1, photoemission feature of liquid water is a robust energetic reference for aqueous liquid microjet photoemission studies.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water

Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.

The importance of nuclear quantum effects in spectral line broadening of optical spectra and electrostatic properties in aromatic chromophores

In this work, we examine the importance of nuclear quantum effects on capturing the line broadening and vibronic structure of optical spectra. We determine the absorption spectra of three aromatic molecules indole, pyridine, and benzene using time dependent density functional theory with several molecular dynamics sampling protocols: force-field based empirical potentials, ab initio simulations, and finally path-integrals for the inclusion of nuclear quantum effects. We show that the absorption spectrum for all these chromophores are similarly broadened in the presence of nuclear quantum effects regardless of the presence of hydrogen bond donor or acceptor groups. We also show that simulations incorporating nuclear quantum effects are able to reproduce the heterogeneous broadening of the absorption spectra even with empirical force fields. The spectral broadening associated with nuclear quantum effects can be accounted for by the broadened distribution of chromophore size as revealed by a particle in the box model. We also highlight the role that nuclear quantum effects have on the underlying electronic structure of aromatic molecules as probed by various electrostatic properties.

Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer

Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.

On the importance of initial conditions for excited-state dynamics

The vast majority of ab initio excited-state simulations are performed within semiclassical, trajectory-based approaches. Apart from the underlying electronic-structure theory, the reliability of the simulations is controlled by a selection of initial conditions for the classical trajectories. We discuss appropriate choices of initial conditions for simulations of different experimental arrangements: dynamics initiated by continuum-wave (CW) laser fields or triggered by ultrashort laser pulses.

Proton transfer in guanine-cytosine base pair analogues studied by NMR spectroscopy and PIMD simulations

It has been hypothesised that proton tunnelling between paired nucleobases significantly enhances the formation of rare tautomeric forms and hence leads to errors in DNA replication. Here, we study nuclear quantum effects (NQEs) using deuterium isotope-induced changes of nitrogen NMR chemical shifts in a model base pair consisting of two tautomers of isocytosine, which form hydrogen-bonded dimers in the same way as the guanine-cytosine base pair. Isotope effects in NMR are consequences of NQEs, because ro-vibrational averaging of different isotopologues gives rise to different magnetic shielding of the nuclei. The experimental deuterium-induced chemical shift changes are compared with those calculated by a combination of path integral molecular dynamics (PIMD) simulations with DFT calculations of nuclear shielding. These calculations can directly link the observable isotope-induced shifts with NQEs. A comparison of the deuterium-induced changes of 15N chemical shifts with those predicted by PIMD simulations shows that inter-base proton transfer reactions do not take place in this system. We demonstrate, however, that NMR isotope shifts provide a unique possibility to study NQEs and to evaluate the accuracy of the computational methods used for modelling quantum effects in molecules. Calculations based on the PBE functional from the general-gradient-approximation family provided significantly worse predictions of deuterium isotope shifts than those with the hybrid B3LYP functional.