三硝基甲烷键离解能和生成焓的理论计算

采用密度泛函(DFT)四种交换/相关函数(B3LYP、B3P86、B3PW91和PBE0)结合不同的基函数,求得了三硝基甲烷C-NO2键的离解能(BDE),并且通过合理选择参考物硝基甲烷,设计等键等电子对反应,计算了气相三硝基甲烷分子的生成焓(HOF).与实验数据进行比较,PBE0/6-31g*计算出的BDE值最好,误差为-2.1 kcal mol-1;PBE0密度泛函结合带极化函数的6-31g基组得到的HOF值与实验值吻合的最好(误差在0.1 kcal mol-1以内).     

Evaluation of the performance of non-local and hybrid density functional theory methods for pi-radical hyperfine structures

A total of 22 substituted benzene radicals are investigated by means of density functional theory (DFT). Two different functional schemes, the B3LYP hybrid DFT-HF functional and the ‘pure’ gradient corrected DFT functional PWP86, previously employed successfully in calculations of radical hyperfine parameters, are employed. The results are compared with experimental data. The same strategy is used for all calculations: initial geometry optimization using a double-ζ basis set, followed by single point hyperfine calculations using the larger IGLO-III, 6-311G(2df, p) and aug-cc-pCVTZ basis sets. It is found that of the two functionals, the PWP86 functional gives the more homogeneous data, with nearly constant underestimations in HFCCs by 5–15%. With the B3LYP approach the smaller couplings, found at the meta positions, generally are overestimated considerably. The best overall performance is observed at the PWP86/6-311G(2df, p) level.     

Pdn(n=1-9)团簇的结构及磁性的第一性原理计算

采用密度泛函理论中的广义梯度近似泛函BPW91和三参数杂化密度泛函B3LYP对Pdn(n=1-9)团簇的结构、稳定性和磁性进行了详细的计算. 两种泛函得到了相同的稳定结构, 除n=3、4外, 两种方法得到的基态结构是完全一致的, 但在n=3、4时用三参数杂化密度泛函B3LYP得到的基态结构与文献[8,9,23]的相同. 两种方法得到的平均配位数和平均键长有相似的变化规律, 总体上随团簇尺寸的增大而增大, n=2-5时, 增幅较大, n=5-9时, 增幅较小. 两种方法得到团簇能量的二阶差分、分裂能在n=4时均有较大的值, 说明相对应的团簇具有较高的稳定性、较低的化学活性. 两种方法得到团簇的平均每原子磁矩随团簇尺寸的增大有逐渐减小的趋势, 个别团簇有振荡. 结果表明两种泛函都可以描述团簇结构、稳定性和磁性的演变规律, 但B3LYP泛函可以更加精确地描述Pdn团簇的结构演化.     

Pdn(n=1~9)团簇的结构及磁性的第一性原理计算简

采用密度泛函理论中的广义梯度近似泛函BPW91和三参数杂化密度泛函B3LYP对Pdn(n=1~9)团簇的结构、稳定性和磁性进行了详细的计算.两种泛函得到了相同的稳定结构,除n=3、4外,两种方法得到的基态结构是完全一致的,但在n=3、4时用三参数杂化密度泛函B3LYP得到的基态结构与文献[8,9,23]的相同.两种方法得到的平均配位数和平均键长有相似的变化规律,总体上随团簇尺寸的增大而增大,n=2~5时,增幅较大,n=5~9时,增幅较小.两种方法得到团簇能量的二阶差分、分裂能在n=4时均有较大的值,说明相对应的团簇具有较高的稳定性、较低的化学活性.两种方法得到团簇的平均每原子磁矩随团簇尺寸的增大有逐渐减小的趋势,个别团簇有振荡.结果表明两种泛函都可以描述团簇结构、稳定性和磁性的演变规律,但B3LYP泛函可以更加精确地描述Pdn团簇的结构演化.     

TiN分子基态(X2Σ)结构和势能函数

