Assessing the Viability of the Methylsulfinyl Radical-Ozone Reaction.

Although integral to remote marine atmospheric sulfur chemistry, the reaction between methylsulfinyl radical (CH3 SO) and ozone poses challenges to theoretical treatments. The lone theoretical study on this reaction reported an unphysically large barrier of 66 kcal mol-1 for abstraction of an oxygen atom from O3 by CH3 SO. Herein, we demonstrate that this result stems from improper use of MP2 with a single-reference, unrestricted Hartree-Fock (UHF) wavefunction. We characterized the potential energy surface using density functional theory (DFT), as well as multireference methodologies employing a complete active-space self-consistent field (CASSCF) reference. Our DFT PES shows, in contrast to previous work, that the reaction proceeds by forming an addition adduct [CH3 S(O3 )O] in a deep potential well of 37 kcal mol-1 . An O-O bond of this adduct dissociates via a flat, low barrier of 1 kcal mol-1 to give CH3 SO2 +O2 . The multireference computations show that the initial addition of CH3 SO+O3 is barrierless. These results provide a more physically intuitive and accurate picture of this reaction than the previous theoretical study. In addition, our results imply that the CH3 SO2 formed in this reaction can readily decompose to give SO2 as a major product, in alignment with the literature on CH3 SO reactions.

Energetics and transition-state dynamics of the F + HOCH3 → HF + OCH3 reaction.

The F + HOCH3 → HF + OCH3 reaction is a system with 15 internal degrees of freedom that can provide a benchmark for the development of theory for increasingly complex chemical reactions. The dynamics of this reaction were studied by photoelectron-photofragment coincidence (PPC) spectroscopy carried out on the F-(HOCH3) anion, aided by a computational study of both the anion and neutral potential energy surfaces, with energies extrapolated to the CCSDT(Q)/CBS level of theory. Photodetachment at 4.80 eV accesses both the reactant and product channels for this reaction. In the product channel (HF + OCH3 + e-) of the neutral potential energy surface, vibrationally excited HF products in addition to the stable product-channel hydrogen-bonded complex (FH-OCH3) are observed in the PPC and photoelectron spectra. In addition, experimental evidence is observed for the reactant-channel van der Waals complex (F-HOCH3), in good agreement with the theoretical predictions. The relative stability of these long-lived complexes was probed by reducing the ion beam energy, increasing the product time-of-flight, indicating lifetimes on the microsecond timescale for the reactant- and product-channel complexes as well as providing evidence for long-lived vibrational Feshbach resonances associated with the HF(v > 0) + OCH3 product states. This system will provide a model for extending full-dimensionality quantum dynamics to larger numbers of degrees of freedom.

Photodissociation of Cerium Oxide Nanocluster Cations.

Cerium oxide cluster cations, CexOy(+), are produced via laser vaporization in a pulsed nozzle source and detected with time-of-flight mass spectrometry. The mass spectrum displays a strongly preferred oxide stoichiometry for each cluster with a specific number of metal atoms x, with x ≤ y. Specifically, the most prominent clusters correspond to the formula CeO(CeO2)n(+). The cluster cations are mass selected and photodissociated with a Nd:YAG laser at either 532 or 355 nm. The prominent clusters dissociate to produce smaller species also having a similar CeO(CeO2)n(+) formula, always with apparent leaving groups of (CeO2). The production of CeO(CeO2)n(+) from the dissociation of many cluster sizes establishes the relative stability of these clusters. Furthermore, the consistent loss of neutral CeO2 shows that the smallest neutral clusters adopt the same oxidation state (IV) as the most common form of bulk cerium oxide. Clusters with higher oxygen content than the CeO(CeO2)n(+) masses are present with much lower abundance. These species dissociate by the loss of O2, leaving surviving clusters with the CeO(CeO2)n(+) formula. Density functional theory calculations on these clusters suggest structures composed of stable CeO(CeO2)n(+) cores with excess oxygen bound to the surface as a superoxide unit (O2(-)).

Molybdenum-molybdenum multiple bonding in homoleptic molybdenum carbonyls: comparison with their chromium analogues.

