ディッケ模型

**ディッケ模型**（ディッケもけい、）とは、量子光学において光と物質の相互作用を扱う基礎的なモデルの一つ。この模型で「光」は単一の量子モードとして、「物質」は2準位系の集合として記述される。ディッケ模型によれば、光と物質の結合の強さが臨界値を超えると平均場相転移を経てが現れる。この相転移はに属し、の実験によって実証されている。超放射相転移はレーザー不安定性といくつかの点で類似しているが、異なるユニバーサリティクラスに属する相転移である。

## 解説
ディッケ模型は単一モードの空洞共振器と  個の2準位系（ 個の1/2スピン自由度としても同じ）の結合を記述する量子力学的なモデルである。1973年にとエリオット・H・リーブによって最初に提案された。その研究は自由空間における超放射光に関するロバート・H・ディッケの先駆的な業績に触発されていた。

量子力学におけるあらゆるモデルと同様に、ディッケ模型は量子状態の集合（ヒルベルト空間）と全エネルギーに関する演算子（ハミルトニアン）を含んでいる。ディッケ模型のヒルベルト空間は共振器の状態と2準位系集団の状態（正確にはそれらのテンソル積）からなる。共振器のヒルベルト空間は光子数  のフォック状態 |n\rangle によって張られる。それぞれの状態は真空状態  |n=0\rangle に正準生成消滅演算子  および  を作用させて作られる。これらの演算子は共振器に光子1個を追加もしくは除去するはたらきがある。一つの2準位系の状態は「アップ」および「ダウン」で呼びわけられ、スピン演算子 \vec\sigma_j = (\sigma^x_j,~\sigma^y_j,~\sigma^z_j) によって定義される。これらはスピン代数 [\sigma^x_j,~\sigma^y_k]=i\hbar\sigma^z_j\delta_{j,k} を満たす。ここで  は換算プランク定数、j = (0,1,2,...,N) は2準位系の番号を表す。

ディッケ模型のハミルトニアンは以下で与えられる。

 (a+a^\dagger) \sum_j \sigma^x_j\;
|}}

第1項は共振器のエネルギーを表しており、共振器内の光子1個分のエネルギー \hbar\omega_\mathrm{c}（\omega_\mathrm{c} は共振器の周波数）と光子数 n_\mathrm{c} =a^\dagger a の積で与えられる。第2項は2準位系のエネルギーを表しており、\hbar\omega_z が準位間のエネルギー差である。最後の項は共振器と2準位系の結合項にあたり、2準位系の個数  の平方根の逆数、および定数  の積に比例することが仮定されている。この仮定により N\to\infty の極限で相転移が起きる（下記参照）。結合項は回転項と反回転項の和として書くこともできる。回転項は a \sigma^+ + a^\dagger \sigma^- に比例し、励起数が保存される過程を表す。反回転項は a \sigma^- + a^\dagger \sigma^+ に比例する。ここで \sigma^\pm = \sigma^x \pm i \sigma^y はスピン昇降演算子である。

式のハミルトニアンはすべてのスピン自由度が同等である（準位間隔や共振器への結合がすべて等しい）と仮定している。それにより全スピン演算子が S^\alpha=\sum_{j=1}^N\sigma_j^\alpha と定義される（ \alpha=x,y,z ）。これはスピン代数 [S^x,S^y]=i\hbar S^z を満たす。これらの演算子を用いて式のハミルトニアンを書き直すと

(a+a^\dagger)S^x
|}}

となる。この記法では S\leq N/2 であるような単一のスピン  を扱えばよく、ヒルベルト空間のサイズが 2S+1 にしかならないため、 個の1/2スピンをサイズ 2^N のヒルベルト空間で取り扱うのに比べて数値解析が簡便になる。

ディッケ模型には以下の変換 \mathcal{P} に関する大域的な対称性がある。

\mathcal{P} の二乗は1となる（2度続けて変換を行うと元の状態に戻る）ため、固有値は 1 もしくは -1 の二通りである。この対称性と関連するとして全励起数のパリティ  P=(-1)^{N_{\mathrm{ex}}} がある。ここで

