Analytical solution of a free-fermion chain with time-dependent ramps

技術要約：時間依存するランプ（傾斜）を受ける自由フェルミオン鎖の解析解

問題提起
本論文は、時間依存する線形ポテンシャルの影響下にある、非相互作用フェルミオン（ジョルダン・ワグナー変換により、相互作用のないボゾンとも等価）の一次元鎖の、非平衡ダイナミクスを扱っている。この系は、時間とともに任意に変化する傾き $\xi(t)$ を持つタイトバインディング・ハミルトニアンによって記述される。静的な線形ポテンシャルについては、ブロッホ振動やワニエール・スターク局在化を含めよく理解されているが、特定の時間依存ケース（急激なクエンチや周期駆動など）の研究はあるものの、任意の時間依存駆動プロトコルに対する一般的な厳密な解析解は欠けていた。著者らは、この空白を埋めるべく、この系の単一粒子シュレディンガー方程式の厳密解を提供し、粒子密度、電流、およびもつれエントロピーといった観測量の厳密なダイナミクスを導出することを目的としている。

手法
著者らは、格子上の第 $k$ モードの単一粒子シュレディンガー方程式 $\psi_k(j, t)$ （ポテンシャルは $j/\xi(t)$）から出発する。静的な系の固有基底（ベッセル関数を含む）に着想を得て、著者らは時間発展した波動関数の厳密なアンザッツとして以下を提案する：
$$\psi_k(j, t) = \exp\left(i(j - k)\theta(t) - ik\phi(t)\right) J_{j-k}[\ell(t)]$$
ここで、$J_\nu$ は第一種ベッセル関数であり、$\ell(t)$、$\theta(t)$、$\phi(t)$ は決定されるべき時間依存関数である。

このアンザッツをシュレディンガー方程式に代入し、ベッセル関数の漸化式を利用することで、著者らは $\ell(t)$ と $\theta(t)$ に関する結合非線形微分方程式のセットを導出し、$\phi(t)$ を逆傾斜の位相積分として特定する。任意の駆動関数 $\xi(t)$ に対してこれらの結合方程式を解くために、著者らは補助変数 $u(t) = \ell(t)\cos\theta(t)$ および $v(t) = \ell(t)\sin\theta(t)$ を導入する。この変換により、問題は二つの線形常微分方程式へとデカップリングされ、これらは駆動プロトコル $\xi(t)$ と累積位相 $\phi(t)$ を含む積分を用いた閉形式の解を持つ。

主要な貢献と結果

1. 厳密な解析解： 主要な貢献は、任意の $\xi(t)$ に対する厳密な時間発展単一粒子波動関数の導出である。これにより、厳密な時間発展二点相関関数（等時刻グリーン関数）の構成が可能となる：
   $$\tilde{C}_{i,j}(t) = e^{i(i-j)\theta(t)} C_{i,j}(\ell(t))$$
   ここで $C_{i,j}(\ell)$ は、傾き $\ell$ を持つ静的系の基底状態相関を表す。この結果は、多体ダイナミクスが自己相似的な構造を保持しており、界面幅 $\ell(t)$ と位相シフト $\theta(t)$ によって進化が完全に特徴付けられることを示している。

2. 流体力学的解釈： 傾きが大きい極限（$\xi \gg 1$）において、著者らは占有関数 $n_p(x,t)$ に関するオイラー方程式を用いた流体力学的記述を導出する。彼らは、厳密解が、フェルミ点 $p_\pm(x,t)$ が時間依存する長さスケール $\ell(t)$ と位相 $\theta(t)$ によって決定される流体力学的状態に対応することを示す。これは、格子解をドメインウォール（領域壁）の融解問題へと結びつけ、静的な時間変数 $t$ を動的なスケール $\ell(t)$ に置き換えるものである。

3. 観測量： 厳密な相関関数を用いて、著者らは以下の明示的な表現を導出している：

   • 粒子密度： $\rho(x,t) \approx \frac{1}{\pi} \arccos\left(\frac{x}{\ell(t)}\right)$
   • 電流： $J(x,t) \approx \frac{\sin\theta(t)}{\pi} \sqrt{1 - \frac{x^2}{\ell(t)^2}}$
   • もつれエントロピー： $S(x,t) \approx \frac{1}{6} \log\left[ \ell(t) \left(1 - \frac{x^2}{\ell(t)^2}\right)^{3/2} \right] + \Upsilon$
     著者らは、もつれエントロピーは $\ell(t)$ にのみ依存し、位相 $\theta(t)$ には依存しないことを指摘しており、計算をスケーリングされた時間を用いたドメインウォール融解問題へと還元している。

4. 特殊な極限：

   • 断熱極限： 駆動が遅い場合（$\dot{\xi} \to 0$）、界面幅 $\ell(t)$ は小さな振動補正を伴いつつ、駆動傾斜 $\xi(t)$ に断熱的に追従する。
   • 量子クエンチ： 傾きが $\xi_0$ から $\xi_1$ へ急激に変化する場合、著者らはドメインウォール融解や文献に見られる特定の極限（例：参考文献 [8, 18]）を回収する。
   • ワニエール・スターク局在化： クエンチ極限において、界面幅 $\ell(t)$ は「呼吸（breathing）」挙動を示し、周期的に振動する。著者らはこれを、連続系で見られる無制限なドリフトを阻止するブリルアンゾーンの離散的な性質による、ワニエール・スターク局在化の実現として解釈している。これは、相関領域の周期的な圧縮と膨張をもたらす。

意義と主張
本論文は、任意の時間依存する傾斜を持つ、駆動されたスターク局在系のダイナミクスに対する最初の一般的な厳密解を提供すると主張している。この厳密解を確立することにより、本研究は流体力学的および数値的なアプローチに対する厳密なベンチマークを提供し、駆動と局在化の相互作用に関する新たな知見を与える。特に、著者らは、クエンチ極限で見られる呼吸する界面が、以前は流体力学的な議論を通じてのみ示唆されていた現象である、非平衡設定におけるワニエール・スターク局在化の具体的な実現を提供することを強調している。これらの結果は、特定のプロトコル（急激なクエンチやドメインウォール融解など）に関する既存の知見を一般化するものであり、相互作用のある系やもつれハミルトニアンに関する将来の研究への基礎を築くものである。著者らは、自らの研究が、導出された単一粒子ダイナミクスを通じてアクセス可能な、駆動された量子多体系の新しい厳密解を構成していることを明示している。