The binuclear molybdenum carbonyls Mo(2)(CO)(n) (n = 11, 10, 9, 8) have been studied by density functional theory using the BP86 and MPW1PW91 functionals. The lowest energy Mo(2)(CO)(11) structure is a singly bridged singlet structure with a Mo-Mo single bond. This structure is essentially thermoneutral toward dissociation into Mo(CO)(6) + Mo(CO)(5), suggesting limited viability similar to the analogous Cr(2)(CO)(11). The lowest energy Mo(2)(CO)(10) structure is a doubly semibridged singlet structure with a Mo═Mo double bond. This structure is essentially thermoneutral toward disproportionation into Mo(2)(CO)(11) + Mo(2)(CO)(9), suggesting limited viability. The lowest energy Mo(2)(CO)(9) structure has three semibridging CO groups and a Mo≡Mo triple bond analogous to the lowest energy Cr(2)(CO)(9) structure. This structure appears to be viable toward CO dissociation, disproportionation into Mo(2)(CO)(10) + Mo(2)(CO)(8), and fragmentation into Mo(CO)(5) + Mo(CO)(4) and thus appears to be a possible synthetic objective. The lowest energy Mo(2)(CO)(8) structure has one semibridging CO group and a Mo≡Mo triple bond similar to that in the lowest energy Mo(2)(CO)(9) structure. This differs from the lowest energy Cr(2)(CO)(8) structure, which is a triply bridged structure. A higher energy unbridged D(2d) Mo(2)(CO)(8) structure was found with a very short Mo-Mo distance of 2.6 Å. This interesting structure has two degenerate imaginary vibrational frequencies. Following the corresponding normal modes leads to a Mo(2)(CO)(8) structure, lying ~5 kcal/mol above the global minimum, with two four-electron donor bridging CO groups and a Mo═Mo distance suggesting a formal double bond. All of the triplet Mo(2)(CO)(n) (n = 10, 9, 8) structures were found to be relatively high energy structures, lying at least 22 kcal/mol above the corresponding global minimum. The singlet-triplet splittings for the Mo(2)(CO)(n) (n = 10, 9, 8) structures are significantly higher than those of the Cr(2)(CO)(n) analogues. The Mo-Mo Wiberg bond indices confirm our assigned bond orders based on predicted bond distances.

Monobridged Si2H4.

The rotational spectrum of a new monobridged isomer of Si(2)H(4), denoted here as H(2)Si(H)SiH, has been detected by Fourier transform microwave spectroscopy of a supersonic molecular beam through the discharge products of silane. On the basis of high-level coupled cluster theory, this isomer is calculated to lie only 7 kcalmol above disilene (H(2)SiSiH(2)), the most stable isomeric arrangement of Si(2)H(4), and to be fairly polar, with a calculated dipole moment of mu = 1.14 D. The rotational spectrum of H(2)Si(H)SiH exhibits closely spaced line doubling, characteristic of a molecule undergoing high-frequency inversion. Transition state calculations indicate that inversion probably occurs in two steps: migration of the bridged hydrogen atom to form silylsilylene, H(3)SiSiH, and then internal rotation of the SiH(3) group, followed by the reverse process. The potential energy surface for this type of inversion is quite shallow, with a barrier height of only 2-3 kcalmol. Searches for the rotational lines of silylsilylene, calculated to be of comparable stability to H(2)Si(H)SiH but about five times less polar (mu = 0.23 D), have also been undertaken, so far without success, even though strong lines of H(2)Si(H)SiH have been detected. The favorable energetics and high polarity of monobridged Si(2)H(4) with respect to either disilene or silylsilylene make it a plausible candidate for radioastronomical detection in sources such as IRC + 10216, where comparably large silicon molecules such as SiS, SiC(3), and SiC(4) have already been discovered.

Analysis of the origin of through-space proton NMR deshielding by selected organic functional groups.

GIAO-HF and IGLO-DFT computations of isotropic magnetic shieldings were used to map the NMR shielding environments of small molecules exemplifying selected organic functional groups. Two different probes were employed: a methane molecule and NICS (nucleus-independent chemical shifts) based on computed absolute isotropic shieldings. The reason for the different results obtained using these two probes is perturbation of the wave function by the proximity of methane to the pi bond, as analyzed by the localized orbital contributions to the shieldings. [structure: see text]

Binuclear homoleptic copper carbonyls Cu(2)(CO)(x) (x = 1-6): remarkable structures contrasting metal-metal multiple bonding with low-dimensional copper bonding manifolds.

