ジェインズ・カミングスモデル

ジェインズ・カミングスモデルの図解。光共振器内の原子は左上に赤い点で示されている。共振器内の場のモードに結合する原子のエネルギー準位は、右下の円で示されている。2つの状態間の遷移により、原子は共振器モード内(外)に光子を放出(吸収)する。

量子光学において、ジェインズ・カミングスモデル(Jaynes–Cummings model 、 JCMと略されることもある)は、光共振器(またはボソン場)の量子化モードと相互作用する2準位原子の系を記述する理論モデルである。このモデルは回転波近似を仮定し、初期段階では散逸を無視し、光(自発放出および吸収を引き起こす可能性のある電磁放射の浴の形態)の有無にかかわらず、単一の場モードと2準位原子のみを扱う。このモデルはもともと、共振器での光子の自発放出および吸収現象を調べるために、量子化電磁場と原子の相互作用を研究するために開発された。1960年代のエドウィン・トンプソン・ジェインズフレッド・カミングスにちなんで命名され、1987年に実験的に確認された。

ジェインズ・カミングス模型は、原子物理学量子光学固体物理学量子情報回路において、実験的にも理論的にも大きな関心を集めています。50周年[ 1 ] 、60周年[2]を記念した特別号が発行されています。[ 2 ]また、コヒーレント制御量子情報処理にも応用されています。

歴史

1963年: ジェインズとカミングス

このモデルは、エドウィン・ジェインズフレッド・カミングスが1963年に発表した論文で、電磁場と相互作用する原子の挙動に完全な量子力学的処理を施すことの効果を明らかにするために考案されました。数学を簡素化し、扱いやすい計算を可能にするために、ジェインズとカミングスは、量子化された電磁場の単一モードと原子との相互作用に焦点を絞りました。 [ 3 ] [ 4 ](数学的な詳細については下記を参照)。

このアプローチは、従来の半古典的手法とは対照的です。半古典的手法では、原子のダイナミクスのみが量子力学的に扱われ、原子が相互作用する場は古典的な電磁気学理論に従って振舞うと仮定されます。ジェインズ・カミングス模型における場の量子力学的扱いは、以下を含む多くの新しい特徴を明らかにします。

  • 量子場と相互作用する二準位系の状態間のラビ振動の存在。これは当初純粋に量子力学的な効果であると考えられていたが、後に線形分散と吸収という観点から半古典的な説明が提示された[ 5 ]。
  • 量子化されたエネルギー準位の梯子はジェインズ・カミングス梯子と呼ばれ、エネルギーは非線形にスケーリングされます。ここで、は結合系における量子の総数です。このエネルギーの量子化と非線形スケーリングは、純粋に量子力学的な性質を持っています。n{\displaystyle {\sqrt {n}}}n{\displaystyle n}
  • 場が当初コヒーレント状態にあるとき、与えられた状態における二準位系を検出する確率の崩壊とその後の回復。崩壊は単純な古典的説明で説明できるが、回復は場の量子的性質に起因するエネルギースペクトルの離散性によってのみ説明できる。 [ 6 ] [ 7 ]

ジェインズ・カミングスモデルによって予測されるダイナミクスを実験的に実現するには、非常に高い品質係数を持つ量子力学共振器が必要である。これにより、2準位系(典型的には原子内の2つのエネルギーサブレベル)の状態間の遷移が、原子と場のモードとの相互作用によって非常に強く結合する。これにより、原子内の他のサブレベル間の結合や場の他のモードとの結合が同時に抑制され、損失はジェインズ・カミングスモデルによって予測されるダイナミクスを観測できるほど小さくなる。このような装置の実現が困難であったため、このモデルは長らく数学的な好奇心の対象となっていた。1985年には、いくつかのグループがマイクロ波空洞内のメーザーリュードベリ原子を用いて、予測されたラビ振動を実証した。[ 8 ] [ 9 ]しかし、前述のように、この効果は後に半古典的な説明ができることが判明した。[ 5 ]

1987年: レンペ、ワルサー、クライン

1987年になってようやく、ゲルハルト・レンペハーバート・ヴァルター、ノルベルト・クラインが単一原子メーザーを用いて、モデルによって予測された確率の復活を実証することができた。[ 10 ]それ以前は、研究グループは、原子と単一モードの結合を増強し、同時に他のモードを抑制できる実験装置を構築することができなかった。実験的には、システムのダイナミクスを単一モードフィールドのダイナミクスと同等と見なすには、共振器のQ値が十分に高くなければならない。量子力学的なフィールドモデルでしか説明できないダイナミクスのこの実証成功は、この研究に使用するための高品質共振器のさらなる開発を促した。

1原子メーザーの出現により、実験的観点から、単一原子(通常はリュードベリ原子)と空洞内の単一共鳴モードの電磁場との相互作用を研究することが可能となり、 [ 11 ] [ 12 ]ジェインズ・カミングスモデルのさまざまな側面を研究することができた。

砂時計形状は、モードが占める体積を最大化すると同時に、結合強度を最大化するために高い品質係数を維持し、モデルのパラメータをより適切に近似できることがわかった。[ 13 ]可視光周波数での強い原子場結合を観測するには、砂時計型の光モードが役立つ可能性がある。これは、その大きなモード体積が最終的にキャビティ内の強い場と一致するためである。[ 13 ]フォトニック結晶ナノキャビティ内の量子ドットも、可視光周波数でのラビサイクルの崩壊と復活を観測するための有望なシステムである。[ 14 ]

さらなる展開

最近の多くの実験は、量子情報処理とコヒーレント制御への応用が期待されるシステムへのモデルの適用に焦点が当てられている。様々な実験により、量子ドットとマイクロキャビティのモードの結合におけるジェインズ・カミングスモデルのダイナミクスが実証されており、はるかに小さいサイズの物理システムへの適用が可能になる可能性がある。 [ 15 ] [ 16 ] [ 17 ] [ 18 ]他の実験は、直接分光観測によってジェインズ・カミングスエネルギーレベルのラダーの非線形性を実証することに焦点を当てている。これらの実験は、超伝導RLC回路の形で非常に高品質の発振器に結合された人工原子を含む超伝導回路と、スピンを介して結合したリュードベリ原子の集合の両方において、場の量子的性質から予測される非線形挙動の直接的な証拠を発見した。[ 19 ] [ 20 ]後者の場合、集団内の集団リュードベリ励起の有無が2レベルシステムの役割を果たし、ボソン場モードの役割は発生するスピン反転の総数によって果たされる。[ 20 ]

理論的研究では、典型的には現象論的アプローチにより、元のモデルを拡張して散逸と減衰の効果を取り入れてきた。[ 21 ] [ 22 ] [ 23 ]提案された拡張では量子場の多重モードも組み込まれており、原子内の追加的なエネルギーレベルとの結合や、同じ場と相互作用する複数の原子の存在が可能になっている。通常用いられるいわゆる回転波近似を超える試みもなされている(下記の数学的導出を参照)[ 24 ] [ 25 ] [ 26 ]単一の量子場モードと複数()の2状態サブシステム(1/2 よりも高いスピンに相当)との結合は、ディッケモデルまたはタヴィス・カミングスモデルとして知られている。例えば、これは空洞共鳴付近の遷移を持つ複数の同一原子を含む高品質共鳴空洞や、超伝導回路上で複数の量子ドットに結合された共振器に適用される。これは、ケース の Jaynes-Cummings モデルに簡約されます。 >1{\displaystyle N>1}1{\displaystyle N=1}

このモデルは、実験環境においていくつかのエキゾチックな理論的可能性を実現する可能性を提供する。例えば、ラビ振動が崩壊する期間、原子-空洞系はマクロスケールで量子重ね合わせ状態にあることが分かった。このような状態は、マクロスケールの系における量子もつれの発現に関する直感に反する効果の探究を可能にするため、シュレーディンガーの猫と呼ばれることもある。 [ 27 ]また、量子場における量子情報伝達のモデル化にも用いることができる。 [ 28 ]

数式1

完全なシステムを記述するハミルトニアンは、 自由場ハミルトニアン、原子励起ハミルトニアン、およびジェインズ・カミングス相互作用ハミルトニアンから構成されます。 H^H^分野+H^原子+H^整数{\displaystyle {\hat {H}}={\hat {H}}_{\text{field}}+{\hat {H}}_{\text{atom}}+{\hat {H}}_{\text{int}}}H^分野ωc1つの^1つの^H^原子ω1つのσ^z2H^整数Ω2E^S^{\displaystyle {\begin{aligned}{\hat {H}}_{\text{field}}&=\hbar \omega _{c}{\hat {a}}^{\dagger}{\hat {a}}\\{\hat {H}}_{\text{atom}}&=\hbar \omega _{a}{\frac {{\hat {\sigma }}_{z}}{2}}\\{\hat {H}}_{\text{int}}&={\frac {\hbar \Omega }{2}}{\hat {E}}{\hat {S}}.\end{aligned}}}

