The Adenine−Thymine Base Pair Radical Anion: Adding an Electron Results in a Major Structural Change

Electron Attachment to Wobble Base Pairs

We have analyzed the low-energy electron attachment to wobble base pairs using the equation of motion coupled cluster method and extended basis sets. A doorway mechanism exists for the attachment of the additional electron to the base pairs, where the initially formed dipole-bound anion captures the incoming electron. The doorway dipole-bound anionic state subsequently leads to the formation of a valence-bound state, and the transfer of extra electron occurs by mixing of electronic and nuclear degrees of freedom. The formation of the valence-bound anion is associated with proton transfer in hypoxanthine-cytosine and hypoxanthine-adenine base pairs, which happens through a concerted electron-proton transfer process. The calculated rate constant for the dipole-bound to valence-bound transition in wobble base pairs is slower than that observed in the Watson-Crick guanine-cytosine base pair.

Efficient implementation of the single-reference algebraic diagrammatic construction theory for charged excitations: Applications to the TEMPO radical and DNA base pairs

We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction (ADC) theory for electron attachment and ionization energies and spectra [EA/IP-ADC(n), n = 2, 3]. Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, and vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in good agreement with the best available results from the experiment and theoretical simulations.

Crystallographic evidence of Watson-Crick connectivity in the base pair of anionic adenine with thymine

Utilizing an ionic liquid strategy, we report crystal structures of salts of free anionic nucleobases and base pairs previously studied only computationally and in the gas phase. Reaction of tetrabutylammonium ([N₄₄₄₄]⁺) or tetrabutylphosphonium ([P₄₄₄₄]⁺) hydroxide with adenine (HAd) and thymine (HThy) led to hydrated salts of deprotonated adenine, [N₄₄₄₄][Ad]·2H₂O, and thymine, [P₄₄₄₄][Thy]·2H₂O, as well as the double salt cocrystal, [P₄₄₄₄]₂[Ad][Thy]·3H₂O·2HThy. The cocrystal includes the anionic [Ad^-(HThy)] base pair which is a stable formation in the solid state that has previously not even been suggested. It exhibits Watson-Crick connectivity as found in DNA but which is unusual for the free neutral base pairs. The stability of the observed anionic bases and their supramolecular formations and hydrates has also been examined by electronic structure calculations, contributing to more insight into how base pairs can bind when a proton is removed and highlighting mechanisms of stabilization or chemical transformation in the DNA chains.

Hydrogen-bonded proton transfer in the hydrated adenine–thymine anion

The proton transfer processes of microhydrated adenine–thymine anions are studied using density functional theory with the B3LYP method and DZP++ basis set. The microhydration effects on the geometrical structures, adsorption site, and the proton transfer reaction of the adenine–thymine anion are investigated. The site N10 atom of the adenine moiety has a larger proton affinity than the site O24 atom of thymine, which facilitates the proton H26 transfers from the N25 site of thymine to the N10 site of adenine. Therefore, the first single-proton transfer pathway (SPT1) is found for the all of the monohydrated adenine–thymine anions (AN4T)−·H2O, (AN13T)−·H2O, (ATO24)−·H2O, and (ATO28)−·H2O and tetrahydrated adenine–thymine anions (AT)−·4H2O. The proton H9 at the N7 site of adenine is also found to transfer to the O24 site of thymine for (AN4T)−·H2O and (AN13T)−·H2O in the gas phase. The double-proton transferred pathway is found when one water molecule interacts with the O28 atom of thymine. The reactant structures before the proton transfer are more stable than the product structures, and the structural changes mainly occur in thymine. The reaction energies of the microhydrated adenine–thymine anion in the gas phase and in the aqueous environment predict that the proton transfer process of the microhydrated adenine–thymine anion are more favorable in the gas phase than in aqueous solution.

Photoelectron spectroscopy and density functional theory studies on the uridine homodimer radical anions

