Electron affinity of liquid water

Predicting electronic screening for fast Koopmans spectral functional calculations

Koopmans spectral functionals are a powerful extension of Kohn-Sham density-functional theory (DFT) that enables the prediction of spectral properties with state-of-the-art accuracy. The success of these functionals relies on capturing the effects of electronic screening through scalar, orbital-dependent parameters. These parameters have to be computed for every calculation, making Koopmans spectral functionals more expensive than their DFT counterparts. In this work, we present a machine-learning model that-with minimal training-can predict these screening parameters directly from orbital densities calculated at the DFT level. We show in two prototypical use cases that using the screening parameters predicted by this model, instead of those calculated from linear response, leads to orbital energies that differ by less than 20 meV on average. Since this approach dramatically reduces run times with minimal loss of accuracy, it will enable the application of Koopmans spectral functionals to classes of problems that previously would have been prohibitively expensive, such as the prediction of temperature-dependent spectral properties. More broadly, this work demonstrates that measuring violations of piecewise linearity (i.e., curvature in total energies with respect to occupancies) can be done efficiently by combining frozen-orbital approximations and machine learning.

Dissecting the hydrogen bond network of water: Charge transfer and nuclear quantum effects

The molecular structure of water is dynamic, with intermolecular hydrogen (H) bond interactions being modified by both electronic charge transfer and nuclear quantum effects (NQEs). Electronic charge transfer and NQEs potentially change under acidic or basic conditions, but such details have not been measured. In this work, we developed correlated vibrational spectroscopy, a symmetry-based method that separates interacting from noninteracting molecules in self- and cross-correlation spectra, giving access to previously inaccessible information. We found that hydroxide (OH-) donated ~8% more negative charge to the H bond network of water, and hydronium (H3O+) accepted ~4% less negative charge from the H bond network of water. Deuterium oxide (D2O) had ~9% more H bonds compared with water (H2O), and acidic solutions displayed more dominant NQEs than basic ones.

The Effect of Electric Fields on Oxidization Processes at the Air-Water Interface

At the air-water interface, many reactions are accelerated, sometimes by several orders of magnitude. This phenomenon has proved to be particularly important in water microdroplets, where the spontaneous oxidation of many species stable in bulk has been experimentally demonstrated. Different theories have been proposed to explain this finding, but it is currently believed that the role of interfacial electric fields is key. In this work, we have carried out a quantum chemistry study aimed at shedding some light on this question. We have studied two prototypical processes in which a hydroxide anion transfers its excess electron to either the water environment or a dioxygen molecule. To model the interface, we use a cluster of 21 water molecules immersed in an electric field, and we examine the energetics of the studied reactions as a function of field magnitude. Our results reveal that electric fields close to those estimated for the neat air-water interface (∼0.15 V ⋅ Å-1) have a moderate effect on the reaction energetics and that much stronger fields (>1 V ⋅ Å-1) are required to get spontaneous electron transfer. Therefore, the study suggests that additional factors such as an excess charge in microdroplets need to be considered for explaining the experimental observations.

Current Status of the MB-pol Data-Driven Many-Body Potential for Predictive Simulations of Water Across Different Phases

Developing a molecular-level understanding of the properties of water is central to numerous scientific and technological applications. However, accurately modeling water through computer simulations has been a significant challenge due to the complex nature of the hydrogen-bonding network that water molecules form under different thermodynamic conditions. This complexity has led to over five decades of research and many modeling attempts. The introduction of the MB-pol data-driven many-body potential energy function marked a significant advancement toward a universal molecular model capable of predicting the structural, thermodynamic, dynamical, and spectroscopic properties of water across all phases. By integrating physics-based and data-driven (i.e., machine-learned) components, which correctly capture the delicate balance among different many-body interactions, MB-pol achieves chemical and spectroscopic accuracy, enabling realistic molecular simulations of water, from gas-phase clusters to liquid water and ice. In this review, we present a comprehensive overview of the data-driven many-body formalism adopted by MB-pol, highlight the main results and predictions made from computer simulations with MB-pol to date, and discuss the prospects for future extensions to data-driven many-body potentials of generic and reactive molecular systems.

Chalcogen Bonds Enable Efficient Photoreduction of Sulfur-Containing Heterocycles

Chalcogen bonding interactions have attracted significant attention in a broad chemistry community, with a particular focus on their ability to stabilise the key transition states in various organic synthetic routes. In this work, we demonstrate that they can also be harnessed in selective photoredox reactions, which cannot be otherwise achieved with alternative approaches to photoreduction. We demonstrate this concept through the photoreduction of the sulfur-containing DNA nucleoside precursor thioanhydrouridine to 2'-deoxy-thiouridine, revealing the previously unrecognized role of bisulfide in this process. Based on quantum chemical simulations, we identify a stable chalcogen-bonding complex of the hydrosulfide anion and thionhydrouridine (HS-⋅⋅⋅S contacts), which enables directional photoinduced electron transfer, resulting in the formation of non-canonical DNA nucleoside. We also disprove the possibility that photoreduction of thioanhydronucleosides could be initiated by hydrated electrons generated from irradiated bisulfide anions which do not interact with the chromophore. Finally, we show that selective photoreduction mediated by chalcogen bonds can only occur for chromophores, which exhibit sufficiently long excited-state lifetimes in the locally-excited states to undergo transition to the productive charge transfer state. These findings can be further used in the design of similar photoredox reactions which can employ the potential of chalcogen bonding interactions.

