Base and phosphate electron detachment energies of deoxyribonucleotide anions

Photoelectron spectra of deoxyribonucleotide anions are interpreted with ab initio, electron propagator calculations. Ground-state structures display hydrogen bonds which are not present in less stable minima that resemble Watson-Crick fragment geometries. For the adenosine and thymidine anions, there are two vertical electron detachment energies (VEDEs) within 0.1 eV of each other that correspond to phosphate- and base-centered Dyson orbitals (DOs). The first VEDE of the cytidine anion belongs to a phosphate-centered DO. The anomalously low VEDE of the guanosine anion is assigned to a base-centered, pi DO. Higher VEDEs of all four anions also are assigned.     

Nonconventional hydrogen bonds: a theoretical study of [uracil-L](-) (L = F, Cl, Br, I, Al, Ga, In) complexes

The interaction of L (-) (L = F, Cl, Br, I, Al, Ga and In) with a uracil molecule has been studied with B3LYP density-functional geometry optimizations and electron-propagator calculations of vertical electron detachment energies. Because the extra electron of the anion is localized on L, nonconventional hydrogen bonds are formed. The interactions of halide anions with uracil are similar to the interactions of uracil with Cu (-), Ag (-) and Au (-) that were reported previously. Whereas halide and transition metal anion complexes with uracils are singlets, the anions formed with Al, Ga and In are triplets. Vertical electron detachment energies (VEDEs) are higher for (uracil-L) (-) than the analogous values for isolated L (-) anions. Predicted VEDEs are assigned to Dyson orbitals that may be localized on L (-) or uracil.     

Study of the valence wave function of thiophene with high resolution electron momentum spectroscopy and advanced Dyson orbital theories

Results of an exhaustive experimental study of the valence electronic structure of thiophene using high resolution electron momentum spectroscopy at impact energies of 1200 and 2400 eV are presented. The measurements were performed using an electron momentum spectrometer of the third generation at Tsinghua University, which enables energy, polar and azimuthal angular resolutions of the order of DeltaE = 0.8 eV, Deltatheta = +/-0.53 degrees and Deltaphi = +/-0.84 degrees . These measurements were interpreted by comparison with Green's function calculations of one-electron and shake-up ionization energies as well as of the related Dyson orbital electron momentum distributions, using the so-called third-order algebraic diagrammatic construction scheme (ADC(3)). Comparison of spherically averaged theoretical electron momentum distributions with experimental results very convincingly confirms the presence of two rather intense pi-2 pi*+1 shake-up lines at electron binding energies of 13.8 and 15.5 eV, with pole strengths equal to 0.18 and 0.13, respectively. Analysis of the electron momentum distributions associated with the two lowest 2A2 (pi3-1) and 2B1 (pi2-1) cationic states provides indirect evidence for a symmetry lowering and nuclear dynamical effects due to vibronic coupling interactions between these two states. ADC(3) Dyson orbital momentum distributions are systematically compared with distributions derived from Kohn-Sham (B3LYP) orbitals, and found to provide most generally superior insights into experiment.     

Probing molecular conformations in momentum space: the case of n-pentane

A comprehensive study, throughout the valence region, of the electronic structure and electron momentum density distributions of the four conformational isomers of n-pentane is presented. Theoretical (e,2e) valence ionization spectra at high electron impact energies (1200 eV+electron binding energy) and at azimuthal angles ranging from 0 degrees to 10 degrees in a noncoplanar symmetric kinematical setup are generated according to the results of large scale one-particle Green's function calculations of Dyson orbitals and related electron binding energies, using the third-order algebraic-diagrammatic construction [ADC(3)] scheme. The results of a focal point analysis (FPA) of relative conformer energies [A. Salam and M. S. Deleuze, J. Chem. Phys. 116, 1296 (2002)] and improved thermodynamical calculations accounting for hindered rotations are also employed in order to quantitatively evaluate the abundance of each conformer in the gas phase at room temperature and reliably predict the outcome of experiments on n-pentane employing high resolution electron momentum spectroscopy. Comparison with available photoelectron measurements confirms the suggestion that, due to entropy effects, the trans-gauche (tg) conformer strongly dominates the conformational mixture characterizing n-pentane at room temperature. Our simulations demonstrate therefore that experimental measurements of (e,2e) valence ionization spectra and electron momentum distributions would very consistently and straightforwardly image the topological changes and energy variations that molecular orbitals undergo due to torsion of the carbon backbone. The strongest fingerprints for the most stable conformer (tt) are found for the electron momentum distributions associated with ionization channels at the top of the inner-valence region, which sensitively image the development of methylenic hyperconjugation in all-staggered n-alkane chains.     

