碳氟键均裂解离能的密度泛函理论研究

运用6种密度泛函方法(B3LYP, B3P86, B3PW91, PBE1PBE, MPW1B95, MPW1K)对15个含氟有机化合物的碳氟键均裂解离能进行理论计算, 得到的理论值与实验值比较, 发现B3P86方法用于碳氟键均裂解离能的计算相对可靠. 使用验证后的理论方法对含氟杂环有机化合物和卤氟烃中的碳氟键均裂解离能进行了预测和分析, 并进一步讨论了α-取代基效应以及Hammett型取代基效应对碳氟键均裂解离能的影响.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.     

邻苯二甲酸单甲酯桥联双核铜(Ⅱ)配合物磁性的理论研究

采用对称性破损态方法结合密度泛函理论,选用典型的强反铁磁双核配合物作为研究对象,通过与实验数据相比较,探讨了不同密度泛函方法与基组对计算铜配合物[Cu2(mMP)4(H2O)2]·H2O交换耦合常数的准确度.结果表明,4种混合密度泛函DFT(B3LYP,B3P86,B3PW91和PBE0)的计算结果都能和实验所观察到的值-324cm-1符号一致,但B3PW91方法得到的结果和实验结果吻合程度最好,同时采用方法B3PW91方法计算所得的交换耦合常数Jab对基组的依赖性较大.研究表明,2个Cu(Ⅱ)离子之间的反铁磁相互作用主要源于单占据分子轨道SOMOs大的能量劈裂和桥联配体O-C-O轨道的重叠.     

间位取代卟啉及其锌卟啉的电子吸收光谱

取代的卟啉类衍生物在气敏传感器方面具有广泛的应用前景.本文采用了密度泛函理论(DFT)和含时密度泛函理论(TD-DFT)研究了四种不同取代基的卟啉衍生物(meso位四硝基苯基/四氨基苯基卟啉(NO2PP,NH2PP)及其相应的锌金属卟啉衍生物(NO2ZnPP,NH2ZnPP))的紫外和近紫外光谱特征.利用两种不同的交换相关泛函(广义梯度近似泛函(PBE)和杂化密度泛函(B3LYP))优化了上述四种物质的结构,并应用TD-DFT计算了相应的电子激发能量和振动强度.结果表明,取代卟啉的吸收光谱与大量的电子跃迁有关;与B3LYP泛函预测的光谱相比,PBE泛函所得B带以及Q带的波长位置与实验值更为接近.另外,计算所得硝基取代基卟啉的B带相对于氨基取代基卟啉的B带发生了红移,这与实验现象也保持一致.由于卟啉衍生物的三重激发态在电子转移中有很重要的应用,因此在PBE/6-31G(d)水平上计算了四种物质的最低三重激发态能量,分别为1.426、1.469、1.608和1.581eV.     

萤光素类似物激发态的TD-DFT计算精度评估

采用含时密度泛函理论(TD-DFT)对8种已知萤光素类似物的垂直激发能和发射能进行了系统的考察.选取10种交换-相关(XC)泛函对8种萤光素类似物的基态和第一单重激发态结构、吸收和发射光谱进行了计算,并将得到的结果与实验数据进行对照.结果表明,该系列物质吸收光谱和发射光谱的计算对XC泛函的选择非常敏感.B3LYP、mPW3PBE、B3PW91方法能够提供较好理论计算结果,吸收光谱的均方根误差(RMS)在0.40eV以内,标准差(SD)在0.27eV以内;发射光谱的RMS在0.24eV以内,SD在0.17eV以内.     

3-甲基-4-对甲基苯基-5-(2-吡啶)-1,2,4-三唑桥联的双核铜配合物的磁性理论研究

采用对称性破损态方法结合密度泛函理论,选用反铁磁双核配合物[Cu2(MMBPT)2Cl4(H2O)2.5](MMBPT=3-甲基-4-对甲基苯基-5-(2-吡啶)-1,2,4-三唑)作为研究对象,通过与实验数据相比较,探讨了不同密度泛函方法与基组对计算铜配合物交换耦合常数的准确度.结果表明,4种混合密度泛函DFT(B3LYP,B3P86,B3PW91和PBE0)的计算结果都能和实验所观察到的值-31cm-1符号一致,但只有B3PW91方法得到的结果和实验结果吻合程度最好,同时采用方法B3PW91方法计算所得的交换耦合常数Jab对基组的依赖性较大.研究表明,2个Cu(Ⅱ)离子之间弱的反铁磁相互作用主要源于单占据分子轨道SOMOs小的能量劈裂.     

