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.

Oxygen Radicals Entrapped between MgO Nanocrystals: Formation, Spectroscopic Fingerprints, and Reactivity toward Water.

Compaction of dehydroxylated MgO nanocrystal powders produces adsorbed oxygen radicals with characteristic UV-vis spectroscopic fingerprints. Identical absorption bands arise upon UV excitation in an oxygen atmosphere but in the absence of uniaxial pressure. Photophysical calculations on MgO gas-phase clusters reveal that the observed optical transitions at 4.4 and 3.0 eV are consistent with adsorbed superoxide (O2·-) and ozonide (O3·-) species, respectively. The presence of these oxygen radicals is corroborated by electron paramagnetic resonance spectroscopy. Upon reaction with interfacial water, oxygen radicals convert into diamagnetic products with no absorptions in the UV-vis range. Since superoxide O2·- and ozonide anions O3·- play a key role in a variety of processes in heterogeneous catalysis, sensing, or as transient species in cold sintering, their UV-vis spectroscopic detection will enable in situ monitoring of transient oxygen radicals inside metal oxide powders.

Resonant Inner-Shell Photofragmentation of Adamantane (C10H16).

Adamantane, the smallest diamondoid molecule with a symmetrical cage, contains two distinct carbon sites, CH and CH2. The ionization/excitation of the molecule leads to the cage opening and strong structural reorganization. While theoretical predictions suggest that the carbon site CH primarily causes the cage opening, the role of the other CH2 site remains unclear. In this study, we used advanced experimental Auger electron-ion coincidence techniques and theoretical calculations to investigate the fragmentation dynamics of adamantane after resonant inner-shell photoexcitation. Our results demonstrate that some fragmentation channels exhibit site-sensitivity of the initial core-hole location, indicating that different carbon site excitations could lead to unique cage opening mechanisms.

Spin-Vibronic Control of Intersystem Crossing in Iodine-Substituted Heptamethine Cyanines.

Spin-orbit coupling between electronic states of different multiplicity can be strongly coupled to molecular vibrations, and this interaction is becoming recognized as an important mechanism for controlling the course of photochemical reactions. Here, we show that the involvement of spin-vibronic coupling is essential for understanding the photophysics and photochemistry of heptamethine cyanines (Cy7), bearing iodine as a heavy atom in the C3' position of the chain and/or a 3H-indolium core, as potential triplet sensitizers and singlet oxygen producers in methanol and aqueous solutions. The sensitization efficiency was found to be an order of magnitude higher for the chain-substituted than the 3H-indolium core-substituted derivatives. Our ab initio calculations demonstrate that while all optimal structures of Cy7 are characterized by negligible spin-orbit coupling (tenths of cm-1) with no dependence on the position of the substituent, molecular vibrations lead to its significant increase (tens of cm-1 for the chain-substituted cyanines), which allowed us to interpret the observed position dependence.

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.

Observation of intermolecular Coulombic decay and shake-up satellites in liquid ammonia.

We report the first nitrogen 1s Auger-Meitner electron spectrum from a liquid ammonia microjet at a temperature of ∼223 K (-50 °C) and compare it with the simultaneously measured spectrum for gas-phase ammonia. The spectra from both phases are interpreted with the assistance of high-level electronic structure and ab initio molecular dynamics calculations. In addition to the regular Auger-Meitner-electron features, we observe electron emission at kinetic energies of 374-388 eV, above the leading Auger-Meitner peak (3a1 2). Based on the electronic structure calculations, we assign this peak to a shake-up satellite in the gas phase, i.e., Auger-Meitner emission from an intermediate state with additional valence excitation present. The high-energy contribution is significantly enhanced in the liquid phase. We consider various mechanisms contributing to this feature. First, in analogy with other hydrogen-bonded liquids (noticeably water), the high-energy signal may be a signature for an ultrafast proton transfer taking place before the electronic decay (proton transfer mediated charge separation). The ab initio dynamical calculations show, however, that such a process is much slower than electronic decay and is, thus, very unlikely. Next, we consider a non-local version of the Auger-Meitner decay, the Intermolecular Coulombic Decay. The electronic structure calculations support an important contribution of this purely electronic mechanism. Finally, we discuss a non-local enhancement of the shake-up processes.

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.

Solvation energies of ions with ensemble cluster-continuum approach.

Solvation free energies can be advantageously estimated by cluster-continuum approaches. They proved useful especially for systems with high charge density. However, the clusters are assumed to be single minimum rigid species. It is an invalid condition for larger clusters and it complicates the assessment of convergence with the system size. We present a new variant of the cluster-continuum approach, "Ensemble Cluster-Continuum" scheme, where the single minima problem is circumvented by a thermodynamic cycle based on vertical quantities (ionization energies, electron affinities). Solvation free energies are calculated for a charged-neutralized system and solvation correction for the vertical quantities is estimated for an ensemble of structures from molecular dynamics simulation. We test the scheme on a set of various types of anions and cations, we study the convergence of the cluster-continuum model and assess various types of errors. The quantitative data depend on the applied continuum solvation model yet the convergence is analogous. We argue that the assessment of convergence provides a measure of the reliability of the calculated solvation energies.

