镱硫属化合物的密度泛函理论研究

用密度泛函理论(DFT)研究镱硫属化合物的电子结构和性质,通过与实验比较考察了现有的几种近似密度泛函公式对镧系元素化合物的适用程度和相对论效应的影响.结果表明,用DFT计算的YbO键长对实验值的偏差约为0.002nm;但得到的键能即使在考虑梯度校正和相对论效应之后,仍比实验值高,在定域密度近似基础上引入交换梯度校正使键能计算值减小,其中PW86x使键能计算值减小稍多些,结果更接近实验值;相关梯度校正使键能计算值升高.相对论效应使键长缩短0.004～0.006nm,键能减小约0.5eV.计算结果的分析表明,Yb的5d轨道和配体的np轨道间形成σ键和π键.在所研究的分子体系中,配体原子从O到Te、Yb原子的5d轨道布后数依次减少,同键能减弱的顺序一致.相对论效应使键能减小的主要原因是在成键过程中发生了Yb的6s电子向5d轨道的转移,而相对论效应使该过程能量增加.偶极矩和电荷分布的计算表明,Yb-L键以共价性为主,相对论效应使共价性成份增加.     

等电子三原子铀化物OUO2+、NUN和NUO+的结构和谐振频率的CCSD(T)计算研究

用小核相对论有效势和CCSD(T)方法计算了三原子铀化物OUO²⁺, NUN和NUO⁺的平衡键长和谐振频率. 计算结果显示U原子内层5s5p5d 电子相关能对这些化合物性质的影响非常小. 除NUN的弯曲振动频率,旋轨耦合效应对这些化合物的结构和频率的影响并不明显. 本文的计算结果与其他研究组的计算结果以及已有的实验值相比符合较好, 这表明作为单参考态方法, CCSD(T)能够对这些体系的键长和频率给出较精确的计算结果. 与此前密度泛函理论(DFT)的计算结果相比, CCSD(T)方法与PBE0泛函的结果吻合最好. 本文的工作有助于在用密度泛函方法研究这些体系时选择合适的交换相关泛函, 也为今后的实验研究提供了新的理论数据.     

N—O键解离焓的密度泛函理论研究

在高精度计算方法G3和G3B3的基础上,比较了密度泛函理论(DFT)十几种方法对N—O键解离焓(BDE)相对于实验值的计算精度,发现用B3P86方法计算15种化合物N—O键的BDE,均方根误差最小,仅为6.36kJ·mol-1,计算值与实验值的线性相关系数为0.991.在此基础上,用该方法分别计算了非芳香化合物及芳香化合物的N—O键BDE.通过自然键轨道分析,发现部分N—O键的BDE与N—O键的键长、原子电荷密度及键级之间存在定量关系.此外,在B3P86方法的基础上预测了几种典型的杂环芳香化合物N—O键BDE值.     

密度泛函理论研究镥二聚体低电子能态的相对稳定性、结构以及成键性质

用密度泛函理论方法研究了镥二聚体（Lu2）低能量电子态的性质,计算了电子态相对能量、平衡键长、振动频率以及基态解离能,考察了密度泛函性质、相对论有效势种类以及Hartree—Fock交换作用大小对计算结果的影响．结果表明,无论采用何种密度泛函和相对论有效势,体系的基态都为三重态,与其他一些基于分子轨道理论的从头计算方法得到的结论是一致的．另外,与分子轨道从头计算结果以及实验结果比较发现,采用杂化密度泛函理论和Stuttgart小核有效势计算得到的结果总体吻合最好．最后,特别分析研究了B3LYP计算中Hartree—Fock交换作用大小对基态键长和基态解离能的影响,发现随着交换作用的增大,键长增长,解离能减小,这是由于5d轨道杂化导致的共价成键作用减弱造成的．     

3,7-二硝亚胺基-2,4,6,8-四硝基-2,4,6,8-四氮杂双环[3.3.0]辛烷结构与性能的量子化学研究

设计了一种新型高能量密度化合物--3,7-二硝亚胺基-2,4,6,8-四硝基-2,4,6,8-四氮杂双环[3.3.0]辛烷, 应用密度泛函理论(DFT)的B3LYP 方法在6-31G(d,p)基组水平上对该化合物进行了结构全优化, 并计算得到其红外(IR)光谱; 通过键级分析获得热解引发键的位置为N7-N22, 同时求得校正后的键离解能为91.47 kJ/mol. 采用Monte-Carlo方法预测该化合物的理论密度为2.16 g/cm³; 基于理论密度并结合等键反应及Kamlet-Jacobs公式预测了生成焓、爆速、爆压和爆热值分别为1219.94 kJ/mol, 10.43 km/s, 53.44 GPa和7407.84 J/g. 以上性能参数显示, 该目标化合物达到了高能量密度化合物的基本要求, 是一种潜在的含能材料. 同时给出了该化合物的逆合成路线.     

