Search

We performed a computational study of positron attachment to hydrated amino acids, namely glycine, alanine, and proline in the zwitterionic form. We combined the sequential quantum mechanics/molecular mechanics (s-QM/MM) method with various levels of any particle molecular orbital (APMO) calculations. Consistent with previous studies, our calculations indicate the formation of energetically stable states for the isolated and microsolvated amino acids, in which the positron localizes around the carboxylate group. However, for the larger clusters, composed of 7 to 40 water molecules, hydrogen bonding between the solute and solvent molecules disfavors positron attachment to the amino acids, giving rise to surface states in which the positron is located around the water–vacuum interface. The analysis of positron binding energies, positronic orbitals, radial probability distributions, and annihilation rates consistently pointed out the change from positron–solute to positron–solvent states. Even with the inclusion of an electrostatic embedding around the aggregates, the positrons did not localize around the solute. Positron attachment to molecules in the gas phase is a well-established fact. The existence of hydrated positronic molecules could also be expected from the analogy with transient anion states, which are believed to participate in radiation damage. Our results indicate that positron attachment to hydrated biomolecules, even to zwitterions with negatively charged carboxylated groups, would not take place. For the larger clusters, in which positron–water interactions are favored, the calculations indicate an unexpectedly large contribution of the core orbitals to the annihilation rates, between 15 and 20%. Finally, we explored correlations between positron binding energies (PBEs) and dipole moments, as well as annihilation rates and PBEs, consistent with previous studies for smaller clusters.

8-oxo-Guanine is a mutagenic lesion produced by reactions involving reactive oxygen species and guanine in DNA. Its production induces mispairing between the canonical nucleobases during DNA replication such that various types of cancers are associated with the DNA lesion. Since radiation therapy is used in some cases, the interaction of low-energy electrons with 8-oxo-guanine can in turn produce other reactive species, which in principle could have either a detrimental or protective effect on the organism. Motivated by these facts, we report a comparative experimental study of electron-induced fragmentation of guanine and 8-oxo-guanine, along with a theoretical study of the π* shape resonances and bound anion states, which may trigger those dissociation reactions. The electron-induced fragmentation of 8-oxo-guanine is remarkably distinct from the native form. More complex reactions were observed for the oxidized species, which may produce several anion fragments at very low energies (∼0 eV). The dehydrogenated parent anion, which is already a minor fragment in guanine, was completely suppressed in 8-oxo-guanine. The calculated thermodynamical thresholds also suggest that NH2elimination in guanine, at sub-excitation energies, proceeds via a complex reaction involving rearrangement steps.

We propose a model for solvated positronium (Ps) atoms in water, based on the sequential quantum mechanics/molecular mechanics (s-QM/MM) protocol. We developed a Lennard-Jones force field to account for Ps–water interactions in the MM step. The repulsive term was obtained from a previously reported model for the solvated electron, while the dispersion constant was derived from the Slater–Kirkwood formula. The force field was employed in classical Monte Carlo (MC) simulations to generate Ps–solvent configurations in the NpTensemble, while the quantum properties were computed with the any-particle molecular orbital method in the subsequent QM step. Our approach is general, as it can be applied to other liquids and materials. One basically needs to describe the solvated electron in the environment of interest to obtain the Ps solvation model. The thermodynamical properties computed from the MC simulations point out similarities between the solvation of Ps and noble gas atoms, hydrophobic solutes that form clathrate structures. We performed convergence tests for the QM step, with particular attention to the choice of basis set and expansion centers for the positronic and electronic subsystems. Our largest model was composed of the Ps atom and 22 water molecules in the QM region, corresponding to the first solvation shell, surrounded by 128 molecules described as point charges. The mean electronic and positronic vertical detachment energies were (4.73 ± 0.04) eV and (5.33 ± 0.04) eV, respectively. The latter estimates were computed with Koopmans’ theorem corrected by second-order self-energies, for a set of statistically uncorrelated MC configurations. While the Hartree–Fock wave functions do not properly account for the annihilation rates, they were useful for numerical tests, pointing out that annihilation is more sensitive to the choice of basis sets and expansion centers than the detachment energies. We further explored a model with reduced solute cavity size by changing the Ps–solvent force field. Although the pick-off annihilation lifetimes were affected by the cavity size, essentially the same conclusions were drawn from both models.

