Benchmarking of density functionals for the kinetics and thermodynamics of the hydrolysis of glycosidic bonds catalyzed by glycosidases

In recent years, there has been an increased interest in understanding the enzymatic mechanism of glycosidases resorting mostly to DFT and DFT/MM calculations. However, the performance of density functionals (DFs) is well known to be system and property dependent. Trends drawn from general studies, despite important to evaluate the quality of the DFs and to pave the way for the development of new DFs, may be misleading when applied to a single specific system/property. To overcome this issue, we carried out a benchmarking study of 40 DFs applied to the geometry optimization and to the electronic barrier height (E _Barrier) and electronic energy of reaction (E _R) of prototypical glycosidase‐catalyzed reactions. Additionally, we report calculations with SCC‐DFTB and four semiempirical MO methods applied to the same problem. We have used a universal molecular model for retaining glycosidases, comprising only a 22‐atoms system that mimics the active site and substrate. High accuracy reference geometries and energies were calculated at the CCSD(T)/CBS//MP2/aug‐cc‐pVTZ level of theory. Most DFs reproduce the reference geometries extremely well, with mean unsigned errors (MUE) smaller than 0.05 Å for bond lengths and 3° for bond angles. Among the DFs, wB97X‐D, CAM‐B3LYP, B3P86, and PBE1PBE have the best performance in geometry optimizations (MUE = 0.02 Å). Conversely, semiempirical MO and SCC‐DFTB methods yielded less accurate geometries (MUE between 0.09 and 0.17 Å). The inclusion of D3 correction has a small, but still relevant, influence in the geometry predicted by some DFs. Regarding E _Barrier, 11 DFs (MPW1B95, CAM‐B3LYP, M06 ‐ 2X, PBE1PBE, wB97X ‐ D, B1B95, BMK, MN12 – SX, M05, M06, and M11) presented errors below 1 kcal.mol^?1, in relation to the reference energy. Most of these functionals belong to the family of hybrid functionals (H‐GGA, HH‐GGA, and HM‐GGA), which shows a positive influence of HF exchange in the determination of E _Barrier. The inclusion of D3 correction has not affected significantly the E _Barrier and E _R. The use of geometries at the accurate but expensive MP2/aug‐cc‐pVTZ level of theory has a small, albeit not insignificant, influence in the E _Barrier when compared with energies calculated with geometries determined with the DFs (usually a few tenths of kcal.mol^?1, with exceptions). In general, semiempirical MO methods and DFTB are associated with larger errors in the determination of E _Barrier, with unsigned errors from 6.9 to 24.7 kcal.mol^?1.     

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     

One‐bond 29Si‐1H spin‐spin coupling constants in the series of halosilanes: benchmark SOPPA and DFT calculations,relativistic effects,and vibrational corrections

A number of most representative second order polarization propagator approach (SOPPA) based wavefunction methods, SOPPA, SOPPA(CC2) and SOPPA(CCSD), and density functional theory (DFT) based methods, B3LYP, PBE0, KT2, and KT3, have been benchmarked in the calculation of the one‐bond ²⁹Si‐¹H spin‐spin coupling constants in the series of halosilanes SiH_nX_4?n (X = F, Cl, Br, I), both at the non‐relativistic and full four‐parameter Dirac's relativistic levels taking into account vibrational corrections. At the non‐relativistic level, the wavefunction methods showed much better results as compared with those of DFT. At the DFT level, out of four tested functionals, the Perdew, Burke, and Ernzerhof's PBE0 showed best performance. Taking into account, relativistic effects and vibrational corrections noticeably improves wavefunction methods results, but generally worsens DFT results. Copyright © 2013 John Wiley & Sons, Ltd.     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

Comparison of aromatic NH···π, OH···π, and CH···π interactions of alanine using MP2, CCSD, and DFT methods

A comparison of the performance of various density functional methods including long‐range corrected and dispersion corrected methods [MPW1PW91, B3LYP, B3PW91, B97‐D, B1B95, MPWB1K, M06‐2X, SVWN5, ωB97XD, long‐range correction (LC)‐ωPBE, and CAM‐B3LYP using 6‐31+G(d,p) basis set] in the study of CH···π, OH···π, and NH···π interactions were done using weak complexes of neutral (A) and cationic (A⁺) forms of alanine with benzene by taking the Møller–Plesset (MP2)/6‐31+G(d,p) results as the reference. Further, the binding energies of the neutral alanine–benzene complexes were assessed at coupled cluster (CCSD)/6‐31G(d,p) method. Analysis of the molecular geometries and interaction energies at density functional theory (DFT), MP2, CCSD methods and CCSD(T) single point level reveal that MP2 is the best overall performer for noncovalent interactions giving accuracy close to CCSD method. MPWB1K fared better in interaction energy calculations than other DFT methods. In the case of M06‐2X, SVWN5, and the dispersion corrected B97‐D, the interaction energies are significantly overrated for neutral systems compared to other methods. However, for cationic systems, B97‐D yields structures and interaction energies similar to MP2 and MPWB1K methods. Among the long‐range corrected methods, LC‐ωPBE and CAM‐B3LYP methods show close agreement with MP2 values while ωB97XD energies are notably higher than MP2 values. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Ab initio reaction path for cisplatin interaction with L‐cysteine and L‐methionine

