Long-time fermionic quantum transport with controlled full-state error using an adaptive reservoir-mode window

技術要約：適応型リザーバー・モード・ウィンドウを用いた、制御された全状態誤差を伴う長時間のフェルミオン量子輸送

問題提起
相互作用するナノ構造がフェルミオン・リザーバーと結合している実時間シミュレーションは、有限のデバイスとマクロな環境との間で生じるエンタングルメントの継続的な増大という、根本的な計算上のボトルネックに直面している。定常的な非相互作用輸送はランダウアー・ビュッティカー形式によって効率的に記述できるが、相互作用、駆動、または強非マルコフ的な系に関する実時間ダイナミクスには、開いた系のダイナミクスの明示的な取り扱いが必要となる。HEOM（階層型方程式法）、TEMPO（時間発展行列積演算子）、あるいは連鎖写像テンソルネットワークといった既存の系統的に収束する手法は、階層の深さ、メモリ長、またはボンド次元の切り捨てに依存している。要求の厳しい領域では、長時間にわたって規定の精度を維持するために、計算リソースの無制限な増大を余儀なくされることが多く、これがシミュレーションを短時間または中間的な時間に限定し、定常状態へのアクセスを困難にしている。

手法：テープレコーダー・コース・グレイン化（Tape-Recorder Coarse Graining）
著者らは、「テープレコーダー・コース・グレイン化」と呼ばれる手法を導入する。これは、各非相互作用リードを、流入（incoming）、アクティブ（active）、および**流出（outgoing）**という3つの明確なカテゴリのモードへと再構成する手法である。このアプローチは、モードがデバイスの進化に対して物理的に関連性を持つ場合にのみ保持される、時間依存のリザーバー基底を利用している。

1. 結合重みによるモード分類:

   • 結合重み: 本手法は、リザーバー・モードをデバイスへの時間積分された結合強度によって特徴付ける。単一粒子ヒルベルト空間において、2つの演算子が定義される：経過結合重み（$I^+$）および、規定されたシミュレーション区間 $[0, T]$ における残存結合重み（$I^-$）。
   • 流入モード: 累積結合重み（$I^+$）が最大重みの規定された割合（$r_{cut}$）を超えると、モードはアクティブな部分空間に進入する。
   • 流出モード: 区間 $[t, T]$ における残存結合重み（$I^-$）が同じ閾値を下回ると、モードはアクティブな部分空間から退出する。
   • アクティブ・ウィンドウ: 「アクティブ」なモードは、有限の移動窓を形成する。シミュレーションが進むにつれて、新しいモードが進入（流入）し、古いモードが退出（流出）することで、明示的に伝播される部分空間の次元が有界に保たれる。

2. フェルミオンへの拡張と粒子・正孔写像:

   • 本手法はフェルミオン輸送へと拡張されており、フェルミオンの反交換関係と有限密度初期状態という2つの課題に対処している。
   • ゼロ温度かつ最大バイアス設定（左リードが充填され、右リードが空の状態）の場合、粒子・正孔変換を用いることで、充填された左リードを準粒子真空へと写像できる。これにより、両方のリザーバーを準粒子表現において空として扱うことが可能になり、ライトコーン構造のモード空間伝播が保持される。

3. ダイナミクスと確率的アンラベリング:

   • システムは、デバイス、アクティブなリザーバー・モード、および基底回転に由来する運動項を含む有効ハミルトニアンの下で進化する。
   • モードが「流出」状態になると、その残留結合は切り捨てられる。このモードを伝播されるフォック空間から除去しつつ、簡約されたダイナミクスを保持するために、著者らはモードの占有数の確率的サンプリングを採用している。これは、流出モードに対する厳密な部分トレースを表しており、これによりデバイスとアクティブなモードが、モンテカルロ・トラジェトリの枠組み内で純粋状態の波動関数として進化を継続することを可能にする。

4. 誤差制御:

   • 流出モードの切り捨てに対する、厳密な状態と切り捨てられた全状態との間のインフィデリティ（不純度）に関する非摂動的な上界が導出されている。
   • この上界は、モードの残存結合重みと、保持されたダイナミクスの有限区間応答因子に依存する。決定的なのは、この上界が、弱い瞬時結合や流出モードの特定の占有数を仮定していない点である。

主要な結果とベンチマーク
本手法は、ゼロ温度・最大バイアス条件下で、半無限タイトバインディング・リードに結合した2サイト相互作用量子点接触（QPC）を用いてベンチマークが行われた。

• ローレンツ型リザーバー: 本手法は、収束したHEOMの計算およびランダウアー・ビュッティカーの定常状態電流と一致する。これにより、HEOMが効率的であることが知られている領域における本アプローチの妥当性が検証された。
• フラットバンド・リザーバー: 代数的に減衰する相関（フラットなスペクトル密度）を持つリザーバーにおいて、本手法は有限サイズの再帰が発生する前の直接的なシュレディンガー発展の結果を再現し、ランダウアー・ビュッティカーの定常電流とも一致した。対照的に、同じ領域でテストされた有限指数型のHEOM展開は収束せず、鋭いスペクトル特徴を指数関数和で表現することの困難さを浮き彫りにした。
• 相互作用系: 本手法は、相互作用するQPCの電流－電圧特性におけるクーロン・ブロックADE・ピーク分裂を正常に捉えた。
• 駆動系: 非相互作用駆動の極限において、結果は厳密なフロケ・グリーン関数計算と一致し、特定の駆動振幅におけるコヒーレントなトンネル抑制（コヒーレント破壊トンネリング）を再現した。この抑制は、有限のクーロン斥力の存在下でも持続する。
• スケーリング: 数値解析により、固定された相対閾値（$r_{cut}$）において、アクティブなモードの数は時間とともに飽和することが示された。アクティブな幅は、閾値を厳しくした場合、時間の経過に対して線形ではなく、対数的（$\ln(1/r_{cut})$）に増大する。これにより、シミュレーションの長い時間窓においても、ピークとなるヒルベルト空間の次元が有界に保たれることが保証される。

意義と主張
本論文は、テープレコーダー・コース・グレイン化が、相互作用するフェルミオン輸送における全デバイス－リザーバー状態の実用的な長時間シミュレーションを可能にすると主張している。リザーバー・モードをその時間的結合プロファイルに基づいて移動窓として再構成することにより、本手法は、実時間シミュレーションを制限するエンタングルメントの継続的な増大を回避する。

著者らは、自身のアプローチが以下の特性を持つことを強調している：

1. 制御された精度: 保持されたモードの占有数に依存しない、全状態に対する厳密な非摂動的誤差上界。
2. 効率性: （固定された閾値において）計算コストが時間ステップに対して線形にスケールし、シミュレーション時間の増大に伴ってリソースが無制限に必要となることはない。
3. 汎用性: 非相互作用および相互作用の両方の領域、および従来の指数関数展開が失敗する複雑なスペクトル密度（例：フラットバンド）を持つ系への適用可能性。

本研究は、さらなる圧縮や伝播技術を適用する前に、物理的なモードを時間的に明示的に分離するという点で、テンソルネットワークや影響汎関数論的アプローチに対する補完的なツールとしての地位を確立している。著者らは、現在の構成は二次的なリザーバーとゼロ温度の初期状態を想定しているものの、本フレームワークは、一般的なガウス状態やマルチターミナル幾何学への将来的な拡張のための基礎を提供するものであると述べている。