Successive attachment of electrons to protonated Guanine: (G+H)* radicals and (G+H)- anions

The structures, energetics, and vibrational frequencies of nine hydrogenated 9H-keto-guanine radicals (G+H)(*) and closed-shell anions (G+H)(-) are predicted using the carefully calibrated (Chem. Rev. 2002, 102, 231) B3LYP density functional method in conjunction with a DZP++ basis set. These radical and anionic species come from consecutive electron attachment to the corresponding protonated (G+H)(+) cations in low pH environments. The (G+H)(+) cations are studied using the same level of theory. The proton affinity (PA) of guanine computed in this research (228.1 kcal/mol) is within 0.7 kcal/mol of the latest experiment value. The radicals range over 41 kcal/mol in relative energy, with radical r1, in which H is attached at the C8 site of guanine, having the lowest energy. The lowest energy anion is a2, derived by hydride ion attachment at the C2 site of guanine. No stable N2-site hydride should exist in the gas phase. Structure a9 was predicted to be dissociative in this research. The theoretical adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies were computed, with AEAs ranging from 0.07 to 3.12 eV for the nine radicals.     

Pi and sigma-phenylethynyl radicals and their isomers o-, m-, and p-ethynylphenyl: structures, energetics, and electron affinities

Molecular structures, energetics, vibrational frequencies, and electron affinities are predicted for the phenylethynyl radical and its isomers. Electron affinities are computed using density functional theory, -namely, the BHLYP, BLYP, B3LYP, BP86, BPW91, and B3PW91 functionals-, employing the double-zeta plus polarization DZP++ basis set; this level of theory is known to perform well for the computation of electron affinities. Furthermore, ab initio computations employing perturbation theory, coupled cluster with single and double excitations [CCSD], and the inclusion of perturbative triples [CCSD(T)] are performed to determine the relative energies of the isomers. These higher level computations are performed with the correlation consistent family of basis sets cc-pVXZ (X = D, T, Q, 5). Three electronic states are probed for the phenylethynyl radical. In C2v symmetry, the out-of-plane (2B1) radical is predicted to lie about 10 kcal/mol below the in-plane (2B2) radical by DFT methods, which becomes 9.4 kcal/mol with the consideration of the CCSD(T) method. The energy difference between the lowest pi and sigma electronic states of the phenylethynyl radical is also about 10 kcal/mol according to DFT; however, CCSD(T) with the cc-pVQZ basis set shows this energy separation to be just 1.8 kcal/mol. The theoretical electron affinities of the phenylethynyl radical are predicted to be 3.00 eV (B3LYP/DZP++) and 3.03 eV (CCSD(T)/DZP++//MP2/DZP++). The adiabatic electron affinities (EAad) of the three isomers of phenylethynyl, that is, the ortho-, meta-, and para-ethynylphenyl, are predicted to be 1.45, 1.40, and 1.43 eV, respectively. Hence, the phenylethynyl radical binds an electron far more effectively than the three other radicals studied. Thermochemical predictions, such as the bond dissociation energies of the aromatic and ethynyl C-H bonds and the proton affinities of the phenylethynyl and ethynylphenyl anions, are also reported.     

Base-pair interactions in the gas-phase proton-bonded complexes of C+G and C+GC

Interactions involved in the formation of gas-phase proton-bonded molecular complexes of cytosine (C) and guanine (G) were theoretically investigated for the case of C(+)G and C(+)GC using B3LYP density functional theory. In this study, particular focus was on the dimeric interaction of proton-bonded C(+)G, where a proton bond and a hydrogen bond are cooperatively involved. The dimer interaction energy in terms of dissociation energy (D(e)) was predicted to be 41.8 kcal/mol. The lowest (frozen) energy structure for the C(+)G dimeric complex was found to be CH(+)...G rather than C...H(+)G in spite of the lower proton affinity of the cytosine moiety, which was more stable by 3.3 kcal/mol. The predicted harmonic vibrational frequencies and bond lengths suggest that the combined contributions of proton and hydrogen bonding may determine the resultant stability of each complex structure. In contrast to the dimer case, in the case of the isolated C(+)GC triplet, the two minimum energy structures of CH(+)...GC and C...H(+)GC were predicted to be almost equivalent in total energy. The dissociation energy (D(e)) for the C(+)G pairing in the C(+)GC triplet was 43.7 kcal/mol. Other energetics are also reported. As for the proton-transfer reaction in the proton-bond axis, the forward proton-transfer barriers for the dimer and trimer complexes were also predicted to be very low, 3.6 and 1.5 kcal/mol (DeltaE(e)(PT)), respectively.     

Microhydration of the guanine-cytosine (GC) base pair in the neutral and anionic radical states: a density functional study

A density functional study of the effects of microhydration on the guanine-cytosine (GC) base pair and its anion radical is presented. Geometries of the GC base pair in the presence of 6 and 11 water molecules were fully optimized in the neutral (GC-nH2O) and anion radical [(GC-nH2O)*-] (n = 6 and 11) states using the B3LYP method and the 6-31+G** basis set. Further, vibrational frequency analysis at the same level of theory (B3LYP/6-31+G**) was also performed to ensure the existence of local minima in these hydrated structures. It was found that water molecules surrounding the GC base pair have significant effects on the geometry of the GC base pair and promote nonplanarity in the GC base pair. The calculated structures were found to be in good agreement with those observed experimentally and obtained in molecular dynamics (MD) simulation studies. The water molecules in neutral GC-nH2O complexes lie near the ring plane of the GC base pair where they undergo hydrogen bonding with both GC and each other. However, in the GC anion radical complexes (GC-nH2O, n = 6, 11), the water molecules are displaced substantially from the GC ring plane. For GC-11H2O*-, a water molecule is hydrogen-bonded with the C6 atom of the cytosine base. We found that the hydration shell initially destabilizes the GC base pair toward electron capture as a transient anion. Energetically unstable diffuse states in the hydration shell are suggested to provide an intermediate state for the excess electron before molecular reorganization of the water molecules and the base pair results in a stable anion formation. The singly occupied molecular orbital (SOMO) in the anion radical complexes clearly shows that an excess electron localizes into a pi orbital of cytosine. The zero-point-energy (ZPE-) corrected adiabatic electron affinities (AEAs) of the GC-6H2O and GC-11H2O complexes, at the B3LYP/6-31+G** level of theory, were found to be 0.74 and 0.95 eV, respectively. However, the incorporation of bulk water as a solvent using the polarized continuum model (PCM) increases the EAs of these complexes to 1.77 eV.     

Ab initio and DFT benchmark study for nucleophilic substitution at carbon (SN2@C) and silicon (SN2@Si)

To obtain a set of consistent benchmark potential energy surfaces (PES) for the two archetypal nucleophilic substitution reactions of the chloride anion at carbon in chloromethane (S(N)2@C) and at silicon in chlorosilane (S(N)2@Si), we have explored these PESes using a hierarchical series of ab initio methods [HF, MP2, MP4SDQ, CCSD, CCSD(T)] in combination with a hierarchical series of six Gaussian-type basis sets, up to g polarization. Relative energies of stationary points are converged to within 0.01 to 0.56 kcal/mol as a function of the basis-set size. Our best estimate, at CCSD(T)/aug-cc-pVQZ, for the relative energies of the [Cl(-), CH(3)Cl] reactant complex, the [Cl-CH(3)-Cl](-) transition state and the stable [Cl-SiH(3)-Cl](-) transition complex is -10.42, +2.52, and -27.10 kcal/mol, respectively. Furthermore, we have investigated the performance for these reactions of four popular density functionals, namely, BP86, BLYP, B3LYP, and OLYP, in combination with a large doubly polarized Slater-type basis set of triple-zeta quality (TZ2P). Best overall agreement with our CCSD(T)/aug-cc-pVQZ benchmark is obtained with OLYP and B3LYP. However, OLYP performs better for the S(N)2@C overall and central barriers, which it underestimates by 2.65 and 4.05 kcal/mol, respectively. The other DFT approaches underestimate these barriers by some 4.8 (B3LYP) to 9.0 kcal/mol (BLYP).     

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.     

Proton affinities of maingroup-element hydrides and noble gases: trends across the periodic table, structural effects, and DFT validation

We have carried out an extensive exploration of the gas-phase basicity of archetypal neutral bases across the periodic system using the generalized gradient approximation (GGA) of the density functional theory (DFT) at BP86/QZ4P//BP86/TZ2P. First, we validate DFT as a reliable tool for computing proton affinities and related thermochemical quantities: BP86/QZ4P//BP86/TZ2P is shown to yield a mean absolute deviation of 2.0 kcal/mol for the proton affinity at 298 K with respect to experiment, and 1.2 kcal/mol with high-level ab initio benchmark data. The main purpose of this work is to provide the proton affinities (and corresponding entropies) at 298 K of the neutral bases constituted by all maingroup-element hydrides of groups 15-17 and the noble gases, that is, group 18, and periods 1-6. We have also studied the effect of step-wise methylation of the protophilic center of the second- and third-period bases.     

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.     

Theoretical prediction of the hydride affinities of various p- and o-quinones in DMSO

The hydride affinities of 80 various p- and o-quinones in DMSO solution were predicted by using B3LYP/6-311++G (2df,p)//B3LYP/6-31+G* and MP2/6-311++G**//B3LYP/6-31+G* methods, combined with the PCM cluster continuum model for the first time. The results show that the hydride affinity scale of the 80 quinones in DMSO ranges from -47.4 kcal/mol for 9,10-anthraquinone to -124.5 kcal/mol for 3,4,5,6-tetracyano-1,2-quinone. Such a long scale of the hydride affinities (-47.4 to -124.5 kcal/mol) indicates that the 80 quinones can form a large and useful library of organic oxidants, which can provide various organic hydride acceptors that the hydride affinities are known for chemists to choose in organic syntheses. By examining the effect of substituent on the hydride affinities of quinones, it is found that the hydride affinities of quinones in DMSO are linearly dependent on the sum of the Hammett substituent parameters sigma: DeltaGH-(Q) approximately -16.0Sigmasigmai - 70.5 (kcal/mol) for p-quinones and DeltaGH-(Q) approximately -16.2Sigmasigmai - 81.5 (kcal/mol) for o-quinones only if the substituents have no large electrostatic inductive effect and large ortho-effect. Study of the effect of the aromatic properties of quinone on the hydride affinities showed that the larger the aromatic system of quinone is, the smaller the hydride affinity of the quinone is, and the decrease of the hydride affinities is linearly to take place with the increase of the number of benzene rings in the molecule of quinones, from which the hydride affinities of aromatic quinones with multiple benzene rings can be predicted. By comparing the hydride affinities of p-quinones and the corresponding o-quinones, it is found that the hydride affinities of o-quinones are generally larger than those of the corresponding p-quinones by ca. 11 kcal/mol. Analyzing the effect of solvent on the hydride affinities of quinones showed that the effects of solvent (DMSO) on the hydride affinities of quinones are mainly dependent on the electrostatic interaction of the charged hydroquinone anions (QH-) with solvent (DMSO). All the information disclosed in this work should provide some valuable clues to chemists to choose suitable quinones or hydroquinones as efficient hydride acceptors or donors in organic syntheses and to predict the thermodynamics of hydride exchange between quinones and hydroquinones in DMSO solution.     

The ground and excited state electron affinities of cytosine and trans-azobenzene

The electron capture detector, reduction potential, electron transfer and photon methods of determining electron affinities are compared. The adiabatic electron affinities are (in eV): t-azobenzene(O₂), 1.578(5); t-azobenzene, 1.378(5); cytosine, 1.043(5) from anion photoelectron spectra. The largest or ground state value for trans-azobenzene and an excited state electron affinity for cytosine, 0.70 eV are also determined by reduction potentials. Other excited state energies are (in eV): t-azobenzene, 0.328(5), 0.589(5), 0.690(5), 0.768(5), 0.954(5), 1.038(5), 1.150(5), 1.275(5) and cytosine, 0.089(5), 0.098(5), 0.198(5), 0.235(5). The cytosine values are consistent with electron transport and radiation damage and repair in DNA.     

Microhydration of cytosine and its radical anion: cytosine.(H2O)n (n=1-5)

Microhydration effects on cytosine and its radical anion have been investigated theoretically, by explicitly considering various structures of cytosine complexes with up to five water molecules. Each successive water molecule (through n=5) is bound by 7-10 kcal mol(-1) to the relevant cytosine complex. The hydration energies are uniformly higher for the analogous anion systems. While the predicted vertical detachment energy (VDE) of the isolated cytosine is only 0.48 eV, it is predicted to increase to 1.27 eV for the lowest-lying pentahydrate of cytosine. The adiabatic electron affinity (AEA) of cytosine was also found to increase from 0.03 to 0.61 eV for the pentahydrate, implying that the cytosine anion, while questionable in the gas phase, is bound in aqueous solution. Both the VDE and AEA values for cytosine are smaller than those of uracil and thymine for a given hydration number. These results are in qualitative agreement with available experimental results from photodetachment-photoelectron spectroscopy studies of Schiedt et al. [Chem. Phys. 239, 511 (1998)].     

The proton affinity of proline analogs using the kinetic method with full entropy analysis

The proton affinity of proline analogs, L-azetidine-2-carboxylic acid (Aze), L-proline (Pro), and L-pipecolic acid (Pip), have been measured using the Armentrout modification of the extended kinetic method in a quadrupole ion trap instrument. Experimental values of 223.0 +/- 1.5, 224.9 +/- 1.6, and 225.6 +/- 1.6 kcal/mol have been determined for the 298K proton affinities of Aze, Pro, and Pip respectively. High level theoretical calculations using both MP2 and B3LYP methods at a variety of basis sets were carried out in order to give theoretical predictions for the 298 K proton affinity and gas phase basicity of all three analogs. Recommended values for the gas phase basicity and proton affinity for proline based on our work and other recent determinations are 216 +/- 2 and 224 +/- 2 kcal/mol.     

Hydrogen atom and hydride anion addition to adenine: structures and energetics.

The radicals and anions derived from the 9H tautomer of adenine by adding a hydrogen atom to one of the four double bonds of the adenine framework have been studied. Computations were carried out using a carefully calibrated density functional (B3LYP) method and basis set (DZP++). Optimized geometries, energies, and vibrational frequencies are predicted for eight radicals and anions. The radicals are found to lie in a range of 22 kcal mol(-1), with the radical derived by addition to the C(8) carbon atom being the lowest lying energetically. The anions are predicted to be bound species in the gas phase with an energetic range of 43 kcal mol(-1). Anions produced by addition of a hydride ion to adenine carbon atoms are found to be the most favorable. Six of the anions are predicted to be stable species with respect to electron detachment. The adiabatic electron affinities, vertical electron affinities, and vertical detachment energies are computed for the first time. Electron affinities for these radicals range from 0.0 to 2.0 eV. Radicals produced by addition to a nitrogen atom have near-zero adiabatic electron affinities, while radicals produced by addition at carbon atoms have considerably higher electron affinities.     

全氟代金刚烷及其自由基的理论研究

采用密度泛函理论(DFT)方法:BHLYP、B3LYP、BP86、BLYP,在全电子的双ζ基组加极化函数和弥散函数(DZP＋)基组下,计算了全氟代金刚烷(C10F16)及其自由基(C10F15)的总能量、优化几何构型、电子亲和势和谐振频率.在B3LYP水平上所得到的可靠绝热电子亲和势(EAad)分别为: C10F16, 1.06 eV; C10F15, 4.11和 3.03 eV.     

Effects of fluorine on the structures and energetics of the propynyl and propargyl radicals and their anions

[reaction: see text] The adiabatic electron affinity (EA(ad)) of the CH(3)-C[triple bond]C(*) radical [experiment = 2.718 +/- 0.008 eV] and the gas-phase basicity of the CH(3)-C[triple bond]C:(-) anion [experiment = 373.4 +/- 2 kcal/mol] have been compared with those of their fluorine derivatives. The latter are studied using theoretical methods. It is found that there are large effects on the electron affinities and gas-phase basicities as the H atoms of the alpha-CH(3) group in the propynyl system are substituted by F atoms. The predicted electron affinities are 3.31 eV (FCH(2)-C[triple bond]C(*)), 3.86 eV (F(2)CH-C[triple bond]C(*)), and 4.24 eV (F(3)C-C[triple bond]C(*)), and the predicted gas-phase basicities of the fluorocarbanion derivatives are 366.4 kcal/mol (FCH(2)-C[triple bond]C:(-)), 356.6 kcal/mol (F(2)CH-C[triple bond]C:(-)), and 349.8 kcal/mol (F(3)C-C[triple bond]C:(-)). It is concluded that the electron affinities of fluoropropynyl radicals increase and the gas-phase basicities decrease as F atoms sequentially replace H atoms of the alpha-CH(3) in the propynyl system. The propargyl radicals, lower in energy than the isomeric propynyl radicals, are also examined and their electron affinities are predicted to be 0.98 eV ((*)CH(2)-C[triple bond]CH), 1.18 eV ((*)CFH-C[triple bond]CH), 1.32 eV ((*)CF(2)-C[triple bond] CH), 1.71 eV ((*)CH(2)-C[triple bond]CF), 2.05 eV ((*)CFH-C[triple bond]CF), and 2.23 eV ((*)CF(2)-C[triple bond]CF).     

Structures and heats of formation of the neutral and ionic PNO, NOP, and NPO systems from electronic structure calculations

High level ab initio electronic structure calculations using the coupled cluster CCSD(T) method with augmented correlation-consistent basis sets extrapolated to the complete basis set limit have been performed on the PNO, NOP, and NPO isomers and their corresponding anions and cations. Geometries for all species were optimized up through the aug-cc-pV(Q+d)Z level and vibrational frequencies were calculated with the aug-cc-pV(T+d)Z basis set. The most stable of the three isomers is NPO and it is predicted to have a heat of formation of 23.3 kcal/mol. PNO is predicted to be only 1.7 kcal/mol higher in energy. The calculated adiabatic ionization potential of NPO is 12.07 eV and the calculated adiabatic electron affinity is 2.34 eV. The calculated adiabatic ionization potential of PNO is 10.27 eV and the calculated adiabatic electron affinity is only 0.24 eV. NOP is predicted to be much higher in energy by 29.9 kcal/mol. The calculated rotational constants for PNO and NPO should allow for these species to be spectroscopically distinguished. The adiabatic bond dissociation energies for the P[Single Bond]N, P[Single Bond]O, and N[Single Bond]O bonds in NPO and PNO are the same within approximately 10 kcal/mol and fall in the range of 72-83 kcal/mol.     

Temporary anion states and dissociative electron attachment in diphenyl disulfide

The temporary anion states of gas-phase diphenyl disulfide are characterized by means of electron transmission (ET) and dissociative electron attachment (DEA) spectroscopies. The measured energies of vertical electron attachment are compared to the virtual orbital energies of the neutral state molecule supplied by MP2 and B3LYP calculations with the 6-31G basis set. The calculated energies, scaled with empirical equations, reproduce satisfactorily the attachment energies measured in the ET spectrum. The first anion state of diphenyl disulfide is stable, thus escaping detection in ETS. The vertical and adiabatic electron affinities, evaluated with B3LYP/6-31+G calculations as the energy difference between the neutral and anion states, are predicted to be 0.37 and 1.38 eV, respectively. The anion current displayed in the DEA spectrum has a sharp and intense peak at zero energy, essentially due to the C6H5S- negative fragment. In agreement, according to the calculations, the localization properties of the first anion state are strongly S-S antibonding, and the energetic requirement for its dissociation along the S-S bond is fulfilled even at zero energy.     

Enthalpies of formation of gas-phase N3, N3-, N5+, and N5- from Ab initio molecular orbital theory, stability predictions for N5(+)N3(-) and N5(+)N5(-), and experimental evidence for the instability of N5(+)N3(-)

Ab initio molecular orbital theory has been used to calculate accurate enthalpies of formation and adiabatic electron affinities or ionization potentials for N3, N3-, N5+, and N5- from total atomization energies. The calculated heats of formation of the gas-phase molecules/ions at 0 K are DeltaHf(N3(2Pi)) = 109.2, DeltaHf(N3-(1sigma+)) = 47.4, DeltaHf(N5-(1A1')) = 62.3, and DeltaHf(N5+(1A1)) = 353.3 kcal/mol with an estimated error bar of +/-1 kcal/mol. For comparison purposes, the error in the calculated bond energy for N2 is 0.72 kcal/mol. Born-Haber cycle calculations, using estimated lattice energies and the adiabatic ionization potentials of the anions and electron affinities of the cations, enable reliable stability predictions for the hypothetical N5(+)N3(-) and N5(+)N5(-) salts. The calculations show that neither salt can be stabilized and that both should decompose spontaneously into N3 radicals and N2. This conclusion was experimentally confirmed for the N5(+)N3(-) salt by low-temperature metathetical reactions between N5SbF6 and alkali metal azides in different solvents, resulting in violent reactions with spontaneous nitrogen evolution. It is emphasized that one needs to use adiabatic ionization potentials and electron affinities instead of vertical potentials and affinities for salt stability predictions when the formed radicals are not vibrationally stable. This is the case for the N5 radicals where the energy difference between vertical and adiabatic potentials amounts to about 100 kcal/mol per N5.     

F+ and F⁻ affinities of simple N(x)F(y) and O(x)F(y) compounds

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the neutral and ionic N(x)F(y) and O(x)F(y) systems using coupled cluster theory with single and double excitations and including a perturbative triples correction (CCSD(T)) method with correlation consistent basis sets extrapolated to the complete basis set (CBS) limit. To achieve near chemical accuracy (±1 kcal/mol), three corrections to the electronic energy were added to the frozen core CCSD(T)/CBS binding energies: corrections for core-valence, scalar relativistic, and first order atomic spin-orbit effects. Vibrational zero point energies were computed at the CCSD(T) level of theory where possible. The calculated heats of formation are in good agreement with the available experimental values, except for FOOF because of the neglect of higher order correlation corrections. The F(+) affinity in the N(x)F(y) series increases from N(2) to N(2)F(4) by 63 kcal/mol, while that in the O(2)F(y) series decreases by 18 kcal/mol from O(2) to O(2)F(2). Neither N(2) nor N(2)F(4) is predicted to bind F(-), and N(2)F(2) is a very weak Lewis acid with an F(-) affinity of about 10 kcal/mol for either the cis or trans isomer. The low F(-) affinities of the nitrogen fluorides explain why, in spite of the fact that many stable nitrogen fluoride cations are known, no nitrogen fluoride anions have been isolated so far. For example, the F(-) affinity of NF is predicted to be only 12.5 kcal/mol which explains the numerous experimental failures to prepare NF(2)(-) salts from the well-known strong acid HNF(2). The F(-) affinity of O(2) is predicted to have a small positive value and increases for O(2)F(2) by 23 kcal/mol, indicating that the O(2)F(3)(-) anion might be marginally stable at subambient temperatures. The calculated adiabatic ionization potentials and electron affinities are in good agreement with experiment considering that many of the experimental values are for vertical processes.     

密度泛函研究细梗胡枝子提取物中新型抗氧化剂的分子结构和活性氧清除机理

The molecular structure and radical scavenging activity of three novel antioxidants from Lespedeza Virgata, lespedezavirgatol, lespedezavirgatal, and lespedezacoumestan, have been studied using density functional theory with the B3LYP and BhandHLYP methods. The optimized geometries of neutral, radical cation, radical and anion forms were obtained at the B3LYP/6-31G(d) level, in which it was found that all the most stable conformations contain intramolecular hydrogen bonds. The same results were obtained from the MP2 method. The homolytic O-H bond dissociation enthalpy and the adiabatic ionization potential of neutral and anion forms for the three new antioxidants and adiabatic electron affinity and H-atom affinity for hydroxyl radical, superoxide anion radical, and hydrogen peroxide radical were determined both in gas phase and in aqueous solution using IEF-PCM and CPCM model with UAHF or Bondi cavity. The antioxidant activities and reactive oxygen species scavenging mechanisms were then discussed, and the results obtained from different methods are consistent. Furthermore, the antioxidant activities are consistent with the experimental findings of the compounds under investigation.