BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

The basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Comprehensive study of the interaction between hydrogen halides and methanol derivatives

The R? CH₂? HO…H? X (R = SCl, Cl, SH, NO₂, OMe, CHO, CN, C₂H₅, CH₃, H; X = F, Cl, Br) complexes are considered here as the interest sample for the consideration of different measures of H‐bond strength. The intermolecular interaction energies are predicted by using MP2/6‐31++G(d,p) and B3LYP/6‐31++G(d,p) methods with basis set superposition error and zero‐point energy corrections. The results showed that intermolecular hydrogen bonds for complexes with HF are stronger than such interactions in complexes with HCl and HBr. Quantum theory of “Atoms in Molecules” and natural bond orbitals method were applied to analyzed H‐bond interactions. The gas phase thermodynamic properties of complexes were predicted using quantum mechanical computations. The obtained results showed a strong influence of the R and X substituents on the thermodynamic properties of complexes. Numerous correlations between topological, geometrical, thermodynamic properties and energetic parameters were also found. © 2011 Wiley Periodicals, Inc.     

Effects of different basis sets and donor‐acceptor groups on linear and second‐order nonlinear optical properties and molecular frontier orbital energies

The calculation of molecular hyperpolarizability, molecular frontier orbital energies of some donor‐acceptor oxadiazoles ( 5a – f , 8a – f , and 9a – f ) have been investigated using ab initio methods and different basis sets. Ab initio optimizations were performed at the Hartree–Fock (HF) and density functional (Beckee‐3–Lee–Yang–Parr; B3LYP) levels of theory with 6‐31G basis set. The polarizability (<α>), anisotropy of polarizability (Δα), and ground‐state dipole moment (μ), first hyperpolarizability (β), and molecular frontier orbital (HOMO, highest occupied molecular orbital and LUMO, lowest unoccupied molecular orbital) energies of 5a – f , 8a – f , and 9a – f have been calculated at the HF and B3LYP methods with 6‐31G, 6‐31G(d), 6‐31+G(d), 6‐31++G(d,p), 6‐311G, 6‐311G(d), 6‐311+G(d), and 6‐311++G(d,p) basis sets. Also, the molecular hardness (η) and electronegativity (χ) parameters have been obtained using molecular frontier orbital energies. The <α>, Δα, μ, β, HOMO, LUMO energies, η and χ parameters have been investigated as dependence on the choice of method and basis set. The variation graphics of <α>, Δα, μ, β, η, and χ parameters using HF and B3LYP methods with different basis sets are presented. We have examined the frontier molecular orbital pictures of 5a – f , 8a – f , and 9a – f using B3LYP/6‐31++G(d,p) level. The 5a – f , 8a – f , and 9a – f display significant linear, second‐order molecular nonlinearity, and molecular parameters and provide the basis for future design of efficient nonlinear optical materials having the 1,3,4‐oxadiazole core. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

DFT-GIAO calculations of 19F NMR chemical shifts for perfluoro compounds

The (19)F NMR shieldings for 53 kinds of perfluoro compounds were calculated by the B3LYP-GIAO method using the 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311++G(d,p), 6-311G(2d,2p), 6-311++G(2d,2p), 6-311++G(2df,2p), 6-311++G(3d,2p), and 6-311++G(3df,2p) basis sets. The diffuse functions markedly reduce the difference between the calculated and experimental chemical shifts. The calculations using the 6-31++G(d,p) basis set give the chemical shifts within 10 ppm deviations from experimental values except for the fluorine nuclei attached to an oxygen atom, a four- and a six-coordinated sulfur atom, and FC(CF(3))(2) attached to a sulfur atom.     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical Studies on Intermolecular Hydrogen-bond Interactions between Hexamethylenetetramine and Nitric Acid

The structures of the complexes generated by hexamethylenetetramine and nitric acid have been fully optimized by B3LYP method at the 6-311++G** and aug-cc-pVTZ levels. The intermolecular hydrogen-bonding interactions have been calculated by the B3LYP/6-311++G**, B3LYP/aug-cc-pVTZ, MP2(full)/6-311++G** and CCSD(T)/6-311++G** methods, respectively. The NBO (nature bond orbital), AIM (atom in molecule), temperature effect and solvation effect have been analyzed to reveal the origin of the interactions. The results indicate that the stable hydrogen-bonded complexes could be generated by hexamethylenetetramine and nitric acid. The interactions follow the order of (a)>>(e)>(b)>(c)>(d)>(f)>(g). The C-N bonds which are adjacent to the methylene involving the hydrogen bonds tend to break in the chemical reaction. Due to the exothermic process, low temperature is conducive to the formation of the composition, which tallies with the experimental result.     

