TensorMixedStates: a Julia library for simulating pure and mixed quantum states using matrix product states

本論文は、行列積状態表現を用いて、純粋状態および混合状態の量子状態（リンドブラッド方程式および非ユニタリゲートによる散逸ダイナミクスを含む）の効率的なシミュレーションを可能にする、ITensor に基づいたユーザーフレンドリーな Julia ライブラリ「TensorMixedStates」を紹介する。

要約

この論文を、平易な言葉と日常的な比喩を用いて解説します。

全体像：漏れのある量子世界のシミュレーション

複雑な機械がどのように機能するかを予測しようとしていると想像してください。量子物理学の世界において、この機械は原子や電子のような微小な粒子で構成されています。通常、科学者たちは、これらの粒子を、何も出入りしない完全な密閉箱の中にあるかのようにシミュレーションしようとします。これを「純粋状態」と呼びます。

しかし、現実世界では、完全に密閉されたものなど存在しません。これらの量子機械は常に環境と衝突し、エネルギーを漏らしたり、「ノイズ」を受けたりしています。これを「混合状態」と呼びます。漏れがありノイズの混じった系をシミュレーションすることは、計算機にとって極めて困難です。なぜなら、数学が複雑になり、非常に急速に計算量が爆発的に増大してしまうからです。

TensorMixedStates は、Julia 言語で書かれた新しいコンピュータプログラム（ライブラリ）であり、科学者たちがこれらの「漏れのある」量子系をシミュレーションするのを助けます。これは、ノイズ、熱、あるいは散逸によって乱されたときに、量子状態がどのように変化するかを追跡できる、専門的なツールキットのような役割を果たします。

中核となるツール：「MPS」というバックパック

このライブラリがどのように機能するかを理解するには、「行列積状態（MPS）」という概念を理解する必要があります。

長い鎖のように人々が手をつながっている様子を想像してください。この鎖全体を記述しようとする場合、すべての人の正確な位置を一度に書き下そうとすると、鎖が長ければ長いほど、そのリストは信じられないほど巨大なものになってしまいます。

その代わりに、MPS 法はこう提案します。「各人が隣接する隣人とどのように手をつないでいるかだけを記述しましょう」。大きな問題を、小さな局所的なつながりに分解することで、情報を圧縮できるのです。これは、長い物語を記述する際に、本全体を毎回書き直すのではなく、登場人物の各ペアの関係性を要約して記述するようなものです。

TensorMixedStates ライブラリは、この「バックパック」方式を採用し、それをアップグレードしました。

• これらのツールの「旧バージョン」は、「純粋状態」（完全な密閉箱）しか運ぶことができませんでした。
• TensorMixedStates は、「混合状態」（漏れがありノイズの混じった箱）を運ぶことができます。それは、複雑で漏れのある情報を、依然として効率的に圧縮・管理できる特殊なベクトルとして扱います。

仕組み：「レゴ」アプローチ

この論文では、このライブラリが ITensor というもう一つの有名なツールの上に構築されていると説明しています。ITensor を、非常にうまく組み合わさる高品質なレゴブロックのセットだと考えてください。

• 問題点: 元のレゴセット（ITensor）は、完璧で剛体な構造（純粋状態）を構築するように設計されていました。ぐらついたり溶けたりする構造（混合状態）には、適切な接続部品が備わっていなかったのです。
• 解決策: 著者たちは、レゴセットの上に載る新しい「アダプターキット」（TensorMixedStates）を構築しました。このキットにより、同じ強力なレゴブロックを土台としながら、ぐらついたり溶けたりする構造を構築できるようになります。

このライブラリは、主に 3 つのスーパーパワーを提供します。

1. 混乱の処理: 純粋状態に使用される効率的な「バックパック」（MPS）方式を用いて、混合状態の数学である密度行列を表現できます。
2. 時間旅行: これらの系が時間とともにどのように変化するかをシミュレートできます。これには以下が含まれます。
   • シュレーディンガー進化: 完全に孤立しているときに系がどのように変化するか。
   • リンドブラッド進化: エネルギーが漏れたり、ノイズの混じった環境と相互作用したりするときに系がどのように変化するか。
   • 量子チャネル: 誤差を導入する可能性のある特定の「ゲート」や操作を適用したときに系がどのように変化するか（ノイズの混じった量子コンピュータなど）。
3. ユーザーフレンドリーなインターフェース: 著者たちは「高レベル」のインターフェースを構築しました。これは、科学者が複雑なシミュレーションを、何千行もの生々しい数式コードを書くのではなく、まるでレシピを書くように数行のコードだけで記述できることを意味します。

論文内の実世界例

この論文は理論について語るだけでなく、6 つの異なる物理シナリオにおいてライブラリが実際に機能していることを示しています。以下に、彼らがテストした内容の簡単な内訳を示します。

