Theoretical study on the spin-state energy splittings and local spin in cationic [Re]-Cn-[Re] complexes

Total spin-state energy splittings are calculated for mono- and dications of the formula {[Re]-Cn-[Re]}z+ where [Re] = eta5-(C5Me5)Re(NO)(PPh3). Cn is an even-numbered carbon chain with n ranging from 4 to 20, and z is 1 or 2. These complexes are experimentally known, and their potential role as molecular electronic devices initiated this work. We have considered the different total spin states monocation/doublet, monocation/quartet, dication/singlet, and dication/triplet. Data obtained for two density functionals BP86 and B3LYP were compared to verify the internal consistency of the results. In both ionization states, the low-spin state is the ground state, but the spin-state splittings decrease as the chain gets longer. For the dications, the splitting reaches a nearly constant value of about 10 kJ/mol with BP86 and about 4 kJ/mol with B3LYP when there are at least 14 carbon atoms in the chain, whereas for the monocations, no constant value appears to be reached asymptotically, not even if 20 carbon atoms are in the chain. For monocations, the splittings range from 138 kJ/mol (n = 4) to 68 kJ/mol (n = 20) with BP86 and from 134 kJ/mol (n = 4) to 73 kJ/mol (n = 20) with B3LYP and are thus considerably higher than those of the dications. The spin-state splittings are qualitatively mirrored by the energy splitting between the highest-occupied molecular orbital with beta spin (HOMObeta) and the lowest-unoccupied molecular orbital with alpha spin (LUMOalpha) as obtained in the low-spin state. Furthermore, the HOMOalpha-LUMOalpha gaps decrease as the carbon chain lengthens. In addition, the local distribution of the ?z expectation value is analyzed for the monocation/doublet, the monocation/quartet, and the dication/triplet state using a modified L?wdin partitioning scheme. In the monocation/doublet and the dication/triplet state, the electron spin is distributed mainly on the metal centers and slightly delocalized onto the carbon chain. In the monocation/quartet state for chain lengths of more than 8 carbon atoms, the electron spin is mainly localized on selected atoms of the chain and not on the metal centers. In all cases, the spin delocalization onto the chain increases as the chain gets longer.     

Dehydrophenylnitrenes: quartet versus doublet states

The geometries and energies of 4-, 3-, and 2-dehydrophenylnitrenes (3, 4, and 5) are investigated using complete active space self-consistent field (CASSCF), multiconfiguration quasi-degenerate second-order perturbation (MCQDPT), and internally contracted multiconfiguration-reference configuration interaction (MRCI) theories in conjunction with a correlation consistent triple-zeta basis set. 4-Dehydrophenylnitrene 3 has a quartet ground state ((4)A(2)). The adiabatic excitation energies to the (2)A(2), (2)B(2), (2)A(1), and (2)B(1) states are 5, 21, 34, and 62 kcal mol(-1), respectively. The (2)B(2) state has pronounced closed-shell carbene/iminyl radical character, while the lowest-energy (2)B(1) state is a combination of a planar allene and a 2-iminylpropa-1,3-diyl. The MCQDPT treatment overestimates the excitation energy to (2)B(2) significantly as compared to CASSCF and MRCI+Q. Among quartet states, (4)A(2)-3 is the most stable one, while those of 4 and 5 (both (4)A') are 3 and 1 kcal mol(-1) higher in energy. 5 also has a quartet ground state and a (2)A' ' state 7 kcal mol(-1) higher in energy. On the other hand, the doublet-quartet energy splitting is -6 kcal mol(-1) for 4 in favor of the doublet state ((2)A'). Hence, (2)A'-4 is the most stable dehydrophenylnitrene, 3.5 kcal mol(-1) below (4)A(2) of 3. The geometry of (2)A'-4 shows the characteristic features of through-bond interaction between the in-plane molecular orbitals at N and at C3. The (2)A' state of 4 resembles the (2)A(1) state of 3 and lies 32 kcal mol(-1) above (4)A'-4. The lowest-energy (2)A' state of 5, on the other hand, resembles the (2)B(2) state of 3 and lies 22 kcal mol(-1) above (4)A'-5.     

Comparative ab initio study of the structures and stabilities of the ethane dication C2H6(2+) and its silicon analogues Si2H6(2+) and CSiH6(2+)

