Hydrodesulfurization of Dibenzothiophene: A Machine Learning Approach.

The hydrodesulfurization (HDS) process is widely used in the industry to eliminate sulfur compounds from fuels. However, removing dibenzothiophene (DBT) and its derivatives is a challenge. Here, the key aspects that affect the efficiency of catalysts in the HDS of DBT were investigated using machine learning (ML) algorithms. The conversion of DBT and selectivity was estimated by applying Lasso, Ridge, and Random Forest regression techniques. For the estimation of conversion of DBT, Random Forest and Lasso offer adequate predictions. At the same time, regularized regressions have similar outcomes, which are suitable for selectivity estimations. According to the regression coefficient, the structural parameters are essential predictors for selectivity, highlighting the pore size, and slab length. These properties can connect with aspects like the availability of active sites. The insights gained through ML techniques about the HDS catalysts agree with the interpretations of previous experimental reports.

Seeking potential anticonvulsant agents that target GABAA receptors using experimental and theoretical procedures.

The aim of this study was to identify compounds that possess anticonvulsant activity by using a pentylenetetrazol (PTZ)-induced seizure model. Theoretical studies of a set of ligands, explored the binding affinities of the ligands for the GABA(A) receptor (GABA(A)R), including some benzodiazepines. The ligands satisfy the Lipinski rules and contain a pharmacophore core that has been previously reported to be a GABA(A)R activator. To select the ligands with the best physicochemical properties, all of the compounds were analyzed by quantum mechanics and the energies of the highest occupied molecular orbital and lowest unoccupied molecular orbital were determined. Docking calculations between the ligands and the GABA(A)R were used to identify the complexes with the highest Gibbs binding energies. The identified compound D1 (dibenzo(b,f)(1,4)diazocine-6,11(5H,12H)-dione) was synthesized, experimentally tested, and the GABA(A)R-D1 complex was submitted to 12-ns-long molecular dynamics (MD) simulations to corroborate the binding conformation obtained by docking techniques. MD simulations were also used to analyze the decomposition of the Gibbs binding energy of the residues involved in the stabilization of the complex. To validate our theoretical results, molecular docking and MD simulations were also performed for three reference compounds that are currently in commercial use: clonazepam (CLZ), zolpidem and eszopiclone. The theoretical results show that the GABA(A)R-D1, and GABA(A)R-CLZ complexes bind to the benzodiazepine binding site, share a similar map of binding residues, and have similar Gibbs binding energies and entropic components. Experimental studies using a PTZ-induced seizure model showed that D1 possesses similar activity to CLZ, which corroborates the predicted binding free energy identified by theoretical calculations.

Fundamental thermochemical properties of amino acids: gas-phase and aqueous acidities and gas-phase heats of formation.

The gas-phase acidities of the 20 L-amino acids have been predicted at the composite G3(MP2) level. A broad range of structures of the neutral and anion were studied to determine the lowest energy conformer. Excellent agreement is found with the available experimental gas-phase deprotonation enthalpies, and the calculated values are within experimental error. We predict that tyrosine is deprotonated at the CO(2)H site. Cysteine is predicted to be deprotonated at the SH but the proton on the CO(2)H is shared with the S(-) site. Self-consistent reaction field (SCRF) calculations with the COSMO parametrization were used to predict the pK(a)'s of the non-zwitterion form in aqueous solution. The differences in the non-zwitterion pK(a) values were used to estimate the free energy difference between the zwitterion and nonzwitterion forms in solution. The heats of formation of the neutral compounds were calculated from atomization energies and isodesmic reactions to provide the first reliable set of these values in the gas phase. Further calculations were performed on five rare amino acids to predict their heats of formation, acidities, and pK(a) values.

Thermodynamic properties of the XO(2), X(2)O, XYO, X(2)O(2), and XYO(2) (X, Y = Cl, Br, and I) isomers.