We report the photoelectron spectrum (PES) of the homogeneous dimer anion radical of uridine, (rU)(2)(●-). It features a broad band consisting of an onset of ∼1.2 eV and a maximum at the electron binding energy (EBE) ranging from 2.0 to 2.5 eV. Calculations performed at the B3LYP∕6-31++G∗∗ level of theory suggest that the PES is dominated by dimeric radical anions in which one uridine nucleoside, hosting the excess charge on the base moiety, forms hydrogen bonds via its O8 atom with hydroxyl of the other neutral nucleoside's ribose. The calculated adiabatic electron affinities (AEAGs) and vertical detachment energies (VDEs) of the most stable homodimers show an excellent agreement with the experimental values. The anionic complexes consisting of two intermolecular uracil-uracil hydrogen bonds appeared to be substantially less stable than the uracil-ribose dimers. Despite the fact that uracil-uracil anionic homodimers are additionally stabilized by barrier-free electron-induced proton transfer, their relative thermodynamic stabilities and the calculated VDEs suggest that they do not contribute to the experimental PES spectrum of (rU)(2)(●-).

Electron attachment to the guanine-cytosine nucleic acid base pair and the effects of monohydration and proton transfer

The guanine-cytosine (GC) radical anion and its interaction with a single water molecule is studied using ab initio and density functional methods. Z-averaged second-order perturbation theory (ZAPT2) was applied to GC radical anion for the first time. Predicted spin densities show that the radical character is localized on cytosine. The Watson-Crick monohydrated GC anion is compared to neutral GC·H2O, as well as to the proton-transferred analogue on the basis of structural and energetic properties. In all three systems, local minima are identified that correspond to water positioned in the major and minor grooves of macromolecular DNA. On the anionic surface, two novel structures have water positioned above or below the GC plane. On the neutral and anionic surfaces, the global minimum can be described as water interacting with the minor groove. These structures are predicted to have hydration energies of 9.7 and 11.8 kcal mol(-1), respectively. Upon interbase proton-transfer (PT), the anionic global minimum has water positioned in the major groove, and the hydration energy increases to 13.4 kcal mol(-1). PT GC·H2O(•-) has distonic character; the radical character resides on cytosine, while the negative charge is localized on guanine. The effects of proton transfer are further investigated through the computed adiabatic electron affinities (AEA) of GC and monohydrated GC, and the vertical detachment energies (VDE) of the corresponding anions. Monohydration increases the AEAs and VDEs by only 0.1 eV, while proton-transfer increases the VDEs substantially (0.8 eV). The molecular charge distribution of monohydrated guanine-cytosine radical anion depends heavily on interbase proton transfer.

Electron attachment to solvated dGpdG: effects of stacking on base-centered and phosphate-centered valence-bound radical anions

To explore the nature of electron attachment to guanine-centered DNA single strands in the presence of a polarizable medium, a theoretical investigation of the DNA oligomer dinucleoside phosphate deoxyguanylyl-3',5'-deoxyguanosine (dGpdG) was performed by using density functional theory. Four different electron-distribution patterns for the radical anions of dGpdG in aqueous solution have been located as local minima on the potential energy surface. The excess electron is found to reside on the proton of the phosphate group (dGp(H-)dG), or on the phosphate group (dGp(.-)dG), or on the nucleobase at the 5' position (dG(.-)pdG), or on the nucleobase at the 3' position (dGpdG(.-)), respectively. These four radical anions are all expected to be electronically viable species under the influence of the polarizable medium. The predicted energetics of the radical anions follows the order dGp(.-)dG>dG(.-)pdG>dGpdG(.-)>dGp(H-)dG. The base-base stacking pattern in DNA single strands seems unaffected by electron attachment. On the contrary, intrastrand H-bonding is greatly influenced by electron attachment, especially in the formation of base-centered radical anions. The intrastrand H-bonding patterns revealed in this study also suggest that intrastrand proton transfer might be possible between successive guanines due to electron attachment to DNA single strands.

Isoguanine formation from adenine

Several possible mechanisms underlying isoguanine formation when OH radical attacks the C(2) position of adenine (A C 2) are investigated theoretically for the first time. Two steps are involved in this process. In the first step, one of two low-lying A C 2⋅⋅⋅OH reactant complexes is formed, leading to C(2)-H(2) bond cleavage. Between the two reactant complexes there is a small isomerization barrier, which lies well below separated adenine plus OH radical. The complex dissociates to free molecular hydrogen and an isoguanine tautomer (isoG 1 or isoG 2). The local and activation barriers for the two pathways are very similar. This evidence suggests that the two pathways are competitive. After dehydrogenation, there are two possible routes for the second step of the reaction. One is direct hydrogen transfer, via enol-keto tautomerization, which has high local barriers for both tautomers and is not favored. The other option is indirect hydrogen transfer involving microsolvation by one water molecule. The water lowers the reaction barrier by over 20 kcal  mol(-1) , indicating that water-mediated hydrogen transfer is much more favorable. Both A+OH(⋅) →isoG+H(⋅) reactions are exothermic and spontaneous. Among four isoguanine tautomers, isoG 1 has the lowest energy. Our findings explain why only the N(1)H and O(2)H tautomers of isolated isoguanine and isoguanosine have been observed experimentally.