Hints of Biological Activity of Xerosydryle: Preliminary Evidence on the Early Stages of Seedling Development

Xerosydryle belongs to a new category of materials resulting from the interaction of water with various hydrophilic polymers. These materials can exhibit different properties depending on the kind of polymer-water interaction. Previous research confirmed the existence of a solid manifestation of water at room temperature. The thermal properties of dissolved xerosydryle in water are similar to those of biological macromolecules during denaturation but with greater stability. This study investigated the biological effect of xerosydryle on a living system for the first time, using a seed germination model. The interaction was evaluated using physiological assays such as chlorophyll shifts, potassium (re)uptake during the onset of germination and a transcriptome approach. Seeds were treated with samples of xerosydryle and distilled water. Transcriptome analysis of germinating seeds highlighted differences (up- and down-regulated genes) between seeds treated with xerosydryle and those treated with distilled water. Overall, the experiments performed indicate that xerosydryle, even at low concentrations, interferes with seedling growth in a manner similar to an osmotic modulator. This work paves the way for a more comprehensive exploration of the active biological role of xerosydryle and similar compounds on living matter and opens up speculation on the interactions at the boundaries between physics, chemistry, and biology.

Nuclear Quantum Effects on the Electronic Structure of Water and Ice

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

Direct Observation of Anionic Yields from the Liquid-Vapor Interface by Electron Irradiation

Dissociative electron attachment (DEA) is widely believed to play a high-profile role in ionizing radiation damages of bioorganic molecules, and its fundamentals are mainly learned from the gas-phase studies. However, the DEA process in aqueous solution is still in debate. Here we provide experimental evidence about the DEA processes of liquid methanol by using electron-impact-time-delayed mass spectrometry. In contrast to the gas- and solid-phase DEAs, methoxide ion CH3O- is the predominant product from the liquid interface. Furthermore, this anion can be produced with both the primary low-energy electrons and the inelastically scattered and secondary low-energy electrons. On the contrary, the primary low-energy electrons in the liquid bulk are more likely to be solvated, rather than directly participating in the DEA process. Our study provides new insights into radiation chemistry, particularly of bioorganic relevance.

Photo-excited charge transfer from adamantane to electronic bound states in water

Aqueous nanodiamonds illuminated by UV light produce free solvated electrons, which may drive high-energy reduction reactions in water. However, the influence of water conformations on the excited-state electron-transfer mechanism are still under debate. In this work, we offer a theoretical study of charge-transfer states in adamantane-water structures obtained by linear-response time-dependent density-functional theory. Small water clusters with broken hydrogen bonds are found to efficiently bind the electron from adamantane. A distinction is made with respect to the nature of the water clusters: some bind the electron in a water cavity, others along a strong permanent total dipole. These two types of bound states are more strongly binding, the higher their electron affinity and their positive electrostatic potential, the latter being dominated by the energy of the lowest unoccupied molecular orbital of the isolated water clusters. Structural sampling in a thermal equilibrium at room temperature via molecular dynamics snapshots confirms under which conditions the underlying waters clusters can occur and verifies that broken hydrogen bonds in the water network close to adamantane can create traps for the solvated electron.

Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections

We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.

Quadratic scaling bosonic path integral molecular dynamics

Bosonic exchange symmetry leads to fascinating quantum phenomena, from exciton condensation in quantum materials to the superfluidity of liquid 4He. Unfortunately, path integral molecular dynamics (PIMD) simulations of bosons are computationally prohibitive beyond ∼100 particles, due to a cubic scaling with the system size. We present an algorithm that reduces the complexity from cubic to quadratic, allowing the first simulations of thousands of bosons using PIMD. Our method is orders of magnitude faster, with a speedup that scales linearly with the number of particles and the number of imaginary time slices (beads). Simulations that would have otherwise taken decades can now be done in days. In practice, the new algorithm eliminates most of the added computational cost of including bosonic exchange effects, making them almost as accessible as PIMD simulations of distinguishable particles.

Use of FLOSIC for understanding anion-solvent interactions

An Achille's heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3-, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth's atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi-Löwdin-Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree-Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Flow through Deformed Carbon Nanotubes Predicted by Rigid and Flexible Water Models

In this study, using nonequilibrium molecular dynamics simulation, the flow of water in deformed carbon nanotubes is studied for two water models TIP4P/2005 and simple point charge/FH (SPC/FH). The results demonstrated a nonuniform dependence of the flow on the tube deformation and the flexibility imposed on the water molecules, leading to an unexpected increase in the flow in some cases. The effects of the tube diameter and pressure gradient are investigated to explain the abnormal flow behavior with different degrees of structural deformation.

The Back Door to the Surface Hydrated Electron

We use a Mg+ metal to extend the size regime of aqueous clusters to extrapolate to the bulk limit of the vertical detachment energy (VDE) of the solvated electron to >3,200, a value between 1 to over 2 orders of magnitude larger than the one previously measured experimentally or computed theoretically. We relate the VDE to the energy difference between the Mg+(H2O)n and Mg2+(H2O)n systems and the metal's second ionization potential. The extrapolated bulk VDEs of the localized surface electron, which moves away from the metal as n increases, are 1.89 ± 0.01 eV for semiempirical (n ∼ 3,200; PM6-D3H4) and 1.73 ± 0.03 eV (n ∼ 150; HF) and 1.83 ± 0.02 eV (n ∼ 150; MP2) for ab initio, in excellent agreement with the 1.6-1.8 eV range of experimental results. The VDEs converge from above (larger values) to the bulk limit, in a manner that is qualitatively opposite from previous studies and experiments lacking a charged metal, a fact justifying the "back door" approach to the solvated electron.

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