High level ab initio electronic structure calculations at the coupled cluster level with a correction for triples extrapolated to the complete basis set limit have been made for the thermodynamics of the BrBrO(2), IIO(2), ClBrO(2), ClIO(2), and BrIO(2) isomers, as well as various molecules involved in the bond dissociation processes. Of the BrBrO(2) isomers, BrOOBr is predicted to be the most stable by 8.5 and 9.3 kcal/mol compared to BrBrO(2) and BrOBrO at 298 K, respectively. The weakest bond in BrOOBr is the O-Br bond with a bond dissociation energy (BDE) of 15.9 kcal/mol, and in BrBrO(2), it is the Br-Br bond of 19.1 kcal/mol. The smallest BDE in BrOBrO is for the central O-Br bond with a BDE of 12.6 kcal/mol. Of the IIO(2) isomers, IIO(2) is predicted to be the most stable by 3.3, 9.4, and 28.9 kcal/mol compared to IOIO, IOOI, and OIIO at 298 K, respectively. The weakest bond in IIO(2) is the I-I bond with a BDE of 22.2 kcal/mol. The smallest BDEs in IOIO and IOOI are the terminal O-I bonds with values of 19.0 and 5.2 kcal/mol, respectively.

Dehydrogenation reactions of cyclic C(2)B(2)N(2)H(12) and C(4)BNH(12) isomers.

The energetics for different dehydrogenation pathways of C(2)B(2)N(2)H(12) and C(4)BNH(12) cycles were calculated at the B3LYP/DGDZVP2 and G3(MP2) levels with additional calculations at the CCSD(T)/complete basis set level. The heats of formation of the different isomers were calculated from the G3(MP2) relative energies and the heats of formation of the most stable isomers of c-C(2)B(2)N(2)H(6), c-C(2)B(2)N(2)H(12), and c-C(4)BNH(12) at the CCSD(T)/CBS including additional corrections together with the previously reported value for c-C(4)BNH(6). Different isomers were analyzed for c-C(2)B(2)N(2)H(x) and c-C(4)BNH(x) (x = 6 and 12), and the most stable cyclic structures were those with C-C-B-N-B-N and C-C-C-C-B-N sequences, respectively. The energetics for the stepwise loss of three H(2) were predicted, and the most feasible thermodynamic pathways were found. Dehydrogenation of the lowest energy c-C(2)B(2)N(2)H(12) isomer (6-H(12)) is almost thermoneutral with DeltaH(3dehydro) = 3.4 kcal/mol at the CCSD(T)/CBS level and -0.6 kcal/mol at the G3(MP2) level at 298 K. Dehydrogenation of the lowest energy c-C(4)BNH(12) isomer (7-H(12)) is endothermic with DeltaH(3dehydro) = 27.9 kcal/mol at the CCSD(T)/CBS level and 23.5 kcal/mol at the G3(MP2) level at 298 K. Dehydrogenation across the B-N bond is more favorable as opposed to dehydrogenation across the B-C, N-C, and C-C bonds. Resonance stabilization energies in relation to that of benzene are reported as are NICS NMR chemical shifts for correlating with the potential aromatic character of the rings.

The effect of the NH2 substituent on NH3: hydrazine as an alternative for ammonia in hydrogen release in the presence of boranes and alanes.

Potential energy surfaces for H(2) release from hydrazine interacting with borane, alane, diborane, dialane and borane-alane were constructed from MP2/aVTZ geometries and zero point energies with single point energies at the CCSD(T)/aug-cc-pVTZ level. With one borane or alane molecule, the energy barrier for H(2)-loss of approximately 38 or 30 kcal mol(-1) does not compete with the B-N or Al-N bond cleavage ( approximately 30 or approximately 28 kcal mol(-1)). The second borane or alane molecule can play the role of a bifunctional catalyst. The barrier energy for H(2)-elimination is reduced from 38 to 23 kcal mol(-1), or 30 to 20 kcal mol(-1) in the presence of diborane or dialane, respectively. The mixed borane-alane dimer reduces the barrier energy for H(2) release from hydrazine to approximately 17 kcal mol(-1). A systematic comparison with the reaction pathways from ammonia borane shows that hydrazine could be an alternative for ammonia in producing borane amine derivatives. The results show a significant effect of the NH(2) substituent on the relevant thermodynamics. The B-N dative bond energy of 31 kcal mol(-1) in NH(2)NH(2)BH(3) is approximately 5 kcal mol(-1) larger than that of the parent BH(3)NH(3). The higher thermodynamic stability could allow hydrazine-borane to be used as a material for certain energetic H(2) storage applications.

