Finite-volume and thermodynamic-limit lower bounds for the coherently driven Bose--Hubbard model via the Lieb--Yamazaki method

技術要約：Lieb–Yamazaki法による、コヒーレント駆動ボース–ハバード模型の有限体積および熱力学的極限における下界

問題提起
本論文は、コヒーレント駆動ボース–ハバード模型の基底状態エネルギーに対する厳密な下界を得るという課題に取り組んでいる。変分試行状態は容易に上界を提供する一方で、四次オンサイト相互作用と（粒子数保存を破る）コヒーレント源項の両方が存在する場合には、厳密な下界を確立することが困難である。単一サイトや結合クラスターによる初等的な推定などの既存の手法は、モデルのパラメータ（デチューニング、ホッピング、駆動）への明示的な依存性を欠いているか、あるいは熱力学的極限における観測量や励起スペクトルの解析のための安定した枠組みを提供できないことが多い。著者らは、熱力学的極限に対して安定であり、明示的なパラメータ依存性を保持し、かつマクロな観測量を導出するための解ける補助系を提供する、厳密な下界を構築することを目的としている。

手法
著者らは、任意の次元 $d \ge 1$ における周期的な超立方格子上の駆動ボース–ハバード・ハミルトニアンに対し、Lieb–Yamazakiの正演算子分解を適用している。手法は以下のステップで進行する：

1. 正演算子分解: 無次元のハミルトニアン $h$ を、正定値演算子の和 $w^2 = \sum B_k^\dagger B_k$ と残差ハミルトニアン $h_1$ に分解する。演算子 $B_k$ は、実変分パラメータ $K$ および $g$ を用いて以下のように定義される：
   $$B_k = a_k^\dagger a_k - K - g(a_k - a_k^\dagger)$$
   これにより、演算子の不等式 $h \ge h_1$ が得られる。

2. 残差の対角化: 残差 $h_1$ は二次形式のハミルトニアンである。これは以下の手順で対角化される：

   • 線形項を除去するために、一様な実変位 $\lambda$ を用いる。
   • 二次項を対角化するために、ボゴリューボフ変換を用いる。
     このプロセスにより、「境界ハミルトニアン」である $h_1$ と、その基底状態であるガウス状態 $|0_B\rangle$ が定義される。

3. 安定性条件: 分解および対角化の妥当性には、パラメータ $K$ および $g$ に関する特定の安定性制約が必要である。具体的には、ボゴリューボフ変換が正準であり、正のスペクトルを与えるためには、すべての運動量モード $\mu$ に対して $\omega_\mu > 2g^2$ である必要がある。ここで $\omega_\mu$ は $K, g$, ホッピング $j$, およびデチューニング $\nu_0$ に依存する。

4. 熱力学的極限: 著者らは、格子和をブリルアンゾーン積分に変換することにより、下界エネルギー密度 $e_{0,l}^{(d)}$ および連続励起スペクトル $\Omega(q)$ の熱力学的極限の存在を証明する。

5. テンプル不等式の解析: ガウス状態における全ハミルトニアンの分散を用いて下界を精緻化するために、テンプル不等式の適用可能性を調査する。これにより、不等式自体は有限体積では厳密であるが、エネルギー分散の規模がスペクトルギャップの強度に対して広範（extensive）であるため、システムサイズが大きくなるにつれて条件がますます制限的になることを確立している。

主要な貢献および結果

• 厳密な下界 (定理 II.3): 著者らは、任意の次元 $d \ge 1$ において、明示的かつ厳密な基底状態エネルギーの下界を導出している。この下界は、安定性制約に従う $K$ と $g$ によってパラメータ化されている。この下界の熱力学的極限は、ブリルアンゾーン上の閉じた形式の積分として表現される。
• 励起スペクトル: 境界ハミルトニアン $h_1$ は、連続的な励起バンド $\Omega(q)$ を持つ。
  • 安定領域において、スペクトルは有限のエネルギーギャップと二次分散（質量を持つ準粒子）を持つ。
  • 許容領域の境界において、ギャップは閉じ、分散は線形（$\Omega(q) \propto |q|$）となり、その傾きはペアリングパラメータ $g$ とホッピング $j$ によって決定される。著者らは、この線形な傾きは「境界ハミルトニアン」を特徴付けるものであり、必ずしもフル駆動モデルを特徴付けるものではないことを強調している。
• マクロな観測量: 境界ハミルトニアンのガウス基底状態 $|0_B\rangle$ 内において、著者らは以下の量の厳密な表現を導出している：
  • 秩序パラメータ（変位 $\lambda$）。
  • 充填率（粒子密度）。
  • エネルギー密度（運動およびオンサイト）。
  • 粒子数分散および二次相関関数 $g^{(2)}(0)$。
  • 異常平均。
    これらすべての量は、明確な熱力学的極限を持つ。
• アンチバンチング: 解析により、十分に強いコヒーレント駆動が、補助的なガウス状態内においてアンチバンチング（$g^{(2)}(0) < 1$）およびサブポアソン統計を誘起することを示す。これは、変位項がゆらぎを克服することによって引き起こされる。
• 数値的評価:
  • 著者らは、Lieb–Yamazakiの下界を、初等的なサイトおよび結合クラスターの下界と比較している。結合クラスターの下界の方が数値的に鋭い場合が多いものの、Lieb–Yamazakiの構成は、解析的に解ける補助状態と明示的なパラメータ依存性を持つという利点があることを指摘している。
  • 変分パラメータ $(K, g)$ の数値的最適化により、変分上界（同じパラメータで評価されたもの）と下界との相対的なギャップは、強い駆動（$f \gtrsim 3$）および中程度のホッピングにおいて最も狭くなることが明らかになった。
  • システムサイズ $M=4$ における有限サイズ効果は、エネルギー・ブラケットの幅に比べれば小さいことが判明しており、これは「下界の緩さ」が有限サイズによるアーティファクトではなく、主にガウス近似の限界に起因することを示している。
  • テンプル不等式の条件（$M^d \Delta_M < \Omega(0)$）は、試行状態が厳密な基底状態に極めて近い場合を除いて、大きなシステムでは成立しないことが示されており、一般的な熱力学的補正としての有用性に制限がある。

意義および主張
本論文は、任意の次元で有効であり、明確な熱力学的極限を持つ、駆動ボース–ハバード模型に対する変分的な厳密下界の族を確立したと主張している。主な意義は、以下の機能を持つ明示的に解ける補助ガウス状態（$h_1$）の構築にある：

1. 基底状態エネルギーに対する、認定された下限閾値。
2. 境界状態内において、観測量および相関関数に対する厳密な解析的表現。
3. 駆動によるアンチバンチングの領域を特定するための枠組み。

著者らは、下界の精度については控えめな姿勢をとっている。ガウス近似は、強い非ガウス相関（例：安定性の境界付近や、対称性の破れにより再現されない未駆動モデルのモット領域における相関）を捉えることができないことを認めている。したがって、これらの下界は一様に最適であるとは主張されていない。実際、いくつかのテストケースでは、初等的な結合クラスターの下界の方が数値的に鋭い。本研究の価値は、高精度な数値計算手法（DMRGやテンソルネットワークなど）に代わるものではなく、物理パラメータに明示的に依存する、解析的に扱いやすい厳密な下界を提供することにある。得られた結果は、境界ハミルトニアンの性質（例：ギャップレスな境界の傾き）を特徴付けるものであり、これらの特徴がフル相互作用モデルに対応するかどうかについては、有限サイズスケーリングやその他の数値的手法を用いた独立した検証が必要であると、著者らは明記している。