[論文レビュー] Hybrid Monte Carlo Simulation of Graphene on the Hexagonal Lattice

Richard C. Brower, C. Rebbi+1 more|arXiv (Cornell University)|Jan 26, 2011

Graphene research and applications参考文献 2被引用数 7

ひとこと要約

本稿では、連続近似を回避するため、六角格子上のグラフェンにおける直接的ハイブリッドモンテカルロ（HMC）シミュレーションを提示する。タイトビンディングハミルトニアンを最近接 hopping とクーロン相互作用を、グラスマン経路積分フレームワークで定式化する。2サイト系に対して正確な結果を再現し、SU(2)対称な連続極限に収束することを示し、電子相関を持つより大きな系に対しても本手法の妥当性を検証する。

ABSTRACT

We present a method for direct hybrid Monte Carlo simulation of graphene on the hexagonal lattice. We compare the results of the simulation with exact results for a unit hexagonal cell system, where the Hamiltonian can be solved analytically.

研究の動機と目的

• 先行の格子ゲージ理論的手法で用いられる連続近似に依存せずに、グラフェンの本来の六角格子上での直接的HMCシミュレーションを開発すること。
• クーロン相互作用をグラスマン変数とハッブル＝ストラトニヴィッチ変換を用いて表現する経路積分定式化を確立すること。
• 調整可能な電子間相互作用を持つ最小の2サイト系に対して、正確な解と照合して手法を検証すること。
• 連続時間極限におけるSU(2)アイソスピン対称性の回復をテストし、スピン簡約状態の場合にフェルミオンの符号問題が生じないことを確認すること。
• 電子相関、格子歪み、外部場を伴うより大きなグラフェン系のシミュレーションの基盤を築くこと。

提案手法

• AおよびBサブラットスに属するフェルミオンの生成/消滅演算子を用いて、最近接 hopping とクーロン相互作用項を含む六角格子上のタイトビンディングハミルトニアンを定式化する。
• 分割関数と相関関数を、二次的フェルミオン作用と四次相互作用項を持つグラスマン経路積分の形で表現する。
• ハッブル＝ストラトニヴィッチ変換を適用し、実の補助ボソン場を介して四次相互作用を二次形式に変換し、フェルミオンに関するガウス積分を可能にする。
• 補助場とフェルミオン行列式の有効作用を導出し、空孔を時間反転粒子として取り扱う際のスピン成分間のキャンセルにより符号問題が生じないことを保証する。
• 分子動力学更新を用いたハイブリッドモンテカルロアルゴリズムを実装し、時間スライスパラメータ δ = β/Nt を用いてユークリッド時間の離散化を制御する。
• L=4 の有限周期的六角格子上でのシミュレーションを実施し、連続時間極限（Nt = 32, 64, 128, 256）における収束性をテストする。

実験結果

リサーチクエスチョン

• RQ1 連続近似に依存せずに、グラフェンの本来の六角格子上にハイブリッドモンテカルロシミュレーションを直接適用可能か？
• RQ2 HMC手法は連続時間極限においてタイトビンディングハミルトニアンのSU(2)アイソスピン対称性を保存するか？
• RQ3 本手法は、電子間相互作用あり・なしの2サイトグラフェン系について、正確な結果を正確に再現できるか？
• RQ4 スピン簡約状態下におけるフェルミオン行列式の振る舞いはいかなるものか？また、符号問題を引き起こすか？
• RQ5 時間格子間隔 δ → 0 のとき、相関関数およびエネルギー準位はどのように正確な連続極限値に収束するか？

主な発見

• HMCシミュレーションは、e=0（自由）およびe=0.5（相互作用あり）の両ケースにおいて、相関関数 C_a(t) が正確な結果に収束し、Nt=128 で優れた一致を示す。
• C_a(t) に対してSU(2)回転で関連付けられる相関関数 C_b(t) も同じ連続極限結果に収束するため、時間連続極限における対称性回復が確認された。
• HMCデータからの線形補外により、lim_{δ→0} E = 1.262 ± 0.004 を得た。これは正確な連続極限エネルギー E_0 = e^2 + √(4 + e^4) - 1 ≈ 1.266 と整合的である。
• 空孔を時間反転粒子として取り扱う際、上スピンと下スピンの行列式間のキャンセルにより、HMC手法に符号問題が生じないことが示された。
• e=0.5 かつ r0=0.5 の場合でも手法は安定かつ正確に機能し、通常順序化則により化学ポテンシャル μ = e²/r0 が適切に設定されると、SU(2)対称性が保たれる。
• 本手法はより大きな系へ拡張可能であり、次近接 hopping（κ′/κ ≈ 0.05）は再重み付けにより取り扱えると予想される。