[論文レビュー] Single-cone real-space finite difference schemes for the Dirac von Neumann equation

Magdalena Schreilechner, W. Pötz|arXiv (Cornell University)|Mar 6, 2015

Topological Materials and Phenomena参考文献 10被引用数 3

ひとこと要約

本稿では、フェルミオン二重化を回避するためのステアガードグリッドを用いて、(2+1)Dディラックハミルトニアンのフォンノイマン方程式を解くための2つの単一コーン実空間有限差分スキームを提示する。ブラケットスキームは、ケットとブラにハミルトニアンを別々に適用することで正値性を保つが、直接スキームは1ステップで密度行列の時間微分を計算する。両スキームとも単一コーン分散を確保し、安定性を有しており、トポロジカル絶縁体表面状態への応用およびリンドブラッドおよびグリーン関数形式への拡張が可能である。

ABSTRACT

Two finite difference schemes for the numerical treatment of the von Neumann equation for the (2+1)D Dirac Hamiltonian are presented. Both utilize a single-cone staggered space-time grid which ensures a single-cone energy dispersion to formulate a numerical treatment of the mixed-state dynamics within the von Neumann equation. The first scheme executes the time-derivative according to the product rule for "bra" and "ket" indices of the density operator. It therefore directly inherits all the favorable properties of the difference scheme for the pure-state Dirac equation and conserves positivity. The second scheme proposed here performs the time-derivative in one sweep. This direct scheme is investigated regarding stability and convergence. Both schemes are tested numerically for elementary simulations using parameters which pertain to topological insulator surface states. Application of the schemes to a Dirac Lindblad equation and real-space-time Green's function formulations are discussed.

研究の動機と目的

• トポロジカル絶縁体表面状態を含む、(2+1)Dディラック系における混合状態ダイナミクスの数値的シミュレーションという課題に取り組む。
• 1次空間時間微分に対する標準的な有限差分スキームに内在するフェルミオン二重化問題を克服する。
• グリッドサイズに線形にスケーリングする、効率的で安定的かつ正値性を保つフォンノイマン方程式の数値スキームを開発する。
• ステアガード格子上でのペイエルス置換を用いて、電磁ポテンシャルのゲージ不変な取り扱いを可能にする。
• リンドブラッド方程式によるオープン量子系への拡張および実空間時間グリーン関数形式への応用を可能にするフレームワークを構築する。

提案手法

• エネルギー運動量分散関係に単一のディラックコーンを保証するため、単一コーンステアガード空間時間グリッドを採用する。
• ブラケットスキームは、密度演算子の左（ブラ）および右（ケット）添え字にハミルトニアンを別々に適用し、交換子に積の法則を適用する。
• 時間微分を1回のスイープで密度行列に計算する直接スキームを開発し、別々のブラ・ケット伝搬を回避する。
• ステアガードグリッド上で1格子間隔の中心差分を用いることで、精度を維持するとともにフェルミオン二重化を回避する。
• ゲージ不変な形で電磁ベクトルポテンシャルを組み込むためにペイエルス置換を適用する。
• トポロジカル絶縁体表面状態に適したパラメータを用いて数値的検証を行い、安定性および収束性を分析する。

実験結果

リサーチクエスチョン

• RQ1 フォンノイマン方程式の有限差分スキームは、(2+1)Dディラック系において正値性を保ち、フェルミオン二重化を回避できるか？
• RQ2 ブラケットスキームの構造は、純粋状態ディラック方程式からどのように安定性と正値性を継承するか？
• RQ3 密度行列の時間発展に対する直接有限差分スキームの安定性および収束性はいかほどか？
• RQ4 提案されたスキームは、リンドブラッド方程式によるオープン量子系ダイナミクスに拡張可能か？
• RQ5 非平衡輸送を扱うために、実空間時間グリーン関数形式にフレームワークを適合できるか？

主な発見

• ブラケットスキームは、積の法則の実装により、純粋状態ディラック方程式から正値性と安定性を継承する。
• 直接スキームは、1ステップあたりの行列演算数を減らすことで計算効率を向上させるが、ブラケットスキームに比べて収束半径が小さい。
• 両スキームともステアガードグリッドを用いることでフェルミオン二重化を排除し、エネルギー分散に単一のディラックコーンを保証する。
• トポロジカル絶縁体表面状態に適したパラメータを用いた数値的テストにおいて、両スキームは安定かつ収束することが確認された。
• ペイエルス置換は、ステアガードグリッド上でゲージ不変な形で適用可能であり、電磁ポテンシャル下でもスキームの安定性を保つ。
• リンドブラッド方程式およびグリーン関数形式への拡張は可能であり、散乱効果および非平衡輸送のモデル化が可能になる。