Mutagenicity associated with O6-methylguanine-DNA damage and mechanism of nucleotide flipping by AGT during repair

Methylated guanine damage at O6 position (i.e. O6MG) is dangerous due to its mutagenic and carcinogenic character that often gives rise to G:C-A:T mutation. However, the reason for this mutagenicity is not known precisely and has been a matter of controversy. Further, although it is known that O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6MG paired with cytosine in DNA, the complete mechanism of target recognition and repair is not known completely. All these aspects of DNA damage and repair have been addressed here by employing high level density functional theory in gas phase and aqueous medium. It is found that the actual cause of O6MG mediated mutation may arise due to the fact that DNA polymerases incorporate thymine opposite to O6MG, misreading the resulting O6MG:T complex as an A:T base pair due to their analogous binding energies and structural alignments. It is further revealed that AGT mediated nucleotide flipping occurs in two successive steps. The intercalation of the finger residue Arg128 into the DNA double helix and its interaction with the O6MG:C base pair followed by rotation of the O6MG nucleotide are found to be crucial for the damage recognition and nucleotide flipping.

Electron attachment to a hydrated DNA duplex: the dinucleoside phosphate deoxyguanylyl-3',5'-deoxycytidine

The minimal essential section of DNA helices, the dinucleoside phosphate deoxyguanylyl-3',5'-deoxycytidine dimer octahydrate, [dGpdC](2), has been constructed, fully optimized, and analyzed by using quantum chemical methods at the B3LYP/6-31+G(d,p) level of theory. Study of the electrons attached to [dGpdC](2) reveals that DNA double strands are capable of capturing low-energy electrons and forming electronically stable radical anions. The relatively large vertical electron affinity (VEA) predicted for [dGpdC](2) (0.38 eV) indicates that the cytosine bases are good electron captors in DNA double strands. The structure, charge distribution, and molecular orbital analysis for the fully optimized radical anion [dGpdC](2)(·-) suggest that the extra electron tends to be redistributed to one of the cytosine base moieties, in an electronically stable structure (with adiabatic electron affinity (AEA) 1.14 eV and vertical detachment energy (VDE) 2.20 eV). The structural features of the optimized radical anion [dGpdC](2)(·-) also suggest the probability of interstrand proton transfer. The interstrand proton transfer leads to a distonic radical anion [d(G-H)pdC:d(C+H)pdG](·-), which contains one deprotonated guanine anion and one protonated cytosine radical. This distonic radical anion is predicted to be more stable than [dGpdC](2)(·-). Therefore, experimental evidence for electron attachment to the DNA double helices should be related to [d(G-H)pdC:d(C+H)pdG](·-) complexes, for which the VDE might be as high as 2.7 eV (in dry conditions) to 3.3 eV (in fully hydrated conditions). Effects of the polarizable medium have been found to be important for increasing the electron capture ability of the dGpdC dimer. The ultimate AEA value for cytosine in DNA duplexes is predicted to be 2.03 eV in aqueous solution.

Electron attachment to hydrated oligonucleotide dimers: guanylyl-3',5'-cytidine and cytidylyl-3',5'-guanosine

