Density functional theory and post-Hartree-Fock studies on molecular structure and harmonic vibrational spectrum of formaldehyde

Density functional theory (DFT) with the Becke's three-parameter exchange correlation functional and the functional of Lee, Yang and Parr, gradient-corrected functionals of Perdew, and Perdew and Wang [the DFT(B3LYP), DFT(B3P86) and DFT(B3PW91) methods, respectively], and several levels of conventional ab initio post-Hartree-Fock theory (second- and fourth-order perturbation theory M?ller-Plesset MP2 and MP4(SDTQ), coupled cluster with the single and double excitations (CCSD), and CCSD with perturbative triple excitation [CCSD(T)], configuration interaction with the single and double excitations [CISD], and quadratic configuration interaction method [QCISD(T)], using several basis sets [ranging from a simple 6-31G(d,p) basis set to a 6-311+ +G(3df, 2pd) one], were applied to study of the molecular structure (geometrical parameters, rotational constants, dipole moment) and harmonized infrared (IR) spectrum of formaldehyde (CH₂O). High-level ab initio methods CCSD(T) and QCISD(T) with the 6-311+ +G(3df, 2pd) predict correctly molecular parameters, vibrational harmonic wavenumbers and the shifts of the harmonic IR spectrum of ¹²CH₂ ¹⁶O upon isotopic substitution. Received: 30 January 1997 / Accepted: 7 May 1997     

CH2＝CHCl与O(3P)反应的理论研究

用量子化学密度泛函理论和QCISD (Quadratic configuration interaction calculation)方法, 对O(³P)与CH₂CHCl的反应进行了理论研究. 在UB3LYP/6-311＋＋G(d,p), UB3LYP/6-31＋＋G(3df, 3pd)计算水平上, 优化了反应物、产物、中间体和过渡态的几何构型, 并在UQCISD(T)/6-311＋＋G(2df,2pd)水平上计算了单点能量. 为了确证过渡态的真实性, 在UB3LYP/6-311＋＋G(3df,3pd)水平上进行了内禀坐标(IRC)计算和频率分析, 并确定了反应机理. 研究结果表明, 反应主要产物为CH₂CHO和Cl.     

Hydride abstraction of methylamine with Cu(S) in the gas phase: A density functional theory study

The gas-phase hydride abstraction of methylamine with Cu⁺(¹S) is theoretically investigated by using density functional theory. Geometries for all the stationary points involved are fully optimized at both the B3LYP/6-311++G(d,p) and B3LYP/6-311++G(3df,2p) levels and the reaction is analyzed in terms of the topology of potential energy surface. Approach of Cu⁺ towards methylamine could form either “classical” N or “nonclassical” η¹-methyl-H attached complex with the former being the global minimum. Both complexes are found to be key intermediates for the hydride abstraction, which could transfer into each other via two parallel routes, i.e., concerted metal movement and stepwise C-H activation-rearrangement. A charge-transfer process is detected for the “nonclassical” complex converting to a precursor species (CuH-NH₂CH₂⁺), which accounts for the final products by a nonreactive dissociation.     

Effects of different GIAO and CSGT models and basis sets on 2-aryl-1,3,4-oxadiazole derivatives

The direct molecular structure implementations of the gage-including atomic orbital (GIAO), individual gages for atoms in molecules (IGAIM) and continuous set of gage transformations (CSGT) methods for calculating nuclear magnetic shielding tensors at both the Hartree-Fock (HF) and density functional (B3LYP) levels of theory with 6-31G(d), 6-311G(d), 6-31++G(d,p), 6-311++G(d,p), and 6-311++G(df,pd) basis sets are presented. Dependence on the ¹H and ¹³C NMR chemical shifts on the choice of method and basis set have been investigated. Also, these chemical shifts of 2-aryl-1,3,4-oxadiazoles 5a–g have been performed related to dihedral angles (C4–C3–C2–O) of two conformers. The optimized molecular geometries and ¹H and ¹³C chemical shift values of 2-aryl-1,3,4-oxadiazoles 5a–g in the ground state have been obtained. The linear correlation coefficients of ¹³C NMR chemical shifts for these molecules were given. The new nuclear magnetic shielding tensors of tetramethylsilane (TMS) were calculated. The data of 2-aryl-1,3,4-oxadiazole derivatives display significant molecular structure and NMR analysis. Also, these provide the basis for future design of efficient materials having the 1,3,4-oxadiazole core.     

DFT Investigation of Nuclear 14N Quadrupole Coupling Constants for Four HMX Polymorphs

Density functional theory investigation of ¹⁴N nuclear quadrupole coupling constants (NQCC) for HMX polymorphs was performed. During the calculation of NQCC and asymmetry parameter (η) of β‐HMX, single molecule model is found to be worse than cluster model. The calculated results are more sensitive to the proper model than to the basis set. The calculation for β‐HMX using cluster model at B3LYP/6‐311++G(d,p) level gives better agreement with experiment. This approach was subsequently applied to α, γ, and δ‐HMXs. The difference of simulated ¹⁴N NQR frequency was found. The NQR spectrum is useful for the study of explosive or propellant on their detection, phase transition, and aging process.     

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.     