Ab initio calculations are used to track the reaction pathway of interaction between cisplatin and the sulfur‐containing amino acids cysteine (Cys) and methionine (Met). Structures of all reactive species as well as thermodynamic and kinetic properties were calculated and discussed based on the role played by the level of theory. Twenty‐three different levels of theory were examined including HF, DFT, and perturbation theory at MP2 and MP4(SDQ) orders. The rate constant for a second‐order associative ligand exchange mechanism (k₂) was calculated by means of transition state theory. This quantity is quite sensitive to small fluctuation of activation free energy, therefore is a good benchmark to assess the performance of different methods of calculations. The k₂ values predicted by DFT methods were in best agreement with experiment, found equal to (10²k₂ in M^?1 s^?1) 3.42 for Met (PBE1PBE) and 1.90 for Cys (B3P86). The experimental values are 3.6 and 2.2 for Met and Cys, respectively. The solvent effect plays a primary role to the kinetic properties, accounting for ～30% of the activation Gibbs free energy. The outcomes from the present study promptly show the adequacy of distinct theoretical approaches to describe the reactivity of cisplatin, thus might be useful for further studies involving other Pt(II) complexes. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Accurate Structure and Bonding Description of the Transition Metal‐Disulfur Monoxide Complexes [(PMe3)2M(S2O)] (M = Ni,Pd, Pt): Grimme Dispersion Corrected DFT Study

Geometry and bonding energy analysis of M–S₂O bonds in the metal‐disulfur monoxide complexes [(PMe₃)₂M(S₂O)] of nickel, palladium, and platinum were investigated at DFT, DFT‐D3, and DFT‐D3(BJ) methods using three different functionals (BP86, PBE, and TPSS). The TPSS/DFT‐D3(BJ) yields better geometry, while the BP86 geometry is least accurate for studied complexes. The geometry of platinum complex optimized at TPSS/DFT‐D3(BJ) level is in excellent agreement with the available experimental values. The M–S bonds are shorter than the M–S(O) bonds. The Mayer bond orders suggest the presence of M–S and M–S(O) single bonds. Both the M–S and M–S(O) bond lengths vary with the density functionals as TPSS‐D3(BJ) < TPSS < PBE < BP86. The Hirshfeld charge distribution indicates that the overall charge flows from metal fragment to [S₂O]. The Ni–S₂O bond has greater degree of covalent character than the ionic. The contribution of dispersion interactions is large in computing accurate bond dissociation energies between the interacting fragments. The BDEs are largest for the functional TPSS and smallest for the functional BP86. The DFT‐D3 dispersion corrections to the BDEs between the metal fragments [(PMe₃)₂M] and ligand fragment [(S₂O)] for the TPSS functional are in the range 7.1–7.3 kcal · mol^–1, which are smaller than the corresponding DFT‐D3(BJ) dispersion corrections (9.4–10.6 kcal · mol^–1).     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Thermochemistry and Kinetics of the Thermal Decomposition of 1‐Chlorohexane

A theoretical kinetic study of the thermal decomposition of 1‐chlorohexane in gas phase between 600 and 1000 K was performed. Transition‐state theory and unimolecular reaction rate theory were combined with molecular information provided by quantum chemical calculations. Particularly, the B3LYP, BMK, M05–2X, and M06–2X formulations of the density functional theory (DFT) and the high‐level ab initio methods G3B3 and G4 were employed. The possible reaction channels for the thermal decomposition of 1‐chlorohexane were investigated, and the reaction takes place through the elimination of HCl with the formation of 1‐hexene. The derived high‐pressure limit rate coefficients are k _∞(600–1000 K) = (8 ± 5) × 10¹³ exp[‐((56.7 ± 0.4) kcal mol⁻¹/RT )] s⁻¹. The pressure effect over the reaction was analyzed from the calculation of the low‐pressure limit rate coefficients and the falloff curves. In addition, the standard enthalpies of formation at 298 K of −46.9 ± 1.5 kcal mol⁻¹ for 1‐chlorohexane and 5.8 ± 1.5 kcal mol⁻¹ for C₆H₁₃ radical were derived from isodesmic and isogiric reactions at high levels of theory.     