核自旋偶合常数2JC-C-F的量子化学计算

用密度泛函B3LYP方法,在6-311G基组上对36个化合物进行构型优化、频率分析,计算了C-C-F的核自旋偶合常数^2JC-C-F。由此拟合出一个基于该方法和相同水平上计算^2JC-C-F的公式,选取四种化合物进行验证,计算值与文献值吻合。     

X-H (X = C, N, O, Si, P, S) 键离解能计算的密度泛函方法研究

使用了不同密度泛函方法计算X-H (X = C, N, O, Si, P, S) 键离解能，并分析不同密度泛函方法的计算精度。研究发现大多数密度泛函方法包括B3LYP, B3P86, B3PW91, G96LYP, PBE1PBE，和BH&HLYP都明显低估键离解能13-25 kJ/mol。该现象与是否使用无限基组无关，因为即使使用无限基组键离解能仍然被低估。因此密度泛函方法不适合用于键离解能的估算。其中B3P86方法的偏差最小。进一步分析表明，使用限制性开壳层计算并无任何优势，在大多数情况下非限制性开壳层计算实际上比限制性开壳层计算要好。最后，我们发现了密度泛函方法对键离解能的低估是系统的，因此建议利用校准后的UDFT/6-311++G(d, p)方法计算化学键离解能。     

脯氨酸的构象及性质

用X3LYP法在6-311++G(d, p)和6-311++G(3df, 3pd)基组水平上对脯氨酸15种构象的几何结构、相对能量、电子结构、红外光谱、偶极矩、极化率等性质进行了研究, 并与PBE1PBE/6-311++G(d, p)的结果和文献相比较, 从而得到: (1) 的脯氨酸的15种构象中能量最低的有4种, 不同构象中存在着强弱不同的5种氢键, 其中以N…H—O氢键最强, 并存在特殊的C—H…O=C氢键. 两种方法计算的几何结构数据相近, 均与实验值吻合; (2) 在构象相对能差计算方面, X3LYP具有明显的优势, 用中等基组就可以得到与高水平从头算法和大基组相同的结果, 而PBE1PBE法计算的相对能值则相差较大; (3) 脯氨酸不同构象中偶极矩最大和极化率最小的是最稳定的构象1和2, 两种方法计算的结果一致.     

结合神经网络方法和扩大训练基组构建新B3LYP泛函

神经网络方法成功地应用于修正密度泛函理论B3LYP方法中的三个参数(a0、ax和ac)以构建新B3LYP交换相关泛函.本文采用包含输入层、隐藏层和输出层的三层式神经网络结构.总电子数、多重度、偶极矩、动能、四极矩和零点能被选为物理描述符.296个能量数据被随机地分成两组,246个能量数据作为训练集以确定神经网络的最优结构和最优突触权重,50个能量数据作为测试集以测试神经网络的预测能力.修正后的三个参数觔0、觔x、觔c从输出层处得到,并用于计算体系的热化学性质如原子化能(AE)、电离势(IP)、质子亲合能(PA)、总原子能(TAE)和标准生成热(ΔfH苓).修正后的计算结果优于传统B3LYP/6-311+G(3df,2p)方法的计算结果.经过神经网络修正后,296个物种的总体均方根偏差从41.0 kJ.mol-1减少到14.2 kJ.mol-1.     

Calculations of ionization energies and electron affinities for atoms and molecules: A comparative study with different methods

In the present work, we examined the performance of 36 density functionals, including the newly developed doubly hybrid density functional XYG3 (Y. Zhang, X. Xu, and W. A. Goddard III, Proc. Natl. Acad. Sci, USA, 2009, 106, 4963), to calculate ionization energies (IEs) and electron affinities (EAs). We used the well-established G2-1 set as reference, which contains 14 atoms and 24 molecules for IE, along with 7 atoms and 18 molecules for EA. XYG3 leads to mean absolute deviations (MADs) of 0.057 and 0.080 eV for IEs and EAs, respectively, using the basis set of 6-311 + G (3df,2p). In comparison with some other functionals, MADs for IEs are 0.109 (B2PLYP), 0.119 (M06-2X), 0.159 (X3LYP), 0.161 (PBE), 0.162 (B3LYP), 0.165 (PBE0), 0.173 (TPSS), 0.200 (BLYP), and 0.215 eV (LC-BLYP). MADs for EAs are 0.090 (X3LYP), 0.090 (B2PLYP), 0.102 (PBE), 0.103 (M06-2X), 0.104 (TPSS), 0.105 (BLYP), 0.106 (B3LYP), 0.126 (LC-BLYP), and 0.128 eV (PBE0).     