Addition of water, methanol, and ammonia to Al3O3- clusters: reaction products, transition states, and electron detachment energies

Products of reactions between the book and kite isomers of Al3O3- and three important molecules are studied with electronic structure calculations. Dissociative adsorption of H2O or CH3OH is highly exothermic and proton-transfer barriers between anion-molecule complexes and the products of these reactions are low. For NH3, the reaction energies are less exothermic and the corresponding barriers are higher. Depending on experimental conditions, Al3O3- (NH3) coordination complexes or products of dissociative adsorption may be prepared. Vertical electron detachment energies of stable anions are predicted with ab initio electron propagator calculations and are in close agreement with experiments on Al3O3- and its products with H2O and CH3OH. Changes in the localization properties of two Al-centered Dyson orbitals account for the differences between the photoelectron spectra of Al3O3- and those of the product anions.     

Interpretation of the photoelectron spectra of superalkali species: Na3O and Na3O-

Recently measured photoelectron spectra of the Na(3)O(-) anion have been interpreted with the aid of ab initio electron propagator calculations. As in the case of the Li(3)O(-), we propose that the photoionization of ground and excited neutral states, in a sequential two photon absorption mechanism, plays a role in the interpretation of the observed spectrum. The lowest vertical electron detachment energy of Na(3)O(-) corresponds to a Dyson orbital that is composed chiefly of diffuse Na s functions and connects a D(3h) singlet anion to an uncharged species with the same point group. Electron binding energies of isomers of the anion with different point groups or multiplicities have been considered. The relative magnitudes of the ionization energies of the neutral Li(3)O and Na(3)O species are also discussed. Whereas the most recent experimental data hold that Na(3)O has the higher ionization energy, this work asserts the opposite trend.     

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.     

Extra electron in (H2O)24- cluster isomers: a theoretical study

The isomers of (H(2)O)(24) (-) tetrakaidecahedral cluster are studied by applying the Becke-3-parameter density functional theory and Lee-Yang-Parr correlation functional (B3LYP) and 6-311++G** basis set. Three isomers are selected on the basis of stabilization energy values. The vertical electron dissociation energies (VDE) of these isomers are 1.353, 0.404, and 0.258 eV, respectively. The experimental VDE value of 1.31 eV [J. Chem. Phys. 92, 3980 (1990)] for this cluster size is in excellent agreement with that calculated for isomer 1, suggesting the dominance of this isomer in the experiment. Four water molecules in this isomer share most of the -1 charge. These four water molecules have non-H-bonding H (NHB H) atoms turned toward the cavity, and the inward turned H atoms exhibit a significant lowering of O-H stretch frequency compared to that of a monomer. Isomers 2 and 3 have all 12 NHB H atoms projected outward and have the -1 charge distributed among 7-8 water molecules on the cluster surface.     

MX3(-) superhalogens (M = Be, Mg, Ca; X = Cl, Br): a photoelectron spectroscopic and ab initio theoretical study

