面向重元素吸收边精确RIXS谱图：一种相对论四分量与精确两分量TDDFT方法

《Journal of Chemical Theory and Computation》：Toward Accurate RIXS Spectra at Heavy Element Edges: A Relativistic Four-Component and Exact Two-Component TDDFT Approach

　　

研究人员提出了一种基于全四分量（4c）Dirac–Coulomb哈密顿量与现代原子平均场精确两分量（amfX2C）哈密顿量模型的相对论含时密度泛函理论（TDDFT）方法，用于模拟共振非弹性X射线散射（RIXS）谱图。该方法建立在伪波函数（pseudo-wavefunction）形式体系与核-价分离（core–valence separation）方案之上，能够高效计算相对于同一基态的两组激发态流形之间的耦合，从而求解RIXS所需的Kramers–Heisenberg方程。相对论表述对标量相对论（scalar relativistic）与自旋-轨道（spin–orbit）效应提供变分描述，这对准确描述RIXS过程涉及的内壳层激发至关重要。通过amfX2C哈密顿量将问题变换至两分量（2c）体系，在计入X2C变换带来的双电子与交换-关联（exchange–correlation）图像变化（picture-change）效应后，显著降低了计算成本，同时给出接近4c质量的结果。除二维RIXS图谱外，该方法还可直接求取高能分辨荧光探测（HERFD）与共振X射线发射谱（RXES）。对选定钌与铀配合物的2p3d与3d4f RIXS图谱的应用表明，amfX2C方法以高精度复现了参考4c结果与实验谱图，捕捉到所有关键谱学特征并给出可靠的峰指认。

研究背景方面，薛定谔在建立波动力学时曾试图纳入相对论但未获成功，狄拉克随后给出相对论性电子方程，解释了光谱精细结构与自旋起源并预言正电子，二人因此共享1933年诺贝尔物理学奖。在多电子体系中引入电子-电子相互作用及Gaunt与Breit等高阶修正后，狄拉克-库仑（-Breit）方程成为相对论量子化学的黄金标准，即全四分量（4c）方法。为降低开销，精确两分量（X2C）哈密顿量通过代数解耦将4c问题缩为2c，其中Knecht等人发展的原子平均场（amfX2C）与扩展原子平均场（eamfX2C）模型在计算效率与精度间取得良好平衡，并能在Kohn–Sham DFT中纳入双电子与交换-关联（XC）的图像变化效应。RIXS作为吸收高能X光子后再发射低能光子的过程，其能量转移反映价层或内壳层电子结构，具备原子选择性且分辨率不受芯空位寿命限制；随着同步辐射与自由电子激光发展，分子气相与溶液相RIXS成为可能，但重元素L、M吸收边因p、d芯轨道受自旋-轨道分裂，必须采用多分量相对论形式。现有高阶波函数方法虽准但昂贵，非相对论TDDFT无法处理自旋-轨道效应，故研究人员开展本研究以实现4c与amfX2C两套相对论TDDFT框架下的RIXS二维谱及HERFD、RXES一维截面的高效高精度计算，论文发表于《Journal of Chemical Theory and Computation》。

关键技术方法上，研究人员采用伪波函数形式体系将LR-TDDFT激发振幅近似为单激子行列式线性组合，并在Tamm–Dancoff近似（B=0）下结合核-价分离方案，避免二阶响应理论的高成本与伪极点问题；在4c层面求解Dirac–Kohn–Sham方程获得四分量分子轨道系数与能量，构造含Dirac单电子、Coulomb双电子与非共线交换-关联贡献的A、B矩阵；进而通过酉解耦矩阵将4c Fock矩阵块对角化，取正能++块得到amfX2C–KS方程，使2c计算继承4c的双电子与XC图像变化修正；样本为文献中已有的钌配合物与铀配合物体系（对应2p3d与3d4f RIXS图谱）。

研究结果部分，引言（Introduction）指出相对论量子化学从4c到X2C的演进及RIXS对重元素内壳层激发的必要性，明确本研究将Nascimento等人的非相对论TDDFT-RIXS框架推广至相对论域。理论（Theory）一节给出各向同性平均下的Kramers–Heisenberg截面公式（式1）、极化依赖跃迁幅F_fn(θ)（式2）及跃迁矩乘积S^uv_fn（式3），由LR-TDDFT本征向量Xⁿ、Yⁿ抽取基态-激发态偶极矩（式4、5、6）；为避免激发态间矩阵元的二阶响应，引入伪波函数|k?（式7、8）并在TDA下化简为单激展开（式9），推导得中间态-末态耦合表达式（式10），在核-价分离下化简为仅含占据-占据块的紧凑形式；随后列出4c框架中A^4c、B^4c矩阵元（式11）、Dirac–KS方程（式12）及Dirac单电子h^D、Coulomb G^4c、非共线XC势V^4c_xc的具体构成（式13a–13c）；最后说明X2C通过酉矩阵U将F^4c解耦为块对角（式14），舍弃负能块后得amfX2C–KS方程，使2c TDDFT自动携带X2C变换后的双电子与XC图像变化积分。应用结果表明，对钌与铀配合物的2p3d、3d4f RIXS图谱，amfX2C在峰位、分裂结构与强度上与参考4c计算及实验高度一致，可可靠指认谱峰。

讨论与结论部分，研究人员总结：所实现的4c与amfX2C相对论TDDFT-RIXS流程，首次在变分意义上统一纳入标量相对论与自旋-轨道耦合，借助伪波函数与核-价分离绕开激发态间矩阵的二阶响应瓶颈；amfX2C在仅约一半维度下复现4c结果，使重元素L、M边RIXS预测在中等尺寸配合物上实用化；该方法同时输出二维RIXS、HERFD与RXES，为多电子强关联的开壳层过渡金属与锕系体系提供兼顾效率与精度的谱学解释工具，验证了X2C图像变化修正对双电子与交换-关联项不可忽略。