ここでは便宜上、真空場エネルギーは に設定されています。 0{\displaystyle 0}

JCM相互作用ハミルトニアンを導出するために、量子化された放射場は、場の演算子 を持つ 単一のボソンモードで構成されるとみなされます。ここで、演算子およびはボソン生成消滅演算子であり、はモードの角周波数です。一方、2レベル原子は、3次元ブロッホベクトルを用いて状態を記述できるスピンハーフと同等です。(ここでの「2レベル原子」は、スピンを持つ実際の原子ではなく、ヒルベルト空間がスピンハーフ同型である一般的な2レベル量子系であることを理解する必要があります。)原子は、分極演算子 を介して場に結合されます。演算子および は、原子の上昇演算子と下降演算子です。演算子 は原子反転演算子であり、は原子遷移周波数です。 E^EZPF1つの^+1つの^{\displaystyle {\hat {E}}=E_{\text{ZPF}}\left({\hat {a}}+{\hat {a}}^{\dagger }\right)}1つの^{\displaystyle {\hat {a}}^{\dagger }}1つの^{\displaystyle {\hat {a}}}ωc{\displaystyle \omega_{c}}S^σ^++σ^{\displaystyle {\hat {S}}={\hat {\sigma }}_{+}+{\hat {\sigma }}_{-}}σ^+|eグラム|{\displaystyle {\hat {\sigma }}_{+}=|e\rangle \langle g|}σ^|グラムe|{\displaystyle {\hat {\sigma }}_{-}=|g\rangle \langle e|}σ^z|ee||グラムグラム|{\displaystyle {\hat {\sigma }}_{z}=|e\rangle \langle e|-|g\rangle \langle g|}ω1つの{\displaystyle \omega _{a}}

ジェインズ・カミングス・ハミルトニアン 1

シュレーディンガー描像から選択によって定義される 相互作用描像(回転座標系ともいう)に移ると、次の式が得られる。 H^0H^分野+H^原子{\displaystyle {\hat {H}}_{0}={\hat {H}}_{\text{field}}+{\hat {H}}_{\text{atom}}}H^整数tΩ21つの^σ^eωc+ω1つのt+1つの^σ^+eωc+ω1つのt+1つの^σ^+eωc+ω1つのt+1つの^σ^eωc+ω1つのt{\displaystyle {\hat {H}}_{\text{int}}(t)={\frac {\hbar \Omega }{2}}\left({\hat {a}}{\hat {\sigma }}_{-}e^{-i(\omega _{c}+\omega _{a})t}+{\hat {a}}^{\dagger }{\hat {\sigma }}_{+}e^{i(\omega _{c}+\omega _{a})t}+{\hat {a}}{\hat {\sigma }}_{+}e^{i(-\omega _{c}+\omega _{a})t}+{\hat {a}}^{\dagger }{\hat {\sigma }}_{-}e^{-i(-\omega _{c}+\omega _{a})t}\right).}

このハミルトニアンには、速く振動する成分と遅く振動する成分の両方が含まれています。解けるモデルを得るために、速く振動する「逆回転」項 は無視されます。これは回転波近似と呼ばれ、速く振動する項が比較的大きなエネルギー差を持つ状態を結合するため有効です。エネルギー差が結合よりもはるかに大きい場合、これらの状態の混合は小さくなります。言い換えれば、結合は状態間のポピュレーション移動をほとんど引き起こしません。シュレーディンガー描像に戻すと、JCMハミルトニアンは次のように書き表されます。 ωc+ω1つの{\displaystyle (\omega _{c}+\omega _{a})}ωcω1つの{\displaystyle (\omega _{c}-\omega _{a})}ωc+ω1つの{\displaystyle (\omega _{c}+\omega _{a})}H^JCωc1つの^1つの^+ω1つのσ^z2+Ω21つの^σ^++1つの^σ^{\displaystyle {\hat {H}}_{\text{JC}}=\hbar \omega _{c}{\hat {a}}^{\dagger }{\hat {a}}+\hbar \omega _{a}{\frac {{\hat {\sigma }}_{z}}{2}}+{\frac {\hbar \Omega }{2}}\left({\hat {a}}{\hat {\sigma }}_{+}+{\hat {a}}^{\dagger }{\hat {\sigma }}_{-}\right).}

固有状態

完全なシステムのハミルトニアンを 2 つの可換部分の合計として記述することが可能であり、多くの場合非常に役立ちます。 ここで 、 はフィールドと 2 レベル システム間の デチューニング(周波数)と呼ばれます。H^JCH^+H^II{\displaystyle {\hat {H}}_{\text{JC}}={\hat {H}}_{\text{I}}+{\hat {H}}_{\text{II}},}H^ωc1つの^1つの^+σ^z2H^IIδσ^z2+Ω21つの^σ^++1つの^σ^{\displaystyle {\begin{aligned}{\hat {H}}_{\text{I}}&=\hbar \omega _{c}\left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {{\hat {\sigma }}_{z}}{2}}\right)\\{\hat {H}}_{\text{II}}&=\hbar \delta {\frac {{\hat {\sigma }}_{z}}{2}}+{\frac {\hbar \Omega }{2}}\left({\hat {a}}{\hat {\sigma }}_{+}+{\hat {a}}^{\dagger }{\hat {\sigma }}_{-}\right)\end{aligned}}}δω1つのωc{\displaystyle \delta =\omega _{a}-\omega _{c}}

の固有状態はテンソル積の形式であるため簡単に解くことができ、 と表されます。ここで、 はモード内の放射量子の数を表します。 H^{\displaystyle {\hat {H}}_{I}}|n+1グラム|ne{\displaystyle |n+1,g\rangle ,|n,e\rangle }n{\displaystyle n\in \mathbb {N} }

状態 とはすべてのに対して に関して退化しているので、の部分空間 において対角化すれば十分である。この部分空間におけるの行列要素は、 |ψ1n:=|ne{\displaystyle |\psi _{1n}\rangle :=|n,e\rangle }|ψ2n:=|n+1グラム{\displaystyle |\psi _{2n}\rangle :=|n+1,g\rangle }H^{\displaystyle {\hat {H}}_{I}}n{\displaystyle n}H^JC{\displaystyle {\hat {H}}_{\text{JC}}}スパン{|ψ1n|ψ2n}{\displaystyle \operatorname {span} \{|\psi _{1n}\rangle ,|\psi _{2n}\rangle \}}H^JC{\displaystyle {\hat {H}}_{\text{JC}}}Hjn:=ψn|H^JC|ψjn{\displaystyle {H}_{ij}^{(n)}:=\langle \psi _{in}|{\hat {H}}_{\text{JC}}|\psi _{jn}\rangle ,}Hnnωc+ω1つの2Ω2n+1Ω2n+1n+1ωcω1つの2{\displaystyle H^{(n)}=\hbar {\begin{pmatrix}n\omega _{c}+{\frac {\omega _{a}}{2}}&{\frac {\Omega }{2}}{\sqrt {n+1}}\\[8pt]{\frac {\Omega }{2}}{\sqrt {n+1}}&(n+1)\omega _{c}-{\frac {\omega _{a}}{2}}\end{pmatrix}}}

与えられた に対して、 のエネルギー固有値はとなる。 ここで は特定の離調パラメータに対するラビ周波数である。エネルギー固有値に関連付けられた固有状態は で与えられる 。 ここで角度は によって定義される。 n{\displaystyle n}Hn{\displaystyle H^{(n)}}E±nωcn+12±12Ωnδ{\displaystyle E_{\pm }(n)=\hbar \omega _{c}\left(n+{\frac {1}{2}}\right)\pm {\frac {1}{2}}\hbar \Omega _{n}(\delta ),}Ωn(δ)=δ2+Ω2(n+1){\textstyle \Omega _{n}(\delta )={\sqrt {\delta ^{2}+\Omega ^{2}(n+1)}}}|n,±{\displaystyle |n,\pm \rangle }|n,+=cos(αn2)|ψ1n+sin(αn2)|ψ2n{\displaystyle |n,+\rangle =\cos \left({\frac {\alpha _{n}}{2}}\right)|\psi _{1n}\rangle +\sin \left({\frac {\alpha _{n}}{2}}\right)|\psi _{2n}\rangle }|n,=sin(αn2)|ψ1ncos(αn2)|ψ2n{\displaystyle |n,-\rangle =\sin \left({\frac {\alpha _{n}}{2}}\right)|\psi _{1n}\rangle -\cos \left({\frac {\alpha _{n}}{2}}\right)|\psi _{2n}\rangle }αn{\displaystyle \alpha _{n}}αn:=tan1(Ωn+1δ).{\displaystyle \alpha _{n}:=\tan ^{-1}\left({\frac {\Omega {\sqrt {n+1}}}{\delta }}\right).}

