Investigating the Effects of Basis Set on Metal-Metal and Metal-Ligand Bond Distances in Stable Transition Metal Carbonyls: Performance of Correlation Consistent Basis Sets with 35 Density Functionals.

Density functional theory (DFT) is a widely used method for predicting equilibrium geometries of organometallic compounds involving transition metals, with a wide choice of functional and basis set combinations. A study of the role of basis set size in predicting the structural parameters can be insightful with respect to the effectiveness of using small basis sets to optimize larger molecular systems. For many organometallic systems, the metal-metal and metal-carbon distances are the most important structural features. In this study, we compare the equilibrium metal-ligand and metal-metal distances of six transition metal carbonyl compounds predicted by the Hood-Pitzer double-ζ polarization (DZP) basis set, against those predicted employing the standard correlation consistent cc-pVXZ (X = D,T,Q) basis sets, for 35 different DFT methods. The effects of systematically increasing the basis set size on the structural parameters are carefully investigated. The Mn-Mn bond distance in Mn2(CO)10 shows a greater dependence on basis set size compared to the other M-M bonds. However, the DZP predictions for re(Mn-Mn) are closer to experiment than those obtained with the much larger cc-pVQZ basis set. Our results show that, in general, DZP basis sets predict structural parameters with an accuracy comparable to the triple and quadruple-ζ basis sets. This finding is very significant, because the quadruple-ζ basis set for Mn2(CO)10 includes 1308 basis functions, while the equally effective double-ζ set (DZP) includes only 366 basis functions. Overall, the DZP M06-L method predicts structures that are very consistent with experiment.

Hydrogen-abstracted adenine-thymine radicals with interesting transferable properties.

The formation of radicals on DNA bases through various pathways can lead to harmful structural alterations. Such processes are of interest for preventing alteration of healthy DNA and, conversely, to develop more refined methods for inhibiting the replication of unwanted mutagenic DNA. In the present work, we explore theoretically the energetic and structural properties of the nine possible neutral radicals formed via hydrogen abstraction from the adenine-thymine base pair. The lowest energy radical is formed by loss of a hydrogen atom from the methyl group of thymine. The next lowest energy radicals, lying 8 and 9 kcal mol-1 higher than the global minimum, are those in which hydrogens are removed from the two nitrogens that would join the base pair to 2-deoxyribose in double-stranded DNA. The other six radicals lie between 16 and 32 kcal mol-1 higher in energy. Unlike the guanine-cytosine base pair, adenine-thymine (A-T) exhibits only minor structural changes upon hydrogen abstraction, with all A-T derived radicals maintaining planarity. Moreover, the energetic ordering for the radicals of the two isolated bases (adenine and thymine) is preserved upon formation of the base pair, though with a wider spread of energies. Even more significantly, the energetic interleaving of the (A-H)*-T and A-(T-H)* radicals is correctly predicted from the X-H bond dissociation energies of the isolated adenine and thymine. This suggests that the addition of the hydrogen-bonded complement base only marginally affects the bond energies.

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne.

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.

The deprotonated guanine-cytosine base pair.

Awareness of the harmful effects of radiation has increased interest in finding the mechanisms of DNA damage. Radical and anion formation among the DNA base pairs are thought to be important steps in such damage [Collins, G. P. (2003) Sci. Am. 289 (3), 26-27]. Energetic properties and optimized geometries of 10 radicals and their respective anions derived through hydrogen abstraction from the Watson-Crick guanine-cytosine (G-C) base pair have been studied using reliable theoretical methods. The most favorable deprotonated structure (dissociation energy 42 kcal x mol(-1), vertical detachment energy 3.79 eV) ejects the proton analogous to the cytosine glycosidic bond in DNA. This structure is a surprisingly large 12 kcal x mol(-1) lower in energy than any of the other nine deprotonated G-C structures. This system retains the qualitative G-C structure but with the H...O2 distance dramatically reduced from 1.88 to 1.58 A, an extremely short hydrogen bond. The most interesting deprotonated G-C structure is a "reverse wobble" incorporating two N-H...N hydrogen bonds. Three different types of relaxation energies (4.3-54 kcal x mol(-1)) are defined and reported to evaluate the energy released via different mechanisms for the preparation of the deprotonated species. Relative energies, adiabatic electron affinities (ranging from 1.93 to 3.65 eV), and pairing energies are determined to discern which radical will most alter the G-C properties. The most stable deprotonated base pair corresponds to the radical with the largest adiabatic electron affinity, 3.65 eV. This value is an enormous increase over the electron affinity (0.60 eV) of the closed-shell G-C base pair.

DNA nucleosides and their radical anions: molecular structures and electron affinities.