Electronic structure calculations based on density functional theory (DFT) have successfully predicted numerous ground-state properties of a variety of molecules and materials. However, exchange and correlation functionals currently used in the literature, including semilocal and hybrid functionals, are often inaccurate to describe the electronic properties of heterogeneous solids, especially systems composed of building blocks with large dielectric mismatch. Here, we present a dielectric-dependent range-separated hybrid functional, screened-exchange range-separated hybrid (SE-RSH), for the investigation of heterogeneous materials. We define a spatially dependent fraction of exact exchange inspired by the static Coulomb-hole and screened-exchange (COHSEX) approximation used in many-body perturbation theory, and we show that the proposed functional accurately predicts the electronic structure of several nonmetallic interfaces, three- and two-dimensional, pristine, and defective solids and nanoparticles.

MBX: A many-body energy and force calculator for data-driven many-body simulations

Many-Body eXpansion (MBX) is a C++ library that implements many-body potential energy functions (PEFs) within the "many-body energy" (MB-nrg) formalism. MB-nrg PEFs integrate an underlying polarizable model with explicit machine-learned representations of many-body interactions to achieve chemical accuracy from the gas to the condensed phases. MBX can be employed either as a stand-alone package or as an energy/force engine that can be integrated with generic software for molecular dynamics and Monte Carlo simulations. MBX is parallelized internally using Open Multi-Processing and can utilize Message Passing Interface when available in interfaced molecular simulation software. MBX enables classical and quantum molecular simulations with MB-nrg PEFs, as well as hybrid simulations that combine conventional force fields and MB-nrg PEFs, for diverse systems ranging from small gas-phase clusters to aqueous solutions and molecular fluids to biomolecular systems and metal-organic frameworks.

The Spontaneous Electron-Mediated Redox Processes on Sprayed Water Microdroplets

Water is considered as an inert environment for the dispersion of many chemical systems. However, by simply spraying bulk water into microsized droplets, the water microdroplets have been shown to possess a large plethora of unique properties, including the ability to accelerate chemical reactions by several orders of magnitude compared to the same reactions in bulk water, and/or to trigger spontaneous reactions that cannot occur in bulk water. A high electric field (∼10⁹ V/m) at the air-water interface of microdroplets has been postulated to be the probable cause of the unique chemistries. This high field can even oxidize electrons out of hydroxide ions or other closed-shell molecules dissolved in water, forming radicals and electrons. Subsequently, the electrons can trigger further reduction processes. In this Perspective, by showing a large number of such electron-mediated redox reactions, and by studying the kinetics of these reactions, we opine that the redox reactions on sprayed water microdroplets are essentially processes using electrons as the charge carriers. The potential impacts of the redox capability of microdroplets are also discussed in a larger context of synthetic chemistry and atmospheric chemistry.

Many-body screening effects in liquid water

The screening arising from many-body excitations is a crucial quantity for describing absorption and inelastic X-ray scattering (IXS) of materials. Similarly, the electron screening plays a critical role in state-of-the-art approaches for determining the fundamental band gap. However, ab initio studies of the screening in liquid water have remained limited. Here, we use a combined analysis based on the Bethe-Salpeter equation and time-dependent density functional theory. We first show that absorption spectra at near-edge energies are insufficient to assess the accuracy by which the screening is described. Next, when the energy range under scrutiny is extended, we instead find that the IXS spectra are highly sensitive and allow for the selection of the optimal theoretical scheme. This leads to good agreement with experiment over a large range of transferred energies and momenta, and enables establishing the elusive fundamental band gap of liquid water at 9.3 eV.

Hydrated electrons as nodes in porous clathrate hydrates

We investigate the structures of hydrated electrons (e^- _aq) in one of water's solid phases, namely, clathrate hydrates (CHs). Using density functional theory (DFT) calculations, DFT-based ab initio molecular dynamics (AIMD), and path-integral AIMD simulations with periodic boundary conditions, we find that the structure of the e^- _aq@node model is in good agreement with the experiment, suggesting that an e^- _aq could form a node in CHs. The node is a H₂O defect in CHs that is supposed to be composed of four unsaturated hydrogen bonds. Since CHs are porous crystals that possess cavities that can accommodate small guest molecules, we expect that these guest molecules can be used to tailor the electronic structure of the e^- _aq@node, and it leads to experimentally observed optical absorption spectra of CHs. Our findings have a general interest and extend the knowledge of e^- _aq into porous aqueous systems.

Surface or Internal Hydration - Does It Really Matter?

The precise location of an ion or electron, whether it is internally solvated or residing on the surface of a water cluster, remains an intriguing question. Subtle differences in the hydrogen bonding network may lead to a preference for one or the other. Here we discuss spectroscopic probes of the structure of gas-phase hydrated ions in combination with quantum chemistry, as well as H/D exchange as a means of structure elucidation. With the help of nanocalorimetry, we look for thermochemical signatures of surface vs internal solvation. Examples of strongly size-dependent reactivity are reviewed which illustrate the influence of surface vs internal solvation on unimolecular rearrangements of the cluster, as well as on the rate and product distribution of ion-molecule reactions.

Absolute Electronic Energetics and Quantitative Work Functions of Liquids from Photoelectron Spectroscopy

