Electrophilic attack on sulfur-sulfur bonds: coordination of lithium cations to sulfur-rich molecules studied by ab initio MO methods

Complex formation between gaseous Li+ ions and sulfur-containing neutral ligands, such as H2S, Me2Sn (n = 1-5; Me = CH3) and various isomers of hexasulfur (S6), has been studied by ab initio MO calculations at the G3X(MP2) level of theory. Generally, the formation of LiS(n) heterocycles and clusters is preferred in these reactions. The binding energies of the cation in the 29 complexes investigated range from -88 kJ mol(-1) for [H2SLi]+ to -189 kJ mol(-1) for the most stable isomer of [Me2S5Li]+ which contains three-coordinate Li+. Of the various S6 ligands (chair, boat, prism, branched ring, and triplet chain structures), two isomeric complexes containing the S5==S ligand have the highest binding energies (-163+/-1 kJ mol(-1)). However, the global minimum structure of [LiS6]+ is of C(3v) symmetry with the six-membered S(6) homocycle in the well-known chair conformation and three Li--S bonds with a length of 256 pm (binding energy: -134 kJ mol(-1)). Relatively unstable isomers of S6 are stabilized by complex formation with Li+. The interaction between the cation and the S6 ligands is mainly attributed to ion-dipole attraction with a little charge transfer, except in cations containing the six sulfur atoms in the form of separated neutral S2, S3, or S4 units, as in [Li(S3)2]+ and [Li(S2)(S4)]+. In the two most stable isomers of the [LiS6]+ complexes, the number of S--S bonds is at maximum and the coordination number of Li+ is either 3 or 4. A topological analysis of all investigated complexes revealed that the Li--S bonds of lengths below 280 pm are characterized by a maximum electron-density path and closed-shell interaction.     

Isomers of cyclo-heptasulfur and their coordination to Li(+): an ab initio molecular orbital study

The potential energy hypersurfaces (PESs) of heptasulfur (S7) and of [LiS7]+ have been investigated by ab initio molecular orbital calculations at the G3X(MP2) level of theory. Besides the chair-like seven-membered ring (1a) as the global minimum structure, eight local minimum structures and one transition state have been located on the PES of S7. The barrier for pseudorotation of 1a is only 5.6 kJ mol(-1). The boat-like S7 ring (1b) is 12.1 kJ mol(-1) less stable than 1a, followed by three isomers of connectivity S6=S and four open-chain isomers. On the basis of multireference calculations at the MRCI(4,4)+Q/6-311G(d) level, the most stable open-chain form of S7 is a triplet of relative energy 133.1 kJ mol(-1). Thus, the reaction energy (deltaE0) for the ring opening of 1a is 133.1 kJ mol(-1), halfway between those of the highly symmetrical rings S6 and S8. Because of their strong multireference characters, the stabilities of the biradicalic singlet chains are significantly overestimated by the single-reference-based G3X(MP2) method. The calculated vibrational spectrum of 1a is in good agreement with experimental data. The various isomers of S7 form stable complexes with Li+ with coordination numbers of 1-4 for the metal atom and binding energies in the range of -93.8 to -165.7 kJ mol(-1). A total of 15 isomeric complexes have been located, with 13 of them containing cyclic ligands. The global minimum structure (2a) is composed of 1a, with the Li+ cation linked to the four negatively charged sulfur atoms (symmetry C(s)). The corresponding complex 2c containing the ligand 1b is by 23.4 kJ mol(-1) less stable than 2a, and a bicyclic crown-shaped LiS7 cation (2e) is by 34.9 kJ mol(-1) less stable than 2a. Even less stable are four complexes with the branched S6=S ligand. SS bond activation by polarization of the valence electrons takes place on coordination of Li+ to cyclo-S7 (1a).     

Chair, boat and twist conformation of dodecamethylcyclohexasilane and undecamethylcyclohexasilane: a combined DFT and Raman spectroscopic study.

