タグ: リウビリアンのKramers縮退

環境とエネルギーや情報をやり取りする量子系では、振動が弱まり、状態が時間とともに変わります。それでも、異なる二つの変化の仕方が、同じ減衰率と同じ振動周波数を持つことはあるのでしょうか。

東京大学物性研究所のKohei Kawabataによる2026年11月2日公開のプレプリントは、開放量子系の時間変化を記述するリウビリアンに、Kramers（クラマース）縮退が生じる条件を示しています。鍵になるのは、一般化した時間反転対称性とパリティ対称性の組み合わせです。

簡単にまとめると？

• リウビリアンの固有値は、量子系の変化のモードがどの速さで減衰し、振動するかを表します。異なるモードが同じ固有値を持つことを、ここでは縮退と呼びます。
• 論文は、一般化した時間反転を二回行う操作がパリティになる場合、その符号が負の演算子部分空間で、モードが対になることを示しています。
• 対になるのは二つの保存可能な量子状態とは限りません。減衰がなくなる保証や、量子メモリーの実機性能を示した研究ではありません。

エネルギーの縮退から、変化の仕方の縮退へ

孤立した量子系のエネルギーを表す基本的な道具が、Hamiltonian（ハミルトニアン）です。異なる状態が同じエネルギーを持つとき、エネルギー準位が縮退しているといいます。通常のクラマースの定理は、半整数スピンを持つ系などで、時間反転対称性によって状態が対になることを説明します。

環境と相互作用する開放量子系では、エネルギーだけを見ても時間変化を記述しきれません。確率的な混合も扱う密度行列に対し、時間変化の生成子として働くのがリウビリアンです。密度行列などの演算子に作用するため、スーパーオペレーターとも呼ばれます。

リウビリアンの固有モードXが固有値λを持つなら、そのモードは時間に対してeのλt乗で変化します。λが複素数の場合、実部が減衰の速さ、虚部が振動の角周波数に対応します。例えばλ＝−γ＋iωなら、振動の振幅はeの−γt乗で小さくなり、角周波数はωです。

この意味での縮退は、異なる独立なモードに同じλが対応することです。同じγとωを持っていても、モードが表す演算子や空間的な変化の形は異なり得ます。また、γが正なら両方とも減衰します。固有値が等しいことと、情報が失われないことは別の性質です。

時間反転を二回行ったときの符号が鍵になる

クラマースの定理で重要なのは、時間反転を表す反ユニタリ演算子Tを二回作用させたとき、Tの二乗が−1になることです。対称性が成り立つと、元の状態に対応する別の独立な状態を作れ、同じ固有値を持つ対が必要になります。

ただし、リウビリアンのような非エルミート演算子では、通常の時間反転と、論文が「時間反転対称性†」と呼ぶ一般化された条件を区別しなければなりません。後者は、リウビリアンLそのものではなく、その随伴L†も使う条件、TL†Tの逆操作＝Lとして定義されています。エルミートな場合には両者は一致しますが、一般には同じ条件ではありません。

「時間反転」という語だけから、散逸する運動を物理的に逆再生できると読んではいけません。ここで使うのは、時間変化を記述する演算子が満たす数学的な関係です。論文のスピン鎖の例でも、通常のスピンの時間反転だけではなく、それに一様な回転を組み合わせた有効な操作を使っています。

全体で同じ符号でなくても、部分空間では負になる

密度行列には左側と右側から演算子が作用します。この二重の構造のため、従来の議論では、対称性ごとに分けたリウビリアンの内部にはクラマース縮退が現れないとされていました。今回の論文は、その議論で置かれていた対称性の表現に関する仮定を見直しています。

用いる構造は、一般化した時間反転を二回作用させた結果が、全体で一律の＋1や−1ではなく、パリティ演算子になるというものです。パリティは、二回作用させると元に戻り、対象を符号＋1と−1の二つの部分に分ける対称性です。粒子数の偶奇に対応する場合もありますが、今回注目するのは演算子が属する部分空間です。

リウビリアンがこのパリティを保存するなら、＋1側と−1側を別々に調べられます。−1側では時間反転の二乗が−1となり、一般化したクラマースの条件が成り立ちます。全体を一つの符号で扱うだけでは見えなかった対が、個々の部分空間の内部に現れる、という仕組みです。

どんなモデルで確かめたのか

論文は、相互作用するフェルミ粒子のモデルと、散逸を伴うスピン鎖を具体的に構成しています。環境との結び付きを表すジャンプ演算子にも、時間反転とエルミート共役に関する条件を課しています。系のハミルトニアンだけに対称性があれば、どのような環境との結合でも同じ結論になる、という話ではありません。

フェルミ粒子の例では、相互作用と位相緩和を含むリウビリアンの複素固有値を計算しています。スピン鎖の数値例は六つのスピンからなり、相互作用の乱れを変えた5000個の実現について、固有値の間隔の統計を調べています。これは理論モデルの数値計算で、六量子ビットの装置による実測ではありません。

間隔の比較にも手順があります。ある固有値から最も近い固有値への距離と、その次に近い固有値への距離を、複素数の比として取ります。縮退した対は同じ位置にあるため、論文は対を一つとして数えたうえで、この統計を計算しています。対を二回数えると、最近接距離がゼロになり、別の比較になってしまいます。

結果は、対応する対称性を持つ非エルミート・ランダム行列の統計とよく一致すると報告されています。ただし、比較するスペクトルの領域を選んでおり、スピン鎖では比の絶対値が1に近いところのずれも記されています。数値の一致は、有限のモデルと選んだ条件の下で評価されたものです。

二つの密度行列を保存できる、という意味ではない

特に注意したいのは、負のパリティ部分空間に、物理的な密度行列そのものは入らない点です。例えばパリティ変換でXが−Xになるなら、その演算子の跡はゼロになります。一方、正規化した密度行列の跡は1です。この違いから、対になった二つのモードを、そのまま二つの保存可能な量子状態と同一視することはできません。

それでも、モードが物理的に無意味になるわけではありません。論文は、負のパリティを持つ演算子の二時刻相関、つまり異なる時刻の量の結び付きから、独立な動的モードを調べられると述べています。状態そのものと、その時間変化や相関を記述する成分を分けると、何を測ればよいのかが見えてきます。

この研究の意義は、散逸する量子系でも、対称性の組み合わせによって縮退するモードが存在し、そのスペクトル統計を分類できると示したことです。応用を考える際には、実際の環境が必要な対称性を保つか、固有値の実部はどれだけ負か、どの相関でモードを読み出せるかを別々に確認する必要があります。縮退の存在を、長寿命な量子メモリーの完成へ直接置き換えることはできません。

出典

• Kohei Kawabata「Liouvillian Kramers Degeneracy」、arXiv:2611.03881v1、2026年11月2日公開。プレプリント。
• 同論文のPDF。一般化した時間反転と部分空間：1〜2ページ、フェルミ粒子とスピン鎖のモデル・数値統計：3〜4ページ、密度行列と二時刻相関の区別：4ページ。7〜8ページには追加のモデルが示されています。