Liquid-jet photoelectron spectroscopy (LJ-PES) enabled a breakthrough in the experimental study of the electronic structure of liquid water, aqueous solutions, and volatile liquids more generally. The novelty of this technique, dating back over 25 years, lies in stabilizing a continuous, micron-diameter LJ in a vacuum environment to enable PES studies. A key quantity in PES is the most probable energy associated with vertical promotion of an electron into vacuum: the vertical ionization energy, VIE, for neutrals and cations, or vertical detachment energy, VDE, for anions. These quantities can be used to identify species, their chemical states and bonding environments, and their structural properties in solution. The ability to accurately measure VIEs and VDEs is correspondingly crucial. An associated principal challenge is the determination of these quantities with respect to well-defined energy references. Only with recently developed methods are such measurements routinely and generally viable for liquids. Practically, these methods involve the application of condensed-matter concepts to the acquisition of photoelectron (PE) spectra from liquid samples, rather than solely relying on molecular-physics treatments that have been commonly implemented since the first LJ-PES experiments. This includes explicit consideration of the traversal of electrons to and through the liquid's surface, prior to free-electron detection. Our approach to measuring VIEs and VDEs with respect to the liquid vacuum level specifically involves detecting the lowest-energy electrons emitted from the sample, which have barely enough energy to surmount the surface potential and accumulate in the low-energy tail of the liquid-phase spectrum. By applying a sufficient bias potential to the liquid sample, this low-energy spectral tail can generally be exposed, with its sharp, low-energy cutoff revealing the genuine kinetic-energy-zero in a measured spectrum, independent of any perturbing intrinsic or extrinsic potentials in the experiment. Together with a precisely known ionizing photon energy, this feature enables the straightforward determination of VIEs or VDEs, with respect to the liquid-phase vacuum level, from any PE feature of interest. Furthermore, by additionally determining solution-phase VIEs and VDEs with respect to the common equilibrated energy level in condensed matter, the Fermi level─the generally implemented reference energy in solid-state PES─solution work functions, eΦ, and liquid-vacuum surface dipole effects can be quantified. With LJs, the Fermi level can only be properly accessed by controlling unwanted surface charging and all other extrinsic potentials, which lead to energy shifts of all PE features and preclude access to accurate electronic energetics. More specifically, conditions must be engineered to minimize all undesirable potentials, while maintaining the equilibrated, intrinsic (contact) potential difference between the sample and apparatus. The establishment of these liquid-phase, accurate energy-referencing protocols importantly enables VIE and VDE determinations from near-arbitrary solutions and the quantitative distinction between bulk electronic structure and interfacial effects. We will review and exemplify these protocols for liquid water and several exemplary aqueous solutions here, with a focus on the lowest-ionization- or lowest-detachment-energy PE peaks, which importantly relate to the oxidative stabilities of aqueous-phase species.

Early dynamics of the emission of solvated electrons from nanodiamonds in water

Solvated electrons are among the most reductive species in an aqueous environment. Diamond materials have been proposed as a promising source of solvated electrons, but the underlying emission process in water remains elusive so far. Here, we show spectroscopic evidence for the emission of solvated electrons from detonation nanodiamonds upon excitation with both deep ultraviolet (225 nm) and visible (400 nm) light using ultrafast transient absorption. The crucial role of surface termination in the emission process is evidenced by comparing hydrogenated, hydroxylated and carboxylated nanodiamonds. In particular, a transient response that we attribute to solvated electrons is observed on hydrogenated nanodiamonds upon visible light excitation, while it shows a sub-ps recombination due to trap states when excited with deep ultraviolet light. The essential role of surface reconstructions on the nanodiamonds in these processes is proposed based on density functional theory calculations. These results open new perspectives for solar-driven emission of solvated electrons in an aqueous phase using nanodiamonds.

Observation of a transient intermediate in the ultrafast relaxation dynamics of the excess electron in strong-field-ionized liquid water

A unified picture of the electronic relaxation dynamics of ionized liquid water has remained elusive despite decades of study. Here, we employ sub-two-cycle visible to short-wave infrared pump-probe spectroscopy and ab initio nonadiabatic molecular dynamics simulations to reveal that the excess electron injected into the conduction band (CB) of ionized liquid water undergoes sequential relaxation to the hydrated electron s ground state via an intermediate state, identified as the elusive p excited state. The measured CB and p-electron lifetimes are 0.26 ± 0.02 ps and 62 ± 10 fs, respectively. Ab initio quantum dynamics yield similar lifetimes and furthermore reveal vibrational modes that participate in the different stages of electronic relaxation, with initial relaxation within the dense CB manifold coupled to hindered translational motions whereas subsequent p-to-s relaxation facilitated by librational and even intramolecular bending modes of water. Finally, energetic considerations suggest that a hitherto unobserved trap state resides ~0.3-eV below the CB edge of liquid water. Our results provide a detailed atomistic picture of the electronic relaxation dynamics of ionized liquid water with unprecedented time resolution.

Hydration Structure of Na+ and K+ Ions in Solution Predicted by Data-Driven Many-Body Potentials

The hydration structure of Na⁺ and K⁺ ions in solution is systematically investigated using a hierarchy of molecular models that progressively include more accurate representations of many-body interactions. We found that a conventional empirical pairwise additive force field that is commonly used in biomolecular simulations is unable to reproduce the extended X-ray absorption fine structure (EXAFS) spectra for both ions. In contrast, progressive inclusion of many-body effects rigorously derived from the many-body expansion of the energy allows the MB-nrg potential energy functions (PEFs) to achieve nearly quantitative agreement with the experimental EXAFS spectra, thus enabling the development of a molecular-level picture of the hydration structure of both Na⁺ and K⁺ in solution. Since the MB-nrg PEFs have already been shown to accurately describe isomeric equilibria and vibrational spectra of small ion-water clusters in the gas phase, the present study demonstrates that the MB-nrg PEFs effectively represent the long-sought-after models able to correctly predict the properties of ionic aqueous systems from the gas to the liquid phase, which has so far remained elusive.