The conformations of dodecamethylcyclohexasilane Si6Me12 and undecamethylcyclohexasilane Si6Me11H have been investigated by ab initio calculations employing the B3LYP density functional with a 6-31+G(d) basis set. Local minima as well as transition structures were calculated with imposed symmetry constraints. For Si6Me12, three unique minima, which correspond to the chair, twist and boat conformations were located with relative zero-point-vibration-corrected energies of 0.0, 7.8 and 11.4 kJ mol(-1). A half-chair conformation with four coplanar silicon atoms connects the chair and twisted minima via an energy barrier of 16.0 and 8.2 kJ mol(-1), respectively. A second transition structure with a barrier of 3.9/0.3 kJ mol(-1) connects the twist with the boat structure. Solution Raman spectra of Si6(CH3)12 and Si6(CD3)12 fully corroborate these results. Below -40 degrees C, the symmetric SiSi ring breathing vibration is a single line, which develops a shoulder (originating from the twist conformer) at longer wavelengths whose intensity increases with increasing temperature. From a Van't Hoff plot, the chair/twist enthalpy difference is 6.6+/-1.5 kJ mol(-1) for Si6(CH3)12 and 6.0+/-1.5 kJ mol(-1) for Si6(CD3)12, which is in reasonable agreement with the ab initio results. Due to the low barrier, the boat conformation cannot be observed, because either the lowest torsional vibration level lies above it or a rapid interconversion between the twist and boat conformations occurs, resulting in averaged Raman spectra. For Si6Me11H, six local minima were located. The chair with the hydrogen atom in the axial position (axial chair) is the global minimum, followed by the equatorial chair (+1.9 kJ mol(-1)) and the three twist conformers (+5.3, +8.0 and +8.1 kJ mol(-1)). The highest local minimum (+11.9 kJ mol(-1)) is a C(s) symmetric boat with the hydrogen atom in the equatorial position. Two possible pathways for the chair-to-chair interconversion with barriers of 13.9 and 14.5 kJ mol(-1) have been investigated. The solution Raman spectra in the SiSi ring breathing region clearly show that below -50 degrees C only the axial and equatorial chairs are present, with an experimental deltaH-value of 0.46 kJ mol(-1). With increasing temperature a shoulder develops which is attributed to the combined twist conformers. The experimental deltaH-value is 6.9 kJ mol(-1), in good agreement with the ab initio results. Due to the low interconversion barriers, the various twist conformers cannot be detected separately.     

The OH* + CH3SH reaction: support for an addition-elimination mechanism from ab initio calculations

Several intermediates for the CH(3)SH + OH(*) --> CH(3)S(*) + H(2)O reaction were identified using MP2(full) 6-311+g(2df,p) ab initio calculations. An adduct, CH(3)S(H)OH(*), I, with electronic energy 13.63 kJ mol(-1) lower than the reactants, and a transition state, II(double dagger), located 5.14 kJ mol(-1) above I, are identified as the entrance channel for an addition-elimination reaction mechanism. After adding zero-point and thermal energies, DeltaH(r,298) ( degrees )(reactants --> I) = -4.85 kJ mol(-1) and DeltaH(298) (double dagger)(I --> II(double dagger)) = +0.10 kJ mol(-1), which indicates that the potential energy surface is broad and flat near the transition state. The calculated imaginary vibrational frequency of the transition state, 62i cm(-1), is also consistent with an addition-elimination mechanism. These calculations are consistent with experimental observations of the OH(*) + CH(3)SH reaction that favored an addition-elimination mechanism rather than direct hydrogen atom abstraction. An alternative reaction, CH(3)SH + OH(*) --> CH(3)SOH + H(*), with DeltaH(r,298) ( degrees ) = +56.94 kJ mol(-1) was also studied, leading to a determination of DeltaH(f,298) ( degrees )(CH(3)SOH) = -149.8 kJ mol(-1).     

Preparation and solid-state characterization of the novel mixed biradical *NSNSC-CNSSN*

Reduction of the radical-cation [*NSSNC-CNSNS][AsF(6)] with ferrocene affords the novel biradical *NSNSC-CNSSN* containing both 1,2,3,5- and 1,3,2,4-isomeric dithiadiazolyl rings. Biradicals form centrosymmetric dimers with pi(*)-pi(*) interactions between different isomeric rings. Biradical *NSNSC-CNSSN* is diamagnetic in the solid state (C = 0.00035, TIP = 6.5 x 10(-)(5) emu/Oe.mol); however, an increase in paramagnetism was observed upon grinding (C = 0.003, TIP = 4.2 x 10(-)(4) emu/Oe.mol).     

Ab initio calculation of boat benzenes in bis(2,4,6-tri-t-butylphenyl) phosphinic chloride

