Theoretical study on structures and stabilities of [H,Ge,C,N

Theoretical investigations are performed for the first time on the simplest hydrogenated germanium cyanide [H,Ge,C,N], whose analogs [H,C(2),N] and [H,Si,C,N] have been detected in space and laboratory, respectively. The detailed potential energy surfaces in both singlet and triplet states are constructed at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-31G(d)+ZPVE level, including 18 minimum isomers and 26 interconversion transition states. The former three low-lying and kinetically stabilized isomers are HGeCN (1)1 (0.0 kcal/mol), HGeNC (1)2 (5.1 kcal/mol), and cyclic cCHNGe(1)7 (11.1 kcal/mol). In addition, five isomers HCNGe (1)3 (33.8), HNCGe (1)5 (29.8), cNHCGe (1)8 (37.9), HGeCN (3)1 (30.1), and HNCGe (3)5 (26.5) each have considerable barriers, despite their high energies. Future laboratory characterization and astrophysical detection of the eight [H,Ge,C,N] isomers, especially the former three low-lying species (1)1, (1)2, and (1)7, are highly recommended. The accurate spectroscopic data at the QCISD/6-311G(d,p) level are provided. For some species, the CBS-QB3 calculations are also performed. Wherever possible, comparisons with the analogous [H,C(2),N] and [H,Si,C,N] are made on the structural, energetic, and bonding properties.     

Monosilicon-substituted cyanoacetylene: a computational study

A detailed theoretical investigation of the [H,Si,C(2),N] potential energy surfaces including 28 minimum isomers and 65 interconversion transition states is reported at the Gaussian-3//B3LYP/6-31G(d) level. Generally, the triplet species lie energetically higher than the singlet ones. The former three low-lying isomers are linear HCCNSi 1 (0.00 kcal/mol), branched SiC(H)CN 12 (7.09 kcal/mol), and bent HNCCSi 7 (14.22 kcal/mol), which are separated by rather high barriers from each other and are kinetically very stable with the least conversion barriers of 32.6-70.5 kcal/mol. Two energetically high-lying isomers HCNCSi 3 (42.99 kcal/mol) and SiC(H)NC 13 (36.05 kcal/mol) are also kinetically stable with a barrier of 49.19 and 21.42 kcal/mol, respectively. Additionally, five high-lying isomers, that is, three chainlike isomers, HCCSiN 2 (55.17), HCSiNC 6 (47.80), HSiNCC 11 (78.83), and one three-membered ring isomer HN-cSiCC 19 (51.21), and one four-membered ring isomer cSiCN(H)C 27 (50.6 kcal/mol), are predicted to each have lower conversion barriers of 12-18 kcal/mol and can be considered as meta-stable species. All of the predicted 10 isomers could exist as stable or meta-stable intermediates under suitable conditions. Finally, the structural and bonding analysis indicate that the [H,Si,C(2),N] molecule contains various properties that are of chemical interest (e.g., silylene, SiC triple bonding, and conjugate SiN triple bonding and CC triple bonding, charge-transfer specie, planar aromatic specie, cumulate double bonding). This is the first detailed theoretical study on the potential energy surfaces of the series of hydrogenated Si,C,C,N-containing molecules. The knowledge of the present monohydrogenated SiC(2)N isomerism could provide useful information for more highly hydrogenated or larger Si,C(2),N-containing species.     

Dicyclobuta[de,ij]naphthalene and dicyclopenta[cd,gh]pentalene: a theoretical study

The structures, energetics, and aromatic character of dicyclobuta[de,ij]naphthalene, 1, dicyclopenta[cd,gh]pentalene, 2, dihydrodicyclobuta[de,ij]naphthalene, 3, and dihydrocyclopenta[cd,gh]pentalene, 4, have been examined at the B3LYP/6-311++G//B3LYP/6-31G level of theory. All molecules are bowl-shaped, and the pentalene isomers, 2 and 4, are most stable. A comparison with other C(12)H(6) and C(12)H(8) isomers indicates that 2 is approximately 25 kcal/mol less stable than 1,5,9-tridehydro[12]annulene and 4 is approximately 100 kcal/mol higher in energy than acenaphthylene, both of which are synthetically accessible. The transition state structure for bowl-to-bowl inversion of 1 is planar (D(2)(h)()) and lies 30.9 kcal/mol higher in energy than the ground state; the transition state for inversion of 2 is C(2)(h)() and lies 46.6 kcal/mol higher in energy. Symmetry considerations, bond length alternations, and NICS values (a magnetic criterion) all indicate that the ground states of 1, 3, and 4 are very aromatic; however, HOMA values (a measure of bond delocalization) indicate that 3S and 4S are aromatic but that 1S is less so. NICS values for the ground state of 2 strongly indicate aromaticity; however, bond localization, symmetry, and HOMA values argue otherwise.     

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.     

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.     

An Ab Initio Molecular Orbital Study of Isophosphinines

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 phosphinine and 13 isophosphinines 7-19 . Isomers 7-11 with an allenic system are calculated to be 8-18 kcal mol m 1 more stable than structures 12-17 with an acetylenic moiety. The calculated energy difference (66.19 kcal mol m 1 ) between phosphinine and the most stable isophosphinine (1-phospha-1,2,4-cyclohexatriene, 10 ) is smaller than the difference (78.96 kcal mol m 1 ) between benzene and the most stable isobenzene (cyclohexa-1,2,4-triene, 2 ). The isophosphinines 18 and 19 , with a butatriene moiety, are calculated to be the least stable isomers.     

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.     

A combined density functional theory and coupled cluster method investigation of the structural properties and stabilities of radical CH(2)CP and its isomers

The doublet potential energy surface of radical system [C(2), H(2), P] is investigated at the UB3LYP/6-311++G(d,p) and UCCSD(T)/6-311++G(2df,2p) (single-point) levels. Eight chainlike and three-membered ring structures are located as energy minima connected by 10 interconversion transition states. At the final UCCSD(T)/6-311++G(2df,2p)//UB3LYP/6-311++G(d,p) level with zero-point vibrational energy correction, species CH(2)CP is found to be thermodynamically the most stable isomer followed by HCCPH, H-cCPC-H, cPCC-H(H), H-cCCP-H, cis-CC(H)PH, trans-CC(H)PH, and CCPH(2) at 11.01, 12.57, 40.07, 43.63, 50.25, 56.82, and 65.36 kcal/mol, respectively. The computed results indicate that the chainlike isomers CH(2)CP and HCCPH and cyclic radical H-cCPC-H possess considerable kinetic stability at extra low pressures and temperatures. Interestingly, radical CCPH(2), whose energy is the highest in all predicted CH(2)CP isomers, can be also regarded as a kinetically stable species with the smallest isomerization barrier of 22.26 kcal/mol at extra low pressures and temperatures. Therefore, considering higher kinetic stability, in addition to the microwave spectroscopy characterized isomer CH(2)CP in previous experiments, the species HCCPH, H-cCPC-H, and CCPH(2) should be considered as excellent candidates for possible experimental observation. Furthermore, the structural nature of stable radical isomers is discussed based on bonding characteristics, single electron spin distribution, and comparison with their analogues.     

Internal degrees of freedom,structural motifs,and conformational energetics of the 5'-deoxyadenosyl radical: implications for function in adenosylcobalamin-dependent enzymes. A computational study

The potential energy surface of the free 5'-deoxyadenosyl radical in the gas phase is explored using density functional and second-order M?ller-Plesset perturbation theories with 6-31G(d) and 6-31++G(d,p) basis sets and interpreted in terms of attractive and repulsive interactions. The 5',8-cyclization is found to be exothermic by approximately 20 kcal/mol but kinetically unfavorable; the lowest cyclization transition state (TS) lies about 7 kcal/mol higher than the highest TS for conversion between most of the open isomers. In open isomers, the two energetically most important attractive interactions are the hydrogen bonds (a) between the 2'-OH group and the N3 adenine center and (b) between the 2'-OH and 3'-OH groups. The relative ribose-adenine rotation about the C1'-N9 glycosyl bond in a certain range changes the energy by as much as 10-15 kcal/mol, the origin being (i) the repulsive 2'-H.H-C8 and O1'.N3 and (ii) the attractive 2'-OH.N3 ribose-adenine interactions. The hypothetical synergy between the glycosyl rotation and the Co-C bond scission may contribute to the experimentally established labilization of the Co-C bond in enzyme-bound adenosylcobalamin. The computational results are not inconsistent with the rotation about the C1'-N9 glycosyl bond being the principal coordinate for long-range radical migration in coenzyme B(12)-dependent enzymes. The effect of the protein environment on the model system results reported here remains an open question.     

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 .     

Bicyclo[3.2.1]octa-3,6-dien-2-yl cation: a bishomoantiaromate

Antiaromatic compounds with a closed loop of 4n p-electrons are relatively unstable and often difficult to study. We report in this article the synthesis of alcohols 2-(4'-fluorophenyl)bicyclo[3.2.1]octan-2-ol 11, 2-(4'-fluorophenyl)bicyclo[3.2.1]oct-3-en-2-ol 12, and 2-(4'-fluorophenyl)bicyclo[3.2.1]octa-3,6-dien-2-ol 13 and their transformations into corresponding carbocations 14-16, respectively, in a superacidic medium (FSO3H/SO2ClF) at -120 degrees C. Cations 14-16 are characterized by NMR analysis (1H, 13C, 19F), and 15 and 16 are further characterized by quenching in NaOCH3/H3COH at -120 degrees C. The relative stabilities of 14-16 are determined experimentally by 19F NMR spectroscopy. Cation 16 is found to be experimentally less stable than cation 15 by 3.7 kcal/mol. DFT calculations (structure and energy: B3LYP/6-31G(d); NMR: B3LYP/6-311+G(2d,p)) are performed for alcohols 11-13 and bicyclo[3.2.1]octyl cations 6, 7, 9, 14-16, 26, 28, and 30. In the case of 11-16, data from DFT calculations is in good agreement with experimental data. Because 6,7-dimethylenebicyclo[3.2.1]oct-3-en-2-yl cation 26 is more stable than cation 7 by 1.69 kcal/mol, the inductive effect of sp(2)-hybridized carbon atoms C6 and C7 in carbocations 6 and 16 cannot be the reason for the destabilization of 6 relative to 7 and 16 relative to 15. Destabilization of 6 relative to 7 and 16 relative to 15 and the calculated NICS of 6 (+4.17 ppm) and 16(+3.3 ppm) document that 6 and 16 are bishomoantiaromates.     

10]Annulene: bond shifting and conformational mechanisms for automerization

We report density-functional and coupled-cluster calculations on conformation change and degenerate bond shifting in [10]annulene isomers 1-5. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, conversion of the twist (1) to the heart (2) has a barrier of 10.1 kcal/mol, compared to Ea = 16.2 kcal/mol for degenerate "two-twist" bond shifting in 1. Pseudorotation in the all-cis boat isomer (3) proceeds with a negligible barrier. The naphthalene-like isomer 4 has a 3.9 kcal/mol barrier to degenerate bond shifting. The azulene-like isomer 5 is the only species for which the nature of the bond-equalized form (5-eq) depends on the method. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, 5-eq is 1.2 kcal/mol more stable than the bond-alternating form 5-alt. Conversion of 5-eq to 4 has a barrier of 12.6 kcal/mol. Despite being significantly nonplanar, both 5-eq and the transition state for bond shifting in 4 are highly aromatic based on magnetic susceptibility exaltations. On the basis of a detailed consideration of these mechanisms and barriers, we can now, with greater confidence, rule out 4 and 5 as candidates to explain the NMR spectra observed by Masamune. Our results support Masamune's original assignments for both isolated isomers.     

Dimers of and tautomerism between 2-pyrimidinethiol and 2(1H)-pyrimidinethione: a density functional theory (DFT) study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.     

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.     

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.     

Favorable pendant-amino metal chelation in VX nerve agent model systems

We have performed DFT computational studies [B3LYP, 6-31+G] to obtain metal ion coordination isomers of VX-Me [MeP(O)(OMe)(SCH2CH2NMe2)], a model of two of the most lethal nerve agents: VX [MeP(O)(OEt)(SCH2CH2N(iPr)2)] and Russian-VX [MeP(O)(OCH2CHMe2)(SCH2CH2N(Et)2)]. Our calculations involved geometry optimizations of the neutral VX-Me model as well as complexes with H+, Li+, Na+, K+, Be2+, Mg2+, and Ca2+ that yielded 2-8 different stable chelation modes for each ion that involved mainly mono- and bidentate binding. Importantly, our studies revealed that the [O(P),N] bidentate binding mode, long thought to be the active mode in differentiating the hydrolytic path of VX from other nerve agents, was the most stable for all ions studied here. Binding energy depended mainly on ionic size as well as charge, with binding energies ranging from 364 kcal mol(-1) for Be2+ to 33 kcal mol(-1) for K+. Furthermore, calculated NMR shifts for VX-Me correlate to experimental values of VX.     

Theoretical study on the structures, isomerization and stability of SiC3H isomers

The calculations of the geometry optimizations, energies, dipole moments, vibrational spectra, rotational constants, and isomerization of doublet SiC₃H species were performed using density functional theory and ab initio methods. Four types of isomers, a total of 18 minima, connected by 16 interconversion transition states, were located on the potential energy surface (PES) at the B3LYP/6-311G (d, p) level. More accurate energies were obtained at the CCSD(T)/6-311G(2df, 2p), and G3(MP2) levels. With the highest isomerization barrier, the lowest lying structure, linear A1 possesses the largest kinetic stability. Besides, the isomerization barriers of A2, A4, C2, F1, F4 and F5 are over 10 kcal/mol, and these isomers are also considered to be higher kinetically stable. Other isomers cannot be kinetically stabilized with considerably low isomerization barriers. Investigation on the bonding properties and the computations of vibrational spectra, dipole moments, and rotational constants for SiC₃H isomers are helpful for understanding their structures and also valuable for their detections in the interstellar space and laboratory.     

Energetic preferences for alpha,beta versus beta,gamma unsaturation

Density functional theory has been applied at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level to examine the energetics of alpha,beta- versus beta,gamma-unsaturation for some common organic functional groups. Specifically, the relative stabilities of allyl-X (H2C=CHCH2X) and 1-propenyl-X (H3CCH=CHX) isomers have been computed for X = methyl, vinyl, phenyl, formyl, acetyl, methoxy, methylthio, methylsulfinyl, methylsulfonyl, sulfamoyl, and methoxysulfonyl, and the results are compared to available experimental data. The intrinsic preference of 3 kcal/mol for the 1-propenyl isomer when X = CH3 is exceeded by 2-4 kcal/mol for first-row conjugating groups, but it is not met for the sulfur-containing groups. In particular, alpha,beta-unsaturation is favored by less than 1 kcal/mol for the sulfone and sulfonamide analogues, while it is preferred by 8 kcal/mol for the vinyl-substituted case. Detailed structural results and torsional energy profiles are also reported.     

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.     

Theoretical investigation on the protonation reactions and products of the stable [N,C,C,S] isomers

A theoretical study on the protonation system of [N,C,C,S], [H,N,C,C,S]+, was performed at the B3LYP/6-311++G(d,p) and CCSD(T)/6-311++G(2df,2p) (single point) levels of theory. On the doublet [H,N,C,C,S]+ surface, 24 species were located as energy minima and 10 of them were considered as kinetically stable species. The species HNCCS+ with 2A' state and a shallow W-shaped skeleton was predicted to be the global minimum and kinetically the most stable species, being in good agreement with previous experimental findings. Furthermore, the protonation reactions of the stable [N,C,C,S] isomers were investigated in detail. The calculation results indicated that the [N,C,C,S] isomers may be significantly stabilized upon protonation. Finally, the possible covalent structures of the [H,N,C,C,S]+ isomers with considerable stability were briefly discussed.