[論文レビュー] 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$正則化エネルギー関数 $\mathcal{J}(\psi) = \frac{1}{\mu}\|\psi\|_1 + \langle\psi|\hat{H}|\psi\rangle$ を直接最小化することで、コンpactにサポートされたWannier関数を計算する新規な手法を提案する。この手法により、1つのパrameter $\mu$ を用いて局在化とエネルギー精度の間で体系的なトレードオフを実現できる。凸最適化によりシフト直交性が保証され、事前のBloch状態計算を経ずにスパースで実数値の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.

研究の動機と目的

• Bloch状態の表現に依存せずに、Wannier関数を直接変分的に計算する手法の開発。
• 1つの正則化パrameter $\mu$ を用いて、Wannier関数の局在化と全エネルギー精度の間の体系的制御を実現すること。
• 非凸な従来の局在化関数の手法で一般的に生じる局所最適解を避けるために、凸最適化フレームワークを用いてWannier関数のシフト直交性を保証すること。
• 既存のDFTコードと互換性のある反復解法を活用して、Wannierモードの効率的数値計算を可能とすること。
• 一般の格子ポテンシャル、金属系を含む周期的格子系に対しても、圧縮平面波フレームワークを一般化すること。

提案手法

• スパarsityと局在化を促進するため、$L_1$正則化エネルギー関数 $\mathcal{J}(\psi) = \frac{1}{\mu}\|\psi\|_1 + \langle\psi|\hat{H}|\psi\rangle$ を最小化する。
• 各バンドに対して、シフト直交性および正規化制約のもとで $\mathcal{J}(\psi)$ を逐次最小化することでWannier関数を計算する。
• シフト直交性は、フーリエ空間の条件 $N\sum_{\mathbf{G}'}|\tilde{\psi}(\mathbf{k}+\mathbf{G}')|^2 = \frac{1}{|\Omega|}$ （$\mathbf{k} \in \text{BZ}$ に対してすべての $\mathbf{k}$ で成立）により強制され、すべての格子並進に対して直交することを保証する。
• ハミルトニアン作用 $\hat{H}\psi$ は、Bloch分解 $\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 Bloch状態の計算を経ずに、変分原理から直接Wannier関数を計算することは可能か？
• RQ2 どのようにして $L_1$ 正則化を用いて、局在化とエネルギー精度を制御可能なコンパクトにサポートされたWannier関数を実現できるか？
• RQ3 周期的系におけるWannier関数の正しい物理的構造を保証するために、シフト直交性が果たす役割は何か？
• RQ4 数値的に安定に、直交性制約のもとで $L_1$-正則化エネルギー関数を効率的に最小化する方法は何か？
• RQ5 この手法は金属系に一般化可能であり、既存のDFTワークフローに統合可能か？

主な発見

• $L_1$ノルムの性質により、$L_1$正則化エネルギー関数は、有限の空間領域内でのみ非ゼロとなるコンパクトにサポートされたWannier関数を生成する。
• パrameter $\mu$ は、エネルギー精度と局在化の間の体系的トレードオフを提供する：$\mu$ が小さいほど、より局在化するがエネルギー精度は低下する。
• シフト直交性は、逆格子ベクトルにわたるBloch振幅の二乗和を含むフーリエ空間の条件により保証される。
• 関数の凸性のおかげで、局所最適解に陥りやすい従来の非凸な局在化関数の手法とは異なり、効率的かつ安定した数値的最小化が可能になる。
• $\hat{H}\psi$ の評価が標準的なブリユアンゾーンサンプリングとフーリエ変換を用いて可能であるため、計算複雑度は標準的なBlochベースのDFT手法と同等である。
• 本手法は、トポロジカルに自明な絶縁体に限定されず、絶縁体および金属系の両方に対応可能である。