Ab initio MP2/6-311G and QCISD(T)/6-311G levels as well as Gaussian-2 theory were used to perform a comparative study of the structures and stabilities of the ethane dication C(2)H(6)(2+) and its silicon analogues Si(2)H(6)(2+) and CSiH(6)(2+). Similar to previous HF/6-31G results, our present calculations also indicate that the two-electron three-center (2e-3c) bonded carbonium-carbenium structure 1 is more stable than the doubly hydrogen bridged diborane-type structure 2 by about 12 kcal/mol. For the silicon analogue Si(2)H(6)(2+) the calculations, however, indicate that the 2e-3c bonded siliconium-silicenium structure 8 is about 9 kcal/mol less stable than doubly hydrogen bridged structure 9. Similar results were also computed for carbon-silicon mixed CSiH(6)(2+) dication structures. These studies are in agreement with the more electropositive character of silicon compared to carbon. Possible dissociation paths of the minimum structures were also calculated.     

Gas-phase ion/molecule reactions in C5F8

Gas-phase ion-molecule reactions in octafluorocyclopentene (C5F8) were studied with a pulsed electron beam mass spectrometer. When a few Torr of major gas, CH4, Ar, or N2, containing approximately 10 mTorr C5F8 was ionized by 2 keV electrons, C5F8+, C5F7+, C4F6+, C4F5+, and C3F3+ were formed as major fragment ions. The interaction between those ions and C5F8 is found to be a weak electrostatic interaction. The cation...C5F8 bonding energies are around 10 kcal/mol, which were reproduced well by (U)B3LYP/6-311+G(d) calculations. The proton affinity of C5F8 (=148.6 kcal/mol by B3LYP/6-311+G(d)) was found to be smaller than that of C2H4 (=162.8 kcal/mol). In the negative mode of operation, the intense signal of C5F8- was observed during the electron pulse. This indicates that C5F8 has a positive electron affinity (1.27 eV by (U)B3LYP/6-311+G(d)). The C5F8- ion was quickly converted to a complex C10F16-. This complex did not react further with C5F8 down to 170 K. The theoretical calculation revealed that a C5F7-F-...C5F8 interaction mode in (C5F8)2- was converted to a C5F7*...C5F9- one via fluoride-ion transfer. The F- ion was found to form a strong covalent bond with C5F8, but the interaction in F-(C5F8)- - -C5F8 is a weak electrostatic interaction due to the charge dispersal in F-(C5F8). The halide ions except F- interact with C5F8 only weakly. Thermochemical stabilities for the cluster ions I-(CH3I)n (n = 1, 2) were also determined.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.     

C7H12(2+): a prototype hexacoordinate carbonium ion

C(7)H(12)(2+) (1), the prototype hexacoordinate carbonium dication was found to be a viable minimum at the MP2/6-31G** and MP2/cc-pVTZ levels. Structure 1 is a propeller shaped molecule resembling a complex involving a C(2+) with three ethylene molecules resulting in the formation of three two-electron, three-center (2e-3c) bonds. Isomeric structure 2 was found to be 21.8 kcal/mol more stable than structure 1. However, conversion of 1 into 2 through transition structure 3 has a barrier of 5.7 kcal/mol. Related structures 4, 5, and 8 were also located as minima for C(7)H(12)(2+). The isoelectronic boron analogue BC(6)H(12)(+) (10) was also computed to be a minimum at the same level of calculations.     

A comparison of the electrophilic reactivities of Zn2+ and acetic acid as catalysts of enolization: imperatives for enzymatic catalysis of proton transfer at carbon

The deprotonation of the alpha-CH3 and alpha-CH2OD groups of hydroxyacetone and the alpha-CH3 groups of acetone in the presence of acetate buffer and zinc chloride in D2O at 25 degrees C was followed by monitoring the incorporation of deuterium by 1H NMR spectroscopy, and the rate laws for catalysis of these reactions by acetate anion and zinc dication were evaluated. Relative to solvent water at a common standard state of 1 M, Zn2+ provides 6.3 and 4.4 kcal/mol stabilizations, respectively, of the transition states for deprotonation of the alpha-CH2OD and alpha-CH3 groups of hydroxyacetone by acetate anion, and a smaller 3.3 kcal/mol stabilization of the transition state for deprotonation of the alpha-CH3 group of acetone. There is only a 1.4 kcal/mol smaller stabilization of the transition state for the acetate-ion-promoted deprotonation of acetone by the Br?nsted acid acetic acid than by Zn2+, which shows that, in the absence of a chelate effect, there is no large advantage to the use of a metal dication rather than a Br?nsted acid to stabilize the transition state for deprotonation of alpha-carbonyl carbon.     

