技術的要約：量子ラビ模型におけるもつれ生成能と対称性

問題提起
量子ラビ模型（QRM）は、量子ビットと単一の調和振動子が結合した系を記述する、光・物質相互作用の標準的な有効ハミルトニアンである。対称なQRM（$\epsilon=0$）は、明示的な離散 $\mathbb{Z}_2$ パリティ対称性を有するが、横バイアス $\epsilon$ によって特徴付けられる非対称量子ラビ模型（AQRM）はこのパリティを破る。しかし、AQRMは整数バイアス値（$\epsilon/\omega \in \mathbb{Z}$）において「隠れた」対称性を備えており、そこではスペクトルに例外的な交差が生じ、ヒルベルト空間が対称セクターへと整理される。励起数（$U(1)$ 対称性）が保存されるジェーンズ・カミングス（JC）モデルとは異なり、AQRMの隠れた対称性はハミルトニアンの構造には現れず、また、固定された結合パラメータ $(g, \Delta)$ において一般的にスペクトルの縮退を生じさせることもない。

本研究が取り組む中心的な問題は、一般的なスペクトルの縮退の兆候を残さないこのような隠れた対称性を、演算子レベルの診断手法、具体的には「時間平均もつれ生成能（time-averaged entangling power）」によって検出できるか否かである。さらに、無限次元の振動子ヒルベルト空間は、通常、有限次元のハール測度に依存する標準的なもつれ診断を困難にする。

手法
著者らは、診断ツールとして時間平均もつれ生成能（$e_p$）を用いる。この量は、積状態の入力からの時間発展演算子 $U(t) = e^{-iHt}$ によって生成される平均的なもつれを、時間および入力状態のアンサンブルの両方について平均化したものである。

無限次元の振動子を扱うために、本論文では2つの異なる有限入力アンサンブルを比較している。

1. 有限ハール平均： 振動子のヒルベルト空間を $N_{tr}$ 個のフォック状態に截断する。その上で、この有限部分空間に対してハール平均を行う。このアンサンブルの平均占有数は、カットオフに応じてスケールする（$\bar{n}_H \approx N_{tr}/2$）。
2. コヒーレント状態平均： 入力アンサンブルを、固定された平均占有数 $\bar{n}$ を持つコヒーレント状態 $|\alpha\rangle$ で構成する。このアプローチにより、物理的な入力スケール（$\bar{n}$）と、数値的な截断カットオフ（$N_{tr}$）を分離しており、後者はコヒーレント状態のサポートを捉えるのに十分な大きさに設定される。

もつれ生成能は、量子ビット部分系の簡約密度行列の線形エントロピーを用いて計算される。時間平均は、エネルギー差に従って固有状態のペアをグループ化する手法を用いて計算され、これにより、厳密な縮退および繰り返されるエネルギーギャップからの寄与を分離している。

主な貢献と結果

1. 隠れた対称性の検出：
   バイアスパラメータ $\epsilon/\omega$ に関する時間平均もつれ生成能の数値スキャンにより、整数値（$\epsilon/\omega \in \mathbb{Z}$）において鋭いピークが現れることが明らかになった。この特徴は、弱結合、強結合、および深強結合領域において頑健であり、有限ハールおよびコヒーレント状態の両方の診断において一貫して観察される。これは、もつれ生成能が、一般的なスペクトルの縮退の兆候がない場合でも、AQRMの隠れた対称性に反応することを裏付けている。

2. ピークのメカニズム（固有ベクトルの再編成）：
   本論文は、整数バイアスのピークが、対称性によるディップ（低下）に見られる典型的なスペクトルメカニズム（すなわち、縮退によるエネルギー差グループの拡大）によって引き起こされるのではないことを示している。代わりに、これらのピークは固有ベクトルの再編成から生じている。

• 変位二重項（displaced-doublet）の描像を用いることで、$\epsilon = k\omega$ において、変位した振動子ラダー間の混合が最大化されることを示している。
• この最大化された混合により、固有状態の簡約量子ビット・ブロッホベクトルの長さが原点に向かって収縮（短縮）し、これが量子ビットと振動子の間のより強いもつれを示唆する。
• 時間平均純度 $\langle \text{Tr}\rho_A^2 \rangle_t$ はこの収縮によって減少するため、もつれ生成能 $e_p = 2(1 - \langle \text{Tr}\rho_A^2 \rangle_t)$ は増加し、結果としてピークとなる。

3. ジェーンズ・カミングス極限との対比：
   AQRMのピークとは対照的に、JC点（$U(1)$ 対称性が顕在化している点）では、もつれ生成能に**弱いディップ（低下）**が見られる。

• JC対称性は、ヒルベルト空間を固定電荷セクターへと分解する。
• 対角成分のチャネル（ブロッホベクトルの長さ）は最小化されるが（ブロッホベクトルが消失する）、$U(1)$ セクターの疎な構造がスペクトル展開における非対角成分を集中させるため、正味のもつれ生成能の減少を招く。
• これは、応答の符号（ピークかディップか）が、対称性がスペクトル展開をどのように再編成するか（固有ベクトル構造か縮退のグループ化か）に依存することを浮き彫りにしている。

4. 診断の頑健性：
   整数バイスのピークの位置は、特定の入力アンサンブルの選択ではなく、AQRMの固有構造によって決定される。2つの診断法（有限ハールとコヒーレント）は振幅が異なり、ピークを解像するために異なるパラメータのスケーリングを必要とするが、ピークの位置については一致しており、この結果がハミルトニアンの動力学の特性であることを検証している。

意義
本論文は、時間平均もつれ生成能が、たとえ「隠れた」対称性であっても、一般的なスペクトルの縮退として現れないシステムにおける、光・物質相互作用系の対称性を捉えるための有効な演算子レベルの診断ツールとなることを主張している。

• 縮退を超えて： 本研究は、もつれ診断の有用性を、明らかな縮退を持つシステム（JCモデルやSU(2)対称性を持つスピン鎖など）から、パラメータ依存的で非顕在的な対称性を持つシステムへと拡張するものである。
• 固有ベクトルへの感度： もつれ生成能が、固有ベクトルの再編成（具体的には、変位二重項におけるブロッホベクトルの収縮）に対して敏感であることを確立し、隠れた対称性演算子 $J_k$ の動力学的シグネチャーを提供している。
• 手法的な有用性： 有限ハールサンプリングとコヒーレント状態サンプリングの比較は、コヒーレント状態平均が、無限次元のボゾン系におけるもつれに対する、物理的に動機付けられた有限エネルギーのプローブを提供し、物理的スケールと数値的アーティファクトを分離できることを示している。

これらの結果は、量子リソース（もつれなど）が、強化された対称性の点において極大化される可能性があるという広範な仮説を支持しており、量子系における基本パラメータの整理に関する新しい視点を提供している。