A theoretical study of thermodynamics and kinetics of nitrosamines: a potential no carrier

In this theoretical study, several hybird DFT functionals and MP2 method are used to investigate the properties and the kinetics of a series of nitrosamines. The results show SN or NS transnitrosation reaction to be more favorable via an S_N2-like pathway. The stability is predicted to be in the order of H₂NNO > cis-MeHNNO > trans-MeHNNO > Me₂NNO > trans-PhHNNO > cis-PhHNNO > cis-MeSNO > Ph₂NNO > N-methylenenitrous amide, in which Ph₂NNO and N-methylenenitrous amide will be potential candidates for the NO donor. For N-methylenenitrous amide, which has the strongest NO donating strength among the titled nitroamines, a nearly perpendicular configuration between H₂C=N and NO can plausibly be rationalized by the fact that lone pair of the nitrogen atom on the fragment H₂CN must be π-type, not σ-type, to form a mesomeric effect with π*N-O of the NO group. Using the polarizable continuum model to consider the water solvent effect, all the barriers and endothermicities of the transnitrosation reactions are decreased and the correlated %N–H and %N–S are decreased and increased.     

Applied chemical thermodynamics and kinetics on pharmaceutical compounds

The important role of thermoanalytical methods in the field of pharmaceutical and galenic research is outlined. The thermodynamic stability of polymorphic forms of a substance is discussed.

---

Zusammenfassung Es wird die wichtige Rolle thermoanalytischer Methoden auf dem Gebiet der pharmazeutischen und galenischen Forschung betont und die thermodynamische Stabilität polymorpher Formen von Substanzen besprochen.

---

. .

     

A photo-ionization study of organosulfur ring compounds: Thiirane,thietane and tetrahydrothiophene

The gas phase heats of formation of several organosulfur cations were determined from thiirane, thietane and tetrahydrothiophene precursor molecules by photoionization mass spectrometry. Heats of formation at 0 K and 298 K are reported for the following ions: [H₂CS]^+˙, [H₃CS]⁺, [C₂H₃S]⁺, [C₂H₄S]^+˙, [C₃H₅S]⁺, [C₃H₆S]^+˙, [C₄H₇S]⁺ and [C₄H₈S]^+˙. The [C₄H₇S]⁺ (m/z 87), [C₂H₄S]^+˙ (m/z 60), [C₂H₃S]⁺ (m/z 59), [C₄H₇]⁺ (m/z 55), [C₄H₆]^+˙ (m/z 54) and [CH₂S]^+˙ (m/z 46) ions are produced from metastable tetrahydrothiophene ions at photon energies between 10.2 and 10.7 eV. Decay rates of internal energy selected parent ions to the m/z 60, 59, 55 and 54 fragments were measured by threshold photoelectron-photoion coincidence, the results of which are compared to statistical theory (RRKM/QET) calculations. The [C₂H₄S]^+˙ ion from tetrahydrothiophene is found to have the thioacetaldehyde structure. From the measured [C₂H₄S]^+˙ onset a ΔH = 50±8 kJ mol^?1 was calculated for the thioacetaldehyde molecule.     

A theoretical study on the mechanism and thermodynamics of ozone-water gas phase reaction

Ozone water reaction including a complex was studied at the MP2/6-311++G(d,p) and CCSD/6-311++G(2df,2p)//MP2/6-311++G(d,p) levels of theory. The interaction between water oxygen and central oxygen of ozone produces stable H₂O-O₃ complex with no barrier. With decomposition of this complex through H-abstraction by O₃ and O-abstraction by H₂O, three possible product channels were found. Intrinsic reaction coordinate, topological analyses of atom in molecule, and vibrational frequency calculation have been used to confirm the preferred mechanism. Thermodynamic data at T = 298.15 K and atmospheric pressure have been calculated. The results show that the production of hydrogen peroxide is the main reaction channel with ΔG = ?21.112 kJ mol^-1.     

Modeling the gas-phase thermochemistry of organosulfur compounds