Difluorophosphoryl Nitrene F2P(O)N: Matrix Isolation and Unexpected Rearrangement to F2PNO

Triplet difluorophosphoryl nitrene F₂P(O)N (X³A′′) was generated on ArF excimer laser irradiation (λ=193 nm) of F₂P(O)N₃ in solid argon matrix at 16 K, and characterized by its matrix IR, UV/Vis, and EPR spectra, in combination with DFT and CBS‐QB3 calculations. On visible light irradiation (λ>420 nm) at 16 K F₂P(O)N reacts with molecular nitrogen and some of the azide is regenerated. UV irradiation (λ=255 nm) of F₂P(O)N (X³A′′) induced a Curtius‐type rearrangement, but instead of a 1,3‐fluorine shift, nitrogen migration to give F₂PON is proposed to be the first step of the photoisomerization of F₂P(O)N into F₂PNO (difluoronitrosophosphine). Formation of novel F₂PNO was confirmed with ¹⁵N‐ and ¹⁸O‐enriched isotopomers by IR spectroscopy and DFT calculations. Theoretical calculations predict a rather long P? N bond of 1.922 Å [B3LYP/6‐311+G(3df)] and low bond‐dissociation energy of 76.3 kJ mol^?1 (CBS‐QB3) for F₂PNO.     

Dihalomethanes as Bent Bifunctional XB/XB‐Donating Building Blocks for Construction of Metal‐involving Halogen Bonded Hexagons

The dihalomethanes CH₂X₂ (X=Cl, Br, I) were co‐crystallized with the isocyanide complexes trans‐[MX^M₂(CNC₆H₄‐4‐X^C)₂] (M=Pd, Pt; X^M=Br, I; X^C=F, Cl, Br) to give an extended series comprising 15 X‐ray structures of isostructural adducts featuring 1D metal‐involving hexagon‐like arrays. In these structures, CH₂X₂ behave as bent bifunctional XB/XB‐donating building blocks, whereas trans‐[MX^M₂(CNC₆H₄‐4‐X^C)₂] act as a linear XB/XB acceptors. Results of DFT calculations indicate that all XCH₂–X???X^M–M contacts are typical noncovalent interactions with estimated strengths in the range of 1.3–3.2 kcal mol^?1. A CCDC search reveals that hexagon‐like arrays are rather common but previously overlooked structural motives for adducts of trans‐bis(halide) complexes and halomethanes.     

MP2, DFT‐D,and PCM study of the HMB–TCNE complex: Thermodynamics,electric properties,and solvent effects

Geometry, thermodynamic, and electric properties of the π‐EDA complex between hexamethylbenzene (HMB) and tetracyanoethylene (TCNE) are investigated at the MP2/6‐31G* and, partly, DFT‐D/6‐31G* levels. Solvent effects on the properties are evaluated using the PCM model. Fully optimized HMB–TCNE geometry in gas phase is a stacking complex with an interplanar distance 2.87 × 10^?10 m and the corresponding BSSE corrected interaction energy is ?51.3 kJ mol^?1. As expected, the interplanar distance is much shorter in comparison with HF and DFT results. However the crystal structures of both (HMB)₂–TCNE and HMB–TCNE complexes have interplanar distances somewhat larger (3.18 and 3.28 × 10^?10 m, respectively) than our MP2 gas phase value. Our estimate of the distance in CCl₄ on the basis of PCM solvent effect study is also larger (3.06–3.16 × 10^?10 m). The calculated enthalpy, entropy, Gibbs energy, and equilibrium constant of HMB–TCNE complex formation in gas phase are: ΔH⁰ = ?61.59 kJ mol^?1, ΔS = ?143 J mol^?1 K^?1, ΔG⁰ = ?18.97 kJ mol^?1, and K = 2,100 dm³ mol^?1. Experimental data, however, measured in CCl₄ are significantly lower: ΔH⁰ = ?34 kJ mol^?1, ΔS = ?70.4 J mol^?1 K^?1, ΔG⁰ = ?13.01 kJ mol^?1, and K = 190 dm³ mol^?1. The differences are caused by solvation effects which stabilize more the isolated components than the complex. The total solvent destabilization of Gibbs energy of the complex relatively to that of components is equal to 5.9 kJ mol^?1 which is very close to our PCM value 6.5 kJ mol^?1. MP2/6‐31G* dipole moment and polarizabilities are in reasonable agreement with experiment (3.56 D versus 2.8 D for dipole moment). The difference here is due to solvent effect which enlarges interplanar distance and thus decreases dipole moment value. The MP2/6‐31G* study supplemented by DFT‐D parameterization for enthalpy calculation, and by the PCM approach to include solvent effect seems to be proper tools to elucidate the properties of π‐EDA complexes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

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.     