Gas-phase alkaline earth halide anions, MgX3(-) and CaX3(-) (X = Cl, Br), were produced using electrospray and investigated using photoelectron spectroscopy at 157 nm. Extremely high electron binding energies were observed for all species and their first vertical detachment energies were measured as 6.60 +/- 0.04 eV for MgCl3(-), 6.00 +/- 0.04 eV for MgBr3(-), 6.62 +/- 0.04 eV for CaCl3(-), and 6.10 +/- 0.04 eV for CaBr3(-). The high electron binding energies indicate these are very stable anions and they belong to a class of anions, called superhalogens. Theoretical calculations at several levels of theory were carried out on these species, as well as the analogous BeX3(-). Vertical detachment energy spectra were predicted to compare with the experimental observations, and good agreement was obtained for all species. The first adiabatic detachment energies were found to be substantially lower (by about 1 eV) than the corresponding vertical detachment energies for all the MX3(-) species, indicating extremely large geometry changes between MX3(-) and MX3. We found that all the MX3(-) anions possess D3h ((1)A1') structures and are extremely stable against dissociation into MX2 and X-. The corresponding neutral species MX3, however, were found to be only weakly bound with respect to dissociation toward MX2 + X. The global minimum structures of all the MX3 neutrals were found to be C2v ((2)B2), which can be described as (X2(-))(MX+) charge-transfer complexes, whereas the MX2...X (C2v, (2)B1) van der Waals complexes were shown to be low-lying isomers.     

OH(3) (-) and O(2)H(5) (-) double Rydberg anions: predictions and comparisons with NH(4) (-) and N(2)H(7) (-)

A low barrier in the reaction pathway between the double Rydberg isomer of OH(3) (-) and a hydride-water complex indicates that the former species is more difficult to isolate and characterize through anion photoelectron spectroscopy than the well known double Rydberg anion (DRA), tetrahedral NH(4) (-). Electron propagator calculations of vertical electron detachment energies (VEDEs) and isosurface plots of the electron localization function disclose that the transition state's electronic structure more closely resembles that of the DRA than that of the hydride-water complex. Possible stabilization of the OH(3) (-) DRA through hydrogen bonding or ion-dipole interactions is examined through calculations on O(2)H(5) (-) species. Three O(2)H(5) (-) minima with H(-)(H(2)O)(2), hydrogen-bridged, and DRA-molecule structures resemble previously discovered N(2)H(7) (-) species and have well separated VEDEs that may be observable in anion photoelectron spectra.     

Anion photoelectron spectroscopy of C(5)H(-)

Anion photoelectron spectroscopy is performed on the C(5)H(-) species. Analogous to C(3)H(-) and C(3)D(-), photodetachment transitions are observed from multiple, energetically close-lying isomers of the anion. A linear and a cyclic structure are found to have electron binding energies of 2.421+/-0.019 eV and 2.857+/-0.028 eV, respectively. A cyclic excited state is also found to be 1.136 eV above the linear (2)Pi C(5)H ground state. Based on our assignments of the observed transitions and previous calculations on the energetics of neutral C(5)H isomers, the cyclic (1)A(1) anion state is found to lie 0.163 eV below the (3)A linear anion.     

Photoabsorption spectra of Ti8C12 metallocarbohedrynes: Theoretical spectroscopy within time-dependent density functional theory

The photoabsorption spectra of several of the most stable isomers of the Ti8C12 metallocarbohedryne are calculated using time-dependent density functional theory. Several ground-state magnitudes have been also calculated, such as cohesive energies, electronic gaps between the highest occupied and lowest unoccupied molecular orbitals, and static polarizabilities. Since significant differences are found among the photoabsorption spectra of the different isomers in the low energy region (0-5 eV), we propose the comparison of experimental and the calculated absorption spectra as a tool to elucidate the isomers that appear to form in the experiments. Between 10 and 13 eV all the spectra show a region of high absorption that we identify as due to collective electronic excitations. The existence of this prominent feature explains the occurrence of delayed ionization and delayed ion emission phenomena observed in previous experiments.     

Cyclopentadiene annulated polycyclic aromatic hydrocarbons: investigations of electron affinities

The adiabatic electron affinities of cyclopentadiene and 10 associated benzannelated derivatives have been predicted with both density functional and Hartree-Fock theory. These systems can also be regarded as benzenoid polycyclic aromatic hydrocarbons (PAHs) augmented with five-membered rings. Like the PAHs, the electron affinities of the present systems generally increase with the number of rings. To unequivocally bind an electron, cyclopentadiene must have at least two conventionally fused benzene rings. 1H-Benz[f]indene, a naphthalene-annulated cyclopentadiene, is predicted to have a zero-point energy corrected adiabatic electron affinity of 0.13 eV. Since the experimental E(A) of naphthalene is negative (-0.19 eV), the five-membered ring appendage contributes to the stability of the naphthalene-derived 1H-benz[f]indene radical anion significantly. The key to binding the electron is a contiguous sequence of fused benzenes, since fluorene, the isomer of 1H-benz[f]indene, with separated six-membered rings, has an electron affinity of -0.07 eV. Each additional benzene ring in the sequence fused to cyclopentadiene increases the electron affinity by 0.15-0.65 eV: the most reliable predictions are cyclopentadiene (-0.63 eV), indene (-0.49 eV), fluorene (-0.07 eV), 1H-benz[f]indene (0.13 eV), 1,2-benzofluorene (0.25 eV), 2,3-benzofluorene (0.26 eV), 12H-dibenzo[b,h]fluorene (0.65 eV), 13H-indeno[1,2-b]anthracene (0.82 eV), and 1H-cyclopenta[b]naphthacene (1.10 eV). In contrast, if the six-membered ring-fusion is across the C(2)-C(3) cyclopentadiene single bond, only a single benzene is needed to bind an electron. The theoretical electron affinity of the resulting molecule, isoindene, is 0.49 eV, and this increases to 1.22 eV for 2H-benz[f]indene. The degree of aromaticity is responsible for this behavior. While the radical anions are stabilized by conjugation, which increases with the size of the system, the regular indenes, like PAHs in general, suffer from the loss of aromatic stabilization in forming their radical anions. While indene is 21 kcal mol(-1) more stable than isoindene, the corresponding radical anion isomers have almost the same energy. Nucleus-independent chemical shift calculations show that the highly aromatic molecules lose almost all aromaticity when an extra electron is present. The radical anions of cyclopentadiene and all of its annulated derivatives have remarkably low C-H bond dissociation energies (only 18-34 kcal mol(-1) for the mono-, bi-, and tricyclics considered). Hydrogen atom loss leads to the restoration of aromaticity in the highly stabilized cyclopentadienyl anion congeners.     

Photoelectron imaging of large anionic methanol clusters: (MeOH)n - (n approximately 70-460)

Electron solvation in methanol anion clusters, (MeOH)(n) (-) (n approximately 70-460), is studied by photoelectron imaging. Two isomers are observed: methanol I, with vertical binding energies (VBE) ranging from 2-2.5 eV, and methanol II, with much lower VBE's between 0.2 and 0.5 eV. The VBE's of the two isomers depend linearly on n(-1/3) with nearly identical slopes. We propose that the excess electron is internally solvated in methanol I clusters, whereas in methanol II it resides in a dipole-bound surface-state. Evidence of an excited state accessible at 1.55 eV is observed for methanol I.     

正、负和中性TiP10团簇结构与电子性质的密度泛函研究

采用密度泛函理论的B3LYP方法研究了正、负和中性TiP10团簇的几何构型和电子结构. 计算结果表明, 中性TiP10团簇的基态构型为金属夹心结构, 正、负离子团簇同样具有该基态稳定结构. 通过对基态稳定结构的分子轨道分析表明, δ键对形成夹心结构起到重要作用. 理论计算得到的中性TiP10团簇的垂直和绝热电离能分别为7.84和7.68 eV, 垂直和绝热电子亲和势分别为3.18和3.35 eV.     

Characterization of solvated electrons in hydrogen cyanide clusters: (HCN)n- (n=3, 4)

Theoretical studies of the solvated electrons (HCN)n- (n=3, 4) reveal a variety of electron trapping possibilities in the (HCN)n (n=3, 4) clusters. Two isomers for (HCN)3- and four isomers for (HCN)4- are obtained at the MP2aug-cc-pVDZ+dBF (diffusive bond functions) level of theory. In view of vertical electron detachment energies (VDEs) at the CCSD(T) level, the excess electron always "prefers" locating in the center of the system, i.e., the isomer with higher coordination number shows larger VDE value. However, the most stable isomers of the solvated electron state (HCN)3- and (HCN)4- are found to be the linear Cinfinitynu and Dinfinityh structures, respectively, but not the fullyl symmetric structures which have the largest VDE values.     

Photoionization of C(4) molecular beam: ab initio calculations

Large computations are performed on the C(4) (+) cation in order to characterize its stable isomers and its lowest electronic excited states using configuration interaction methods and large basis sets. Several stable isomers are found including a linear C(4) (+)(l-C(4) (+)), a rhombic C(4) (+)(r-C(4) (+)) (or cyclic), and a branched (d-C(4) (+)) structure. Our calculations show a high density of electronic states for all of these isomers favoring their interactions. By combining the present ab initio data and those on neutral C(4), the l-C(4)(X)+hnu-->l-C(4) (+)(X(+))+e(-), d-C(4)(X)+hnu-->d-C(4) (+)(X(+))+e(-), and r-C(4)(X)+hnu-->r-C(4) (+)(X(+))+e(-) vertical photoionization transition energies are computed at 10.87, 10.92, and 10.77 eV, respectively. Photoionizing a C(4) molecular beam results on an onset at 10.4-10.5 eV and then to a linear increase of the signal due to the opening of several ionization channels involving most of the C(4) and C(4) (+) isomers and electronic states.     

An ab initio MO study on structures and energetics of CH2Si2, CH2Si2, and CH2Si2

The geometric structures and isomeric stabilities of various stationary points in CH₂Si₂ neutral, cation and anion are investigated at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory. For the geometrical survey, the basis sets used are of the cc-pVTZ for the neutral and cation. The final energies are calculated by the use of the CCSD(T) level of theory with the aug-cc-pVTZ basis set at their optimized geometries. To the competitive two-anion isomers, the aug-cc-pVTZ basis sets are applied. The global minimum (N-1) of the CH₂Si₂ neutral has a quite different framework from those of the C₃H₂ (cyclopropenylidene) and Si₃H₂ (trisilacyclopropenylidene) neutrals. No competitive low-lying isomers are found in the CH₂Si₂ neutral. The attractive conformer (C-1) is predicted for the most stable cation, where its framework is quite different from that of the neutral N-1. Both H atoms are connected to the same C atom, but each C–H bond length is different from each other. Two competitive anion isomers with positive (real) electron affinities are predicted. The framework of the most stable anion A-1 is quite similar to that of the cation C-1, whereas both H atoms are equally connected to the same C atom. The framework of the anion isomer A-2 is the same as that in the neutral N-1. The vertical and adiabatic ionization potentials from the most stable neutral N-1 are 9.02 and 8.71 eV, respectively. The adiabatic electron affinity of the lowest lying isomer N-1 is only 0.43 eV and the vertical electron detachment energy form the global minimum anion (A-1) is 2.02 eV. The multi-centered Si–H–Si bonds are found in the neutral, cation, and anion.     

Structure and stability of (CuNO) isomers

Quantum-mechanical calculations have been performed on various isomers of the (CuNO)⁺ system. A ²Π ground state is found for the linear CuNO⁺ and CuON⁺ isomers and a ²A′ state for the bent CuNO⁺ and CuON⁺ isomers. Energy calculations indicate that the linear CuNO⁺ structure is the most stable, the bent CuNO⁺ and CuON⁺ and the linear CuON⁺ structures are at 0.86 eV, 0.99 eV and 1.04 eV above this respectively. In the CuNO⁺ → CuON⁺ interconversion between the linear isomers, three transition states are involved, whereas the bent CuNO⁺ isomer is found to be an intermediate species. The isomerization barriers, dissociation energies, equilibrium geometries and vibration frequencies are given for all isomers in their ground and first excited states.     

Accuracy and limitations of second-order many-body perturbation theory for predicting vertical detachment energies of solvated-electron clusters

Vertical electron detachment energies (VDEs) are calculated for a variety of (H(2)O)(n)(-) and (HF)(n)(-) isomers, using different electronic structure methodologies but focusing in particular on a comparison between second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory with noniterative triples, CCSD(T). For the surface-bound electrons that characterize small (H(2)O)(n)(-) clusters (n< or = 7), the correlation energy associated with the unpaired electron grows linearly as a function of the VDE but is unrelated to the number of monomers, n. In every example considered here, including strongly-bound "cavity" isomers of (H(2)O)(24)(-), the correlation energy associated with the unpaired electron is significantly smaller than that associated with typical valence electrons. As a result, the error in the MP2 detachment energy, as a fraction of the CCSD(T) value, approaches a limit of about -7% for (H(2)O)(n)(-) clusters with VDEs larger than about 0.4 eV. CCSD(T) detachment energies are bounded from below by MP2 values and from above by VDEs calculated using second-order many-body perturbation theory with molecular orbitals obtained from density functional theory. For a variety of both strongly- and weakly-bound isomers of (H(2)O)(20)(-) and (H(2)O)(24)(-), including both surface states and cavity states, these bounds afford typical error bars of +/-0.1 eV. We have found only one case where the Hartree-Fock and density functional orbitals differ qualitatively; in this case the aforementioned bounds lie 0.4 eV apart, and second-order perturbation theory may not be reliable.