Theoretical study of the valence ionization energies and electron affinities of linear C2n+1 (n=1–6) clusters

Electron Propagator Theory of Vertical Electron Detachment Energies of Anions: Benchmarks and Applications to Nucleotides

A new generation of ab initio electron-propagator self-energy approximations that are free of adjustable parameters is tested on a benchmark set of 55 vertical electron detachment energies of closed-shell anions. Comparisons with older self-energy approximations indicate that several new methods that make the diagonal self-energy approximation in the canonical Hartree-Fock orbital basis provide superior accuracy and computational efficiency. These methods and their acronyms, mean absolute errors (in eV), and arithmetic bottlenecks expressed in terms of occupied (O) and virtual (V) orbitals are the opposite-spin, non-Dyson, diagonal second-order method (os-nD-D2, 0.2, OV²), the approximately renormalized quasiparticle third-order method (Q3+, 0.15, O²V³) and the approximately renormalized, non-Dyson, linear, third-order method (nD-L3+, 0.1, OV⁴). The Brueckner doubles with triple field operators (BD-T1) nondiagonal electron-propagator method provides such close agreement with coupled-cluster single, double, and perturbative triple replacement total energy differences that it may be used as an alternative means of obtaining standard data. The new methods with diagonal self-energy matrices are the foundation of a composite procedure for estimating basis-set effects. This model produces accurate predictions and clear interpretations based on Dyson orbitals for the photoelectron spectra of the nucleotides found in DNA.

Understanding the antioxidant behavior of some vitamin molecules: a first-principles density functional approach