污泥热解中HCN与CaO的反应机理:密度泛函理论研究

采用密度泛函理论对污泥热解中CaO与HCN在低温段的反应进行了研究。在B3LYP/6-311++(3df,2p)水平上计算得到了反应路径上各驻点的几何构型与频率,并在此构型上使用CCSD(T)/cc-pVQZ进行单点能计算。结果表明,两个HCN分子吸附于CaO后,质子发生转移时出现反应路径中最大能垒(310.33 kJ/mol)。使用经典过渡态理论拟合了反应中各步骤的阿累尼乌斯公式,计算了三种典型温度下各步骤的反应速率,发现质子转移为该反应的决速步骤,且温度越高CaO对HCN的作用效果越好。     

Density functional theory studies on the dissociation energies of metallic salts: relationship between lattice and dissociation energies

The formation and physicochemical properties of polymer electrolytes strongly depend on the lattice energy of metal salts. An indirect but efficient way to estimate the lattice energy through the relationship between the heterolytic bond dissociation and lattice energies is proposed in this work. The heterolytic bond dissociation energies for alkali metal compounds were calculated theoretically using the Density Functional Theory (DFT) of B3LYP level with 6‐311+G(d,p) and 6‐311+G(2df,p) basis sets. For transition metal compounds, the same method was employed except for using the effective core potential (ECP) of LANL2DZ and SDD on transition metals for 6‐311+G(d,p) and 6‐311+G(2df,p) calculations, respectively. The dissociation energies calculated by 6‐311+G(2df,p) basis set combined with SDD basis set were better correlated with the experimental values with average error of ca. ±1.0% than those by 6‐311+G* combined with the LANL2DZ basis set. The relationship between dissociation and lattice energies was found to be fairly linear (r>0.98). Thus, this method can be used to estimate the lattice energy of an unknown ionic compound with reasonably high accuracy. We also found that the dissociation energies of transition metal salts were relatively larger than those of alkaline metal salts for comparable ionic radii. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 827–834, 2001     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Ab initio and density functional theory study of the enthalpies of formation of F2SOx and FClSOx (x=1, 2)

Enthalpies of formation of F₂SO, F₂SO₂, FClSO and FClSO₂ molecules have been determined using ab initio molecular orbital theory and density functional theory (DFT) calculations. Different DFT approaches and levels of the Gaussian-3 and the complete basis set (CBS) ab initio model chemistries have been employed to calculate enthalpies of formation from both total atomization energies and isodesmic reaction schemes. The best values at 298 K for F₂SO, F₂SO₂, FClSO and FClSO₂ as derived from an average of G3, G3B3, CBS-Q and CBS-QB3 isodesmic energies are −140.6, −181.1, −92.6 and −132.3 kcal mol⁻¹, respectively. The results obtained suggest that the accumulated small component errors found in the DFT-based methods are significantly reduced at the ab initio levels employed. Structural properties, harmonic vibrational frequencies, mode assignations and infrared intensities derived from B3LYP and mPW1PW91 functional with the 6-311+G(3df) basis set are presented.     

Gas phase and solution state stability of complexes L…M, where M = Cu, Ni, or Zn and L = R−C(O)NHOH (R = H, NH2, CH3, CF3, or Phenyl)

The structure and gas-phase metal affinities (M = Cu²⁺, Ni²⁺, and Zn²⁺) of formohydroxamic acid derivatives R–C(O)NHOH (R = H, NH₂, CH₃, CF₃ and Phenyl) were studied using the B3LYP/6-311+G(d,p) method of DFT theory. In order to evaluate the conformational behavior of these systems in water, we carried out CPCM-SCRF optimization calculations at the B3LYP/6-311+G(d,p) levels of theory. The obtained optimized geometries and interaction affinities of the gas and solution phase were compared. The following order of stability was found for ionic complexes of the transition metals: Cu²⁺ > Ni(t)²⁺ > Zn²⁺. The same stability order would be expected according to the Irving–Williams order of stability constants. The high-spin complexes of the Ni²⁺ were more stable than the low-spin complexes. The solvent effect reduced the observed relative stability of individual metallic complexes of substituted hydroxamic acids.     

Density functional theory study of a mechanism of formation of a spiro heterocycle with the Ge spiroatom from dimethylgermylidene germylene (Me2Ge=Ge:) and ethylene

The B3LYP/6-311++G** study of the mechanism of the cycloaddition of singlet 2,2-dimethyl-1,2-digermavinylidene (Me₂Ge=Ge:) to ethylene was performed. [2+2] Cycloaddition of the reactants first produced four-membered cyclic germylene, then the interaction of unoccupied 4p orbital of the Ge atom with the π orbital of another ethylene molecule yielded intermediate with the π → p donor-acceptor bond. Isomerization of this intermediate via low-laying transition state resulted in spirocyclic compound with the sp³ hybridized Ge spiroatom.     

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^?.     