Scaling factors for vibrational frequencies and zero-point vibrational energies of some recently developed exchange-correlation functionals

The scaling factors for the vibrational frequencies and zero-point vibrational energies evaluated at various combinations of recently developed exchange-correlation functionals and various basis sets are reported. The exchange-correlation functionals considered are B972, B98, HCTH, OLYP, O3LYP, G96LYP, PBE0 and VSXC functionals; the basis sets employed are 3-21G, 6-31G*, 6-31G**, 6-31+G, 6-311G*, 6-311G**, 6-311G(df,p), 6-311+G(df,p), cc-pVDZ and aug-cc-pVDZ. The experimental harmonic frequencies of 122 small molecules and the zero-point vibrational energies of 39 small molecules are used to determine the scaling factors through the least-square fitting procedure. It was found that the scaling factors do not depend significantly on the basis sets considered. The vibrational frequency scaling factors evaluated by using the B98 and PBE0 functionals are recommended on the basis of smallest root mean square error. The zero-point vibrational energy scaling factors evaluated from the B972 functional with Pople's double-zeta basis set and the HCTH functional with Pople's triple-zeta basis set are recommended on the basis of smallest root mean square error.     

Intermolecular potentials of the silane dimer calculated with Hartree-Fock theory, Møller-Plesset perturbation theory, and density functional theory

