[论文解读] On the pure state $v$-representability of density matrix embedding theory

Fabian M. Faulstich, Raehyun Kim+4 more|arXiv (Cornell University)|Oct 18, 2021

Physics of Superconductivity and Magnetism参考文献 57被引用 4

一句话总结

该论文提出了一种增广拉格朗日方法（alm-DMET），以放松密度矩阵嵌入理论（DMET）中非相互作用纯态 v-可表示性（NI-PS-V）假设，从而在阿夫堡原理失效时，仍能实现高阶和低阶密度矩阵块的精确匹配。通过允许任意轨道占据分布并保持幂等性，该方法避免了低阶哈密顿量中的虚假无隙性，显著提升了在二维 Hubbard 模型和氢链模型中的精度。

ABSTRACT

Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy.

研究动机与目标

• 解决由于非可行的 NI-PS-V 假设导致的 DMET 中密度矩阵精确匹配失败的问题。
• 解决当阿夫堡原理被违反时，低阶能带隙消失的问题。
• 开发一种稳健的方法，使密度矩阵块的精确匹配不依赖于阿夫堡原理。
• 在标准 DMET 失效的强关联体系中，提升数值稳定性和精度。
• 在挑战性条件下，展示所提方法在二维 Hubbard 模型和氢链模型中的有效性。

提出的方法

• 使用增广拉格朗日方法将相关势拟合表述为约束优化问题，以强制实现密度矩阵的精确匹配。
• 通过允许低阶密度矩阵中任意占据分布，放松 NI-PS-V 条件，同时保持幂等性（占据数被限制为 0 或 1）。
• 使用投影梯度法求解约束优化问题，引入惩罚参数和线性搜索以实现收敛。
• 引入超参数（如初始惩罚值、最大内层迭代次数）以控制收敛性并降低计算成本。
• 限制内层迭代次数，以减少完整体系对角化的次数，降低计算开销。
• 每 10 次外层迭代采用自适应超参数更新策略，以改善收敛行为。

实验结果

研究问题

• RQ1 当高阶密度矩阵块不满足 NI-PS-V 条件时，能否在 DMET 中实现密度矩阵的精确匹配？
• RQ2 在低阶密度矩阵构建中放松阿夫堡原理，是否能消除低阶哈密顿量中的虚假无隙性？
• RQ3 增广拉格朗日方法能否提升强关联体系中 DMET 的数值鲁棒性和收敛性？
• RQ4 与标准最小二乘法 DMET 相比，所提 alm-DMET 方法在精度和收敛性方面的表现如何？
• RQ5 超参数配置对 alm-DMET 的收敛性和计算成本有何影响？

主要发现

• 在具有挑战性的体系中，如空穴掺杂的二维 Hubbard 模型（n = 32/36），阿夫堡原理的放松实现了密度矩阵块的精确匹配。
• alm-DMET 方法成功避免了标准 DMET 在 NI-PS-V 条件失效时出现的低阶能带隙消失问题。
• 数值结果表明，与标准最小二乘法 DMET 相比，alm-DMET 在能量和密度矩阵匹配方面均表现出更高的精度，尤其在强关联区域。
• alm-DMET 的收敛性与初始 RHF 初始猜测无关，这一特性在 H36 链体系中通过不同初始化设置得到验证。
• 限制内层迭代次数可降低计算成本并稳定收敛过程，但最优超参数调优对性能仍至关重要。
• 不同超参数设置下，能量轨迹完全一致，表明该方法在能量收敛方面具有高度鲁棒性，即使迭代次数不同。