应用群论及原子分子反应静力学方法推导了TiN分子基态(X²Σ)的离解极限. 采用不同的密度泛函方法,包括BP86, B3P86, B3LYP, B3PW91, 分别选用不同的基组对TiN分子基态进行结构优化计算.通过比较得出使用BP86方法, 对N原子使用D95V++(d,P)基组和Ti原子使用6-311++G^**基组时,计算得到的平衡几何结构、分子离解能和谐振频率与实验值符合得最好. 并采用最小二乘法拟合改进的Murrell-Sorbie函数得到了相应电子态的完整势能函数. 计算得到的光谱常数与实验光谱数据符合得很好. 关键词： BP86 TiN分子基态 势能函数 光谱常数     

On the convergence of zero-point vibrational corrections to nuclear shieldings and shielding anisotropies towards the complete basis set limit in water

ABSTRACT

The method and basis set dependence of zero-point vibrational corrections (ZPVCs) to nuclear magnetic resonance shielding constants and anisotropies has been investigated using water as a test system. A systematic comparison has been made using the Hartree–Fock, second-order Møller–Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) and Kohn–Sham density functional theory with the B3LYP exchange-correlation functional methods in combination with the second-order vibrational perturbation theory (VPT2) approach for the vibrational corrections. As basis sets, the correlation consistent basis sets cc-pVXZ, aug-cc-pVXZ, cc-pCVXZ and aug-cc-pCVXZ with X = D, T, Q, 5, 6 and the polarisation consistent basis sets aug-pc-n and aug-pcS-n with n = 1, 2, 3, 4 were employed. Our results show that basis set convergence of the vibrational corrections is not monotonic and that very large basis sets are needed before a reasonable extrapolation to the basis set limit can be performed. Furthermore, our results suggest that coupled cluster methods and a decent basis set are required before the error of the electronic structure approach is lower than the inherent error of the VPT2 approximation.     

Sextic centrifugal distortion constants: interplay of density functional and basis set for accurate yet feasible computations

Quantum-chemical calculations assist the analysis of laboratory spectra, and often provide the only means to determine spectroscopic data that cannot be accessed experimentally. Accurate predictions of vibrational and rotational spectroscopic parameters are required for applications in the field of high-resolution molecular spectroscopy. While the accuracy issue of the quantum-chemical calculation of vibrational properties and of equilibrium structures has been addressed in the literature, the same is not true for centrifugal distortion constants that however play an essential role for the interpretation of remote sensing data. In this work, the performance of several model chemistries, rooted mainly in density functional theory, in computing sextic centrifugal distortion constants is assessed employing a benchmark set of molecules of both atmospheric and astrochemical relevance. The Jensen’s (aug-)pcs-n basis sets, different flavours of Dunning’s triple-ζ basis sets and the SNSD basis set, are employed in conjunction with different functionals, and their predictions are benchmarked against experimental and theoretical data at the coupled cluster level of theory. This study also demonstrates the reliability of the calculation of sextic centrifugal distortion constants within the Gaussian16 rev. B.01 program package. Reliable predictions of the sextic centrifugal distortion constants for the gauche- and trans-conformers of ethyl-mercaptan are also presented.     

Gaussian basis set of triple zeta valence quality for the atoms from K to Kr: Application in DFT and CCSD(T) calculations of molecular properties

A contracted basis set of triple zeta (TZ) valence quality for the atoms from K to Kr was constructed from fully-optimized Gaussian basis sets generated in this work. Gaussian polarization functions (d, f, and g symmetries), which were optimized at the second-order Mφller–Plesset level, were added to the TZ set. This extends earlier work on segmented contracted TZ basis set for atoms H-Ar. This set along with the BP86 non-hybrid and B3LYP hybrid functionals were used to calculate geometric parameters, dissociation energy, harmonic vibrational frequency, and electric dipole moment of a sample of molecules and, then, comparison with results obtained with other basis sets and with experimental data reported in the literature is done. CCSD(T) atomic excitation energies and bond lengths, dissociation energies, and harmonic vibrational frequencies of some diatomics were also evaluated. Using density functional theory and gauge-including atomic orbitals, ⁵⁷Fe and ⁷⁷Se nuclear magnetic resonance chemical shifts in Fe(C₅H₅)₂, H₂Se, (CH₃)SeH, CSe₂, SeCO, H₂CSe, and SeF₆ were calculated. Comparison with theoretical and experimental values previously published in the literature was done. It is verified that in general these results give good agreement with experimental and benchmark values.     