High-level ab initio studies of unimolecular dissociation of the ground-state N3 radical

A comprehensive study of the unimolecular dissociation of the N(3) radical on the ground doublet and excited quartet potential energy surfaces has been carried out with multireference single and double excitation configuration interaction and second-order multireference perturbation methods. Two forms of the N(3) radical have been located in the linear and cyclic region of the lowest doublet potential energy surface with an isomerization barrier of 62.2 kcal/mol above the linear N(3). Three equivalent C(2v) minima of cyclic N(3) are connected by low barrier, meaning the molecule is free to undergo pseudorotation. The cyclic N(3) is metastable with respect to ground state products, N((4)S)+N(2), and dissociation must occur via intersystem crossing to a quartet potential energy surface. Minima on the seams of crossing between the doublet and quartet potential surfaces are found to lie substantially higher in energy than the cyclic N(3) minima. This strongly suggests that cyclic N(3) possesses a long collision-free lifetime even if formed with substantial internal excitation.     

Experimental and theoretical studies on the thermal decomposition of heterocyclic nitrosimines

A series of substituted 2-nitrosiminobenzothiazolines (2) were synthesized by the nitrosation of the corresponding 2-iminobenzothiazolines (6). Thermal decomposition of 2a--f and of the seleno analogue 7 in methanol and of 3-methyl-2-nitrosobenzothiazoline (2a) in acetonitrile, 1,4-dioxane, and cyclohexane followed first-order kinetics. The activation parameters for thermal deazetization of 2a were measured in cyclohexane (Delta H(++) = 25.3 +/- 0.5 kcal/mol, Delta S(++) = 1.3 +/- 1.5 eu) and in methanol (Delta H(++) = 22.5 +/- 0.7 kcal/mol, Delta S(++) = -12.9 +/- 2.1 eu). These results indicate a unimolecular decomposition and are consistent with a proposed stepwise mechanism involving cyclization of the nitrosimine followed by loss of N(2). The ground-state conformations of the parent nitrosiminothiazoline (9a) and transition states for rotation around the exocyclic C==N bond, electrocyclic ring closure, and loss of N(2) were calculated using ab initio molecular orbital theory at the MP2/6-31G* level. The calculated gas-phase barrier height for the loss of N(2) from 9a (25.2 kcal/mol, MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE) compares favorably with the experimental barrier for 2a of 25.3 kcal/mol in cyclohexane. The potential energy surface is unusual; the rotational transition state 9a-rot-ts connects directly to the orthogonal transition state for ring-closure 9aTS. The decoupling of rotational and pseudopericyclic bond-forming transition states is contrasted with the single pericyclic transition state (15TS) for the electrocyclic ring-opening of oxetene (15) to acrolein (16). For comparison, the calculated homolytic strength of the N--NO bond is 40.0 kcal/mol (MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE).     

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.     

What is the structure of FeC5H6?

Calculations performed using both density functional theory (DFT) and the modified coupled-pair functional (MCPF) approach show that FeC₅H₆⁺ is more stable than HFeC₅H₅⁺ by about 10 kcal/mol. The ground state of FeC₅H₆⁺ is a quartet state derived from the 3d⁷ occupation of Fe⁺. For HFeC₅H₅⁺, the MCPF approach yields a sextet ground state while DFT yields a quartet; however, these two states are close in energy at both levels of theory.     

Stable ion NMR and GIAO-DFT study of novel cations from 8,16-dicyano[2.2]metacyclophanedienes and from strategically substituted/benzannelated dihydropyrenes: charge-induced tropicity modulation and pi-switching