The dinucleoside phosphate deoxycytidylyl-3',5'-deoxyguanosine (dCpdG) and deoxyguanylyl-3',5'-deoxycytidine (dGpdC) systems are among the largest to be studied by reliable theoretical methods. Exploring electron attachment to these subunits of DNA single strands provides significant progress toward definitive predictions of the electron affinities of DNA single strands. The adiabatic electron affinities of the oligonucleotides are found to be sequence dependent. Deoxycytidine (dC) on the 5' end, dCpdG, has larger adiabatic electron affinity (AEA, 0.90 eV) than dC on the 3' end of the oligomer (dGpdC, 0.66 eV). The geometric features, molecular orbital analyses, and charge distribution studies for the radical anions of the cytidine-containing oligonucleotides demonstrate that the excess electron in these anionic systems is dominantly located on the cytosine nucleobase moiety. The pi-stacking interaction between nucleobases G and C seems unlikely to improve the electron-capturing ability of the oligonucleotide dimers. The influence of the neighboring base on the electron-capturing ability of cytosine should be attributed to the intensified proton accepting-donating interaction between the bases. The present investigation demonstrates that the vertical detachment energies (VDEs) of the radical anions of the oligonucleotides dGpdC and dCpdG are significantly larger than those of the corresponding nucleotides. Consequently, reactions with low activation barriers, such as those for O-C sigma bond and N-glycosidic bond breakage, might be expected for the radical anions of the guanosine-cytosine mixed oligonucleotides.

Stability of the valence anion of cytosine is governed by nucleobases sequence in the double stranded DNA pi-stack: A computational study

The stabilities of the valence anion of cytosine (C(-)) in model trimers of complementary base pairs that possess the B-DNA geometry but differ in base sequence are reported. In order to estimate the energetics of electron attachment to the middle cytosine incorporated in the trimer, a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated cytosine were developed. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. We have demonstrated that contrary to the literature reports, concerning single stranded DNA, the sequence of nucleic bases has a profound effect on the stability of the cytosine valence anion. The anionic 3(')-CCC-5(') complex is the most stable configuration (EA=0.399 eV) and the 3(')-GCG-5(') trimer anion is the most unstable species (EA=-0.193 eV). Moreover, with the energetic correction for the presence of sugar-phosphate backbone all possible double stranded DNA sequences lead to the stable C(-). The predicted electron affinities of the cytosine anion have been compared to the results of analogous studies on the thymine anion published recently [M. Kobyłecka et al., J. Am. Chem. Soc. 130, 15683 (2008)]. The consequences of low-energy barrier proton transfer in the GC anion have been discussed in the context of induced by electrons DNA single strand breaks. The DNA sequences that should dramatically differ in their vulnerability to be damaged by low energy electrons have been proposed.

Stabilization of radical anion states of nucleobases in DNA

Trapping of an electron by DNA leads to the formation of radical anion states of pyrimidine bases. Because these states play an important role in biological and chemical processes, their computational treatment is of particular interest. We show that simple electrostatic and quantum chemical models can accurately reproduce the adiabatic electron affinities (EAs) of short DNA stacks recently derived from high-level ab initio calculations (M. Kobylecka, J. Leszczynski, and J. Rak, J. Am. Chem. Soc., 2008, 130, 15683). The electrostatic interaction of an excess electron localized on cytosine or thymine with intra- and inter-strand adjacent nucleobases is found to strongly affect the energy of the radical anions. This interaction is the main origin of the dependence of EA of nucleobases on the nature of neighboring base pairs. In particular, the states XT(-)Y and XC(-)Y, where X and Y = C, T, are, by ca. 0.7 eV, more stable than radical anions GT(-)G and GC(-)G. We find that second-neighbor effects can also significantly modulate EAs, although being smaller than the effects of adjacent bases. The strongest destabilizing effect is found for 5'-GC and 3'-GC, while the 5'-AT base pair stabilizes the radical anion states. Using a combined QM/MD approach, we consider how structural fluctuations of DNA influence the stability of the radical anion states. Despite large dispersions of the stabilization energies due to conformational dynamics of DNA, there are only few thermally accessible structures where GT(-)G and GC(-)G are energetically more favorable than the corresponding pyrimidine triplets. Although stabilization energies calculated for stacks of regular structure are in qualitative agreement with the QM/MD results, structural fluctuations of pi stacks should be taken into account for more accurate description of the excess electron trapped by DNA. The results obtained in this study suggest that simple electrostatic models, in combination with MD simulations, can be very helpful to explore the long time scale behavior of radical anions in DNA.

Elementary lesions in DNA subunits: electron, hydrogen atom, proton, and hydride transfers