First-principles redox energy estimates under the condition of satisfying the general form of Koopmans’ theorem: An atomistic study of aqueous iron

In this work, we explore the relative accuracy to which a hybrid functional, in the context of density functional theory, may predict redox properties under the constraint of satisfying the general form of Koopmans’ theorem. Taking aqueous iron as our model system within the framework of first-principles molecular dynamics, direct comparison between computed single-particle energies and experimental ionization data is assessed by both (1) tuning the degree of hybrid exchange, to satisfy the general form of Koopmans’ theorem, and (2) ensuring the application of finite-size corrections. These finite-size corrections are benchmarked through classical molecular dynamics calculations, extended to large atomic ensembles, for which good convergence is obtained in the large supercell limit. Our first-principles findings indicate that while precise quantitative agreement with experimental ionization data cannot always be attained for solvated systems, when satisfying the general form of Koopmans’ theorem via hybrid functionals, theoretically robust estimates of single-particle redox energies are most often arrived at by employing a total energy difference approach. That is, when seeking to employ a value of exact exchange that does not satisfy the general form of Koopmans’ theorem, but some other physical metric, the single-particle energy estimate that would most closely align with the general form of Koopmans’ theorem is obtained from a total energy difference approach. In this respect, these findings provide important guidance for the more general comparison of redox energies computed via hybrid functionals with experimental data.

An Algorithmic Approach Based on Data Trees and Genetic Algorithms to Understanding Charged and Neutral Metal Nanocluster Growth

Metal nanocluster growth pathways are challenging to study due to the number of possible species involved and the branching nature of possible mechanisms. We present a data tree-based approach to mapping these reaction pathways based on structures and energies obtained from Density Functional Theory (DFT) computations and including positive, negative, and neutral clusters in a continuum solvent of water. We also develop a genetic algorithm to study the relative stability of the clusters, which we combine with the data tree method to determine the effects that favorable decomposition reactions have on the reaction pathways. Introducing data tree pruning based on the exothermicity of each reaction, and including the first and second most important paths for each cluster up to five atoms including positive, neutral, and negative clusters, is then implemented to determine the cluster growth paths. These most favorable Ag-Ag reaction pathways are in agreement with more limited prior theoretical and experimental results, but they provide more systematic results that include predictions about the importance of clusters not previously identified. A key feature of the data tree approach is that it provides a comprehensive sampling of possible clusters, but without needing to generate the entire reaction tree. Additionally, the data tree-based approach allows for flexibility in the analysis based on restricting reactant charge or size, selecting a starting species to mimic an experimental precursor, and choosing a maximum allowable final cluster size or charge of interest.

GPU Acceleration of Large-Scale Full-Frequency GW Calculations

Many-body perturbation theory is a powerful method to simulate electronic excitations in molecules and materials starting from the output of density functional theory calculations. By implementing the theory efficiently so as to run at scale on the latest leadership high-performance computing systems it is possible to extend the scope of GW calculations. We present a GPU acceleration study of the full-frequency GW method as implemented in the WEST code. Excellent performance is achieved through the use of (i) optimized GPU libraries, e.g., cuFFT and cuBLAS, (ii) a hierarchical parallelization strategy that minimizes CPU-CPU, CPU-GPU, and GPU-GPU data transfer operations, (iii) nonblocking MPI communications that overlap with GPU computations, and (iv) mixed precision in selected portions of the code. A series of performance benchmarks has been carried out on leadership high-performance computing systems, showing a substantial speedup of the GPU-accelerated version of WEST with respect to its CPU version. Good strong and weak scaling is demonstrated using up to 25 920 GPUs. Finally, we showcase the capability of the GPU version of WEST for large-scale, full-frequency GW calculations of realistic systems, e.g., a nanostructure, an interface, and a defect, comprising up to 10 368 valence electrons.

Toward High-level Machine Learning Potential for Water Based on Quantum Fragmentation and Neural Networks

Accurate and efficient simulation of liquids, such as water and salt solutions, using high-level wave function theories is still a formidable task for computational chemists owing to the high computational costs. In this study, we develop a deep machine learning potential based on fragment-based second-order Møller-Plesset perturbation theory (DP-MP2) for water through neural networks. We show that the DP-MP2 potential predicts the structural, dynamical, and thermodynamic properties of liquid water in better agreement with the experimental data than previous studies based on density functional theory (DFT). The nuclear quantum effects (NQEs) on the properties of liquid water are also examined, which are noticeable in affecting the structural and dynamical properties of liquid water under ambient conditions. This work provides a general framework for quantitative predictions of the properties of condensed-phase systems with the accuracy of high-level wave function theory while achieving significant computational savings compared to ab initio simulations.

Photoionization of the aqueous phase: clusters, droplets and liquid jets