The dicyanometacyclophanediene 1 is diprotonated at the cyano groups (1H2 2+) in various superacid media. Upon quenching, intact 1 and the ring-closed CPD 2 were obtained in a 3:2 or 3:1 ratio, depending on the superacid system. Compound 2 undergoes ring opening in the superacid to give the ipso-monoprotonated 2H+, which on quenching furnishes 1-cyanopyrene as a major product together with 2 and 1. The dication 3 2+, with strongly deshielded internal methyls, was generated from the epoxyannulene 3. Ketones 4-6 and ester 7 are O/C diprotonated to give paratropic carboxonium-annulenium dications (4H2 2+, 5H2 2+, 6H2 2+, and 7H2 2+, respectively). Ester 8 gives a trication by two-electron oxidation and O-protonation. Conjugated carboxylic acid 9 gives a mixture of two dications by CO and ring protonation. The dibromo derivatives 10 and 11 form carboxonium ions, whereas the monobromo derivative 12 is O/C diprotonated to give an oxonium-annulenium dication. Charge delocalization modes and tropicity in the resulting species are evaluated by NMR and GIAO-DFT. Facile formation of 2 from 1 in quenching experiments indicates that thermal closing can be achieved with the diprotonated dinitrile, without imposing skeletal rearrangement.     

MINDO-SCF-MO-berechnung der cis-trans-isomerisierung von N-methyl- und N-phenylazomethin

MINDO/2-SCF-MO calculations for the ground state properties of N-methyl- and N-phenyl-azomethin have been carried out. The calculated rotation barrier for the methyl group in N-methyl-azomethin was 0·8 kcal/mol, the eclipsed conformation being most stable. The calculated rotation barrier about the CN bond in the protonated methylazomethin was 27·9 kcal/mol. MINDO/1-SCF-MO treatment for the N-inversion barrier of the unprotonated species yielded 13·00 kcal/mol. Similar MINDO/2 calculations for N-phenylazomethin yielded 4·0 kcal/mol for the rotation barrier of the phenyl ring around the CN= bond, the perpendicular conformation of the ring to the CNC plane being most stable. For the corresponding N protonated derivative the value 27·3 kcal/mol was calculated for the rotation barrier around the CN bond. MINDO/1 treatment yielded an inversion barrier of 14·0 kcal/mol for N-phenylazomethin.     

Theoretical study of isomerization and dissociation of acetylene dication in the ground and excited electronic states

Ab initio calculations employing the configuration interaction method including Davidson's corrections for quadruple excitations have been carried out to unravel the dissociation mechanism of acetylene dication in various electronic states and to elucidate ultrafast acetylene-vinylidene isomerization recently observed experimentally. Both in the ground triplet and the lowest singlet electronic states of C2H2(2+) the proton migration barrier is shown to remain high, in the range of 50 kcal/mol. On the other hand, the barrier in the excited 2 3A" and 1 3A' states decreases to about 15 and 34 kcal/mol, respectively, indicating that the ultrafast proton migration is possible in these states, especially, in 2 3A", even at relatively low available vibrational energies. Rice-Ramsperger-Kassel-Marcus calculations of individual reaction-rate constants and product branching ratios indicate that if C2H(2)2+ dissociates from the ground triplet state, the major reaction products should be CCH+(3Sigma-)+H+ followed by CH+(3Pi)+CH+(1Sigma+) and with a minor contribution (approximately 1%) of C2H+(2A1)+C+(2P). In the lowest singlet state, C2H+(2A1)+C+(2P) are the major dissociation products at low available energies when the other channels are closed, whereas at Eint>5 eV, the CCH+(1A')+H+ products have the largest branching ratio, up to 70% and higher, that of CH+(1Sigma+)+CH+(1Sigma+) is in the range of 25%-27%, and the yield of C2H++C+ is only 2%-3%. The calculated product branching ratios at Eint approximately 17 eV are in qualitative agreement with the available experimental data. The appearance thresholds calculated for the CCH++H+, CH++CH+, and C2H++C+ products are 34.25, 35.12, and 34.55 eV. The results of calculations in the presence of strong electric field show that the field can make the vinylidene isomer unstable and the proton elimination spontaneous, but is unlikely to significantly reduce the barrier for the acetylene-vinylidene isomerization and to render the acetylene configuration unstable or metastable with respect to proton migration.     

Structure,stability, and electronic and NMR properties of various oxo- and nitrido-derivatives of [L(Salen)Mn(III)]+, where L = none and imidazole. A density functional study

