Interpretation of quantum mechanics
確率量子力学は 、様々な外力に加え、固有のランダム過程を受ける粒子のダイナミクスを記述するための枠組みである。この枠組みは、これらの確率粒子に関連する 拡散方程式の導出を提供する。この枠組みは、ある種の保存的(またはユニタリー)拡散に対する シュレーディンガー方程式を コルモゴロフ方程式 として 導出したことで最もよく知られている 。
この導出は、作用の極値化と 量子化の 規定の組み合わせ に基づくことができる 。 この量子化規定は、 正準量子化 や 経路積分の定式化 と比較することができ、ネルソンの確率量子化または確率化と呼ばれることが多い。 この理論はシュレーディンガー方程式の導出を可能にするため、量子力学の確率的解釈を生み出した。この解釈は、確率力学理論の発展の主な動機となっている。
1930年代には エルヴィン・シュレーディンガー と ラインホルト・フュルトがともに古典的 拡散 方程式 と量子論の形式主義との類似性を認識していたが、量子力学における最初の比較的首尾一貫した確率論は1946年にハンガリーの物理学者 イムレ・フェニェス によって提唱された。 ルイ・ド・ブロイ パイロット波 から別のパイロット波に切り替えるために、量子力学の根底にある確率過程を取り入れる必要があると感じた 。 確率量子力学の理論は、この枠組みの中でシュレーディンガー方程式の導出を独立に発見した エドワード・ネルソン に帰せられる。 ゲラ 、ルッジェーロ、パヴォンら によっても発展させられた。
量子力学の確率論的解釈 確率的解釈は、量子力学 の 経路積分定式 における経路を、 確率過程 のサンプル経路として解釈する 。 量子粒子はこれらの経路のいずれかに局在するが、観測者は粒子がどこに局在するかを確実に予測することはできないと仮定する。粒子の位置を特定する唯一の方法は、測定を行うことである。観測者は、過去の測定結果と粒子に作用する力に関する知識に基づいて、そのような測定結果の確率を予測することしかできない。
この解釈は統計力学 、 特に ブラウン運動 の文脈ではよく知られている 。したがって、確率論的解釈によれば、量子力学はブラウン運動に似た方法で解釈されるべきである。 しかし、ブラウン運動の場合、統計的経路積分を定義する確率測度( ウィーナー測度 と呼ばれる)の存在は確立しており、この測度は ウィーナー過程 と呼ばれる確率過程によって生成することができる。 一方、量子力学的経路積分を定義する確率測度の存在を証明することは困難であり、 そのような確率測度が確率過程によって生成できることは保証されていない。確率力学は、量子力学の確率測度を生成するそのような確率過程の構築に関係する枠組みである。
ブラウン運動 においては 、ブラウン運動粒子の統計的変動は、粒子と多数の微小粒子との相互作用によって引き起こされることが多いことが知られている。 ウィーナー過程 によるブラウン運動の記述は 近似としてのみ用いられ、背景にある個々の粒子のダイナミクスは無視される。その代わりに、ウィーナー過程は、これらの背景粒子の影響をそれらの統計的挙動によって記述する。
量子力学の確率論的解釈は、量子粒子の量子ゆらぎの起源については不可知論的である。この解釈は、量子ゆらぎを、背景仮説と呼ばれる新たな確率論的自然法則の結果として導入する。 この仮説は、「神はサイコロを振る」という命題の厳密な実装として解釈できるが、このサイコロゲームがブラウン運動理論のような隠れた変数理論に置き換えられる可能性も残している。
本稿の残りの部分では、このような過程の定義と、 この過程に関連する 拡散方程式の導出について述べる。これは、 ブラウン運動 と 量子力学を 特別な限界とする一般的な設定で行われ、それぞれ 熱方程式 と シュレーディンガー方程式 が得られる。導出は、 ラグランジアン力学 と 確率計算 のツールに大きく依存している。
確率的量子化 この理論の公理は、ネルソンによって定式化された確率的量子化条件に要約することができる。 この条件 に関する非相対論的理論は次のように述べている。 R n {\displaystyle \mathbb {R} ^{n}}
量子粒子の軌道は、複素 セミマルチンゲール の実射影によって記述される。ここ で 、 は連続 有限変化 過程であり、 は複素 マルチンゲール である。 X ( t ) = R e [ Z ( t ) ] {\displaystyle X(t)={\rm {Re}}[Z(t)]} Z ( t ) = C ( t ) + M ( t ) {\displaystyle Z(t)=C(t)+M(t)} C ( t ) {\displaystyle C(t)} M ( t ) {\displaystyle M(t)} 軌道は確率的に行動を極端化する 。 S = E [ ∫ L d t ] {\displaystyle S=\mathbb {E} \left[\int L\,dt\right]} マル チンゲールは、 独立した増分 と有限 モーメント を持つ連続過程である 。さらに、その 二次変化 は構造関係によって固定される。 ここで 、は粒子の質量、 換算 プランク定数 、 は無次元定数、は クロネッカーデルタ である 。 M ( t ) {\displaystyle M(t)} d [ M i , M j ] = α ℏ m δ i j d t {\displaystyle d[M^{i},M^{j}]={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}\,dt} m {\displaystyle m} ℏ {\displaystyle \hbar } α = | α | e i ϕ ∈ C {\displaystyle \alpha =|\alpha |e^{{\rm {i}}\phi }\in \mathbb {C} } δ i j {\displaystyle \delta ^{ij}} 時間を逆転させたプロセスが存在し、同じ力学法則に従います。 分解 と有限変化を持つという事実を用いると、 と の2次変化は 次 のように与えられる ことがわかる。 Z = X + i Y {\displaystyle Z=X+{\rm {i}}\,Y} C {\displaystyle C} X {\displaystyle X} Y {\displaystyle Y}
( d [ X i , X j ] d [ X i , Y j ] d [ Y i , X j ] d [ Y i , Y j ] ) = | α | ℏ 2 m δ i j ( 1 + cos ϕ sin ϕ sin ϕ 1 − cos ϕ ) d t {\displaystyle {\begin{pmatrix}d[X^{i},X^{j}]&d[X^{i},Y^{j}]\\d[Y^{i},X^{j}]&d[Y^{i},Y^{j}]\end{pmatrix}}={\frac {|\alpha |\,\hbar }{2\,m}}\,\delta ^{ij}\,{\begin{pmatrix}1+\cos \phi &\sin \phi \\\sin \phi &1-\cos \phi \end{pmatrix}}\,dt} したがって、 レヴィの ブラウン運動の特徴付けにより、 およびは、 有限変化過程 によって記述されるドリフト 、 に比例する拡散定数 、および角度 に依存する相関を伴う 、相関のある2つのウィーナー過程を記述する。これらの過程は、量子極限において最大相関し、 に関連し 、 に対応する 。一方、ブラウン運動極限においては、 に関連し 、 に対応する無相関となる 。 X ( t ) {\displaystyle X(t)} Y ( t ) {\displaystyle Y(t)} C ( t ) {\displaystyle C(t)} | α | ∈ [ 0 , ∞ ) {\displaystyle |\alpha |\in [0,\infty )} ϕ {\displaystyle \phi } ϕ ∈ { − π 2 , π 2 } {\displaystyle \phi \in \left\{-{\frac {\pi }{2}},{\frac {\pi }{2}}\right\}} α ∈ i × R {\displaystyle \alpha \in {\rm {i}}\times \mathbb {R} } ϕ ∈ { 0 , π } {\displaystyle \phi \in \{0,\pi \}} α ∈ R {\displaystyle \alpha \in \mathbb {R} }
この量子化手順を説明するために確率的量子化という用語が導入されたのは1970年代です。 現在では、 確率的量子化は、1981年に パリシ とウーによって開発された枠組みを指すことが一般的です。 そのため、確率力学で開発された量子化手順は、ネルソンの確率的量子化または確率化と呼ばれることもあります。
プロセスの速度 確率過程は ほぼ確実に 微分不可能であり 、そのため過程に沿った速度は 明確に定義されない。しかし、条件付き期待値を用いて定義される速度場が存在する。これらは次のように与えられる。 Z ( t ) {\displaystyle Z(t)} Z ˙ ( t ) = d Z ( t ) d t {\displaystyle {\dot {Z}}(t)={\frac {dZ(t)}{dt}}}
w + ( x , t ) = lim d t → 0 E [ Z ( t + d t ) − Z ( t ) d t | X ( t ) = x ] , {\displaystyle w_{+}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {Z(t+dt)-Z(t)}{dt}}\,{\Big |}\,X(t)=x\right],} w − ( x , t ) = lim d t → 0 E [ Z ( t ) − Z ( t − d t ) d t | X ( t ) = x ] , {\displaystyle w_{-}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {Z(t)-Z(t-dt)}{dt}}\,{\Big |}\,X(t)=x\right],} は、過程 に沿った イトー積分 に関連付けることができます 。この過程は微分不可能であるため、これらの速度は一般に互いに等しくありません。この事実の物理的な解釈は次のとおりです。粒子は、任意の時点でその速度を から に瞬間的に変化させるランダムな力を受けます 。2つの速度場は等しくないため、過程 に対する速度の唯一の概念は存在しません 。実際、 Z ( t ) {\displaystyle Z(t)} t {\displaystyle t} w − {\displaystyle w_{-}} w + {\displaystyle w_{+}} Z ( t ) {\displaystyle Z(t)} w a = a w + + ( 1 − a ) w − {\displaystyle w_{a}=a\ w_{+}+(1-a)\ w_{-}} は、 過程の速度 の有効な選択肢を表します。これは特に、 で示される 特殊なケースに当てはまり、 に沿った ストラトノビッチ積分 に関連付けることができます 。 a ∈ [ 0 , 1 ] {\displaystyle a\in [0,1]} Z ( t ) {\displaystyle Z(t)} a = 1 2 {\displaystyle a={\frac {1}{2}}} w ∘ = w + + w − 2 {\displaystyle w_{\circ }={\frac {w_{+}+w_{-}}{2}}} Z ( t ) {\displaystyle Z(t)}
は2次非 零変分 を持つので 、次式で表される2次速度場をさらに定義することができる。 Z ( t ) {\displaystyle Z(t)}
w 2 , + ( x , t ) = lim d t → 0 E [ [ Z ( t + d t ) − Z ( t ) ] ⊗ [ Z ( t + d t ) − Z ( t ) ] d t | X ( t ) = x ] , {\displaystyle w_{2,+}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {[Z(t+dt)-Z(t)]\otimes [Z(t+dt)-Z(t)]}{dt}}\ {\Big |}\ X(t)=x\right],} w 2 , − ( x , t ) = lim d t → 0 E [ [ Z ( t ) − Z ( t − d t ) ] ⊗ [ Z ( t ) − Z ( t − d t ) ] d t | X ( t ) = x ] . {\displaystyle w_{2,-}(x,t)=\lim _{dt\rightarrow 0}\mathbb {E} \left[{\frac {[Z(t)-Z(t-dt)]\otimes [Z(t)-Z(t-dt)]}{dt}}\ {\Big |}\ X(t)=x\right].} 時間可逆性の公理は、これら2つの場に という関係を課す 。さらに、二次変化を固定する構造関係を用いると、 となる 。したがって、ストラトノビッチの定式化では、速度の2次成分はゼロ、すなわち となる 。 w 2 , ± = ± w 2 {\displaystyle w_{2,\pm }=\pm \,w_{2}} w 2 i j ( x , t ) = α ℏ m δ i j {\displaystyle w_{2}^{ij}(x,t)={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}} w 2 , ∘ = 0 {\displaystyle w_{2,\circ }=0}
速度の実数部と虚数部は次のように表される。
v = R e ( w ) a n d u = I m ( w ) . {\displaystyle v={\rm {Re}}(w)\qquad {\rm {and}}\qquad u={\rm {Im}}(w)\,.} これらの速度場の存在を利用して、速度過程を 伊藤積分 によって正式に定義することができます 。同様に、過程を ストラトノビッチ積分 によって正式に定義し 、2 次速度過程を スティルチェス積分 によって正式に定義することができます 。構造関係 を用いて、2 次速度過程が によって与えられることがわかります 。しかし、過程 と は 明確に定義されていません。1 次モーメントは存在し によって与えられます が、2 次モーメントは発散します。つまり です 。この発散の物理的な解釈は、位置表現では位置は正確に分かっていますが、速度には無限の不確実性があるということです。 W ± ( t ) {\displaystyle W_{\pm }(t)} ∫ W ± ( t ) d t = ∫ d ± Z ( t ) {\displaystyle \int W_{\pm }(t)\ dt=\int d_{\pm }Z(t)} W ∘ ( t ) {\displaystyle W_{\circ }(t)} ∫ W ∘ ( t ) d t = ∫ d ∘ Z ( t ) {\displaystyle \int W_{\circ }(t)\ dt=\int d_{\circ }Z(t)} W 2 ( t ) {\displaystyle W_{2}(t)} ∫ W 2 ( t ) d t = ∫ d [ Z , Z ] ( t ) {\displaystyle \int W_{2}(t)\ dt=\int d[Z,Z](t)} W 2 i j ( t ) = α ℏ m δ i j {\displaystyle W_{2}^{ij}(t)={\frac {\alpha \ \hbar }{m}}\ \delta ^{ij}} W ± ( t ) {\displaystyle W_{\pm }(t)} W ∘ ( t ) {\displaystyle W_{\circ }(t)} E [ W ± ( t ) | X ( t ) ] = w ± ( X ( t ) , t ) {\displaystyle \mathbb {E} [W_{\pm }(t)\ |\ X(t)]=w_{\pm }(X(t),t)} E [ | | W ± ( t ) | | 2 | X ( t ) ] = ∞ {\displaystyle \mathbb {E} [||W_{\pm }(t)||^{2}\ |\ X(t)]=\infty }
確率的作用 確率的量子化条件は、確率的軌道が確率的作用を極値化しなければならないことを規定している が、確率的ラグランジアンを規定していない。このラグランジアンは、標準的な手順を用いて古典的 ラグランジアン から得ることができる 。ここでは、以下の形式の古典的ラグランジアンを考える。 S = E [ ∫ L d t ] {\displaystyle S=\mathbb {E} \left[\int L\,dt\right]} L {\displaystyle L}
L ( x , v , t ) = m 2 δ i j v i v j + q A i ( x , t ) v i − U ( x , t ) . {\displaystyle L(x,v,t)={\frac {m}{2}}\delta _{ij}v^{i}v^{j}+q\,A_{i}(x,t)\,v^{i}-{\mathfrak {U}}(x,t).} ここで、は 位相空間 ( 接線束 ) の座標、 は 上の 計量 を記述する クロネッカーのデルタ 、 は粒子の質量、 ベクトルポテンシャル における 電荷 、 はスカラーポテンシャルである。さらに、アインシュタイン の 総和法 が仮定されている。 ( x , v ) {\displaystyle (x,v)} δ i j {\displaystyle \delta _{ij}} R n {\displaystyle \mathbb {R} ^{n}} m {\displaystyle m} q {\displaystyle q} A {\displaystyle A} U {\displaystyle {\mathfrak {U}}}
このラグランジアンの重要な性質はゲージ不変性 の原理である。これは 、元の作用に全微分項を追加することで 新たな作用を定義することで明示的に表すことができる。 S ~ {\displaystyle {\tilde {S}}}
S ~ [ x ( t ) ] = S [ x ( t ) ] + ∫ d F ( x , t ) = ∫ [ m 2 δ i j v i v j + q A i v i − U + ∂ t F + v i ∂ i F ] d t = ∫ [ m 2 δ i j v i v j + q A ~ i v i − U ~ ] d t , {\displaystyle {\tilde {S}}[x(t)]=S[x(t)]+\int dF(x,t)=\int \left[{\frac {m}{2}}\delta _{ij}v^{i}v^{j}+qA_{i}v^{i}-{\mathfrak {U}}+\partial _{t}F+v^{i}\partial _{i}F\right]dt=\int \left[{\frac {m}{2}}\delta _{ij}v^{i}v^{j}+q{\tilde {A}}_{i}v^{i}-{\tilde {\mathfrak {U}}}\right]dt,} ここで、 およびである 。したがって、作用に全微分を加えてもダイナミクスは影響を受けないため、任意の微分可能関数に対するポテンシャルの上記の再定義の下では、作用はゲージ不変である 。 A ~ i = A i + q − 1 ∂ i F {\displaystyle {\tilde {A}}_{i}=A_{i}+q^{-1}\partial _{i}F} U ~ = U − ∂ t F {\displaystyle {\tilde {\mathfrak {U}}}={\mathfrak {U}}-\partial _{t}F} F {\displaystyle F}
この古典的ラグランジアンに対応する確率的ラグランジアンを構成するためには、このゲージ不変性を尊重する上記のラグランジアンの最小拡張を探す必要がある。 この理論のストラトノビッチ定式化では、これは簡単に行うことができる。なぜなら、ストラトノビッチ定式化における微分演算子は次のように与えられるからである。
∫ d ∘ F ( x , t ) = ∫ ( ∂ t F + v ∘ i ∂ i F ) d t . {\displaystyle \int d_{\circ }F(x,t)=\int \left(\partial _{t}F+v_{\circ }^{i}\partial _{i}F\right)dt\,.} したがって、ストラトノビッチラグランジアンは、古典的な速度を複素速度に 置き換えることによって得られる 。 v {\displaystyle v} w ∘ {\displaystyle w_{\circ }}
L ∘ ( x , w ∘ , t ) = m 2 δ i j w ∘ i w ∘ j + q A i w ∘ i − U . {\displaystyle L_{\circ }(x,w_{\circ },t)={\frac {m}{2}}\delta _{ij}w_{\circ }^{i}w_{\circ }^{j}+qA_{i}w_{\circ }^{i}-{\mathfrak {U}}\,.} 伊藤の定式化では、全微分が 伊藤の補題 によって与えられるため、事態はより複雑になる。 ∫ d ± F ( x , t ) = ∫ ( ∂ t F + v ± i ∂ i F ± 1 2 v 2 i j ∂ j ∂ i F ) d t . {\displaystyle \int d_{\pm }F(x,t)=\int \left(\partial _{t}F+v_{\pm }^{i}\partial _{i}F\pm {\frac {1}{2}}v_{2}^{ij}\partial _{j}\partial _{i}F\right)dt\,.}
2階微分項の存在によりゲージ不変性は破れる。しかし、ベクトルポテンシャルの微分をラグランジアンに加えることでゲージ不変性は回復できる。したがって、確率ラグランジアンは次のようなラグランジアンで与えられる。
L ± ( x , w ± , w 2 , t ) = m 2 δ i j w ± i w ± j + q A i w ± i ± q 2 ∂ j A i w 2 i j − U . {\displaystyle L_{\pm }(x,w_{\pm },w_{2},t)={\frac {m}{2}}\delta _{ij}w_{\pm }^{i}w_{\pm }^{j}+q\,A_{i}w_{\pm }^{i}\pm {\frac {q}{2}}\partial _{j}A_{i}w_{2}^{ij}-{\mathfrak {U}}\,.} 確率的作用はストラトノビッチラグランジアンを用いて定義することができ、これは発散項を除けば伊藤ラグランジアンによって定義された作用に等しい:
S = E [ ∫ L ∘ d t ] = E [ ∫ L ± d t ] ± E [ ∫ L ∞ d t ] . {\displaystyle S=\mathbb {E} \left[\int L_{\circ }\,dt\right]=\mathbb {E} \left[\int L_{\pm }dt\right]\pm \mathbb {E} \left[\int L_{\infty }dt\right]\,.} 発散項は計算可能であり、次のように与えられる。
E [ ∫ L ∞ d t ] = m 2 ∮ γ δ i j w 2 i j t d t = α ℏ π i ∑ i = 1 n k i , {\displaystyle \mathbb {E} \left[\int L_{\infty }dt\right]={\frac {m}{2}}\oint _{\gamma }{\frac {\delta _{ij}w_{2}^{ij}}{t}}dt=\alpha \,\hbar \,\pi \,{\rm {i}}\,\sum _{i=1}^{n}k_{i},} ここで、 は 、 における極の周り の経路の巻き数を数える巻き数である 。 k i ∈ Z {\displaystyle k_{i}\in \mathbb {Z} } γ ( t ) {\displaystyle \gamma (t)} t = 0 {\displaystyle t=0}
発散項は定数であるため、運動方程式には寄与しない。このため、この項は確率力学の初期の研究では無視されてきた 。しかし、この項を無視すると、確率力学は量子力学における離散スペクトルの出現を説明できなくなる。この問題はヴァルストロム批判 として知られており、発散項を適切に考慮することで解決できる。
確率力学のハミルトン 定式化も存在する 。 これは正準運動量の定義から始まる。
p ∘ , i = ∂ L ∘ ∂ w ∘ i = m δ i j w ∘ j + q A i , {\displaystyle p_{\circ ,i}={\frac {\partial L_{\circ }}{\partial w_{\circ }^{i}}}=m\,\delta _{ij}w_{\circ }^{j}+q\,A_{i}\,,} p ± , i = ∂ L ± ∂ w ± i = m δ i j w ± j + q A i . {\displaystyle p_{\pm ,i}={\frac {\partial L_{\pm }}{\partial w_{\pm }^{i}}}=m\,\delta _{ij}w_{\pm }^{j}+q\,A_{i}\,.} ストラトノビッチ定式化におけるハミルトニアンは、1次 ルジャンドル変換 によって得られる。
H ∘ ( x , p ∘ , t ) = p ∘ , i v ∘ i − L ( x , v ∘ , t ) . {\displaystyle H_{\circ }(x,p_{\circ },t)=p_{\circ ,i}v_{\circ }^{i}-L(x,v_{\circ },t)\,.} 一方、伊藤の定式化では、ハミルトニアンは2次のルジャンドル変換によって得られる。
H ± ( x , p ± , ∂ p ± , t ) = p i ± w ± i ± 1 2 w 2 i j ∂ p i j ± − L ( x , w ± , w 2 , t ) . {\displaystyle H_{\pm }(x,p^{\pm },\partial p^{\pm },t)=p_{i}^{\pm }w_{\pm }^{i}\pm {\frac {1}{2}}w_{2}^{ij}\partial p_{ij}^{\pm }-L(x,w_{\pm },w_{2},t)\,.}
オイラー・ラグランジュ方程式 確率的作用は極限まで極限化することができ、オイラー・ラグランジュ方程式 の確率的バージョンが得られる 。ストラトノビッチの定式化では、これらは次のように与えられる。
∫ d ∘ ( ∂ L ∘ ∂ w ∘ i ) = ∫ ( ∂ L ∘ ∂ x i ) d t . {\displaystyle \int d_{\circ }\left({\frac {\partial L_{\circ }}{\partial w_{\circ }^{i}}}\right)=\int \left({\frac {\partial L_{\circ }}{\partial x^{i}}}\right)dt\,.} 前のセクションで説明したラグランジアンの場合、これは ストラトノビッチ の意味での次の2次 確率微分方程式 につながります。
m δ i j d ∘ 2 Z j ( t ) = q F i j ( X ( t ) , t ) d ∘ Z j ( t ) d t − q ∂ t A i ( X ( t ) , t ) d t 2 − ∂ i U ( X ( t ) , t ) d t 2 , {\displaystyle m\,\delta _{ij}\,d_{\circ }^{2}Z^{j}(t)=q\,F_{ij}(X(t),t)\,d_{\circ }Z^{j}(t)\,dt-q\,\partial _{t}A_{i}(X(t),t)\,dt^{2}-\partial _{i}{\mathfrak {U}}(X(t),t)\,dt^{2}\,,} ここで、場の強度は で与えられる。この式は ニュートンの第二法則 の確率論的バージョンとして機能する 。 F i j := ∂ i A j − ∂ j A i {\displaystyle F_{ij}:=\partial _{i}A_{j}-\partial _{j}A_{i}}
伊藤の定式化では、確率 オイラー・ラグランジュ方程式は 次のように与えられる。
∫ d ± ( ∂ L ± ∂ w ± i ) = ∫ ( ∂ L ± ∂ x i ) d t . {\displaystyle \int d_{\pm }\left({\frac {\partial L_{\pm }}{\partial w_{\pm }^{i}}}\right)=\int \left({\frac {\partial L_{\pm }}{\partial x^{i}}}\right)dt\,.} これは伊藤 の意味で 2階の 確率微分方程式につながり、 ニュートンの第2法則 の確率バージョンとして次のように与え られる
。
m δ i j d ± 2 Z j ( t ) = q F i j ( X ( t ) , t ) d ± Z j ( t ) d t ± α ℏ q 2 m δ j k ∂ k F i j ( X ( t ) , t ) d t 2 − q ∂ t A i ( X ( t ) , t ) d t 2 − ∂ i U ( X ( t ) , t ) d t 2 . {\displaystyle m\,\delta _{ij}\,d_{\pm }^{2}Z^{j}(t)=q\,F_{ij}(X(t),t)\,d_{\pm }Z^{j}(t)\,dt\pm {\frac {\alpha \,\hbar \,q}{2\,m}}\,\delta ^{jk}\,\partial _{k}F_{ij}(X(t),t)\,dt^{2}-q\,\partial _{t}A_{i}(X(t),t)\,dt^{2}-\partial _{i}{\mathfrak {U}}(X(t),t)\,dt^{2}\,.}
ハミルトン・ヤコビ方程式 運動方程式は、古典力学の ハミルトン・ヤコビ 定式化の確率論的一般化によっても得られる。この場合、まずハミルトン主関数を定義する。ラグランジアンの場合 、この関数は次のように定義される
。 L + {\displaystyle L_{+}}
S + ( x , t ; x f , t f ) := − E [ ∫ t t f L + ( X ( s ) , W + ( s ) , W 2 ( s ) , s ) d s | X ( t ) = x , X ( t f ) = x f ] , {\displaystyle S_{+}(x,t;x_{f},t_{f}):=-\mathbb {E} \left[\int _{t}^{t_{f}}L_{+}(X(s),W_{+}(s),W_{2}(s),s)\,ds\,{\Big |}\,X(t)=x,X(t_{f})=x_{f}\right],} ここで、過程は確率オイラー・ラグランジュ方程式に従うと仮定する 。同様に、ラグランジアンに対して 、ハミルトンの主関数は次のように定義される。 { X ( s ) : s ∈ [ t , t f ] } {\displaystyle \{X(s):s\in [t,t_{f}]\}} L − {\displaystyle L_{-}}
S − ( x , t ; x 0 , t 0 ) := E [ ∫ t 0 t L − ( X ( s ) , W − ( s ) , W 2 ( s ) , s ) d s | X ( t ) = x , X ( t 0 ) = x 0 ] , {\displaystyle S_{-}(x,t;x_{0},t_{0}):=\mathbb {E} \left[\int _{t_{0}}^{t}L_{-}(X(s),W_{-}(s),W_{2}(s),s)\,ds\,{\Big |}\,X(t)=x,X(t_{0})=x_{0}\right],} ここで、過程は確率オイラー・ラグランジュ方程式に従うと仮定する 。作用の発散部分により、これらの主関数は同値関係に従う。 { X ( s ) : s ∈ [ t 0 , t ] } {\displaystyle \{X(s):s\in [t_{0},t]\}}
S ~ ± ≡ S ± i f ∃ k ∈ Z n , s u c h t h a t S ~ ± = S ± ± α π i ℏ ∑ i = 1 n k i . {\displaystyle {\tilde {S}}_{\pm }\equiv S_{\pm }\;{\rm {if}}\;\exists k\in \mathbb {Z} ^{n},\;{\rm {such\,that}}\;{\tilde {S}}_{\pm }=S_{\pm }\pm \alpha \,\pi \,{\rm {i}}\,\hbar \,\sum _{i=1}^{n}k_{i}\,.} 主関数を点に関して変化させることで、 ハミルトン・ヤコビ方程式が得られる。これは次のように与えられる。 ( x , t ) {\displaystyle (x,t)}
{ ∂ ∂ x i S ± ( x , t ) = p i ± ( x , t ) , ∂ ∂ t S ± ( x , t ) = − H ± ( x , p ± ( x , t ) , ∂ p ± ( x , t ) , t ) . {\displaystyle {\begin{cases}{\frac {\partial }{\partial x^{i}}}S_{\pm }(x,t)&=p_{i}^{\pm }(x,t)\,,\\{\frac {\partial }{\partial t}}S_{\pm }(x,t)&=-H_{\pm }(x,p_{\pm }(x,t),\partial p_{\pm }(x,t),t)\,.\end{cases}}} これらは古典的な場合と同じように見えることに注意してください。ただし、2番目のハミルトン・ヤコビ方程式のハミルトニアンは、2次のルジャンドル変換を用いて得られます。さらに、作用の発散部分により、非自明な積分制約の形をとる3番目のハミルトン・ヤコビ方程式が存在します。
∮ ( p i ± v ± i ± 1 2 v 2 i j ∂ i p j ) d t = ± α ℏ π i ∑ i = 1 n k i . {\displaystyle \oint \left(p_{i}^{\pm }\,v_{\pm }^{i}\pm {\frac {1}{2}}v_{2}^{ij}\,\partial _{i}p_{j}\right)dt=\pm \alpha \,\hbar \,\pi \,{\rm {i}}\,\sum _{i=1}^{n}k_{i}\,.} 与えられたラグランジアンに対して、最初の2つのハミルトン・ヤコビ方程式は
{ ∂ i S = m δ i j w ± j + q A i , ∂ t S = − m 2 δ i j w ± i w ± j ∓ m 2 δ i j w 2 i k ∂ k w ± j − U . {\displaystyle {\begin{cases}\partial _{i}S&=m\,\delta _{ij}w_{\pm }^{j}+q\,A_{i}\,,\\\partial _{t}S&=-{\frac {m}{2}}\,\delta _{ij}w_{\pm }^{i}w_{\pm }^{j}\mp {\frac {m}{2}}\,\delta _{ij}w_{2}^{ik}\partial _{k}w_{\pm }^{j}-{\mathfrak {U}}\,.\end{cases}}} これら2つの式を組み合わせると、
[ m δ i j ( ∂ t + w ± k ∂ k ± 1 2 w 2 k l ∂ l ∂ k ) − q F i j ] w ± j = ± q 2 w 2 j k ∂ k F i j − q ∂ t A i − ∂ i U . {\displaystyle \left[m\,\delta _{ij}\left(\partial _{t}+w_{\pm }^{k}\partial _{k}\pm {\frac {1}{2}}\,w_{2}^{kl}\partial _{l}\partial _{k}\right)-q\,F_{ij}\right]w_{\pm }^{j}=\pm {\frac {q}{2}}\,w_{2}^{jk}\partial _{k}F_{ij}-q\,\partial _{t}A_{i}-\partial _{i}{\mathfrak {U}}\,.} を用いて 、積分条件と初期条件 または終了条件を課したこの方程式 を について解くことができる 。そして、その解を伊藤方程式に代入することができる。 w 2 i j = α ℏ m δ i j {\displaystyle w_{2}^{ij}={\frac {\alpha \,\hbar }{m}}\,\delta ^{ij}} w + ( x , t 0 ) = w 0 ( x ) {\displaystyle w_{+}(x,t_{0})=w_{0}(x)} w − ( x , t f ) = w f ( x ) {\displaystyle w_{-}(x,t_{f})=w_{f}(x)} w ± ( x , t ) {\displaystyle w_{\pm }(x,t)}
{ d ± Z i ( t ) = w ± i ( x , t ) d t + d M i ( t ) , d [ M i , M j ] ( t ) = α ℏ m δ i j d t , {\displaystyle {\begin{cases}d_{\pm }Z^{i}(t)&=w_{\pm }^{i}(x,t)\,dt+dM^{i}(t)\,,\\d[M^{i},M^{j}](t)&={\frac {\alpha \,\hbar }{m}}\delta ^{ij}\,dt\,,\end{cases}}} これは過程 について解くことができる。したがって、 ( でラベル付けされた未来指向方程式の 場合)初期条件 または ( でラベル付けされた過去指向方程式の場合)終了条件を指定すると、 粒子の軌跡を記述する 唯一の確率過程が見つかる。 { Z ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{Z(t):t\in [t_{0},t_{f}]\}} X ( t 0 ) = x 0 {\displaystyle X(t_{0})=x_{0}} + {\displaystyle +} X ( t f ) = x f {\displaystyle X(t_{f})=x_{f}} − {\displaystyle -} { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}
拡散方程式 ネルソン確率力学の重要な結果は、仮定された確率過程からシュレーディンガー方程式を導出することである。この導出において、ハミルトン・ヤコビ方程式は
{ ∂ ∂ x i S ± ( x , t ) = p i ± ( x , t ) ∂ ∂ t S ± ( x , t ) = − H ± ( x , p ± ( x , t ) , ∂ p ± ( x , t ) , t ) {\displaystyle {\begin{cases}{\frac {\partial }{\partial x^{i}}}S_{\pm }(x,t)&=p_{i}^{\pm }(x,t)\\{\frac {\partial }{\partial t}}S_{\pm }(x,t)&=-H_{\pm }(x,p_{\pm }(x,t),\partial p_{\pm }(x,t),t)\end{cases}}} を組み合わせると、次の式が得られる。
2 m ( ∂ t S ± + U ) + δ i j ( ∂ i S ± ∂ j S ± ± α ℏ ∂ j ∂ i S ± − 2 q A i ∂ j S ± ∓ α ℏ q ∂ j A i + q 2 A i A j ) = 0 . {\displaystyle 2\,m\left(\partial _{t}S_{\pm }+{\mathfrak {U}}\right)+\delta ^{ij}\left(\partial _{i}S_{\pm }\partial _{j}S_{\pm }\pm \alpha \,\hbar \,\partial _{j}\partial _{i}S_{\pm }-2\,q\,A_{i}\partial _{j}S_{\pm }\mp \alpha \,\hbar \,q\,\partial _{j}A_{i}+q^{2}\,A_{i}A_{j}\right)=0\,.} 次に波動関数を定義する。
Ψ ± ( x , t ) = exp ( ± S ± ( x , t ) α ℏ ) . {\displaystyle \Psi _{\pm }(x,t)=\exp \left(\pm {\frac {S_{\pm }(x,t)}{\alpha \,\hbar }}\right).} ハミルトンの主関数は多価関数であるため、波動関数は同値関係に従うことがわかる。
Ψ ~ + ≡ Ψ + i f Ψ ~ + = ± Ψ + a n d Ψ ~ − ≡ Ψ − i f Ψ ~ − = ± Ψ − . {\displaystyle {\tilde {\Psi }}_{+}\equiv \Psi _{+}\quad {\rm {if}}\quad {\tilde {\Psi }}_{+}=\pm \Psi _{+}\qquad {\rm {and}}\qquad {\tilde {\Psi }}_{-}\equiv \Psi _{-}\quad {\rm {if}}\quad {\tilde {\Psi }}_{-}=\pm \Psi _{-}\,.} さらに、波動関数は複素 拡散方程式に従う。
− α ℏ ∂ ∂ t Ψ + = [ δ i j 2 m ( α ℏ ∂ ∂ x i + q A i ) ( α ℏ ∂ ∂ x j + q A j ) + U ] Ψ + , {\displaystyle -\alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{+}=\left[{\frac {\delta ^{ij}}{2\,m}}\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}+q\,A_{i}\right)\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{j}}}+q\,A_{j}\right)+{\mathfrak {U}}\right]\Psi _{+}\,,} α ℏ ∂ ∂ t Ψ − = [ δ i j 2 m ( α ℏ ∂ ∂ x i + q A i ) ( α ℏ ∂ ∂ x j + q A j ) + U ] Ψ − . {\displaystyle \alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{-}=\left[{\frac {\delta ^{ij}}{2\,m}}\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}+q\,A_{i}\right)\left(\alpha \,\hbar \,{\frac {\partial }{\partial x^{j}}}+q\,A_{j}\right)+{\mathfrak {U}}\right]\Psi _{-}\,.} したがって、確率力学の公理を解く任意の過程に対して、これらの拡散方程式に従う波動関数を構築することができます。ハミルトンの主関数における同値関係により、逆のことも成り立ちます。つまり、これらの複素拡散方程式の任意の解に対して、確率力学の公理を解く確率過程を構築することができます。同様の結果は 、ファインマン-カックの定理 によって確立されています 。 { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}
最後に、確率密度を構築することができる。
ρ ± ( x , t ) := | Ψ ± ( x , t ) | 2 ∫ | Ψ ± ( y , t ) | 2 d n y , {\displaystyle \rho _{\pm }(x,t):={\frac {|\Psi _{\pm }(x,t)|^{2}}{\int |\Psi _{\pm }(y,t)|^{2}d^{n}y}}\,,} これは、プロセス の遷移確率を表します 。より正確には、 は、 システムが最終的に状態 になった場合に、状態 になる確率を表します 。したがって、 の拡散方程式は、 プロセス の コルモゴロフ後退方程式 として解釈できます 。同様に、 は 、時間を逆に展開した場合に、 システムが最終的に状態 になった場合に、 状態 になる確率を表します 。したがって、 の拡散方程式は、プロセス を過去に向けて展開した場合 の コルモゴロフ後退方程式 として解釈できます。時間方向を反転すると、 は、 システムが状態で開始し、時間を順方向に展開した場合に、 状態 になる確率を表すこと がわかります。したがって、 の拡散方程式は、 プロセス を未来に向けて展開した場合
の コルモゴロフ前進方程式 としても解釈できます。 { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}} ρ + {\displaystyle \rho _{+}} ( x , t ) {\displaystyle (x,t)} ( x f , t f ) {\displaystyle (x_{f},t_{f})} Ψ + {\displaystyle \Psi _{+}} { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}} ρ − {\displaystyle \rho _{-}} ( x , t ) {\displaystyle (x,t)} ( x 0 , t 0 ) {\displaystyle (x_{0},t_{0})} Ψ − {\displaystyle \Psi _{-}} { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}} ρ − {\displaystyle \rho _{-}} ( x , t ) {\displaystyle (x,t)} ( x 0 , t 0 ) {\displaystyle (x_{0},t_{0})} Ψ − {\displaystyle \Psi _{-}} { X ( t ) : t ∈ [ t 0 , t f ] } {\displaystyle \{X(t):t\in [t_{0},t_{f}]\}}
数学的側面
限定的なケース この理論にはさまざまな特別な制限が含まれています。
の古典的極限 。この場合、プロセス と補助プロセスは 2つの分離された決定論的軌道を記述します。 α = 0 {\displaystyle \alpha =0} X {\displaystyle X} Y {\displaystyle Y} のブラウン運動極限 。この場合、過程は ウィーナー過程 ( ブラウン運動 とも呼ばれる)を表し 、上記の結果は ファインマン-カック定理 によって確立されますが、補助過程は 決定論的過程を表します。 α ∈ ( 0 , ∞ ) {\displaystyle \alpha \in (0,\infty )} X {\displaystyle X} Y {\displaystyle Y} の量子極限 。この場合、過程 と補助過程は 2つの正の相関を持つウィーナー過程を記述します。 α ∈ i × ( 0 , ∞ ) {\displaystyle \alpha \in {\rm {i}}\times (0,\infty )} X {\displaystyle X} Y {\displaystyle Y} の時間反転ブラウン運動極限 。この場合、過程 は 決定論的過程を記述し、補助過程 は ウィーナー過程 を記述します 。 α ∈ ( − ∞ , 0 ) {\displaystyle \alpha \in (-\infty ,0)} X {\displaystyle X} Y {\displaystyle Y} を伴う時間反転量子極限 。この場合、過程 と補助過程は 2つの負の相関を持つウィーナー過程を記述する。 α ∈ i × ( − ∞ , 0 ) {\displaystyle \alpha \in {\rm {i}}\times (-\infty ,0)} X {\displaystyle X} Y {\displaystyle Y} 初期条件または終了条件 ( を意味する) を伴うブラウン極限では 、過程 と は 分離されているため、補助過程のダイナミクスは 無視でき、 は 実 ウィーナー過程 によって記述されます。 を伴うその他のすべてのケースでは 、過程は互いに結合しているため、 のダイナミクスを導出する際には補助過程を考慮する必要があります 。 u ( x , t 0 ) = 0 {\displaystyle u(x,t_{0})=0} u ( x , t f ) = 0 {\displaystyle u(x,t_{f})=0} u ( x , t ) = 0 ∀ t {\displaystyle u(x,t)=0\quad \forall t} X {\displaystyle X} Y {\displaystyle Y} Y {\displaystyle Y} X {\displaystyle X} α ≠ 0 {\displaystyle \alpha \neq 0} Y {\displaystyle Y} X {\displaystyle X}
時間反転対称性 この理論は時間反転操作の下では対称的である 。 ( t , α , q ) ↔ ( − t , − α , − q ) {\displaystyle (t,\alpha ,q)\leftrightarrow (-t,-\alpha ,-q)}
ブラウン運動の極限では理論は最大 散逸的で あるが、量子の極限では ユニタリ的 であり、
d d t ∫ | Ψ ± ( y , t ) | 2 d n y | α ∈ i × R = 0 . {\displaystyle {\frac {d}{dt}}\int |\Psi _{\pm }(y,t)|^{2}d^{n}y{\Big |}_{\alpha \in {\rm {i}}\times \mathbb {R} }=0\,.}
標準的な交換関係 拡散方程式は次のように書き直すことができる。
∓ α ℏ ∂ ∂ t Ψ ± = H ^ ( x ^ , p ^ ± , t ) Ψ ± , {\displaystyle \mp \alpha \,\hbar \,{\frac {\partial }{\partial t}}\Psi _{\pm }={\hat {H}}({\hat {x}},{\hat {p}}^{\pm },t)\,\Psi _{\pm }\,,} ここで ハミルトン 演算子である 。これにより、位置演算子と運動量演算子を次のように導入することができる。 H ^ {\displaystyle {\hat {H}}}
x ^ i = x i a n d p ^ i ± = ± α ℏ ∂ ∂ x i , {\displaystyle {\hat {x}}^{i}=x^{i}\qquad {\rm {and}}\qquad {\hat {p}}_{i}^{\pm }=\pm \alpha \,\hbar \,{\frac {\partial }{\partial x^{i}}}\,,} ハミルトニアンはおなじみの形をとる
H ( x ^ , p ^ , t ) = δ i j 2 m [ p ^ i − q A i ( x ^ , t ) ] [ p ^ j − q A j ( x ^ , t ) ] + U ( x ^ , t ) . {\displaystyle H({\hat {x}},{\hat {p}},t)={\frac {\delta ^{ij}}{2\,m}}\,{\Big [}{\hat {p}}_{i}-q\,A_{i}({\hat {x}},t){\Big ]}{\Big [}{\hat {p}}_{j}-q\,A_{j}({\hat {x}},t){\Big ]}+{\mathfrak {U}}({\hat {x}},t)\,.} これらの演算子は 標準的な交換関係に従う
[ x ^ i , p ^ j ± ] = ∓ α ℏ δ j i . {\displaystyle [{\hat {x}}^{i},{\hat {p}}_{j}^{\pm }]=\mp \alpha \,\hbar \,\delta _{j}^{i}\,.}
参照
注記
参考文献
論文 ド・ブロイ、L. (1967)。 「Le Mouvement Brownien d'une Particule Dans Son Onde」。 CRアカデミー。科学 。 B264 :1041。 キャメロン, RH (1960). 「ウィーナー積分とファインマン積分を繋ぐ積分族」. 数学物理学ジャーナル . 39 ( 1–4 ): 126– 140. doi :10.1002/sapm1960391126. ISSN 0097-1421. Daletskii, Yu L (1962-10-31). 「作用素進化方程式と関連した関数積分」. ロシア数学概論 . 17 (5): 1– 107. Bibcode :1962RuMaS..17....1D. doi :10.1070/RM1962v017n05ABEH004121. ISSN 0036-0279. デイビッドソン, MP (1979). 「量子力学の確率論的定式化における量子作用素代数の起源」. 数理物理学論文集 . 3 (5): 367– 376. arXiv : quant-ph/0112099 . Bibcode :1979LMaPh...3..367D. doi :10.1007/BF00397209. ISSN 0377-9017. S2CID 6416365. Doob, JL (1942). 「ブラウン運動と確率方程式」. Annals of Mathematics . 43 (2): 351– 369. doi :10.2307/1968873. JSTOR 1968873. アインシュタイン、A. (1905)。 「分子動力学理論を理解するために、完全な理論を確立する必要があります。」 アンナレン・デア・フィジーク 。 322 (8): 549–560 。 書誌コード :1905AnP...322..549E。 土井 : 10.1002/andp.19053220806 。 ISSN 0003-3804。 フェニエス、I. (1946)。 「シュレディンガー方程式の演繹」。 アクタ・ボリャイアナ 。 1 (5):ch. 2. フェニエス、I. (1952)。 「Eine wahrscheinlichkeitstheoretische Begründung und Interpretation der Quantenmechanik」。 物理学の時代 。 132 (1): 81–106 。 ビブコード :1952ZPhy..132...81F。 土井 :10.1007/BF01338578。 ISSN 1434-6001。 S2CID 119581427。 ゲラ, F. (1981). 「確率力学と確率場理論の構造的側面」 . 物理学レポート . 77 (3): 263– 312. 書誌コード :1981PhR....77..263G. doi :10.1016/0370-1573(81)90078-8. Kac, M. (1949). 「ある種のウィーナー関数の超関数について」. アメリカ数学会誌 . 65 (1): 1– 13. doi : 10.2307/1990512 . JSTOR 1990512. Lindgren, J.; Liukkonen, J. (2019). 「量子力学は時空上の確率的最適化を通して理解できる」. Scientific Reports . 9 (1): 19984. Bibcode :2019NatSR...919984L. doi : 10.1038/s41598-019-56357-3 . PMC 6934697. PMID 31882809 . ネルソン, E. (1966). 「ニュートン力学からのシュレーディンガー方程式の導出」 . フィジカル・レビュー . 150 (4): 1079– 1085. 書誌コード :1966PhRv..150.1079N. doi :10.1103/PhysRev.150.1079. ネルソン, E. (1986). 「場の理論と確率力学の未来」. アルベベリオ, S.; カサティ, G.; メルリーニ, D. (編). 古典系および量子系における確率過程 . 物理学講義ノート. 第262巻. ベルリン: シュプリンガー出版社. pp. 438– 469. doi :10.1007/3-540-17166-5. ISBN 978-3-662-13589-1 . OCLC 864657129。 ネルソン, E. (2014). 「相対論的場の確率力学」. Journal of Physics: Conference Series . 504 (1) 012013. Bibcode :2014JPhCS.504a2013N. doi : 10.1088/1742-6596/504/1/012013 . S2CID 123706792. Pavon, M. (1995). 「確率力学の新しい定式化」 . Physics Letters A. 209 ( 3–4 ) : 143–149 . Bibcode :1995PhLA..209..143P. doi :10.1016/0375-9601(95)00847-4. Perrin, J. (1911). 「ブラウン運動と分子的実在性」. Nature . 86 (2160): 105. Bibcode :1911Natur..86..105R. doi : 10.1038/086105a0 . hdl : 2027/mdp.39015010952425 . ISSN 0028-0836. ウォールストロム, TC (1989). 「確率力学からのシュレーディンガー方程式の導出について」 . 物理学基礎レターズ . 2 (2): 113– 126. 書誌コード :1989FoPhL...2..113W. doi :10.1007/BF00696108. Wallstrom, TC (1994). 「シュレーディンガー方程式とマーデルング流体力学方程式の不等式」 . Physical Review A. 49 ( 3): 1613– 1617. Bibcode :1994PhRvA..49.1613W. doi :10.1103/PhysRevA.49.1613. PMID :9910408. ウィーナー、ノーバート (1923). 「微分空間」. 数学物理学ジャーナル . 2 ( 1–4 ): 131–174 . doi :10.1002/sapm192321131. ISSN 0097-1421. 安江 憲一 (1979). 「確率量子化:レビュー」 . 国際理論物理学ジャーナル . 18 (12): 861– 913. Bibcode :1979IJTP...18..861Y. doi :10.1007/BF00669566. S2CID 120858652. ザンブリニ, JC (1985). 「確率力学:確率変分法のレビュー」 . 国際理論物理学ジャーナル . 24 (3): 277– 327. Bibcode :1985IJTP...24..277Z. doi :10.1007/BF00669792. S2CID 122076991. Huang, Q.; Zambrini, JC (2023). 「二次微分幾何学から確率幾何力学へ」. Journal of Nonlinear Science . 33 (67): 1– 127. arXiv : 2201.03706 . Bibcode :2023JNS....33...67H. doi : 10.1007/s00332-023-09917-x .
本