The [12]annulene global minimum

A new global minimum for [12]annulene has been computationally located. This mono-trans minimum 5 (CCCCCT) is computed to be 1.5 kcal/mol more stable (CCSD(T)/cc-pVDZ//BHHLYP/6-311+G**) than the known tri-trans isomer 1 (CTCTCT) and 2.4 kcal/mol lower than the di-trans isomer 4 (CCTCCT), for which there is indirect evidence. The barriers for several rearrangements of 5 were all computed to be above 15 kcal/mol, indicating that direct experimental characterization of 5 should be possible. The computed barriers for the dynamic processes (including conformational automerization) coupled with computed 1H NMR shift values should aid in the future characterization of this [12]annulene isomer.     

Dehydro[12]annulenes: structures, energetics, and dynamic processes

Density functional and coupled cluster calculations on neutral monodehydro[12]annulenes (C(12)H(10)) reveal a global minimum that should be kinetically stable. At the CCSD(T)/cc-pVDZ//BHLYP/6-31G* level, the unsymmetrical CTCTC conformer 1a lies at least 3 kcal/mol below all other isomers studied. The two isomers closest in energy to 1a are Mo?bius structure 5a (CCTCC) and all-cis 6a. Isomer 1a can undergo conformational automerization with E(a) = 3.9 kcal/mol, implying that this process would be rapid on the NMR time scale, and computed (1)H NMR parameters (GIAO-B3LYP/6-311+G**//RHF/6-31G*) are presented. Cumulenic dehydro[12]annulene isomers, with 1,2,3-butatriene subunits, were found to be reactive intermediates in the interconversion of different configurations of the alkyne forms. Pathways for configuration change of 1a, and for subsequent rearrangement to biphenyl, were investigated. The 28 kcal/mol overall barrier for the lowest energy pathway connecting 1a to biphenyl suggests that 1a is kinetically stable with respect to valence isomerization.     

Theoretical study on stability and properties of NC2O isomers

The structures, energetics, spectroscopies, and stabilities of the doublet NC(2)O radical are explored at density functional theory and ab initio levels. Nine minimum isomers are located connected by 22 interconversion transition states. At the CCSD(T)/6-311+G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is bent NCCO 1 (0.0 kcal/mol) with (2)A' state followed by bent isomer CNCO 2 (16.7). Two isomers (1 and 2) and another high-lying species CCNO 4 (99.4) with bent structure are considerably stabilized by a barrier of at least 20 kcal/mol. All of the three isomers should be experimentally or astrophysically observable. This result is consistent with their indication of neutralization-reionization mass spectrometry experiments. Also, the calculated spectroscopic properties and bond distances of known NCCO 1 are consistent with recent experimental observations and theoretical studies. The bonding natures of the isomers 1, 2, and 4 are analyzed. Their molecular properties including the heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities are calculated at the higher levels G3//B3LYP, G3(MP2)//B3LYP, QCISD, and CCSD(T) (single-point). Possible formation strategies of the isomers 1, 2, and 4 in laboratory and space are also discussed in detail.     

Dynamic processes in [16]annulene: Möbius bond-shifting routes to configuration change

Density functional and ab initio methods have been used to study the mechanisms for key dynamic processes of the experimentally known S4-symmetric [16]annulene (1a). Using BH&HLYP/6-311+G** and B3LYP/6-311+G**, we located two viable stepwise pathways with computed energy barriers (Ea = 8-10 kcal/mol) for conformational automerization of 1a, in agreement with experimental data. The transition states connecting these conformational minima have M?bius topology and serve as starting points for non-degenerate pi-bond shifting (configuration change) via M?bius aromatic transition states. The key transition state, TS1-2, that connects the two isomers of [16]annulene (CTCTCTCT, 1 --> CTCTTCTT, 2) has an energy, relative to the S4 isomer, that ranged from 6.9 kcal/mol (B3LYP/6-311+G**) to 16.7 kcal/mol (BH&HLYP/6-311+G**), bracketing the experimental barrier. At our best level of theory, CCSD(T)/cc-pVDZ(est), this barrier is 13.7 kcal/mol. Several other M?bius bond-shifting transition states, as well as M?bius topology conformational minima, were found with BH&HLYP energies within 22 kcal/mol of 1a, indicating that many possibilities exist for facile thermal configuration change in [16]annulene. This bond-shifting mechanism and the corresponding low barriers contrast sharply with those observed for cis/trans isomerization in acyclic polyenes, which occurs via singlet diradical transition states. All M?bius bond-shifting transition states located in [16]- and [12]annulene were found to have RHF --> UHF instabilities with the BH&HLYP method but not with B3LYP. This result appears to be an artifact of the BH&HLYP method. These findings support the idea that facile thermal configuration change in [4n]annulenes can be accounted for by mechanisms involving twist-coupled bond shifting.     

Möbius aromaticity in [12]annulene: cis-trans isomerization via twist-coupled bond shifting

Density functional and coupled cluster calculations show that facile thermal configuration change in [12]annulene occurs via a twist-coupled bond-shifting mechanism. The transition state for this process is highly aromatic with M?bius topology. At the CCSD(T)/cc-pVDZ//BH&HLYP/6-311+G** level, the isomerization of tri-trans-[12]annulene 1a (CTCTCT) to its di-trans isomer 2 (CCCTCT) via such a mechanism has a barrier of 18.0 kcal/mol, in good agreement with earlier experiments. Two other aromatic M?bius bond-shifting transition states were located that result in configuration change for other [12]annulene conformers. This mechanism contrasts sharply with diradical configuration change for acyclic polyenes and with planar bond-shifting mechanisms generally assumed for annulenes. This constitutes evidence that neutral M?bius aromatic annulenes play a role in the dynamic processes of neutral [4n]annulenes.     

Theoretical study on structures and stability of C4P isomers

The structures, energetics, spectroscopies, and stabilities of doublet C(4)P isomeric species are explored at the DFT/B3LYP, QCISD, and CCSD(T) (singlet-point) levels. A total of 12 minimum isomers and 27 interconversion transition states are located. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the lowest-lying isomer is a floppy CCCCP 1 (0.0 kcal/mol) mainly featuring a cumulenic structure |C=C=C=C=P*|, which differs much from the analogous C4N radical (|*C-C[triple bond]C-C[triple bond]N|). The quasi-linearity and the low bending mode of 1 are in contrast to the previous prediction. The second energetically followed isomer PC-cCCC 3 (14.9 kcal/mol) possesses a CCC ring-bonded to CP. The two low-lying isomers are separated by a high-energy ring-closure/open transition state (26.5 kcal/mol) and thus are very promising candidates for future laboratory and astrophysical detection. Furthermore, four high-energy isomers, that is, two bent isomers CCPCC 2 (68.4 kcal/mol) and CCPCC 2' (68.5 kcal/mol) and two cagelike species 10 (56.0 kcal/mol) and 11 (67.9 kcal/mol), are also stabilized by considerable barriers. The present work is the first detailed potential energy survey of CnP clusters and can provide useful information for the investigation of larger CnP radicals and for understanding the isomerism of P-doped C vaporization processes.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

12]Annulene radical anions revisited: evaluation of structure assignments based on computed energetic and electron spin resonance data

