DFT and ab initio calculations on the reaction between fluorine atoms and the fire suppressant, 2-H heptafluoropropane

Benchmark values for the reaction enthalpy (298 K) and the barrier height (0 K) of the reaction, CF₃CHFCF₃ + F → CF₃CFCF₃ + HF, have been calculated at state-of-the-art ab initio level to be −34.7 ± 1.0 and −0.9 ± 0.9 kcal/mole, respectively. The B3LYP, BH&HLYP, BB1K, MPW1K, MPWB1K and TPSS1KCIS functionals and the model method, the integrated molecular orbital + molecular orbital (IMOMO) method, have also been used to study this reaction.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

Kinetic mechanism of the hydrogen abstraction reactions of the chlorine atoms with CH3CF2Cl and CH3CFCl2: a dual level direct dynamics study

The mechanisms of the reactions: CH(3)CFCl(2) + Cl (R1) and CH(3)CF(2)Cl + Cl (R2) are studied over a wide temperature range (200-3000 K) using the dual-level direct dynamics method. The minimum energy path calculation is carried out at the MP2/6-311G(d,p) and B3LYP/6-311G(d,p) levels, and energetic information is further refined by the G3(MP2) theory. The H-abstraction from the out-of-plane for (R1) is the major reaction channel, while the in-plane H-abstraction is the predominant route of (R2). The canonical variational transition-state theory (CVT) with the small-curvature tunneling (SCT) correction method is used to calculate the rate constants. Using group-balanced isodesmic reactions and hydrogenation reactions as working chemical reactions, the standard enthalpies of formation for CH(3)CFCl(2), CH(3)CF(2)Cl, CH(2)CFCl(2), and CH(2)CF(2)Cl are evaluated at the CCSD(T)/6-311 + G(3df,2p)//MP2/6-311G(d,p) level of theory. The results indicate that the substitution of fluorine atom for the chlorine atom leads to a decrease in the C-H bond reactivity with a small increase in reaction enthalpies. Also, for all reaction pathways the variational effect is small and the SCT effect is only important in the lower temperature range on the rate constants.     

The singlet-triplet separation in CF2: state-of-the-art ab initio calculations and Franck-Condon simulations including anharmonicity.

Geometrical parameters, vibrational frequencies and relative electronic energies of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been calculated. Core-electron effects on the computed minimum-energy geometries and relative electronic energies have been investigated, and relativistic contributions to the computed relative electronic energies calculated. Potential energy functions of the X2B1 state of CF2- and the X1A1 and ?3B1 states of CF2 have been determined, and anharmonic vibrational wavefunctions of these states calculated variationally. Franck-Condon factors including anharmonicity and Duschinsky rotation have been computed and used to simulate the ?-X emission spectrum of CF2 determined by S. Koda [Chem. Phys. Lett. 1978, 55, 353] and the 364 nm laser photodetachment spectrum of CF2- obtained by R. L. Schwartz et al. [J. Phys. Chem. A 1999, 103, 8213]. Comparison between theory and experiment shows that the theoretical approach benchmarked in the present study is able to give highly reliable positions for the CF2(X1A1) + e <-- CF2-(X2B1) and CF2(?3B1) + e <-- CF2-(X2B1) bands in the photoelectron spectrum of CF2- and a reliable singlet-triplet gap for CF2. It is therefore concluded that the same theoretical approach should give reliable simulated CCl2(X1A1) + e <-- CCl2-(X2B1) and CCl2(?3B1) + e <-- CCl2-(X2B1) bands in the photodetachment spectrum of CCl2- and a reliable singlet-triplet gap for CCl2.     

A theoretical ab initio study on the H(2)NO + O(3) reaction

The deviation of the NH(2) pseudo-first-order decay Arrhenius plots of the NH(2) + O(3) reaction at high ozone pressures measured by experimentalists, has been attributed to the regeneration of NH(2) radicals due to the subsequent reactions of the products of this reaction with ozone. Although these products have not yet been characterized experimentally, the radical H(2)NO has been postulated, because it can regenerate NH(2) radicals through the reactions: H(2)NO + O(3) --> NH(2) + O(2) and H(2)NO + O(3) --> HNO + OH + O(2). With the purpose of providing a reasonable explanation from a theoretical point of view to the kinetic observed behaviour of the NH(2) + O(3) system, we have carried ab initio electronic structure calculations on both H(2)NO + O(3) possible reactions. The results obtained in this article, however, predict that of both reactions proposed, only the H(2)NO + O(3) --> NH(2) + O(2) reaction would regenerate indeed NH(2) radicals, explaining thus the deviation of the NH(2) pseudo-first-order decay observed experimentally.     

