量子系の時間発展を逆にたどれるとしても、その一部分を観測した結果に、最初の状態の違いがいつまでも現れるとは限りません。この二つが両立する条件を、無限に並ぶ振動子の模型で調べた研究があります。
UCLAのRoman Geiko氏とJake Gerenraich氏による「Thermalization with Gaussian Quantum Cellular Automata」は、ガウス量子セル・オートマトン（GQCA）と呼ぶ更新規則について、特定の局所観測量が無限温度状態の値へ近づく十分条件を示しています。初稿は2026年6月4日、参照するv2は11月3日公開の小改訂です。実機の温度測定ではなく、観測量と初期状態の範囲を明確にした数学的な研究です。
簡単にまとめると？
- 対象は、各格子点に一つのボソンモードを置いた無限格子。有限の範囲にだけ作用するWeyl演算子の期待値で、熱化を捉えます。
- 影響する格子点の数が増えても、観測量を指定する係数が十分速く伸びれば、非自明なWeyl演算子の期待値はゼロへ近づきます。
- 「伸びる方向がある」だけでは不十分です。局所的な観測量が縮む方向だけに残らないことなど、更新規則にも条件が要ります。
セルを更新するのは、二値のビットではなく振動子
セル・オートマトンは、近くのセルに関係する規則を繰り返し適用して、全体を動かす模型です。この論文ではセルを整数格子上に置き、各点に位置と運動量を持つ一つのボソンモードを対応させます。二状態だけの量子ビットとは異なり、一つのセルの状態空間も無限次元です。
ガウス量子セル・オートマトンでは、位置・運動量に対応する位相空間の変換を使って、観測量の時間発展を記述できます。同じ規則を格子のどこでも用い、1ステップで影響が届く範囲は有限。逆向きの更新にも同じ有限範囲の条件を課します。「ガウス」はこの変換の構造を指し、初期状態までガウス状態に限るという意味ではありません。
何を見たときに「熱化した」と言えるのか
論文が追うのは、有限個のセルの位置・運動量を係数付きでまとめ、それを指数の中に入れたWeyl演算子です。係数の組をξと書くと、W(ξ)は量子状態の特徴を調べるための一つの「問い」に当たります。その期待値は一般に複素数なので、実部と虚部に対応する二つの実数の観測結果としても読めます。
定理の結論は、どの固定した非ゼロの局所係数ξについても、時間発展後のWeyl演算子の期待値がゼロへ近づくことです。一方、ξ＝0ならW(0)は恒等演算子で、期待値は常に1。この対応が、論文でいう局所Weyl代数上の「無限温度状態」です。
ここで無限温度と呼ぶのは、この観測量の集合に対する状態の性質です。実際の装置を無限の温度に加熱した、有限の箱の密度行列が一様分布へ収束した、という結果ではありません。位置や粒子数を含むあらゆる物理量が同じように収束するとも述べていません。
初期状態には、有限の領域ごとに通常の密度演算子で表せる「局所正規性」と、粒子密度が一様に有界であることを仮定します。後者は、どの有限領域を選んでも、その平均粒子数が「一定値ν×格子点数」を超えない条件です。場所ごとに粒子数が際限なく増える状態は除きますが、初期状態の並進対称性やガウス性、熱平衡状態であることは要求しません。
観測する範囲の広がりより、係数の伸びが速い
仕組みの中心は、Weyl演算子の期待値の絶対値に対する上限です。係数が非ゼロの格子点の集合をΓとすると、論文はおおよそ「粒子密度と格子点数から決まる平方根」を「係数の大きさ」で割った量で抑えられることを示します。正確には、π√(2ν|Γ|＋1)／‖ξ‖₂という上限です。
観測量を時間発展させると、関係するセルも増えます。ただし1回の更新が有限範囲なので、d次元格子でその点数はステップ数kの多項式程度、最大で定数倍の(1＋k)dに抑えられます。それに対して係数のノルム‖ξ‖₂が指数関数的に増えるなら、上限全体はゼロへ向かいます。
大きくなるのは、Weyl演算子を指定する位相空間の係数です。Weyl演算子そのものの作用素ノルムは1のままなので、「観測装置の出力が指数的に増幅される」と読み替えることはできません。また、観測量を時間発展させて期待値を計算するこの見方と、状態を時間発展させて固定した観測量を測る見方は対応しています。
伸びる方向があるだけでは、違いは消えない
更新によって係数が伸びる方向と縮む方向に分かれる性質を、双曲性と呼びます。論文は二組の十分条件を示しました。一つは、すべての波数で一様な双曲性があり、非ゼロの局所係数が縮む方向だけに属さないこと。もう一つは、一部の波数領域に双曲性があり、局所係数が「同じ形を定数倍して平行移動するだけ」で更新される例外を持たないことです。後者の例外を持たない性質を、論文では正則性と呼びます。
なぜ追加条件が要るのでしょうか。各セルで独立に同じスクイージングを行うだけなら、ある方向が伸びる一方で、局所的な係数が縮む方向にそのまま残れます。その係数に対するWeyl期待値はゼロではなく1へ近づく場合があります。双曲性があるという一言だけでは、論文の熱化を保証できません。
一次元の例で、条件の違いを見る
論文の第5.1節には、実数パラメーターaと左右の隣接セルへの結合で定める一次元の例があります。波数θごとの変換行列のトレースはa＋2cosθ、行列式は1。その固有値が単位円から離れる条件は、|a＋2cosθ|＞2です。次の表はこの例だけに当てはまります。
| パラメーター | 双曲性と正則性 | この論文の十分条件 |
|---|---|---|
| a＝0 | 局所双曲性も正則性も満たさない | 適用できない。形を保って移動する局所係数がある |
| a＝1 | 一部の波数で双曲的、かつ正則 | 局所双曲性と正則性の組で適用できる |
| a＝5 | 全波数で一様に双曲的、かつ正則 | より強い条件の組も満たす |
この族ではa≠0なら局所双曲性と正則性を満たし、|a|＞4なら一様双曲性も満たします。全波数で伸縮の条件をそろえなくても、局所的に指定した観測量に伸びる成分が必ず含まれれば、熱化を示せるところが見どころです。
ただし、格子点あたり一つのボソンモードという条件は大切です。論文自身が二モードの場合の反例を挙げており、正則性と他の条件の関係を、そのまま多モードの模型へ移すことはできません。収束の速さについても上限評価はありますが、実際の最適な熱化時間を一般に求めたわけではありません。
この研究から得られるのは、可逆な更新と、局所観測から初期状態の違いが見えなくなることを結ぶ、具体的な数学の道筋です。熱化という言葉を読むときには、系の種類に加えて「どの観測量が、どの初期状態から、何へ収束するのか」をそろえると、結果の意味がはっきりします。
参考資料
- Roman Geiko, Jake Gerenraich, Thermalization with Gaussian Quantum Cellular Automata, arXiv:2606.05542v2。初稿2026年6月4日、v2改訂2026年11月3日。定理1.1・1.2・1.5、第2～5節、特に系3.19と第5.1節、第6節。原論文PDF。