=a^\dagger a + \sum_{j=1}^N\sigma^z_j\;
|}}

である。このパリティ保存性は、ディッケ・ハミルトニアンの反回転項を除く各項が励起数を変化させず、反回転項も励起数を \pm 2 ずつしか変化させない事実からも確かめられる。ディッケ模型の状態はこの対称性が保たれている場合「通常相」と、対称性が自発的に敗れている場合「超放射相」と呼ばれる。

### 関連するモデル
ディッケ模型は量子光学におけるほかのモデルと密接な関係がある。特に、ディッケ模型における2準位系の個数を N=1 としたモデルはラビ模型と呼ばれる。またディッケ模型から反回転項を削除したモデルは、N=1 ならば、N>1 ならばと呼ばれる。これら二つのモデルは U(1) 対称性を持ち、励起数 N_{\mathrm{ex}} が保存される。この対称性が自発的に破れるとレーザー状態が生じる（後段参照）。

ディッケ模型とほかのモデルとの関係を以下の表に示す。

## 超放射相転移

ディッケ模型を用いた初期の研究は平衡状態の性質を対象にしていた。それらの研究は N\to\infty の極限（「熱力学的極限」とも呼ばれる）を考え、熱力学的な分配関数  Z=\exp(-H/k_\mathrm{B} T) を仮定していた。ここで k_\mathrm{B} はボルツマン定数、T は温度である。このとき、結合定数 \lambda が臨界値 \lambda_\mathrm{c} を超えるとディッケ模型が超放射相転移として知られる二次相転移を起こすことが発見された。ヘップとリーブによる最初の導出では反回転項が無視されていたため、扱われていたモデルは実際にはタヴィス＝カミングズ模型であった（前節参照）。その後、完全なディッケ模型の研究により、反回転項が含まれる場合にも相転移が起きるものの、臨界結合定数は異なることが判明した。

超放射相転移では式で定義されるパリティ対称性 \mathcal{P} が自発的に破れる。相転移の秩序変数は \langle{a}\rangle/\sqrt{N} である。熱力学的極限においては、秩序変数は通常相でゼロに近づき、超放射相では二つの取りうる値のどちらかに近づく。二つの値は互いに逆位相の共振器場状態に対応している（式に表されているように、互いに逆の  成分を持つスピン状態に対応していとも言える）。超放射相転移の近傍では、秩序変数は \lambda に \langle{a}\rangle/\sqrt{N}\sim(\lambda_\mathrm{c}-\lambda)^{-1/2} のように依存する。この依存性は平均場理論でいう臨界指数 \beta = 1/2 に対応する。

### 超放射相転移の平均場的な取り扱い
超放射相転移を理論的に扱うもっとも単純な方法は平均場近似である。共振器場の演算子をその期待値で置き換える近似であり、熱力学的極限においては厳密に正しい。このとき式のディッケ・ハミルトニアンはそれぞれ異なる2準位系に作用する独立な項の総和となり、それらを個別に対角化することができる。熱的平衡（前段参照）において2準位系ごとの自由エネルギー  は

=\alpha\right) = \omega_\mathrm{c} \alpha^2 - k_\mathrm{B} T \ln\left(2 \cosh\left(\frac{\sqrt{\omega_z^2+16\lambda^2\alpha^2}}{2 k_\mathrm{B} T}\right)\right)
|}}

となる。相転移の臨界結合定数は dF/d\alpha(\alpha=0)=0 の条件から

と求められる。\lambda の場合  の極小は一つで、\lambda>\lambda_\mathrm{c} ならば極小が二つに分かれる。T\to0 の極限を取ることで、ゼロ温度における超放射相転移の臨界結合定数が \lambda_\mathrm{c}=\sqrt{\omega_\mathrm{c} \omega_z}/2 と求められる。

## 半古典極限とカオス
### 半古典極限
原子のヒルベルト空間が対称な部分空間に限定されているディッケ模型では、以下に示す光子状態と原子状態のテンソル積によって相空間を構成することができる。光子はグラウバーのコヒーレント状態

