Ab initio studies on the reactivity of the CF3OCH2O radical: Thermal decomposition vs. reaction with O2

Hydrofluoroethers are being considered as potential candidates for third generation refrigerants. The present investigation involves the ab initio quantum mechanical study of the decomposition mechanism of CF₃OCH₂O radical formed from a hydrofluoroether, CF₃OCH₃ (HFE-143a) in the atmosphere. The geometries of the reactant, products and transition states involved in the decomposition pathways are optimized and characterized at the DFT (B3LYP) level of theory using 6-311G(d,p) basis set. Energy calculations have been performed at the G2(MP2) and G2M(CC,MP2) level of theory. Two prominent decomposition channels, C-O bond scission and reaction with atmospheric O₂ have been considered for detailed investigation. Studies performed at the G2(MP2) level reveals that the decomposition channel involving C-O bond scission occurs with a barrier height of 23.8 kcal mol⁻¹ whereas the oxidative pathway occurring with O₂ proceeds with an energy barrier of 7.2 kcal mol⁻¹. On the other hand the corresponding values at G2M(CC,MP2) are 24.5 and 5.9 kcal mol⁻¹ respectively. Using canonical transition state theory (CTST) rate constants for the two pathways considered are calculated at 298 K and 1 atm pressure and found to be 5.9 × 10⁻⁶ s⁻¹ and 2.3 × 10⁻⁵ s⁻¹ respectively. The present study concludes that reaction with O₂ is the dominant path for the consumption of CF₃OCH₂O in the atmosphere. Transition states are searched and characterized on the potential energy surfaces involved in both of the reaction channels. The existence of transition state on the corresponding potential energy surface is ascertained by performing intrinsic reaction coordinate (IRC) calculation.     

Hydrogen atom abstraction from HOOOCl by chlorine atom and OH radical

Gas-phase reactions of HOOOCl with both Cl atom and OH radical are investigated using ab initio methods. The structures of all reactants, products, intermediates, and transition states have been optimized and characterized with the quadratic configuration interaction (QCISD) method. The overall mechanism for the Cl + HOOOCl and OH + HOOOCl reaction is the formation of HCl + O₂ + ClO and H₂O + O₂ + ClO, respectively. The rate-limiting step in each reaction is the abstraction of hydrogen from HOOOCl by either Cl or OH radicals and the barrier height is predicted to be 1.9 kcal mol⁻¹ and 8.1 kcal mol⁻¹ for abstraction by Cl atom and OH radical, respectively. Since both barriers for hydrogen abstraction are high, the reaction is suggested to be slow. These results also suggest that an atmospheric removal mechanism for HOOOCl may result from reaction with Cl atoms rather than with OH radicals, and that photolysis of HOOOCl may be the major removal mechanism for the intermediate.     

A Comparative Semiempirical, Ab initio and DFT Study of Decarbonylation of Cyclic Ketones Part 1: Thermal Fragmentation of Cyclopropanone