シュレーディンガー描像力学

一般状態のダイナミクスは、これを前述の固有状態へ展開することで得られる。場の初期状態として、状態の重ね合わせを考え、励起状態の原子が場に注入されると仮定する。系の初期状態は |ψfield(0)=nCn|n{\textstyle |\psi _{\text{field}}(0)\rangle =\sum _{n}{C_{n}|n\rangle }}|ψtot(0)=nCn|n,e=nCn[cos(αn2)|n,++sin(αn2)|n,].{\displaystyle |\psi _{\text{tot}}(0)\rangle =\sum _{n}{C_{n}|n,e\rangle }=\sum _{n}C_{n}\left[\cos \left({\frac {\alpha _{n}}{2}}\right)|n,+\rangle +\sin \left({\frac {\alpha _{n}}{2}}\right)|n,-\rangle \right].}

は場-原子系の定常状態な ので、時間に対する状態ベクトルは次のように与えられる。 |n,±{\displaystyle |n,\pm \rangle }t>0{\displaystyle t>0}|ψtot(t)=eiH^JCt/|ψtot(0)=nCn[cos(αn2)|n,+eiE+(n)t/+sin(αn2)|n,eiE(n)t/].{\displaystyle |\psi _{\text{tot}}(t)\rangle =e^{-i{\hat {H}}_{\text{JC}}t/\hbar }|\psi _{\text{tot}}(0)\rangle =\sum _{n}C_{n}\left[\cos \left({\frac {\alpha _{n}}{2}}\right)|n,+\rangle e^{-iE_{+}(n)t/\hbar }+\sin \left({\frac {\alpha _{n}}{2}}\right)|n,-\rangle e^{-iE_{-}(n)t/\hbar }\right].}

ラビ振動は、状態ベクトルのsin関数とcos関数で容易に観察できます。光子の状態数によって周期は異なります。実験で観測されるのは、非常に広範囲に振動し、ある瞬間には相殺的に和がゼロになる多くの周期関数の和です。しかし、後の瞬間には再びゼロではなくなります。この瞬間の有限性は、周期性の議論の離散性から生じます。もし電場の振幅が連続的であれば、有限時間でこの復活は決して起こらなかったでしょう。

ハイゼンベルクの描像力学

ハイゼンベルク記法では、ハミルトニアンからユニタリー発展演算子を直接決定することができる。[ 29 ] ここで演算子は次のように定義され 、次のように与えられる。 U^(t)=eiH^JCt/=(eiωct(a^a^+12)(costφ^+g2iδ/2sintφ^+g2φ^+g2)igeiωct(a^a^+12)sintφ^+g2φ^+g2a^igeiωct(a^a^12)sintφ^φ^a^eiωct(a^a^12)(costφ^+iδ/2sintφ^φ^)){\displaystyle {\begin{matrix}{\begin{aligned}{\hat {U}}(t)&=e^{-i{\hat {H}}_{\text{JC}}t/\hbar }\\&={\begin{pmatrix}e^{-i\omega _{c}t\left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right)}\left(\cos t{\sqrt {{\hat {\varphi }}+g^{2}}}-i\delta /2{\frac {\sin t{\sqrt {{\hat {\varphi }}+g^{2}}}}{\sqrt {{\hat {\varphi }}+g^{2}}}}\right)&-ige^{-i\omega _{c}t\left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right)}{\frac {\sin t{\sqrt {{\hat {\varphi }}+g^{2}}}}{\sqrt {{\hat {\varphi }}+g^{2}}}}\,{\hat {a}}\\-ige^{-i\omega _{c}t\left({\hat {a}}^{\dagger }{\hat {a}}-{\frac {1}{2}}\right)}{\frac {\sin t{\sqrt {\hat {\varphi }}}}{\sqrt {\hat {\varphi }}}}{\hat {a}}^{\dagger }&e^{-i\omega _{c}t\left({\hat {a}}^{\dagger }{\hat {a}}-{\frac {1}{2}}\right)}\left(\cos t{\sqrt {\hat {\varphi }}}+i\delta /2{\frac {\sin t{\sqrt {\hat {\varphi }}}}{\sqrt {\hat {\varphi }}}}\right)\end{pmatrix}}\end{aligned}}\end{matrix}}}φ^{\displaystyle {\hat {\varphi }}}φ^=g2a^a^+δ2/4{\displaystyle {\hat {\varphi }}=g^{2}{\hat {a}}^{\dagger }{\hat {a}}+\delta ^{2}/4}g{\displaystyle g}g=Ω{\displaystyle g={\frac {\Omega }{\hbar }}}

のユニタリー性は、恒等式 とそれらのエルミート共役によって保証されます。 U^{\displaystyle {\hat {U}}}sintφ^+g2φ^+g2a^=a^sintφ^φ^,costφ^+g2a^=a^costφ^,{\displaystyle {\begin{aligned}{\frac {\sin t\,{\sqrt {{\hat {\varphi }}+g^{2}}}}{\sqrt {{\hat {\varphi }}+g^{2}}}}\;{\hat {a}}&={\hat {a}}\;{\frac {\sin t\,{\sqrt {\hat {\varphi }}}}{\sqrt {\hat {\varphi }}}},\\\cos t\,{\sqrt {{\hat {\varphi }}+g^{2}}}\;{\hat {a}}&={\hat {a}}\;\cos t{\sqrt {\hat {\varphi }}},\end{aligned}}}

ユニタリ進化演算子を使用すると、密度行列によって記述されるシステムの状態の時間発展を計算し、そこから初期状態が与えられた場合の任意の観測可能な値の期待値を計算することが可能です。 ρ^(t){\displaystyle {\hat {\rho }}(t)}ρ^(t)=U^(t)ρ^(0)U^(t){\displaystyle {\hat {\rho }}(t)={\hat {U}}^{\dagger }(t){\hat {\rho }}(0){\hat {U}}(t)}Θ^t=Tr[ρ^(t)Θ^]{\displaystyle \langle {\hat {\Theta }}\rangle _{t}={\text{Tr}}[{\hat {\rho }}(t){\hat {\Theta }}]}

システムの初期状態は で表され、は観測可能量を表す演算子です。 ρ^(0){\displaystyle {\hat {\rho }}(0)}Θ^{\displaystyle {\hat {\Theta }}}

数式2