Deciphering the Structure-Property Relations in Substituted Heptamethine Cyanines.

Heptamethine cyanines (Cy7) are fluorophores essential for modern bioimaging techniques and chemistry. Here, we systematically evaluated the photochemical and photophysical properties of a library of Cy7 derivatives containing diverse substituents in different positions of the heptamethine chain. A single substitution allows modulation of their absorption maxima in the range of 693-805 nm and photophysical properties, such as quantum yields of singlet-oxygen formation, decomposition, and fluorescence or affinity to singlet oxygen, within 2-3 orders of magnitude. The same substituent in different positions of the chain often exhibits distinctly contradictory effects, demonstrating that both the type and position of the substituent are pivotal for the design of Cy7-based applications. The combination of experimental results with quantum-chemical calculations provides insights into the structure-property relationship, the elucidation of which will accelerate the development of cyanines with properties tailored for specific applications, such as fluorescent probes and sensors, photouncaging, photodynamic therapy, or singlet-oxygen detection.

Ionization energies in solution with the QM:QM approach.

We discuss a fragment-based QM:QM scheme as a practical way to access the energetics of vertical electronic processes in the condensed phase. In the QM:QM scheme, we decompose the large molecular system into small fragments, which interact solely electrostatically. The energies of the fragments are calculated in a self-consistent field generated by the other fragments and the total energy of the system is calculated as a sum of the fragment energies. We show on two test cases (cytosine and a sodium cation) that the method allows one to accurately simulate the shift of vertical ionization energies (VIE) while going from the gas phase to the bulk. For both examples, the predicted solvent shifts and peak widths estimated at the DFT level agree well with the experimental observations. We argue that the QM:QM approach is more suitable than either an electrostatic embedding based QM/MM approach, a full quantum description at the DFT level with a generally used functional or a combination of both. We also discuss the potential scope of the applicability for other electronic processes such as Auger decay.

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.

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.

Optimal Tuning of Range-Separated Hybrids for Solvated Molecules with Time-Dependent Density Functional Theory.

The applicability range of density functional theory (DFT) can be improved with no additional parametrization by imposing some exact conditions. Enforcing equality between the orbital energy of the highest occupied Kohn-Sham orbital and ionization energy, determined from the total energy difference between neutral and ionized states (ΔKS), leads to the concept of optimally tuned range-separated hybrid functionals. Here, we present an alternative tuning scheme for range-separated hybrid functionals based on enforcing the equality between the ΔKS ionization energy and the ionization energy calculated by means of the time-dependent DFT using the concept of ionization as an excitation to the distant center (OT-IEDC scheme). The scheme can be naturally applied to solvated systems described either within the explicit solvation or dielectric continuum models. We test the performance of the scheme on a benchmark set of molecules. We further show that the scheme allows for reliably modeling liquid phase photoemission spectra.

Observation of electron-transfer-mediated decay in aqueous solution.

Photoionization is at the heart of X-ray photoelectron spectroscopy (XPS), which gives access to important information on a sample's local chemical environment. Local and non-local electronic decay after photoionization-in which the refilling of core holes results in electron emission from either the initially ionized species or a neighbour, respectively-have been well studied. However, electron-transfer-mediated decay (ETMD), which involves the refilling of a core hole by an electron from a neighbouring species, has not yet been observed in condensed phase. Here we report the experimental observation of ETMD in an aqueous LiCl solution by detecting characteristic secondary low-energy electrons using liquid-microjet soft XPS. Experimental results are interpreted using molecular dynamics and high-level ab initio calculations. We show that both solvent molecules and counterions participate in the ETMD processes, and different ion associations have distinctive spectral fingerprints. Furthermore, ETMD spectra are sensitive to coordination numbers, ion-solvent distances and solvent arrangement.

Photoswitching of Azobenzene-Based Reverse Micelles above and at Subzero Temperatures As Studied by NMR and Molecular Dynamics Simulations.

We designed and studied the structure, dynamics, and photochemistry of photoswitchable reverse micelles (RMs) composed of azobenzene-containing ammonium amphiphile 1 and water in chloroform at room and subzero temperatures by NMR spectroscopy and molecular dynamics simulations. The NMR and diffusion coefficient analyses showed that micelles containing either the E or Z configuration of 1 are stable at room temperature. Depending on the water-to-surfactant molar ratio, the size of the RMs remains unchanged or is slightly reduced because of the partial loss of water from the micellar cores upon extensive E → Z or Z → E photoisomerization of the azobenzene group in 1. Upon freezing at 253 or 233 K, E-1 RMs partially precipitate from the solution but are redissolved upon warming whereas Z-1 RMs remain fully dissolved at all temperatures. Light-induced isomerization of 1 at low temperatures does not lead to the disintegration of RMs remaining in the solution; however, its scope is influenced by a precipitation process. To obtain a deeper molecular view of RMs, their structure was characterized by MD simulations. It is shown that RMs allow for amphiphile isomerization without causing any immediate significant structural changes in the micelles.

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

