Eva Pluhařová

Chasing charge localization and chemical reactivity following photoionization in liquid water

The ultrafast dynamics of the cationic hole formed in bulk liquid water following ionization is investigated by ab initio molecular dynamics simulations and an experimentally accessible signature is suggested that might be tracked by femtosecond pump-probe spectroscopy. This is one of the fastest fundamental processes occurring in radiation-induced chemistry in aqueous systems and biological tissue. However, unlike the excess electron formed in the same process, the nature and time evolution of the cationic hole has been hitherto little studied. Simulations show that an initially partially delocalized cationic hole localizes within ~30 fs after which proton transfer to a neighboring water molecule proceeds practically immediately, leading to the formation of the OH radical and the hydronium cation in a reaction which can be formally written as H(2)O(+) + H(2)O → OH + H(3)O(+). The exact amount of initial spin delocalization is, however, somewhat method dependent, being realistically described by approximate density functional theory methods corrected for the self-interaction error. Localization, and then the evolving separation of spin and charge, changes the electronic structure of the radical center. This is manifested in the spectrum of electronic excitations which is calculated for the ensemble of ab initio molecular dynamics trajectories using a quantum mechanics/molecular mechanics (QM∕MM) formalism applying the equation of motion coupled-clusters method to the radical core. A clear spectroscopic signature is predicted by the theoretical model: as the hole transforms into a hydroxyl radical, a transient electronic absorption in the visible shifts to the blue, growing toward the near ultraviolet. Experimental evidence for this primary radiation-induced process is sought using femtosecond photoionization of liquid water excited with two photons at 11 eV. Transient absorption measurements carried out with ~40 fs time resolution and broadband spectral probing across the near-UV and visible are presented and direct comparisons with the theoretical simulations are made. Within the sensitivity and time resolution of the current measurement, a matching spectral signature is not detected. This result is used to place an upper limit on the absorption strength and/or lifetime of the localized H(2)O(+) ((aq)) species.

Ion Pairing in Aqueous Lithium Salt Solutions with Monovalent and Divalent Counter-Anions

Molecular dynamics simulations of concentrated aqueous solutions of LiCl and Li2SO4 were conducted in order to provide molecular insight into recent neutron scattering data. The structures predicted from the molecular dynamics simulations using standard nonpolarizable force fields provided a very poor fit to the experiment; therefore, refinement was needed. The electronic polarizability of the medium was effectively accounted for by implementing the electronic continuum correction, which practically means rescaling the ionic charges. Consistent with previous studies, we found that this approach in each case provided a significantly improved fit to the experimental data, which was further enhanced by slightly adjusting the radius of the lithium ion. The polarization effect was particularly pronounced in the Li2SO4 solution where the ions in the nonpolarizable simulations tended to cluster unphysically. With the above alterations, the employed force field displayed an excellent fit to the neutron scattering data and provided a useful interpretative framework for the experimental measurements. At the same time, the present study underlines the importance of solvent polarization effects in hydration of ions with high charge density.

Ionization of purine tautomers in nucleobases, nucleosides, and nucleotides: from the gas phase to the aqueous environment.

We have simulated ionization of purine nucleic acid components in the gas phase and in a water environment. The vertical and adiabatic ionization processes were calculated at the PMP2/aug-cc-pVDZ level with the TDDFT method applied to obtain ionization from the deeper lying orbitals. The water environment was modeled via microsolvation approach and using a nonequilibrium polarizable continuum model. We have characterized a set of guanine tautomers and investigated nucleosides and nucleotides in different conformations. The results for guanine, i.e., the nucleic acid base with the lowest vertical ionization potential, were also compared to those for the other purine base, adenine. The main findings of our study are the following: (i) Guanine remains clearly the base with the lowest ionization energy even upon aqueous solvation. (ii) Water solvent has a strong effect on the ionization energetics of guanine and adenine and their derivatives; the vertical ionization potential (VIP) is lowered by about 1 eV for guanine while it is ∼1.5 eV higher in the nucleotides, overall resulting in similar VIPs for GMP(-), guanosine and guanine in water. (iii) Water efficiently screens the electrostatic interactions between nucleic acid components. Consequently, ionization in water always originates from the base unit of the nucleic acid and all the information about conformational state is lost in the ionization energetics. (iv) The energy splitting between ionization of the two least bound electrons increases upon solvation. (v) Tautomerism does not contribute to the width of the photoelectron spectra in water. (vi) The effect of specific short-range interactions with individual solvent molecules is negligible for purine bases, compared to the long-range dielectric effects of the aqueous medium.