Thermochemistry for the dehydrogenation of methyl-substituted ammonia borane compounds.

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for (CH3)H2N-BH3, (CH3)HN=BH2, (BH3)HN=CH2, (CH3)H2B-NH3, (CH3)HB=NH2, and (NH3)HB=CH2, as well as various molecules involved in the different bond-breaking processes, from coupled cluster theory (CCSD(T)) calculations. In order to achieve near-chemical accuracy (+/-1 kcal/mol), three corrections were added to the complete basis set binding energies based on frozen core CCSD(T) energies, corrections for core-valence, scalar relativistic, and first-order atomic spin-orbit effects. Scaled vibrational zero-point energies were computed with the MP2 method. The heats of formation were predicted for the respective dimethyl- and trimethyl-substituted ammonia boranes, their dehydrogenated derivatives, and the various molecules involved in the different bond breaking processes, based on isodesmic reaction schemes calculated at the G3(MP2) level. Thermodynamics for dehydrogenation pathways in the monomethyl-substituted molecules were predicted. Dehydrogenation across the B-N bond is more favorable as opposed to dehydrogenation across the B-C and N-C bonds. Methylation at N reduces the exothermocity of the dehydrogenation reaction and makes the reaction more thermoneutral, while methylation at B moves it away from thermoneutral. Various mixtures of CH3NH2BH3 and NH3BH3 were made, and their melting points were measured. The lowest melting mixture contained approximately 35% NH3BH3 by weight and melted at 35-37 degrees C.

Thermochemistry and electronic structure of small boron and boron oxide clusters and their anions.