ab initio Calculation of the benzene skeletons of the titled compound indicated that the molecule suffers 61.4 kcal/mol higher energy than that of flat benzene rings and that the boat benzene form is even more stable than the envelope or chair form benzene.     

Structures and vibrational spectra of the sulfur-rich oxides SnO (n = 4-9): the importance of pi*-pi* interactions

The structures of a large number of isomers of the sulfur oxides S(n)O with n = 4-9 have been calculated at the G3X(MP2) level of theory. In most cases, homocyclic molecules with exocyclic oxygen atoms in an axial position are the global minimum structures. Perfect agreement is obtained with experimentally determined structures of S(7)O and S(8)O. The most stable S(4)O isomer as well as some less stable isomers of S(5)O and S(6)O are characterized by a strong pi*-pi* interaction between S==O and S==S groups, which results in relatively long S--S bonds with internuclear distances of 244-262 pm. Heterocyclic isomers are less stable than the global minimum structures, and this energy difference approximately increases with the ring size: 17 (S(4)O), 40 (S(5)O), 32 (S(6)O), 28 (S(7)O), 45 (S(8)O), and 54 kJ mol(-1) (S(9)O). Owing to a favorable pi*-pi* interaction, preference for an axial (or endo) conformation is calculated for the global energy minima of S(7)O, S(8)O, and S(9)O. Vapor-phase decomposition of S(n)O molecules to SO(2) and S(8) is strongly exothermic, whereas the formation of S(2)O and S(8) is exothermic if n<7, but slightly endothermic for S(7)O, S(8)O, and S(9)O. The calculated vibrational spectra of the most stable isomers of S(6)O, S(7)O, and S(8)O are in excellent agreement with the observed data.     

Ab initio molecular orbital study of dications of 1,5‐dichalcogenacyclooctanes

The molecular and electronic structures of the dications of three homonuclear and three heteronuclear dichalcogenacyclooctanes (chalcogen = S, Se, or Te) were investigated by ab initio molecular orbital calculations. Four energy‐minimum structures were located for each dication. Three of those (chair‐chair, boat‐boat, and boat‐chair) have the cis configuration with respect to the chalcogen lone pairs, and the remaining one has the trans configuration. The cis isomers were found to be much more stable than the trans isomer. Among the three cis structures, the stability is in the order of boat‐chair > boat‐boat > chair‐chair for all dications. This order can be explained by considering the nonbonding H···H interactions. The chair‐chair structure (C_2v symmetry) of the 1,5‐dithiacyclooctane dication has a very low vibrational frequency of a₂ symmetry, and its LUMO energy is lower than those of boat‐boat and boat‐chair. These can rationalize the fact that in the crystalline state the dication adopts a distorted C₂ chair‐chair conformation. The transition states between the three conformers of the homonuclear dications were also located. The corresponding energy barriers are relatively low, which is consistent with their NMR spectra. The relative stabilities of the homonuclear and heteronuclear dications were elucidated on the basis of their energies and those of the corresponding neutral compounds. © 2000 John Wiley & Sons, Inc. Heteroatom Chem 11:31–41, 2000     

从头计算法研究沙蚕毒和杀虫环分子的几何构型与电子结构

用Hartree-Fock/6-31G^*从头算确定了沙蚕毒和杀虫环分子的几何构型,在全局优化中发现杀虫环分子的椅式和船式两种稳定构象,在二级Moller-plesset微扰理论MP2/6-31^*水平下,椅式较船式稳定27.06kJ/mol.用MP2/6-31G^*波函数计算电子相关校正的分子静电势,以此为基础讨论生物活性与静电势的关系。发现对此二分子,Mulliken布居分析获得的原子净电荷存在问题,本文用Breneman提出的从静电势导出原子净电荷的CHELPG方法计算了原子净电荷。     

S3O2--an unusual structure with pi*-pi* interaction

The structures and relative stabilities of 15 S3O2 isomers have been investigated by G3X(MP2), CCSD(T)/aug-cc-pVTZ and MRCI/CASSCF calculations. The global energy minimum is a three-membered sulfur ring with two adjacent sulfoxide groups in a trans conformation, i.e. a vic-disulfoxide of C2 symmetry. The SS bond lengths are 2.136 (2x) and 2.354 angstroms at the CCSD(T)/cc-pVTZ level of theory. There is a strong interaction between the pi* orbitals of the two S=O moieties both in the trans and in the almost degenerate cis conformer. The corresponding chain-like singlet and triplet isomers of connectivity OSSSO lie close in energy (ca. 67 kJ mol(-1)) while five-membered and branched four-membered rings are significantly less stable. The structure of S3O2 is in contrast to that of the isoelectronic analogue S5, which exists as a five-membered twisted heterocycle.     