BSSE‐free description of the formamide dimers

The different configurations (linear, zig‐zag, and cyclic) of formamide dimers have been studied at the level of both Hartree–Fock (HF) and second order Møller–Plesset perturbation theory (MP2). The widely used a posteriori Boys–Bernardi “counterpoise” (CP) correction scheme has been compared with our a priori methods utilizing the “chemical Hamiltonian approach” (CHA). The appropriate interaction energies have been calculated in six different basis sets (6‐31G, 6‐31G**, DZV, DZP, TZV, and cc‐pVDZ). © 2001 John Wiley & Sons, Inc. Int J Quant Chem, 2001     

Theoretical analysis of vibrational spectra and scaling‐factor of 2‐aryl‐1,3,4‐oxadiazole derivatives

In this study, the direct molecular structure implementations for calculating vibrational spectra and scaling factors, and infrared intensities 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), and 6‐311++G(d,p) basis sets are presented. Also, vibrational frequencies have been investigated as dependence on the choice of method and basis set. The parameters of molecular geometry and vibrational frequencies values of 2‐aryl‐1,3,4‐oxadiazoles 5a–g in the ground state have been calculated. Theoretical determination of vibrational frequencies is quite useful both in understanding the relationship between the molecular structures and scaling factor. The data of 2‐aryl‐1,3,4‐oxadiazoles 5a–g display significant electronic properties provide the basis for future design of efficient materials having the oxadiazole core and theoretical IR studies. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Theoretical calculations on the mechanisms of the gas‐phase elimination kinetics of 2‐chloro‐1‐phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane,and their corresponding chloroalkanes: The effect of the phenyl ring

The kinetics and mechanisms of the dehydrochlorination of 2‐chloro‐1‐ phenylethane, 3‐chloro‐1‐phenylpropane, 4‐chloro‐1‐phenylbutane, 5‐chloro‐1‐phenylpentane, and their corresponding chloroalkanes were examined by means of electronic structure calculation using density functional theory methods B3LYP/6–31G(d,p), B3LYP/6–31++G(d,p), MPW1PW91/6–31G(d,p), MPW1PW91/6–31++G(d,p), PBEPBE/6–31G(d,p), and PBEPBE/6–31++G(d,p). The potential energy surface was investigated for the minimum energy path. Calculated enthalpies and energies of activation are in good agreement with experimental values using the MPW1PW91 and B3LYP methods. The transition state of these reactions is a four‐centered cyclic structure. The reported experimental results proposing neighboring group participation by the phenyl group was not supported by theoretical calculations. The rate‐determining process in these reactions is the breaking of Cl? C bond. The reactions are described as concerted moderately polar and nonsynchronous. © 2011 Wiley Peiodicals, Inc. Int J Chem Kinet 43: 292–302, 2011     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Theoretical study of the N-H…O red-shifted and blue-shifted hydrogen bonds

Theoretical calculations are performed to study the nature of the hydrogen bonds in complexes HCHO…HNO, HCOOH…HNO, HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F. The geometric structures and vibrational frequencies of these six complexes at the MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) levels are calculated by standard and counterpoise-corrected methods, respectively. The results indicate that in complexes HCHO…HNO and HCOOH…HNO the N-H bond is strongly contracted and N-H…O blue-shifted hydrogen bonds are observed. While in complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F, the N-H bond is elongated and N-H…O red-shifted hydrogen bonds are found. From the natural bond orbital analysis it can be seen that the X-H bond length in the X-H…Y hydrogen bond is controlled by a balance of four main factors in the opposite directions: hyperconjugation, electron density redistribution, rehybridization and structural reorganization. Among them hyperconjugation has the effect of elongating the X-H bond, and the other three factors belong to the bond shortening effects. In complexes HCHO…HNO and HCOOH…HNO, the shortening effects dominate which lead to the blue shift of the N-H stretching frequencies. In complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F where elongating effects are dominant, the N-H…O hydrogen bonds are red-shifted.     

A combined density functional theory and coupled cluster method investigation of the structural properties and stabilities of radical CH(2)CP and its isomers