We report density functional and coupled cluster calculations on numerous monocyclic and bicyclic (CH)12(*-) isomers. At the RCCSD(T)/cc-pVDZ//UB3LYP/6-31+G* level, a nearly planar, bond-equalized radical anion of 1,7-di-trans-[12]annulene (4a(*-)) is lowest in energy; several other isomers and conformations lie within 3 kcal/mol of 4a(*-). RCCSD(T)/AUG-cc-pVDZ//UB3LYP/6-31+G* results place the all-cis isomer 3(*-) slightly below 4a(*-) in energy. Validation studies on the heptalene radical anion, [16]annulene radical anion, and tri-trans-[12]annulene radical anion indicate that electron spin resonance (ESR) hyperfine coupling constants (aH values) computed at the BLYP/EPR-III level on DFT geometries give much better agreement with experimental values than those computed using B3LYP/6-31G*. We were unable to locate any C12H12(*-) isomer that could account for the ESR spectrum previously attributed to a highly twisted structure for the 1,7-di-trans-[12]annulene radical anion. Our computed energetic and ESR data for [12]annulene radical anions and their valence isomers suggest that 4a(*-) may have been made, yet its ESR spectrum was incorrectly assigned to the bicyclic isomer 6b(*-). Finally, the computed (1)H NMR shift values of the dianion of 4 reveal a distinct diatropic ring current that should aid in its characterization.     