This perspective article reviews specific challenges associated with photoemission spectroscopy of bulk liquid water, aqueous solutions, water droplets and water clusters. The main focus lies on retrieving accurate energetics and photoelectron angular information from measured photoemission spectra, and on the question how these quantities differ in different aqueous environments. Measured photoelectron band shapes, vertical binding energies (ionization energies), and photoelectron angular distributions are influenced by various phenomena. We discuss the influences of multiple energy-dependent electron scattering in aqueous environments, and we discuss different energy referencing methods, including the application of a bias voltage to access absolute energetics of solvent and solute. Recommendations how to account for or minimize the influence of electron scattering are provided. The example of the hydrated electron in different aqueous environments illustrates how one can account for electron scattering, while reliable methods addressing parasitic potentials and proper energy referencing are demonstrated for ionization from the outermost valence orbital of neat liquid water.

Review-Hysteresis in Carbon Nano-Structure Field Effect Transistor

In recent decades, the research of nano-structure devices (e.g., carbon nanotube and graphene) has experienced rapid growth. These materials have supreme electronic, thermal, optical and mechanical properties and have received widespread concern in different fields. It is worth noting that gate hysteresis behavior of field effect transistors can always be found in ambient conditions, which may influence the transmission appearance. Many researchers have put forward various views on this question. Here, we summarize and discuss the mechanisms behind hysteresis, different influencing factors and improvement methods which help decrease or eliminate unevenness and asymmetry.

Uptake of N2O5 by aqueous aerosol unveiled using chemically accurate many-body potentials

The reactive uptake of N2O5 to aqueous aerosol is a major loss channel for nitrogen oxides in the troposphere. Despite its importance, a quantitative picture of the uptake mechanism is missing. Here we use molecular dynamics simulations with a data-driven many-body model of coupled-cluster accuracy to quantify thermodynamics and kinetics of solvation and adsorption of N2O5 in water. The free energy profile highlights that N2O5 is selectively adsorbed to the liquid-vapor interface and weakly solvated. Accommodation into bulk water occurs slowly, competing with evaporation upon adsorption from gas phase. Leveraging the quantitative accuracy of the model, we parameterize and solve a reaction-diffusion equation to determine hydrolysis rates consistent with experimental observations. We find a short reaction-diffusion length, indicating that the uptake is dominated by interfacial features. The parameters deduced here, including solubility, accommodation coefficient, and hydrolysis rate, afford a foundation for which to consider the reactive loss of N2O5 in more complex solutions.

Low-energy electron distributions from the photoionization of liquid water: a sensitive test of electron mean free paths

The availability of accurate mean free paths for slow electrons (<50 eV) in water is central to the understanding of many electron-driven processes in aqueous solutions, but their determination poses major challenges to experiment and theory alike. Here, we describe a joint experimental and theoretical study demonstrating a novel approach for testing, and, in the future, refining such mean free paths. We report the development of Monte-Carlo electron-trajectory simulations including elastic and inelastic electron scattering, as well as energy loss and secondary-electron production to predict complete photoelectron spectra of liquid water. These simulations are compared to a new set of photoelectron spectra of a liquid-water microjet recorded over a broad range of photon energies in the extreme ultraviolet (20-57 eV). Several previously published sets of scattering parameters are investigated, providing direct and intuitive insights on how they influence the shape of the low-energy electron spectra. A pronounced sensitivity to the escape barrier is also demonstrated. These simulations considerably advance our understanding of the origin of the prominent low-energy electron distributions in photoelectron spectra of liquid water and clarify the influence of scattering parameters and the escape barrier on their shape. They moreover describe the reshaping and displacement of low-energy photoelectron bands caused by vibrationally inelastic scattering. Our work provides a quantitative basis for the interpretation of the complete photoelectron spectra of liquids and opens the path to fully predictive simulations of low-energy scattering in liquid water.

Electronic structure of the solvated benzene radical anion

The benzene radical anion is a molecular ion pertinent to several organic reactions, including the Birch reduction of benzene in liquid ammonia. The species exhibits a dynamic Jahn-Teller effect due to its open-shell nature and undergoes pseudorotation of its geometry. Here, we characterize the complex electronic structure of this condensed-phase system based on ab initio molecular dynamics simulations and GW calculations of the benzene radical anion solvated in liquid ammonia. Using detailed analysis of the molecular and electronic structure, we find that the spatial character of the excess electron of the solvated radical anion follows the underlying Jahn-Teller distortions of the molecular geometry. We decompose the electronic density of states to isolate the contribution of the solute and to examine the response of the solvent to its presence. Our findings show the correspondence between instantaneous molecular structure and spin density; provide important insights into the electronic stability of the species, revealing that it is, indeed, a bound state in the condensed phase; and offer electronic densities of states that aid in the interpretation of experimental photoelectron spectra.

Shallow and deep trap states of solvated electrons in methanol and their formation, electronic excitation, and relaxation dynamics

We present condensed-phase first-principles molecular dynamics simulations to elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons (e_met⁻) in methanol. Excess electrons injected into liquid methanol are most likely trapped by methyl groups, but rapidly diffuse to more stable trapping sites with dangling OH bonds. After localization at the sites with one free OH bond (1OH trapping sites), reorientation of other methanol molecules increases the OH coordination number and the trap depth, and ultimately four OH bonds become coordinated with the excess electrons under thermal conditions. The simulation identified four distinct trapping states with different OH coordination numbers. The simulation results also revealed that electronic transitions of e_met⁻ are primarily due to charge transfer between electron trapping sites (cavities) formed by OH and methyl groups, and that these transitions differ from hydrogenic electronic transitions involving aqueous solvated electrons (e_aq⁻). Such charge transfer also explains the alkyl-chain-length dependence of the photoabsorption peak wavelength and the excited-state lifetime of solvated electrons in primary alcohols.

