Lewis-Ermakov approach for the time-dependent two-level system

技術要約：時間依存型二準位系におけるルイス–エルマコフ・アプローチ

問題提起
時間依存型二準位系（TLS）は量子制御の基本であるが、通常、時間依存型調和振動子（THO）とは別個のケースとして分析されてきた。ルイス–エルマコフ法は、THOに対しては、動的不変量、閉形式の進化演算子、および任意の時間依存パラメータに対する厳密なプロパゲータ（伝播関数）を導出するという堅牢な枠組みを提供しているが、TLSに対する統一的かつ十分に発展したルイス–エルマコフの枠組みは欠如していた。具体的には、任意の結合 $a(t)$ およびデチューニング $b(t)$ に対して有効な、動的不変量、因子分解された時間依存ハミルトニアン、および閉形式のプロパゲータを備えた定式化が、THOで使用されるものと同じ統一的な言語を用いて確立されていなかった。このギャップは、特に、例外点（exceptional points）や状態依存の利得が断熱追従の概念を複雑にする、開放系または実質的に非エルミート的な設定において顕著である。

手法
著者らは、TLSとTHOの間の代数的対応関係を利用することで、TLSのための明示的なルイス–エルマコフ型の動的不変量を構築する。この手法は、両方のシステムが同一の複素化されたリー代数 $sl(2, \mathbb{C})$ の表現を実現しているという事実に依拠している：

• THOは、$\{\hat{x}^2, \hat{p}^2, \hat{x}\hat{p} + \hat{p}\hat{x}\}$ によって生成される非コンパクトな $su(1, 1)$ 代数を利用する。
• TLSは、$\{\hat{\sigma}_-, \hat{\sigma}_+, \hat{\sigma}_z\}$ によって生成されるコンパクトな $su(2)$ 代数を利用する。

コアとなる手順は、時間依存ハミルトニアン $\hat{H}_{TLS}(t) = a(t)(\hat{\sigma}_+ + \hat{\sigma}_-) + b(t)\hat{\sigma}_z$ を、一対の時間依存変換 $\hat{G}(t)$ および $\hat{T}(t)$ を用いて定常形式へと写像することである。

1. 変換： 変換 $\hat{G}(t)$ はデチューニングとスケーリングを組み込み、$\hat{T}(t)$ は残差項を除去する。注目すべきは、$\hat{\sigma}_-$ の冪零性により $\hat{G}(t)$ が非ユニタリであるため、中間的な複素化された状態空間が必要となる点である。
2. エルマコフ–ピニー方程式： この写像には、補助的なスケーリング関数 $\rho(t)$ が以下の複素エルマコフ–ピニー方程式を満たす必要がある：
   $$\ddot{\rho} + \theta^2(t)\rho = \frac{\theta_0^2}{\rho^3}$$
   ここで、有効周波数 $\theta^2(t)$ は物理パラメータ $a(t), b(t)$ およびそれらの時間微分を符号化している：
   $$\theta^2(t) = a^2 + b^2 + i\dot{b} - i\frac{\dot{a}b}{a} + \frac{\ddot{a}}{2a} - \frac{3\dot{a}^2}{4a^2}$$
3. 不変量の構築： 動的不変量 $\hat{I}_{TLS}(t)$ は、定常演算子 $(\hat{\sigma}_+ + \hat{\sigma}_-)$ に対して複合変換 $\hat{S}_{TLS}(t) = \hat{G}(t)\hat{T}(t)$ を作用させることで構築される。この不変量は、ルイスの条件 $\frac{d\hat{I}}{dt} = \frac{\partial \hat{I}}{\partial t} - i[\hat{I}, \hat{H}] = 0$ を満たす。
4. 進化演算子： 変換されたハミルトニアンが常に自身と可換であるため、時間順序演算子は自明となる。厳密なプロパゲータは、スケーリング関数 $\rho(t)$ と累積位相 $f(t) = \int dt'/\rho^2(t')$ を用いて、閉形式で導出される。

主な貢献と結果
本論文は、以下の具体的な結果を提示している：

