$\texttt{Symdyn}$: an automated algebraic solution for high-order quantum systems

技術要約：Symdyn：高次量子系のための自動代数解法

問題提起
多くの重要な量子物理系は、ハミルトニアンが閉じたリー代数の時間独立な生成子の線形結合として $\hat{H}(t) = \sum_{l=1}^L \eta_l(t)\hat{g}_l$ と表されるという特徴を持っています。このような系の時間発展演算子（TEO）は、ウェイ・ニマン法（WNM）を介して、因子分解された形式 $\hat{U}(t) = \prod_{l=1}^L e^{\Lambda_l(t)\hat{g}_l}$ で表現できます。しかし、高次量子系（$L \ge 6$）に対してWNMを適用することは、依然として大きな課題となっています。主な困難は、群の次元（例：$SU(N)$は$N^2-1$ 個の生成子を持つ）が増加するにつれて代数的次数が急速に増大することであり、これにより、手計算での導出が困難な、複雑に結合した非線形微分方程式の増殖を招きます。さらに、指数の順序や基底を変更する場合、これらの複雑な関係を最初から再計算する必要があります。現在、このプロセスを自動化するオープンソースソフトウェアは存在しません。

手法
著者らは、ウェイ・ニマン法の適用を自動化するために設計されたPythonライブラリである symdyn を導入しています。コアとなる手法は、因子分解されたリー群要素の微分を計算するための系統的な代数的アプローチに基づいています。主な技術的構成要素は以下の通りです：

1. 構造テンソルとBCH行列： 本ライブラリは、リー代数の構造定数を含む構造テンソル $\gamma$ を利用します。これは、代数の生成子に関連する横行列（transverse matrices）$\Upsilon_i$ の行列指数関数から導かれるBCH行列（$b_i$）を定義することにより、単一および入れ子状の相似変換（BCH関係）を計算します。
2. 結合行列 ($\xi$)： BCH行列の積を計算することで、ライブラリは結合行列 $\xi$ を構築します。この行列は、TEOの係数の時間微分（$\dot{\Lambda}$）とハミルトニアンの係数（$\eta$）を、$\eta = i \xi^T \dot{\Lambda}$ という方程式を通じて関連付けます。
3. 基底選択（カルタン・ワイル）： 解の全域的な妥当性は基底の選択に依存することを論文では強調しています。本ライブラリは、半単純リー代数に対するカルタン・ワイル基底（CWB）の使用を実装しています。この基底において、横行列は厳密な上/下三角行列（べき零行列）またはブロック対角行列となり、これにより結合行列の逆行列が存在することを保証し（特定の特異点を除いて）、微分方程式のシステムをブロック分離された行列リーカチ方程式の階層へと簡略化します。
4. 自動化： ライブラリは、ユーザー定義の構造テンソルまたは有限行列生成子を受け取り、交換子、BCH行列、結合行列、および結果として得られる微分方程式のシステムを自動的に計算します。また、構造テンソルを再配置することで、指数生成子の順序を変更することもサポートしています。

主な貢献と結果

• symdynの開発： 主要な貢献は、高次量子系に対するWNMの適用を自動化するための初のオープンソースPythonライブラリのリリースです。
• 解析的回収： ライブラリは、低次の代数（$su(1,1)$, $su(2)$, $sl(2)$, $so(2,1)$）における既知の解析的結果を正常に回収し、その実装を検証しました。また、$su(2)$ においてパウリ基底から昇降演算子基底（CWB）へ変更することで、結合行列の特異点が解消され、大域的な解が得られることを示しています。
• 高次系 ($L=11$)： 著者らは、時間依存する結合を持つ2つの時間依存結合調和振動子の系に対してsymdynを適用しました。ライブラリは、11個の結合した非線形微分方程式を正常に導出し、これまで文献では計算されていなかった、ブロック分離されたリーカチ方程式の階層構造を明らかにしました。
• **$SU(N)$への一般化：** ライブラリは、汎用的なカルタン・ワイル基底を提供することで、リー群$SU(N)$に特化しています。これにより、任意の$N$ の処理が可能になります。
  • $SU(2)$： HadamardゲートおよびTゲートを導出するために使用され、ユニバーサル量子計算への有用性を確認しました。
  • $SU(3)$： ライブラリは8つの生成子（ゲルマン行列）を扱い、それらをCWBに変換することで、因子分解されたTEOと対応するブロック分離された微分方程式を取得します。
  • $SU(4)$： フレームワークは、CNOTゲートを構成するために必要な係数を導出するために適用され、ユニバーサル量子計算に必要なゲートのセットを完結させました。

意義と主張
論文は、symdyn が高次元量子系における手計算による導出の限界を克服し、これまで困難であった $L \ge 6$ の系の研究を可能にすると主張しています。相似変換と結果としての微分方程式の計算を自動化することで、本ライブラリは、量子光学、開放系ダイナミクス、量子コンピューティングに関連する複雑な量子系の時間発展の解析を容易にします。

著者らは、全域的な解を得て微分方程式を管理可能な階層へと簡略化するためには、カルタン・ワイル基底の使用が極めて重要であると強調しています。彼らは、彼らの $SU(N)$ に関する結果が、リー代数の次元が増加することに伴う計算の時間発展のスケールアップ問題に対処することに寄与すると主張しています。論文は、symdynを量子ダイナミクスや最適制御の新たな領域を探索するためのツールとして位置付けており、今後は擬似直交群への拡張や、結果として得られる微分方程式をより効率的に解くための「symdyn II」の開発が計画されています。