When DNA is damaged by ionizing radiation, the genes in a cell may acquire mutations or the cell could die. The smallest known DNA-damaging unit is an electron, often low-energy secondary electrons. Additional electrons and transfers involving hydrogen atoms, protons, and hydride anions can damage DNA subunits, including individual nucleobases and nucleoside pairs. Researchers would like to better understand the molecular mechanisms involved in DNA damage from ionizing radiation. In this Account, we highlight our theoretical investigations of the molecular mechanisms of DNA damage using quantum mechanical models. Our investigations use robust theoretical methods with computations conducted in the gas phase and with solution models. We calculate adiabatic electron affinities (AEAs), which describe the energetics of electronic attachment to closed-shell DNA subunits, for the free bases, nucleosides, nucleotides, base pairs, and single and double DNA strand units. Electron affinities for free nucleobases yield the order uracil > thymine > cytosine > guanine > adenine and the same order for the DNA nucleosides, mononucleotides, and nucleoside 3',5'-diphosphates. AEA values increase steadily with the size and complexity of the system because of greater hydration, glycosylation, nucleotide formation, and base pairing. We predict and experimental results partially confirm that most of the more complex and hydrated species are observable as radical anions. Our modeling studies indicate that depyrimidination reactions of radical anion nucleosides release cytosine more often than thymine. Recent experimental results support those findings. Our theoretical studies of DNA base-pair radical anions predict increases in electron affinity accompanying H bonding and solvation. Electron addition facilitates some proton transfers in these pairs, which results in strongly perturbed pairing configurations. Of all nucleobase moieties within the more complex radical anion systems, thymine is best able to retain a negative charge. Charge and spin are well-separated in some of these systems. Radical species derived via hydrogen abstraction from DNA subunits yield large AEA values because they form closed-shell anions. Our studies predict single-strand breaks following H abstraction from nucleotides. Some H-abstraction processes in the DNA base pairs lead to severe distortions in pairing configuration based on our calculations. This body of systematic theoretical studies provides realistic descriptions of some events that lead to elementary DNA lesions, while providing rationalizations for many observed phenomena. Such approaches can lead to the design of new experiments, which would contribute to our understanding of the chemical physics of nucleic acids.

Ab initio determination of the electron affinities of DNA and RNA nucleobases

High-level quantum-chemical ab initio coupled-cluster and multiconfigurational perturbation methods have been used to compute the vertical and adiabatic electron affinities of the five canonical DNA and RNA nucleobases: uracil, thymine, cytosine, adenine, and guanine. The present results aim for the accurate determination of the intrinsic electron acceptor properties of the isolated nucleic acid bases as described by their electron affinities, establishing an overall set of theoretical reference values at a level not reported before and helping to rule out less reliable theoretical and experimental data and to calibrate theoretical strategies.

On the hydrogen loss from protonated nucleobases after electronic excitation or collisional electron capture

In this work, we have subjected protonated nucleobases MH(+) (M = guanine, adenine, thymine, uracil and cytosine) to a range of experiments that involve high-energy (50 keV) collision induced dissociation and electron capture induced dissociation. In the latter case, both neutralisation reionisation and charge reversal were done. For the collision induced dissociation experiments, the ions interacted with O(2). In neutral reionisation, caesium atoms were used as the target gas and the protonated nucleobases captured electrons to give neutrals. These were reionised to cations a microsecond later in collisions with O(2). In choosing Cs as the target gas, we have assured that the first electron transfer process is favourable (by about 0.1-0.8 eV depending on the base). In the case of protonated adenine, charge reversal experiments (two Cs collisions) were also carried out, with the results corroborating those from the neutralisation reionisation experiments. We find that while collisional excitation of protonated nucleobases in O(2) may lead to hydrogen loss with limited probabilities, this channel becomes dominant for electron capture events. Indeed, when sampling reionised neutrals on a microsecond timescale, we see that the ratio between MH(+) and M(+) is 0.2-0.4 when one electron is captured from Cs. There are differences in these ratios between the bases but no obvious correlation with recombination energies was found.

Valence anion of thymine in the DNA pi-stack