A comparative study of DFT and traditional ab initio methodologies on the OsO4 molecule

The performance of different conventional ab initio methodologies and density functional procedures is compared through its application to the theoretical calculation of the bond distance and harmonic vibrational frequencies of the OsO₄ molecule. The problem of the basis set is first considered, with up to nine different basis sets being tested in calculations using the hybrid Becke3LYP density functional, and the most appropriate basis set is used in the comparison of Hartree–Fock, post‐Hartree–Fock, and density functional methods. The post‐Hartree–Fock methods analyzed are MP2, CISD, and CCSD(T), and the density functionals tested are SVWN, BLYP, BPW91, and Becke3LYP. The results show that for this particular system density functional methods perform better than do HF‐based methods with the exception of CCSD(T), which gives the best overall results. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 544–551, 2000     

DFT and ab initio direct dynamics studies on the hydrogen abstraction reactions of chlorine atoms with CH(4-n)F(n) (n = 1-3)

A direct dynamics study is carried out for the hydrogen abstraction reactions Cl + CH(4-n)F(n) (n = 1-3) in the temperature range of 200-1,000 K. The minimum energy paths (MEPs) of these reactions are calculated at the BH&H-LYP/6-311G(d,p) level, and the energies along the MEPs are further refined at the QCISD(T)/6-311+G(2df,2p) and QCISD(T)/6-311+G(d,p) (single-point) level. The rate constants obtained by using the improved canonical variational transition state theory incorporating small-curvature tunneling correction (ICVT/SCT) are in good agreement with the available experimental results. It is shown that the vibrational adiabatic potential energy curves for these reactions have two barriers, a situation similar to the analogous reactions CH(3)X+Cl (X=Cl, Br). The theoretical results show that for the title reactions the variational effect should not be neglected over the whole considered temperature range, while the small-curvature tunneling effect is only important in the lower temperature range. The effects of fluorine substitution on the rate of this kind of reactions are also examined.     

Rationalization of the behavior of M2(CH3CS2)4I (M = Ni,Pt) chains at room temperature from periodic density functional theory and ab initio cluster calculations

The electrical conductivities and plausible charge‐ordering states in the room temperature (r.t.) phase for MMX chains [Ni₂(dta)₄I]_∞ and [Pt₂(dta)₄I]_∞ (dta = CH₃CS) have been analyzed with periodic density functional theory (DFT) and correlated ab initio calculations combined with the effective Hamiltonian theory. Periodic DFT calculations show a more delocalized nature of the ground state in [Pt₂(dta)₄I]_∞ compared to [Ni₂(dta)₄I]_∞, which features a rather large energy gap between the occupied and empty bands, and charge polarized dimer units. A larger electrical conductivity for the Pt chain can be expected, especially because the Fermi level lies within a band with contributions from Pt and I orbitals. Electronic structure parameters extracted from ab initio cluster calculations show that the large difference between the observed conductivities at 300 K for Ni and Pt compounds, of 3 orders of magnitude, cannot be explained from the parameters extracted from an embedded M₂(dta)₄I₂ dimer fragment alone. When tetramer fragments are considered, we observe that the interdimer transfer integral (t) between neighboring M₂ units connected by an iodine atom at correlated level is comparable in both chains. On the other hand, the energy to transfer an electron from a dimer to the neighboring one (Coulomb repulsion U) is three times larger in the Ni compound with respect to the Pt chain, in line with the poor conductivity of the former. The electronic structure of the M₄(dta)₈I₃ fragment points to an alternate charge‐polarization state for Ni and an average valence state for Pt when the r.t. X‐ray structure is considered. © 2012 Wiley Periodicals, Inc.     

MP2, DFT and ab initio calculations on thioxanthone

Density functional theory (DFT), HF and MP2 calculations have been carried out to investigate thioxanthone molecule using the standard 6-31+G(d,p) basis set. The results of MP2 calculations show a butterfly structure for thioxanthone. The calculated results show that the predicted geometry can well reproduce the structural parameters. The predicted vibrational frequencies were assigned and compared with experimental IR spectra. A good harmony between theory and experiment is found. The theoretical electronic absorption spectra have been calculated using CIS method. ¹³C and ¹H NMR of the title compound have been calculated by means of B3LYP density functional method with 6-31+G(d,p) basis set. The comparison of the experimental and the theoretical results indicate that density functional B3LYP method is able to provide satisfactory results for predicting NMR properties.     

