ウィグナー関数

ウィグナー関数（ウィグナーかんすう、）とは、ユージン・ウィグナーにより1932年に導入された、古典統計力学を量子補正するための関数である。その目標は、シュレーディンガー方程式に表われる波動関数を位相空間上の確率分布と結びつけることであった。ウィグナーの擬確率分布関数（）、ウィグナー・ビレ分布 () とも。

ウィグナー関数は量子力学的波動関数のすべての空間的自己相関の母関数である。従って、ウィグナー関数と密度行列との間の写像により、実位相空間上の関数とヘルマン・ワイルが1927年に導入した; ; エルミート演算子とを表現論的な文脈で対応づけられる()。ウィグナー関数は密度行列をしたものとみなすことができ、よって密度行列の位相空間上での表現とみなせる。1948年、によって独立にスペクトログラムの一種、信号エネルギーの局所時間・周波数表示方法として再導入された。

1949年、は量子化された運動量の母関数としてウィグナー関数を再導入し、これを用いて全ての量子期待値を計算する方法を確立し、位相空間上における量子力学の基礎を築いた（を参照）。統計力学、量子化学、量子光学、古典光学、および電子工学、地震学、音楽の時間周波数解析、生物学のスペクトログラム、音声合成、エンジンの設計など、信号処理を伴う幅広い分野で応用されている。

古典力学との関係古典力学的には、粒子は決まった位置と運動量を持ち、その運動状態は位相空間上の一点により表現される。多数の粒子の集合体が与えられたとき、位相空間内の特定の領域に粒子をみいだす確率はリウビル確率密度と呼ばれる確率密度関数に従う。しかし、このような決定論的な取扱いは量子力学的な粒子に対しては不確定性原理のために不可能である。ウィグナー関数は古典的な確率密度分布と同様に取り扱うことができるが、ウィグナー関数は古典的な確率密度関数の満すべき条件を全て満たしてはいない。そのかわり、古典的な分布が必ずしも満たさない有界性を満たしている。

たとえば、ウィグナー関数は古典的分布ではありえない負値をとることがよくある。そして、ウィグナー関数が負値をとることは量子干渉が起きていることを示す指標である。ウィグナー関数にディラック定数よりも小さな位相空間体積における構造を無視するような処理(たとえば伏見表示（後述)を得るために位相空間上のガウス関数で畳み込むなど）を施すと、半正定値関数となり半古典形式に粗視化できるこの畳み込みは可逆であるため、情報は全く失われておらず、量子エントロピーは増加していない。。

負の値をとる領域が存在しても、（幅の小さいガウス関数と畳み込んだ場合）多くの場合その領域は「小さく」なる。つまり、その領域はの数倍より大きくなることはなく、そのためにおいては消滅する。これは、位相空間上でよりも小さな体積をもつ領域に粒子の運動状態を特定することはできないとする不確定性原理による遮蔽であり、「負の確率」という概念の矛盾を軽減している。

定義と意味を波動関数とし、をそれぞれ位置および運動量、または他の正準共役量（例えば電磁場の実部および虚部、もしくは信号における時間と周波数など)とすると、ウィグナー関数は以下のように定義される。

ここで、ウィグナー関数はが上に台を持たない領域でも台を持つことがある。

ウィグナー関数の定義はおよびについて対称である。

ここで、はのフーリエ変換である。

三次元系では、以下のようになる。

混合状態を含む一般の場合には、密度行列のウィグナー変換を用いて以下のように定義される。

ここで、である。このウィグナー変換（または写像）は、位相空間上の関数をヒルベルト空間上の作用素へと移すの逆になっている。

よって、ウィグナー関数は位相空間上の量子力学における基礎となっている。

1949年、ホセ・エンリケ・モヤルはウィグナー関数が確率密度関数と同様に、位相空間に測度を与えていることを明らかにした。つまり、古典確率論と同様に、c-数を返す位相空間上の一価の関数とワイル変換によって関係づけられる作用素の期待値をウィグナー関数を使って定義することができる（後述のウィグナー・ワイル変換の性質を参照）。

具体的に書き下せば、作用素の期待値は作用素をウィグナー変換して得られる関数の「位相空間上の平均値」として以下のように定義される。

数学的特徴 thumb|図1: それぞれ a) 真空、 b) のフォック状態 (例：単一光子)、 c) のフォック状態、のウィグナー関数。は実関数である。およびの確率密度関数は次の周辺確率により与えられる。系が純粋状態ならば系が純粋状態ならば

通常密度行列のトレースは1である。は次の鏡映対称性をもつ。時間対称性: 空間対称性: はガリレイ共変（ガリレイ変換に対して不変）である。

ローレンツ共変ではない。位相空間上の各点における運動方程式は力のない古典力学の方程式である。実際、調和力が働いている場合も古典的である。状態の重なり積分は以下のように計算される。

作用素の期待値（平均値）はウィグナー変換したのちに位相空間上の平均値をとることにより与えられる。

が物理的な（正の）密度行列を持つためには、全ての純粋状態に対して以下を満たす必要がある。

コーシー・シュワルツ不等式を用い、純粋状態においては以下のように有界である。

古典極限においては非有界となる。このことから、は座標空間においては確率密度関数に帰着し、通常は非常に局在化した、運動量方向にデルタ関数のかかった分布になる。つまり、古典極限は「尖って」いる。このことから、この有界性は不確定性原理を反映し、ウィグナー関数が位相空間上で完全に局在化した関数になることを防いできると言える; 。

ウィグナー関数の時間発展方程式

right|thumb|220x220px|図2: 単純な調和振動子の、位相空間の原点からずれた基底状態（コヒーレント状態）を表すウィグナー関数。（クリックして動画を表示）古典力学の場合と同じように、剛体回転している。これは単純な調和振動子に特異な性質である。一般教育ウェブサイトより。ウィグナー変換は、ヒルベルト空間上の作用素を位相空間上の関数へと写す可逆な変換であり、以下のように定義される。

