[論文レビュー] A phase-space method for the Bose-Hubbard model

Poonam Jain, C. W. Gardiner|arXiv (Cornell University)|Apr 27, 2004

Cold Atom Physics and Bose-Einstein Condensates参考文献 14被引用数 3

ひとこと要約

本稿では、Bose-Hubbardモデルの基底状態を近似するため、Q関数表現に基づく位相空間法を提案する。1サイトおよび2サイトの平均場近似を用いる。量子的制約、特に不確定性関係が位相図を支配しており、正確な解と良好な定性的一致を示すが、パラメータ化の制限により2サイトモデルはMott絶縁体転移を正しく予測できない。

ABSTRACT

We present a phase-space method for the Bose-Hubbard model based on the Q-function representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard "one site" model where quantum tunneling is approximated entirely using mean-field terms; the second "two site" model explicitly includes tunneling between two adjacent sites while treating tunneling with other neighbouring sites using the mean-field approximation. The ground state is determined by minimizing the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose-Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime.

研究の動機と目的

• 超流動相およびMott絶縁体相の両方において、Bose-Hubbardモデルの基底状態を計算的に効率よく近似する手法を開発すること。
• 量子力学的制約、特に不確定性関係が、モデルの位相構造を決定する役割を調査すること。
• 平均場近似のもとで、1サイトおよび2サイトモデルに対してQ関数位相空間法と正確な数値解を比較すること。
• 2サイトQ関数アプローチがMott絶縁体相転移を正しく捉えられない理由を評価すること。
• 2サイトを超えるサイト数に拡張することで、DMRG や QMC などの数値的正確な結果との一致がどのように向上するかを検討すること。

提案手法

• Q関数表現を用いて、Bose-Hubbardハミルトニアンをガウス変数の形に再パラメータ化し、古典的位相空間記述を可能にする。
• 量子的制約（数の分散に関する1つと、数と位相のフラクチュエーションを結ぶ1つ）から導かれる不確定性関係に基づき、古典的エネルギー関数を最小化することで基底状態を求める。
• 2つのモデルハミルトニアンを検討する：1サイトモデル（トンネル項を完全に平均場的に扱う）、2サイトモデル（相互作用サイト間のトンネルを明示的に含み、他の近接サイトについては平均場近似を適用）。
• 1サイトモデルは、数の分散と平均場振幅の変分最適化により正確に解け、Gutzwillerアンザッツに等価である。
• 2サイトモデルはQ関数法と、格子の対称性を仮定した正確な数値的最小化の両方で解かれる。
• 結果は、特に1次元格子において、DMRG や量子モンテカルロ法などの正確な数値的手法と比較される。

実験結果

リサーチクエスチョン

• RQ1 Q関数に基づく位相空間法は、Bose-Hubbardモデルにおける超流動からMott絶縁体への転移を正確に記述できるか？
• RQ2 量子力学的不確定性制約は、1サイトおよび2サイトの平均場近似において、基底状態構造にどのように影響を与えるか？
• RQ3 2サイトQ関数法は、良好な定性的一致を示すにもかかわらず、なぜMott絶縁体相を正しく予測できないのか？
• RQ4 Q関数法と正確な数値的手法との間で、転移の臨界相互作用強度 $\overline{U}_c$ はどう異なるか？
• RQ5 2サイトモデルをより多くのサイトに拡張することで、DMRG や QMC などの正確な結果との一致はどの程度向上するか？

主な発見

• Q関数法は、超流動からMott絶縁体への転移を定性的に正確に記述しており、1サイトモデルでは通常、正確な解と5%以内の一致を示す。
• 1サイトQ関数法は正確な数値的解と良好に一致し、計算時間のわずか数分の1で同等の結果を得られる。
• 2サイトQ関数法は、Mott絶縁体相を正しく捉えていない。正確な2サイト解とは異なり、転移点で平均場量がゼロにならないためである。
• 正確な2サイト解では、整数占有状態の下で明確な位相転移が観察され、平均場量がゼロとなる。これはMott絶縁体相と整合的である。
• 転移の臨界相互作用強度 $\overline{U}_c$ は、2サイト正確解において1サイトモデルよりも低く、特に1次元系で顕著なシフトを示す。
• 結果から、特にMott絶縁体相の始まりに近い段階で、量子的制約（特に数の分散の下限）が基底状態を決定する上で支配的役割を果たしていることが示唆される。