[論文レビュー] Hybrid Monte-Carlo simulation of interacting tight-binding model of graphene

Dominik Smith, Lorenz von Smekal|arXiv (Cornell University)|Nov 5, 2013

Graphene research and applications参考文献 12被引用数 3

ひとこと要約

本研究では、長距離クーロン相互作用をモデル化するために Hubbard-Stratonovich 変換を用いた、グラフェンの相互作用的タイトバインディング模型に対するハイブリッドモンテカルロシミュレーションを実施する。有効微細構造定数 α_C ≈ 4 における金属絶縁体相転移を同定し、先行研究と整合的であることを示し、コンパクトな Hubbard 細胞と2次精度離散化スキームが、高次の補正項を含んでも同一の連続極限をもたらすことを示した。

ABSTRACT

In this work, results are presented of Hybrid-Monte-Carlo simulations of the tight-binding Hamiltonian of graphene, coupled to an instantaneous long-range two-body potential which is modeled by a Hubbard-Stratonovich auxiliary field. We present an investigation of the spontaneous breaking of the sublattice symmetry, which corresponds to a phase transition from a conducting to an insulating phase and which occurs when the effective fine-structure constant $α$ of the system crosses above a certain threshold $α_C$. Qualitative comparisons to earlier works on the subject (which used larger system sizes and higher statistics) are made and it is established that $α_C$ is of a plausible magnitude in our simulations. Also, we discuss differences between simulations using compact and non-compact variants of the Hubbard field and present a quantitative comparison of distinct discretization schemes of the Euclidean time-like dimension in the Fermion operator.

研究の動機と目的

• 非摂動的格子シミュレーションを用いて、グラフェンの相互作用的タイトバインディング模型における自発的サブラットスイッチ対称性の破れを調査すること。
• 金属絶縁体相転移が発生する有効微細構造定数 α_C の臨界値を特定すること。
• 異なる離散化スキームおよび Hubbard 細胞のバリエーション（コンパクト対非コンパクト）が、シミュレーションの精度および収束性に与える影響を評価すること。
• 一貫したポテンシャルパラメータと境界条件を用いて、先行結果の α_C を再現できるか、シミュレーションフレームワークを検証すること。

提案手法

• タイトバインディングハミルトニアンの経路積分表現を導出し、Hubbard-Stratonovich補助場を用いて長距離2体ポテンシャルを組み込む。
• cRPA 値を用いて、局所項、最近接項、準最近接項、ヘキサゴン跨ぎ項に分けて非局所的で区分的なポテンシャルを構築し、スクリーニングのないクーロン尾部を含む。
• GPUアクセラレーションを用いたハイブリッドモンテカルロ法を用い、ユークリッド時間離散化に基づく 6×6 の長方形グラフェン格子上でシミュレーションを実施する。
• 複数の配置においてオーダーパラメータ ⟨Δ_N⟩ を測定し、自発的対称性の破れを検出するとともに、ゼロ質量極限 (m→0) に外挿する。
• 標準的および改良型2次フェルミオン作用の両方に対して、⟨Δ_N⟩ を 1/N_t および 1/N_t² 項にフィットさせることで、離散化誤差を定量化する。
• 収束性および非ゼロのオーダーパラメータ検出の観点から、非コンパクトな場合に比べてコンパクトな Hubbard 細胞が好ましい。

実験結果

リサーチクエスチョン

• RQ1 グラフェンの相互作用的タイトバインディング模型は、物理的に妥当な α_C において金属から絶縁体への相転移を示すか？
• RQ2 コンパクトと非コンパクトな Hubbard 細胞の定式化は、サブラットスイッチ対称性の破れのシミュレーションにどのように影響を与えるか？
• RQ3 1次精度と2次精度の離散化スキームが、フェルミオン作用の連続極限に与える定量的影響は何か？
• RQ4 一貫したポテンシャルパラメータと境界条件を用いた場合、シミュレーションフレームワークは先行研究の α_C の結果を再現できるか？

主な発見

• シミュレーションにより、α_C ≈ 4 における金属絶縁体相転移が同定され、先行研究が報告した α_C ≈ 3.12 と定性的に整合的である。
• 非コンパクトな Hubbard 細胞では非ゼロのオーダーパラメータが得られず、数値的不安定性および収束性の悪さが示された。
• 1次および2次精度の両方の離散化スキームが、同じ連続極限に収束し、係数 c1 および c2 の差が無視できるほど小さいことが確認された。
• 2次精度スキームは標準的スキームに比べて顕著な改善をもたらさず、この設定では高次項の利点がないことが示唆された。
• 結果は、今後の境界処理の改善やポテンシャルマッチングを想定した大規模な研究に向けた信頼性の高いプラットフォームとしてのシミュレーションフレームワークの妥当性を裏付けた。
• FFT を用いたポテンシャル計算の採用により、将来的にはより大きなスケールのシミュレーションと一貫した境界条件の実現が可能になる。