14N nuclear quadrupole coupling in solid cinnolin-4-ones. A study by nuclear quadrupole resonance,x-ray crystallography and ab initio calculations

The ¹⁴N nuclear quadrupole coupling constants (NQCCs) for 1H-cinnolin-4-one and its 1- and 2-methyl derivatives have been obtained by nuclear quadrupole double resonance. The X-ray crystal structures of the 1H and 1Me compounds have been solved, and show lengths and angles consistent with related heterocycles; attempts to obtain an X-ray structure for the 2Me compound have so far been unsuccessful owing to disorder. Computed equilibrium structures by ab initio SCF methods for the Me derivatives and also the experimental 1H structures yielded electric field gradients, which were converted into NQCCs in reasonable agreement with experiment. The 2Me compound has unusual NQCC, both experimentally and in the theoretical study.     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

A theoretical study of the structure and vibrations of 2,4,6-trinitrotolune

Theoretical calculations of the structure, internal rotations and vibrations of 2,4,6-trinitrotolune, TNT, in the gas phase were performed at the B3LYP/6-31G* and B3LYP/6-311+G** levels of theory. Two genuine energy minimum structures were found. In both structures the 4-nitro group is planar to the phenyl ring, while the 2,6-nitro groups are slightly out of plane with the phenyl ring due to steric interaction with the methyl group. The two structures are related by internal rotations of the methyl and 2, or 6-nitro group. The lowest energy route for interconversion between them is a concerted motion of the methyl group and 2 or 6 nitro group in a ‘cog wheel’ type of mechanism. The geometry of the low energy structure A is closest to that observed in the crystal structures of TNT, where all three nitro groups are out of plane with the phenyl ring. FTIR and Raman spectra of solid TNT and ¹³C, ¹⁵N enriched TNT are presented and assigned with the help of the B3LYP/6-311+G** calculations on A. The lower level B3LYP/6-31G* calculation fails to predict the correct vibrational coupling between the nitro and phenyl groups. The B3LYP/6-311+G** calculation gives a good prediction of the nitro vibrations and the isotopic shifts observed for TNT isotopomers.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

A theoretical study of MgH2 ambient and high-pressure phases using NQCC parameters

Quadrupolar parameters of nuclei can be used as a tool to understand the electronic structure of the compounds. Magnesium hydride (MgH₂) is a potential hydrogen storage material due to its outstanding hydrogen capacity, however, its high thermodynamic stability is unfavorable for dehydrogenation processes. Understanding the bonding nature of Mg and H is essential for improving its dehydrogenation performance. In this work the charge density distribution in MgH₂ is studied. For this purpose, using calculated NQCCs of hydrogen atoms, the electronic structure of α-MgH₂ with several high pressure forms of MgH₂ were compared. The results show that in the high pressure phases (β, γ, and δ) some hydrogens have very small NQCC and therefore these hydrogens form weaker bond with Mg. In other words, easier condition for dehydrogenation in pressure-induced forms is expected. The electric field gradient (EFG) at the site of quadrupolar nuclei were calculated to obtain NQCC parameters using Gaussian 03 at B3LYP/6-31G level of theory. The selected level and basis set give the rather acceptable qualitative NQCCs of hydrogen atoms.     

Further evidence for the inseparability of the theoretical enthalpic terminal and penultimate unit effects in the radical copolymerization of styrene with acrylonitrile

The theoretical enthalpies of propagation reactions at 0 K without zero‐point vibrational energy corrections according to terminal and penultimate models of the radical copolymerization of styrene with acrylonitrile are reported from molecular orbital calculations at the following levels of theory and basis sets: HF/6‐31G(d); B3‐LYP/6‐31G(d); B3‐LYP/6‐311G(d,p) and B3‐LYP/6‐311+G(3df)//6‐311G(d,p). Both the enthalpic terminal and penultimate unit effects, determined according to the theoretical thermochemistry, depend on the level of theory and basis set used for the molecular orbital calculations. The best performing B3LYP/6‐311+G(3df)//B3LYP/6‐311G(d,p) procedure gives theoretical enthalpies for the addition of styrene and acrylonitrile to CH that differ from experimental values by 0.6 and 1.6 kcal mol^?1, respectively. An analysis of the results obtained here leads to the conclusion that at least for the styrene–acrylonitrile monomer system, that is, a monomer system known as one of the few systems that do not conform to terminal model composition and microstructure equations, the enthalpic terminal unit effects seem to depend on the penultimate units of the growing radical. This finding, together with the outcome from our previous work on the dependence of the penultimate effects on the terminal units in a growing macroradical, indicates the inseparability of the enthalpic terminal (implicit) and explicit penultimate unit effects on the radical copolymerization. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 1778–1787, 2003