Condensed-phase first-principles molecular dynamics simulations elucidate the presence of different electron trapping sites in liquid methanol and their roles in the formation, electronic transitions, and relaxation of solvated electrons.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

Data-driven many-body models enable a quantitative description of chloride hydration from clusters to bulk

We present a new data-driven potential energy function (PEF) describing chloride-water interactions, which is developed within the many-body-energy (MB-nrg) theoretical framework. Besides quantitatively reproducing low-order many-body energy contributions, the new MB-nrg PEF is able to correctly predict the interaction energies of small chloride-water clusters calculated at the coupled cluster level of theory. Importantly, classical and quantum molecular dynamics simulations of a single chloride ion in water demonstrate that the new MB-nrg PEF predicts x-ray spectra in close agreement with the experimental results. Comparisons with an popular empirical model and a polarizable PEF emphasize the importance of an accurate representation of short-range many-body effect while demonstrating that pairwise additive representations of chloride-water and water-water interactions are inadequate for correctly representing the hydration structure of chloride in both gas-phase clusters and solution. We believe that the analyses presented in this study provide additional evidence for the accuracy and predictive ability of the MB-nrg PEFs, which can then enable more realistic simulations of ionic aqueous systems in different environments.

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.

The Fermi level as an energy reference in liquid jet X-ray photoelectron spectroscopy studies of aqueous solutions

To advance the understanding of key electrochemical and photocatalytic processes that depend on the electronic structure of aqueous solutions, X-ray photoemission spectroscopy has become an invaluable tool, especially when practiced with liquid microjet setups. Determining vertical ionization energies referenced to the vacuum level, and binding energies referenced to the Fermi level, including the much-coveted reorganization energy of the oxidized species of a redox couple, requires that energy levels be properly defined. The present paper addresses specifically how the vacuum level "just outside the surface" can be known through the energy position of the rising edge of the secondary electrons, and how the Fermi level reference is uniquely determined via the introduction of a redox couple. Taking the case of the ferricyanide/ferrocyanide and ferric/ferrous couples, this study also tackles issues related to the electrokinetic effects inherent to the production of a liquid jet in a vacuum, which has become the standard water sample environment for photoemission experiments.

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.

Relationship between Hydrogen-Bonding Motifs and the 1b1 Splitting in the X-ray Emission Spectrum of Liquid Water

The split of the 1b₁ peak observed in the X-ray emission (XE) spectrum of liquid water has been the focus of intense research. Although several hypotheses have been proposed to explain the origin of this split, a consensus has not yet been reached. Here, we introduce a novel theoretical/computation approach which, combining path-integral molecular dynamics simulations with the MB-pol model and time-dependent density functional theory calculations, predicts the 1b₁ splitting in liquid water and not in crystalline ice, in agreement with the experimental observations. A systematic analysis of the underlying local structure of liquid water at ambient conditions indicates that several different hydrogen-bonding motifs contribute to the overall XE line shape in the energy range corresponding to emissions from the 1b₁ orbitals. This suggests that it is not possible to unambiguously attribute the split of the 1b₁ peak to only two specific structural arrangements of the underlying hydrogen-bonding network.

In situ monitoring of the influence of water on DNA radiation damage by near-ambient pressure X-ray photoelectron spectroscopy

Ionizing radiation damage to DNA plays a fundamental role in cancer therapy. X-ray photoelectron-spectroscopy (XPS) allows simultaneous irradiation and damage monitoring. Although water radiolysis is essential for radiation damage, all previous XPS studies were performed in vacuum. Here we present near-ambient-pressure XPS experiments to directly measure DNA damage under water atmosphere. They permit in-situ monitoring of the effects of radicals on fully hydrated double-stranded DNA. The results allow us to distinguish direct damage, by photons and secondary low-energy electrons (LEE), from damage by hydroxyl radicals or hydration induced modifications of damage pathways. The exposure of dry DNA to x-rays leads to strand-breaks at the sugar-phosphate backbone, while deoxyribose and nucleobases are less affected. In contrast, a strong increase of DNA damage is observed in water, where OH-radicals are produced. In consequence, base damage and base release become predominant, even though the number of strand-breaks increases further.

Low-energy constraints on photoelectron spectra measured from liquid water and aqueous solutions

We report on the effects of electron collision and indirect ionization processes, occurring at photoexcitation and electron kinetic energies well below 30 eV, on the photoemission spectra of liquid water. We show that the nascent photoelectron spectrum and, hence, the inferred electron binding energy can only be accurately determined if electron energies are large enough that cross sections for quasi-elastic scattering processes, such as vibrational excitation, are negligible. Otherwise, quasi-elastic scattering leads to strong, down-to-few-meV kinetic energy scattering losses from the direct photoelectron features, which manifest in severely distorted intrinsic photoelectron peak shapes. The associated cross-over point from predominant (known) electronically inelastic to quasi-elastic scattering seems to arise at surprisingly large electron kinetic energies, of approximately 10-14 eV. Concomitantly, we present evidence for the onset of indirect, autoionization phenomena (occurring via superexcited states) within a few eV of the primary and secondary ionization thresholds. These processes are inferred to compete with the direct ionization channels and primarily produce low-energy photoelectrons at photon and electron impact excitation energies below ∼15 eV. Our results highlight that vibrational inelastic electron scattering processes and neutral photoexcitation and autoionization channels become increasingly important when photon and electron kinetic energies are decreased towards the ionization threshold. Correspondingly, we show that for neat water and aqueous solutions, great care must be taken when quantitatively analyzing photoelectron spectra measured too close to the ionization threshold. Such care is essential for the accurate determination of solvent and solute ionization energies as well as photoelectron branching ratios and peak magnitudes.