• ノイズのあるフェルミオン鎖: ワイヤーを伝ってホッピングする電子の列を想像してください。研究者たちは、電子がどのように広がり広がるかを見るために、「位相ノイズ」（ラジオの雑音のようなもの）を追加しました。ライブラリの結果は、正確な数学的答えと完全に一致しました。
• 漏れのあるスピン鎖: 小さな磁石（スピン）の列を想像してください。この列の両端は、磁石を反転させようとする「リザーバー（熱浴）」に接続されています。ライブラリは、磁気性が鎖全体をどのように流れるかを正常にシミュレートしました。
• ボソン源: 空のスポットの列に粒子を注入するパイプを想像してください。ライブラリは、粒子が局所的なスペースに制限されている場合でも、時間とともに粒子が列をどのように満たしていくかを追跡しました。
• グラフ状態のデコヒーレンス: 絡み合った量子ビット（量子ビット）の複雑な網を想像してください。研究者たちは、この網がノイズにさらされたときにどのようにほどけ（デコヒーレンス）、崩れていくかを観察しました。このライブラリは、このタイプの計算としては非常に巨大な 512 量子ビットの系について、このシミュレーションを実行することができました。
• ノイズのある回路: ゲート（スイッチ）が時々誤りを犯す量子コンピュータの回路を想像してください。ライブラリは、ゲートと誤りの「レンガ壁」パターンをシミュレートし、ノイズによって「もつれ」（部分間の量子接続）がどのように成長し、その後破壊されるかを示しました。

なぜこれが重要なのか（論文によると）

この論文は、このライブラリが欠落部分を埋めていると主張しています。完璧な量子系をシミュレートするための優れたツールは存在するものの、現実的でノイズの混じった系をシミュレートするためのツールは不足していたか、使用が困難でした。

• 効率性: 計算を高速かつ正確に保つために、TDVP や DMRG などの最良のアルゴリズムを使用しています。
• 精度: シミュレーションが「いい加減」になっているかどうか（例えば、数学が物理的な現実から逸れ始めている場合など）をユーザーに知らせるための組み込みチェック機能が含まれています。
• アクセシビリティ: 研究者が数行のコードで高度なシミュレーションを設定できるようにし、現実のノイズの混じった世界における量子系の挙動を研究しやすくしています。

要約すると、TensorMixedStates は、科学者たちが理論の滑らかで完璧な道路だけでなく、現実世界の荒々しくノイズの混じった地形を、量子シミュレーションの車で行き来できるようにする、新しくユーザーフレンドリーなエンジンです。

技術概要

「TensorMixedStates: 行列積状態を用いた純粋状態および混合量子状態のシミュレーションのための Julia ライブラリ」の技術的概要

問題提起
開放量子多体系のシミュレーションは、冷たい原子、トラップイオン、超伝導回路における実験の進展に駆動され、現代物理学における重要な課題です。閉じた系とは異なり、開放系は環境と相互作用し、ノイズ、散逸、およびデコヒーレンスを引き起こします。これらは通常、リンドブラッド方程式によって記述されます。テンソルネットワーク手法、特に行列積状態（MPS）は、純粋状態のシミュレーション（ITensor や TenPy などのライブラリを使用）の標準となっていますが、MPS 枠組み内で混合状態（密度行列）をシミュレートするためのソフトウェアは、はるかに限られています。既存の解決策は、柔軟性に欠け、特定の局所ヒルベルト空間（例えば、量子ビットのみ）に制限されているか、最先端のアルゴリズムとシームレスに統合されていないことが多くあります。さらに、混合状態の表現には、密度行列を倍増されたヒルベルト空間内のベクトルとして扱う必要があり、演算子の作用（例えば、$O\rho$、$O\rho O^\dagger$、交換子）における複雑さが生じます。これは、標準的な純粋状態の枠組みではネイティブにサポートされていません。

手法
著者らは、ITensor ライブラリを基盤として構築された TensorMixedStates (TMS) という Julia ライブラリを導入します。TMS は、MPS 表現を用いて純粋状態および混合量子状態の両方をシミュレートするように設計されています。