When low-energy electrons interact with molecules, they can give rise to transient anion states commonly known as resonances. These states are formed through vertical electron attachment processes and have the potential to induce various forms of DNA lesions, including base damage, single- and double-strand breaks, cross-links, and clustered lesions that are challenging to repair. So far, most experimental and theoretical studies have investigated the formation of resonances of (bio)molecules in the gas phase or in microsolvated environments. Since cellular environments are mainly composed of water molecules, it is crucial to understand how bulk water affects the resonances of (bio)molecules. Given the existing gap in studies on resonances of bulk-solvated molecules, we propose a novel theoretical-computational approach to address this void. Our approach combines the multibasis-set (time-dependent-)density functional theory and self-consistent sequential quantum mechanics/molecular mechanics polarizable electrostatic embedding methods. We apply this combined methodology to predict the vertical electron attachment energies of 1-methyl-5-nitroimidazole (1M5NI), a well-known radiosensitizer model, in bulk water. In addition, we analyze the rapid mutual polarization between the resonances (both shape- and core-excited) of 1M5NI and the surrounding bulk water environment. For comparison, we also studied the isolated and microsolvated 1M5NI. Overall, while the polarization of the environment is clearly sensitive to the solute charge, causing a significant impact on the vertical electron affinity and consequently on the attachment electron energies, it does not have a significant impact on the excitation energies of the anion.

We report a computational study of the anion states of the resveratrol (RV) and resorcinol (RS) molecules, also investigating dissociative electron attachment (DEA) pathways. RV has well-known beneficial effects in human health, and its antioxidant activity was previously associated with DEA reactions producing H2. Our calculations indicate a valence bound state (π1*)and four resonances (π2*to π5*) for that system. While the computed thermodynamic thresholds are compatible with DEA reactions producing H2at 0 eV, the well-known mechanism involving vibrational Feshbach resonances built on a dipole bound state should not be present in RV. Our results suggest that the shallow π1*valence bound state is expected to account for H2elimination, probably involving π1*/σOH*couplings along the vibration dynamics. The RS molecule is also an oxidant and a subunit of RV. Because two close-lying hydroxyl groups are found in the RS moiety, the H2-elimination reaction in RV should take place at the RS site. Our calculations point out a correspondence between the anion states of RV and RS and even between the thresholds. Nevertheless, the absence of bound anion states in RS, indicated by our calculations, is expected to suppress the H2-formation channel at 0 eV. One is led to conclude that the ethene and phenol subunits in RV stabilize the π1*state, thus switching on the DEA mechanism producing H2.

We have characterized the size, intensity, density, and distribution of charge-transfer (CT) excitons as a function of the acceptor–donor architecture of prototypical organic interfaces. This characterization was done by computational analysis of 17 models of varying numbers, positions, and orientations of the donor and acceptor molecules. The models’ building blocks were phenyl-C61-butyric acid methyl ester (PCBM) fullerene acceptors and dual-band donor polymers composed of thiophene, benzothiadiazole, and benzotriazole subunits. The electronic structure of the donor–acceptor complexes was computed with the time-dependent long-range-corrected density-functional tight-binding method and analyzed with the fragment-based one-electron transition density matrix. In all models, the complexes with edge-on orientation have denser spectra of low-energy CT states lying below the absorption bands compared to the complexes with face-on orientation. This CT-state distribution in edge-on complexes provides a gate to efficiently populate cold CT excitons. Moreover, the cold CT excitons have a higher degree of charge separation in the edge-on than in the face-on complexes. The CT amount and the CT exciton size generally increase with the energy of the CT states, although the electron remains localized on a single molecule in cold CT states. Delocalization over two PCBM molecules was observed for high-energy CT states. The exciton size also depends on the orientation. Larger excitons are produced by the delocalization of the electrons perpendicularly to the interface. When the delocalization is parallel, the smaller electron–hole distances yield moderately sized CT excitons.