含能材料中键离解能的密度泛函理论计算

通过计算一个包括臭氧,硝基甲烷以及1,3,5 -三硝基-1,3,5-三氮杂环已烷(RDX)在内的典型系统的键离解能,对由四种交换/相关函数(BLYP,B3LYP,B3PW91和B3P86)加上不同的基函数组合而成的多种密度泛函方法的准确性进行了比较研究.结果表明:B3P86/6-31G**是计算该系统C-NO2,O-O和N-NO2键离解能的最可靠的方法.     

On the performance of density functional schemes for computing the static dipole polarizability of 4d transition-metal monohalides

This article presents the static dipole polarizabilities of 4d transition-metal monohalides in the framework of density functional theory. The behavior of a large number of density functionals, including various types of the local spin density approximation (LSDA), the generalized gradient approximation (GGA), meta GGA (MGGA), hybrid GGA (HGGA), and hybrid meta GGA (HMGGA), has been assessed in polarizability calculations. The highly accurate CCSD(T) methodology is used to benchmark a variety of functionals. Our results indicate that there is a large variation in performance of the various functionals of each type in the calculation of the dipole polarizability. It turned out that, in the LSDA class, SVWN5 gives better results than SVWN3 for the reference values. Of the GGA methods, BP86 and BPW91 give the best results for our data set, followed by G96LYP. The TPSS functional is the best method of the MGGA class, followed by the VSXC and BB95 functionals. Moreover, the B98 and mPW1PW91 functionals are the best methods of the HGGA lineage, predicting the most accurate static dipole polarizabilities. Finally, TPSS1KCIS and PBE1KCIS are the best methods of the HMGGA functionals for that purpose. In general, it was found that the HGGA and HMGGA functionals provide the best performance.     

Theoretical Study of the C-Cl Bond Dissociation Enthalpy and Electronic Structure of Substituted Chlorobenzene Compounds

