[論文レビュー] Monte Carlo Simulations of Quantum Spin Systems in the Valence Bond Basis

Anders W. Sandvik, K. S. D. Beach|ArXiv.org|Apr 11, 2007

Physics of Superconductivity and Magnetism参考文献 1被引用数 6

ひとこと要約

本稿では、結合状態（VB）基底におけるプロジェクター・モンテカルロ法を提案し、量子スピン系の高効率な基底状態シミュレーションと有限運動量の三重項励起状態への直接的アクセスを可能にする。シミュレーション中に試行状態を繰り返し最適化し、VB基底の性質を活用することで、スイング・トリプレットギャップと三重項分散関係を高精度に計算する。16×16格子においては正確な対角化結果と統計誤差の範囲で一致し、既知のスピン波速度 c = 1.58 と良好に一致する。

ABSTRACT

We discuss a projector Monte Carlo method for quantum spin models formulated in the valence bond basis, using the S=1/2 Heisenberg antiferromagnet as an example. Its singlet ground state can be projected out of an arbitrary basis state as the trial state, but a more rapid convergence can be obtained using a good variational state. As an alternative to first carrying out a time consuming variational Monte Carlo calculation, we show that a very good trial state can be generated in an iterative fashion in the course of the simulation itself. We also show how the properties of the valence bond basis enable calculations of quantities that are difficult to obtain with the standard basis of Sz eigenstates. In particular, we discuss quantities involving finite-momentum states in the triplet sector, such as the dispersion relation and the spectral weight of the lowest triplet.

研究の動機と目的

• 標準的な S^z 固有状態の定式化に起因する制限を克服し、量子スピン系に対してバリエンス・ボンド基底におけるプロジェクター・モンテカルロ法を開発すること。
• 従来の基底定式化では取りにくい、有限運動量の三重項励起状態（分散関係やスペクトル重みなど）を直接計算可能にすること。
• 高価な変分モンテカルロの事前計算を回避するため、試行状態の反復的自己最適化スキームを導入すること。
• スイング・トリプレットギャップとスピン波速度の高精度推定を、スイング状態と三重項状態の共同サンプリングによる誤差相殺によって達成すること。

提案手法

• ハミルトニアンの負の冪 (-H)^m を用いて、バリエンス・ボンド基底で試行状態から基底状態をプロジェクションする。m を増加させることで収束を達成する。
• サブラティス A と B 間のスイング状態に制限を加えた制限付き VB 基底を用い、計算コストを低減しながらも本質的な物理を保持する。
• VB 構造のサンプリングと S^z_q 演算子を含む行列要素の計算により、E(q) やスペクトル重みといった運動量空間の性質にアクセスする。
• エネルギーおよび行列要素の式を再定式化し、初期状態にのみ依存するサンプリング重みを分離することで、運動量依存観測量の効率的測定を可能にする。
• VB 状態間の重なりルール（ループ数による）を活用し、シミュレーション中に行列要素を効率的に計算する。
• スイング状態と三重項状態の共同サンプリングにより、スイング・トリプレットギャップ Δ と三重項分散 ω(q) を計算し、誤差相殺を実現する。

実験結果

リサーチクエスチョン

• RQ1 標準的な変分最適化を経ずに、S=1/2 ハイゼンベルグ反強磁性体の基底状態をバリエンス・ボンド基底におけるプロジェクター・モンテカルロ法で効率的に計算できるか？
• RQ2 有限運動量の三重項励起状態（分散関係やスペクトル重みなど）を、VB 基底で直接計算できるか？
• RQ3 スイング状態と三重項状態の共同サンプリングにより、スイング・トリプレットギャップとスピン波速度を高精度に推定できるか？
• RQ4 VB 基底は、動的構造因子行列要素 ⟨τ(q)|S^z_q|σ(0)⟩ のような物理量にどのようにアクセスを可能にするか？
• RQ5 三重項結合長分布は、励起スペクトルにおけるスピンオンの脱コンfinement やトリプレオン結合の特徴をどのように示すか？

主な発見

• 投影長 N/m における指数的収束が達成され、L=16 格子において大きな m で明確な収束が観察され、L=4 時に正確な対角化結果と一致する。
• スイング・トリプレットギャップ Δ が高精度に計算され、q₁ = (π−2π/L, π) における三重項エネルギーからスピン波速度 c = 1.58 が得られ、既知の値と良好に一致する。
• VB 基底により、運動量空間の行列要素（特に (π,π) 近傍での単一磁気励起重み）の直接計算が可能となり、動的構造因子に寄与する。
• スイング状態と三重項状態の共同サンプリングにより、誤差相殺が顕著に発生し、独立計算と比較してギャップ Δ の統計的精度が桁違いに向上する。
• 三重項結合長分布は、スピンオンの脱コンfinement に一致する特徴を示しており、拡散的分布である可能性を示唆しているが、さらなる解析が進行中である。