• コアアーキテクチャ: TMS は、ITensor の堅牢な低レベルのテンソル操作と高性能アルゴリズム（DMRG、TDVP）を活用します。しかし、ITensor のネイティブな枠組みが純粋状態に最適化されていることを認識し、TMS は混合状態に対応するために、新しい状態および演算子枠組みを一から実装しています。
• 混合状態の表現: このライブラリは、密度行列 $\rho$ を、次元 $d^2$ のより大きな空間におけるベクトル化された状態 $|\rho\rangle\rangle$ として表現します。これにより、密度行列を 4 指標テンソルの積（またはベクトル化された基底では 3 指標テンソル）に分解する MPS 手法を利用することが可能になります。
• 演算子枠組み: 重要な設計上の選択として、演算子が密度行列に作用する 4 つの異なる方法（$O\rho$、$O\rho O^\dagger$、$[O, \rho]$、および $\{O^\dagger O, \rho\}$）を処理できる柔軟な演算子枠組みの開発が含まれます。これには、リンドブラッド散逸項および一般的な量子チャネル（クラウス演算子）のサポートが含まれます。
• アルゴリズム:
  • 時間発展: TMS は、リンドブラッド方程式を介した連続的な時間発展を、2 つの主要なアルゴリズムである TDVP（時間依存変分原理）および WI/WII 近似（指数写像の MPO 近似）を使用してサポートします。WI/WII 法は、特に非エルミートな時間発展において TDVP を補完するものとして注目されています。
  • 定常状態: このライブラリは、長時間の積分を必要とせずに、エルミート演算子 $L^\dagger L$ の基底状態を DMRG によって見つけることで定常状態を計算します。
  • 離散時間発展: 非ユニタリゲートおよび量子チャネルの適用をサポートしています。
• インターフェース: TMS は、微細な制御を可能にする低レベル API と、数行のコードで複雑なシミュレーションワークフロー（状態の作成、時間発展、測定）を定義できる高レベルインターフェース（runTMS）という、二重のインターフェースを提供します。

主な貢献

1. 一般化された混合状態シミュレーション: TMS は、量子ビットを超えて、ボソン、フェルミオン、スピン、および $t-J$ モデルを含む一般的な局所ヒルベルト空間へと MPS の能力を拡張し、開放量子系のソフトウェア生態系におけるギャップを埋めています。
2. 柔軟な演算子枠組み: ハミルトニアン、リンドブラッド散逸項、および一般的な量子チャネルを構築するための包括的なツールセットを提供し、混合状態操作の代数的複雑さを処理します。
3. アルゴリズムスイート: このライブラリは、時間発展のための TDVP および WI/WII 近似、およびリンドブラッド形式に特化して適応された定常状態計算のための DMRG を統合しています。
4. ユーザーフレンドリーな設計: 高レベルインターフェースは、洗練されたシミュレーションのセットアップを簡素化し、基盤となる Julia/ITensor アーキテクチャはマルチスレッド処理と効率的なテンソル縮約を確保します。
5. 検証: このライブラリは、いくつかの準自由モデル（位相ずれまたは注入を伴うフェルミオンおよびボソン鎖、境界散逸を伴う XX スピン鎖）の厳密な解析的解、およびグラフ状態のデコヒーレンスおよびノイズのある量子回路に関する既知の結果に対して検証されています。

結果
論文は、6 つの例を通じてライブラリの能力を実証しています。

• 位相ずれを伴うフェルミオン鎖: フェルミオン密度および演算子空間エンタングルメントエントロピー（OSEE）のシミュレーションは、さまざまな散逸強度に対して厳密な漸近結果と優れた一致を示します。
• 境界散逸を伴う XX スピン鎖: 磁化およびスピン電流の時間発展は、以前の文献からの結果と一致します。
• 局所ソースを伴う自由ボソン: このライブラリは、粒子注入を伴うボソン鎖のダイナミクスを成功裡にモデル化し、密度プロファイルおよび全粒子数の増加を捉えますが、計算コストは局所ヒルベルト空間の次元とともに増加します。
• 局所ソースを伴う自由フェルミオン: TDVP と WI/WII アルゴリズムの比較は、このモデルに対して両者が同様の精度を提供することを示していますが、TDVP ははるかに多くの CPU 時間を必要とします。
• 完全グラフ状態のデコヒーレンス: 状態の低いエンタングルメントにより小さな結合次元で済むため、このライブラリは最大 $N=512$ の大きなシステムサイズを効率的に処理します。
• ノイズのある量子回路: このライブラリは、ユニタリゲートおよび脱分極ノイズを伴うレンガ壁型量子回路をシミュレートし、ノイズによる OSEE の増加およびその後の減少を追跡します。

意義と主張
著者らは、MPS を用いた開放量子系のシミュレーションにおけるソフトウェアのギャップを埋めるために TMS を必要なツールとして位置づけています。彼らは、このライブラリが ITensor を介して「最先端のアルゴリズム」へのアクセスを維持しながら、「数行のコード」で洗練されたシミュレーションを可能にする「ユーザーフレンドリーなインターフェース」を提供すると主張しています。論文は、このライブラリが軌道アンラベリングを実装するのではなく、密度行列の直接時間発展に焦点を当てていることを強調しています。その意義は、先行する Lindbladmpo とは異なり、一般的な局所ヒルベルト空間を処理できる能力と、現代の Julia 生態系との統合にあります。著者らは控えめに、ライブラリは強力であるものの、局所次元が大きい系（ボソンのような系）の計算コストは依然として高く、今後の研究では効率をさらに向上させるために保存量子数（QNS）の導入と自動 MPO 圧縮に焦点を当てることを述べています。