The deoxyribonucleosides have been studied to determine the properties of combinations of 2-deoxyribose with each of the isolated DNA bases for both neutral and anionic species. We have used a carefully calibrated theoretical method [Chem. Rev. 2002, 102, 231], employing the B3LYP hybrid Hartree-Fock/DFT functional with the DZP++ basis set. Predictions are made of the geometric parameters, adiabatic electron affinities, charge distributions based on natural population analysis, and decomposition enthalpy for the neutral and anionic forms of the four 2'-deoxyribonucleosides in DNA: 2'-deoxyriboadenosine (dA), 2'-deoxyribocytidine (dC), 2'-deoxyriboguanosine (dG), and 2'-deoxyribothymidine (dT). Geometric changes in the anions show that the glycosidic bond exhibits little change with excess charge for the guanosine but significant shortening for the adenosine and for the pyrimidines. The zero-point corrected adiabatic electron affinities in eV for each of the 2'-deoxyribonucleosides are as follows: 0.06, dA; 0.09, dG; 0.33, dC; and 0.44, dT. These values are uniformly greater than those of the corresponding isolated bases (-0.28, A; -0.07, G; 0.03, C; 0.20, T) and the isolated 2-deoxyribose (-0.38) at the same level of theory. The vertical detachment energies of dT and dC are substantial, 0.72 and 0.94 eV, and these anions should be observable. A high VDE, 0.91 eV, is also found for dA but its anion is unlikely to be stable due to the small AEA of 0.06 eV. The high VDE reflects the fact that the molecular structures of the anions and the corresponding neutral species are quite different. Valence character is displayed for the SOMOs of dA, dC, and dT, while some component of diffuse character is visible in the SOMO of dG. Further analysis of electronic changes upon electron attachment include an examination of the NPA charges, which show that in the neutral 2'-deoxyribonucleosides the sum of NPA charges for every base is the same, -0.28 with the sum of 2-deoxyribose charges being positive, +0.28. In the anions, the trend in charge division varies based on the nature of the excess electron in the anions. Thermodynamically, the overall enthalpy change for the reaction of water with the neutral nucleosides to give bases and ribose is approximately zero. The analogous decomposition is exothermic by 8 to 11 kcal mol-1 for the anions, indicating possible challenges for anionic gas-phase nucleoside exploration in the presence of water.

The rule breaking Cr2(CO)10. A 17 electron Cr system or a Cr=Cr double bond?

Density functional theory (DFT) has been used to investigate the conformations and thermochemistry on the singlet and triplet potential energy surfaces (PES) of Cr2(CO)10. The global minimum energy structure for the lowest singlet state of C2h symmetry is consistent with a model of two interacting Cr(CO)5 fragments in which one carbonyl in each fragment acts as an asymmetric four-electron donor bridging carbonyl, with chromium-chromium distances of 2.93 A (B3LYP) or 2.83 A (BP86). Avoiding a Cr...Cr bond by incorporating four-electron donor CO groups in this way allows each chromium atom in singlet Cr(CO)10 to attain the favored 18-electron configuration by using, in a simple picture of the bonding, only the six octahedral sp3d2 hybrids. The dissociation energy to two Cr(CO)5 fragments or to Cr(CO)6 + Cr(CO)4 fragments is predicted to be 10 kcal mol(-1). The lowest triplet state of Cr2(CO)10 is predicted to lie approximately 10 kcal mol(-1) above the singlet global minimum. In the case of triplet Cr2(CO)10 the lowest energy minima were found to be of C2 and C2h symmetry, with similar energies. The chromium-chromium distances in triplet Cr2(CO)10 were found to be shorter than those in the corresponding singlet structures, namely 2.81 (B3LYP) or 2.68 A (BP86) suggesting a sigma + 2(1/2) pi Cr=Cr double bond similar to the O=O bond in O2 or the Fe=Fe bond in the experimentally observed triplet state (Me5C5)2Fe2(mu-CO)3.

Electron affinity of the guanine-cytosine base pair and structural perturbations upon anion formation.

The adiabatic electron affinity (AEA) for the Watson-Crick guanine-cytosine (GC) DNA base pair is predicted using a range of density functional methods with double- and triple-zeta plus polarization plus diffuse (DZP++ and TZ2P++) basis sets in an effort to bracket the true electron affinity. The methods used have been calibrated against a comprehensive tabulation of experimental electron affinities (Chem.Rev. 2002, 102, 231). Optimized structures for GC and the GC anion are compared to the neutral and anionic forms of the individual bases as well as Rich's 1976 X-ray structure for sodium guanylyl-3',5'-cytidine nonahydrate, GpC.9H(2)O. Structural distortions and natural population (NPA) charge distributions of the GC anion indicate that the unpaired electron is localized primarily on the cytosine moiety. Unlike treatments using second-order perturbation theory (MP2), density functional theory consistently predicts a substantial positive adiabatic electron affinity for the GC pair (e.g., TZ2P++/B3LYP: +0.48 eV). The stabilization of C(-) via three hydrogen bonds to guanine is sufficient to facilitate adiabatic binding of an electron to GC and is also consistent with the positive experimental electron affinities obtained by photoelectron spectroscopy of cytosine anions incrementally microsolvated with water molecules. The pairing (dissociation) energy for GC(-) (35.6 kcal/mol) is determined with inclusion of electron correlation and shows the anion to have greater thermodynamic stability; the pairing energy for neutral GC (TZ2P++/B3LYP 23.9 kcal/mol) compares favorably to previous MP2/6-31G (23.4 kcal/mol) results and a debated experiment (21.0 kcal/mol).

Examination of the Stabilities of Group 14 (C, Si, Ge, Sn, Pb) Congeners of Dihydroxycarbene and Dioxirane. Comparison to Formic Acid and Hydroperoxycarbene Congeners.

The relative energetics of four XH(2)O(2) (X = C, Si, Ge, Sn, Pb) isomers, dihydroxycarbene, formic acid, dioxirane, and hydroperoxycarbene, were determined using the BLYP and B3LYP density functionals with DZP and TZ2P basis sets, as well as CCSD and CCSD(T) single-point energies at the BLYP/TZ2P optimized geometries. Relative to dihydroxycarbene, formic acid was 41.8 kcal/mol lower in energy while dioxirane and hydroperoxycarbene were 51.3 and 63.6 kcal/mol higher, respectively, with CCSD(T). Furthermore, using an effective core potential (ECP) the dihydroxy congener was shown to be the most stable isomer for X = Si-Pb. The formic acid and dioxirane congeners become increasingly less stable as one descends group 14. Our results show that divalency is preferred for Si-Pb (dihydroxy congeners are the most stable) but the tetravalent formic acid congeners remain more stable than the hydroperoxy congeners, showing that divalency is not universally preferred among these isomers.