Real-axis least-squares bath discretization for real-time transport in quantum-dot Josephson junctions

技術要約：量子ドット・ジョセフソン接合におけるリアルタイム輸送のための実軸最小二乗浴離散化

問題提起
駆動される超伝導ナノ構造の効率的なリアルタイム・シミュレーションは、超伝導リード（電極）の表現によって阻害されている。そこには、以下の3つの相反する要件によるボトルネックが存在する：

1. 回帰性（Recurrence）： 有限の離散的な浴はポアンカレ回帰（Poincaré recurrences）に苦しみ、これが遅い観測量に影響を与える（$t_{rec} \approx \pi N_b/D$）。これを抑制するには、モード数 $N_b$ をシミュレーション時間 $T$ に対して線形に増加させる必要があり、密な相関行列伝播において計算コストが $T^4$ に近いスケールで増大する。
2. スペクトル分解能： BCSギャップ端の特異性とサブギャップ・コヒーレンス（アンドレエフ物理において重要）は、標準的な離散化が最も困難な領域において高いスペクトル分解能を要求する。
3. 減衰のトレードオフ： グローバルな減衰は回帰を抑制するが、ギャップ端を広げ、人工的にコヒーレントな寿命を短縮させてしまう。単一のグローバルな減衰率では、長時間の忠実度と定常状態の精度を同時に維持することはできない。既存の手法（Floquet、NRG、固定されたDynes広がりを持つNEGF、HEOM）は、過渡現象に対して失敗するか、計算コストが長時間のダイナミクスに対して高すぎるか、あるいは広がりスケール以下の減衰率を覆い隠してしまう分解能の底（resolution floor）を導入してしまう。

手法
著者らは、駆動型リウヴィル・フォン・ノイマン（DLvN）アプローチを超伝導リザーバー（Bogoliubov–de Gennesリード）へと拡張している。核心となる革新は、浴の離散化を、あらかじめ規定されたグリッドとしてではなく、制約付き最適化問題として扱うことである。

• 定式化： 系は、BCSリードに結合した非相互作用量子ドットとしてモデル化される。DLvN方程式はナンブ（Nambu）基底へと適応されており、これにより、各モードの粒子およびホールの成分に対して同じ減衰 $\gamma_j$ を割り当てることで、減衰構造が粒子・ホール対称性を保持することを保証している。
• ゲージの規律（Gauge Discipline）： 静的な参照状態 $C_{ref}$ との一貫性を維持するため、時間依存のバイアス/駆動はホッピング位相（Peierls）ゲージを介して導入される。これにより、減衰のあるサイト上で回転項が不整合を起こすことを防ぎ、浴の軌道を固定したままにすることができる。
• 実軸最小二乗（LS）フレームワーク： 連続的なリードのハイブリダイゼーション自己エネルギーは、$N_b$ 個の減衰したBdGモードの和として近似される。モードエネルギー ($\varepsilon_j$)、スペクトル重み ($w_j$)、およびモード分解された減衰 ($\gamma_j$) は、ターゲットとなる自己エネルギーに対する加重最小二乗誤差を最小化することによって、共同で決定される。
• 4つの設計制約：
  1. ギャップ・クリーンな重み付け： 寄生的オーム抵抗による短絡を防ぎ、アンドレエフ束縛状態（ABS）を保護するために、ギャップ内部のスペクトル重みに対して高いペナルティを課す。
  2. モードごとの $\gamma$： バンド内のモードには大きな減衰（吸収）を与え、ギャップ端のモードには小さな減衰（保護）を与える。
  3. 密度底（Density floor）： ギャップ端付近に高密度で等間隔なコアを配置することで、離散モードが主信号よりも長く生き残ることによって生じる「エッジ・ビート」のアーティファクトを防ぐ。
  4. 常伝導リードの制約： 常伝導リードの場合、フェルミ占有数を正しく固定するために輸送窓内に高密度なコアを保持し、ウィング部分にのみ圧縮を適用する。

主な貢献