Most of theoretical data on the stability of radical anions supported by nucleic acid bases have been obtained for anions of isolated nucleobases, their nucleosides, or nucleotides. This approach ignores the hallmark forces of DNA, namely, hydrogen bonding and pi-stacking interactions. Since these interactions might be crucial for the electron affinities of nucleobases bound in DNA, we report for the first time on the stability of the thymine valence anion in trimers of complementary bases possessing the regular B-DNA geometry but differing in base sequence. In order to estimate the energetics of electron attachment to a trimer, we developed a thermodynamic cycle employing all possible two-body interaction energies in the neutral and anionic duplex as well as the adiabatic electron affinity of isolated thymine. All calculations were carried out at the MP2 level of theory with the aug-cc-pVDZ basis set. The two-body interaction energies were corrected for the basis set superposition error, and in benchmark systems, they were extrapolated to the basis set limit and supplemented with correction for higher order correlation terms calculated at the CCSD(T) level. We have demonstrated that the sequence of nucleic bases has a profound effect on the stability of the thymine valence anion: the anionic 5'-CTC-3' (6.0 kcal/mol) sequence is the most stable configuration, and the 5'-GTG-3' (-8.0 kcal/mol) trimer anion is the most unstable species. On the basis of obtained results, one can propose DNA sequences that are different in their vulnerability to damage by low energy electron.

Formation of spectral intermediate G-C and A-T anion complex in duplex DNA studied by pulse radiolysis

The dynamics of electron adducts of 2'-deoxynucleotides and oligonucelotides (ODNs) were measured spectroscopically by nanosecond pulse radiolysis. The radical anions of the nucleotides were produced within 10 ns by the reaction of hydrated electrons (e(aq)(-)) and were protonated to form the corresponding neutral radicals. At pH 7.0, the radical anion of deoxythymidine (dT(*-)) was protonated to form the neutral radical dT(H)(*) in the time range of microseconds. The rate constant for the protonation was determined as 1.8 x 10(10) M(-1) s(-1). In contrast, the neutral radical of dC(H)(*) was formed immediately after the pulse, suggesting that the protonation occurs within 10 ns. The transient spectra of excess electrons of the double-stranded ODNs 5'-TAATTTAATAT-3' (AT) and 5'-CGGCCCGGCGC-3' (GC) differed from those of pyrimidine radicals (C and T) and their composite. In contrast, the spectra of the electron adducts of the single-stranded ODNs GC and AT exhibited characteristics of C and T, respectively. These results suggest that, in duplex ODNs, the spectral intermediates of G-C and A-T anions complex were formed. On the microsecond time scale, the subsequent changes in absorbance of the ODN AT had a first-order rate constant of 4 x 10(4) s(-1), reflecting the protonation of T.

Electron attachment to the hydrogenated Watson-Crick guanine cytosine base pair (GC+H): conventional and proton-transferred structures

The anionic species resulting from hydride addition to the Watson-Crick guanine-cytosine (GC) DNA base pair are investigated theoretically. Proton-transferred structures of GC hydride, in which proton H1 of guanine or proton H4 of cytosine migrates to the complementary base-pair side, have been studied also. All optimized geometrical structures are confirmed to be minima via vibrational frequency analyses. The lowest energy structure places the additional hydride on the C6 position of cytosine coupled with proton transfer, resulting in the closed-shell anion designated 1T (G(-)C(C6)). Energetically, the major groove side of the GC pair has a greater propensity toward hydride/hydrogen addition than does the minor grove side. The pairing (dissociation) energy and electron-attracting ability of each anionic structure are predicted and compared with those of the neutral GC and the hydrogenated GC base pairs. Anion 8T (G(O6)C(-)) is a water-extracting complex and has the largest dissociation energy. Anion 2 (GC(C4)(-)) and the corresponding open-shell radical GC(C4) have the largest vertical electron detachment energy and adiabatic electron affinity, respectively. From the difference between the dissociation energy and electron-removal ability of the normal GC anion and the most favorable structure of GC hydride, it is clear that one may dissociate the GC anion and maintain the integrity of the GC hydride.

Species with negative electron affinity and standard DFT methods. Finding the valence anions

Recently, we have shown that traditional bound-electron DFT models are reliable enough to reproduce negative electron affinities (EA) within a few meV, as long as the valence anion state is found, but they seem to fail in predicting the lowest EA when the ground anion state obtained is non-valence, which holds the extra electron in a diffuse orbital around the molecule; here we propose an alternative approach for finding the valence anion state, based on the stabilization exerted by a polar solvent; the methodology yields correct EA values (i.e. beyond the Koopman's theorem approximation) by gradually decreasing the dielectric constant of the medium.