説明を簡単にするために、原子の2つのエネルギーサブレベルと量子化された電磁場との相互作用を考えてみましょう。ボソン場と結合した他の任意の2状態系の挙動は、このダイナミクスと同型になります。この場合、原子-場系のハミルトニアンは[ 30 ] です 。ここで、以下の定義を行いました。 H^=H^A+H^F+H^AF{\displaystyle {\hat {H}}={\hat {H}}_{A}+{\hat {H}}_{F}+{\hat {H}}_{AF}}

  • H^A=Eg|gg|+Ee|ee|{\displaystyle {\hat {H}}_{A}=E_{g}|g\rangle \langle g|+E_{e}|e\rangle \langle e|}は原子のハミルトニアンであり、文字はそれぞれ励起状態と基底状態を表します。エネルギーのゼロを原子の基底状態エネルギーに設定すると、これは次のように簡略化されます。ここで、は原子のサブレベル間の遷移の共鳴周波数です。e,g{\displaystyle e,g}H^A=Ee|ee|=ωeg|ee|{\displaystyle {\hat {H}}_{A}=E_{e}|e\rangle \langle e|=\hbar \omega _{eg}|e\rangle \langle e|}ωeg{\displaystyle \omega _{eg}}
  • H^F=k,λωk(a^k,λa^k,λ+12){\displaystyle {\hat {H}}_{F}=\sum _{\mathbf {k} ,\lambda }\hbar \omega _{\mathbf {k} }\left({\hat {a}}_{\mathbf {k} ,\lambda }^{\dagger }{\hat {a}}_{\mathbf {k} ,\lambda }+{\frac {1}{2}}\right)}は量子化された電磁場のハミルトニアンです。すべての可能な波数ベクトルと2つの可能な直交偏光状態にわたる無限和 に注目してください。演算子と は、場の各インデックス付きモードに対する光子生成および消滅演算子です。ジェインズ・カミングスモデルの単純さは、この一般的な和を抑制し、場の単一モードのみを考慮することで、 と書くことができることにあります。ここで、添え字は空洞の共鳴モードのみを考慮していることを示しています。k{\displaystyle \mathbf {k} }λ{\displaystyle \lambda }a^k,λ{\displaystyle {\hat {a}}_{\mathbf {k} ,\lambda }^{\dagger }}a^k,λ{\displaystyle {\hat {a}}_{\mathbf {k} ,\lambda }}H^F=ωc(a^ca^c+12){\textstyle {\hat {H}}_{F}=\hbar \omega _{c}\left({\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}+{\frac {1}{2}}\right)}c{\displaystyle c}
  • H^AF=d^E^(R){\displaystyle {\hat {H}}_{AF}=-{\hat {\mathbf {d} }}\cdot {\hat {\mathbf {E} }}(\mathbf {R} )}は双極子原子場相互作用ハミルトニアンです(ここで原子の位置は です)。量子化された電磁場の電場演算子は で与えられ、双極子演算子は で与えられます。 を設定して定義すると、は直交場モードであり、とは原子の部分空間に作用する上昇演算子と下降演算子です。ジェインズ・カミングス模型を適用することで、この和を抑制し、場の単一モードにのみ注目することができます。したがって、原子場ハミルトニアンは となります。R{\displaystyle \mathbf {R} }E^(R)=ik,λ2πωkVuk,λ(a^k,λeikRa^k,λeikR){\displaystyle {\hat {\mathbf {E} }}(\mathbf {R} )=i\sum _{\mathbf {k} ,\lambda }{\sqrt {\frac {2\pi \hbar \omega _{\mathbf {k} }}{V}}}\mathbf {u} _{\mathbf {k} ,\lambda }\left({\hat {a}}_{\mathbf {k} ,\lambda }e^{i\mathbf {k} \cdot \mathbf {R} }-{\hat {a}}_{\mathbf {k} ,\lambda }^{\dagger }e^{-i\mathbf {k} \cdot \mathbf {R} }\right)}d^=σ^+e|d^|g+σ^g|d^|e{\displaystyle {\hat {\mathbf {d} }}={\hat {\sigma }}_{+}\langle e|{\hat {\mathbf {d} }}|g\rangle +{\hat {\sigma }}_{-}\langle g|{\hat {\mathbf {d} }}|e\rangle }R=0{\displaystyle \mathbf {R} =\mathbf {0} }gk,λ=i2πωkVe|d^|guk,λ,{\displaystyle \hbar g_{\mathbf {k} ,\lambda }=i{\sqrt {\frac {2\pi \hbar \omega _{\mathbf {k} }}{V}}}\langle e|{\hat {\mathbf {d} }}|g\rangle \cdot \mathbf {u} _{\mathbf {k} ,\lambda },}uk,λ{\displaystyle \mathbf {u} _{\mathbf {k} ,\lambda }}H^AF=k,λ(gk,λσ^+a^k,λgk,λσ^a^k,λgk,λσ^+a^k,λ+gk,λσ^a^k,λ),{\displaystyle {\hat {H}}_{AF}=-\sum _{\mathbf {k} ,\lambda }\hbar \left(g_{\mathbf {k} ,\lambda }{\hat {\sigma }}_{+}{\hat {a}}_{\mathbf {k} ,\lambda }-g_{\mathbf {k} ,\lambda }^{*}{\hat {\sigma }}_{-}{\hat {a}}_{\mathbf {k} ,\lambda }^{\dagger }-g_{\mathbf {k} ,\lambda }{\hat {\sigma }}_{+}{\hat {a}}_{\mathbf {k} ,\lambda }^{\dagger }+g_{\mathbf {k} ,\lambda }^{*}{\hat {\sigma }}_{-}{\hat {a}}_{\mathbf {k} ,\lambda }\right),}σ^+=|eg|{\displaystyle {\hat {\sigma }}_{+}=|e\rangle \langle g|}σ^=|ge|{\displaystyle {\hat {\sigma }}_{-}=|g\rangle \langle e|}{|e,|g}{\displaystyle \{|e\rangle ,|g\rangle \}}H^AF=[(gcσ^+a^cgcσ^a^c)+(gcσ^+a^c+gcσ^a^c)]{\displaystyle {\hat {H}}_{AF}=\hbar \left[\left(g_{c}{\hat {\sigma }}_{+}{\hat {a}}_{c}-g_{c}^{*}{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }\right)+\left(-g_{c}{\hat {\sigma }}_{+}{\hat {a}}_{c}^{\dagger }+g_{c}^{*}{\hat {\sigma }}_{-}{\hat {a}}_{c}\right)\right]}

回転フレームと回転波近似

次に、いわゆる「共回転」系への受動的な変換を行うことで、解析を簡略化できます。これを行うには、相互作用描像を使用します。 とします 。すると、相互作用ハミルトニアンは 次のようになります。ここで、空洞の共鳴周波数が原子の遷移周波数に近い、つまり であると仮定します。この条件下では、 で振動する指数項はほぼ共鳴し、 で振動する他の指数項はほぼ反共鳴します。共鳴項が1回の完全な振動を完了するのにかかる時間で、反共鳴項は多くの完全なサイクルを完了します。反共鳴振動の各完全なサイクルにおいて、共鳴挙動を解析したい時間スケールにおいて、高速で振動する反共鳴項の正味の効果は平均して0になる傾向があるため、反共鳴項の値はほぼ共鳴項の値と比較して無視できるため、反共鳴項を完全に無視することができます。この近似は回転波近似として知られており、エネルギーは保存されるという直感と一致します。すると、相互作用ハミルトニアン(簡潔にするために実数とします)は次のようになります。 H^0=H^A+H^F{\displaystyle {\hat {H}}_{0}={\hat {H}}_{A}+{\hat {H}}_{F}}H^AF(t)=eiH^0t/H^AFeiH^0t/=(gcσ^+a^cei(ωc+ωeg)t+gcσ^a^cei(ωc+ωeg)tgcσ^a^cei(ωegωc)tgcσ^+a^cei(ωegωc)t){\displaystyle {\hat {H}}_{AF}(t)=e^{i{\hat {H}}_{0}t/\hbar }{\hat {H}}_{AF}e^{-i{\hat {H}}_{0}t/\hbar }=\hbar \left(g_{c}{\hat {\sigma }}_{+}{\hat {a}}_{c}^{\dagger }e^{i(\omega _{c}+\omega _{eg})t}+g_{c}^{*}{\hat {\sigma }}_{-}{\hat {a}}_{c}e^{-i(\omega _{c}+\omega _{eg})t}-g_{c}^{*}{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }e^{-i(\omega _{eg}-\omega _{c})t}-g_{c}{\hat {\sigma }}_{+}{\hat {a}}_{c}e^{i(\omega _{eg}-\omega _{c})t}\right)}|ωegωc|ωeg+ωc{\displaystyle |\omega _{eg}-\omega _{c}|\ll \omega _{eg}+\omega _{c}}ωegωc0{\displaystyle \omega _{eg}-\omega _{c}\simeq 0}ωeg+ωc2ωc{\displaystyle \omega _{eg}+\omega _{c}\simeq 2\omega _{c}}τ=2πΔ,Δωegωc{\displaystyle \tau ={\frac {2\pi }{\Delta }},\Delta \equiv \omega _{eg}-\omega _{c}}2π2ωcτ{\displaystyle {\frac {2\pi }{2\omega _{c}}}\ll \tau }gc{\displaystyle g_{c}}

H^AF(t)=gc(σ^+a^cei(ωegωc)t+σ^a^cei(ωegωc)t){\displaystyle {\hat {H}}_{AF}(t)=-\hbar g_{c}\left({\hat {\sigma }}_{+}{\hat {a}}_{c}e^{i(\omega _{eg}-\omega _{c})t}+{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }e^{-i(\omega _{eg}-\omega _{c})t}\right)}

この近似値を使って(そして負の符号を に吸収して)、シュレーディンガー像に戻すことができます。 gc{\displaystyle g_{c}}

H^AF=eiH^0t/H^AF(t)eiH^0t/=gc(σ^+a^c+σ^a^c){\displaystyle {\hat {H}}_{AF}=e^{-i{\hat {H}}_{0}t/\hbar }{\hat {H}}_{AF}(t)e^{i{\hat {H}}_{0}t/\hbar }=\hbar g_{c}\left({\hat {\sigma }}_{+}{\hat {a}}_{c}+{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }\right)}

ジェインズ・カミングス・ハミルトニアン 2