元素电负性和硬度的密度泛函理论研究

应用密度泛函理论的DFT LDA、DFT LDA/NL和改进的Slater过渡态方法,把元素的电离能和电子亲合能的计算扩展到周期表的103种元素.并用有限差分方法计算了这103种元素的电负性和硬度.计算中考虑了相对论效应.计算结果比以前Robles等用密度泛函理论的XGL和Xα近似的交换相关泛函的计算结果有所改进,更接近实验值.     

N6H6结构和性质的理论研究

采用密度泛函理论的b3lyp方法在6-311＋＋G**基组上对15种分子式为N₆H₆的氮氢化合物进行了理论计算, 并且应用了自然键轨道理论(Nature Bond Orbital, NBO)和分子中的原子理论(Atoms In Molecules, AIM)分析了这些化合物的成键特征和相对稳定性. NBO分析表明N原子孤对电子到相邻的氮氮键的超共轭作用是影响氮氮键长变化的主要因素, AIM计算的氮氮键的键临界点电荷密度与键长呈反比关系. 而且, NBO的立体和超共轭分析表明立体交换排斥能和超共轭作用对这些分子的相对稳定性起了重要作用. G3MP2计算结果表明氮氢化合物的生成热均为正, 并且环状分子的能量和生成热都高于链状分子.     

N6H6结构和性质的理论研究

采用密度泛函理论的b3lyp方法在6-311＋＋G**基组上对15种分子式为N₆H₆的氮氢化合物进行了理论计算, 并且应用了自然键轨道理论(Nature Bond Orbital, NBO)和分子中的原子理论(Atoms In Molecules, AIM)分析了这些化合物的成键特征和相对稳定性. NBO分析表明N原子孤对电子到相邻的氮氮键的超共轭作用是影响氮氮键长变化的主要因素, AIM计算的氮氮键的键临界点电荷密度与键长呈反比关系. 而且, NBO的立体和超共轭分析表明立体交换排斥能和超共轭作用对这些分子的相对稳定性起了重要作用. G3MP2计算结果表明氮氢化合物的生成热均为正, 并且环状分子的能量和生成热都高于链状分子.     

密度泛函计算条件对结果精度的影响

通过对有代表性的分子BCl₃,SO₂,ZnO,TiCl₄,LuF₃的一系列计算,从三个方面考察密度泛函理论的Kohn-Sham方法中采用的计算条件对结果精度的影响.发现基组完备性对计算结果影响最大,数值积分取点数目对计算结果影响略小.电荷多极投影收敛较快.分子几何构型和振动基频对计算条件不太敏感,总能量和键能对计算结果较敏感.在计算键能时对分子及其组成原子采用相同的基组和积分布点方式可以减少误差.选择适当的计算条件,用较小的计算量可以得到与现有近似能量密度泛函精度相匹配的结果.     

Au(II)化合物[Au(CH_2)_2PH_2]_2X_2(X=F,Cl,Br,I)的量子化学理论研究

采用从头计算HF,MP2方法和密度泛函理论,对Au(II)系列化合物[Au(CH2)2PH2]2X2(X=F,Cl,Br,I)的几何结构、电子结构和振动频率进行了研究.研究表明Au的5d和6s电子参与Au—Au以及Au—X之间的成键.Au—Au,Au—X键强烈的电子相关作用使HF方法不适于该体系的研究,BP86和B3LYP两种泛函给出较大的Au—Au和Au—X键长,而MP2方法和局域的密度泛函方法则给出了合理的结构参数.局域密度泛函方法计算得到的Au—Au键和Au—X键振动频率也与实验数据符合较好.还运用含时密度泛函理论计算了[Au(CH2)2PH2]2X2的电子激发能,对分子在紫外-可见光谱范围内的电子跃迁进行了分析,考察了卤素配体对激发能的影响,并结合分子轨道能级的变化对此给予了解释.     

On low‐lying excited states of extended nanographenes