Machine learning dielectric screening for the simulation of excited state properties of molecules and materials

Accurate and efficient calculations of absorption spectra of molecules and materials are essential for the understanding and rational design of broad classes of systems. Solving the Bethe-Salpeter equation (BSE) for electron-hole pairs usually yields accurate predictions of absorption spectra, but it is computationally expensive, especially if thermal averages of spectra computed for multiple configurations are required. We present a method based on machine learning to evaluate a key quantity entering the definition of absorption spectra: the dielectric screening. We show that our approach yields a model for the screening that is transferable between multiple configurations sampled during first principles molecular dynamics simulations; hence it leads to a substantial improvement in the efficiency of calculations of finite temperature spectra. We obtained computational gains of one to two orders of magnitude for systems with 50 to 500 atoms, including liquids, solids, nanostructures, and solid/liquid interfaces. Importantly, the models of dielectric screening derived here may be used not only in the solution of the BSE but also in developing functionals for time-dependent density functional theory (TDDFT) calculations of homogeneous and heterogeneous systems. Overall, our work provides a strategy to combine machine learning with electronic structure calculations to accelerate first principles simulations of excited-state properties.

Molecular Dynamics Characterization of Dielectron Hydration in Liquid Water with Unique Double Proton Transfers

Radiation chemistry of water and aqueous solutions has always been an interesting scientific issue owing to involving electronic excitations, ionization of solvated species, and formation of radiolytic species and many elementary reactions, but the underlying mechanisms are still poorly understood. Here, we for the first time molecular dynamics characterize the hydration dynamics of two correlated electrons and their triggered unique phenomena in liquid water associated with radiolysis of water using the combined hybrid functional and nonlocal dispersion functional. Hydration of two electrons may experience two distinctly different mechanisms, one forming a spin-paired closed-shell unicaged dielectron hydrate (e₂^2-_aq) and the other forming a spin-paired metastable open-shell bicaged hydrated electron pair (e^-_aq···e^-_aq) which exhibits intriguing antiferromagnetic spin coupling dynamics (in a range of -40 cm^-1 to -500 cm^-1). e^-_aq···e^-_aq can recombine to e₂^2-_aq through a unique solvent fluctuation-controlled gradual-flowing mechanism, and enlarging fluctuation can promote the conversion. Interestingly, we directly observe that e₂^2-_aq as the precursor can trigger hydrogen evolution via unique continuous spontaneous double proton transfer to the dielectron with a short-lived H^-_aq intermediate, but e^-_aq···e^-_aq does not directly. This is the first direct observation for the connection between e₂^2-_aq and spontaneous hydrogen evolution including participation of H^-_aq in aqueous solution, bridging relevant experimental phenomena. This work also evidences an unnoticed process, the double proton transfer mediated charge separation, and presents the first detailed analysis regarding the evolution dynamics of e₂^2-_aq for the understanding of the radiolysis reactions in aqueous solutions.

Identification of novel potential cyclooxygenase-2 inhibitors using ligand- and structure-based virtual screening approaches

Cyclooxygenase 2 (COX-2) is a well-established target for the design of anti-inflammatory intermediates. Celecoxib was selected as a template molecule to perform ligand-based virtual screening, i.e. to search for structures with similarity in shape and electrostatic potential, with a gradual increase in accuracy through the combined fitting of several steps using eight commercial databases. The molecules ZINC408709 and ZINC2090319 reproduced values within the limits established in an initial study of absorption and distribution in the body. No alert was fired for possible toxic groups when these molecules were subjected to toxicity prediction. Molecular docking results with these compounds showed a higher binding affinity in comparison to rofecoxib for the COX-2 target. Additionally, ZINC408709 and ZINC2090319 were predicted to be potentially biologically active. In in silico prediction of endocrine disruption potential, it was established that the molecule ZINC2090319 binds strongly to the target related to cardiovascular risk in a desirable way as a non-steroidal antagonist and ZINC408709 binds strongly to the target that is associated with the treatment of inflammatory pathologies and similar to celecoxib. Metabolites generated from these compounds are less likely to have side effects. Simulations were used to evaluate the interaction of compounds with COX-1 and COX-2 during 200 ns. Despite the differences, ZINC408709 molecule showed better stability for COX-2 during molecular dynamics simulation. In the calculations of free energy MM/PBSA, the molecule ZINC408709 ΔGbind value has a higher affinity to celecoxib and rofecoxib COX-2. This demonstrates that the selected substances can be considered as promising COX-2 inhibitors. Communicated by Ramaswamy H. Sarma.

Below Band Gap Formation of Solvated Electrons in Neutral Water Clusters?

Below band gap formation of solvated electrons in neutral water clusters using pump-probe photoelectron imaging is compared with recent data for liquid water and with above band gap excitation studies in liquid and clusters. Similar relaxation times on the order of 200 fs and 1-2 ps are retrieved for below and above band gap excitation, in both clusters and liquid. The independence of the relaxation times from the generation process indicates that these times are dominated by the solvent response, which is significantly slower than the various solvated electron formation processes. The analysis of the temporal evolution of the vertical electron binding energy and the electron binding energy at half-maximum suggests a dependence of the solvation time on the binding energy.