Key to understanding the involvement of organosulfur compounds in a variety of radical chemistries, such as atmospheric chemistry, polymerization, pyrolysis, and so forth, is knowledge of their thermochemical properties. For organosulfur compounds and radicals, thermochemical data are, however, much less well documented than for hydrocarbons. The traditional recourse to the Benson group additivity method offers no solace since only a very limited number of group additivity values (GAVs) is available. In this work, CBS‐QB3 calculations augmented with 1D hindered rotor corrections for 122 organosulfur compounds and 45 organosulfur radicals were used to derive 93 Benson group additivity values, 18 ring‐strain corrections, 2 non‐nearest‐neighbor interactions, and 3 resonance corrections for standard enthalpies of formation, standard molar entropies, and heat capacities for organosulfur compounds and organosulfur radicals. The reported GAVs are consistent with previously reported GAVs for hydrocarbons and hydrocarbon radicals and include 77 contributions, among which 26 radical contributions, which, to the best of our knowledge, have not been reported before. The GAVs allow one to estimate the standard enthalpies of formation at 298 K, the standard entropies at 298 K, and standard heat capacities in the temperature range 300–1500 K for a large set of organosulfur compounds, that is, thiols, thioketons, polysulfides, alkylsulfides, thials, dithioates, and cyclic sulfur compounds. For a validation set of 26 organosulfur compounds, the mean absolute deviation between experimental and group additively modeled enthalpies of formation amounts to 1.9 kJ mol^?1. For an additional set of 14 organosulfur compounds, it was shown that the mean absolute deviations between calculated and group additively modeled standard entropies and heat capacities are restricted to 4 and 2 J mol^?1 K^?1, respectively. As an alternative to Benson GAVs, 26 new hydrogen‐bond increments are reported, which can also be useful for the prediction of radical thermochemistry.     

A theoretical study of the kinetics of the benzylperoxy radical isomerization

The rate constant of the benzylperoxy isomerization reaction has been computed using 54 different levels of theory and has been compared to the experimental value reported at 773 K. The aim of this methodology work is to demonstrate that standard theoretical methods are not adequate to obtain quantitative rate constants for the reaction under study. The use of the elaborated CASPT2 method is essential to estimate a quantitative rate constant. Geometry optimizations and vibrational frequency calculations are performed using three different methods (B3LYP, MPW1K, and MP2) and six different basis sets (6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), and cc-pVDZ). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the 6-31G(d,p) basis set, and with the CASPT2 level of theory with the ANO-L-VDZP basis set. Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. We recommend the use of the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory to compute the temperature dependence of the rate constant of the four-center isomerization of the benzylperoxy radical. It is given by the following relation: k(600-2000 K) (in s (-1)) = (1.29 x 10 (10)) T (0.79) exp[(-133.1 in kJ mol (-1))/ RT]. These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. This computational procedure can be extended to predict rate constants for other similar reactions where no available experimental data exist.     

A theoretical study on the aromaticity of thiadiazoles and related compounds

Three different approaches have been considered to determinate the aromatic character of isomeric thiadiazoles. Several indices have been obtained from calculations made in the context of nonlocal (B3LYP) density functional calculations. The results were compared with the existing experimental evidence for analog heterocycles. It is concluded that the thiadiazoles are aromatic compounds, with a larger aromatic character than pyrrole, thiophene and furan. The monoxide and dioxide derivatives are shown to be non-aromatic or antiaromatic molecules.     

Study of polymorphism of organosulfur and organoselenium compounds

The thermal behavior of a series of organosulfur and organoselenium compounds has been studied by means of differential scanning calorimetry, X-ray diffraction, and thermomicroscopy in order to investigate their polymorphism. In this study, the polymorphism of some of the products has been established. The results of the experiments show that there are four types of thermal behavior for compounds studied. The physicochemical characterization of sulfur and selenium compounds showed differences in structural parameters and fusion temperatures among polymorphic forms.     

Synthesis of highly reactive organosulfur compounds

Application of novel bowl‐type and dendrimer‐type steric protection groups to the first synthesis of stable aromatic S‐nitrosothiols is described. These compounds showed remarkable thermal stability whereas they easily reacted with appropriate reagents. X‐ray crystallographic analysis established their structures, where the C S N O linkage adopts only the syn conformation. Synthesis of a stable sulfenic acid by taking advantage of the bowl‐type substituent is also delineated. © 2002 Wiley Periodicals, Inc. Heteroatom Chem 13:414–418, 2002; Published online in Wiley Interscience (www.interscience.wiley.com). DOI 10.1002/hc.10068     

A theoretical study of malononitrile addition to carbonyl compounds

The nucleophilic addition of the malononitrile anion (MN⁻) to formaldehyde was studied theoretically by the AM1 semiempirical MO method. The addition is found to be endothermic with a late productlike transition state on the reaction coordinate. Additions of MN⁻ to a series of carbonyl compounds were studied in order to investigate the substituent effect on the energetics of the title addition and the nucleophilic attack reactivity. The solvent effect was stimulated by hydrogen bonding a single molecule of water to the formaldehyde oxygen and/or to the MN⁻ anion. Its influence on the energetics and the transition-state geometry was estimated. The Hammond postulate was satisfied for the studied additions. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62 : 419–426, 1997     

Chemical thermodynamics and theoretical models

This paper presents a brief discussion of the relationship between classical thermodynamics and theoretical models of chemical systems. Applications of thermodynamics to real chemical systems and to related theoretical models concerning “complexes in nonelectrolyte solutions” and “substituent and solvent effects in organic chemistry” are described, with particular emphasis on research that my colleagues and I have done in these areas. This paper concludes with a brief autobiographical section on “people and places”.     