で表される。ここでD(q,p)=e^{-\frac{S}{4}\left (q^2 + p^2\right)}\exp\left(\sqrt{\frac{S}{2}}\left(q+ip\right) a^\dagger\right)は変位演算子、\left\vert 0\right \rangle はフォック真空状態である。原子はSU(2)コヒーレント状態

で表される。ここで R(Q,P)=\left(1-\frac{Q^2+P^2}{4}\right)^S \exp\left(\sqrt{\frac{1}{4-Q^2-P^2}}\left(Q+iP\right) \frac{S^+}{\hbar}\right) はブロッホ球におけるであり、Q^2+P^2 \leq 4 の条件がある。\left \vert {S,-S}\right \rangle はすべての原子が基底状態にある状態である。これらのテンソル積により正準座標 (q,p) および (Q,P) を持つ4次元相空間が得られる。

この状態に対して式で与えられるディッケ・ハミルトニアンの期待値を取ることで、古典的なハミルトニアン

が得られる。N\to \infty の極限では、式の量子ハミルトニアンで表される量子的ダイナミクスと、式で与えられる古典的ダイナミクスが一致する。系のサイズが有限なら、古典的挙動と量子的挙動との対応は、N に反比例するある長さの時間（エーレンフェスト時間）を超えると破綻する。

### 量子カオス

ディッケ模型は量子描像と古典描像の対応や量子カオスを研究するための理想的な系となる。

式で与えられる古典系はパラメータ \lambda、\omega_\mathrm{c}、\omega_z やエネルギー E の値によってカオス的にも規則的にもふるまう。通常領域でも超放射領域でもカオスが生じうることには注意が必要である。

近年、カオス的領域において、ならびに規則的領域の不安定点において、非時間順序相関関数の指数関数的な成長速度が古典的なリアプノフ指数と一致することが発見された。さらに、極度に非局在化されたエネルギー固有基底を持つ初期コヒーレント状態の生存確率（すなわち、時間経過後の自分自身との）の時間発展はランダム行列理論によって良く記述される一方で、の影響を強く受ける初期コヒーレント状態ではエルゴード性が破られることが分かった。

## 開放ディッケ模型
式のディッケ模型は共振器モードと2準位系が外部環境から完全に孤立していることを前提としていた。実際の実験では、自由空間の光モードとの結合によって共振器光子の損失や2準位系の緩和（散逸チャネル）が生じるため、この前提は成り立たない。この種の実験では共振器モードと2準位系の間に結合を作るために駆動場（レーザー場）を用いていることにも注意が必要である。そのような散逸チャネルは環境自由度との結合を新たに導入することで記述される。それらの外部自由度のダイナミクスを平均することでの運動方程式が得られる。標準的なボルン＝マルコフ近似に従うなら、系のダイナミクスはの

によって記述される。ここで \rho はその系の密度行列、L_\alpha は緩和チャネル \alpha に関するリンドブラッド演算子、\gamma_\alpha は緩和レートである。ハミルトニアン H が式で与えられる場合、このモデルは開放ディッケ模型と呼ばれる。

実験において一般的に見られる緩和過程を以下に示す。

このモデルを用いた理論的な取り扱いでは d\rho/dt = 0 となる定常状態を考えることが多い。N\to\infty の極限において、解放ディッケ模型の定常状態は連続相転移を起こす。これは一般に非平衡超放射相転移と呼ばれる。この相転移の臨界指数は有限温度における平衡超放射相転移と等しい（したがってゼロ温度での超放射相転移とは異なる）。

### 超放射相転移とディッケ超放射

開放ディッケ模型における超放射相転移と関連しているが異なる現象にディッケ超放射がある。

ディッケ超放射とは自由空間に置かれた多数の2準位系がコヒーレントに光子を放出する集団的現象である。2準位系が初期状態として励起状態にあり、2準位系間の距離が関連する光子の波長よりも十分に小さい場合に起きる。これらの条件の下では2準位系の自発的な崩壊が非常に早くなり、大振幅の短い光パルスが放出される。理想的な条件においてパルスの持続時間は2準位系の数 N に反比例し、放出される光の最大強度は N^2 に比例する。これに対して、N 個の独立な2準位系からの自発放出では崩壊時間は N に依存せず、パルス強度は N に比例する。