对氨基苯甲酸热力学性质的量热和热分析研究

在80～400 K温区,用高精度全自动绝热量热仪测定了对氨基苯甲酸摩尔热容,得到摩尔热容随温度的变化的关系式为：     

Thermochemical properties for n‐propyl,iso‐propyl,and tert‐butyl nitroalkanes,alkyl nitrites,and their carbon‐centered radicals

Density functional theory (DFT) based calculations are performed on a series of alkyl nitrites and nitroalkanes representing large‐scale primary, secondary, and tertiary nitro compounds and their radicals resulting from the loss of their skeletal hydrogen atoms. Geometries, vibration frequencies, and thermochemical properties [S°(T) and C°_p(T) (10 K ? T ? 5000 K)] are calculated at the B3LYP/6‐31G(d,p) DFT level. Δ_fH°₂₉₈ values are from B3LYP/6‐31G(d,p), B3LYP/6‐31+G(2d,2p), and the composite CBS‐QB3 levels. Potential energy barriers for the internal rotations have been computed at the B3LYP/6‐31G(d,p) level of theory, and the lower barrier contributions are incorporated into entropy and heat capacity data. The standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Recommended values derived from the most stable conformers of respective nitro‐ and nitrite isomers include ?30.57 and ?28.44 kcal mol^?1 for n‐propane‐, ?33.89 and ?32.32 kcal mol^?1 for iso‐propane‐, ?42.78 and ?41.36 kcal mol^?1 for tert‐butane‐nitro compounds and nitrites, respectively. Entropy and heat capacity values are also reported for the lower homologues: nitromethane, nitroethane, and corresponding nitrites. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 181–199, 2010     

Structures and stabilization energies of methyl anions with main group substituents from the first five periods

The effects of substituents (X) on the structures and stabilities of CH₂X^? anions for groups comprised of fourth- and fifth-period main group elements (X = K, CaH, GaH₂, GeH₃, AsH₂, SeH, Br, Rb, SrH, InH₂, SnH₃, SbH₂, TeH, and I) have been investigated by ab initio pseudopotential calculations. Full geometry optimizations have been carried out on the CH₂X^? anions and the corresponding neutral parent molecules, CH₃X, at HF/DZP + and MP2/DZP + levels. Results for substituents from the second (X = Li? F) and third (X = Na? Cl) periods provide comparisons of substituent effects of the main group elements of the first four rows of the periodic table on methyl anions. Frequency calculations characterize the nature of stationary points and show pyramidal CH₂X^? anion structures to be the most stable unless π acceptor interactions (e.g., with BH₂, AlH₂, GaH₂, and InH₂ favor planar geometries. The CH₂X^? stabilization energies [at QCISD(T)/DZP + /MP2/DZP + + ZPE level for X = K? I and QCISD(T)/6?31 + G*/MP2/6?31 + G* + ZPE level] for X = Li? Cl) also show strong π-stabilizing effects for the same substituents. With the exception of CH₃ and NH₂, all substituents stabilize methyl anions, although the σ stabilization by OH and F is small. The SiH₃? PH₂? SH? Cl, GeH₃? AsH₂? SeH? Br, and SnH₃? SbH₂? TeH? I sets of substituents give stabilization energies between 19 and 30 kcal/mol. The stability of methyl anions substituted by the halogens and the chalcogens (X = OH, SH, SeH, and TeH) increases down a group in accord with the increasing substituent polarizability, while for π acceptors (BH₂, AlH₂, GaH₂, and InH₂) the stability decreases down a group in line with their π-accepting ability. © 1994 by John Wiley & Sons, Inc.     

4- (4,6-二甲基嘧啶-2-基)-3-硫代脲酸甲酯的合成晶体结构及量子化学研究

4-（4,6-Dimethylpyrimidin-2-yl）-3-thio-allophanic acid methyl ester was synthesized with mixing 2-amino-4,6- dimethylpyrimidine, potassium thiocyanate and methyl chloroformate in ethyl acetate. Single crystals suitable for X-ray diffraction measurement were obtained by recrystallization from dimethylformamide at room temperature. The crystal belongs to monoclinic symmetry with space group C2/m, and crystal parameters of a= 1.7537（5） nm, b= 0.6759（2） nm, c=1.1148（3） nm, β=118.557（4）°, V=1.1605（6） nm^3, Z=4, De= 1.375 g/cm^3,μ=0.271 mm^-1, F（000）=504, and 1519 [1〉2σ（I）] observable independent reflections were used for the determination and refmement of the crystal structures with final R1 of 0.0372 and wR2 of 0.0992. The theoretical investigation of the title compound was carried out with DRT-B3LYP/6-311G, HF/6-311G and MP2/6-311G methods, and the atomic net charges and the population were discussed.     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.