Binuclear homoleptic copper carbonyls Cu(2)(CO)(x) (x = 1-6) have been studied using four different density functional theory methods (DFT) in conjunction with a basis set of extended double-zeta plus polarization quality, labeled as DZP. For each homoleptic binuclear copper carbonyl compound, several stationary point structures are presented, and these structures are characterized in terms of their geometries, thermochemistry, and vibrational frequencies. The optimal unsaturated Cu(2)(CO)(x) (x = 1-6) structures are generated by joining 18-electron tetrahedral, 16-electron trigonal, 14-electron linear copper carbonyl building blocks, and/or bare copper atoms with copper-copper single bonds rather than by joining 18-electron copper carbonyl units with multiple copper-copper bonds. For Cu(2)(CO)(6) the eclipsed and staggered ethane-like structure are virtually degenerate and lie significantly lower in energy than other possible structures. The eclipsed Cu-Cu single bond distance is predicted to be 2.61 A, while that for the staggered structure is 2.65 A. The lowest energy structure for Cu(2)(CO)(5) is the eclipsed ethyl radical-like structure, with r(e)(Cu-Cu) = 2.51 A. The staggered ethyl radical-like structure lies only 0.1 kcal/mol higher in energy, with a Cu-Cu distance shorter by only approximately 0.001 A. For Cu(2)(CO)(4) a methylcarbene-like structure is predicted to lie lowest, with Cu-Cu distance 2.40 A. However, twisted and planar ethylene-like structure lie only 3-5 kcal/mol higher. For Cu(2)(CO)(3) a surprising methylcarbyne-like structure with r(e)(Cu-Cu) = 2.38 A is predicted to lie lowest with all four DFT methods. However, a classical vinyl radical-like lies only 2-4 kcal/mol higher. For Cu(2)(CO)(2) theory predicts a vinylidene-like structure with r(e)(Cu-Cu) = 2.34 A to be essentially degenerate with cis and trans bent acetylene structures with copper-copper distances 2.33 A. Finally, and consistent with earlier theoretical studies, the linear end on Cu-Cu-CO structure with r(e)(Cu-Cu) = 2.27 A is the predicted global minimum for Cu(2)(CO).

Electron affinities of the DNA and RNA bases.

Adiabatic electron affinities (AEAs) for the DNA and RNA bases are predicted by using a range of density functionals with a double-zeta plus polarization plus diffuse (DZP++) basis set in an effort to bracket the true EAs. Although the AEAs exhibit moderate fluctuations with respect to the choice of functional, systematic trends show that the covalent uracil (U) and thymine (T) anions are bound by 0.05-0.25 eV while the adenine (A) anion is clearly unbound. The computed AEAs for cytosine (C) and guanine (G) oscillate between small positive and negative values for the three most reliable functional combinations (BP86, B3LYP, and BLYP), and it remains unclear if either covalent anion is bound. AEAs with B3LYP/TZ2P++ single points are 0.19 (U), 0.16 (T), 0.07 (G), -0.02 (C), and -0.17 eV (A). Favorable comparisons are made to experimental estimates extrapolated from photoelectron spectra data for the complexes of the nucleobases with water. However, experimental values scaled from liquid-phase reduction potentials are shown to overestimate the AEAs by as much as 1.5 eV. Because the uracil and thymine covalent EAs are in energy ranges near those of their dipole-bound counterparts, preparation and precise experimental measurement of the thermodynamically stable covalent anions may prove challenging.

The almost bottleable triplet carbene: 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene.

Computations on 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene (1) using ab initio and density functional theory methods underscore the unusual stability of the triplet over the singlet state. At the B3LYP/6-311G(d,p) level, the triplet state had a slightly bent central C-C-C bond angle of 167 degrees, whereas this angle in the singlet was 134 degrees. The B3LYP singlet-triplet splitting (12.2 kcal/mol) was larger than that of the parent molecule (5.8 kcal/mol), diphenylcarbene (2), which also has a triplet ground state. The energy of a suitable isodesmic reaction showed the triplet and singlet states of (1) to be destabilized, by 6.3 and 12.5 kcal/mol, respectively, due to the combined effects of the CF3, Br, and alkyl substituents. The linear-coplanar form of (3)(1), which might facilitate dimerization or electrophilic attack at the more exposed diradical center, was prohibitively (35.9 kcal/mol) higher in energy. Our results confirm Tomioka's conclusion that the triplet diarylcarbene, ortho-substituted with bulky CF3 and Br substituents, is persistent due to steric protection of the diradical center. Dimerization and other possible reaction pathways are inhibited, not only by the bulky ortho substituents but also by the para alkyl groups. The increase in stability of the triplet ((3)(1)) state relative to the singlet ((1)(1)) state does not influence the reactivity directly.

The photohydration of N-alkylpyridinium salts: theory and experiment.

Bicyclic aziridines formed by the irradiation of pyridinium salts in basic solution have recently been recognized to have great synthetic potential. We have undertaken a joint computational and experimental investigation of the mechanism of this photoreaction. We have computationally determined the structures and relative energies of the relevant stationary points on the lowest potential energy surface (PES) of the pyridinium and methylpyridinium ions. Two important intermediates are shown to be bound minima on the ground-state PES: azoniabenzvalene and a 6-aza[3.1.0]bicyclic ion with an exo-oriented substituent (analogous to prefulvene). We advance a mechanism which involves initial formation of this exo-bicyclic ion, followed by nitrogen migration around the ring via the azoniabenzvalene intermediate. Thus, the barrier separating the two intermediates is the factor that determines the degree of scrambling observed in the photoproducts when the carbon atoms are labeled with deuterium or substituted with additional methyl groups. For N-methylpyridinium, the exo-methyl bicyclic ion was computed to be approximately 1 kcalmol(-1) lower in energy than N-methyl-azoniabenzvalene. The transition state was computed to lie several kcal mol(-1) above the exo-methyl bicyclic ion (+8.4kcalmol(-1), 6-31G* RHF; +3.7kcalmol(-1), 6-31G* B3LYP), but still well below the energy available from the 254 nm excitation of the N-methylpyridinium ion. The computed relative energies correspond splendidly with several experimental findings which include the preference for exo products, the results of deuterium labeling, and the impact of additional substituent methyl groups on the product distribution.