前の2つの節で得られた結果を用いて、完全なジェインズ・カミングス・ハミルトニアンを書き表すことができる。[ 30 ] 定数項は場の零点エネルギーを表す。これは力学には寄与しないので無視することができ、以下のようになる。 H^JC=ωc(a^ca^c+12)+ωeg|ee|+gc(σ^+a^c+σ^a^c){\displaystyle {\hat {H}}_{JC}=\hbar \omega _{c}\left({\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}+{\frac {1}{2}}\right)+\hbar \omega _{eg}|e\rangle \langle e|+\hbar g_{c}\left({\hat {\sigma }}_{+}{\hat {a}}_{c}+{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }\right)}12ωc{\displaystyle {\frac {1}{2}}\hbar \omega _{c}}H^JC=ωca^ca^c+ωeg|ee|+gc(σ^+a^c+σ^a^c){\displaystyle {\hat {H}}_{JC}=\hbar \omega _{c}{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}+\hbar \omega _{eg}|e\rangle \langle e|+\hbar g_{c}\left({\hat {\sigma }}_{+}{\hat {a}}_{c}+{\hat {\sigma }}_{-}{\hat {a}}_{c}^{\dagger }\right)}

次に、いわゆる数演算子を次のように 定義します。この演算子と原子場ハミルトニアンとの 交換子を考えます。N^=|ee|+a^ca^c{\displaystyle {\hat {N}}=|e\rangle \langle e|+{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}}[H^AF,N^]=gc([a^cσ^+,|ee|+a^ca^c]+[a^cσ^,|ee|+a^ca^c])=gc(a^c[σ^+,|ee|]+[a^c,a^ca^c]σ^++a^c[σ^,|ee|]+[a^c,a^ca^c]σ^)=gc(a^cσ^++a^cσ^++a^cσ^a^cσ^)=0{\displaystyle {\begin{aligned}\left[{\hat {H}}_{AF},{\hat {N}}\right]&=\hbar g_{c}\left(\left[{\hat {a}}_{c}{\hat {\sigma }}_{+},|e\rangle \langle e|+{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}\right]+\left[{\hat {a}}_{c}^{\dagger }{\hat {\sigma }}_{-},|e\rangle \langle e|+{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}\right]\right)\\&=\hbar g_{c}\left({\hat {a}}_{c}\left[{\hat {\sigma }}_{+},|e\rangle \langle e|\right]+\left[{\hat {a}}_{c},{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}\right]{\hat {\sigma }}_{+}+{\hat {a}}_{c}^{\dagger }\left[{\hat {\sigma }}_{-},|e\rangle \langle e|\right]+\left[{\hat {a}}_{c}^{\dagger },{\hat {a}}_{c}^{\dagger }{\hat {a}}_{c}\right]{\hat {\sigma }}_{-}\right)\\&=\hbar g_{c}\left(-{\hat {a}}_{c}{\hat {\sigma }}_{+}+{\hat {a}}_{c}{\hat {\sigma }}_{+}+{\hat {a}}_{c}^{\dagger }{\hat {\sigma }}_{-}-{\hat {a}}_{c}^{\dagger }{\hat {\sigma }}_{-}\right)\\&=0\end{aligned}}}

このように、数演算子は原子場ハミルトニアンと可換である。数演算子の固有状態は、場の状態が特定の数の光子を持つ状態であるテンソル積状態 の基礎となる。数演算子は、原子場系における量子の 総数を数える。{|g,0;|e,0,|g,1;;|e,n1,|g,n}{\displaystyle \left\{|g,0\rangle ;|e,0\rangle ,|g,1\rangle ;\cdots ;|e,n-1\rangle ,|g,n\rangle \right\}}{|n}{\displaystyle \left\{|n\rangle \right\}}n{\displaystyle n}N^{\displaystyle {\hat {N}}}n{\displaystyle n}

この(全状態数)の固有状態の基底において、ハミルトニアンはブロック対角構造をとる:[ 30 ]N^{\displaystyle {\hat {N}}}H^JC=[H00000H^10000H^200H^n0]{\displaystyle {\hat {H}}_{JC}={\begin{bmatrix}H_{0}&0&0&0&\cdots &\cdots &\cdots \\0&{\hat {H}}_{1}&0&0&\ddots &\ddots &\ddots \\0&0&{\hat {H}}_{2}&0&\ddots &\ddots &\ddots \\\vdots &\ddots &\ddots &\ddots &\ddots &\ddots &\ddots \\\vdots &\ddots &\ddots &0&{\hat {H}}_{n}&0&\ddots \\\vdots &\ddots &\ddots &\ddots &\ddots &\ddots &\ddots \\\end{bmatrix}}}

スカラーを除いて、対角線上の各 は、それ自体が形式の行列です。 H0{\displaystyle H_{0}}H^n{\displaystyle {\hat {H}}_{n}}2×2{\displaystyle 2\times 2}H^n=[ωc(n1)+ωege,n1|H^JC|g,ng,n|H^JC|e,n1nωc]{\displaystyle {\hat {H}}_{n}={\begin{bmatrix}\hbar \omega _{c}(n-1)+\hbar \omega _{eg}&\langle e,n-1|{\hat {H}}_{JC}|g,n\rangle \\\langle g,n|{\hat {H}}_{JC}|e,n-1\rangle &n\hbar \omega _{c}\\\end{bmatrix}}}

ここで、次の関係を使用します。 g,n|H^JC|e,n1=gcg,n|a^cσ^|e,n1+gcg,n|a^cσ^+|e,n1=ngc{\displaystyle \langle g,n|{\hat {H}}_{JC}|e,n-1\rangle =\hbar g_{c}\langle g,n|{\hat {a}}_{c}^{\dagger }{\hat {\sigma }}_{-}|e,n-1\rangle +\hbar g_{c}\langle g,n|{\hat {a}}_{c}{\hat {\sigma }}_{+}|e,n-1\rangle ={\sqrt {n}}\hbar g_{c}}

n番目の部分空間で作用するハミルトニアンの部分は次のようになります。 H^n=[nωcΔnΩ2nΩ2nωc]{\displaystyle {\hat {H}}_{n}={\begin{bmatrix}n\hbar \omega _{c}-\hbar \Delta &{\frac {{\sqrt {n}}\hbar \Omega }{2}}\\{\frac {{\sqrt {n}}\hbar \Omega }{2}}&n\hbar \omega _{c}\\\end{bmatrix}}}

エネルギーを から に の量だけシフトすると、[ 30 ]得られる。|e{\displaystyle |e\rangle }|g{\displaystyle |g\rangle }12Δ{\displaystyle {\frac {1}{2}}\hbar \Delta }H^n=[nωc12ΔnΩ2nΩ2nωc+12Δ]=nωcI^(n)Δ2σ^z(n)+12nΩσ^x(n){\displaystyle {\hat {H}}_{n}={\begin{bmatrix}n\hbar \omega _{c}-{\frac {1}{2}}\hbar \Delta &{\frac {{\sqrt {n}}\hbar \Omega }{2}}\\{\frac {{\sqrt {n}}\hbar \Omega }{2}}&n\hbar \omega _{c}+{\frac {1}{2}}\hbar \Delta \\\end{bmatrix}}=n\hbar \omega _{c}{\hat {I}}^{(n)}-{\frac {\hbar \Delta }{2}}{\hat {\sigma }}_{z}^{(n)}+{\frac {1}{2}}{\sqrt {n}}\hbar \Omega {\hat {\sigma }}_{x}^{(n)}}

ここで、 は系のラビ周波数、 は共振と原子遷移の周波数間の いわゆる「デチューニング」である。また、以下の演算子も定義した。2gc=Ω{\displaystyle 2g_{c}=\Omega }Δ=ωcωeg{\displaystyle \Delta =\omega _{c}-\omega _{eg}}I^(n)=|e,n1e,n1|+|g,ng,n|σ^z(n)=|e,n1e,n1||g,ng,n|σ^x(n)=|e,n1g,n|+|g,ne,n1|.{\displaystyle {\begin{aligned}{\hat {I}}^{(n)}&=\left|e,n-1\right\rangle \left\langle e,n-1\right|+\left|g,n\right\rangle \left\langle g,n\right|\\[1ex]{\hat {\sigma }}_{z}^{(n)}&=\left|e,n-1\right\rangle \left\langle e,n-1\right|-\left|g,n\right\rangle \left\langle g,n\right|\\[1ex]{\hat {\sigma }}_{x}^{(n)}&=\left|e,n-1\right\rangle \left\langle g,n\right|+\left|g,n\right\rangle \left\langle e,n-1\right|.\\[-1ex]\,\end{aligned}}}