Transforming Anion Instability into Stability: Contrasting Photoionization of Three Protonation Forms of the Phosphate Ion upon Moving into Water

We use photoelectron emission spectroscopy with vacuum microjet technique and quantum chemistry calculations to investigate electronic structure and stability of aqueous phosphate anions. On the basis of the measured photoelectron spectra of sodium phosphates at different pH, we report the lowest vertical ionization energies of monobasic (9.5 eV), dibasic (8.9 eV), and tribasic (8.4 eV) anions. Electron binding energies were in tandem modeled with ab initio methods, using a mixed dielectric solvation model together with up to 64 explicitly solvating water molecules. We demonstrate that two solvation layers of explicit water molecules are needed to obtain converged values of vertical ionization energies (VIEs) within this mixed solvation model, leading to very good agreement with experiment. We also show that the highly charged PO(4)(3-) anion, which is electronically unstable in the gas phase, gains the electronic stability with about 16 water molecules, while only 2-3 water molecules are sufficient to stabilize the doubly charged phosphate anion. We also investigate the effect of ion pairing on the vertical ionization energy. In contrast to protonation (leading to a formation of covalent O-H bond), sodiation (leading to an anion···Na(+) ion pair) has only a weak effect on the electron binding energy.

Oxidation Half-Reaction of Aqueous Nucleosides and Nucleotides via Photoelectron Spectroscopy Augmented by ab Initio Calculations

Oxidative damage to DNA and hole transport between nucleobases in oxidized DNA are important processes in lesion formation for which surprisingly poor thermodynamic data exist, the relative ease of oxidizing the four nucleobases being one such example. Theoretical simulations of radiation damage and charge transport in DNA depend on accurate values for vertical ionization energies (VIEs), reorganization energies, and standard reduction potentials. Liquid-jet photoelectron spectroscopy can be used to directly study the oxidation half-reaction. The VIEs of nucleic acid building blocks are measured in their native buffered aqueous environment. The experimental investigation of purine and pyrimidine nucleotides, nucleosides, pentose sugars, and inorganic phosphate demonstrates that photoelectron spectra of nucleotides arise as a spectral sum over their individual chemical components; that is, the electronic interactions between each component are effectively screened from one another by water. Electronic structure theory affords the assignment of the lowest energy photoelectron band in all investigated nucleosides and nucleotides to a single ionizing transition centered solely on the nucleobase. Thus, combining the measured VIEs with theoretically determined reorganization energies allows for the spectroscopic determination of the one-electron redox potentials that have been difficult to establish via electrochemistry.

Exploring Ion–Ion Interactions in Aqueous Solutions by a Combination of Molecular Dynamics and Neutron Scattering

Recent advances in computational and experimental techniques have allowed for accurate description of ion pairing in aqueous solutions. Free energy methods based on ab initio molecular dynamics, as well as on force fields accounting effectively for electronic polarization, can provide quantitative information about the structures and occurrences of individual types of ion pairs. When properly benchmarked against electronic structure calculations for model systems and against structural experiments, in particular neutron scattering, such force field simulations represent a powerful tool for elucidating interactions of salt ions in complex biological aqueous environments.

Modeling photoionization of aqueous DNA and its components.