Calculations using different quantum mechanical methods including semiempirical (MNDO,AM1 and PM3), ab initio (RHF and MP2 calculations using the 6-311G and 6-311++G^** basis sets), and density functional theory (LSDA, BP, MIXBP and B3LYP, i.e., B3LYP/6-311+G**//B3LYP/6-31G*) have been performed on the thermal fragmentation of cyclopropanone to ethylene and carbon monoxide. All RHF calculations predict a concerted single step mechanism for this conversion. The estimated activation energies vary from 34.4 to 54.6 kcal·mol^-1, mainly localized around 37±2 kcal·mol^-1, depending on the method. Whereas the calculated RHF reaction energies also varied from 14.5 to -33.3 kcal·mol^-1, the B3LYP/6-311+G**//B3LYP/6-31G* method predicts the experimental value (-17.7 kcal·mol^-1) within experimental uncertainties. Remarkably, semiempirical AM1 and PM3 methods and simple DFT calculations, LSDA, predict comparable results to the more advanced methods. UHF ab initio calculations predict the same single step mechanism, whereas a multistep biradical mechanism with an unrealistically low activation energy is favored by the semiempirical methods. Structures of the activated complex of the single step mechanism, estimated by different methods, are very similar and consistent with a nonlinear cheletropic [_2s + _2a] reaction, as predicted by the orbital symmetry rules and earlier EHT calculations.Electronic Supplementary Material available.     

Ab initio studies on the decomposition kinetics of CF<Subscript>3</Subscript>OCF<Subscript>2</Subscript>O radical

The present study deals with the decomposition of CF₃OCF₂O radical formed from a hydrofluoroether, CF₃OCHF₂ (HFE-125), in the atmosphere. The study is performed using ab initio quantum mechanical methods. Two plausible pathways of decomposition of the titled species have been considered, one involving C-O bond scission and the other occurring via F atom elimination. The geometries of the reactant, products and transition states involved in the decomposition pathways are optimized and characterized at DFT (B3LYP) level of theory using 6-311G(d,p) basis set. Single point energy calculations have been performed at G2M(CC,MP2) level of theory. Out of the two prominent decomposition channels considered, the C-O bond scission is found to be dominant involving a barrier height of 15.3 kcal mol⁻¹ whereas the F-elimination path proceeds with a barrier of 26.1 kcal mol⁻¹. The thermal rate constants for the above two decomposition pathways are evaluated using canonical transition state theory (CTST) and these are found to be 1.78 × 10⁶ s⁻¹ and 2.83 × 10⁻⁷ s⁻¹ for C-O bond scission and F-elimination respectively at 298 K and 1 atm pressure. Transition states are searched on the potential energy surfaces involved during the decomposition channels and each of the transition states is characterized. The existence of transition states on the corresponding potential energy surface is ascertained by performing intrinsic reaction coordinate (IRC) calculation.     

Quantum chemical studies on substitution effects within silyl group in the silylative coupling of olefins

In this study quantum chemical calculations based on the density functional theory (DFT) have been carried out to examine the effects of methoxy substituent attached to a silicon atom on the reaction of silylative coupling of olefins. It has been shown, that substituted substrate undergoes the reaction according to the recently proposed insertion-rotation-elimination mechanism. During the rotation around C-C single bond additional stabilization by oxygen-ruthenium interaction was observed. Similarly to the (trimethylsilyl)ethene the rate determining step of the reaction is the insertion of the alkene into Ru-Si single bond. The substitution of SiMe₃ by Si(OMe)₃ decreases the energy span of the reaction by almost 3 kcal mol^-1 that is from 21 kcal mol^-1 to 18 kcal mol^-1. The decrease of the energy barrier of the reaction seems to be the result of the increase of point charge differences between the Ru and Si atoms which increases electrostatic attraction between these atoms. Moreover, for Si(OMe)₃ the rate-determining transition state is closer to the alkene interacting with the Ru centre side of the reaction.     

Cα torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

A multipolar, polarizable electrostatic method for future use in a novel force field is described. Quantum Chemical Topology (QCT) is used to partition the electron density of a chemical system into atoms, then the machine learning method Kriging is used to build models that relate the multipole moments of the atoms to the positions of their surrounding nuclei. The pilot system serine is used to study both the influence of the level of theory and the set of data generator methods used. The latter consists of: (i) sampling of protein structures deposited in the Protein Data Bank (PDB), or (ii) normal mode distortion along either (a) Cartesian coordinates, or (b) redundant internal coordinates. Wavefunctions for the sampled geometries were obtained at the HF/6-31G(d,p), B3LYP/apc-1, and MP2/cc-pVDZ levels of theory, prior to calculation of the atomic multipole moments by volume integration. The average absolute error (over an independent test set of conformations) in the total atom-atom electrostatic interaction energy of serine, using Kriging models built with the three data generator methods is 11.3 kJ mol^-1 (PDB), 8.2 kJ mol^-1 (Cartesian distortion), and 10.1 kJ mol^-1 (redundant internal distortion) at the HF/6-31G(d,p) level. At the B3LYP/apc-1 level, the respective errors are 7.7 kJ mol^-1, 6.7 kJ mol^-1, and 4.9 kJmol^-1, while at the MP2/cc-pVDZ level they are 6.5 kJ mol^-1, 5.3 kJ mol^-1, and 4.0 kJmol^-1. The ranges of geometries generated by the redundant internal coordinate distortion and by extraction from the PDB are much wider than the range generated by Cartesian distortion. The atomic multipole moment and electrostatic interaction energy predictions for the B3LYP/apc-1 and MP2/cc-pVDZ levels are similar, and both are better than the corresponding predictions at the HF/6-31G(d,p) level.     

Reductive elimination pathway for homocysteine to methionine conversion in cobalamin-dependent methionine synthase

Density functional theory has been applied to investigate the methyl transfer from methylcobalamin (MeCbl) cofactor to homocysteine (Hcy) as catalyzed by methionine synthase (MetH). Specifically, the S_N2 and the reductive elimination pathways have been probed as the possible mechanistic pathways for the methyl transfer reaction. The calculations indicate that the activation barrier for the reductive elimination reaction (24.4 kcal mol⁻¹) is almost four times higher than that for the S_N2 reaction (7.3 kcal mol⁻¹). This high energy demand of the reductive elimination pathway is rooted in the structural distortion of the corrin ring that is induced en route to the formation of the triangular transition state. Furthermore, the reductive elimination reaction demands the syn accommodation of the methyl group and the substrate over the upper face of the corrin ring, which also accounts for the high energy demand of the reaction. Consequently, the reductive elimination pathway for MetH-catalyzed methyl transfer from MeCbl to Hcy cannot be considered as one of the possible mechanistic routes.     

Direct dynamics simulations of the hydrogen abstraction reaction Cl + CF3CF2CH2OH

The mechanism and kinetics of 2,2,3,3,3-pentafluoropropanol (CF₃CF₂CH₂OH) reaction with Chlorine atom (Cl) is investigated in this work. Two hydrogen abstraction channels of the title reaction are identified. The geometries of all the stationary points in the potential energy surface are obtained at the BHandHLYP/6-311G** level, and the energies of the selected points along the minimum energy path (MEP) are improved by the CCSD(T) method. A dual-level direct dynamics method is employed to study the kinetic nature of the hydrogen-abstraction reaction channels. The calculated rate coefficients show that the hydrogen abstraction from the CH₂ group is the primary channel. The calculated total rate coefficients are in best agreement with the experimental values. The four-parameter rate coefficients expression of the title reaction between the temperatures 200 K and 1000 K is provided.     

Molecular prediction of oseltamivir efficiency against probable influenza A (H1N1-2009) mutants: molecular modeling approach

To predict the susceptibility of the probable 2009 influenza A (H1N1-2009) mutant strains to oseltamivir, MD/LIE approach was applied to oseltamivir complexed with the most frequent drug-resistant strains of neuraminidase subtypes N1 and N2: two mutations on the framework residues (N294S and H274Y) and the two others on the direct-binding residues (E119V and R292K) of oseltamivir. Relative to those of the wild type (WT), loss of drug–target interaction energies, especially in terms of electrostatic contributions and hydrogen bonds were dominantly established in the E119V and R292K mutated systems. The inhibitory potencies of oseltamivir towards the WT and mutants were predicted according to the ordering of binding-free energies: WT (−12.3 kcal mol⁻¹) > N294S (−10.4 kcal mol⁻¹) > H274Y (−9.8 kcal mol⁻¹) > E119 V (−9.3 kcal mol⁻¹) > R292K (−7.7 kcal mol⁻¹), suggesting that the H1N1-2009 influenza with R292K substitution, perhaps, conferred a high level of oseltamivir resistance, while the other mutants revealed moderate resistance levels. This result calls for an urgent need to develop new potent anti-influenza agents against the next pandemic of potentially higher oseltamivir-resistant H1N1-2009 influenza.     

Theoretical investigation of the ring opening process of verdoheme to biliverdin in the presence of dioxygen

The conversion of ferrous verdoheme to ferric biliverdin in the presence of O₂ was investigated using the B3LYP method. Both 6-31G and 6-31G (d) basis sets were employed for geometry optimization calculation as well as energy stabilization estimation. Three possible pathways for the conversion of iron verdoheme to iron biliverdin were considered. In the first route oxygen and reducing electron were employed. In this path formation of ferrous verdoheme-O₂ complex was followed by the addition of one electron to the ferrous-oxycomplex to produce ferric peroxide intermediate. The ferric peroxide intermediate experienced an intramolecular nucleophilic attack to the most positive position at 5-oxo carbons on the ring to form a closed ring biliverdin. Subsequently the ring opening process took place and the iron (III) biliverdin complex was formed. Closed ring iron biliverdin intermediate and open ring iron biliverdin formed as a product of verdoheme cleavage were respectively 13.20 and 32.70 kcal mol⁻¹ more stable than ferric peroxide intermediate. Barrier energy for conversion of ferric peroxide to closed ring Fe (III) biliverdin and from the latter to Fe (III) biliverdin were respectively 8.67 and 3.35 kcal mol⁻¹. In this path spin ground states are doublet except for iron (III) biliverdin in which spin state is quartet. In the second path a ferrous-O₂ complex was formed and, without going to a one electron reduction process, nucleophilic attack of iron superoxide complex took place followed by the formation of iron (III) biliverdin. This path is thermodynamically and kinetically less favorable than the first one. In addition, iron hydro peroxy complex or direct attack of O₂ to macrocycle to form an isoporphyrin type intermediate have shown energy surfaces less favorable than aforementioned routes.     

Instability of 2,2-di(pyridin-2-yl)acetic acid. Tautomerization versus decarboxylation

The DFT calculations at the B3LYP level with 6-311G** basis set were carried out in order to reveal whether tautomerization or decarboxylation is responsible for the instability of 2,2-di(pyridin-2-yl)acetic (DPA) and 1,8-diazafluorene-9-carboxylic (DAF) acids. The carboxyl protons in both compounds are involved in the intramolecular hydrogen bonds (the pyridine nitrogen atoms are the hydrogen bond acceptors). Although formation of two intramolecular OH···N hydrogen bonds in the enols of both carboxylic acids enables effective electron delocalization within the quasi rings (···HO − C = C − C = N), only ene-1,1-diol of DAF has somewhat lower energy than DAF itself (ΔE is ca. 7 kcal mol^-1). DPA and its enediol have comparable energies. Migration of the methine proton toward the carbonyl oxygen atom (to form enediols) requires overstepping the energy barriers of 55-57 kcal mol^-1 for both DPA and DAF. The enaminone tautomers of the acids, formed by migration of this proton toward the pyridine nitrogen atom, are thermodynamically somewhat more stable than the respective enediols. The energy barriers of these processes are equal to ca. 44 and 62 kcal mol^-1 for DPA and DAF, respectively. Thus, such tautomerization of the acids is not likely to proceed. On the other hand, the distinct energetic effects (ca. 15 kcal mol^-1) favor decarboxylation. This process involves formation of (E)-2-(pyridin-2(1H)-ylidenemethyl)pyridine and its cyclic analogue followed by their tautomerization to (dipyridin-2-yl)methane and 1,8-diazafluorene, respectively. Although the later compound was found to be somewhat thermodynamically more stable, kinetic control of tautomerization of the former is more distinct.     

Ab initio and DFT study on the electrophilic addition reaction of bromine to tetracyclo[5.3.0.02,6.03,10]deca–4,8–diene

The electronic and geometric structures of tetracyclo[5.3.0.0^2,6.0^3,10]deca-4,8-diene (hypostrophene) have been investigated by ab initio and DFT/B3LYP methods using the 6-31G* and 6-311G* basis sets. The double bonds of hypostrophene are endo-pyramidalized. The cationic intermediates and products formed in the addition reaction have been investigated using the HF/6-311G*, HF/6-311G**, and B3LYP/6-311G* methods. The bridged bromonium cation was more stable than the U-type cation. Considering that the bridged cation does not isomerize to the less stable U-type cation, it is not possible for the U-type product to be obtained in the reaction. The bridged bromonium cation transformed into the more stable N-type cation and the N-type product was obtained via this cation. The thermodynamic stability of the exo, exo and exo, endo isomers of the N-type dibromide molecule were almost identical. The N-type product was 16.6 kcal mol⁻¹ more stable than the U-type product. Figure General energy diagram of the hypostrophene–bromine (HS–Br₂) system (kcal mol⁻¹) (MP2/6-311G*//HF/6-311G*)     

<Emphasis Type="Italic">Ab initio</Emphasis> and DFT investigation of electrophilic addition reaction of bromine to <Emphasis Type="Italic">endo,endo</Emphasis>-tetracyclo[4.2.1.1<Superscript>3,6</Superscript>.0<Superscript>2,7</Superscript>]dodeca-4,9-diene

Full geometric optimization of endo,endo-tetracyclo[4.2.1.1^3,6.0^2,7]dodeca-4,9-diene (TTDD) has been carried out by ab initio and DFT/B3LYP methods and the structure of the molecule investigated. The double bonds of TTDD molecule are endo pyramidalized. The structure of π-orbitals and their mutual interactions for TTDD molecule were investigated. The cationic intermediates and products obtained as a result of the addition reaction have been studied using the HF/6-311G(d), HF/6-311G(d,p) and B3LYP/6-311G(d) methods. The bridged bromonium cation isomerized into the more stable N- and U-type cations and the difference between the stability of these cations is small. The N- and U-type reaction products are obtained as a result of the reaction, which takes place via the cations in question. The stability of exo, exo and exo, endo isomers of N-type product are nearly the same and the formation of both isomers is feasible. The U-type product basically formed from the exo, exo-isomer. Although the U-type cation was 0.68 kcal mol⁻¹ more stable than the N-type cation, the U-type product was 4.79 kcal mol⁻¹ less stable than the N-type product. Figure The energy diagram of TTDD–Br₂ system (kcal mol⁻¹)(MP2/6-311G*//HF/6-311G*)     

A B3LYP and MP2 theoretical investigation into host-guest interaction between calix[4]arene and Li+ or Na+

The DFT-B3LYP and MP2 methods with 6-311G** and 6-311++G** basis sets have been applied to study the complexation energies of the host-guest complexes between the cone calix[4]arene and Li⁺ or Na⁺ on the B3LYP optimized geometries. A comparison of the complexation energies obtained from the MP2(full) with those from MP2(fc) method is also carried out. The result shows that it is essential to introduce the diffuse basis set into the geometry optimizations and complexation energy calculations of the alkali-metal cation-π interaction complexes of calix[4]arene, and the D _e values show a maximum of 21.13 kJ mol⁻¹ (14.45% of relative error) between the MP2(full)/6-311++G** and MP2(fc)/6-311++G** method. For Li⁺ cation, the complexation is mainly energetically stabilized by the lower rim/cation (namely O–Li⁺) interaction. However, binding energies and NBO analyses confirm that Na⁺ cation prefers to enter the calix[4]arene cavity and the cation-π interaction is predominant, which contradicts the previous low-level theoretical studies. Furthermore, the complexation with Li⁺ is preferred over that with Na⁺ by at least 12.70 kJ mol⁻¹ at MP2(full)/6-311++G**//B3LYP/6-311++G** level.      

Theoretical <Emphasis Type="Italic">ab initio </Emphasis>Study of TiCl<Subscript>4</Subscript> Ammonolysis: Gas Phase Reactions of TiN Chemical Vapor Deposition

The mechanism of TiCl₄ ammonolysis has been studied theoretically at ab initio Hartree-Fock, B3LYP, MP2 and CCSD(T)//B3LYP levels using effective core potentials for Ti and Cl and 6-31G* basis sets for N and H. TiCl₄ and products of its ammonolysis form five- and six-coordinated complexes with ammonia, which intermediate substitution of Cl atoms by NH₂ groups. Transition state energies for the subsequent steps of ammonolysis decrease with increasing number of NH₂ groups bound to Ti. The energy of the transition state for the first step of ammonolysis is 19 kcal mol^-1 above the energy of the reactants (TiCl₄+NH₃) and 8 kcal mol^-1 above the products (TiCl₃NH₂ + HCl). The following steps have transition states energetically located below the products, indicating weak hydrogen bonded complex formation as intermediate between transition state and product. A thermodynamic estimation shows the last step of ammonolysis to be endothermic, while the first three steps are exothermic if the adduct formation energy is taken into account.     

EVALUATION OF SEVERAL ECONOMICAL COMPUTATIONAL METHODS FOR GEOMETRY OPTIMISATION OF PHOSPHORUS ACID DERIVATIVES

Several economical methods for geometry optimisation, applicable to larger molecules, have been evaluated for phosphorus acid derivatives. MP2/cc-pVDZ and B3LYP/6-31+G(d) geometry optimisations are used as reference points, results from geometry optimisations for other methods and their subsequent single point energy calculations are compared to these references. The geometries from HF/MIDI! optimisations were close to those of the references and subsequent single point energies with B3LYP/6-31+G(d,p) or EDF1/6-31+G(d) gave a mean average deviation (MAD) of less than 0.5 kcal mol^?1 from those obtained with the reference geometries.     

Pragmatic ab initio prediction of enthalpies of formation for large molecules: accuracy of MP2 geometries and frequencies using CCSD(T) correlation energies

We have addressed the accuracy of calculating the enthalpy of formation of an arbitrary single reference molecule using practical ab initio methodologies. It is known that MP2 geometries with a triple zeta basis set are almost as reliable as CCSD(T) geometries. It is also known that CCSD(T) correlation energies, with basis extrapolation, feature chemical accuracy for single-reference molecules. We investigate what accuracy one might expect in enthalpies of formation from a MP2 geometry, MP2 harmonic vibrational frequencies, a CCSD(T) correlation energy using triple zeta basis sets. It is far from obvious, a priori, as to which error source contributes most significantly. We observe that the accuracy in calculating enthalpies of formation of single-reference molecules with this protocol is 4 kcal mol^-1; our error analysis shows this comes almost exclusively from the correlation energy basis extrapolation, rather than errors intrinsic to MP2.     

Insight into the interaction of benzothiazole tethered triazole analogues with human serum albumin: Spectroscopy and molecular docking approaches

The interaction of four benzothiazole tethered triazole analogues (MS43, MS70, MS71, and MS78) with human serum albumin (HSA) was investigated using various spectroscopic techniques (ultraviolet–visible (UV–vis) light absorption, fluorescence, circular dichroism (CD), molecular docking and density functional theory (DFT) studies). Fluorescence quenching constants (~10¹²) revealed a static mode of quenching and binding constants (K_b ~10⁴) indicating the strong affinity of these analogues for HSA. Further alteration in the secondary structure of HSA in the presence of these analogues was also confirmed by far UV–CD spectroscopy. The intensity loss in CD studied at 222 nm indicated an increase in random coil/β‐sheet conformations in the protein. Binding energy values (MS71 (?9.3 kcal mol^?1), MS78 (?8.02 kcal mol^?1), MS70 (?7.16 kcal mol^?1) and MS43 (?6.81 kcal mol^?1)) obtained from molecular docking revealed binding of these analogues with HSA. Molecular docking and DFT studies validated the experimental results, as these four analogues bind with HSA at site II through hydrogen bonding and hydrophobic interactions.     

Theoretical investigation on the kinetics and branching ratio of the gas phase reaction of sevoflurane with Cl atom

The present work deals with the theoretical investigation on the Cl initiated H-atom abstraction reaction of sevoflurane, (CF₃)₂CHOCH₂F. A dual-level procedure has been adopted for studying the kinetics of the reaction. Geometrical optimization and frequency calculation were performed at DFT(BHandHLYP)/6-311G(d,p) while single-point energy calculation was made at CCSD(T)/6-311G(d,p) level of theory. The intrinsic reaction coordinate (IRC) calculation has also been performed to confirm the smooth transition from the reactant to product through the respective transition state. The rate constants were calculated using conventional transition state theory (TST). It has been found that 99 % of the reaction proceeded via the H-atom abstraction from the –CH₂F end of the sevoflurane. The rate constant of the dominant path is found to be 1.13 × 10^?13 cm³ molecule^?1 s^?1. This is in excellent agreement with the reported experimental rate constant of 1.10 × 10^?13 cm³ molecule^?1 s^?1 obtained by relative rate method using FTIR/Smog chamber and LP/LIF techniques.     

Theoretical study of N<Subscript>4</Subscript>X (X = O,S, Se) systems

A series of N₄X (X = O, S, Se) compounds have been examined with ab initio and density functional theory (DFT) methods. To our knowledge, these compounds, except for the C_2v ring and the C_3v towerlike isomers of N₄O, are first reported here. The ring structures are the most energetically favored for N₄X (X = O and S) systems. For N₄Se, the cagelike structure is the most energetically favored. Several decomposition and isomerization pathways for the N₄X species have been investigated. The dissociation of C_2v ring N₄O and N₄S structures via ring breaking and the barrier height are only 1.1 and −0.2 kcal mol⁻¹ at the CCSD(T)/6-311+G*//MP2/6-311+G* level of theory. The dissociation of the cagelike N₄X species is at a cost of 12.1–16.2 kcal mol⁻¹. As for the towerlike and triangle bipyramidal isomers, their decomposition or isomerization barrier heights are all lower than 10.0 kcal mol⁻¹. Although the C_S cagelike N₄S isomer has a moderate isomerization barrier (18.3–29.1 kcal mol⁻¹), the low dissociation barrier (−1.0 kcal mol⁻¹) indicates that it will disappear when going to the higher CCSD(T) level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.