利用密度泛函(DFT)三种交换/相关函数(B3LYP, B3PW91,B3P86)结合6-31G**和6-311G**基组,计算了13个取代氯苯化合物的键离解能. 结果表明B3P86/6-311G**方法是计算取代氯苯化合物键离解能的可信方法,研究发现C-Cl键的键离解能与所使用的基组和计算方法密切相关,取代基对C-Cl键的键离解能的影响不明显. 研究了目标化合物的前线轨道能级差,并对取代氯苯化合物的热稳定性做了评估.     

Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies(BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory(DFT)(B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set(CBS-Q) method in conjunction with the 6311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.     

密度泛函理论研究MN@H_2O (M=Ga,Ge, In,Sn, Sb; N=M,Al)团簇的特性

采用密度泛函理论的B3LYP, B3P86, B1B95, P3PW91和PBE1PBE方法结合SDD, LANL2DZ和CEP-121G基组计算了d~(10)组态二聚物MN(M=Ga, Ge, In, Sn和Sb; N=M和Al)的几何结构.采用B3P86/SDD进一步研究了MN@H_2O团簇的几何结构及吸附能.结果表明,水分子结合在二聚物M_2上时,对二聚物影响较大,对水分子自身影响较小.将M_2中Ga, Ge, In, Sn或Sb替换一个原子为Al时,水分子在GeAl和SnAl上的吸附能变化较大,而在GaAl, InAl和SbAl上吸附能变化较小.另外, H_2O吸附在Ga, Ge, In, Sn和Sb上时,与吸附在Al上时,吸附能的变化不大.     

On the possibility of interpretation of the spectra of tetrapyrrole compounds and their radical ions on the basis of quantum-chemical calculations by the density functional method

Using the time-dependent formalism of the density functional theory (time-dependent density functional theory (TDDFT)), the energies and intensities of the electronic transitions of radical anions of Mg porphyrin and Mg phthalocyanine are calculated quantum chemically. Based on these calculations, an interpretation of the electronic spectra of these compounds is proposed. It is shown that, for the neutral Mg porphin molecule, the results of the TDDFT calculations are consistent with the generally accepted notions and with the experimental data only if hybrid density functionals are applied. In this case, the best results are obtained for an exchange-correlation functional whose parameters are close to the standard parameters of the B3LYP functional. The spectra of the electronic transitions of radical anions of the considered tetrapyrrole molecules calculated by the TDDFT (B3LYP) method are consistent both with previous semiempirical calculations and with the experimental data.     

Ab initio study of superoxide anion—water clusters O2 − (H2O)n=1-5

The structures, energies, Spectra, and charge transfer to solvent (CTTS) energies of the hydrated superoxide anion clusters have been investigated using both density functional calculations (DFT) with Becke-3-parameters employing Lee-Yang-Parr functionals (B3LYP) and second-order M?ller Plesset perturbation (MP2) calculations employing highly diffuse basis sets. Given the good agreement of our predicted results with the experimental spectra and enthalpy, it would be interesting if the CTTS energies could be experimentally determined.     

Difficulties in applying pure Kohn-Sham density functional theory electronic structure methods to protein molecules

Self-consistency-based Kohn-Sham density functional theory (KS-DFT) electronic structure calculations with Gaussian basis sets are reported for a set of 17 protein-like molecules with geometries obtained from the Protein Data Bank. It is found that in many cases such calculations do not converge due to vanishing HOMO-LUMO gaps. A sequence of polyproline I helix molecules is also studied and it is found that self-consistency calculations using pure functionals fail to converge for helices longer than six proline units. Since the computed gap is strongly correlated to the fraction of Hartree-Fock exchange, test calculations using both pure and hybrid density functionals are reported. The tested methods include the pure functionals BLYP, PBE and LDA, as well as Hartree-Fock and the hybrid functionals BHandHLYP, B3LYP and PBE0. The effect of including solvent molecules in the calculations is studied, and it is found that the inclusion of explicit solvent molecules around the protein fragment in many cases gives a larger gap, but that convergence problems due to vanishing gaps still occur in calculations with pure functionals. In order to achieve converged results, some modeling of the charge distribution of solvent water molecules outside the electronic structure calculation is needed. Representing solvent water molecules by a simple point charge distribution is found to give non-vanishing HOMO-LUMO gaps for the tested protein-like systems also for pure functionals.     

Density functional theory calculations of nuclear quadrupole coupling constants with calibrated 14N quadrupole moments

Density functional calculations of the electric field gradient tensor at the nitrogen nucleus in 13 test molecules, containing 14 nitrogen sites, have been performed using the linear combination of Gaussian-type orbital Kohn-Sham density functional theory (LCGTO-KSDFT) approach. Local and gradient corrected functionals were used for all-electron calculations. All the molecular structures were optimized at their respective levels of theory with extended basis sets. Calibrated ¹⁴N nuclear quadrupole moments were obtained through a fitting procedure between calculated electric field gradients and experimental nuclear quadrupole coupling constants of the test set of molecules for each basis set and functional considered. With these calibrated ¹⁴N nuclear quadrupole moments, the nuclear quadrupole coupling constants of the following selected systems were determined: fluoromethylisonitrile, pyridine, pyrrole, imadazole, pyrazole, 1,8-bis(dimethyl-amino)naphthalene, cyclotetramethylenetetranitramine, cocaine and heroin.