Complexation of Phenol and Thiophenol by Amine N‐Oxides: Isothermal Titration Calorimetry and ab Initio Calculations

To develop a new solvent‐impregnated resin (SIR) system for removal of phenols from water, the complex formation of dimethyldodecylamine N‐oxide (DMDAO), trioctylamine N‐oxide (TOAO), and tris(2‐ethylhexyl)amine N‐oxide (TEHAO) with phenol (PhOH) and thiophenol (PhSH) is studied. To this end we use isothermal titration calorimetry (ITC) and quantum chemical modeling (on B3LYP/6‐311G(d,p)‐optimized geometries: B3LYP/6‐311+G(d,p), B3LYP/6‐311++G(2d,2p), MP2/6‐311+G(d,p), and spin component scaled (SCS) MP2/6‐311+G(d,p); M06‐2X/6‐311+G(d,p)//M06‐2X/6‐311G(d,p), MP2 with an extrapolation to the complete basis set limit (MP2/CBS), as well as CBS‐Q). The complexes are analyzed in terms of structural (e.g., bond lengths) and electronic elements (e.g., charges). Furthermore, complexation and solvent effects (in benzene, toluene, and mesitylene) are investigated by ITC measurements, yielding binding constants K, enthalpies ΔH⁰, Gibbs fre energies ΔG⁰, and entropies ΔS⁰ of complex formation, and stoichiometry N. The ITC measurements revealed strong 1:1 complex formation between both DMDAO–PhOH and TOAO–PhOH. The binding constant (K=1.7–5.7×10⁴ M ^?1) drops markedly when water‐saturated toluene was used (K=5.8×10³ M ^?1), and π–π interaction with the solvent is shown to be relevant. Quantum mechanical modeling confirms formation of stable 1:1 complexes with linear hydrogen bonds that weaken on attachment of electron‐withdrawing groups to the amine N‐oxide moiety. Modeling also showed that complexes with PhSH are much weaker than those with PhOH, and in fact too weak for ITC determination. CBS‐Q incorrectly predicts equal or even higher binding enthalpies for PhSH than for PhOH, which invalidates it as a benchmark for other calculations. Data from the straightforward SCS‐MP2 method without counterpoise correction show very good agreement with the MP2/CBS values.     

Preparing Gold(I) for Interactions with Proton Donors: The Elusive [Au]⋅⋅⋅HO Hydrogen Bond

MP2 and DFT calculations with correlation consistent basis sets indicate that isolated linear anionic dialkylgold(I) complexes form moderately strong (ca. 10 kcal mol^?1) Au???H hydrogen bonds with single H₂O molecules as donors in the absence of sterically demanding substituents. Relativistic effects are critically important in the attraction. Such bonds are significantly weaker in neutral, strong σ‐donor N‐heterocyclic carbene (NHC) complexes (ca. 5 kcal mol^?1). The overall association (>11 kcal mol^?1), however, is strengthened by co‐operative, synergistic classical hydrogen bonding when the NHC ligands bear NH units. Further manipulation of the interaction by ligands positioned trans to the carbene, is possible.     

Exploring new species on the [H,S, Se,Cl] potential energy surface

This investigation is concerned with the characterization of seleno‐sulfide‐halogen model systems, the isomerization processes, and the dissociation into diatomic fragment channels on the [H, S, Se, Cl] potential energy surface. Structural, energetic, and vibrational data were obtained at the CCSD(T) and MP2 levels of theory with the series of correlation consistent basis sets and extrapolated to the complete basis set (CBS) limit. For the frequencies, additional computations were performed to include the contribution of anharmonic effects, and for the determination of the heats of formation, important corrections incorporating core‐valence correlation effects and relativistic effects (scalar and spin‐orbit) were taken into account. CCSD(T)/CBS relative stability (kcal mol^?1) follows the order: HSSeCl (0.0), HSeSCl (8.80), SSeHCl (23.52), and SeSHCl (25.87). The cis‐rotational barrier for the two lowest isomers is practically identical (10.14 and 10.09 kcal mol^?1), whereas for the trans barrier, we obtained 9.25 (HSSeCl) and 8.45 (HSeSCl) kcal mol^?1. Dissociation of HSSeCl (HSeSCl) into HS (HSe) + SeCl (SCl) requires 59.70 (56.30) kcal mol^?1. For the most stable isomer, we predict a value of the heat of formation at 298.15 K of 2.53 kcal mol^?1. One of the outcomes of this research is that the MP2 results are consistent with those of CCSD(T). The MP2 method turns out to be a reliable alternative for a first exploration of larger catenated species, although it accounts for a lesser fraction of correlation effects. © 2012 Wiley Periodicals, Inc.     

Accurate reaction barrier heights of pericyclic reactions: Surprisingly large deviations for the CBS‐QB3 composite method and their consequences in DFT benchmark studies

Accurate barrier heights are obtained for the 26 pericyclic reactions in the BHPERI dataset by means of the high‐level Wn‐F12 thermochemical protocols. Very often, the complete basis set (CBS)‐type composite methods are used in similar situations, but herein it is shown that they in fact result in surprisingly large errors with root mean square deviations (RMSDs) of about 2.5 kcal mol^?1. In comparison, other composite methods, particularly G4‐type and estimated coupled cluster with singles, doubles, and quasiperturbative triple excitations [CCSD(T)/CBS] approaches, show deviations well below the chemical‐accuracy threshold of 1 kcal mol^?1. With the exception of SCS‐MP2 and the herein newly introduced MP3.5 approach, all other tested Møller‐Plesset perturbative procedures give poor performance with RMSDs of up to 8.0 kcal mol^?1. The finding that CBS‐type methods fail for barrier heights of these reactions is unexpected and it is particularly troublesome given that they are often used to obtain reference values for benchmark studies. Significant differences are identified in the interpretation and final ranking of density functional theory (DFT) methods when using the original CBS‐QB3 rather than the new Wn‐F12 reference values for BHPERI. In particular, it is observed that the more accurate Wn‐F12 benchmark results in lower statistical errors for those methods that are generally considered to be robust and accurate. Two examples are the PW6B95‐D3(BJ) hybrid‐meta‐general‐gradient approximation and the PWPB95‐D3(BJ) double‐hybrid functionals, which result in the lowest RMSDs of the entire DFT study (1.3 and 1.0 kcal mol^?1, respectively). These results indicate that CBS‐QB3 should be applied with caution in computational modeling and benchmark studies involving related systems. © 2015 Wiley Periodicals, Inc.     

Formic acid catalyzed gas-phase reaction of H2O with SO3 and the reverse reaction: a theoretical study

The formic acid catalyzed gas‐phase reaction between H₂O and SO₃ and its reverse reaction are respectively investigated by means of quantum chemical calculations at the CCSD(T)//B3LYP/cc‐pv(T+d)z and CCSD(T)//MP2/aug‐cc‐pv(T+d)z levels of theory. Remarkably, the activation energy relative to the reactants for the reaction of H₂O with SO₃ is lowered through formic acid catalysis from 15.97 kcal mol^?1 to ?15.12 and ?14.83 kcal mol^?1 for the formed H₂O ??? SO₃ complex plus HCOOH and the formed H₂O ??? HCOOH complex plus SO₃, respectively, at the CCSD(T)//MP2/aug‐cc‐pv(T+d)z level. For the reverse reaction, the energy barrier for decomposition of sulfuric acid is reduced to ?3.07 kcal mol^?1 from 35.82 kcal mol^?1 with the aid of formic acid. The results show that formic acid plays a strong catalytic role in facilitating the formation and decomposition of sulfuric acid. The rate constant of the SO₃+H₂O reaction with formic acid is 10⁵ times greater than that of the corresponding reaction with water dimer. The calculated rate constant for the HCOOH+H₂SO₄ reaction is about 10^?13 cm³ molecule^?1 s^?1 in the temperature range 200–280 K. The results of the present investigation show that formic acid plays a crucial role in the cycle between SO₃ and H₂SO₄ in atmospheric chemistry.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

Oxidative Dehydrogenation of Propane over a VO2‐Exchanged MCM‐22 Zeolite: A DFT Study

The adsorption and the mechanism of the oxidative dehydrogenation (ODH) of propane over VO₂‐exchanged MCM‐22 are investigated by DFT calculations using the M06‐L functional, which takes into account dispersion contributions to the energy. The adsorption energies of propane are in good agreement with those from computationally much more demanding MP2 calculations and with experimental results. In contrast, B3LYP binding energies are too small. The reaction begins with the movement of a methylene hydrogen atom to the oxygen atom of the VO₂ group, which leads to an isopropyl radical bound to a HO? V? O intermediate. This step is rate determining with the apparent activation energy of 30.9 kcal mol^?1, a value within the range of experimental results for ODH over other silica supports. In the propene formation step, the hydroxyl group is the more reactive group requiring an apparent activation energy of 27.7 kcal mol^?1 compared to that of the oxy group of 40.8 kcal mol^?1. To take the effect of the extended framework into account, single‐point calculations on 120T structures at the same level of theory are performed. The apparent activation energy is reduced to 28.5 kcal mol^?1 by a stabilizing effect caused by the framework. Reoxidation of the catalyst is found to be important for the product release at the end of the reaction.     

C?X⋅⋅⋅O Halogen Bonding: Interactions of Trifluoromethyl Halides with Dimethyl Ether

The formation of weakly bound molecular complexes between dimethyl ether (DME) and the trifluoromethyl halides CF₃Cl, CF₃Br and CF₃I dissolved in liquid argon and in liquid krypton is investigated, using Raman and FTIR spectroscopy. For all halides evidence is found for the formation of C? X???O halogen‐bonded 1:1 complexes. At higher concentrations of CF₃Br, a weak absorption due to a 1:2 complex is also observed. Using spectra recorded at temperatures between 87 and 125 K, the complexation enthalpies for the complexes are determined to be ?6.8(3) kJ mol^?1 (DME?CF₃Cl), ?10.2(1) kJ mol^?1 (DME?CF₃Br), ?15.5(1) kJ mol^?1 (DME?CF₃I), and ?17.8(5) kJ mol^?1 [DME(?CF₃Br)₂]. Structural and spectral information on the complexes is obtained from ab initio calculations at the MP2/ 6‐311++G(d,p) and MP2/6‐311++G(d,p)+LanL2DZ* levels. By applying Monte Carlo free energy perturbation calculations to account for the solvent influences, and statistical thermodynamics to estimate the zero‐point vibrational and thermal influences, the ab initio complexation energies are converted into complexation enthalpies for the solutions in liquid argon. The resulting values are compared with the experimental data deduced from the cryosolutions.     

The thermodynamic stability of P8, A CBS‐Q study

The thermodynamic stabilities of P₂, P₄, and three P₈ cage structure were investigated through high‐precision CBS‐Q calculations. The CBS‐Q values for the bond energy of P₂ (ΔE_o: +115.7 kcal mol⁻¹) and the formation of P₄ from P₂ (Δ E_o:‐56.6 kcal mol⁻¹) were in excellent agreement with the experimental values (E_o: +117 and ‐56.4 kcal mol⁻¹ respectively). Among the P₈ cages, the cubane structure was the least stable (Δ E_o +37 kcal vs. 2×P₄). The most stable P₈ isomer adopts a cuneane structure resembling S₄N₄, and is more stable than white phosphorus at T = 0 K (Δ E_o −3.3 kcal mol⁻¹), but still unstable under standard conditions for entropic reasons (Δ G^o of +8.1 kcal mol⁻¹ vs. 2×P₄). The CBS‐Q energies represent significant revisions (6–20 kcal mol⁻¹) of previous computational predictions obtained by high‐level single method calculations. © 2005 Wiley Periodicals, Inc. Heteroatom Chem 16:453–457, 2005; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20119     

Anionic CH⋅⋅⋅X− Hydrogen Bonds: Origin of Their Strength,Geometry, and Other Properties

CF₃H as a proton donor was paired with a variety of anions, and its properties were assessed by MP2/aug‐cc‐pVDZ calculations. The binding energy of monoanions halide, NO₃^?, formate, acetate, HSO₄^?, and H₂PO₄^? lie in the 12–17 kcal mol^?1 range, although F^? is more strongly bound, by 26 kcal mol^?1. Dianions SO₄^2? and HPO₄^2? are bound by 27 kcal mol^?1, and trianion PO₄^3? by 45 kcal mol^?1. When two O atoms are available on the anion, the CH???O^? H‐bond (HB) is usually bifurcated, although asymmetrically. The CH bond is elongated and its stretching frequency redshifted in these ionic HBs, but the shift is reduced in the bifurcated structures. Slightly more than half of the binding energy is attributed to Coulombic attraction, with smaller contributions from induction and dispersion. The amount of charge transfer from the anions to the σ*(CH) orbital correlates with many of the other indicators of bond strength, such as binding energy, CH bond stretch, CH redshift, downfield NMR spectroscopic chemical shift of the bridging proton, and density at bond critical points.     

Enthalpy of formation and bond energies on unsaturated oxygenated hydrocarbons using G3MP2B3 calculation methods

Enthalpies of unsaturated oxygenated hydrocarbons and radicals corresponding to the loss of hydrogen atoms from the parent molecules are intermediates and decomposition products in the oxidation and combustion of aromatic and polyaromatic species. Enthalpies (Δ_fH⁰₂₉₈) are calculated for a set of 27 oxygenated and nonoxygenated, unsaturated hydrocarbons and 12 radicals at the G3MP2B3 level of theory and with the commonly used B3LYP/6‐311g(d,p) density functional theory (DFT) method. Standard enthalpies of formation (Δ_fH⁰₂₉₈) are determined from the calculated enthalpy of reaction (ΔH⁰_rxn,298) using isodesmic work reactions with reference species that have accurately known Δ_fH⁰₂₉₈ values. The deviation between G3MP2B3 and B3LYP methods is under ±0.5 kcal mol^?1 for 9 species, 18 other species differs by less than ±1 kcal mol^?1 , and 11 species differ by about 1.5 kcal mol^?1. Under them are 11 radicals derived from the above‐oxygenated hydrocarbons that show good agreement between G3MP2B3 and B3LYP methods. G3 calculations have been performed to further validate enthalpy values, where a discrepancy of more than 2.5 kcal mol^?1 exists between the G3MP3B3 and density functional results. Surprisingly the G3 calculations support the density functional calculations in these several nonagreement cases. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 633–648, 2005     

The Affinity of Some Lewis Bases for Hexafluoroisopropanol as a Reference Lewis Acid: An ITC/DFT Study

To figure out the possible role of 1,1,1,3,3,3-hexafluoropropan-2-ol (HFIP) as well as to provide reference thermochemical data in solution, the formation of Lewis acid-base complexes between HFIP (Lewis acid) and a series of 8 different Lewis bases (3 sulfoxides, 3 Nsp² pyridine derivatives, 1 aromatic amine, 1 cyclic aliphatic ether) was examined by isothermal titration calorimetry (ITC) experiments and static density functional theory augmented with Dispersion (DFT−D) calculations. Measured ITC association enthalpy values (ΔH_a) lie between −9.3 and −14 kcal mol⁻¹. Computations including a PCM implicit solvation model produced similar exothermicity of association of all studied systems compared to the ITC data with ΔH_a values ranging from −8.5 to −12.7 kcal mol⁻¹. An additional set of calculations combining implicit and explicit solvation by chlorobenzene of the reactants, pointed out the relatively low interference of the solvent with the HFIP-base complexation: its main effect is to slightly enhance the Gibbs energy of the HFIP-Lewis base association. It is speculated that the interactions of bulk HFIP with Lewis bases therefore may significantly intervene in catalytic processes not only via the dynamic microstructuring of the medium but also more explicitly by affecting bonds’ polarization at the Lewis bases.     

Modulation of Stacking Interactions by Transition‐Metal Coordination: Ab Initio Benchmark Studies

A series of ab initio calculations are used to determine the C? H???π and π???π‐stacking interactions of aromatic rings coordinated to transition‐metal centres. Two model complexes have been employed, namely, ferrocene and chromium benzene tricarbonyl. Benchmark data obtained from extrapolation of MP2 energies to the basis set limit, coupled with CCSD(T) correction, indicate that coordinated aromatic rings are slightly weaker hydrogen‐bond acceptors but are significantly stronger hydrogen‐bond donors than uncomplexed rings. It is found that π???π stacking to a second benzene is stronger than in the free benzene dimer, especially in the chromium case. This is assigned, by use of energy partitioning in the local correlation method, to dispersion interactions between metal d and benzene π orbitals. The benchmark data is also used to test the performance of more efficient theoretical methods, indicating that spin‐component scaling of MP2 energies performs well in all cases, whereas various density functionals describe some complexes well, but others with errors of more than 1 kcal mol^?1.     

First Principles (DFT) Characterization of RhI/dppp‐Catalyzed CH Activation by Tandem 1,2‐Addition/1,4‐Rh Shift Reactions of Norbornene to Phenylboronic Acid

The C?H activation in the tandem, “merry‐go‐round”, [(dppp)Rh]‐catalyzed (dppp=1,3‐bis(diphenylphosphino)propane), four‐fold addition of norborene to PhB(OH)₂ has been postulated to occur by a C(alkyl)?H oxidative addition to square‐pyramidal Rh^III?H species, which in turn undergoes a C(aryl)?H reductive elimination. Our DFT calculations confirm the Rh^I/Rh^III mechanism. At the IEFPCM(toluene, 373.15 K)/PBE0/DGDZVP level of theory, the oxidative addition barrier was calculated to be 12.9 kcal mol^?1, and that of reductive elimination was 5.0 kcal mol^?1. The observed selectivity of the reaction correlates well with the relative energy barriers of the cycle steps. The higher barrier (20.9 kcal mol^?1) for norbornyl–Rh protonation ensures that the reaction is steered towards the 1,4‐shift (total barrier of 16.3 kcal mol^?1), acting as an equilibration shuttle. The carborhodation (13.2 kcal mol^?1) proceeds through a lower barrier than the protonation (16.7 kcal mol^?1) of the rearranged aryl–Rh species in the absence of o‐ or m‐substituents, ensuring multiple carborhodations take place. However, for 2,5‐dimethylphenyl, which was used as a model substrate, the barrier for carborhodation is increased to 19.4 kcal mol^?1, explaining the observed termination of the reaction at 1,2,3,4‐tetra(exo‐norborn‐2‐yl)benzene. Finally, calculations with (Z)‐2‐butene gave a carborhodation barrier of 20.2 kcal mol^?1, suggesting that carborhodation of non‐strained, open‐chain substrates would be disfavored relative to protonation.     

Protonation of phenylboronic acid: Comparison of G3B3 and G2MP2 methods

G3B3 and G2MP2 calculations using Gaussian 03 have been carried out to investigate the protonation preferences for phenylboronic acid. All nine heavy atoms have been protonated in turn. With both methodologies, the two lowest protonation energies are obtained with the proton located either at the ipso carbon atom or at a hydroxyl oxygen atom. Within the G3B3 formalism, the lowest‐energy configuration by 4.3 kcal · mol^?1 is found when the proton is located at the ipso carbon, rather than at the electronegative oxygen atom. In the resulting structure, the phenyl ring has lost a significant amount of aromaticity. By contrast, calculations with G2MP2 show that protonation at the hydroxyl oxygen atom is favored by 7.7 kcal · mol^?1. Calculations using the polarizable continuum model (PCM) solvent method also give preference to protonation at the oxygen atom when water is used as the solvent. The preference for protonation at the ipso carbon found by the more accurate G3B3 method is unexpected and its implications in Suzuki coupling are discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Dimers of N‐Heterocyclic Carbene Copper,Silver, and Gold Halides: Probing Metallophilic Interactions through Electron Density Based Concepts

Homobimetallic metallophilic interactions between copper, silver, and gold‐based [(NHC)MX]‐type complexes (NHC=N‐heterocyclic carbene, i.e, 1,3,4‐trimethyl‐4,5‐dihydro‐1H‐1,2,4‐triazol‐5‐ylidene; X=F, Cl, Br, I) were investigated by means of ab initio interaction energies, Ziegler–Rauk‐type energy‐decomposition analysis, the natural orbital for chemical valence (NOCV) framework, and the noncovalent interaction (NCI) index. It was found that the dimers of these complexes predominantly adopt a head‐to‐tail arrangement with typical M ??? M distance of 3.04–3.64 Å, in good agreement with the experimental X‐ray structure determined for [{(NHC)AuCl}₂], which has an Au ??? Au distance of 3.33 Å. The interaction energies between silver‐ and gold‐based monomers are calculated to be about ?25 kcal mol^?1, whereas that for the Cu congener is significantly lower (?19.7 kcal mol^?1). With the inclusion of thermal and solvent contributions, both of which are destabilizing, by about 15 and 8 kcal mol^?1, respectively, an equilibrium process is predicted for the formation of dimer complexes. Energy‐decomposition analysis revealed a dominant electrostatic contribution to the interaction energy, besides significantly stabilizing dispersion and orbital interactions. This electrostatic contribution is rationalized by NHC(δ⁺) ??? halogen(δ^?) interactions between monomers, as demonstrated by electrostatic potentials and derived charges. The dominant NOCV orbital indicates weakening of the π backdonation in the monomers on dimer formation, whereas the second most dominant NOCV represents an electron‐density deformation according to the formation of a very weak M ??? M bond. One of the characteristic signals found in the reduced density gradient versus electron density diagram corresponds to the noncovalent interactions between the metal centers of the monomers in the NCI plots, which is the manifestation of metallophilic interaction.     

A Highly Oxidizing and Isolable Oxoruthenium(V) Complex [RuV(N4O)(O)]2+: Electronic Structure,Redox Properties,and Oxidation Reactions Investigated by DFT Calculations

The electronic structure and redox properties of the highly oxidizing, isolable Ru^V?O complex [Ru^V(N₄O)(O)]²⁺, its oxidation reactions with saturated alkanes (cyclohexane and methane) and inorganic substrates (hydrochloric acid and water), and its intermolecular coupling reaction have been examined by DFT calculations. The oxidation reactions with cyclohexane and methane proceed through hydrogen atom transfer in a transition state with a calculated free energy barrier of 10.8 and 23.8 kcal mol^?1, respectively. The overall free energy activation barrier (ΔG^≠=25.5 kcal mol^?1) of oxidation of hydrochloric acid can be decomposed into two parts: the formation of [Ru^III(N₄O)(HOCl)]²⁺ (ΔG=15.0 kcal mol^?1) and the substitution of HOCl by a water molecule (ΔG^≠=10.5 kcal mol^?1). For water oxidation, nucleophilic attack on Ru^V?O by water, leading to O? O bond formation, has a free energy barrier of 24.0 kcal mol^?1, the major component of which comes from the cleavage of the H? OH bond of water. Intermolecular self‐coupling of two molecules of [Ru^V(N₄O)(O)]²⁺ leads to the [(N₄O)Ru^IV? O₂? Ru^III(N₄O)]⁴⁺ complex with a calculated free energy barrier of 12.0 kcal mol^?1.     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Stabilization of HHeF by Complexation: Is it a Really Viable Strategy?

Ab initio calculations at the MP2 and CCSD(T) levels of theory have disclosed the conceivable existence of fluorine‐coordinated complexes of HHeF with alkali‐metal ions and molecules M⁺ (M⁺=Li⁺–Cs⁺), M⁺–OH₂, M⁺–NH₃ (M⁺=Li⁺, Na⁺), and MX (M=Li, Na; X=F, Cl, Br). All these ligands L induce a shortening of the H? He distance and a lengthening of the He? F distance accompanied by consistent blue‐ and redshifts, respectively, of the H? He and He? F stretching modes. These structural effects are qualitatively similar to those predicted for other investigated complexes of the noble gas hydrides HNgY, but are quantitatively more pronounced. For example, the blueshifts of the H? He stretching mode are exceptionally large, ranging between around 750 and 1000 cm^?1. The interactions of HHeF with the ligands investigated herein also enhance the (HHe)⁺F^? dipole character and produce large complexation energies of around 20–60 kcal mol^?1. Most of the HHeF–L complexes are indeed so stable that the three‐body dissociation of HHeF into H+He+F, exothermic by around 25–30 kcal mol^?1, becomes endothermic. This effect is, however, accompanied by a strong decrease in the H? He? F bending barrier. The complexation energies, ΔE, and the bending barriers, E*, are, in particular, related by the inverse relationship E*(kcal mol^?1)=6.9exp[?0.041ΔE(kcal mol^?1)]. Therefore the HHeF? L complexes, which are definitely stable with respect to H+He+F+L (ΔE≈25–30 kcal mol^?1), are predicted to have bending barriers of only 0.5–2 kcal mol^?1. Overall, our calculations cast doubt on the conceivable stabilization of HHeF by complexation.     

A Theoretical Study of the Kinetics of the Hydrogen Atom Abstraction Reactions from Cyclopropane by H,O (3P), and Cl (2P3/2) Atoms and OH Radicals

The rate constants of the H‐abstraction reactions from cyclopropane by H, O (³P), Cl (²P_3/2), and OH radicals have been calculated over the temperature range of 250?2500 K using two different levels of theory. Calculations of optimized geometrical parameters and vibrational frequencies are performed using the MP2 method combined with the cc‐pVTZ basis set and the 6–311++G(d,p) basis set. Single‐point energy calculations have been carried out with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (perturbatively) electron excitations CCSD(T) using either the cc‐pVTZ, aug‐cc‐pVTZ, and aug‐cc‐pVQZ basis sets or the 6–311++G(3df,3pd) basis set. The CCSD(T) calculated potential energies have been extrapolated to the complete basis limit (CBS) limit. The Full Configuration Interaction (FCI) energies have been also estimated using the continued‐fraction approximation as proposed by Goodson (J. Chem. Phys., 2002, 116, 6948–6956). Canonical transition‐state theory combined with an Eckart tunneling correction has been used to predict the rate constants as a function of temperature using two kinetic models (direct abstraction or complex mechanism) at two levels of theory (CCSD(T)‐cf/CBS//MP2/cc‐pVTZ and CCSD(T)‐cf/6–311++G(3df,3pd)//MP2/6–311++G(d,p)). The calculated kinetic parameters are in reasonable agreement with their literature counterparts for all reactions. In the light of these trends, the use of the Pople‐style basis sets for studying the reactivity of other systems such as larger cycloalkanes or halogenated cycloalkanes is recommended because the 6–311++G(3df,3pd) basis set is less time consuming than the aug‐cc‐pVQZ basis set. Based on our calculations performed at the CCSD(T)‐cf/CBS//MP2/cc‐pVTZ level of theory, the standard enthalpy of formation at 298 K for the cyclopropyl radical has been reassessed and its value is (290.5 ± 1.6) kJ mol^?1.