Structure, stability, and electronic and NMR properties of [(Salen)Mn(III)](+)-derived intermediates/reactants in the epoxidation/amination of unfunctionalized olefins, namely [(Salen)Mn(V)O](+) (1-oxo), [(Salen)Mn(IV)O] (2-oxo), and [(Salen)Mn(V)N] (3), have been studied with the B3LYP density functional method. It has been shown that the (1)A, (3)A, and (5)A states of cationic 1-oxo species are virtually degenerate, while for the neutral 2-oxo species the ground (4)A state lies 6.4 kcal/mol lower than (2)A. In the nitrido species 3, the (1)A state has been shown to be the ground state in agreement with experiment. We have investigated isomerization of 1-oxo and 2-oxo species into unusual [(OSalen)Mn(III)](+) (1-N-oxo and 1-peroxo) and [(OSalen)Mn(II)] (2-N-oxo and 2-peroxo) species, respectively. For cationic species 1, the 1-N-oxo isomers are more stable (by 10-12 kcal/mol) than the 1-oxo isomer and are separated from the latter by 21-22 kcal/mol barriers. On the other hand, 1-peroxo isomers are calculated to be 14-16 kcal/mol higher than the 1-oxo isomer. For neutral species 2, however, both 2-N-oxo and 2-peroxo isomers lie significantly higher in energy than the 2-oxo isomer. It has been shown that coordination of axial imidazole ligand alters relative energies of spin states for 1- and 2-oxo species, destabilizing low-spin states. For singlet states of H(2)Salen, 1-oxo, and 3, we have calculated (1)H, (13)C, (15)N, and (17)O NMR chemical shifts using the gauge-independent-atomic orbital (GIAO) approach.     