Radiation damage to DNA is usually considered in terms of UVA and UVB radiation. These ultraviolet rays, which are part of the solar spectrum, can indeed cause chemical lesions in DNA, triggered by photoexcitation particularly in the UVB range. Damage can, however, be also caused by higher energy radiation, which can ionize directly the DNA or its immediate surroundings, leading to indirect damage. Thanks to absorption in the atmosphere, the intensity of such ionizing radiation is negligible in the solar spectrum at the surface of Earth. Nevertheless, such an ionizing scenario can become dangerously plausible for astronauts or flight personnel, as well as for persons present at nuclear power plant accidents. On the beneficial side, ionizing radiation is employed as means for destroying the DNA of cancer cells during radiation therapy. Quantitative information about ionization of DNA and its components is important not only for DNA radiation damage, but also for understanding redox properties of DNA in redox sensing or labeling, as well as charge migration along the double helix in nanoelectronics applications. Until recently, the vast majority of experimental and computational data on DNA ionization was pertinent to its components in the gas phase, which is far from its native aqueous environment. The situation has, however, changed for the better due to the advent of photoelectron spectroscopy in liquid microjets and its most recent application to photoionization of aqueous nucleosides, nucleotides, and larger DNA fragments. Here, we present a consistent and efficient computational methodology, which allows to accurately evaluate ionization energies and model photoelectron spectra of aqueous DNA and its individual components. After careful benchmarking, the method based on density functional theory and its time-dependent variant with properly chosen hybrid functionals and polarizable continuum solvent model provides ionization energies with accuracy of 0.2-0.3 eV, allowing for faithful modeling and interpretation of DNA photoionization. The key finding is that the aqueous medium is remarkably efficient in screening the interactions within DNA such that, unlike in the gas phase, ionization of a base, nucleoside, or nucleotide depends only very weakly on the particular DNA context. An exception is the electronic interaction between neighboring bases which can lead to sequence-specific effects, such as a partial delocalization of the cationic hole upon ionization enabled by presence of adjacent bases of the same type.

Aqueous Cation-Amide Binding: Free Energies and IR Spectral Signatures by Ab Initio Molecular Dynamics.

Understanding specific ion effects on proteins remains a considerable challenge. N-methylacetamide serves as a useful proxy for the protein backbone that can be well characterized both experimentally and theoretically. The spectroscopic signatures in the amide I band reflecting the strength of the interaction of alkali cations and alkaline earth dications with the carbonyl group remain difficult to assign and controversial to interpret. Herein, we directly compute the infrared (IR) shifts corresponding to the binding of either sodium or calcium to aqueous N-methylacetamide using ab initio molecular dynamics simulations. We show that the two cations interact with aqueous N-methylacetamide with different affinities and in different geometries. Because sodium exhibits a weak interaction with the carbonyl group, the resulting amide I band is similar to an unperturbed carbonyl group undergoing aqueous solvation. In contrast, the stronger calcium binding results in a clear IR shift with respect to N-methylacetamide in pure water.

Hydration of the chloride ion in concentrated aqueous solutions using neutron scattering and molecular dynamics

Neutron scattering experiments were performed on 6 m LiCl solutions in order to obtain the solvation structure around the chloride ion. Molecular dynamics simulations on systems mirroring the concentrated electrolyte conditions of the experiment were carried out with a variety of chloride force-fields. In each case the simulations were run with both full ionic charges and employing the electronic continuum correction (implemented through charge scaling) to account effectively for electronic polarisation. The experimental data were then used to assess the successes and shortcomings of the investigated force-fields. We found that due to the very good signal-to-noise ratio in the experimental data, they provide a very narrow window for the position of the first hydration shell of the chloride ion. This allowed us to establish the importance of effectively accounting for electronic polarisation, as well as adjusting the ionic size, for obtaining a force-field which compares quantitatively to the experimental data. The present results emphasise the utility of performing neutron diffraction with isotopic substitution as a powerful tool in gaining insight and examining the validity of force-fields in concentrated electrolyte solutions.

Peptide salt bridge stability:From gas phase via microhydration to bulk water simulations

The salt bridge formation and stability in the terminated lysine-glutamate dipeptide is investigated in water clusters of increasing size up to the limit of bulk water. Proton transfer dynamics between the acidic and basic side chains is described by DFT-based Born-Oppenheimer molecular dynamics simulations. While the desolvated peptide prefers to be in its neutral state, already the addition of a single water molecule can trigger proton transfer from the glutamate side chain to the lysine side chain, leading to a zwitterionic salt bridge state. Upon adding more water molecules we find that stabilization of the zwitterionic state critically depends on the number of hydrogen bonds between side chain termini, the water molecules, and the peptidic backbone. Employing classical molecular dynamics simulations for larger clusters, we observed that the salt bridge is weakened upon additional hydration. Consequently, long-lived solvent shared ion pairs are observed for about 30 water molecules while solvent separated ion pairs are found when at least 40 or more water molecules hydrate the dipeptide. These results have implications for the formation and stability of salt bridges at partially dehydrated surfaces of aqueous proteins.