A combined crossed molecular beam and ab initio investigation of C2 and C3 elementary reactions with unsaturated hydrocarbons--pathways to hydrogen deficient hydrocarbon radicals in combustion flames.

Crossed molecular beam experiments on dicarbon and tricarbon reactions with unsaturated hydrocarbons acetylene, methylacetylene, and ethylene were performed to investigate the dynamics of channels leading to hydrogen-deficient hydrocarbon radicals. In the light of the results of new ab initio calculations, the experimental data suggest that these reactions are governed by an initial addition of C2/C3 to the pi molecular orbitals forming highly unsaturated cyclic structures. These intermediates are connected via various transition states and are suggested to ring open to chain isomers which decompose predominantly by displacement of atomic hydrogen, forming C4H, C5H, HCCCCCH2, HCCCCCCH3, H2CCCCH and H2CCCCCH. The C2(1 sigma g+) + C2H4 reaction has no entrance barrier and the channel leading to the H2CCCCH product is strongly exothermic. This is in strong contrast with the C3(1 sigma g+) + C2H4 reaction as this is characterized by a 26.4 kJ mol-1 threshold to form a HCCCCCH2 isomer. Analogous to the behavior with ethylene, preliminary results on the reactions of C2 and C3 with C2H2 and CH3CCH showed the H-displacement channels of these systems to share many similarities such as the absence/presence of an entrance barrier and the reaction mechanism. The explicit identification of the C2/C3 vs. hydrogen displacement demonstrates that hydrogen-deficient hydrocarbon radicals can be formed easily in environments like those of combustion processes. Our work is a first step towards a systematic database of the intermediates and the reaction products which are involved in this important class of reactions. These findings should be included in future models of PAH and soot formation in combustion flames.

Methylene: a paradigm for computational quantum chemistry.

The year 1970 has been suggested as a starting date for the "third age of quantum chemistry," in which theory takes on not only qualitative but also quantitative value. In fact, each of the years 1960, 1970, 1972, and 1977 is of historical value in the unraveling of the structure and energetics of the CH(2) molecule, methylene. What took place for methylene, namely the establishment of credibility for theory, has subsequently taken place for many other molecules. Three important roles for quantitative theory are outlined: (i) theory precedes experiment; (ii) theory overturns experiment, as resolved by later experiments; and (iii) theory and experiment work together to gain insight that is afforded independently to neither. Several examples from each of the three classes are given.

Theoretical studies of metal-phosphate interactions: interaction of Li+, Na+, K+, Be++, Mg++, and Ca++ with H2PO4- and (CH3O)2PO2-: implications for nucleic acid solvation.

Model phosphate-metal solvation complexes have been studied by ab-initio self-consistent-field techniques. The complexes studied include (RO)2PO2-(R = H or CH3) with Li+, Na+, K+, Be++, Mg++, Ca++, H2O, and Cl-. The geometries of the complexes were chosen to approximate reasonable model solvation complexes for phosphate groups in a nucleic acid environment. Calculated energies of formation vary as Be++ greater than Mg++ greater than Ca++ greater than Li+ greater than Na+ greater than K+ for all isostructural complexes, consistent with experimental binding trends. These results suggest that site binding of this type can successfully account for the relative specificities of ion binding in polynucleotides and other phosphate-containing molecules.

Barrier Height for the Exchange Reaction F + HF --> FH + F.

There exists a body of conflicting data as to the existence or nonexistence of FHF, ClHCl, BrHBr, and IHI as chemically bound molecular species. Ab initio quantum mechanical electronic structure calculations are presented which predict linear symmetric FHF to be unstable. The barrier height for the F + HF exchange reaction is suggested to be no less than 18 kcal/mol, much larger than expected either intuitively or on the basis of certain experiments on related systems. The expected reliability of the calculations is based upon comparable results for diatomic molecules and the F + H(2) and H + F(2) potential energy surfaces.

Potential Energy Surface Including Electron Correlation for F + H2 rarr FH + H: Refined Linear Surface.

A priori quantum mechanical calculations have been carried out at about 150 linear geometries for the fluorine plus hydrogen molecule system. An extended basis set of Gaussian functions was used, and electron correlation was treated explicitly by configuration interaction. Comparison with the experimental activation energy and exothermicity suggests that the theoretical potential surface is quite realistic.