原子場系のn番目のエネルギー準位のヒルベルト空間における恒等演算子とパウリのx, z演算子である。この単純なハミルトニアンは、ラビ問題で見られるものと同じ形である。対角化により、エネルギー固有値固有状態は以下のようになる。[ 30 ] [ 31 ] ここで、角度は関係式で定義される。 2×2{\displaystyle 2\times 2}En,±=(nωc12Δ)±12Δ2+nΩ2|n,+=cos(θn2)|e,n1+sin(θn2)|g,n|n,=cos(θn2)|g,nsin(θn2)|e,n1{\displaystyle {\begin{aligned}E_{n,\pm }&=\left(n\hbar \omega _{c}-{\frac {1}{2}}\hbar \Delta \right)\pm {\frac {1}{2}}\hbar {\sqrt {\Delta ^{2}+n\Omega ^{2}}}\\|n,+\rangle &=\cos \left({\frac {\theta _{n}}{2}}\right)|e,n-1\rangle +\sin \left({\frac {\theta _{n}}{2}}\right)|g,n\rangle \\|n,-\rangle &=\cos \left({\frac {\theta _{n}}{2}}\right)|g,n\rangle -\sin \left({\frac {\theta _{n}}{2}}\right)|e,n-1\rangle \\\end{aligned}}}θn{\displaystyle \theta _{n}}tanθn=nΩΔ{\displaystyle \tan \theta _{n}=-{\frac {{\sqrt {n}}\Omega }{\Delta }}}

真空ラビ振動

原子が励起状態で共振器に進入し、共振器が真空状態にある場合を考える。さらに、モードの角周波数は、 を含む原子遷移周波数に近似できると仮定する。このとき、原子場系の状態は時間の関数として次のように表される。 Δ0{\displaystyle \Delta \approx 0}|ψ(t)=cos(Ωt2)|e,0isin(Ωt2)|g,1{\displaystyle |\psi (t)\rangle =\cos \left({\frac {\Omega t}{2}}\right)|e,0\rangle -i\sin \left({\frac {\Omega t}{2}}\right)|g,1\rangle }

したがって、一定時間空洞と相互作用した後にシステムが基底状態または励起状態にある確率は次のようになる:[ 32 ]t{\displaystyle t}Pe(t)=|e,0|ψ(t)|2=cos2(Ωt2)Pg(t)=|g,1|ψ(t)|2=sin2(Ωt2){\displaystyle {\begin{aligned}P_{e}(t)&=|\langle e,0|\psi (t)\rangle |^{2}=\cos ^{2}\left({\frac {\Omega t}{2}}\right)\\P_{g}(t)&=|\langle g,1|\psi (t)\rangle |^{2}=\sin ^{2}\left({\frac {\Omega t}{2}}\right)\\\end{aligned}}}

したがって、原子がどちらの状態にあるかの確率振幅は振動する。これは真空ラビ振動現象の量子力学的説明である。この場合、原子場系には最初に励起された原子によって運ばれた単一の量子しか存在しなかった。一般に、量子の原子場系に関連するラビ振動は周波数を持つ。後述するように、この離散的な周波数スペクトルが、モデルにおける崩壊とそれに続く再生の確率の根本的な理由である。 n{\displaystyle n}Ωn=nΩ2{\displaystyle \Omega _{n}={\frac {{\sqrt {n}}\Omega }{2}}}

ジェインズ・カミングス・ラダー

前の節で示したように、原子空洞系の初期状態がまたは である場合(最初は特定の状態(基底状態または励起状態)にある原子が、既知の数の光子を含む空洞に入る場合)、その後の原子空洞系の状態は、原子空洞系の 新しい固有状態の重ね合わせになります。|e,n1{\displaystyle |e,n-1\rangle }|g,n{\displaystyle |g,n\rangle }|n,+=cos(θn2)|e,n1+sin(θn2)|g,n|n,=cos(θn2)|g,nsin(θn2)|e,n1{\displaystyle {\begin{aligned}|n,+\rangle &=\cos \left({\frac {\theta _{n}}{2}}\right)|e,n-1\rangle +\sin \left({\frac {\theta _{n}}{2}}\right)|g,n\rangle \\|n,-\rangle &=\cos \left({\frac {\theta _{n}}{2}}\right)|g,n\rangle -\sin \left({\frac {\theta _{n}}{2}}\right)|e,n-1\rangle \\\end{aligned}}}

原子-場相互作用によって引き起こされるハミルトニアンの変化によるこの固有状態の変化は、原子の「ドレッシング」と呼ばれることもあり、新しい固有状態はドレッシング状態と呼ばれます。[ 30 ] ドレッシング状態間のエネルギー差は: 特に興味深いのは、空洞周波数が原子の遷移周波数と完全に共鳴する場合です。したがって、共鳴の場合、ドレッシング状態は:[ 31 ]δE=E+E=Δ2+nΩ2{\displaystyle \delta E=E_{+}-E_{-}=\hbar {\sqrt {\Delta ^{2}+n\Omega ^{2}}}}ωeg=ωcΔ=0{\displaystyle \omega _{eg}=\omega _{c}\implies \Delta =0}|n,±=12(|g,n|e,n1){\displaystyle |n,\pm \rangle ={\frac {1}{\sqrt {2}}}\left(|g,n\rangle \mp |e,n-1\rangle \right)}

エネルギー差 の場合。したがって、原子と場の相互作用により、状態と状態の縮退はによって 分割されます。このエネルギー準位の非線形階層は、 にスケーリングされ、ジェインズ・カミングス・ラダーとして知られています。この非線形分裂効果は純粋に量子力学的なものであり、いかなる半古典的モデルでも説明できません。[ 20 ]δE=nΩ{\displaystyle \delta E={\sqrt {n}}\hbar \Omega }|e,n1{\displaystyle |e,n-1\rangle }|g,n{\displaystyle |g,n\rangle }nΩ{\displaystyle {\sqrt {n}}\hbar \Omega }n{\displaystyle {\sqrt {n}}}

確率の崩壊と復活

最初は基底状態にあった原子が、最初はコヒーレント状態に準備された場のモードと相互作用すると仮定します。この場合、原子場システムの初期状態は次のようになります。 |ψ(0)=|g,α=n=0e|α|2/2αnn!|g,n{\displaystyle |\psi (0)\rangle =|g,\alpha \rangle =\sum _{n=0}^{\infty }e^{-|\alpha |^{2}/2}{\frac {\alpha ^{n}}{\sqrt {n!}}}|g,n\rangle }

簡単にするために共鳴ケース( )を取ると、n番目の数の部分空間のハミルトニアンは次のようになります。 Δ=0{\displaystyle \Delta =0}H^n=(n+12)I^(n)+nΩ2σ^x(n){\displaystyle {\hat {H}}_{n}=\left(n+{\frac {1}{2}}\right){\hat {I}}^{(n)}+{\frac {\hbar {\sqrt {n}}\Omega }{2}}{\hat {\sigma }}_{x}^{(n)}}

これを用いると、原子場システムの時間発展は次のようになります。 定数因子もも、零点エネルギーを表すため、全体的な位相を超えてダイナミクスに寄与しないことに注意してください。この場合、原子が後の時点で励起状態に反転している確率は次のようになります。 ここで、 はコヒーレント状態における平均光子数であることを確認しました。平均光子数が大きい場合、コヒーレント状態の統計はポアソン分布に従うため、分散と平均の比は となります。この結果を用いて、を の最小の非ゼロ次数まで展開すると、次のようになります。 これを和に代入すると、複雑な指数積が得られます。 |ψ(t)=eiH^nt/|ψ(0)=e|α|2/2|g,0+n=1e|α|2/2αnn!einωct(cos(nΩt/2)I^(n)isin(nΩt/2)σ^x(n))|g,n=e|α|2/2|g,0+n=1e|α|2/2αnn!einωct(cos(nΩt/2)|g,nisin(nΩt/2)|e,n1){\displaystyle {\begin{aligned}|\psi (t)\rangle &=e^{-i{\hat {H}}_{n}t/\hbar }|\psi (0)\rangle \\&=e^{-|\alpha |^{2}/2}|g,0\rangle +\sum _{n=1}^{\infty }e^{-|\alpha |^{2}/2}{\frac {\alpha ^{n}}{\sqrt {n!}}}e^{-in\omega _{c}t}\left(\cos {({\sqrt {n}}\Omega t/2)}{\hat {I}}^{(n)}-i\sin {({\sqrt {n}}\Omega t/2)}{\hat {\sigma }}_{x}^{(n)}\right)|g,n\rangle \\&=e^{-|\alpha |^{2}/2}|g,0\rangle +\sum _{n=1}^{\infty }e^{-|\alpha |^{2}/2}{\frac {\alpha ^{n}}{\sqrt {n!}}}e^{-in\omega _{c}t}\left(\cos {({\sqrt {n}}\Omega t/2)}|g,n\rangle -i\sin {({\sqrt {n}}\Omega t/2)}|e,n-1\rangle \right)\end{aligned}}}ωc2I^(n){\displaystyle {\frac {\hbar \omega _{c}}{2}}{\hat {I}}^{(n)}}H^0{\displaystyle {\hat {H}}_{0}}t{\displaystyle t}Pe(t)=|e|ψ(t)|2=n=1e|α|2n!|α|2nsin2(12nΩt)=n=1ennnn!sin2(12nΩt)=n=1ennnn!sin2(Ωnt){\displaystyle {\begin{aligned}P_{e}(t)=\left|\langle e|\psi (t)\rangle \right|^{2}&=\sum _{n=1}^{\infty }{\frac {e^{-|\alpha |^{2}}}{n!}}|\alpha |^{2n}\sin ^{2}\left({\tfrac {1}{2}}{\sqrt {n}}\Omega t\right)\\[2ex]&=\sum _{n=1}^{\infty }{\frac {e^{-\langle n\rangle }\langle n\rangle ^{n}}{n!}}\sin ^{2}\left({\tfrac {1}{2}}{\sqrt {n}}\Omega t\right)\\[2ex]&=\sum _{n=1}^{\infty }{\frac {e^{-\langle n\rangle }\langle n\rangle ^{n}}{n!}}\sin ^{2}(\Omega _{n}t)\\{}\end{aligned}}}n=|α|2{\displaystyle \langle n\rangle =|\alpha |^{2}}(Δn)2/n21/n{\displaystyle \langle (\Delta n)^{2}\rangle /\langle n\rangle ^{2}\simeq 1/\langle n\rangle }Ωn{\displaystyle \Omega _{n}}n{\displaystyle \langle n\rangle }n{\displaystyle n}ΩnΩ2n(1+12nnn){\displaystyle \Omega _{n}\simeq {\frac {\Omega }{2}}{\sqrt {\langle n\rangle }}\left(1+{\frac {1}{2}}{\frac {n-\langle n\rangle }{\langle n\rangle }}\right)}Pe(t)12en4(einΩt/2exp[nexp(iΩt2n)]+einΩt/2exp[nexp(iΩt2n)]){\displaystyle P_{e}(t)\simeq {\frac {1}{2}}-{\frac {e^{-\langle n\rangle }}{4}}\cdot \left(e^{-i{\sqrt {\langle n\rangle }}\Omega t/2}\exp \left[\langle n\rangle \exp \left(-{\frac {i\Omega t}{2{\sqrt {\langle n\rangle }}}}\right)\right]+e^{i{\sqrt {\langle n\rangle }}\Omega t/2}\exp \left[\langle n\rangle \exp \left({\frac {i\Omega t}{2{\sqrt {\langle n\rangle }}}}\right)\right]\right)}