High‐level ab initio study of the N+(3P)+SH2 reactions in the gas phase: Role of spin‐forbidden pathways

The singlet and triplet potential energy surfaces involved in N⁺+SH₂ reactions have been explored using high‐level ab initio techniques. The geometries of the stationary points were optimized at the QCISD/6‐311G(df,p) level. The final energies were obtained in CCSD(T)/6‐311+G(3df,2p) single‐point calculations. The results obtained show that, although the N⁺(¹D)+SH₂ entrance channel is higher in energy than the N⁺(³P)+SH₂ one, most of the [H₂, S, N]⁺ singlet state cations are lower in energy than the corresponding triplets, due to their different bonding characteristics. Both singlet and triplet potential energy surfaces are quite close each other, and crossover between them can occur. The minimum energy crossing points were located by means of CASSCF(6,5) calculations. The spin‐orbit couplings show that the transition probability from the triplet to the singlet potential energy surface is significantly large. One of the most important consequences is that some of the products of the reaction, such as SH⁺, can be formed in typical spin‐forbidden processes. Since all the relevant structures along these pathways are much lower in energy than the reactants, this mechanism should be accessible even at low impact energies and therefore could be important in processes taking place in interstellar media. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Dissociation of sulfur oxoacids by two water molecules studied using ab initio and density functional theory calculations

Using ab initio [SCS‐MP2 and CCSD(T)] and density functional theory (M062X) calculations, we have studied the geometries and energies of sulfur oxoacids H₂S_mO₆ (m = 2–4) and their monohydrated and dihydrated clusters. When including the results from previously reported disulfuric acid (H₂S₂O₇) cases, the gas phase acidity is ordered as H₂S₂O₆ < H₂S₃O₆ < H₂S₂O₇ < H₂S₄O₆. The intramolecular H‐bonding, which may indicate the degree of structural flexibility in this molecular series, is an important factor for the order of the gas phase acidity. All these sulfur oxoacids show dissociated (or deprotonated) geometries with only two water molecules, although the energies of the dissociated conformers are ranked differently. All of the dissociated conformers form a unique H‐bonding network structure in which the protonated first water (H₃O⁺) is triply H‐bonded to each oxygen atom of two SO₃ moieties as well as the second water, which in turn is H‐bonded to a SO₃ moiety. H₂S₃O₆ has the best molecular flexibility for adopting such an H‐bonding network structure, and thereby all the low‐lying conformers of H₂S₃O₆(H₂O)₂ are dissociated. In contrast, the least flexible H₂S₂O₆ forms such a structure with a high strain, and dissociation of H₂S₂O₆(H₂O)₂ is found from the third lowest conformer. Although the gas phase acidity of H₂S₄O₆ is the highest in this series, the lowest dissociated conformer and the lowest undissociated conformer of H₂S₄O₆(H₂O)₂ are very close in energy. This is because forming the H‐bonding network structure is somewhat difficult due to the large distance between the two SO₃ moieties.     

On the kinetic mechanism of the hydrogen abstraction reactions of the hydroxyl radical with CH3CF2Cl and CH3CFCl2: A dual level direct dynamics study

By means of the dual‐level direct dynamics method, the mechanisms of the reactions, CH₃CF₂Cl + OH → products (R1) and CH₃CFCl₂ + OH → products (R2), are studied over a wide temperature range 200–2000 K. The optimized geometries and frequencies of the stationary points are calculated at the MP2/6‐311G(d,p) level, and then the energy profiles of the reactions are refined with the interpolated single‐point energy method at the G3(MP2) level. The canonical variational transition‐state theory with the small‐curvature tunneling (SCT) correction method is used to calculate the rate constants. For the title reactions, three reaction channels are identified and the H‐abstraction channel is the major pathway. The results indicate that F substitution has a significant (reductive) effect on hydrochlorofluorocarbon reactivity. Also, for all H‐abstraction reaction channels the variational effect is small and the SCT effect is only important in the lower temperature range on the rate constants calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

DNA base trimers: empirical and quantum chemical ab initio calculations versus experiment in vacuo