1. 明示的なハミルトニアン・フレームワーク： 本手法は、最適化されたパラメータを持つ浴の明示的なハミルトニアン表現を構築し、マルチターミナル常伝導・超伝導幾何学へ自然に一般化するとともに、相互作用のあるインパリティ・ソルバーへの道筋を提供する。
2. 広がり（Broadening）の独立した較正： 減衰パラメータ $\gamma$ は物理的なパラメータとしてではなく、独立して較正される数値的な広がりとして扱われる。その長時間の主要な影響は、パラメータフリーのスケール則 $\kappa \simeq w_{bath}\gamma_{eff}$ に従うトラップされた準粒子の崩壊である。
3. 共同精緻化パス： 本フレームワークは、制御可能な精緻化パス $(\gamma \to 0, N_b \to \infty)$ を定義しており、従来の有限広がり計算の分解能スケル以下の減衰率を数値的に解決することを可能にする。
4. 計算効率： 浴を最適化することにより、精度を高い解像度のリファレンスに対して維持しながら、計算コストを1桁以上削減することに成功した（例：7.7時間のシミュレーションを20分に短縮）。

結果と検証
本フレームワークは、以下の4つのシナリオにおいて、確立された結果および高解像度のリファレンスを用いて検証された：

1. 静的な電流・位相関係（CPR）： 本手法は、連続体自由エネルギーCPRのベンチマークを再現する。静的なCPRの収束は、虚軸を用いたバリアントよりも実軸LSの方が遅いが、ダイナミクスには十分である。
2. 多重アンドレエフ反射（MAR）： 本手法は、サブギャップの階段状構造、パリティ効果（奇数次での増強）、および分散関係を再現する。これは、110個のモードを用いた高解像度ガウス・リファレンスに対し、わずか32個のモードで（相対偏差3〜4%）一致している。
3. 準粒子捕捉ダイナミクス： 本手法は、バイアス・クエンチ後の捕捉された準粒子の振動周波数を再現する。決定的なことに、本手法はスケール則 $\kappa \simeq w_{bath}\gamma_{eff}$ を検証し、従来の有限広がり計算の分解能スケル以下の減衰率 ($\kappa \approx 5.7 \times 10^{-4}$) を解決できることを示した。
4. シャピロ・ロッキング： 本手法は、整数および分数シャピロ・ステップを捉える。また、選択則（例：$m=2$ ロックにおける $\pi$ 周期性）を裏付け、周波数領域のFloquet法が計算上困難な領域（無理数バイアス点や高分母の有理数）においても有効であることを示した。

意義と境界
本論文は、本フレームワークが、回帰抑制と寿命の忠実度の間のトレードオフを克服し、超伝導接合におけるリアルタイム輸送のための堅牢でスケーラブルな経路を提供することを主張している。これにより、広範な時間およびパラメータ範囲にわたる過渡的なコヒーレンスと非線形輸送の研究が可能になる。

明示された制限事項：

• 相互作用： 現在の研究は、厳密に非相互作用ドット（$U=0$）に限定されている。
• 有限ギャップ・強駆動振幅： シャピロ・ステップの位置と選択則は堅牢であるが、有限ギャップ領域における強駆動の絶対振幅は、周波数領域のFloquet定常状態から $O(1)$ レベルで逸脱する。これらの振幅は定性的であると報告されている。
• 分解能要件： 極端に小さい $\kappa$ の領域では、必要なモード密度は $1/\gamma$ に比例してスケールする。係数の圧縮は最適化を行うが、この根本的な分解能要件を回避することはできない。
• 高周波観測量： 長時間分光のために使用される特定の「密度底」を持つ浴は、遠方の裾（far-tail）のアンカーが存在するため、高周波スペクトルモーメントには適していない。

著者らは、本手法が確立された物理を再現しつつ、パラメータマッピングや過渡ダイナミクスにおける新たな能力を可能にすることに成功したと結論付けており、その定量的精度の境界についても今後の発展のために明示的にマッピングされている。