### 精确模型化学中MPn微扰外推的一个问题

1. 西北工业大学化工系 西安 710072|西北大学现代物理所 西安 710068
• 收稿日期:2000-02-21 修回日期:2000-05-09 发布日期:2000-10-15
• 通讯作者: 苏克和 E-mail:sukehe@nwpu.edu.cn

### Single-reference M?ller-Plesset Perturbation Extrapolation Technique in the Accurate Model Chemistry is Unreliable for Non-equilibrium Molecules

Su Ke-He, Wei Jun, Hu Xiao-Ling, Yue Hong, Lv Ling, Wang Yu-Bin, Wen Zhen-Yi

1. Department of Chemical Engineering,Northwestern Polytechnical University,Xi'an 710072|Institute of Modern Physice,Northwest University,Xi'an Shanxi 710068
• Received:2000-02-21 Revised:2000-05-09 Published:2000-10-15
• Contact: Su Ke-He E-mail:sukehe@nwpu.edu.cn

Abstract:

Potential energy curves of BH(1Σ ) and AlH(1Σ ) molecules in their ground electronic states were calculated at RQCISD(T), UQCISD(T), RCCSD(T), UCCSD(T) and CASSCF (with four active electrons in five active orbitals) theoretical levels. Basis sets examined were 6-311＋ G(3df,2p),G3 large ,aug-cc-pVDZ and aug-cc-pVTZ and the bond lengths tested were up to 1.0 nm. Potential energy curves of BF(1Σ ) and AlF(1Σ ) were also examined at MP4/aug-cc-pVDZ level. It was shown that the MPn (n=2－ 4) curves for BH, AlH and BF were seriously in error. In accordance with CASSCF wavefunctions, the single Hartree-Fock reference state was insufficient to describe the electronic structure of those molecules. In this case, the Hartree-Fock energy gaps between the LUMO and HOMO orbitals for these molecules were decreased when the bond distances were increased. By examination of the molecular orbital integrals, it was found that the interaction (integral values) between HOMO and LUMO orbitals increased greatly. As in the case of second-order MP perturbation calculations, the "stronger interaction" and the "smaller energy gap" resulted in too large correlation energy corrections and therefore led to serious errors in the total electronic energy. This predicted clearly that the MPn extrapolation technique in the widely used G1,G2,G3 and CBS series accurate model chemistry was not reliable in general if the molecule is not at its equilibrium geometry, as what in the study of chemical reaction paths. However, QCISD(T) calculations reflected approximately the correct potential energy curves and more reliable dissociation limits. The maximum error,40.6 kJ• mol－ 1,of dissociation energy at the bond length of 1.0 nm was found in AlH molecule by the G2 basis sets 6-311＋ G(3df,2p).The CCSD(T) was also found less reliable than the QCISD(T) method in describing the dissociation energies. The same error was 101.9 kJ• mol－ 1.All of the results in this work suggested the non-extrapolation G2(QCI) be more reliable a method in the energy calculations of non-equilibrium structure of molecules and be more valuable in further developments and modifications to the methods of accurate model chemistry.