Departamento de Química, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
J Chem Phys. 2013 Feb 21;138(7):074102. doi: 10.1063/1.4790166.
At the end of the lanthanide series, 4f → 5d and other interconfigurational transitions, in which one electron is excited from a tight 4f orbital to a much more diffuse one, occur with a break of many f-f pairs, which make the electron correlation effects dominant. For instance, the large energy gap of 25 000 cm(-1) (∼29 500 cm(-1) without spin-orbit coupling) above the 4f(14) ground state of the SrCl2:Yb(2+) material is mostly due to electron correlation. In effect, a minimal multiconfigurational restricted active space (RASSCF) calculation that includes only the 4f(14) ground and 4f(13)5d and 4f(13)6s open-shell excited configurations gives a very small gap (5400 cm(-1)), whereas the correlation corrections to the 4f(14) → 4f(13)5d(eg) transition energies at the second order perturbation theory (RASPT2) level are very large: 35 599 ± 439 cm(-1), in average, for all excited states. These corrections are too large to be accurate at second order perturbation level. When a second f-shell is also included in the active space and single and double excitations to the 5d, 6s, and 5f shells are treated variationally, the (extended) RASSCF energy gap above the ground state and the electronic transitions increase by 22 038 ± 120 cm(-1) and the RASPT2 correlation energy corrections become small (-721 ± 571 cm(-1)), as it is desirable for a second order perturbation. A comparative analysis of both RASPT2 results reveals that the lack of the second f-shell accounts for 12 700 cm(-1) of the 14 223 ± 80 cm(-1) overestimation of interconfigurational transitions energies by the minimal RASPT2 calculation, which indicates an inaccurate calculation of the differential radial correlation between the 4f(14) and 4f(13)5d configurations by second order perturbation theory. In order to establish practical and accurate procedures for the calculation of 4f → 5d and other interconfigurational transitions at the end of the lanthanide series, the above and other RASSCF/RASPT2 calculations on the ionization potential of Yb(2+) in gas phase and in SrCl2 have been benchmarked in this paper against coupled cluster (coupled cluster singles and doubles and triples ) calculations, and RASSCF/RASPT2 calculations on the absorption spectrum of SrCl2:Yb(2+) have been compared with experiment. The results support that variational calculation of SD 4f → 5f excitations prior to RASPT2 calculations can be a realistic, accurate, and feasible choice to model radial correlation effects at the end of the lanthanide series.
在镧系元素系列的末端,4f→5d 和其他组态间跃迁,其中一个电子从紧密的 4f 轨道被激发到更扩散的轨道,伴随着许多 f-f 对的断裂,这使得电子相关效应占主导地位。例如,SrCl2:Yb(2+)材料的 4f(14)基态上方的 25000 cm(-1)(无自旋-轨道耦合时约为 29500 cm(-1))的大能隙主要归因于电子相关。实际上,仅包含 4f(14)基态和 4f(13)5d 和 4f(13)6s 开壳激发态的最小多组态受限活性空间(RASSCF)计算给出了一个非常小的能隙(5400 cm(-1)),而二阶微扰理论(RASPT2)水平上对 4f(14)→4f(13)5d(eg)跃迁能的相关修正非常大:所有激发态的平均为 35599±439 cm(-1)。这些修正在二阶微扰水平上太大而无法准确。当第二个 f 壳层也包含在活性空间中并且对 5d、6s 和 5f 壳层进行单和双激发的变分处理时,基态以上的(扩展)RASSCF 能隙和电子跃迁增加 22038±120 cm(-1),而 RASPT2 相关能量修正变得很小(-721±571 cm(-1)),这是二阶微扰所期望的。对两种 RASPT2 结果的比较分析表明,最小 RASPT2 计算中对组态间跃迁能量的 14223±80 cm(-1)的过度估计中,缺少第二个 f 壳层占 12700 cm(-1),这表明二阶微扰理论对 4f(14)和 4f(13)5d 构型之间的微分径向相关性的计算不准确。为了建立在镧系元素末端计算 4f→5d 和其他组态间跃迁的实用和准确程序,本文对气相中的 Yb(2+)的电离势和 SrCl2 进行了 RASSCF/RASPT2 计算,并与耦合簇(耦合簇单双和三重)计算进行了基准测试,并将 SrCl2:Yb(2+)的吸收光谱的 RASSCF/RASPT2 计算与实验进行了比较。结果表明,在进行 RASPT2 计算之前,对 SD 4f→5f 激发的变分计算可以是一种现实、准确和可行的选择,用于模拟镧系元素末端的径向相关效应。