The doublet potential energy surface of radical system [C(2), H(2), P] is investigated at the UB3LYP/6-311++G(d,p) and UCCSD(T)/6-311++G(2df,2p) (single-point) levels. Eight chainlike and three-membered ring structures are located as energy minima connected by 10 interconversion transition states. At the final UCCSD(T)/6-311++G(2df,2p)//UB3LYP/6-311++G(d,p) level with zero-point vibrational energy correction, species CH(2)CP is found to be thermodynamically the most stable isomer followed by HCCPH, H-cCPC-H, cPCC-H(H), H-cCCP-H, cis-CC(H)PH, trans-CC(H)PH, and CCPH(2) at 11.01, 12.57, 40.07, 43.63, 50.25, 56.82, and 65.36 kcal/mol, respectively. The computed results indicate that the chainlike isomers CH(2)CP and HCCPH and cyclic radical H-cCPC-H possess considerable kinetic stability at extra low pressures and temperatures. Interestingly, radical CCPH(2), whose energy is the highest in all predicted CH(2)CP isomers, can be also regarded as a kinetically stable species with the smallest isomerization barrier of 22.26 kcal/mol at extra low pressures and temperatures. Therefore, considering higher kinetic stability, in addition to the microwave spectroscopy characterized isomer CH(2)CP in previous experiments, the species HCCPH, H-cCPC-H, and CCPH(2) should be considered as excellent candidates for possible experimental observation. Furthermore, the structural nature of stable radical isomers is discussed based on bonding characteristics, single electron spin distribution, and comparison with their analogues.     

复合物HNO···H2O2内N-H···O蓝移氢键的理论研究

A theoretical study on the blue-shifted H-bond N-H…O and red-shifted H-bond O-H…O in the complexHNO…H_2O_2 was conducted by employment of both standard and counterpoise-corrected methods to calculate thegeometric structures and vibrational frequencies at the MP2/6-31G(d),MP2/6-31 G(d,p),MP2/6-311 q G(d,p),B3LYP/6-31G(d),B3LYP/6-31 G(d,p) and B3LYP/6-311 G(d,p) levels.In the H-bond N-H…O,the calcu-lated blue shift of N-H stretching frequency is in the vicinity of 120 cm~(-1) and this is indeed the largest theoreticalestimate of a blue shift in the X-H…Y H-bond ever reported in the literature.From the natural bond orbital analy-sis,the red-shifted H-bond O-H…O can be explained on the basis of the dominant role of the hyperconjugation.For the blue-shifted H-bond N-H…O,the hyperconjugation was inhibited due to the existence of significant elec-tron density redistribution effect,and the large blue shift of the N-H stretching frequency was prominently due tothe rehybridization of sp~n N-H hybrid orbital.     

Assessment of density functionals with long‐range and/or empirical dispersion corrections for conformational energy calculations of peptides

Density functionals with long‐range and/or empirical dispersion corrections, including LC‐ωPBE, B97‐D, ωB97X‐D, M06‐2X, B2PLYP‐D, and mPW2PLYP‐D functionals, are assessed for their ability to describe the conformational preferences of Ac‐Ala‐NHMe (the alanine dipeptide) and Ac‐Pro‐NHMe (the proline dipeptide) in the gas phase and in water, which have been used as prototypes for amino acid residues of peptides. For both dipeptides, the mean absolute deviation (MAD) is estimated to be 0.22–0.40 kcal/mol in conformational energy and 2.0–3.2° in torsion angles ? and ψ using these functionals with the 6‐311++G(d,p) basis set against the reference values calculated at the MP2/aug‐cc‐pVTZ//MP2/aug‐cc‐pVDZ level of theory in the gas phase. The overall performance is obtained in the order B2PLYP‐D ≈ mPW2PLYP‐D > ωB97X‐D ≈ M06‐2X > MP2 > LC‐ωPBE > B3LYP with the 6–311++G(d,p) basis set. The SMD model at the M06‐2X/6‐31+G(d) level of theory well reproduced experimental hydration free energies of the model compounds for backbone and side chains of peptides with MADs of 0.47 and 4.3 kcal/mol for 20 neutral and 5 charged molecules, respectively. The B2PLYP‐D/6‐311++G(d,p)//SMD M06‐2X/6‐31+G(d) level of theory provides the populations of backbone and/or prolyl peptide bond for the alanine and proline dipeptides in water that are consistent with the observed values. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Theory Studies on the Intermolecular Interactions of Urea Nitrate with RDX