Theoretical study of the thermodynamics and kinetics of hydrogen abstractions from hydrocarbons

Thermochemical and kinetic data were calculated at four cost-effective levels of theory for a set consisting of five hydrogen abstraction reactions between hydrocarbons for which experimental data are available. The selection of a reliable, yet cost-effective method to study this type of reactions for a broad range of applications was done on the basis of comparison with experimental data or with results obtained from computationally demanding high level of theory calculations. For this benchmark study two composite methods (CBS-QB3 and G3B3) and two density functional theory (DFT) methods, MPW1PW91/6-311G(2d,d,p) and BMK/6-311G(2d,d,p), were selected. All four methods succeeded well in describing the thermochemical properties of the five studied hydrogen abstraction reactions. High-level Weizmann-1 (W1) calculations indicated that CBS-QB3 succeeds in predicting the most accurate reaction barrier for the hydrogen abstraction of methane by methyl but tends to underestimate the reaction barriers for reactions where spin contamination is observed in the transition state. Experimental rate coefficients were most accurately predicted with CBS-QB3. Therefore, CBS-QB3 was selected to investigate the influence of both the 1D hindered internal rotor treatment about the forming bond (1D-HR) and tunneling on the rate coefficients for a set of 21 hydrogen abstraction reactions. Three zero curvature tunneling (ZCT) methods were evaluated (Wigner, Skodje & Truhlar, Eckart). As the computationally more demanding centrifugal dominant small curvature semiclassical (CD-SCS) tunneling method did not yield significantly better agreement with experiment compared to the ZCT methods, CD-SCS tunneling contributions were only assessed for the hydrogen abstractions by methyl from methane and ethane. The best agreement with experimental rate coefficients was found when Eckart tunneling and 1D-HR corrections were applied. A mean deviation of a factor 6 on the rate coefficients is found for the complete set of 21 reactions at temperatures ranging from 298 to 1000 K. Tunneling corrections play a critical role in obtaining accurate rate coefficients, especially at lower temperatures, whereas the hindered rotor treatment only improves the agreement with experiment in the high-temperature range.     

Synthesis of new types of organosulfur nitro compounds

Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 5, pp. 1237–1238, May, 1991.     

Dynamics and kinetics of the S + HO2 reaction: A theoretical study

We report a quasi‐classical trajectory study of the S + HO₂ reaction using a previously reported global potential energy surface for the ground electronic state of HSO₂. Zero‐point energy leakage is approximately accounted for by using the vibrational energy quantum mechanical threshold method. Calculations are carried out both for specific ro‐vibrational states of the reactants and thermalized ones, with rate constants being reported as a function of temperature. The results suggest that the title reaction is capture type, with OH and SO showing as the most favorable products. The internal energy distribution of such products and the reaction mechanism are also investigated. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 40: 533–540, 2008     

Estimation of the activity of catalysts and the reactivity of organosulfur compounds

The kinetics of hydrogenolysis of various classes of organosulfur compounds were studied in the presence of MeAlCl ₄ catalysts (Me=Li, Na, K). The reactions were conducted at 100–250°C without introduction of hydrogen into the reaction zone, but in the presence of hydrocarbons capable of acting as hydride—ion donors; these donors were necessary for H ₂ S formation. A series of catalytic activities and optimal reaction conditions for each catalyst were determined from the rate constants of hydrogenolysis at the C-S bond. Reactivities of mercaptans, sulfides, disulfides, and thiophene were estimated.Institute of Chemistry, Bashkir Scientific Center, Ural Branch, Russian Academy of Sciences, Ufa 450054. Translated from Izvestiya Akademii Nauk, Seriya Khimicheskaya, No. 3, pp. 562–564, March, 1992.     

A theoretical study of the mechanism and kinetics of F+N3 reactions.

Both the singlet(1A') and triplet(3A') potential energy surfaces (PESs) of F+N(3) reactions are investigated using the complete-active-space self-consistent field (CASSCF) and the multireference configuration interaction (MRCI) methods with a proper active space. The minimum energy crossing point (MECP) at the intersection seam between the 1A' and 3A' PESs is located and used to clarify the reaction mechanisms. Two triplet transition states are found, with one in the cis form and the other one in the trans form. Further kinetic calculations are performed with the canonical unified statistical (CUS) theory on the singlet PES and the improved canonical variational transition-state (ICVT) method on the triplet PES. The rate constants are also reported. At 298 K, the calculated rate constant is in reasonably good agreement with experimental values, and spin-orbit coupling effects lower it by 28 %. The spectroscopic constants derived from the fitted potential-energy curves for the singlet and triplet states of NF are in very good agreement with experimental values. Our calculations indicate that the adiabatic reaction on the singlet PES leading to NF(a(1)Delta)+N(2) is the major channel, whereas the nonadiabatic reaction through the MECP, which leads to NF(X(3)Sigma(-))+N(2), is a minor channel.