Synthesis and characterization of the non-Kekulé, singlet biradicaloid Ar'Ge(micro-NSiMe(3))(2)GeAr' (Ar' = 2,6-Dipp(2)C(6)H(3), Dipp = 2,6-i-Pr(2)C(6)H(3))

Reaction of Ar'GeGeAr' (1) with an excess of Me3SiN3 gives the non-Kekulé, biradicaloid Ar'Ge(mu-NSiMe3)2GeAr' (3, Ar' = 2,6-Dipp2C6H3, Dipp = 2,6-i-Pr2C6H3) which has a planar Ge2N2Si2 array and pyramidal geometry at the germaniums. DFT calculations for the model MeGe(mu-NSiH3)2GeMe indicate no Ge-Ge bonding and a singlet ground state. The calculated energy difference between the optimized singlet and triplet states is 17.51 kcal/mol.     

Force field parameters for sulfates and sulfamates based on ab initio calculations: Extensions of AMBER and CHARMm fields

Ab initio self-consistent field (SCF) Hartree-Fock calculations of sulfates R? O? SO₃(?1) (R = Me, Et, i-Pr) and sulfamates R? NHSO₃(?1) (R = H, Me, Et, i-Pr) were performed at the 4-31G(*S*N) //3-21G(*S*N) basis set levels, where asterisks indicate d functions on sulfur and nitrogen atoms. These standard levels were determined by comparing calculation results with several basis sets up to MP2/6-31G*//6-31G*. Several conformations per compound were studied to obtain molecular geometries, rotational barriers, and potential derived point charges. In methyl sulfate, the rotational barrier around the C? O bond is 1.6 kcal/mol at the MP2 level and 1.4 kcal/mol at the standard level. Its ground state has one of three HCOS torsion angles trans and one of three COSO torsion angles trans. Rotation over 60° around the single O? S bond in the sulfate group costs 2.5 kcal/mol at the MP2 and 2.1 kcal/mol at the standard level. For ethyl sulfate, the calculated rotational barrier in going from the ground state, which has its CCOS torsion angle trans, to the syn-periplanar conformation (CCOS torsion angle cis) is 4.8 kcal/mol. However, a much lower barrier of 0.7 kcal/mol leads to a secondary gauchelike conformation about 0.4 kcal/mol above the ground state, with the CCOS torsion angle at 87.6°. Again, one of the COSO torsion angles is trans in the ground state, and the rotational barrier for a 60° rotation of the sulfate group amounts to 1.8 kcal/mol. For methyl sulfamate, the rotational barriers are 2.5 kcal/mol around the C? N bond and 3.3 kcal/mol around the N? S bond. This is noteworthy because sulfamate itself has a calculated rotational barrier around the N? S bond of only 1.7 kcal/mol. These and other data were used to parameterize the well-known empirical force fields AMBER and CHARMm. When the new fields were tested by means of vibrational frequency calculations at the 6-31G*//6-31G* level for methyl sulfate, sulfamate, and methyl sulfamate ground states, the frequencies compared favorably with the AMBER and CHARMm calculated frequencies. The transferability of the force parameters to β-D -glucose-6-sulfate and isopropyl sulfate appears to be better than to isopropyl sulfamate. © 1995 by John Wiley & Sons, Inc.     

Ring opening of 2,5-didehydrothiophene: structures and rearrangements of C4H2S isomers

Electronic structures and rearrangement pathways of several C4H2S isomers are computationally investigated by methods based on coupled cluster theory and density functional theory. Six singlet C4H2S isomers lie within ca. 30 kcal/mol above butatrienethione (6), the apparent global minimum. Ethynylthioketene (7) lies only 2 kcal/mol higher in energy than cumulene 6. Two open-chain isomers, butadiynylthiol (8) and diethynyl sulfide (9), reside ca. 9 and 24 kcal/mol above 6, respectively. Lying 30 kcal/mol above 6, two cyclic singlet isomers, ethynylthiirene (10) and cyclopropenylidenemethanthione (11), are nearly degenerate in energy. Thiophene-2,5-diyl (12) lies substantially higher in energy than 6 (ca. 45 kcal/mol) and is predicted to rearrange preferentially by C-S bond cleavage, leading to thioketene 7, rather than by C-C bond cleavage, leading to diethynyl sulfide (9; retro-Bergman cyclization). Accurate spectroscopic properties of these C4H2S isomers, as well as an understanding of their rearrangement pathways, should facilitate the detection and characterization of these isomers in the laboratory and the interstellar medium.     

Protonated tert-pentyl dication (C5H122+, isopentane dication)

Structures of the tert-pentyl cation (C(5)H(11)(+)) and its protonated dication (C(5)H(12)(2+), isopentane dication) were studied using ab initio methods at the MP2/cc-pVTZ level. Both C-C and C-H hyperconjugatively stabilized structures 1 and 2 , respectively, were found to be minima on the potential energy surface (PES) of the tert-pentyl cation. Structure 1 was computed to be about as stable as structure 2 (slightly more stable by 0.5 kcal mol(-1)). Inter-conversion between 1 and 2 through transition state 3 has a kinetic barrier of only 1.5 kcal mol(-1). The C-H protonated form (H(3)C)(2)C(+)CH(2)CH(4)(+)4 was found to be the global minimum for the protonated tert-pentyl dication. Charges and (13)C NMR chemical shifts of the dication 4 were calculated and compared to those of monocation 1 to study the effect of the additional charge in the dication.     

Density functional study of the retrocyclization of norbornadiene and norbornene catalyzed by Fe(+)

In the presence of Fe(+) catalyst, the retro Diels-Alder reaction of norbornadiene (NBD) is predicted to be stepwise with an activation barrier of 18.8 kcal/mol, which is 3.1 kcal/mol lower than the concerted retro reaction. For norbornene (NBN), the Fe(+)-catalyzed retro reaction is also calculated to be stepwise with an activation barrier of 24.9 kcal/mol, which is 8.5 kcal/mol lower than the uncatalyzed stepwise reaction but 3.8 kcal/mol higher than the concerted reaction. The intermediates from the NBD and NBN retro Diels-Alder reactions, C(5)H(6)FeC(2)H(2)(+) and C(5)H(6)FeC(2)H(4)(+), are predicted to have low activation barriers for ligand-to-ligand hydrogen transfers (through an iron-hydrido intermediate) to form CpFeC(2)H(3)(+) and CpFeC(2)H(5)(+) and, ultimately, vinyl- and ethyl-substituted cyclopentadiene-iron complexes, respectively. In contrast to FeC(2)H(2)(+) and FeC(2)H(4)(+), the lowest-energy pathways on the C(5)H(6)FeC(2)H(2)(+) and C(5)H(6)FeC(2)H(4)(+) potential energy surfaces involve only one multiplicity (quartet). The C(2)H(2) and C(2)H(4) complexes of CpFe(+) and C(5)H(6)Fe(+) are compared.