The structures, energetics, vertical and adiabatic ionization potentials, electron affinities, and global reactivity descriptors of antioxidant vitamins (both water- and fat-soluble) in neutral, positively charged, and negatively charged states were investigated theoretically. We worked within the framework of first-principles density functional theory (DFT), using the B3LYP functional and both localized (6-311G+(d,p) and plane-wave basis sets. Solvent effects were modeled via the polarizable continuum model (PCM), using the integral equation formalism variant (IEFPCM). From the computed structural parameters, ionization potentials, electron affinities, and spin densities, we deduced that these vitamins prefer to lose electrons to neutral reactive oxygen species (·OH and ·OOH), making them good antioxidants. Conceptual DFT was used to determine global chemical reactivity parameters. The computed chemical hardnesses showed that these antioxidant vitamins are more reactive than neutral reactive oxygen species (ROS), thus supporting their antioxidant character towards these species. However, in the neutral state, these vitamins do not act as antioxidants for [Formula: see text]. The reactivity of vitamins towards ROS depends on the nature of the solvent. Amongst the ROS, ·OH has the greatest propensity to attract electrons from a generic donor. The reactivities of fat-soluble vitamins towards neutral reactive oxygen species were found to be larger than those of water-soluble vitamins towards these species, showing that the former are better antioxidants.

Formation and emission of gold and silver carbide cluster ions in a single C60- surface impact at keV energies: experiment and calculations

Impact of fullerene ions (C(60)(-)) on a metallic surface at keV kinetic energies and under single collision conditions is used as an efficient way for generating gas phase carbide cluster ions of gold and silver, which were rarely explored before. Positively and negatively charged cluster ions, Au(n)C(m)(+) (n = 1-5, 1 ≤ m ≤ 12), Ag(n)C(m)(+) (n = 1-7, 1 ≤ m ≤ 7), Au(n)C(m)(-) (n = 1-5, 1 ≤ m ≤ 10), and Ag(n)C(m)(-) (n = 1-3, 1 ≤ m ≤ 6), were observed. The Au(3)C(2)(+) and Ag(3)C(2)(+) clusters are the most abundant cations in the corresponding mass spectra. Pronounced odd/even intensity alternations were observed for nearly all Au(n)C(m)(+/-) and Ag(n)C(m)(+/-) series. The time dependence of signal intensity for selected positive ions was measured over a broad range of C(60)(-) impact energies and fluxes. A few orders of magnitude immediate signal jump instantaneous with the C(60)(-) ion beam opening was observed, followed by a nearly constant plateau. It is concluded that the overall process of the fullerene collision and formation∕ejection of the carbidic species can be described as a single impact event where the shattering of the incoming C(60)(-) ion into small C(m) fragments occurs nearly instantaneously with the (multiple) pickup of metal atoms and resulting emission of the carbide clusters. Density functional theory calculations showed that the most stable configuration of the Au(n)C(m)(+) (n = 1, 2) clusters is a linear carbon chain with one or two terminal gold atoms correspondingly (except for a bent configuration of Au(2)C(+)). The calculated AuC(m) adiabatic ionization energies showed parity alternations in agreement with the measured intensity alternations of the corresponding ions. The Au(3)C(2)(+) ion possesses a basic Au(2)C(2) acetylide structure with a π-coordinated third gold atom, forming a π-complex structure of the type [Au(π-Au(2)C(2))](+). The calculation shows meaningful contributions of direct gold-gold bonding to the overall stability of the Au(3)C(2)(+) complex.

Energetics of linear carbon chains in one-dimensional restricted environment

The energetics of even and odd linear C(n) carbon chain clusters are investigated by hybrid density functional theory (DFT) calculations. These molecular species are especially interesting due to their recent observation inside carbon nanotubes by polarized resonant Raman spectroscopy and high-resolution transmission electron microscopy (HRTEM) by different research groups. Neutral, anionic and dianionic carbon chains were studied with sizes up to n=75, although most presented calculations are limited to n<or= 24. Aggregation into longer chains is favored for neutral and anionic chains of any size. The barrier to aggregation of 2C(n)<-->C(2n) is of the order of 40-20 kcal mol(-1), which gradually decreases with increasing chain size, n. These barriers can be overcome during the high temperature synthesis or annealing conditions, but not when cooled down for the HRTEM and Raman experiments. Therefore, in addition to the already observed long chains also shorter chains should be observable under appropriate conditions inside carbon nanotubes.

Ionization Thresholds of Small Carbon Clusters: Tunable VUV Experiments and Theory

Small carbon clusters (Cn, n = 2-15) are produced in a molecular beam by pulsed laser vaporization and studied with vacuum ultraviolet (VUV) photoionization mass spectrometry. The required VUV radiation in the 8-12 eV range is provided by the Advanced Light Source (ALS) at the Lawrence Berkeley National Laboratory. Mass spectra at various ionization energies reveal the qualitative relative abundances of the neutral carbon clusters produced. By far the most abundant species is C3. Using the tunability of the ALS, ionization threshold spectra are recorded for the clusters up to 15 atoms in size. The ionization thresholds are compared to those measured previously with charge-transfer bracketing methods. To interpret the ionization thresholds for different cluster sizes, new ab initio calculations are carried out on the clusters for n = 4-10. Geometric structures are optimized at the CCSD(T) level with cc-pVTZ (or cc-pVDZ) basis sets, and focal point extrapolations are applied to both neutral and cation species to determine adiabatic and vertical ionization potentials. The comparison of computed and measured ionization potentials makes it possible to investigate the isomeric structures of the neutral clusters produced in this experiment. The measurements are inconclusive for the n = 4-6 species because of unquenched excited electronic states. However, the data provide evidence for the prominence of linear structures for the n = 7, 9, 11, 13 species and the presence of cyclic C10.

Higher-order equation-of-motion coupled-cluster methods for electron attachment

High-order equation-of-motion coupled-cluster methods for electron attachment (EA-EOM-CC) have been implemented with the aid of the symbolic algebra program TCE into parallel computer programs. Two types of size-extensive truncation have been applied to the electron-attachment and cluster excitation operators: (1) the electron-attachment operator truncated after the 2p-1h, 3p-2h, or 4p-3h level in combination with the cluster excitation operator after doubles, triples, or quadruples, respectively, defining EA-EOM-CCSD, EA-EOM-CCSDT, or EA-EOM-CCSDTQ; (2) the combination of up to the 3p-2h electron-attachment operator and up to the double cluster excitation operator [EA-EOM-CCSD(3p-2h)] or up to 4p-3h and triples [EA-EOM-CCSDT(4p-3h)]. These methods, capable of handling electron attachment to open-shell molecules, have been applied to the electron affinities of NH and C2, the excitation energies of CH, and the spectroscopic constants of all these molecules with the errors due to basis sets of finite sizes removed by extrapolation. The differences in the electron affinities or excitation energies between EA-EOM-CCSD and experiment are frequently in excess of 2 eV for these molecules, which have severe multideterminant wave functions. Including higher-order operators, the EA-EOM-CC methods predict these quantities accurate to within 0.01 eV of experimental values. In particular, the 3p-2h electron-attachment and triple cluster excitation operators are significant for achieving this accuracy.

Molecular and Electronic Structures by Design: Tuning Symmetrical and Unsymmetrical Linear Trichromium Chains

The preparation, properties, and crystal structures of 12 trichromium extended metal atom chain (EMAC) compounds of the type Cr(3)(L)(4)X(2) (L = equatorial ligands dipyridylamide (dpa) or di-4,4'-ethyl-2,2'-pyridylamide (depa), and X = axial ligands, e.g., halide or pseudohalide ions) with large variations in metal-metal distances are reported here. These complexes, which belong to a broad class of fundamentally interesting trinuclear molecules over which the electrons may or may not be delocalized, pose significant theoretical and experimental challenges which are dealt with in this report. Complexes with strongly donating axial or equatorial ligands tend to favor a symmetrical (D(4)) molecular structure, while more weakly donating ligands give rise to unsymmetrical (C(4)) structures; the physical properties of these two classes of compounds are discussed fully, and important comparisons with a reported DFT model of the electronic structures of the compounds are made.

Performance of time-dependent density functional and Green functions methods for calculations of excitation energies in radicals and for Rydberg electronic states

Time-dependent density functional (TD-DFT) and perturbation theory-based outer valence Green functions (OVGF) methods have been tested for calculations of excitation energies for a set of radicals, molecules, and model clusters simulating points defects in silica. The results show that the TD-DFT approach may give unreliable results not only for diffuse Rydberg states, but also for electronic states involving transitions between MOs localized in two remote from each other spatial regions, for example, for charge-transfer excitations. For the. O-SiX(3) clusters, where X is a single-valence group, TD-DFT predicts reasonable excitation energies but incorrect sequence of electronic transitions. For a number of cases where TD-DFT is shown to be unreliable, the OVGF approach can provide better estimates of excitation energies, but this method also is not expected to perform universally well. The OVGF performance is demonstrated to be satisfactory for excitations with predominantly single-determinant wave functions where the deviations of the calculated energies from experiment should not exceed 0.1-0.3 eV. However, for more complicated transitions involving multiple bonds or for excited states with multireference wave functions the OVGF approach is less reliable and error in the computed energies can reach 0.5-1 eV.

What is the nature of polyacetylene neutral and anionic chains HC(2n)H and HC(2n)H(-) (n = 6-12) that have recently been observed?

The optimized geometries, adiabatic electron affinities, and IR-active vibrational frequencies have been predicted for the long linear carbon chains HC(2n)H. The B3LYP density functional combined with the DZP basis set was used in this theoretical study. The computed physical properties are discussed. The predicted electron affinities form a remarkably regular sequence: 1.78 (HC(12)H), 2.08 (HC(14)H), 2.32 (HC(16)H), 2.53 (HC(18)H), 2.69 (HC(20)H), 2.83 (HC(22)H), and 2.95 eV (HC(24)H). The predicted structures display an alternating triple and very short single bond pattern, with the degree of bond alternation significantly less for the radical anions.