Low‐lying excited states of planarly extended nanographenes are investigated using the long‐range corrected (LC) density functional theory (DFT) and the spin‐flip (SF) time‐dependent density functional theory (TDDFT) by exploring the long‐range exchange and double‐excitation correlation effects on the excitation energies, band gaps, and exciton binding energies. Optimizing the geometries of the nanographenes indicates that the long‐range exchange interaction significantly improves the C C bond lengths and amplify their bond length alternations with overall shortening the bond lengths. The calculated TDDFT excitation energies show that long‐range exchange interaction is crucial to provide accurate excitation energies of small nanographenes and dominate the exciton binding energies in the excited states of nanographenes. It is, however, also found that the present long‐range correction may cause the overestimation of the excitation energy for the infinitely wide graphene due to the discrepancy between the calculated band gaps and vertical ionization potential (IP) minus electron affinity (EA) values. Contrasting to the long‐range exchange effects, the SF‐TDDFT calculations show that the double‐excitation correlation effects are negligible in the low‐lying excitations of nanographenes, although this effect is large in the lowest excitation of benzene molecule. It is, therefore, concluded that long‐range exchange interactions should be incorporated in TDDFT calculations to quantitatively investigate the excited states of graphenes, although TDDFT using a present LC functional may provide a considerable excitation energy for the infinitely wide graphene mainly due to the discrepancy between the calculated band gaps and IP–EA values. © 2017 Wiley Periodicals, Inc.     

Numerical examination of performance of some exchange-correlation functionals for molecules containing heavy elements

The performance of 17 exchange-correlation functionals for molecules containing heavy elements are numerically examined through four-component relativistic density DFT calculations. The examined functionals show the similar accuracy as they do for the molecules containing light elements only except for bond lengths. LDA and OP86 produce good results for bond lengths and frequencies but bad bond energies. Different functionals do not show much different performance for bond energies except LDA. BP86 and GP86 produce results with average accuracy while LYP does not perform well. Although encouraging results are obtained with functional B97GGA-1, other heavily parameterized and meta-GGA functionals do not produce impressive results.     

A systematic sequence of relativistic approximations

An approach to the development of a systematic sequence of relativistic approximations is reviewed. The approach depends on the atomically localized nature of relativistic effects, and is based on the normalized elimination of the small component in the matrix modified Dirac equation. Errors in the approximations are assessed relative to four-component Dirac-Hartree-Fock calculations or other reference points. Projection onto the positive energy states of the isolated atoms provides an approximation in which the energy-dependent parts of the matrices can be evaluated in separate atomic calculations and implemented in terms of two sets of contraction coefficients. The errors in this approximation are extremely small, of the order of 0.001 pm in bond lengths and tens of microhartrees in absolute energies. From this approximation it is possible to partition the atoms into relativistic and nonrelativistic groups and to treat the latter with the standard operators of nonrelativistic quantum mechanics. This partitioning is shared with the relativistic effective core potential approximation. For atoms in the second period, errors in the approximation are of the order of a few hundredths of a picometer in bond lengths and less than 1 kJ mol(-1) in dissociation energies; for atoms in the third period, errors are a few tenths of a picometer and a few kilojoule/mole, respectively. A third approximation for scalar relativistic effects replaces the relativistic two-electron integrals with the nonrelativistic integrals evaluated with the atomic Foldy-Wouthuysen coefficients as contraction coefficients. It is similar to the Douglas-Kroll-Hess approximation, and is accurate to about 0.1 pm and a few tenths of a kilojoule/mole. The integrals in all the approximations are no more complicated than the integrals in the full relativistic methods, and their derivatives are correspondingly easy to formulate and evaluate.     

密度泛函理论研究温度和取代基对4-亚芳基亚氨基-1,2,4-三唑-3(2H)-酮气相热分解的影响

采用密度泛函理论(DFT)方法, 在B3LYP/6-31＋＋G(d,p)水平上研究了4-亚芳基亚氨基-1,2,4-三唑-3(2H)-酮及其类似物的气相热分解反应. 从热力学性质、几何结构参数、自然电荷分布、温度效应等角度探讨不同取代基对标题物热分解反应的影响. 结果表明: R (R＝H, CH2CH2CN)取代基对反应影响很小; 而苯环上Y (Y＝MeO, Me, H, Cl和NO2)取代基对活化自由能(ΔG≠)、键长(ΔR≠)、键角(Δα≠)和自然电荷(Δq≠)变化的影响与Hammet常数(σ)呈现线性关系ΔP≠＝ρσ＋C (P＝G, R, α和q, r＞0.95). 给电子基团使得相应的活化自由能降低, 键长变短, 键角变大, 自然电荷增加; 吸电子反之.     