Thermochemical properties of a set of small boron (B(n)) and boron oxide (B(n)O(m)) clusters, with n = 1-4 and m = 0-3, their anions, and the B(4)(2-) dianion, were calculated by using coupled-cluster theory CCSD(T) calculations with the aug-cc-pVnZ (n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. Enthalpies of formation, bond dissociation energies, singlet-triplet or doublet-quartet separation gaps, adiabatic electron affinities (EA), and both vertical electron attachment and detachment energies were evaluated. The predicted heats of formation show agreement close to the error bars of the literature results for boron oxides with the largest error for OBO. Our calculated adiabatic EAs are in good agreement with recent experiments: B (calc, 0.26 eV; exptl, 0.28 eV), B(2) (1.95, 1.80), B(3) (2.88, 2.820 +/- 0.020), B(4) (1.68, 1.60 +/- 0.10), BO (2.50, 2.51), BO(2) (4.48, 4.51), BOB (0.07), B(2)O(2) (0.37), B(3)O (2.05), B(3)O(2) (2.94, 2.94), B(4)O (2.58), and B(4)O(2) (3.14, 3.160 +/- 0.015). The BO bond is strong, so this moiety is maintained in most of the clusters. Thermochemical parameters of clusters are not linearly additive with respect to the number of B atoms. The EA tends to be larger in the dioxides. The growth mechanism of small boron oxides should be determined by a number of factors: (i) formation of BO bonds, (ii) when possible, formation of a cyclic B(3) or B(4), and (iii) combination of a boron cycle and a BO bond. When these factors compete, the strength of the BO bonds tends to compensate the destabilization arising from a loss of binding in the cyclic boron clusters, in such a way that a linear boron oxide prevails. When the B(2) moiety is present in these linear clusters, the oxide derivatives prefer a high spin state.

Coordination of aminoborane, NH(2)BH(2), dictates selectivity and extent of H(2) release in metal-catalysed ammonia borane dehydrogenation.

In situ(11)B NMR monitoring, computational modeling, and external trapping studies show that selectivity and extent of H(2) release in metal-catalysed dehydrogenation of ammonia borane, NH(3)BH(3), are determined by coordination of reactive aminoborane, NH(2)BH(2), to the metal center.

ClClO2 is the most stable isomer of Cl2O2. Accurate coupled cluster energetics and electronic spectra of Cl2O2 isomers.

High level ab initio electronic structure calculations at the coupled cluster level with a correction for triples (CCSD(T)) extrapolated to the complete basis set limit have been made for the thermodynamics of the Cl2O2 isomers: ClClO2, ClOOCl, and ClOClO. The ClClO2 isomer is predicted to be the most stable isomer and is more stable than ClOOCl by 3.1 kcal/mol at 298 K. The ClOClO isomer is less stable than ClOOCl by 8.3 kcal/mol at 298 K. The weakest bond in ClClO2 is the Cl-Cl bond with a bond dissociation energy (BDE) of 24.4 kcal/mol, and the smallest BDE in ClOOCl is the O-O bond with a value of 18.0 kcal/mol. The smallest BDE in ClOClO is for the central O-Cl bond with a BDE of 9.7 kcal/mol. Electronic transitions were calculated with the equations of motion EOM-CCSD method. The calculations clearly demonstrate that singlet states of ClClO2 absorb to longer wavelengths in the visible than do the singlet states of ClOOCl as does ClOClO.

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation.

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.

Reactions of diborane with ammonia and ammonia borane: catalytic effects for multiple pathways for hydrogen release.

High-level electronic structure calculations have been used to construct portions of the potential energy surfaces related to the reaction of diborane with ammonia and ammonia borane (B2H6 + NH3 and B2H6 + BH3NH3)to probe the molecular mechanism of H2 release. Geometries of stationary points were optimized at the MP2/aug-cc-pVTZ level. Total energies were computed at the coupled-cluster CCSD(T) theory level with the correlation-consistent basis sets. The results show a wide range of reaction pathways for H2 elimination. The initial interaction of B2H6 + NH3 leads to a weak preassociation complex, from which a B-H-B bridge bond is broken giving rise to a more stable H3BHBH2NH3 adduct. This intermediate, which is also formed from BH3NH3 + BH3, is connected with at least six transition states for H2 release with energies 18-93 kal/mol above the separated reactants. The lowest-lying transition state is a six-member cycle, in which BH3exerts a bifunctional catalytic effect accelerating H2 generation within a B-H-H-N framework. Diborane also induces a catalytic effect for H2 elimination from BH3NH3 via a three-step pathway with cyclic transition states. Following conformational changes, the rate-determining transition state for H2 release is approximately 27 kcal/mol above the B2H6 + BH3NH3 reactants, as compared with an energy barrier of approximately 37 kcal/mol for H2 release from BH3NH3. The behavior of two separated BH3 molecules is more complex and involves multiple reaction pathways. Channels from diborane or borane initially converge to a complex comprising the H3BHBH2NH3adduct plus BH3. The interaction of free BH3 with the BH3 moiety of BH3NH3 via a six-member transition state with diborane type of bonding leads to a lower-energy transition state. The corresponding energy barrier is approximately 8 kcal/mol, relative to the reference point H3BHBH2NH3 adduct + BH3. These transition states are 27-36 kcal/mol above BH3NH3 + B2H6, but 1-9 kcal/mol below the separated reactants BH3NH3 + 2 BH3. Upon chemical activation of B2H6 by forming 2 BH3, there should be sufficient internal energy to undergo spontaneous H2 release. Proceeding in the opposite direction, the H2 regeneration of the products of the B2H6 + BH3NH3reaction should be a feasible process under mild thermal conditions.

Thermochemical parameters of CHFO and CF2O.

Thermochemical properties of CHFO and CF 2O and their derivatives were calculated by using coupled-cluster theory (U)CCSD(T) calculations with the aug-cc-pV nZ ( n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. The predicted properties include the following. Enthalpies of formation (298 K, kcal/mol): Delta H f (CF 2O) = -144.7, Delta H f(CHFO) = -91.1, Delta H f (CFO (*)) = -41.6. Bond dissociation energy (0 K, kcal/mol): BDE(CFO-F) = 120.7, BDE(CHO-F) = 119.1, BDE(CFO-H) = 100.2. Ionization potential (eV): IP 1(CF 2O) = 13.04, IP 2(CF 2O) = 14.09, IP 1(CHFO) = 12.41, IP 2(CHFO) = 13.99, IP 1(CFO (*)) = 9.34. Proton affinity (298 K, kcal/mol), PA O(CF 2O) = 148.8, PA O(CHFO) = 156.7, PA F(CHFO) = 154.5 kcal/mol. Electron affinity: EA(CFO (*)) = 2.38 eV. Triplet-singlet separation gap (eV): Delta E T1-S0(CF 2O) = 4.47, Delta E T1-S0(CHFO) = 4.36. Triplet-triplet transition energy (eV): Delta E T2-T1(CF 2O) = 0.44. The new calculated values contribute to solving some persistent discrepancies in the literature. The effects of F-atoms on thermochemical parameters are not linearly additive, and the changes are largely dominated by the first F-substitution. On the basis of the calculated proton affinities of CF 2O and CF 3OH, the nucleophilicities of the oxygen atoms are, within computational errors, the same in both compounds.

Bond dissociation energies in second-row compounds.

Heats of formation at 0 and 298 K are predicted for PF3, PF5, PF3O, SF2, SF4, SF6, SF2O, SF2O2, and SF4O as well as a number of radicals derived from these stable compounds on the basis of coupled cluster theory [CCSD(T)] calculations extrapolated to the complete basis set limit. In order to achieve near chemical accuracy (+/-1 kcal/mol), additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and vibrational zero-point energies. The calculated values substantially reduce the error limits for these species. A detailed comparison of adiabatic and diabatic bond dissociation energies (BDEs) is made and used to explain trends in the BDEs. Because the adiabatic BDEs of polyatomic molecules represent not only the energy required for breaking a specific bond but also contain any reorganization energies of the bonds in the resulting products, these BDEs can be quite different for each step in the stepwise loss of ligands in binary compounds. For example, the adiabatic BDE for the removal of one fluorine ligand from the very stable closed-shell SF6 molecule to give the unstable SF5 radical is 2.8 times the BDE needed for the removal of one fluorine ligand from the unstable SF5 radical to give the stable closed-shell SF4 molecule. Similarly, the BDE for the removal of one fluorine ligand from the stable closed-shell PF3O molecule to give the unstable PF2O radical is higher than the BDE needed to remove the oxygen atom to give the stable closed-shell PF3 molecule. The same principles govern the BDEs of the phosphorus fluorides and the sulfur oxofluorides. In polyatomic molecules, care must be exercised not to equate BDEs with the bond strengths of given bonds. The measurement of the bond strength or stiffness of a given bond represented by its force constant involves only a small displacement of the atoms near equilibrium and, therefore, does not involve any reorganization energies, i.e., it may be more appropriate to correlate with the diabatic product states.

Energetics and mechanism of the decomposition of trifluoromethanol.

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.

Heats of formation of triplet ethylene, ethylidene, and acetylene.

Heats of formation of the lowest triplet state of ethylene and the ground triplet state of ethylidene have been predicted by high level electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T) energies extrapolated to the complete basis set limit using correlation consistent basis sets (CBS), plus additional corrections predict the following heats of formation in kcal/mol: DeltaH0r(C2H4,3A1) = 80.1 at 0 K and 78.5 at 298 K, and DeltaH0t(CH3CH,3A' ') = 86.8 at 0 K and 85.1 at 298 K, with an error of less than +/-1.0 kcal/mol. The vertical and adiabatic singlet-triplet separation energies of ethylene were calculated as DeltaES-T,vert = 104.1 and DeltaES-T,adia = 65.8 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo (DMC) values of 103.5 +/- 0.3 and 66.4 +/- 0.3 kcal/mol. Both sets of computational values differ from the experimental estimate of 58 +/- 3 kcal/mol for the adiabatic splitting. The computed singlet-triplet gap at 0 K for acetylene is DeltaES-T,adia(C2H2) = 90.5 kcal/mol, which is in notable disagreement with the experimental value of 82.6 kcal/mol. The heat of formation of the triplet is DeltaH0tC2H2,3B2) = 145.3 kcal/mol. There is a systematic underestimation of the singlet-triplet gaps in recent photodecomposition experiments by approximately 7 to 8 kcal/mol. For vinylidene, we predict DeltaH0t(H2CC,1A1) = 98.8 kcal/mol at 298 K (exptl. 100.3 +/- 4.0), DeltaH0t(H2CC,3B2) = 146.2 at 298 K, and an energy gap DeltaES-T-adia(H2CC) = 47.7 kcal/mol.