エルミート演算子は実関数に写される。位相空間からヒルベルト空間への逆変換はワイル変換と呼ばれる。

（別の定義のワイル変換も存在することに注意。）

thumb|230x230px|図3: 量子フリップフロップ: 単純調和振動子の基底状態と第一励起状態の重ね合わせ状態を表すウィグナー関数。位相空間での剛体回転は座標空間での振動を表わす。（クリックして動画を表示）この項で取り扱ってきたウィグナー関数は、密度行列をウィグナー変換したものと捉えることができる。よって、ある作用素と密度行列をかけたもののトレースは、その作用素をウィグナー変換したものと、ウィグナー関数との位相空間上の重なり積分と等しい。

シュレーディンガー描像における密度行列の時間発展を記述する、フォン・ノイマン方程式のウィグナー変換はウィグナー関数に対するモヤル方程式

に帰着する。ここで、はハミルトニアン、はを表わす。古典極限では、モヤル括弧はポアソン括弧に帰着し、従ってこの時間発展方程式は古典統計力学におけるリウビル方程式に帰着する。

Quantum characteristicsの記法を用いて、上の方程式の形式的な厳密解は以下のように書ける。ここでとはいわゆる量子ハミルトン方程式の解で、初期条件及びに従い、積の合成は全ての関数について成り立つものとする。合成は完全に非局所（モヤルが指摘したように、「量子確率流体」は拡散する）であるため、通常はウィグナー関数の発展につれて、局所的な軌道のなごりはほとんど確認できなくなるQuantum characteristics とファインマンの経路積分やド・ブロイ-ボーム理論におけるトラジェクトリを混同しないように注意。この三つ巴の曖昧さから、ニールス・ボーアの立場をより理解できるだろう。彼は熱心に原子物理学におけるトラジェクトリへの言及に反対していた。例えば、1948年のポコノ会議において彼はリチャード・ファインマンに対して「…原子内の電子のトラジェクトリについて議論することはできない。なぜなら観測できないからである」 ("The Beat of a Different Drum: The Life and Science of Richard Feynman", by Jagdish Mehra (Oxford, 1994, pp. 245-248))と言っている。このような議論はエルンスト・マッハがかつて原子論を批判した際や、1960年代にやトゥーリオ・レッジェらが局所的量子場理論をS行列で置き換えようとした際などに広く用いられた論法である。今日では、完全に原子論的概念に基いた統計物理学が標準的に教えられているし、S行列理論は時代遅れになっているのに対してファインマンの経路積分法はゲージ理論において最も効率的な手法であるとみなされている。。積の積分表示においては、積を連続的に位相空間経路積分に適用することで、このウィグナー関数の発展方程式を解くことができる（以下も参照特に86-89ページの5節 "Path integral for the propagator" を参照。）。

古典極限ウィグナー関数により古典極限を記述することで、位相空間上の古典動力学と量子動力学とを対応づけることができる。

近年、ウィグナー関数法は1932年にとフォン・ノイマンによって導入された、古典力学の演算子表式の量子的アナロジーになっていることが示唆されている。の極限では、ウィグナー関数の時間発展はクープマン・フォンノイマン波動関数の時間発展に漸近するDenys I. Bondar, Renan Cabrera, Dmitry V. Zhdanov, Herschel A. Rabitz: Wigner Function's Negativity Demystified (submitted February 2012, version of 3 November 2012)Renan Cabrera, Denys I. Bondar, Herschel A. Rabitz: Relativistic Wigner function and consistent classical limit for spin 1/2 particles, (submitted on 26 July 2011, version of 22 August 2012)。

ウィグナー関数と他の量子力学の表現との関係ウィグナー関数は、ドブロイ・ボームアンサンブルを表わす位相空間分布関数の -とみなせることが示されている (submitted 26 August 2002)。は、ウィグナー関数は位相空間上の「セル」における平均座標と平均運動量で密度行列を表わしたものと見ることができ、ドブロイ・ボーム表式はその「セル」の中心が従うダイナミクスを表わしていることを示したB. Hiley: Moyal's characteristic function, the density matrix and von Neumann's idempotent (preprint) 。

ウィグナー関数による量子状態の表現は、による量子状態の再構成と密接な関係があるF.C. Khanna, P.A. Mello, M. Revzen, Classical and Quantum Mechanical State Reconstruction, [quant-ph] (submitted December 14, 2011)。

量子力学以外でのウィグナー関数の利用

望遠鏡や光ファイバー通信機器の設計において、ウィグナー関数は単純なと波形解析とのギャップを埋めるために用いられる。ここで近軸近似の下では、はと置き換えられる。この文脈では、ウィグナー関数は干渉の影響をとりこんだまま光線の位置と角度で系を取り扱う最善の方法である。ウィグナー関数がいずれかの点で負になった場合、単純なレイトレーシングでは系をモデル化するのに不十分であることを示している。

信号解析においては、時間変化する電気信号、機械的振動、音波などをウィグナー分布により解析することがある。ここで、は時間、は角振動数（は振動数）に置き換えられる。

超高速光学において、短レーザーパルスは上と同じようにおよびで置換されたウィグナー関数により特徴づけられる。チャープ（周波数の時間依存性）などのパルスの乱れをウィグナー関数により可視化することができる。図7を参照。

量子光学では、およびを及び成分、電磁場の実部と虚部に置き換える（コヒーレント状態も参照）。図1は光の量子状態を示している。

ウィグナー関数の測定トモグラフィー

(Frequency-resolved optical gating)

ウィグナー関数

Źródło i licencja