Crown ether inclusion complexes of the early actinide elements, [AnO2(18-crown-6)]n+, An = U, Np, Pu and n = 1, 2: a relativistic density functional study

The title compounds, [AnO2(18-crown-6)]n+, An = U, Np, and Pu and n = 1 and 2, as well as the related (experimentally observed) complex [UO2(dicyclohexyl-18-crown-6)]2+ are studied using relativistic density functional theory (DFT). Different relativistic methods (large-core and small-core effective core potentials, all-electron scalar four-component) and two flavors of approximate DFT (B3LYP and PBE) are used. Calculated bond lengths agree well with the available experimental data for the NpV complex, while larger differences for the UVI complexes appear to be related to the large uncertainties in the experimental data. The axial AnO bonds are found to be weaker and longer than in the corresponding penta-aquo complexes, though still of partial triple-bond character. The AnO bond lengths and strengths decrease along the actinide series, consistent with the actinide contraction. Gas-phase binding energies calculated for the penta-aquo complexes and crown-ether complexes of the actinides studied, as well as ligand-exchange energies, show that there is no intrinsic preference, or "better fit", for actinyl(V) cations as compared to actinyl(VI) ones. Rather, the ability of NpO2+ (NpV) to form in-cavity 18-crown-6 complexes in water, which is impossible for UO22+, is traced to solvation effects in polar solvents. Thus, the experimentally observed stabilization of the pentavalent oxidation state as compared to the hexavalent one is due to the effective screening of the charge provided by the macrocycle, and this leads to destabilization of the AnVI crown complexes relative to their AnV counterparts.     

外电场对甲醛分子的结构与电子光谱的影响

采用密度泛函理论(DFT)在B3LYP/6-311++G(d,p)基组水平上,计算了不同外加电场(-8.22×10~9~8.22×10~9 V/m)下甲醛分子基态稳定构型、分子键长、电荷分布、能级分布、能隙、红外光谱、拉曼光谱和分子的总能量.在此基础上利用TDDFT/B3LYP/6-311++G(d,p)方法研究了甲醛分子由基态跃迁到前25个激发态的激发能E、谐振强度f、吸收波长λ受外电场的影响.结果表明:随着C=O连线方向外电场的增加,C=O键键长、氢原子电荷、偶极矩和能隙递增;C—H键键长、C,O原子电荷递减,总能量降低.振动频率与红外强度及拉曼强度由于不同振动有不同变化.甲醛分子UV-Vis光谱随外电场的增加,不同的吸收峰发生了不同程度的蓝移或者红移;外电场对甲醛分子的激发能、谐振强度和吸收波长的强度有一定影响,但随电场变化比较复杂.     

Meta-generalized gradient approximation: explanation of a realistic nonempirical density functional

Tao, Perdew, Staroverov, and Scuseria (TPSS) have constructed a nonempirical meta-generalized gradient approximation (meta-GGA) [Phys. Rev. Lett. 91, 146401 (2003)] for the exchange-correlation energy, imposing exact constraints relevant to the paradigm densities of condensed matter physics and quantum chemistry. Results of their extensive tests on molecules, solids, and solid surfaces are encouraging, suggesting that this density functional achieves uniform accuracy for diverse properties and systems. In the present work, this functional is explained and details of its construction are presented. In particular, the functional is constructed to yield accurate energies under uniform coordinate scaling to the low-density or strong-interaction limit. Its nonlocality is displayed by plotting the factor F(xc) that gives the enhancement relative to the local density approximation for exchange. We also discuss an apparently harmless order-of-limits problem in the meta-GGA. The performance of this functional is investigated for exchange and correlation energies and shell-removal energies of atoms and ions. Non-self-consistent molecular atomization energies and bond lengths of the TPSS meta-GGA, calculated with GGA orbitals and densities, agree well with those calculated self-consistently. We suggest that satisfaction of additional exact constraints on higher rungs of a ladder of density functional approximations can lead to further progress.     

Local exchange-correlation approximations and first-row molecular dissociation energies

Recent Xα calculations of bond energies and other related properties of first-row diatomic molecules show very encouraging agreement with experiment. In the worst cases, however, the Xα dissociation energies overestimate the experimental values by almost 2 eV. Therefore, we have examined several refinements of the Xα theory and their effects on molecular bond lengths, bond energies, and vibrational frequencies. Among them, gradient corrections to the Xα exchange energy and also some variations of the local spin-density correlation energy approximation are considered. We find that a local exchange-correlation functional with gradient corrections gives dissociation energies in significantly better agreement with experiment than the Xα approximation.