[论文解读] Compressed Wannier modes found from an $L_1$ regularized energy functional

Farzin Barekat, Ke Yin+4 more|arXiv (Cornell University)|Mar 26, 2014

Optical Polarization and Ellipsometry参考文献 2被引用 3

一句话总结

本文提出了一种新颖的方法，通过$L_1$-正则化能量泛函直接计算紧支撑的Wannier函数，利用单一参数$\mu$系统性地权衡局域化与能量精度。该方法通过凸优化确保平移正交性，生成稀疏的实值Wannier模态，无需预先计算布洛赫态，从而实现高效且精确的电子结构计算。

ABSTRACT

We propose a method for calculating Wannier functions of periodic solids directly from a modified variational principle for the energy, subject to the requirement that the Wannier functions are orthogonal to all their translations ("shift-orthogonality"). Localization is achieved by adding an $L_1$ regularization term to the energy functional. This approach results in "compressed" Wannier modes with compact support, where one parameter $μ$ controls the trade-off between the accuracy of the total energy and the size of the support of the Wannier modes. Efficient algorithms for shift-orthogonalization and solution of the variational minimization problem are demonstrated.

研究动机与目标

• 开发一种不依赖布洛赫态表示的直接变分方法来计算Wannier函数。
• 通过单一正则化参数$\mu$，系统性地控制Wannier函数局域化与总能量精度之间的权衡。
• 通过凸优化框架确保Wannier函数的平移正交性，避免传统局域化泛函中常见的局部极小值问题。
• 利用与现有DFT代码兼容的迭代求解器，实现Wannier模态的高效数值计算。
• 将压缩平面波框架推广至具有通用晶体势的周期性固体，包括金属体系。

提出的方法

• 通过最小化$L_1$-正则化能量泛函$\mathcal{J}(\psi) = \frac{1}{\mu}\|\psi\|_1 + \langle\psi|\hat{H}|\psi\rangle$来促进稀疏性和局域化。
• 通过在每个能带下对$\mathcal{J}(\psi)$施加平移正交性和归一化约束，递归地计算Wannier函数。
• 通过基于傅里叶的条件强制实现平移正交性：对所有$\mathbf{k} \in \text{BZ}$，有$N\sum_{\mathbf{G}'}|\tilde{\psi}(\mathbf{k}+\mathbf{G}')|^2 = \frac{1}{|\Omega|}$，确保与所有晶格平移的正交性。
• 哈密顿算符作用$\hat{H}\psi$通过布洛赫分解$\hat{H}\psi = \sum_{\mathbf{k}} e^{i\mathbf{k}\mathbf{r}} \hat{H}_{\mathbf{k}} u_{\mathbf{k}}(\mathbf{r})$计算，实现与标准DFT代码的兼容性。
• 使用高效的迭代求解器求解凸最小化问题，利用现有$\hat{H}_{\mathbf{k}}u_{\mathbf{k}}$的计算例程。
• 该方法将压缩平面波基推广至周期性体系，实现紧支撑性与通过$\mu$的系统性控制。

实验结果

研究问题

• RQ1 Wannier函数能否在不先计算布洛赫态的情况下，直接从变分原理中计算得出？
• RQ2 如何利用$L_1$正则化实现具有可控局域化与能量精度的紧支撑Wannier函数？
• RQ3 在周期性体系中，平移正交性在确保Wannier函数正确物理结构方面起什么作用？
• RQ4 如何在数值上稳定地、高效地在正交性约束下最小化$L_1$-正则化能量泛函？
• RQ5 该方法能否推广至金属体系，并集成到现有的DFT工作流中？

主要发现

• $L_1$-正则化能量泛函生成的Wannier函数具有紧支撑性，即仅在有限空间区域内非零，这是$L_1$范数特性的结果。
• 参数$\mu$提供了能量精度与局域化之间的系统性权衡：$\mu$越小，Wannier函数越局域化，但能量精度越低。
• 该方法通过涉及倒易晶格矢量上平方布洛赫振幅之和的傅里叶空间条件，确保了Wannier函数的平移正交性。
• 泛函的凸性使得其能够实现高效且鲁棒的数值最小化，避免了传统非凸局域化泛函中常见的局部极小值问题。
• 计算复杂度与标准布洛赫基DFT方法相当，因为$\hat{H}\psi$可通过标准布里渊区采样和傅里叶变换进行评估。
• 该方法适用于绝缘体和金属体系，突破了以往仅限于拓扑平凡绝缘体的方法限制。