We have calculated the intermolecular interaction potentials of the silane dimer at the D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with 108 functionals chosen from the combinations of 9 exchange and 12 correlation functionals. Single-point coupled cluster [CCSD(T)] calculations have also been carried out to calibrate the correlation effect. The HF calculations yield unbound potentials largely because of the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater type orbitals fitted with Gaussian functions (STO-nG, n = 3 approximately 6), Pople's medium size basis sets [up to 6-311++G(3df,3pd)], to Dunning's correlation consistent basis sets (cc-pVXZ and aug-cc-pVXZ, X = D, T, Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(3d,3p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy ( approximately 0.05 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the expected C6 value from molecular polarizability calculations. We attribute the slow convergence partly to the inefficacy of using the MP2 calculations with Gaussian type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the expected potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energies calculated using the OPTXHCTH147, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals and the equilibrium bond lengths calculated using the MPWHCTH93, TPSSHCTH, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals are close to the MP2 results using the 6-311++G(3df,3pd) basis set. A correlation between the calculated DFT potentials and the exchange and correlation enhancement factors at the low-density region has been elucidated. The asymptotic behaviors of the DFT potentials are also analyzed.     

Dimers of boroglycine and methylamine boronic acid: a computational comparison of the relative importance of dative versus hydrogen bonding

Boronic acids are widely used in materials science, pharmacology, and the synthesis of biologically active compounds. In this Article, geometrical structures and relative energies of dimers of boroglycine, H2N-CH2-B(OH)2, and its constitutional isomer H3C-NH-B(OH)2, were computed using second-order M?ller-Plesset perturbation theory and density functional theory; Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the MP2 calculations, and the Pople 6-311++G(d,p) basis set was employed for a majority of the DFT calculations. Effects of an aqueous environment were incorporated into the results using PCM and COSMO-RS methodology. The lowest-energy conformer of the H2N-CH2-B(OH)2 dimer was a six-membered ring structure (chair conformation; Ci symmetry) with two intermolecular B:N dative-bonds; it was 14.0 kcal/mol lower in energy at the MP2/aug-cc-pVDZ computational level than a conformer with the classic eight-centered ring structure (Ci symmetry) in which the boroglycine monomers are linked by a pair of H-O...H bonds. Compared to the results of MP2 calculations with correlation-consistent basis sets, DFT calculations using the PBE1PBE and TPSS functionals with the 6-311++G(d,p) basis set were significantly better at predicting relative conformational energies of the H2N-CH2-B(OH)2 and H3C-NH-B(OH)2 dimers than corresponding calculations using the BLYP, B3LYP, OLYP, and O3LYP functionals, particularly with respect to dative-bonded structures.     

Reaction energy benchmarks of hydrocarbon combustion by Gaussian basis and plane wave basis approaches

The reaction energies of 275 elementary reactions from the hydrocarbon combustion model GRI-Mech 3.0 were evaluated by electronic structure calculations using both localized Gaussian basis and plane wave basis sets. In the Gaussian basis calculations, the d-polarization function on C, N, and O elements reduces the mean absolute deviation (MAD) from the experimental value by 53%, a significant improvement in computational accuracy. In the plane wave basis calculation using different exchange-correlation (XC) functionals, the MAD values were 0.316–0.426 eV when non-hybrid type XC functionals such as RPBE, PBE, PW91, revPBE, and PBEsol were used. On the other hand, hybrid functionals like B3LYP and HSE06 reduced the MAD values significantly down to 0.182 and 0.233 eV, respectively. The B3LYP results have 49% less MAD compared to the PBE results. These demonstrated the strong advantage of the hybrid functional for calculating gas-phase reaction energies. The present comprehensive benchmarks will be crucial for future microkinetics as well as machine learning studies on the catalytic reactions. © 2019 Wiley Periodicals, Inc.     

Ab initio procedure for aqueous-phase pKa calculation: the acidity of nitrous acid

We present an ab initio procedure for accurately calculating aqueous-phase pKa values and apply it to study the acidity of nitrous acid (HNO2, or HONO). The aqueous-phase pK(a) of nitrous acid was obtained from calculated gas-phase acidities and solvation free energies via a thermodynamic cycle and the solvation model chemistry of Barone et al. (J. Chem. Phys. 1997, 107, 3210). Solvation free energies were calculated at the HF/6-31G(d) level using the dielectric-polarizable continuum and the integral equation formalism-polarizable continuum solvent models (D-PCM and IEF-PCM, respectively), with the D-PCM model yielding the most accurate pKa values. For HF free energies of solvation, significant improvements in accuracy could be made by moving to the larger 6-311++G(3df,3pd) and aug-cc-pVQZ basis sets. Solvation free energies were also calculated using the density functional theory (DFT) methods B3LYP, TPSS, PBE0, B1B95, VSXC, B98 and O3LYP, with the most accurate methods being TPSS and VSXC, which provided average errors of less than 0.11 pKa units. Solvation free energies calculated with the different DFT methods were relatively insensitive to the basis set used. Our theoretical calculations are compared with experimental results obtained using stopped flow spectrophotometry. The pKa of nitrous acid was measured as 3.16 at 25 degrees C, and the enthalpy and entropy of nitrous acid dissociation were calculated from measurements as 6.7 kJ mol(-1) and -38.4 J mol(-1) K(-1), respectively, between 25 and 45 degrees C. The UV/visible absorption spectra of the nitrite ion and nitrous acid were also examined, and molar extinction coefficients were obtained for each.     

Assessing alkyl-, silyl-, and halo-substituent effects on the electron affinities of silyl radicals

Neutral anion energy differences for a large class of alpha-substituted silyl radicals have been computed to determine the effect of alkyl, silyl, and halo substituents on their electron affinities. In particular, we report theoretical predictions of the adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) for a series of methyl-, silyl-, and halo-substituted silyl radical compounds. This work utilizes the carefully calibrated DZP++ basis set, in conjunction with the pure BLYP and OLYP functionals, as well as with the hybrid B3LYP, BHLYP, PBE1PBE, MPW1K, and O3LYP functionals. Bromine has the largest effect in stabilizing the anions, and the BLYP/DZP++ AEA for SiBr(3) is 3.29 eV. The other predicted electron affinities are for SiH(3) (1.37 eV), SiH(2)CH(3) (1.09 eV), SiH(2)F (1.54 eV), SiH(2)Cl (1.94 eV), SiH(2)Br (2.05 eV), SiH(2)(SiH(3)) (1.77 eV), SiH(CH(3))(2) (0.92 eV), SiHF(2) (1.86 eV), SiHCl(2) (2.53 eV), SiHBr(2) (2.67 eV), Si(CH(3))(3) (0.86 eV), SiF(3) (2.66 eV), SiCl(3) (3.21 eV), Si(SiH(3))(3) (2.25 eV), and SiFClBr (3.13 eV). For the five silyl radicals where experimental data are available, the BLYP functional gives the most accurate determination of AEAs; the average absolute error is 0.04(1) eV, whereas the corresponding errors for the O3LYP, MPW1K, PBE1PBE, B3LYP, OLYP, and BHLYP functionals are 0.05(8), 0.06(0), 0.06(3), 0.08(5), 0.11(5), and 0.15(3) eV, respectively.