前記の通り、開放ディッケ模型が対象としているのは量子化された空洞共振器と結合し、かつ外部からポンピングを受けている2準位系である。通常相において共振器場の強度は原子数 N に比例することはないが、超放射相においては \langle a^\dagger a \rangle \sim N が成り立つ。

ディッケ超放射、およびディッケ模型における超放射相転移のスケーリング則を以下の表に示す。

## 実験的実現

共振器との間に双極子相互作用がはたらく2準位原子はディッケ模型をもっとも単純に実現する。しかしこの系で超放射相転移を観測するには障害が二つ考えられる。(1) 通常、原子と共振器の結合は素のままでは弱く、式の臨界値を得るには不十分である。(2) この物理系を正確にモデル化するには A^2 項を取り入れなければいけないが、この項は（ある「no-go定理（何らかの物理現象の禁止に関する法則）」によれば）超放射相転移を阻害する可能性がある。これら二つの制約は、原子に外部ポンピングを施すとともに適切な回転座標系において有効なディッケ模型を構成することによって回避できる。

2010年、開放ディッケ模型における超放射相転移が光共振器中にトラップされた中性ルビジウム原子を用いて実験的に観測された。この種の実験では、原子系と共振器系の間の結合は、それらの間の直接的な双極子相互作用にはよらない。その代わりに原子は外部ポンピング光によって誘導ラマン遷移を誘起される。誘導ラマン遷移は2光子過程で、2準位系の状態がダウンからアップに（またはその逆に）遷移すると同時に共振器に光子が放出（または吸収）される。実験ではポンピング強度がある臨界値を超えると共振器中の光子数が急激に上昇することが確認され、しきい値がディッケ模型の臨界結合定数と関連付けられた。

この種の実験では、ダウン／アップ状態を表す物理的状態の選び方が2種類あった。いくつかの実験では、異なる運動量を持つ原子が二つの状態に対応しており、ダウン状態は運動量ゼロのボース＝アインシュタイン凝縮状態、アップ状態は共振器光子1個とポンピング光子1個の運動量の和に等しい運動量を持つ状態にあたる。後期の実験では異なる超微細構造準位を持つ磁場中ルビジウム原子が利用され、それにより一般化ディッケ模型（下記参照）の研究が可能になった。いずれの方式でも系は時間依存性を持ち、（一般化）ディッケ・ハミルトニアンはポンピング周波数で回転する座標系において実現されている。

## 一般化ディッケ模型とレーザー不安定性
ディッケ模型は式のハミルトニアンへの追加項を考えることで一般化できる。例として、2018年には独立に調整可能な回転項と反回転項を持つ開放ディッケ模型が実験的に実現された。この一般化ディッケ模型は超放射相転移に加えて "inverted lasing" や "counter-lasing" と呼ばれる特殊なレーザー不安定性を示す。この遷移はディッケ模型の反回転項によって引き起こされるもので、回転項よりも反回転項の方が大きい場合に顕著となる。

このレーザー不安定性は非平衡超放射相転移との間に類似点と相違点の両方を持っている。どちらも平均場的な転移であって1自由度のダイナミクスによって理解できるが、超放射相転移が超臨界ピッチフォーク分岐にあたるのに対し、レーザー不安定性はホップ分岐にあたる。二つの分岐の最大の違いは、前者では転移後に二つの安定解が生じる一方で、後者は周期的な解（リミットサイクル）に落ち込む点である。その結果、共振器場は超放射相においては（ポンプ場の座標系から見れば）静的であるが、レーザー相においては周期的な振動を示す。

## 脚注
この記事は以下の出典からCC BY 4.0ライセンスのもとで2020年に転載され（査読コメント）、2025年に翻訳された。
* Mor M Roses; Emanuele Dalla Torre (4 September 2020). . *[[PLOS ONE|PLOS One]]*. **15** (9): e0235197. doi:. ISSN&nbsp;1932-6203. PMID&nbsp;32886669. Wikidata&nbsp;.

### 注釈

### 出典

## 関連項目
*