Six fully optimized geometries of urea nitrate cation and RDX complexes have been obtained with DFT-B3LYP and MP2 methods at the 6-311++G** level. The intermolecular interaction energies have been calculated with basis set superposition error (BSSE) and zero point energy (ZPE) correction. The nature of intermolecular interaction has been revealed by the analysis of AIM and NBO. The results indicate that the greatest binding energy of urea nitrate with RDX is –82.47kJ/mol. The O–H…O and N–H…O hydrogen bonds are important intermolecular interactions of urea nitrate cation with RDX, and the origin of hydrogen bonds is the oxygen atom offering its lone-pair electrons to the σ(O-H)* or σ(O-H)* antibonding orbital. The intermolecular interactions strengthen the N–NO2 bond, leading to the reduced sensitivity of urea nitrate and RDX mixture explosive.     

Comparison of basis set superposition error corrected perturbation theories for calculating intermolecular interaction energies

Recently, two different but conceptually similar basis set superposition error (BSSE) free second‐order perturbation theoretical schemes were developed by us that are being based on the chemical Hamiltonian approach (CHA). Using these CHA‐MP2 and CHA‐PT2 methods, a comparison is made between the a posteriori and a priori BSSE correction schemes at the correlated level. Sample calculations are presented for four hydrogen bonded complexes (HFH₃N, HFH₂O, H₂SHF, and H₂OHCl) in nine different basis sets (from 6–31G to TZV**++). The results show that the BSSE content is very significant in the interaction energy if electron correlation is accounted for, so removing the BSSE is very important. The differences of the two perturbational theories discussed are connected solely with the different one electron orbital sets used for building up the unperturbed single determinant wave function. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 274–283, 1999     

BeH2与HX(X=F,Cl, Br,I)形成的二氢键复合物的结构特征与本质

对BeH2与HX(X=F, Cl, Br, I)形成的二氢键复合物的结构特征及本质进行了探讨. 在MP2/6-311++G(3d,3p)水平优化、频率验证, 得到复合物的分子结构, 用分子间距离及电子密度拓扑理论确认BeH2与卤化氢已形成了二氢键型复合物. 在MP2/6-311++G(3d, 3p)水平下进行基函数重叠误差(BSSE)校正后的结合能在-14.468 kJ·mol-1到-5.464 kJ·mol-1之间.用对称匹配微扰理论(SAPT)对复合物的结合能进行分解, 结果表明, BeH2…HX二氢键复合物中静电能对总吸引能的贡献都是最主要的, 但交换排斥能、诱导能、色散能对总结合能的贡献也很重要. 从BeH2…HF到BeH2…HI, 诱导能对总吸引能的贡献从37.8%逐渐减小到24.0%. 而色散能对总吸引能的贡献从BeH2…HF体系中的16.0%逐渐增加到BeH2…HI体系中的33.8%.     

A theoretical study on decomposition and rearrangement reaction mechanism of trichloroacetyl chloride (CCl3COCl)

In the present study, the possible decomposition and rearrangement reaction profile of trichloroacetyl chloride have been studied using UMP2/6‐311++G (2d, 2p) level of ab initio and UB3LYP/6‐311++G (2d, 2p) level of density functional theory methods. The harmonic vibrational frequencies were calculated at the same level of theory used for the characterization of stationary points and zero‐point vibrational energy corrections. The potential energy barrier and activation energy between each step of the reaction have been calculated for the seven possible reaction pathways (Ia–c, IIa–b, IIIa–b). The trichloroacetyl chloride is an asymmetric ketone where the two α bonds of acetyl chloride, the C? C and C? Cl bonds are strong with dissociation energy of 72 kcal/mol. The phosgene (COCl₂), dichloroketene (CCl₂CO), carbon dichloride (CCl₂), carbon tetrachloride (CCl₄), and carbon monoxide (CO) are the major dominant products on the decomposition of the trichloroacetyl chloride. These resultant products are more hazards than the parent trichloroacetyl chloride molecules. The positive value of the reaction energy indicate that the overall reaction profile is found to be endothermic at the UMP2 and UBLYP/6‐311++G(2d, 2p) levels of theory, respectively, at UMP2/6‐311G** optimized geometry. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011