Theoretical study on structures and stability of Si2CP isomers

The structures, energetics, spectroscopies, and isomerization of various doublet Si2CP species are explored theoretically. In contrast to the previously studied SiC2N and SiC2P radicals that have linear SiCCN and SiCCP ground states, the title Si2CP radical has a four-membered-ring form cSiSiPC 1 (0.0 kcal/mol) with Si-C cross-bonding as the ground-state isomer at the CCSD(T)/6-311G(2df)//B3LYP/6-311G(d)+ZPVE level, similar to the Si2CN radical. The second low-lying isomer 2 at 11.6 kcal/mol has a SiCSiP four-membered ring with C-P cross-bonding, yet it is kinetically quite unstable toward conversion to 1 with a barrier of 3.5 kcal/mol. In addition, three cyclic species with divalent carbene character, i.e., cSiSiCP 7, 7' with C-P cross-bonding and cSiCSiP 8 with Si-Si cross-bonding, are found to possess considerable kinetic stability, although they are energetically high lying at 44.4, 46.5, and 41.4 kcal/mol, respectively. Moreover, a linear isomer SiCSiP 5 at 44.3 kcal/mol also has considerable kinetic stability and predominantly features the interesting cumulenic /Si=C=Si=P/* form with a slight contribution from the silicon-phosphorus triply bonded form /Si=C*-Si[triple bond]P/. The silicon-carbon triply bonded form *Si[triple bond]C-Si[triple bond]P/ has negligible contribution. All five isomers are expected to be observable in low-temperature environments. Their bonding nature and possible formation strategies are discussed. For relevant species, the QCISD/6-311G(d) and CCSD(T)/6-311+G(2df) (single-point) calculations are performed to provide more reliable results. The calculated results are compared to those of the analogous C3N, C3P, SiC2N, and Si2CN radicals with 17 valence electrons. Implications in interstellar space and P-doped SiC vaporization processes are also discussed.     

SiCNN--a new stable isomer with Si(triple bond)C triple bonding

To predict potentially stable molecules with Si(triple bond)C triple bonding, theoretical calculations at the B3LYP/ 6-311G(d) and CCSD(T)/6-311G(2df) (single-point) levels were employed to study the structures, energetics, and isomerization of various SiCN2 isomers. A schematic potential energy surface (PES) of SiCN2 was established to discuss the kinetic stability of the isomers. A new isomer SiCNN was found to possess a typical Si(triple bond)C triple bond, as confirmed by comparative calculations at the B3LYP, QCISD, QCISD(T), CCSD, and CCSD(T) levels on the bond lengths of SiCNN and other experimentally or theoretically known species of RSiCH (R = H, F, Cl, OH). Moreover, SiCNN resides in a very deep potential, the stabilization barrier is at least 53.2 kcal mol(-1). Thus, SiCNN may be considered as the most kinetically stable isomer with Si(triple bond)C triple bonding known to date, and it may represent a very promising molecule for future experimental characterization. In addition, the stability of the other isomers, such as the four linear species SiNCN, SiNNC, NSiCN and NSiNC, a three-membered NNC ring isomer with exocyclic C-Si bonding, and a four-membered SiCNN ring isomer is discussed and compared with SiCNN.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

Theoretical study on the isomerization behavior between alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals

[reaction: see text] Ab initio calculations using 6-311G**, cc-pVDZ, aug-cc-pVDZ, and a (valence) double-zeta pseudopotential (DZP) basis set, with (QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional methods (BHandHLYP, B3LYP) predict that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals exist as isomers. At the CCSD(T)/cc-pVDZ//BHandHLY/cc-pVDZ level of theory, energy barriers of 15.1 and 17.7-21.7 kJ mol(-)(1) are calculated for the isomerization of s-trans-propenoyl and s-trans-crotonoyl radical to ketenylmethyl and 1-ketenylethyl radical, respectively. Similar results are obtained for the reactions of s-trans isomers involving silyl, germyl, and stannyl groups with energy barriers (DeltaE++) of 12.2-12.4, 13.1-13.9, and 12.9-18.2 kJ mol(-)(1) at the CCSD(T)/DZP//BHandHLYP/DZP calculation, respectively. These results suggest that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals are not canonical forms but are isomeric species that can rapidly interconvert.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

Computational evaluation of the evidence for tri-trans-[12]annulene

[reaction: see text] Automerization in tri-trans-[12]annulene (1) was investigated by DFT, MP2, and coupled-cluster methods. Using the highest level of theory employed here, CCSD(T)/cc-pVDZ//BHandHLYP/6-311+G(d,p), we located two low-energy pathways for degenerate conformational change from the lowest-energy conformer of 1 (1a): one with E(a) = 4.5 kcal/mol that interconverts the three inner trans hydrogens with the three outer trans hydrogens and one with E(a) = 2.7 kcal/mol that interconverts the three inner hydrogens with each other. These results are consistent with the experimental results of Oth and co-workers on [12]annulene 1a (Oth, J. F. M.; R?ttele, H.; Schr?der, G. Tetrahedron Lett. 1970, 61). The conformational exchange of the inner trans hydrogens with the outer ones is predicted to occur via a one-step process involving a C(2)-symmetric transition state and not via the D(3)-symmetric transition state (1b) that was postulated earlier. Conformer 1b was found to be a shallow minimum 6.7 kcal/mol above 1a with a barrier of 0.4 kcal/mol for conversion to 1a. Finally, GIAO-B3LYP/6-311+G(d,p) and BHandHLYP/6-311+G(d,p) computed (1)H NMR chemical shifts of 1a and three other low-lying isomers support Oth's original assignment of observed (1)H NMR peaks to 1a at both low and high temperature.     

Determining the energy gap between the cis and trans isomers of HO3- using geometry optimization within the anti-Hermitian contracted Schrödinger and coupled cluster methods

The cis and trans isomers of the HO3- anion, which are important in proposed mechanisms for ozonization, are studied computationally. Relative energies, geometries, and normal-mode frequencies are calculated with anti-Hermitian contracted Schr?dinger equation (ACSE) and coupled cluster methods. Both the ACSE method and the coupled cluster method with single and double excitations (CCSD) are applied in a correlation-consistent polarized double-zeta basis set (cc-pVDZ). Using coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), we treat the problem with larger basis sets than those in previous work, including correlation-consistent polarized quadruple-zeta basis sets with (aug-cc-pVQZ) and without (cc-pVQZ) diffuse functions, which permit extrapolation of the cis and trans energies to the complete-basis-set limit. The cis isomer is found to be lower in energy than the trans isomer by -3.5 kcal/mol, which is 50% larger in magnitude than the best previous result of -2.2 kcal/mol. The bond lengths between the O2 and OH fragments of the cis- and trans-HO3 are calculated to be 1.713 and 1.857 A, respectively, where both bond lengths are significantly longer than the 1.464 A O-O bond in hydrogen peroxide. In this paper, we extend the ACSE method [Mazziotti, D. A. J. Chem. Phys. 2007, 126, 184101], which computes the two-electron reduced density matrix directly, to include geometry optimization by a Newton's method with numerical derivatives. Calculation of the cis- and trans-HO3- isomers by the ACSE yields energies, geometries, and frequencies that are closer to those from CCSD(T) than those from CCSD.     

Carbene Rearrangements Unsurpassed: Details of the C(7)H(6) Potential Energy Surface Revealed