平均光子数 の系において、系が励起状態にある確率を単位のないパラメータの関数としてプロットしたグラフ。短時間で初期崩壊が起こり、その後長時間で回復する様子に注目してください。この挙動は、場の量子化によって生じる周波数の離散スペクトルに起因します。gt{\displaystyle gt}n=25{\displaystyle \langle n\rangle =25}

となるような「小さい」時間の場合、最後の項の二重指数関数内の内部指数関数は2次まで展開して次式を得ることができます。 Ωt2n{\displaystyle {\frac {\Omega t}{2}}\ll {\sqrt {\langle n\rangle }}}

Pe(t)1212cos[nΩt]eΩ2t2/8{\displaystyle P_{e}(t)\simeq {\frac {1}{2}}-{\frac {1}{2}}\cdot \cos \left[{\sqrt {\langle n\rangle }}\Omega t\right]e^{-\Omega ^{2}t^{2}/8}}

この結果は、励起状態の占有確率が有効周波数で振動することを示している。また、それが特性時間とともに減衰することも示している:[ 6 ] [ 7 ] [ 31 ]Ωeff=nΩ{\textstyle \Omega _{\text{eff}}={\sqrt {\langle n\rangle }}\Omega }τc=2Ω{\displaystyle \tau _{c}={\frac {\sqrt {2}}{\Omega }}}

この崩壊は、異なる周波数成分が位相をずらし、時間の経過とともに破壊的に干渉し始めるため、それらの成分間の破壊的干渉の結果として簡単に理解できます。[ 31 ] [ 32 ]しかし、周波数が離散的なスペクトルを持つという事実は、より長い時間領域で別の興味深い結果につながります。その場合、ゆっくりと変化する二重指数関数の周期的な性質は、時間で確率の復活もあるはずであると予測します。 τr=4πΩn.{\displaystyle \tau _{r}={\frac {4\pi }{\Omega }}{\sqrt {\langle n\rangle }}.}

確率の復活は、様々な離散周波数の再位相によるものである。もし場が古典的なものであれば、周波数は連続スペクトルを持ち、そのような再位相は有限時間内には決して起こり得ない。[ 7 ] [ 31 ] [ 32 ]

右は、コヒーレント状態で準備された共振器と相互作用した後に、基底状態にあった原子が励起状態へ遷移する確率を、単位のないパラメータに対してプロットしたものです。初期の崩壊と、それに続く長時間にわたる明確な回復に注目してください。 gt=Ωt/2{\displaystyle gt=\Omega t/2}

量子振動の崩壊と復活

ColRev3a40

場が最初にコヒーレント状態に準備されると、ジェインズ・カミングス模型は原子反転の崩壊と復活を予測する。当初、コヒーレントラビ振動は異なる光子数成分間の干渉により減衰するが、後に周期的に復活する。この非古典的な効果は半古典的な放射理論には見当たらず、この模型の最も特徴的な予測の一つとされている。[ 33 ] [ 34 ]

この原子反転の量子振動のグラフは、二次スケールのデチューニングパラメータ (デチューニングパラメータ )について、AA KaratsubaとEA Karatsubaによって得られた公式に基づいて作成されたものである。[ 35 ]a=(δ/2g)2=40{\displaystyle a=(\delta /2g)^{2}=40}δ{\displaystyle \delta }

参照