A complete scan of the potential‐energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller–Plesset second‐order perturbation theory (RI‐MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T‐shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol^?1. The different behaviors of nonmethylated and methylated trimers is also discussed. The high‐level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple‐zeta basis set predicts 11 uracil trimer interaction energies with a root‐mean‐square error of 2.3 kcal mol^?1 relative to highly correlated ab initio theoretical calculations.     

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

Thioformaldehyde S-methylide and thioacetone S-methylide: an ab initio MO study of structure and cycloaddition reactivity

The mechanisms of cycloaddition of thioformaldehyde S-methylide and thioacetone S-methylide, as models for an alkyl-substituted ylide, to thioformaldehyde and thioacetone, as well as to ethene as a model for a C=C double bond have been studied by ab initio calculations. Restricted and unrestricted B3LYP/6-31G* calculations were performed for the geometries of ground states, transition structures, and intermediates. Although basis sets with more polarization functions were tested, the 6-31G* basis set was applied throughout. Single-point CASPT2 calculations are reported for analysis of the unsubstituted system. The stabilities of structures with high biradical character seem to be overestimated by DFT methods in comparison to CASPT2. The general trends of the results are independent of the level of theory. Thioformaldehyde adds to thioformaldehyde S-methylide without activation energy, and the activation energies for two-step biradical pathways to 1,3-dithiolane are low. C,S biradicals are more stable than C,C biradicals. The two-step cycloaddition is not competitive with the concerted cycloaddition. Methyl substitution in the 1,3-dipole and the dipolarophile does not change the mechanistic relationships. TSs for the concerted formation of the regioisomeric cycloadducts of thioacetone Smethylide and thioacetone were located. Concerted addition remains the preferred reaction. The reactivity of the C=S double bond is high relative to that of the C=C double bond.     

HNNH3, a new possible isomer of N2H4: An ab initio study

Ab initio studies applying the 3-21G, 6-31G, and 6-31G** basis sets and also including the MP2 correction were carried out on H₂NNH₂, HNNH₃, and the transition state of the reaction H₂NNH₂(DOUBLE BOND)HNNH₃. First, the geometries of molecules were optimized using the theoretical methods mentioned in the restricted Hartree–Fock (RHF) scheme. The energies of the molecules corresponding to RHF/6-31G** geometries were subsequently calculated including electron-correlation effects at the level of the second-order Møller–Plesset (MP2) perturbation theory. The vibrational frequencies, net charges, and dipole moments were obtained from the theoretical calculations. The results of our calculations indicate unambiguously that H₂NNH₂ is thermodynamically more stable than is HNNH₃. On the other hand, an isolated HNNH₃ molecule once created would be stable since barriers for its unimolecular isomerization and decomposition are relatively high. But HNNH₃ is unlikely to be isolated in measurable amounts because of bimolecular tautomerization. Nevertheless, HNNH₃ can be considered as an intermediate in chemical processes involving N₂H₄. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 64 : 447–452, 1997     

Energy barriers of proton transfer reactions between amino acid side chain analogs and water from ab initio calculations

Proton transfer reactions were studied in all titratable pairs of amino acid side chains where, under physiologically reasonable conditions, one amino acid may function as a donor and the other one as an acceptor. Energy barriers for shifting the proton from donor to acceptor atom were calculated by electronic structure methods at the MP2/6-31++G(d,p) level, and the well-known double-well potentials were characterized. The energy difference between both minima can be expressed by a parabola using as argument the donor-acceptor distance R(DA). In this work, the fit parameters of the quadratic expression are determined for each donor-acceptor pair. Moreover, it was found previously that the energy barriers of the reactions can be expressed by an analytical expression depending on the distance between donor and acceptor and the energy difference between donor and acceptor bound states. The validity of this approach is supported by the extensive new data set. This new parameterization of proton transfer barriers between titratable amino acid side chains allows us to very efficiently estimate proton transfer probabilities in molecular modelling studies or during classical molecular dynamics simulation of biomolecular systems.     

Molecular structure and vibrational spectra of phenobaraitone by density functional theory and ab initio hartree-Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra at RHF/6-31++G** and B3LYP/6-31++G** levels for phenobarbitone (C₁₂H₁₂N₂O₃) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. Theoretical vibrational spectra of the title compound were interpreted by means of potential energy distributions (PEDs) using MOLVIB program. A detailed interpretation of the infrared spectra of the title compound is reported. On the basis of the agreement between the calculated and observed results, the assignments of fundamental vibrational modes of phenobarbitone were examined and some assignments were proposed. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title compound have been constructed.