The rearrangement of phenylcarbene (1) to 1,2,4,6-cycloheptatetraene (3) has been studied theoretically, using SCF, CASSCF, CASPT2N, DFT (B3LYP), CISD, CCSD, and CCSD(T) methods in conjunction with the 6-31G, 6-311+G, 6-311G(2d,p), cc-pVDZ, and DZd basis sets. Stationary points were characterized by vibrational frequency analyses at CASSCF(8,8)/6-31G and B3LYP/6-31G. Phenylcarbene (1) has a triplet ground state ((3)A") with a singlet-triplet separation (DeltaE(ST)) of 3-5 kcal mol(-)(1). In agreement with experiment, chiral 3 is the lowest lying structure on this part of the C(7)H(6) potential energy surface. Bicyclo[4.1.0]hepta-2,4,6-triene (2) is an intermediate in the rearrangement of 1 into 3, but it is unlikely to be observable experimentally due to a barrier height of only 1-2 kcal mol(-)(1). The enantiomers of 3 interconvert via the (1)A(2) state of cycloheptatrienylidene (4) with an activation energy of 20 kcal mol(-)(1). The "aromatic" (1)A(1) state, previously believed to be the lowest singlet state of 4, is roughly 10 kcal mol(-)(1) higher in energy than the (1)A(2) state, and, in violation of Hund's rule, (3)A(2) is also calculated to lie above (1)A(2) in energy. Thus, even if (3)A(2) were populated, it is likely to undergo rapid intersystem crossing to (1)A(2). We suggest (3)B(1)-4 is the metastable triplet observed by EPR.     

Thinking out of the black box: accurate barrier heights of 1,3-dipolar cycloadditions of ozone with acetylene and ethylene

Accurate barriers for the 1,3-dipolar cycloadditions of ozone with acetylene and ethylene have been determined via the systematic extrapolation of ab initio energies within the focal point approach of Allen and co-workers. Electron correlation has been accounted for primarily via coupled cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [CCSD, CCSD(T), CCSDT, and CCSDT(Q)]. For the concerted [4 + 2] cycloadditions, the final recommended barriers are DeltaH(0K) = 9.4 +/- 0.2 and 5.3 +/- 0.2 kcal mol(-1) for ozone adding to acetylene and ethylene, respectively. These agree with recent results of Cremer et al. and Anglada et al., respectively. The reaction energy for O3 + C2H2 exhibits a protracted convergence with respect to inclusion of electron correlation, with the CCSDT/cc-pVDZ and CCSDT(Q)/cc-pVDZ values differing by 2.3 kcal mol-1. Recommended enthalpies of formation (298 K) for cycloadducts 1,2,3-trioxole and 1,2,3-trioxolane are +32.8 and -1.6 kcal mol(-1), respectively. Popular composite ab initio approaches [CBS-QB3, CBS-APNO, G3, G3B3, G3(MP2)B3, G4, G4(MP3), and G4(MP2)] predict a range of barrier heights for these systems. The CBS-QB3 computed barrier for ozone and acetylene, DeltaH(0K) = 4.4 kcal mol(-1), deviates by 5 kcal mol(-1) from the focal point value. CBS-QB3 similarly underestimates the barrier for the reaction of ozone and ethylene, yielding a prediction of only 0.7 kcal mol(-1). The errors in the CBS-QB3 results are significantly larger than mean errors observed in application to the G2 test set. The problem is traced to the nontransferability of MP2 basis set effects in the case of these reaction barriers. The recently published G4 and G4(MP2) approaches perform substantially better for O3 + C2H2, predicting enthalpy barriers of 9.0 and 8.4 kcal mol(-1), respectively. For the prediction of these reaction barriers, the additive corrections applied in the majority of the composite approaches considered lead to worse agreement with the reference focal point values than would be obtained relying only on single point energies evaluated at the highest level of theory utilized within each composite method.     

Structure and Stability of Isomers of the Promising Interstellar Molecule PC3O

DFT/B3LYP/6-311G(d) and CCSD(T)/6-311G(2d) single-point calculations are carried out for exploring the doublet potential energy surface (PES) of PC₃O, a molecule of potential interest in interstellar chemistry. A total of 29 minima connected by 65 interconversion transition states are located. The structures of the most relevant isomers and transition states are further optimized at the QCISD level followed by CCSD(T) single-point energy calculations. At the CCSD(T)/6-311G(2df)//QCISD/6-311G(d)+ZPVE level, the global minimum is the quasi-linear structure PCCCO 1 (0.0 kcal/mol) with a great kinetic stability of 47.9 kcal/mol, and the cumulenic form features largely in its resonance structures. Moreover, the chainlike isomer OPCCC 3 (64.5) and five-membered-ring species cPCCCO 19 (77.8) possess considerable kinetic stability of about 18.0 kcal/mol. All these three isomers are very promising candidates for future experimental and astrophysical detection. Additionally, a three-membered-ring isomer CC-cCOP 10 (69.6) has slightly lower kinetic stability of around 15 kcal/mol and may also be experimentally observable. Possible formation mechanisms of the four stable isomers in interstellar space are discussed. The present research is the first attempt to study the isomerization and dissociation mechanisms of PC _n O series. The predicted spectroscopic properties, including harmonic vibrational frequencies, dipole moments and rotational constants for the relevant isomers, are expected to be informative for the identification of PC₃O in laboratory and interstellar medium.     

A theoretical survey on the structures, energetics, and isomerization pathways of the B(5)O radical

Boron-centered radicals have received growing interest. Recently, two groups reported density functional theory investigations (GGA-PW91 and B3LYP) on a hexa-atomic boron-oxide radical, B(5)O, which has led to great discrepancies on the type of low-lying structures. In this article, we not only explore the energetics of doublet and quartet B(5)O isomers at high electron-correlated levels (CCSD(T)/6-311+G(2df), CCSD(T)/aug-cc-pVTZ, and G3B3) but also investigate the isomerization and fragmentation stability of the low-lying B(5)O isomers. All the high-level studies consistently show that the B(5)O radical possesses a belt-like ground structure (2)01 in doublet electronic state followed by isomer (2)02 with an exocyclic - BO moiety at around 3.0 kcal/mol. Kinetically, (2)01 and (2)02 are separated by a considerable barrier of about 20 kcal/mol. Thus, the two isomeric forms of B(5)O radical should be very promising for isolation in laboratory. However, the other four isomers reported recently are all kinetically unstable toward conversion to (2)01 and (2)02. The high thermodynamic and kinetic stability of (2)01 and (2)02 might make them as important building cores in the growth of boron-oxide clusters. This results would also help deeply understand the oxidation and doping mechanism of pure boron clusters.     

Ab initio studies of isomerization and dissociation reactions of methyl peroxynitrate

Using the CCSD(T)/cc-pVDZ//B3LYP/6-311G(2d,2p) method, we calculated the detailed potential energy surfaces (PESs) for the unimolecular isomerization and decomposition of methyl peroxynitrate (CH₃O₂NO₂). The results show that there are the two most stable isomers, IS1a and IS1b, which are a pair of mirror image isomers. From IS1a and IS1b, different isomerization and unimolecular decomposition reaction channels have been studied and discussed. Among them, the predominant thermal decomposition pathways are those leading to CH₃O₂ + NO₂ and cis-CH₃ONO + O₂. The former is the lowest-energy path through the direct O–N bond rupture in IS1a or IS1b. The PES along the O–N bond in IS1a has been scanned, where the energy of IS1a reaches maximum value of 23.5 kcal/mol when the O–N bond is stretched to about 2.8 Å. This energy is 2.7 kcal/mol larger than the O–N bond dissociation energy (BDE) and 2.8 kcal/mol larger than the experimental active energy. In addition, because the energy barriers of IS1a isomerization to IS2a are 23.8 kcal/mol, close to the 20.8 kcal/mol O–N BDE in IS1a or IS1b, the isomerization reaction may compete with the direct bond rupture dissociation reaction.