参考文献

  1. ^ Greentree, Andrew D; Koch, Jens; Larson, Jonas (2013年11月1日). 「ジェインズ・カミングス物理学の50年」. Journal of Physics B: Atomic, Molecular and Optical Physics . 46 (22) 220201. doi : 10.1088/0953-4075/46/22/220201 .
  2. ^ Larson, Jonas; Mavrogordatos, Themistoklis; Parkins, Scott; Vidiella-Barranco, Antonio (2024年8月1日). 「ジェインズ・カミングスモデル:60年を経てなお進化を続ける」アメリカ光学会誌 B . 41 (8): JCM1. doi : 10.1364/JOSAB.536847 .
  3. ^ ET Jaynes; FW Cummings (1963). 「ビームメーザーへの応用による量子および半古典的放射理論の比較」Proc. IEEE . 51 (1): 89– 109. doi : 10.1109/PROC.1963.1664 .
  4. ^ Frederick W. Cummings (2013). 「1950年代のスタンフォード大学におけるE.T. Jaynesとの卒業研究を回想する」 . Journal of Physics B: Atomic, Molecular and Optical Physics . 46 (22): 220202 (3pp). Bibcode : 2013JPhB...46v0202C . doi : 10.1088/0953-4075/46/22/220202 .
  5. ^ a b Zhu, Yifu、Gauthier, Daniel J.、Morin, SE、Wu, Qilin、Carmichael, HJ、Mossberg, TW (1990). 「線形分散理論の特徴としての真空ラビ分裂:解析と実験的観察」. Phys. Rev. Lett . 64 (21). American Physical Society: 2499– 2502. Bibcode : 1990PhRvL..64.2499Z . doi : 10.1103/PhysRevLett.64.2499 . PMID 10041728 . {{cite journal}}: CS1 maint: multiple names: authors list (link)
  6. ^ a b F.W. Cummings (1965). 「単一モードにおける誘導放出」. Phys. Rev. 140 ( 4A): A1051– A1056. Bibcode : 1965PhRv..140.1051C . doi : 10.1103/PhysRev.140.A1051 . S2CID 117003321 . 
  7. ^ a b c J.H. Eberly; NB Narozhny; JJ Sanchez-Mondragon (1980). 「単純な量子モデルにおける周期的な自発的崩壊と再生」. Phys. Rev. Lett . 44 (20): 1323– 1326. Bibcode : 1980PhRvL..44.1323E . doi : 10.1103/PhysRevLett.44.1323 .
  8. ^ S. Haroche; JM Raimond (1985). 「共鳴空洞におけるリュードベリ状態の放射特性」.原子・分子物理学の進歩. 20 : 347–411 . Bibcode : 1985AdAMP..20..347H . doi : 10.1016/S0065-2199(08)60271-7 . ISBN 978-0-12-003820-6
  9. ^ JAC Gallas; G. Leuchs; H. Walther; H. Figger (1985). 「リュードベリ原子:高解像度分光法と放射線相互作用-リュードベリ分子」.原子・分子物理学の進歩. 20 : 413–466 . Bibcode : 1985AdAMP..20..413G . doi : 10.1016/S0065-2199(08)60272-9 . ISBN 978-0-12-003820-6
  10. ^ G. Rempe; H. Walther; N. Klein (1987). 「一原子メーザーにおける量子崩壊と再生の観測」. Phys. Rev. Lett . 58 (4): 353– 356. Bibcode : 1987PhRvL..58..353R . doi : 10.1103/PhysRevLett.58.353 . PMID 10034912 . 
  11. ^ S. Haroche; JM Raimond (1985). 「共鳴空洞におけるリュードベリ状態の放射特性」 D. Bates; B. Bederson (編) 『原子・分子物理学の進歩』第20巻 ニューヨーク: アカデミック・プレス p. 350.
  12. ^ JAC Gallas; G. Leuchs; H. Walther; H. Figger (1985). 「リュードベリ原子:高解像度分光法と放射線相互作用-リュードベリ分子」 D. Bates; B. Bederson (編). 『原子・分子物理学の進歩』 第20巻. ニューヨーク: Academic Press. p. 414.
  13. ^ a b S.E. Morin; CC Yu; TW Mossberg (1994). 「大規模体積における強い原子-空洞結合と空洞内原子線幅の観測」. Phys. Rev. Lett . 73 (11): 1489– 1492. Bibcode : 1994PhRvL..73.1489M . doi : 10.1103/PhysRevLett.73.1489 . PMID 10056806 . 
  14. ^吉江 孝文; 他 (2004). 「フォトニック結晶ナノ共振器における単一量子ドットによる真空ラビ分裂」. Nature . 432 ( 7014): 200– 203. Bibcode : 2004Natur.432..200Y . doi : 10.1038/nature03119 . PMID 15538363. S2CID 4391322 .  
  15. ^ Florian Meier; David D. Awschalom (2004). 「2モード共振器における量子ドットのスピン光子ダイナミクス」 . Phys. Rev. B. 70 ( 20) 205329. arXiv : cond-mat/0405342 . Bibcode : 2004PhRvB..70t5329M . doi : 10.1103/PhysRevB.70.205329 . S2CID 119454297 . 
  16. ^ T. Yoshie; et al. (2004). 「フォトニック結晶ナノ共振器における単一量子ドットによる真空ラビ分裂」Nature . 432 (7014): 200– 203. Bibcode : 2004Natur.432..200Y . doi : 10.1038 / nature03119 . PMID 15538363. S2CID 4391322 .  
  17. ^ J. Kasprzak; et al. (2010). 「量子ドット/マイクロキャビティシステムのジェインズ・カミングス・ラダーの上へ」. Nature Materials . 9 (4): 304– 308. Bibcode : 2010NatMa...9..304K . doi : 10.1038/nmat2717 . PMID 20208523 . 
  18. ^ J. Basset; et al. (2013). 「単一電子二重量子ドットの双極子結合による単一光子モード」 . Phys. Rev. B. 88 ( 12) 125312. arXiv : 1304.5141 . Bibcode : 2013PhRvB..88l5312B . doi : 10.1103/PhysRevB.88.125312 . S2CID 56422948 . 
  19. ^ JM Fink; et al. (2008). 「ジェインズ・カミングス・ラダーを登り、空洞QEDシステムにおけるその非線形性を観察する」Nature . 454 (7202): 315– 318. arXiv : 0902.1827 . Bibcode : 2008Natur.454..315F . doi : 10.1038/nature07112 . PMID 18633413 . S2CID 205214018 .  
  20. ^ a b c J. Lee; et al. (2017). 「リュードベリドレス原子を用いたジェインズ・カミングスラダーの実証」 . Phys. Rev. A. 95 ( 4) 041801. arXiv : 1609.03940 . Bibcode : 2017PhRvA..95d1801L . doi : 10.1103/PhysRevA.95.041801 . S2CID 118360147 . 
  21. ^ Kukliński, J.; Madajczyk, J. (1988). 「ジェインズ・カミングス模型における強いスクイージング」. Physical Review A. 37 ( 8): 3175– 3178. Bibcode : 1988PhRvA..37.3175K . doi : 10.1103/PhysRevA.37.3175 . PMID 9900058 . 
  22. ^ Gea-Banacloche, J. (1993). 「準古典場を伴うジェインズ・カミングス模型:散逸の影響」. Physical Review A. 47 ( 3): 2221– 2234. Bibcode : 1993PhRvA..47.2221G . doi : 10.1103/PhysRevA.47.2221 . PMID 9909174 . 
  23. ^ Rodríguez-Lara, B.; Moya-Cessa, H.; Klimov, A. (2005). 「レーザー駆動トラップイオンシステムにおけるJaynes-Cummingsダイナミクスと反Jaynes-Cummingsダイナミクスの結合」. Physical Review A. 71 ( 2) 023811. Bibcode : 2005PhRvA..71b3811R . doi : 10.1103/PhysRevA.71.023811 .
  24. ^ P. Kochanski; Z. Bialynicka-Birula; I. Bialynicki-Birula (2001). 「トロヤ状態にある電子による空洞内の電磁場の圧縮」. Phys. Rev. A . 63 (1): 013811–013811–8. arXiv : quant-ph/0007033 . Bibcode : 2000PhRvA..63a3811K . doi : 10.1103/PhysRevA.63.013811 . S2CID 36895794 . 
  25. ^ Kundu, A. (2005). 「回転波近似を用いた場合と用いない場合の量子可積分多原子物質放射モデル」.理論数理物理学. 144 (1): 975– 984. arXiv : nlin/0409032 . Bibcode : 2005TMP...144..975K . doi : 10.1007/s11232-005-0125-7 . S2CID 119379204 . 
  26. ^ Hussin, V.; Nieto, LM (2005). 「回転波近似におけるJaynes-Cummings模型のラダー演算子とコヒーレント状態」. Journal of Mathematical Physics . 46 (12): 122102. Bibcode : 2005JMP....46l2102H . doi : 10.1063/1.2137718 .
  27. ^ Shore, Bruce W.; Knight, Peter L. (1993). 「ジェインズ・カミングスモデル」. J​​ournal of Modern Optics . 40 (7): 1195– 1238. Bibcode : 1993JMOp...40.1195S . doi : 10.1080/09500349314551321 .
  28. ^ Ellinas, D.; Smyrnakis, I. (2005). 「光共振器によって駆動される量子ランダムウォークの漸近解析」. J. Opt. B. 7 ( 7): S152. arXiv : quant-ph/0510112 . Bibcode : 2005JOptB...7S.152E . doi : 10.1088/1464-4266/7/7/004 . S2CID 15489395 . 
  29. ^ S. Stenholm、「原子および分子と相互作用する電磁場の量子論」、 Physics Reports 6 (1)、1–121 (1973)。
  30. ^ a b c d e fバーマン, ポール・R.; マリノフスキー, ウラジミール・S. (2011).レーザー分光法と量子光学の原理. プリンストン大学出版局. pp.  360– 370. ISBN 978-0-691-14056-8
  31. ^ a b c d e Vogel、W.;ウェルシュ、DG. (2006)。量子光学(第 3 版)。ワイリー-VCH。ページ 408–423。ISBN 978-3-527-40507-7
  32. ^ a b cピエール・メイストル、サージェントIII、マレー (2007). 『量子光学の要素』(第4版). シュプリンガー. pp.  333– 338. ISBN 978-3-540-74209-8
  33. ^ Jaynes, E. T; Cummings, F. W (1963).量子および半古典的放射理論の比較とビームメーザーへの応用”. Proceedings of the IEEE .
  34. ^ Rempe, G; Walther, H. 1原子メーザーにおける量子崩壊と再生の観測。Physical Review Letters
  35. ^ AA Karatsuba; EA Karatsuba (2009). 「ジェインズ・カミングス模型における崩壊と再生の再総和式」. J. Phys. A: Math. Theor . 42 (19): 195304, 16. Bibcode : 2009JPhA...42s5304K . doi : 10.1088/1751-8113/42/19/195304 . S2CID 120269208 . 

さらに読む

  • バーマン, PR; マリオフスキー, VS (2011). 『レーザー分光法と量子光学の原理』プリンストン大学出版局. ISBN 978-0-691-14056-8
  • ジェリー, CC; ナイト, PL (2005). 『量子光学入門』 ケンブリッジ: ケンブリッジ大学出版局. ISBN 0-521-52735-X
  • Scully, MO; Zubairy, MS (1997).量子光学. ケンブリッジ: ケンブリッジ大学出版局. ISBN 0-521-43595-1
  • フォーゲル、W.ウェルシュ、DG (2006)。量子光学(第 3 版)。ワイリー-VCH。ISBN 978-3-527-40507-7
  • ウォールズ, DF; ミルバーン, GJ (1995).量子光学. シュプリンガー・フェアラーク. ISBN 3-540-58831-0