A PM3 Molecular Orbital Study on the Conformational Equilibrium of Cyclohexane upon α-Cyclodextrin Inclusion Complexation

The effect ofα-cyclodextrin (α-CD) inclusion complexation on the conformational equilibrium of cyclo- hexane wasstudied with thesemiempirical PM3 molecularorbital calculations. The calculation results indicated that the chair form of cyclohexane is 18.5 kJ·mol-1lower than that of boat one in energy, however, theα-cyclodextrin inclusion complex of boat cycl ohexane is 4.4 kJ·mol-1more stable than the complex of chair form. It demon-strated that the conformational equilibrium of cyclohexane was influenced by theα-CD inclusion complexation. Hence, caution should be given when extrapolating the conformational behaviors of the guest compounds in the supramolecular systems totheir free forms, since the interactionsbetween the host and guest significantly affect the conformation of the guest compounds.     

CASSCF and CASPT2 calculations on the cleavage and ring inversion of bicyclo[2.2.0]hexane find that these reactions involve formation of a common twist-boat diradical intermediate.

(6/6)CASSCF and CASPT2/6-31G calculations have been performed to understand the experimental finding of Goldstein and Benzon (J. Am. Chem. Soc. 1972, 94, 5119) that exo-bicyclo[2.2.0]hexane-d(4) (1b) undergoes ring inversion to form endo-bicyclo[2.2.0]hexane-d(4) (4b) faster than it undergoes cleavage to form cis,trans-1,5-hexadiene-d(4) (3b). Goldstein and Benzon also found that the latter reaction, which must occur via a chairlike transition structure (TS), is much faster than cleavage of 1b to trans,trans-1,5-hexadiene-d(4) (2b) via a boatlike TS. Our calculations reveal that all three of these reactions involve ring opening of 1, through a boat diradical TS (BDTS), to form a twist-boat diradical intermediate (TBDI). TBDI can reclose to 4 via a stereoisomeric boat diradical TS (BDTS'), or TBDI can cleave, either via a half-chair diradical TS (HCDTS) to form 3 or via a boat TS (BTS) to form 2. The calculated values of DeltaH(++) = 34.6 kcal/mol, DeltaS(++) = -1.6 eu, and DeltaH(++) = 35.2 kcal/mol, DeltaS(++) = 2.0 eu for ring inversion of 1 to 4 and cleavage of 1 to 3, respectively, are in excellent agreement with the values measured by Goldstein and Benzon. The higher value of DeltaH(++) = 37.6 kcal/mol, computed for cleavage of TBDI to 2, is consistent with the experimental finding that very little 2b is formed when 1b is pyrolyzed. The relationships between BDTS, HCDTS, and BTS and the chair and boat Cope rearrangement TSs (CCTS and BCTS) are discussed.     

New insights into an old reaction. High-resolution x-ray powder diffraction of Wiberg's aminoalane intermediate

In accordance with the procedure described by E. Wiberg, Me(3)Al-NH(3) was heated as a bulk material in inert atmosphere to give a colorless liquid which slowly loses methane. Close to the end of this elimination reaction, the melt crystallized to give a microcrystalline powder of (Me(2)AlNH(2))(x)(). The structure of this intermediate has been solved by the method of high-resolution X-ray powder diffraction. The compound crystallizes in the monoclinic space group C2/c with the cell parameters of a = 15.0047(6) A, b = 8.7500(2) A, c = 24.4702(8) A, and beta = 107.290(2) degrees, with eight trimers (Me(2)AlNH(2))(3) per unit cell. These trimers crystallize in a boat conformation in contrast to the known trimers of the same composition where a twist-boat conformation had been found by single crystal determination. Different conformers of (Me(2)AlNH(2))(3) have been investigated by theoretical methods (HF/6-31G(d), B3LYP/6-31G(d), B3LYP/6-311G(d,p), MP2(fc)/6-31G(d), and MP2(fc)/6-311G(d,p)). The twist-boat and the chair conformer correspond to minima at the potential energy surface, whereas the boat conformer corresponds to a first-order transition state (relative energies of 0.45-2.56 kJ/mol (boat) and 6.66-11.91 kJ/mol (chair)). Relaxed scans of the potential energy surface at the HF/6-31G(d) and B3LYP/6-31G(d) levels have shown that the boat conformer (C(s)() symmetry) connects two enantiomers of the twist-boat form (C(2) symmetry).     

Structure and dynamics of binary and ternary lanthanide(III) and actinide(III) tris[4,4,4-trifluoro-1-(2-thienyl)-1,3-butanedione] (TTA)-tributylphosphate (TBP) complexes. Part 3, the structure, thermodynamics and reaction mechanisms of 8- and 9-coordinated binary and ternary Y-TTA-TBP complexes studied by quantum chemical methods

Possible mechanisms for intermolecular exchange between coordinated and solvent water in the complexes Y(TTA)(3)(OH(2))(2) and Y(TTA)(3)(TBP)(OH(2)) and intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(OH(2))(HTTA) and Y(TTA)(3)(TBP)(HTTA) have been investigated using ab initio quantum chemical methods. The calculations comprise both structures and energies of isomers, intermediates and transition states. Based on these data and experimental NMR data (Part 2) we have suggested intimate reaction mechanisms for water exchange, intramolecular exchange between structure isomers and intermolecular exchange between free HTTA and coordinated TTA. A large number of isomers are possible for the complexes investigated, but only some of them have been investigated, in all of them the most stable geometry is a more or less distorted square anti-prism or bicapped trigonal prism; the energy differences between the various isomers are in general small, less than 10 kJ mol(-1). 9-coordinated intermediates play an important role in all reactions. Y(TTA)(3)(OH(2))(3) has three non-equivalent water ligands that can participate in ligand exchange reactions. The fastest of these exchanging sites has a QM activation energy of 18.1 kJ mol(-1), in good agreement with the experimental activation enthalpy of 19.6 kJ mol(-1). The mechanism for the intramolecular exchange between structure isomers in Y(TTA)(3)(OH(2))(2) involves the opening of a TTA-ring as the rate determining step as suggested by the good agreement between the QM activation energy and the experimental activation enthalpy 47.8 and 58.3 J mol(-1), respectively. The mechanism for the intermolecular exchange between free and coordinated HTTA in Y(TTA)(3)(HTTA) and Y(TTA)(3)(TBP)(HTTA) involves the opening of the intramolecular hydrogen bond in coordinated HTTA followed by proton transfer to coordinated TTA. This mechanism is supported by the good agreement between experimental activation enthalpies (within parenthesis) and calculated activation energies 68.7 (71.8) and 35.3 (38.8) kJ mol(-1). The main reason for the difference between the two systems is the much lower energy required to open the intramolecular hydrogen bond in the latter. The accuracy of the QM methods and chemical models used is discussed.     

Synthetic and computational studies of the palladium(iv) system Pd(alkyl)(aryl)(alkynyl)(bidentate)(triflate) exhibiting selectivity in C-C reductive elimination

Synthetic routes to methyl(aryl)alkynylpalladium(iv) motifs are presented, together with studies of selectivity in carbon-carbon coupling by reductive elimination from Pd(IV) centres. The iodonium reagents IPh(C[triple bond, length as m-dash]CR)(OTf) (R = SiMe(3), Bu(t), OTf = O(3)SCF(3)) oxidise Pd(II)Me(p-Tol)(L(2)) (1-3) [L(2) = 1,2-bis(dimethylphosphino)ethane (dmpe) (1), 2,2'-bipyridine (bpy) (2), 1,10-phenanthroline (phen) (3)] in acetone-d(6) or toluene-d(9) at -80 °C to form complexes Pd(IV)(OTf)Me(p-Tol)(C[triple bond, length as m-dash]CR)(L(2)) [R = SiMe(3), L(2) = dmpe (4), bpy (5), phen (6); R = Bu(t), L(2) = dmpe (7), bpy (8), phen (9)] which reductively eliminate predominantly (>90%) p-Tol-C[triple bond, length as m-dash]CR above ～-50 °C. NMR spectra show that isomeric mixtures are present for the Pd(IV) complexes: three for dmpe complexes (4, 7), and two for bpy and phen complexes (5, 6, 8, 9), with reversible reduction in the number of isomers to two occurring between -80 °C and -60 °C observed for the dmpe complex 4 in toluene-d(8). Kinetic data for reductive elimination from Pd(IV)(OTf)Me(p-Tol)(C[triple bond, length as m-dash]CSiMe(3))(dmpe) (4) yield similar activation parameters in acetone-d(6) (66 ± 2 kJ mol(-1), ΔH(?) 64 ± 2 kJ mol(-1), ΔS(?)-67 ± 2 J K(-1) mol(-1)) and toluene-d(8) (E(a) 68 ± 3 kJ mol(-1), ΔH(?) 66 ± 3 kJ mol(-1), ΔS(?)-74 ± 3 J K(-1) mol(-1)). The reaction rate in acetone-d(6) is unaffected by addition of sodium triflate, indicative of reductive elimination without prior dissociation of triflate. DFT computational studies at the B97-D level show that the energy difference between the three isomers of 4 is small (12.6 kJ mol(-1)), and is similar to the energy difference encompassing the six potential transition state structures from these isomers leading to three feasible C-C coupling products (13.0 kJ mol(-1)). The calculations are supportive of reductive elimination occurring directly from two of the three NMR observed isomers of 4, involving lower activation energies to form p-TolC[triple bond, length as m-dash]CSiMe(3) and earlier transition states than for other products, and involving coupling of carbon atoms with higher s character of σ-bonds (sp(2) for p-Tol, sp for C[triple bond, length as m-dash]C-SiMe(3)) to form the product with the strongest C-C bond energy of the potential coupling products. Reductive elimination occurs predominantly from the isomer with Me(3)SiC[triple bond, length as m-dash]C trans to OTf. Crystal structure analyses are presented for Pd(II)Me(p-Tol)(dmpe) (1), Pd(II)Me(p-Tol)(bpy) (2), and the acetonyl complex Pd(II)Me(CH(2)COMe)(bpy) (11).     

Experimental and theoretical study of the photodissociation reaction of thiophenol at 243 nm: intramolecular orbital alignment of the phenylthiyl radical

The photoinduced hydrogen (or deuterium) detachment reaction of thiophenol (C(6)H(5)SH) or thiophenol-d(1) (C(6)H(5)SD) pumped at 243 nm has been investigated using the H (D) ion velocity map imaging technique. Photodissociation products, corresponding to the two distinct and anisotropic rings observed in the H (or D) ion images, are identified as the two lowest electronic states of phenylthiyl radical (C(6)H(5)S). Ab initio calculations show that the singly occupied molecular orbital of the phenylthiyl radical is localized on the sulfur atom and it is oriented either perpendicular or parallel to the molecular plane for the ground (B(1)) and the first excited state (B(2)) species, respectively. The experimental energy separation between these two states is 2600+/-200 cm(-1) in excellent agreement with the authors' theoretical prediction of 2674 cm(-1) at the CASPT2 level. The experimental anisotropy parameter (beta) of -1.0+/-0.05 at the large translational energy of D from the C(6)H(5)SD dissociation indicates that the transition dipole moment associated with this optical transition at 243 nm is perpendicular to the dissociating S-D bond, which in turn suggests an ultrafast D+C(6)H(5)S(B(1)) dissociation channel on a repulsive potential energy surface. The reduced anisotropy parameter of -0.76+/-0.04 observed at the smaller translational energy of D suggests that the D+C(6)H(5)S(B(2)) channel may proceed on adiabatic reaction paths resulting from the coupling of the initially excited state to other low-lying electronic states encountered along the reaction coordinate. Detailed high level ab initio calculations adopting multireference wave functions reveal that the C(6)H(5)S(B(1)) channel may be directly accessed via a (1)(n(pi),sigma(*)) photoexcitation at 243 nm while the key feature of the photodissociation dynamics of the C(6)H(5)S(B(2)) channel is the involvement of the (3)(n(pi),pi(*))-->(3)(n(sigma),sigma(*)) profile as well as the spin-orbit induced avoided crossing between the ground and the (3)(n(pi),sigma(*)) state. The S-D bond dissociation energy of thiophenol-d(1) is accurately estimated to be D(0)=79.6+/-0.3 kcalmol. The S-H bond dissociation energy is also estimated to give D(0)=76.8+/-0.3 kcalmol, which is smaller than previously reported ones by at least 2 kcalmol. The C-H bond of the benzene moiety is found to give rise to the H fragment. Ring opening reactions induced by the pi-pi(*)n(pi)-pi(*) transitions followed by internal conversion may be responsible for the isotropic broad translational energy distribution of fragments.     

β‐Amyrin Biosynthesis: The Methyl‐30 Group of (3S)‐2,3‐Oxidosqualene Is More Critical to Its Correct Folding To Generate the Pentacyclic Scaffold than the Methyl‐24 Group

Oxidosqualene cyclases catalyze the transformation of oxidosqualene ( 1 ) into numerous cyclic triterpenes. Enzymatic reactions of 24‐noroxidosqualene ( 8 ) and 30‐noroxidosqualene ( 9 ) with Euphorbia tirucalli β‐amyrin synthase were conducted to examine the role of the branched methyl groups of compound 1 in the β‐amyrin biosynthesis. Substrate 8 almost exclusively afforded 30‐nor‐β‐amyrin (>95.5 %), which was produced through a normal cyclization pathway, along with minor products (<4.5 %). However, a lack of the Me‐30 group (analogue 9 ) resulted in significantly high production of premature cyclization products, including 6/6/6/5‐fused tetracyclic and 6/6/6/6/5‐fused pentacyclic skeletons (64.6 %). In addition, the fully cyclized product (35.4 %) having the 6/6/6/6/6‐fused pentacycle was produced; however, the normally cyclized product, 29‐nor‐β‐amyrin was present in only 18.6 % of these products. The conversion yield of substrate 8 possessing a Z‐Me group at the terminus was approximately twofold greater than that of compound 9 with an E‐Me group. Thus, the Me‐30 group is essential for the correct folding of a chair–chair–chair–boat–boat conformation of compound 1 for the production of the β‐amyrin scaffold, whereas the Me‐24 group exerts little influence on the normal polycyclization cascade. Here, we show that the Me‐30 group plays critical roles in constructing the ordered architecture of a chair–chair–chair–boat–boat structure, in facilitating the ring‐expansion reactions, and in performing the final deprotonation reaction at the correct position.     

The conversion among various B4C clusters: a density functional theoretical study

Geometry optimizations and vibration frequencies of B4C clusters were performed with Becke-3LYP method using 6-31G(d) basis set. We have found 14 stable isomers, and the most stable structure among them is the five-member ring containing two three-member boron rings. We also analyzed these stable isomers in detail, and the results show that the structures containing three-member boron rings are predominant in energy for B4C clusters. In terms of MO and NBO analysis, the three-centered bond and the pi-electron delocalization play an important role in stabilizing the planar five-member rings of these B4C clusters. Our calculations suggest that isomer4 can be converted into isomer7 with only an energy barrier of 0.31 kJ mol(-1) at the B3LYP/6-311G+(3df) level. Although the planar structures of the five-member rings (isomers12-14) can be converted with each other, the conversions of isomer14 to isomer13 and isomer13 to isomer12 have high-energy barriers of 70.99 and 68.51 kJ mol(-1) at the B3LYP/6-31G(d) level, respectively.     

S2O2分子稳定结构的密度泛函理论研究

采用密度泛函理论中的B3LYP方法,在6-311+G(3df)和Aug-cc-pVTZ水平上,研究了单态S2O2分子的各种可能的异构体及其相对稳定性。结果表明,具有C2v对称性的三角平面分叉异构体的热力学稳定性要高于目前实验上唯一被发现的具有C2v对称性的cis-OSSO异构体,同时trans-OSSO的稳定性与顺式异构体十分接近,这2种异构体应该可以在实验上被观察到。同时本文还讨论了3个最稳定构型的前线分子轨道和链型OSSO的内扭转势能。     

密度泛函方法研究HS2X的稳定构型

采用量子化学中的密度泛函理论(DFT),在B3LYP／6-311 G(3df,2p)水平上全优化得到了单卤代二硫烷HSSX(X-H,F,Cl,Br)链型和分叉型两种异构体的平衡结构。计算结果表明．在热力学上,所有链型的HSSX为稳定构型．其能量分别较分叉型构型低109．8、60．2、74．8和73．1kJ／mol。同时,采用统计热力学方法,研究了两种平衡结构之间相互转化的化学热力学性质,发现异构化反应的平衡常数很小。