• 一般解： $\rho(t)$ に対する複素エルマコフ–ピニー方程式が解を持つ限り、任意の $a(t)$ および $b(t)$ に対して、厳密な閉形式の進化演算子およびプロパゲータを提供する。
• ランダウ–ゼネー–マヨラナ–シュテュッケルベルク（LZMS）遷移： 本フレームワークを標準的なLZMSモデル（$a(t)=g, b(t)=\alpha t$）に適用している。この手法は、既知の非断熱ダイナミクスと遷移確率を回収する。その際、複素周波数 $\theta^2(t)$ が、虚数成分 $i\alpha$ を通じてスイープレートの効果を捉える。
• 非エルミート力学： 虚数デチューニング $b(t) = i\alpha t$ を持つ散逸的LZMSモデルへと拡張している。著者らは、静的な例外点（固有値が合体する点）は瞬時スペクトルによって決定される一方で、「動的な」例外点（指数関数的な増幅が始まる点）はシフトすることを実証している。このシフトは、有効周波数 $\theta^2(t)$ のゼロ点によって支配されており、有限のスイープレートが対称性の破れの開始点を変化させることを浮き彫りにしている。
• 断熱急速通過（ARP）： アレン–エベリーモデル（双曲線正弦結合および双曲線正接デチューニング）に適用している。不変量に基づくアプローチは、人口移動（ポピュレーション・トランスファー）の近似なしの厳密な記述を提供し、有限のスイープレートによる完全な断熱性からの逸脱を正確に捉えている。
• 逆設計（断熱への近道）： 残留励起（量子摩擦）なしに、指定された軌道に沿ってシステムを駆動する制御場 $a(t)$ と $b(t)$ を設計するためのプロトコルを開発している。
  • 積分定数 $\theta_0$ を自由な幾何学的パラメータ $\vartheta_0$ に格上げすることで、不変量がプロトコルの開始時と終了時にハミルトニアンと可換であることを保証する $\rho(t)$ の境界条件を導出している。
  • 実数のスケーリング関数 $\rho(t)$ については、逆設計によって完全に明示的な閉形式の制御場 $a(t) = \vartheta_0 \cos\phi / \rho^2(t)$ および $b(t) = \vartheta_0 \sin\phi / \rho^2(t)$ （ここで $\phi$ は固定角）が得られる。これは、結合とデチューニングの比例的な変調に限定されるが、解析的な解法可能性を保証する。
  • 複素の $\rho(t)$ については、制御場は結合された非線形積分微分方程式によって決定される。これは、$a(t)$ と $b(t)$ を独立に成形するためのより大きな柔軟性を提供するが、数値解を必要とする。

意義と主張
本論文は、その主要な貢献は、時間依存型二準位系に対するルイス–エルマコフ動的不変量の明示的な構築とその実践的な応用にあると主張している。その意義は以下の点にある：

1. 統一： 連続変数（振動子）と離散変数（二準位系）の $sl(2, \mathbb{C})$ 代数の実現の間の概念的および数学的リンクを確立し、因子分解、補助方程式、および不変量構築という同一のプロセスが両方に適用されることを示している。
2. 計算上の有用性： 問題固有の積分を単一の補助方程式（複素エルマコフ–ピニー方程式）に置き換え、非ユニタリな領域を含む任意の駆動場に対して厳密なプロパゲータを提供する。
3. 制御設計： 不変量を理論的なツールから、逆設計のための実践的なリソースへと昇華させている。これは、「断熱への近道」を実現するための制御場を設計するための体系的な方法を提供しており（実数 $\rho$ の場合は完全な明示的制御場、複素 $\rho$ の場合は数値的に決定される場）、駆動場に対する恣意的な制約を排除している。
4. 汎用性： 著者らは、代数的フレームワークが任意の $N$ レベル系（スピン $j$ 表現）へ自然に拡張可能であることに触れている。そこでのダイナミクスもまた、係数が修正されたエルマコフ–ピニー方程式によって支配されるため、マルチレベルの量子ダイナミクスに対する普遍的なアプローチを示唆している。

本研究は新しい実験セットアップを提案するものではなく、むしろ既存のシナリオ（捕捉イオン、超伝導回路など）を分析し、それらのための制御プロトコルを設計するための理論的メカニズムを提供するものである。