We present a computational protocol for modeling valence photoemission spectra of liquids. We use water as an experimentally well-characterized model system, and we represent its liquid state by larger finite-sized droplets. The photoemission spectrum is evaluated for an ensemble of structures along molecular dynamics simulations. The nuclear quantum effects are accounted for by ab initio based path-integral molecular dynamics simulations that are greatly accelerated with the so-called colored noise thermostat (PI+GLE) method. The ionization energies for the valence electrons are evaluated as orbital energies of optimally tuned range-separated hybrid functionals (OT-RSH). This approach provides Koopmans-type ionization energies including relaxation energy. We show that the present protocol can quantitatively describe the valence photoemission spectrum of liquid water, i.e., the positions, shapes, and widths of the photoemission peaks. With the PI+GLE simulations, even the subtle isotope effects that have been recently observed experimentally can be modeled. The electronic properties of finite-sized droplets are shown to converge rapidly to those of liquids. We discuss the importance of proper tuning of the range-separation parameter in OT-RSH as well as possible sources of error in our simulations. The present approach seems to be a viable route to modeling photoemission spectra of liquids, especially in conjunction with efficient implementation of density functional methods on graphical processing units.

Nature of CTAB/Water/Chloroform Reverse Micelles at Above- and Subzero Temperatures Studied by NMR and Molecular Dynamics Simulations.

The nature and stability of cetyltrimethylammonium bromide (CTAB) reverse micelles in chloroform formed above the critical micellar concentration at above- and subzero temperatures were examined by NMR and molecular dynamics simulations. The experiments showed that the supercooled micellar water pool becomes unstable upon cooling to relatively high temperatures (253 K), and smaller micelles are formed. Upon freezing at lower temperatures (233 K), micelles become completely frozen and remain intact in the solution. With an average hydrodynamic radius of approximately 1.3 nm, we estimate that the water pool contains approximately 50 water molecules, which is well below the onset of ice crystal formation. To support the experimental results, molecular dynamics simulations were used to model the structure of CTAB/water/chloroform reverse micelles of different sizes. The MD simulations show that the reverse micelles contain a water pool with bromide anions residing on its surface and their shape is nonspherical, especially in the case of larger water pools. Upon fast freezing, the mobility of the water molecules is suppressed, and the pool becomes more spherical.

The ice-vapor interface and the melting point of ice I(h) for the polarizable POL3 water model.

We use molecular dynamics simulations to determine the melting point of ice I(h) for the polarizable POL3 water force field (Dang, L. X. J. Chem. Phys.1992, 97, 2659). Simulations are performed on a slab of ice I(h) with two free surfaces at several different temperatures. The analysis of the time evolution of the total energy in the course of the simulations at the set of temperatures yields the melting point of the POL3 model to be T(m) = 180 ± 10 K. Moreover, the results of the simulations show that the degree of hydrogen-bond disorder occurring in the bulk of POL3 ice is larger (at the corresponding degree of undercooling) than in ice modeled by nonpolarizable water models. These results demonstrate that the POL3 water force field is rather a poor model for studying ice and ice-liquid or ice-vapor interfaces. While a number of polarizable water models have been developed over the past years, little is known about their performance in simulations of supercooled water and ice. This study thus highlights the need for testing of the existing polarizable water models over a broad range of temperatures, pressures, and phases, and developing a new polarizable water force field, reliable over larger areas of the phase diagram.

Glycine in an electronically excited state: ab initio electronic structure and dynamical calculations.

The goal of this study is to explore the photochemical processes following optical excitation of the glycine molecule into its two low-lying excited states. We employed electronic structure methods at various levels to map the PES of the ground state and the two low-lying excited states of glycine. It follows from our calculations that the photochemistry of glycine can be regarded as a combination of photochemical behavior of amines and carboxylic acid. The first channel (connected to the presence of amino group) results in ultrafast decay, while the channels characteristic for the carboxylic group occur on a longer time scale. Dynamical calculations provided the branching ratio for these channels. We also addressed the question whether conformationally dependent photochemistry can be observed for glycine. While electronic structure calculations favor this possibility, the ab initio multiple spawning (AIMS) calculations showed only minor relevance of the reaction path resulting in conformationally dependent dynamics.

Theoretical study of photoacidity of HCN: the effect of complexation with water.

The character of the hydrogen bonding and the excited state proton transfer (ESPT) in the model system HCN...H(2)O is investigated. The PES of the two lowest excited states of the H(2)O...HCN complex was calculated using the CASPT2 method. The nonadiabatic coupling of the two states of the (pi-->pi*) and (pi-->sigma*) character is responsible for the excited state proton/hydrogen transfer. Compared to the ground state, the barrier for this process is significantly smaller. An increased number of water molecules in the complex with cyclic hydrogen-bonded network causes a large blue shift of the state of the (pi-->sigma*) character. The question of the dissociation of the complex in its excited state is also addressed.