Coupled-cluster study of the electronic structure and energetics of tetrasulfur, S4.

Ab initio electronic structure calculations are reported for S4. Geometric and energetic parameters are calculated using the singles and doubles coupled-cluster method, including a perturbutional correction for connected triple excitation, CCSD(T), together with systematic sequences of correlation consistent basis sets extrapolated to the complete basis set limit. The geometry for the ground state singlet C2v structure of S4 is in good agreement with the microwave structure determined for S4. There is a low-lying D2h transition state at 1.6 kcal/mol which interchanges the long S-S bond. S4 has a low-lying triplet state (3B 1u) in D2h symmetry which is 10.8 kcal/mol above the C2v singlet ground state. The S-S bond dissociation energy for S4 into two S2(3Sigma*g) molecules is predicted to be 22.8 kcal mol(-1). The S-S bond energy to form S3+S(3P) is predicted to be 64 kcal/mol.

Heats of formation of boron hydride anions and dianions and their ammonium salts [BnHmy-][NH4+]y with y=1-2.

Thermochemical parameters of the closo boron hydride BnHn2- dianions, with n=5-12, the B3H8- and B11H14- anions, and the B5H9 and B10H14 neutral species were predicted by high-level ab initio electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T)/complete basis set (CBS) extrapolated energies, plus additional corrections were used to predict the heats of formation of the simplest BnHmy- species in the gas phase in kcal/mol at 298 K: DeltaHf(B3H8-)=-23.1+/-1.0; DeltaHf(B5H52-)=119.4+/-1.5; DeltaHf(B6H62-)=64.1+/-1.5; and DeltaHf(B5H9)=24.1+/-1.5. The heats of formation of the larger species were evaluated by the G3 method from hydrogenation reactions (values at 298 K, in kcal/mol with estimated error bars of+/-3 kcal/mol): DeltaHf(B7H72-)=51.8; DeltaHf(B8H82-)=46.1; DeltaHf(B9H92-)=24.4; DeltaHf(B10H102-)=-12.5; DeltaHf(B11H112-)=-11.8; DeltaHf(B12H122-)=-86.3; DeltaHf(B11H14-)=-57.3; and DeltaHf(B10H14)=18.7. A linear correlation between atomization energies of the dianions and energies of the BH units was found. The heats of formation of the ammonium salts of the anions and dianions were predicted using lattice energies (UL) calculated from an empirical expression based on ionic volumes. The UL values (0 K) of the BnHn2- dianions range from 319 to 372 kcal/mol. The values of UL for the B3H8- and B11H14- anions are 113 and 135 kcal/mol, respectively. The calculated lattice energies and gas-phase heats of formation of the constituent ions were used to predict the heats of formation of the ammonium crystal salts [BnHmy-][NH4+]y. These results were used to evaluate the thermodynamics of the H2 release reactions from the ammonium hydro-borate salts.