この記事の主題は、ランダムに動く粒子の平均的な振る舞いを、関数から関数への操作として扱う方法です。たとえば粒子が位置 x x x から出発して時間 t t t 後にいる場所での、ある量 f f f の平均を、出発点 x x x の関数として並べると、熱が広がるときと同じ形の操作になります。さらに粒子の通った道に沿って重みを掛けて平均すると、量子力学の Schrödinger 作用素が時間発展させる様子を、道についての平均として書けます。読み終えると、その公式がどこまで特異なポテンシャルで意味をもつか、また基底状態を使うとなぜ再びランダムな動きの話に戻せるかを説明できます。
Brown 運動の Markov 性 (定理 2.2 )は、時刻 s s s 以降の道が「位置 X s X_s X s から新しく出発した Brown 運動」と同じ法則をもつことです。ここから、関数に作用する半群 P t f ( x ) = E x [ f ( X t ) ] P_tf(x)=E_x[f(X_t)] P t f ( x ) = E x [ f ( X t )] が得られます。
Feller 半群 は C 0 ( E ) C_0(E) C 0 ( E ) 上の正値・縮小・強連続な半群です(定義 3.1 )。生成子は正の最大値原理 を満たし(命題 3.2 )、半群は生成子で一意に決まります(定理 3.4 )。Brown 運動の生成子は 1 2 Δ \frac12\Delta 2 1 Δ です(定理 3.5 )。
Feynman–Kac 公式 T t V f ( x ) = E x [ e − ∫ 0 t V ( X s ) d s f ( X t ) ] T^V_tf(x)=E_x\bigl[e^{-\int_0^tV(X_s)ds}f(X_t)\bigr] T t V f ( x ) = E x [ e − ∫ 0 t V ( X s ) d s f ( X t ) ] は、有界一様連続な V V V では 1 2 Δ − V \frac12\Delta-V 2 1 Δ − V の生成する半群です(定理 4.2 )。Kato クラス の V V V (Coulomb ポテンシャルを含む)では、Khasminskii の補題 (補題 4.5 )が指数モーメントを抑え、半群が意味をもちます(定理 4.6 )。
Doob の h 変換 は、Q t h = e λ t h Q_th=e^{\lambda t}h Q t h = e λ t h を満たす正の関数 h h h で核を h ( y ) / h ( x ) h(y)/h(x) h ( y ) / h ( x ) 倍して Markov 半群に直す操作です(定理 5.2 )。Schrödinger 作用素の基底状態 ψ \psi ψ で行うと、ドリフト ∇ log ψ \nabla\log\psi ∇ log ψ をもつ拡散が得られます(基底状態変換 、定理 5.4 )。
Stratonovich 積分 は台形則の Riemann 和の極限で、伊藤積分との差は 1 2 ∫ f ′ ( B s ) d s \frac12\int f'(B_s)ds 2 1 ∫ f ′ ( B s ) d s です(定理 6.2 )。その代わり通常の連鎖律が成り立ちます(系 6.3 )。
この記事で使う事前知識
強連続半群と生成子 : 時間 t t t でパラメータ付けされた作用素の族で、t → 0 t\to0 t → 0 で恒等作用素に近づくもの。生成子はその t = 0 t=0 t = 0 での微分です。→ 1 径数ユニタリ群と強連続半群
期待値 : 確率変数の値を確率で重み付けして平均したもの。→ 期待値
マルチンゲール : 現在までの情報のもとで、将来の値の期待値が現在の値に等しい確率過程。→ マルチンゲール
Wiener 測度 : Brown 運動の道の空間に与える確率測度。→ Kolmogorov の拡張定理と連続性定理
伊藤積分 : Brown 運動で積分する確率積分。各分割区間の左端点の値を使った和の極限として定まります。→ 確率微分方程式
自己共役作用素 : 対称で、定義域も共役作用素のものと一致する作用素。Schrödinger 作用素がその例です。→ 自己共役作用素
量子力学の Schrödinger 作用素 H = − 1 2 Δ + V H=-\frac12\Delta+V H = − 2 1 Δ + V の熱半群 e − t H e^{-tH} e − t H は、Brown 運動の道についての期待値で書けます。これが Feynman–Kac 公式で、構成的場の量子論では作用素の問題(自己共役性、基底状態の一意性、相関関数の評価)を確率論の問題に翻訳する基本の道具です(Feynman–Kac 公式と場の模型 )。その記事では有界連続な V V V について Trotter 積公式 から L 2 L^2 L 2 の等式を示しました。本稿は、同じ公式を確率過程の側から組み立て直します。
組み立ての材料は三つあります。第一は半群と生成子 です。Markov 過程の情報は P t f ( x ) = E x [ f ( X t ) ] P_tf(x)=E_x[f(X_t)] P t f ( x ) = E x [ f ( X t )] という作用素の族に集約され、その微分 A = d d t P t ∣ t = 0 A=\frac{d}{dt}P_t|_{t=0} A = d t d P t ∣ t = 0 が過程の「無限小の動き」を表します。Brown 運動なら A = 1 2 Δ A=\frac12\Delta A = 2 1 Δ です。第二はポテンシャルによる重み です。道に沿って e − ∫ V e^{-\int V} e − ∫ V の重みを掛けると、生成子は A − V A-V A − V に変わります。ただし水素原子の Coulomb ポテンシャル − 1 / ∣ x ∣ -1/|x| − 1/∣ x ∣ のように V V V が特異だと、重みの期待値が有限かどうかから問題になります。これを解決するのが Kato クラスと Khasminskii の補題です。第三はh 変換 です。重み付きの半群は確率を保存しませんが、基底状態で割り直すと再び Markov 過程になり、量子力学の基底状態が「平衡分布 ψ 2 d x \psi^2dx ψ 2 d x をもつ拡散過程」として見えます。
最後の Stratonovich 積分は、ドリフトをもつ拡散を幾何的に(座標によらず)書くときに現れる積分です。伊藤積分と 1 2 \frac12 2 1 の補正だけ異なるという関係を、Riemann 和から直接証明します。
Kolmogorov の拡張定理と連続性定理 の 定理 7.2[Kolmogorov の拡張定理と連続性定理] で作った道の空間を使います。Ω = W d = C ( [ 0 , ∞ ) , R d ) \Omega=W^d=C([0,\infty),\mathbb{R}^d) Ω = W d = C ([ 0 , ∞ ) , R d ) 、X t ( w ) = w ( t ) X_t(w)=w(t) X t ( w ) = w ( t ) を座標過程、F = B ( W d ) \mathcal{F}=\mathcal{B}(W^d) F = B ( W d ) 、F t = σ ( X s : s ≤ t ) \mathcal{F}_t=\sigma(X_s:s\le t) F t = σ ( X s : s ≤ t ) 、P x = W x P_x=\mathbb{W}_x P x = W x (x x x から出発する d d d 次元 Brown 運動の法則)、E x E_x E x をその期待値とします。シフト θ s : W d → W d \theta_s:W^d\to W^d θ s : W d → W d を ( θ s w ) ( t ) = w ( s + t ) (\theta_sw)(t)=w(s+t) ( θ s w ) ( t ) = w ( s + t ) で定めると、X t ∘ θ s = X s + t X_t\circ\theta_s=X_{s+t} X t ∘ θ s = X s + t なので θ s \theta_s θ s は可測です(補題 7.1[Kolmogorov の拡張定理と連続性定理] )。( t , w ) ↦ w ( t ) (t,w)\mapsto w(t) ( t , w ) ↦ w ( t ) は連続なので、ボレル関数 V V V について ( s , w ) ↦ V ( X s ( w ) ) (s,w)\mapsto V(X_s(w)) ( s , w ) ↦ V ( X s ( w )) は B ( [ 0 , ∞ ) ) ⊗ F \mathcal{B}([0,\infty))\otimes\mathcal{F} B ([ 0 , ∞ )) ⊗ F -可測で、Tonelli の定理が使えます。
d d d 次元 Brown 運動の増分 X t 1 − X t 0 , … , X t k − X t k − 1 X_{t_1}-X_{t_0},\ldots,X_{t_k}-X_{t_{k-1}} X t 1 − X t 0 , … , X t k − X t k − 1 (t 0 < ⋯ < t k t_0<\cdots<t_k t 0 < ⋯ < t k )は独立で X t − X s ∼ N ( 0 , ( t − s ) I d ) X_t-X_s\sim N(0,(t-s)I_d) X t − X s ∼ N ( 0 , ( t − s ) I d ) です。成分ごとの増分ベクトルは独立な Gauss ベクトル なので全体が Gauss ベクトルで、異なる時間区間どうしの共分散が 0 0 0 だからです(命題 5.3[特性関数と Gauss ベクトル] )。
次の補題は、過去の情報で決まる変数 Y Y Y と、過去と独立な変数 Z Z Z があるとき、Y Y Y を固定した値とみなして Z Z Z だけについて平均してよい、という主張です。直後の Markov 性の証明で使います。
補題 2.1 (独立な変数の凍結 ) ✓
G ⊂ F \mathcal{G}\subset\mathcal{F} G ⊂ F を部分 σ-加法族、Y Y Y を G \mathcal{G} G -可測な R d \mathbb{R}^d R d 値確率変数、Z Z Z を G \mathcal{G} G と独立な可測空間 ( M , M ) (M,\mathcal{M}) ( M , M ) 値確率変数、Φ : R d × M → [ 0 , ∞ ] \Phi:\mathbb{R}^d\times M\to[0,\infty] Φ : R d × M → [ 0 , ∞ ] を B ( R d ) ⊗ M \mathcal{B}(\mathbb{R}^d)\otimes\mathcal{M} B ( R d ) ⊗ M -可測とする。g ( y ) = E [ Φ ( y , Z ) ] g(y)=E[\Phi(y,Z)] g ( y ) = E [ Φ ( y , Z )] とおくと g g g は可測で、E [ Φ ( Y , Z ) ∣ G ] = g ( Y ) E[\Phi(Y,Z)\mid\mathcal{G}]=g(Y) E [ Φ ( Y , Z ) ∣ G ] = g ( Y ) 。Φ \Phi Φ が有界な場合も同じ。
証明(補題 2.1) Φ = 1 A × B \Phi=\mathbf 1_{A\times B} Φ = 1 A × B のとき g ( y ) = 1 A ( y ) P ( Z ∈ B ) g(y)=\mathbf 1_A(y)P(Z\in B) g ( y ) = 1 A ( y ) P ( Z ∈ B ) で、G ∈ G G\in\mathcal{G} G ∈ G について 1 G 1 A ( Y ) \mathbf 1_G\mathbf 1_A(Y) 1 G 1 A ( Y ) は G \mathcal{G} G -可測、1 B ( Z ) \mathbf 1_B(Z) 1 B ( Z ) はそれと独立なので E [ 1 G 1 A ( Y ) 1 B ( Z ) ] = E [ 1 G 1 A ( Y ) ] P ( Z ∈ B ) = E [ 1 G g ( Y ) ] E[\mathbf 1_G\mathbf 1_A(Y)\mathbf 1_B(Z)]=E[\mathbf 1_G\mathbf 1_A(Y)]P(Z\in B)=E[\mathbf 1_Gg(Y)] E [ 1 G 1 A ( Y ) 1 B ( Z )] = E [ 1 G 1 A ( Y )] P ( Z ∈ B ) = E [ 1 G g ( Y )] です。等式 E [ 1 G 1 C ( Y , Z ) ] = E [ 1 G g C ( Y ) ] E[\mathbf 1_G\mathbf 1_C(Y,Z)]=E[\mathbf 1_Gg_C(Y)] E [ 1 G 1 C ( Y , Z )] = E [ 1 G g C ( Y )] (g C ( y ) = P ( ( y , Z ) ∈ C ) g_C(y)=P((y,Z)\in C) g C ( y ) = P (( y , Z ) ∈ C ) )が成り立つ C C C の全体は λ-系で、長方形(π-系)を含むので、定理 4.4[測度空間と Carathéodory の拡張定理] によりすべての C ∈ B ( R d ) ⊗ M C\in\mathcal{B}(\mathbb{R}^d)\otimes\mathcal{M} C ∈ B ( R d ) ⊗ M で成り立ちます。g C g_C g C の可測性は 補題 3.1[直積測度と Fubini–Tonelli の定理] から従います。線形性で非負単関数へ、単調収束定理で非負可測関数へ広がります。有界な Φ \Phi Φ は正部分と負部分に分けます。
∎
定理 2.2 (Brown 運動の Markov 性 ) ✓
F : W d → [ 0 , ∞ ] F:W^d\to[0,\infty] F : W d → [ 0 , ∞ ] (または有界)を可測関数、s ≥ 0 s\ge0 s ≥ 0 、x ∈ R d x\in\mathbb{R}^d x ∈ R d とする。P x P_x P x -ほとんど確実に
E x [ F ∘ θ s ∣ F s ] = E X s [ F ] , E_x\bigl[F\circ\theta_s\mid\mathcal{F}_s\bigr]=E_{X_s}[F], E x [ F ∘ θ s ∣ F s ] = E X s [ F ] , ここで右辺は関数 y ↦ E y [ F ] y\mapsto E_y[F] y ↦ E y [ F ] に y = X s y=X_s y = X s を代入したものである。
証明(定理 2.2) P x P_x P x の下で W t = X s + t − X s W_t=X_{s+t}-X_s W t = X s + t − X s とおくと、θ s X = X s + W \theta_sX=X_s+W θ s X = X s + W です。W W W は連続で W 0 = 0 W_0=0 W 0 = 0 、増分は X X X の時刻 s s s 以降の増分なので独立かつ N ( 0 , ( t − u ) I d ) N(0,(t-u)I_d) N ( 0 , ( t − u ) I d ) で、d d d 次元標準 Brown 運動です。定理 7.2[Kolmogorov の拡張定理と連続性定理] の一意性から W W W の法則は W \mathbb{W} W です。
W W W が F s \mathcal{F}_s F s と独立であることを示します。r 1 < ⋯ < r k ≤ s r_1<\cdots<r_k\le s r 1 < ⋯ < r k ≤ s と 0 < u 1 < ⋯ < u m 0<u_1<\cdots<u_m 0 < u 1 < ⋯ < u m について、( X r 1 , … , X r k ) (X_{r_1},\ldots,X_{r_k}) ( X r 1 , … , X r k ) は X 0 = x X_0=x X 0 = x と時刻 s s s 以前の増分の関数、( W u 1 , … , W u m ) (W_{u_1},\ldots,W_{u_m}) ( W u 1 , … , W u m ) は時刻 s s s 以降の増分の関数なので、両者は独立です。{ ( X r i ) i ∈ B } \{(X_{r_i})_i\in B\} {( X r i ) i ∈ B } の形の事象全体と { ( W u j ) j ∈ B ′ } \{(W_{u_j})_j\in B'\} {( W u j ) j ∈ B ′ } の形の事象全体はそれぞれ π-系で、F s \mathcal{F}_s F s と σ ( W ) \sigma(W) σ ( W ) を生成します。一方の π-系の元 A A A を固定すると P ( A ∩ B ) = P ( A ) P ( B ) P(A\cap B)=P(A)P(B) P ( A ∩ B ) = P ( A ) P ( B ) を満たす B B B は λ-系なので、定理 4.4[測度空間と Carathéodory の拡張定理] を二回使えば F s \mathcal{F}_s F s と σ ( W ) \sigma(W) σ ( W ) は独立です。
Φ ( y , w ) = F ( y + w ) \Phi(y,w)=F(y+w) Φ ( y , w ) = F ( y + w ) は連続写像 ( y , w ) ↦ y + w (y,w)\mapsto y+w ( y , w ) ↦ y + w と F F F の合成なので可測です(R d × W d \mathbb{R}^d\times W^d R d × W d は可分距離空間なので、ボレル集合族は直積 σ-加法族に等しい。補題 2.2[Kolmogorov の拡張定理と連続性定理] と同じ議論)。補題 2.1 を Y = X s Y=X_s Y = X s 、Z = W Z=W Z = W に使うと E x [ F ( X s + W ) ∣ F s ] = g ( X s ) E_x[F(X_s+W)\mid\mathcal{F}_s]=g(X_s) E x [ F ( X s + W ) ∣ F s ] = g ( X s ) 、g ( y ) = ∫ F ( y + w ) W ( d w ) = E y [ F ] g(y)=\int F(y+w)\,\mathbb{W}(dw)=E_y[F] g ( y ) = ∫ F ( y + w ) W ( d w ) = E y [ F ] です。
∎
以下 E E E は局所コンパクト な可分距離空間、C 0 ( E ) C_0(E) C 0 ( E ) は無限遠で消える連続関数の空間(一様ノルム )です(Riesz–Markov–角谷の表現定理 の記号)。
Feller 半群は、P t f ( x ) = E x [ f ( X t ) ] P_tf(x)=E_x[f(X_t)] P t f ( x ) = E x [ f ( X t )] (出発点 x x x から時間 t t t 後の f f f の平均)を作る操作 f ↦ P t f f\mapsto P_tf f ↦ P t f の性質を、過程そのものから切り離して取り出したものです。
定義 3.1 (Feller 半群 ) ↗
C 0 ( E ) C_0(E) C 0 ( E ) 上の線形作用素の族 ( T t ) t ≥ 0 (T_t)_{t\ge0} ( T t ) t ≥ 0 が Feller 半群 であるとは、次を満たすことをいう。
T 0 = I T_0=I T 0 = I 、T t + s = T t T s T_{t+s}=T_tT_s T t + s = T t T s (s , t ≥ 0 s,t\ge0 s , t ≥ 0 )。
(正値性・劣 Markov 性)0 ≤ f ≤ 1 0\le f\le1 0 ≤ f ≤ 1 ならば 0 ≤ T t f ≤ 1 0\le T_tf\le1 0 ≤ T t f ≤ 1 。
(強連続性)すべての f ∈ C 0 ( E ) f\in C_0(E) f ∈ C 0 ( E ) で ∥ T t f − f ∥ → 0 \|T_tf-f\|\to0 ∥ T t f − f ∥ → 0 (t ↓ 0 t\downarrow0 t ↓ 0 )。
その生成子 A A A を 定義 4.1[1 径数ユニタリ群と強連続半群] と同じく、D ( A ) = { f ∣ lim t ↓ 0 ( T t f − f ) / t D(A)=\{f\mid\lim_{t\downarrow0}(T_tf-f)/t D ( A ) = { f ∣ lim t ↓ 0 ( T t f − f ) / t が C 0 ( E ) C_0(E) C 0 ( E ) で存在} \} } 、A f Af A f をその極限で定める。
Riesz–Markov–角谷の定理(定理 4.1[Radon 測度と Riesz–Markov–角谷の表現定理] )を f ↦ T t f ( x ) f\mapsto T_tf(x) f ↦ T t f ( x ) に使うと、全質量 1 1 1 以下の測度 P t ( x , d y ) P_t(x,dy) P t ( x , d y ) で T t f ( x ) = ∫ f ( y ) P t ( x , d y ) T_tf(x)=\int f(y)P_t(x,dy) T t f ( x ) = ∫ f ( y ) P t ( x , d y ) と書けます。これが推移確率で、P t ( x , E ) < 1 P_t(x,E)<1 P t ( x , E ) < 1 は「時刻 t t t までに消滅(無限遠へ脱出)する確率」です。
次の命題の (3) の直感は、f f f が最大になる点では、平均をとっても値は上がらない(T t f ( x 0 ) ≤ f ( x 0 ) T_tf(x_0)\le f(x_0) T t f ( x 0 ) ≤ f ( x 0 ) )ので、その点での時間変化の割合 A f ( x 0 ) Af(x_0) A f ( x 0 ) は正になれない、ということです。
命題 3.2 (Feller 半群の基本性質 ) ✓
( T t ) (T_t) ( T t ) を Feller 半群、A A A をその生成子とする。
∥ T t ∥ ≤ 1 \|T_t\|\le1 ∥ T t ∥ ≤ 1 で、t ↦ T t f t\mapsto T_tf t ↦ T t f はすべての t ≥ 0 t\ge0 t ≥ 0 で連続である。
D ( A ) D(A) D ( A ) は稠密、A A A は閉作用素で、f ∈ D ( A ) f\in D(A) f ∈ D ( A ) なら T t f ∈ D ( A ) T_tf\in D(A) T t f ∈ D ( A ) 、d d t T t f = A T t f = T t A f \frac{d}{dt}T_tf=AT_tf=T_tAf d t d T t f = A T t f = T t A f 、T t f − f = ∫ 0 t T s A f d s T_tf-f=\int_0^tT_sAf\,ds T t f − f = ∫ 0 t T s A f d s 。λ > 0 \lambda>0 λ > 0 なら ( λ − A ) − 1 f = ∫ 0 ∞ e − λ t T t f d t (\lambda-A)^{-1}f=\int_0^\infty e^{-\lambda t}T_tf\,dt ( λ − A ) − 1 f = ∫ 0 ∞ e − λ t T t f d t 、∥ ( λ − A ) − 1 ∥ ≤ 1 / λ \|(\lambda-A)^{-1}\|\le1/\lambda ∥ ( λ − A ) − 1 ∥ ≤ 1/ λ 。
(正の最大値原理)f ∈ D ( A ) f\in D(A) f ∈ D ( A ) が x 0 x_0 x 0 で最大値 f ( x 0 ) ≥ 0 f(x_0)\ge0 f ( x 0 ) ≥ 0 をとるならば A f ( x 0 ) ≤ 0 Af(x_0)\le0 A f ( x 0 ) ≤ 0 。
(散逸性)f ∈ D ( A ) f\in D(A) f ∈ D ( A ) 、λ > 0 \lambda>0 λ > 0 ならば ∥ ( λ − A ) f ∥ ≥ λ ∥ f ∥ \|(\lambda-A)f\|\ge\lambda\|f\| ∥ ( λ − A ) f ∥ ≥ λ ∥ f ∥ 。
証明(命題 3.2) (1) ∥ f ∥ ≤ 1 \|f\|\le1 ∥ f ∥ ≤ 1 とし f = f + − f − f=f^+-f^- f = f + − f − と分けると、f ± f^\pm f ± は C 0 ( E ) C_0(E) C 0 ( E ) の元で値は [ 0 , 1 ] [0,1] [ 0 , 1 ] です。条件 2 から ∣ T t f ∣ ≤ T t f + + T t f − = T t ∣ f ∣ ≤ 1 |T_tf|\le T_tf^++T_tf^-=T_t|f|\le1 ∣ T t f ∣ ≤ T t f + + T t f − = T t ∣ f ∣ ≤ 1 。連続性は ∥ T t + h f − T t f ∥ = ∥ T t ( T h f − f ) ∥ ≤ ∥ T h f − f ∥ \|T_{t+h}f-T_tf\|=\|T_t(T_hf-f)\|\le\|T_hf-f\| ∥ T t + h f − T t f ∥ = ∥ T t ( T h f − f ) ∥ ≤ ∥ T h f − f ∥ と、h < 0 h<0 h < 0 の側も ∥ T t f − T t − ∣ h ∣ f ∥ ≤ ∥ T ∣ h ∣ f − f ∥ \|T_{t}f-T_{t-|h|}f\|\le\|T_{|h|}f-f\| ∥ T t f − T t − ∣ h ∣ f ∥ ≤ ∥ T ∣ h ∣ f − f ∥ から従います。
(2) ( T t ) (T_t) ( T t ) は縮小半群なので、命題 4.2[1 径数ユニタリ群と強連続半群] を M = 1 M=1 M = 1 、ω = 0 \omega=0 ω = 0 で使えます(その (3) の証明は ∥ T ( t ) ∥ ≤ M e ω t \|T(t)\|\le Me^{\omega t} ∥ T ( t ) ∥ ≤ M e ω t を満たす任意の組で通用します)。∥ ∫ 0 ∞ e − λ t T t f d t ∥ ≤ ∫ 0 ∞ e − λ t d t ∥ f ∥ = ∥ f ∥ / λ \|\int_0^\infty e^{-\lambda t}T_tf\,dt\|\le\int_0^\infty e^{-\lambda t}dt\,\|f\|=\|f\|/\lambda ∥ ∫ 0 ∞ e − λ t T t f d t ∥ ≤ ∫ 0 ∞ e − λ t d t ∥ f ∥ = ∥ f ∥/ λ です。
(3) 正値性から T t f ≤ T t f + T_tf\le T_tf^+ T t f ≤ T t f + で、f + / f ( x 0 ) f^+/f(x_0) f + / f ( x 0 ) (f ( x 0 ) > 0 f(x_0)>0 f ( x 0 ) > 0 のとき)は値が [ 0 , 1 ] [0,1] [ 0 , 1 ] なので T t f + ≤ f ( x 0 ) T_tf^+\le f(x_0) T t f + ≤ f ( x 0 ) 。f ( x 0 ) = 0 f(x_0)=0 f ( x 0 ) = 0 なら f + = 0 f^+=0 f + = 0 です。いずれも T t f ( x 0 ) ≤ f ( x 0 ) T_tf(x_0)\le f(x_0) T t f ( x 0 ) ≤ f ( x 0 ) で、( T t f ( x 0 ) − f ( x 0 ) ) / t ≤ 0 (T_tf(x_0)-f(x_0))/t\le0 ( T t f ( x 0 ) − f ( x 0 )) / t ≤ 0 の極限をとれば A f ( x 0 ) ≤ 0 Af(x_0)\le0 A f ( x 0 ) ≤ 0 。
(4) f ≠ 0 f\ne0 f = 0 とし、∣ f ∣ |f| ∣ f ∣ は C 0 ( E ) C_0(E) C 0 ( E ) の元なので最大値をとる点 x 0 x_0 x 0 があります。− f -f − f に取り替えて f ( x 0 ) = ∥ f ∥ f(x_0)=\|f\| f ( x 0 ) = ∥ f ∥ とすると、(3) から ( λ f − A f ) ( x 0 ) ≥ λ ∥ f ∥ (\lambda f-Af)(x_0)\ge\lambda\|f\| ( λ f − A f ) ( x 0 ) ≥ λ ∥ f ∥ 。
∎
逆に、稠密な定義域をもち正の最大値原理を満たす作用素 A A A について、ある λ > 0 \lambda>0 λ > 0 で λ − A \lambda-A λ − A の値域が稠密なら、A A A の閉包は Feller 半群を生成します(Hille–吉田の定理 定理 4.3[1 径数ユニタリ群と強連続半群] の Feller 版。Ethier–Kurtz の第 4 章。本稿では使いません)。Feller 半群からは、右連続で左極限をもつ道の Markov 過程が作れることも知られています。本稿で扱う過程はすべて Brown 運動から作るので、この一般論には立ち入りません。
半群が生成子で決まることは、次の初等的な補題で示せます。
補題 3.3 (右微分が 0 の連続関数 ) ✓
F : [ 0 , t ] → R F:[0,t]\to\mathbb{R} F : [ 0 , t ] → R が連続で、各 s ∈ [ 0 , t ) s\in[0,t) s ∈ [ 0 , t ) で右微分 D + F ( s ) = lim h ↓ 0 ( F ( s + h ) − F ( s ) ) / h D^+F(s)=\lim_{h\downarrow0}(F(s+h)-F(s))/h D + F ( s ) = lim h ↓ 0 ( F ( s + h ) − F ( s )) / h が存在して 0 0 0 ならば、F F F は定数である。
証明(補題 3.3) ε > 0 \varepsilon>0 ε > 0 を固定し、ϕ ( s ) = F ( s ) − F ( 0 ) − ε s \phi(s)=F(s)-F(0)-\varepsilon s ϕ ( s ) = F ( s ) − F ( 0 ) − ε s 、s ∗ = sup { s ∈ [ 0 , t ] ∣ ϕ ≤ 0 on [ 0 , s ] } s^*=\sup\{s\in[0,t]\mid\phi\le0\text{ on }[0,s]\} s ∗ = sup { s ∈ [ 0 , t ] ∣ ϕ ≤ 0 on [ 0 , s ]} とします。ϕ \phi ϕ は連続なので ϕ ( s ∗ ) ≤ 0 \phi(s^*)\le0 ϕ ( s ∗ ) ≤ 0 。s ∗ < t s^*<t s ∗ < t なら、D + ϕ ( s ∗ ) = − ε < 0 D^+\phi(s^*)=-\varepsilon<0 D + ϕ ( s ∗ ) = − ε < 0 から s ∗ s^* s ∗ の右の近くで ϕ ( s ) < ϕ ( s ∗ ) ≤ 0 \phi(s)<\phi(s^*)\le0 ϕ ( s ) < ϕ ( s ∗ ) ≤ 0 となり、s ∗ s^* s ∗ の定義に反します。よって F ( s ) − F ( 0 ) ≤ ε s F(s)-F(0)\le\varepsilon s F ( s ) − F ( 0 ) ≤ ε s 。− F -F − F にも同じ議論を使い、ε → 0 \varepsilon\to0 ε → 0 とすれば F ( s ) = F ( 0 ) F(s)=F(0) F ( s ) = F ( 0 ) です。
∎
定理 3.4 (半群は生成子で決まる ) ✓
( T t ) (T_t) ( T t ) 、( S t ) (S_t) ( S t ) を同じ生成子 A A A (定義域を含めて同じ)をもつ C 0 ( E ) C_0(E) C 0 ( E ) 上の Feller 半群とすると、すべての t t t で T t = S t T_t=S_t T t = S t 。
証明(定理 3.4) f ∈ D ( A ) f\in D(A) f ∈ D ( A ) 、t > 0 t>0 t > 0 、x ∈ E x\in E x ∈ E を固定し、F ( s ) = ( T t − s S s f ) ( x ) F(s)=(T_{t-s}S_sf)(x) F ( s ) = ( T t − s S s f ) ( x ) (0 ≤ s ≤ t 0\le s\le t 0 ≤ s ≤ t )とします。F F F は 命題 3.2 (1) と縮小性から連続です。0 < h < t − s 0<h<t-s 0 < h < t − s について
T t − s − h S s + h f − T t − s S s f h = T t − s − h ( S s + h f − S s f h − A S s f ) + T t − s − h A S s f − T t − s S s f − T t − s − h S s f h . \frac{T_{t-s-h}S_{s+h}f-T_{t-s}S_sf}{h}=T_{t-s-h}\Bigl(\frac{S_{s+h}f-S_sf}{h}-AS_sf\Bigr)+T_{t-s-h}AS_sf-\frac{T_{t-s}S_sf-T_{t-s-h}S_sf}{h}. h T t − s − h S s + h f − T t − s S s f = T t − s − h ( h S s + h f − S s f − A S s f ) + T t − s − h A S s f − h T t − s S s f − T t − s − h S s f . 第 1 項はノルムが括弧の中のノルム以下で 0 0 0 に収束します。S s f ∈ D ( A ) S_sf\in D(A) S s f ∈ D ( A ) なので 命題 3.2 (2) から r ↦ T r S s f r\mapsto T_rS_sf r ↦ T r S s f は微分 T r A S s f T_rAS_sf T r A S s f をもち、第 3 項は T t − s A S s f T_{t-s}AS_sf T t − s A S s f に収束します。第 2 項も同じ極限に収束するので、D + F ( s ) = 0 D^+F(s)=0 D + F ( s ) = 0 です。補題 3.3 から T t f ( x ) = F ( 0 ) = F ( t ) = S t f ( x ) T_tf(x)=F(0)=F(t)=S_tf(x) T t f ( x ) = F ( 0 ) = F ( t ) = S t f ( x ) 。D ( A ) D(A) D ( A ) は稠密で T t T_t T t 、S t S_t S t は有界なので T t = S t T_t=S_t T t = S t です。
∎
定理 3.5 (Brown 運動の Feller 半群 ) ✓
P t f ( x ) = E x [ f ( X t ) ] = ∫ R d p t ( x − y ) f ( y ) d y P_tf(x)=E_x[f(X_t)]=\int_{\mathbb{R}^d}p_t(x-y)f(y)\,dy P t f ( x ) = E x [ f ( X t )] = ∫ R d p t ( x − y ) f ( y ) d y 、p t ( z ) = ( 2 π t ) − d / 2 e − ∣ z ∣ 2 / ( 2 t ) p_t(z)=(2\pi t)^{-d/2}e^{-|z|^2/(2t)} p t ( z ) = ( 2 π t ) − d /2 e − ∣ z ∣ 2 / ( 2 t ) は C 0 ( R d ) C_0(\mathbb{R}^d) C 0 ( R d ) 上の Feller 半群である。その生成子 A A A は、f f f と 2 階までの偏導関数がすべて C 0 ( R d ) C_0(\mathbb{R}^d) C 0 ( R d ) に属する関数の空間 C 0 2 ( R d ) C_0^2(\mathbb{R}^d) C 0 2 ( R d ) を定義域に含み、そこで A f = 1 2 Δ f Af=\frac12\Delta f A f = 2 1 Δ f 。
証明(定理 3.5) C 0 C_0 C 0 を保つこと。 B = X − x B=X-x B = X − x は P x P_x P x の下で標準 Brown 運動なので P t f ( x ) = E 0 [ f ( x + X t ) ] P_tf(x)=E_0[f(x+X_t)] P t f ( x ) = E 0 [ f ( x + X t )] 。x n → x x_n\to x x n → x なら f ( x n + X t ) → f ( x + X t ) f(x_n+X_t)\to f(x+X_t) f ( x n + X t ) → f ( x + X t ) 、∣ x ∣ → ∞ |x|\to\infty ∣ x ∣ → ∞ なら f ( x + X t ) → 0 f(x+X_t)\to0 f ( x + X t ) → 0 で、どちらも ∥ f ∥ \|f\| ∥ f ∥ で押さえられるので、優収束定理から P t f P_tf P t f は連続で無限遠で消えます。
半群性。 Markov 性(定理 2.2 )を F = f ( X s ) F=f(X_s) F = f ( X s ) に使うと、P t + s f ( x ) = E x [ E x [ f ( X s ) ∘ θ t ∣ F t ] ] = E x [ P s f ( X t ) ] = P t P s f ( x ) P_{t+s}f(x)=E_x\bigl[E_x[f(X_s)\circ\theta_t\mid\mathcal{F}_t]\bigr]=E_x[P_sf(X_t)]=P_tP_sf(x) P t + s f ( x ) = E x [ E x [ f ( X s ) ∘ θ t ∣ F t ] ] = E x [ P s f ( X t )] = P t P s f ( x ) 。正値性と P t 1 = 1 P_t1=1 P t 1 = 1 は明らかです。
強連続性。 f f f は一様連続(C 0 C_0 C 0 の関数)なので、連続度 ω f ( δ ) = sup ∣ y ∣ ≤ δ sup x ∣ f ( x + y ) − f ( x ) ∣ → 0 \omega_f(\delta)=\sup_{|y|\le\delta}\sup_x|f(x+y)-f(x)|\to0 ω f ( δ ) = sup ∣ y ∣ ≤ δ sup x ∣ f ( x + y ) − f ( x ) ∣ → 0 。∣ P t f ( x ) − f ( x ) ∣ ≤ E 0 ∣ f ( x + X t ) − f ( x ) ∣ ≤ ω f ( δ ) + 2 ∥ f ∥ P 0 ( ∣ X t ∣ > δ ) |P_tf(x)-f(x)|\le E_0|f(x+X_t)-f(x)|\le\omega_f(\delta)+2\|f\|P_0(|X_t|>\delta) ∣ P t f ( x ) − f ( x ) ∣ ≤ E 0 ∣ f ( x + X t ) − f ( x ) ∣ ≤ ω f ( δ ) + 2∥ f ∥ P 0 ( ∣ X t ∣ > δ ) で、右辺は x x x によらず t → 0 t\to0 t → 0 で ω f ( δ ) \omega_f(\delta) ω f ( δ ) に収束し、δ \delta δ は任意です。
生成子。 f ∈ C 0 2 f\in C_0^2 f ∈ C 0 2 とし、Hesse 行列 ∇ 2 f \nabla^2f ∇ 2 f の連続度を ω ( δ ) \omega(\delta) ω ( δ ) とします(各成分が C 0 C_0 C 0 なので一様連続)。Taylor の定理から f ( x + y ) = f ( x ) + ∇ f ( x ) ⋅ y + 1 2 y T ∇ 2 f ( x ) y + r ( x , y ) f(x+y)=f(x)+\nabla f(x)\cdot y+\frac12y^{\mathsf T}\nabla^2f(x)y+r(x,y) f ( x + y ) = f ( x ) + ∇ f ( x ) ⋅ y + 2 1 y T ∇ 2 f ( x ) y + r ( x , y ) 、∣ r ( x , y ) ∣ ≤ 1 2 ∣ y ∣ 2 ω ( ∣ y ∣ ) ⋅ d |r(x,y)|\le\frac12|y|^2\omega(|y|)\cdot d ∣ r ( x , y ) ∣ ≤ 2 1 ∣ y ∣ 2 ω ( ∣ y ∣ ) ⋅ d です(剰余は中間点の Hesse 行列との差で、行列ノルムを成分の最大値の d d d 倍で押さえた)。y = X t y=X_t y = X t として E 0 E_0 E 0 をとると、E 0 [ X t ] = 0 E_0[X_t]=0 E 0 [ X t ] = 0 、E 0 [ X t i X t j ] = t δ i j E_0[X^i_tX^j_t]=t\delta_{ij} E 0 [ X t i X t j ] = t δ ij から
∣ P t f ( x ) − f ( x ) t − 1 2 Δ f ( x ) ∣ ≤ d 2 t E 0 [ ∣ X t ∣ 2 ω ( ∣ X t ∣ ) ] ≤ d 2 ( ω ( δ ) d + 2 ∥ ∇ 2 f ∥ max t E 0 [ ∣ X t ∣ 2 ; ∣ X t ∣ > δ ] ) . \Bigl|\frac{P_tf(x)-f(x)}{t}-\frac12\Delta f(x)\Bigr|\le\frac{d}{2t}E_0\bigl[|X_t|^2\omega(|X_t|)\bigr]\le\frac d2\Bigl(\omega(\delta)\,d+\frac{2\|\nabla^2f\|_{\max}}{t}E_0\bigl[|X_t|^2;|X_t|>\delta\bigr]\Bigr). t P t f ( x ) − f ( x ) − 2 1 Δ f ( x ) ≤ 2 t d E 0 [ ∣ X t ∣ 2 ω ( ∣ X t ∣ ) ] ≤ 2 d ( ω ( δ ) d + t 2∥ ∇ 2 f ∥ m a x E 0 [ ∣ X t ∣ 2 ; ∣ X t ∣ > δ ] ) . Cauchy–Schwarz の不等式で E 0 [ ∣ X t ∣ 2 ; ∣ X t ∣ > δ ] ≤ ( E 0 ∣ X t ∣ 4 ) 1 / 2 P 0 ( ∣ X t ∣ > δ ) 1 / 2 ≤ C t ⋅ ( d t / δ 2 ) 1 / 2 E_0[|X_t|^2;|X_t|>\delta]\le(E_0|X_t|^4)^{1/2}P_0(|X_t|>\delta)^{1/2}\le Ct\cdot(dt/\delta^2)^{1/2} E 0 [ ∣ X t ∣ 2 ; ∣ X t ∣ > δ ] ≤ ( E 0 ∣ X t ∣ 4 ) 1/2 P 0 ( ∣ X t ∣ > δ ) 1/2 ≤ C t ⋅ ( d t / δ 2 ) 1/2 なので、最後の項は t → 0 t\to0 t → 0 で 0 0 0 です。δ \delta δ は任意なので、x x x について一様に ( P t f − f ) / t → 1 2 Δ f (P_tf-f)/t\to\frac12\Delta f ( P t f − f ) / t → 2 1 Δ f です。
∎
命題 3.6 (Dynkin のマルチンゲール ) ✓
f ∈ D ( A ) f\in D(A) f ∈ D ( A ) (A A A は Brown 運動の生成子)ならば、すべての x x x で
M t = f ( X t ) − f ( X 0 ) − ∫ 0 t A f ( X u ) d u M_t=f(X_t)-f(X_0)-\int_0^tAf(X_u)\,du M t = f ( X t ) − f ( X 0 ) − ∫ 0 t A f ( X u ) d u は P x P_x P x の下で ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールである。
証明(命題 3.6) M t M_t M t は有界な適合過程です。s < t s<t s < t 、G ∈ F s G\in\mathcal{F}_s G ∈ F s とします。Markov 性から E x [ f ( X t ) ; G ] = E x [ P t − s f ( X s ) ; G ] E_x[f(X_t);G]=E_x[P_{t-s}f(X_s);G] E x [ f ( X t ) ; G ] = E x [ P t − s f ( X s ) ; G ] 、Fubini の定理と Markov 性から E x [ ∫ s t A f ( X u ) d u ; G ] = ∫ s t E x [ P u − s A f ( X s ) ; G ] d u = E x [ ∫ 0 t − s P r A f ( X s ) d r ; G ] E_x[\int_s^tAf(X_u)du;G]=\int_s^tE_x[P_{u-s}Af(X_s);G]\,du=E_x[\int_0^{t-s}P_rAf(X_s)dr;G] E x [ ∫ s t A f ( X u ) d u ; G ] = ∫ s t E x [ P u − s A f ( X s ) ; G ] d u = E x [ ∫ 0 t − s P r A f ( X s ) d r ; G ] 。命題 3.2 (2) の P t − s f − f = ∫ 0 t − s P r A f d r P_{t-s}f-f=\int_0^{t-s}P_rAf\,dr P t − s f − f = ∫ 0 t − s P r A f d r を点 X s X_s X s で評価する(点での評価は C 0 C_0 C 0 上の連続線形汎関数)と、E x [ M t − M s ; G ] = 0 E_x[M_t-M_s;G]=0 E x [ M t − M s ; G ] = 0 を得ます。
∎
f ∈ C 0 2 f\in C^2_0 f ∈ C 0 2 のとき、これは伊藤の公式(定理 5.2[確率微分方程式] )の確率積分の項がマルチンゲール であることと同じ内容です。生成子は「f ( X t ) f(X_t) f ( X t ) からマルチンゲールを引いた残りのドリフト」として、過程の無限小の動きを表します。
可測関数 V : R d → R V:\mathbb{R}^d\to\mathbb{R} V : R d → R と有界可測関数 f f f に対し、期待値が意味をもつとき
T t V f ( x ) = E x [ e − ∫ 0 t V ( X s ) d s f ( X t ) ] T^V_tf(x)=E_x\Bigl[e^{-\int_0^tV(X_s)\,ds}f(X_t)\Bigr] T t V f ( x ) = E x [ e − ∫ 0 t V ( X s ) d s f ( X t ) ]
と書きます。これを Feynman–Kac 半群 と呼びます。まず半群性は Markov 性だけから出ます。
補題 4.1 (Feynman–Kac 半群の半群性 ) ✓
可測関数 V V V が、すべての t ≥ 0 t\ge0 t ≥ 0 で c V ( t ) = sup x E x [ e ∫ 0 t ∣ V ( X s ) ∣ d s ] < ∞ c_V(t)=\sup_xE_x\bigl[e^{\int_0^t|V(X_s)|ds}\bigr]<\infty c V ( t ) = sup x E x [ e ∫ 0 t ∣ V ( X s ) ∣ d s ] < ∞ を満たすとする。このとき有界可測な f f f について T t V f T^V_tf T t V f は well-defined な有界可測関数で ∣ T t V f ∣ ≤ c V ( t ) ∥ f ∥ |T^V_tf|\le c_V(t)\|f\| ∣ T t V f ∣ ≤ c V ( t ) ∥ f ∥ 、さらに T t + s V f = T t V ( T s V f ) T^V_{t+s}f=T^V_t(T^V_sf) T t + s V f = T t V ( T s V f ) 。
証明(補題 4.1) 評価は ∣ e − ∫ V f ( X t ) ∣ ≤ e ∫ ∣ V ∣ ∥ f ∥ |e^{-\int V}f(X_t)|\le e^{\int|V|}\|f\| ∣ e − ∫ V f ( X t ) ∣ ≤ e ∫ ∣ V ∣ ∥ f ∥ から従います。Φ ( w ) = e − ∫ 0 t V ( w ( s ) ) d s f ( w ( t ) ) \Phi(w)=e^{-\int_0^tV(w(s))ds}f(w(t)) Φ ( w ) = e − ∫ 0 t V ( w ( s )) d s f ( w ( t )) とおくと T t V f ( x ) = ∫ Φ ( x + w ) W ( d w ) T_t^Vf(x)=\int\Phi(x+w)\,\mathbb{W}(dw) T t V f ( x ) = ∫ Φ ( x + w ) W ( d w ) で、x x x についての可測性は Fubini の定理によります。f ≥ 0 f\ge0 f ≥ 0 として示せば十分です(一般には f = f + − f − f=f^+-f^- f = f + − f − )。F = e − ∫ 0 s V ( X r ) d r f ( X s ) F=e^{-\int_0^sV(X_r)dr}f(X_s) F = e − ∫ 0 s V ( X r ) d r f ( X s ) とおくと F ∘ θ t = e − ∫ t t + s V ( X r ) d r f ( X t + s ) F\circ\theta_t=e^{-\int_t^{t+s}V(X_r)dr}f(X_{t+s}) F ∘ θ t = e − ∫ t t + s V ( X r ) d r f ( X t + s ) です。e − ∫ 0 t V ( X r ) d r e^{-\int_0^tV(X_r)dr} e − ∫ 0 t V ( X r ) d r は非負の F t \mathcal{F}_t F t -可測変数なので、非負変数についての取り出し(定理 4.4[条件付き期待値] を単調収束定理で非負の場合に広げたもの)と 定理 2.2 から
T t + s V f ( x ) = E x [ e − ∫ 0 t V ( X r ) d r E x [ F ∘ θ t ∣ F t ] ] = E x [ e − ∫ 0 t V ( X r ) d r T s V f ( X t ) ] = T t V ( T s V f ) ( x ) . T^V_{t+s}f(x)=E_x\Bigl[e^{-\int_0^tV(X_r)dr}\,E_x[F\circ\theta_t\mid\mathcal{F}_t]\Bigr]=E_x\Bigl[e^{-\int_0^tV(X_r)dr}\,T^V_sf(X_t)\Bigr]=T^V_t(T^V_sf)(x). T t + s V f ( x ) = E x [ e − ∫ 0 t V ( X r ) d r E x [ F ∘ θ t ∣ F t ] ] = E x [ e − ∫ 0 t V ( X r ) d r T s V f ( X t ) ] = T t V ( T s V f ) ( x ) . ∎
定理 4.2 (Feynman–Kac 公式(有界一様連続なポテンシャル) ) ✓
V V V を R d \mathbb{R}^d R d 上の有界かつ一様連続な実数値関数とし、A A A を Brown 運動の生成子(定理 3.5 )とする。( T t V ) (T^V_t) ( T t V ) は C 0 ( R d ) C_0(\mathbb{R}^d) C 0 ( R d ) 上の強連続半群で ∥ T t V ∥ ≤ e t ∥ V ∥ \|T^V_t\|\le e^{t\|V\|} ∥ T t V ∥ ≤ e t ∥ V ∥ を満たし、その生成子は A − V A-V A − V (定義域 D ( A ) D(A) D ( A ) )である。すなわち T t V = e t ( A − V ) T^V_t=e^{t(A-V)} T t V = e t ( A − V ) であり、f ∈ C 0 2 f\in C^2_0 f ∈ C 0 2 なら u ( t , x ) = T t V f ( x ) u(t,x)=T^V_tf(x) u ( t , x ) = T t V f ( x ) は ∂ t u = 1 2 Δ u − V u \partial_tu=\frac12\Delta u-Vu ∂ t u = 2 1 Δ u − V u の解である。
証明は三段です。まず ( T t V ) (T^V_t) ( T t V ) が強連続半群であることを確かめ、次にその生成子が A − V A-V A − V を拡張することを計算し、最後に定義域がちょうど D ( A ) D(A) D ( A ) に一致することを、V V V を掛ける作用素が小さい摂動であることから示します。
証明(定理 4.2) 半群であること。 c V ( t ) ≤ e t ∥ V ∥ c_V(t)\le e^{t\|V\|} c V ( t ) ≤ e t ∥ V ∥ なので 補題 4.1 が使えます。f ∈ C 0 f\in C_0 f ∈ C 0 なら、V V V 、f f f の連続性と優収束定理から T t V f T^V_tf T t V f は連続で、∣ T t V f ∣ ≤ e t ∥ V ∥ P t ∣ f ∣ |T^V_tf|\le e^{t\|V\|}P_t|f| ∣ T t V f ∣ ≤ e t ∥ V ∥ P t ∣ f ∣ から無限遠で消えます。∥ T t V f − P t f ∥ ≤ ∥ f ∥ ( e t ∥ V ∥ − 1 ) \|T^V_tf-P_tf\|\le\|f\|(e^{t\|V\|}-1) ∥ T t V f − P t f ∥ ≤ ∥ f ∥ ( e t ∥ V ∥ − 1 ) と 定理 3.5 から強連続です。
生成子が A − V A-V A − V を拡張すること。 f ∈ D ( A ) f\in D(A) f ∈ D ( A ) とします。a = ∫ 0 t V ( X s ) d s a=\int_0^tV(X_s)ds a = ∫ 0 t V ( X s ) d s とおくと ∣ a ∣ ≤ t ∥ V ∥ |a|\le t\|V\| ∣ a ∣ ≤ t ∥ V ∥ で、∣ e − a − 1 + a ∣ ≤ 1 2 a 2 e ∣ a ∣ |e^{-a}-1+a|\le\frac12a^2e^{|a|} ∣ e − a − 1 + a ∣ ≤ 2 1 a 2 e ∣ a ∣ なので
T t V f − f t − P t f − f t = − E x [ a t f ( X t ) ] + O ( t ) ( x について一様 ) . \frac{T^V_tf-f}{t}-\frac{P_tf-f}{t}=-E_x\Bigl[\frac at\,f(X_t)\Bigr]+O(t)\quad(\text{x について一様}). t T t V f − f − t P t f − f = − E x [ t a f ( X t ) ] + O ( t ) ( x について一様 ) . ∣ E x [ V ( X s ) f ( X t ) ] − V ( x ) f ( x ) ∣ ≤ ∥ f ∥ E x ∣ V ( X s ) − V ( x ) ∣ + ∥ V ∥ E x ∣ f ( X t ) − f ( x ) ∣ \bigl|E_x[V(X_s)f(X_t)]-V(x)f(x)\bigr|\le\|f\|E_x|V(X_s)-V(x)|+\|V\|E_x|f(X_t)-f(x)| E x [ V ( X s ) f ( X t )] − V ( x ) f ( x ) ≤ ∥ f ∥ E x ∣ V ( X s ) − V ( x ) ∣ + ∥ V ∥ E x ∣ f ( X t ) − f ( x ) ∣ で、定理 3.5 の強連続性の証明と同じく、V V V と f f f の一様連続性から右辺は s , t → 0 s,t\to0 s , t → 0 で x x x について一様に 0 0 0 に収束します。よって E x [ a t f ( X t ) ] = 1 t ∫ 0 t E x [ V ( X s ) f ( X t ) ] d s → V ( x ) f ( x ) E_x[\frac at f(X_t)]=\frac1t\int_0^tE_x[V(X_s)f(X_t)]ds\to V(x)f(x) E x [ t a f ( X t )] = t 1 ∫ 0 t E x [ V ( X s ) f ( X t )] d s → V ( x ) f ( x ) が一様に成り立ち、f ∈ D ( A V ) f\in D(A^V) f ∈ D ( A V ) 、A V f = A f − V f A^Vf=Af-Vf A V f = A f − V f です(A V A^V A V は ( T t V ) (T^V_t) ( T t V ) の生成子)。
定義域の一致。 λ > ∥ V ∥ \lambda>\|V\| λ > ∥ V ∥ とします。∥ ( λ − A ) − 1 ∥ ≤ 1 / λ \|(\lambda-A)^{-1}\|\le1/\lambda ∥ ( λ − A ) − 1 ∥ ≤ 1/ λ (命題 3.2 (2))なので、掛け算作用素 V V V との積は ∥ V ( λ − A ) − 1 ∥ < 1 \|V(\lambda-A)^{-1}\|<1 ∥ V ( λ − A ) − 1 ∥ < 1 で、I + V ( λ − A ) − 1 I+V(\lambda-A)^{-1} I + V ( λ − A ) − 1 は Neumann 級数で可逆です。よって λ − A + V = ( I + V ( λ − A ) − 1 ) ( λ − A ) \lambda-A+V=(I+V(\lambda-A)^{-1})(\lambda-A) λ − A + V = ( I + V ( λ − A ) − 1 ) ( λ − A ) は D ( A ) D(A) D ( A ) から C 0 C_0 C 0 への全単射です。一方 ∥ T t V ∥ ≤ e t ∥ V ∥ \|T^V_t\|\le e^{t\|V\|} ∥ T t V ∥ ≤ e t ∥ V ∥ なので、命題 4.2[1 径数ユニタリ群と強連続半群] (3) から λ − A V \lambda-A^V λ − A V は D ( A V ) D(A^V) D ( A V ) 上で単射です。g ∈ D ( A V ) g\in D(A^V) g ∈ D ( A V ) に対し ( λ − A + V ) g ′ = ( λ − A V ) g (\lambda-A+V)g'=(\lambda-A^V)g ( λ − A + V ) g ′ = ( λ − A V ) g となる g ′ ∈ D ( A ) g'\in D(A) g ′ ∈ D ( A ) を取ると、A V ⊃ A − V A^V\supset A-V A V ⊃ A − V から ( λ − A V ) ( g − g ′ ) = 0 (\lambda-A^V)(g-g')=0 ( λ − A V ) ( g − g ′ ) = 0 、よって g = g ′ ∈ D ( A ) g=g'\in D(A) g = g ′ ∈ D ( A ) です。
∎
V V V が有界連続なら、定理 2.1[Feynman–Kac 表示と場の模型の汎関数積分] は L 2 ( R d ) L^2(\mathbb{R}^d) L 2 ( R d ) 上で e − t H f = T t V f e^{-tH}f=T^V_tf e − t H f = T t V f (H = − 1 2 Δ + V H=-\frac12\Delta+V H = − 2 1 Δ + V )を示しています。定理 4.2 はその C 0 C_0 C 0 版で、Trotter 積公式を使わず、生成子の計算と定義域の一致だけで同じ等式を得ています。
V ( x ) = − 1 / ∣ x ∣ V(x)=-1/|x| V ( x ) = − 1/∣ x ∣ (d = 3 d=3 d = 3 )のような特異なポテンシャルでは、∫ 0 t V ( X s ) d s \int_0^tV(X_s)ds ∫ 0 t V ( X s ) d s が有限か、e − ∫ V e^{-\int V} e − ∫ V が可積分かがまず問題です。
定義 4.3 (Kato クラス(確率論的な定義) ) ↗
可測関数 V : R d → R V:\mathbb{R}^d\to\mathbb{R} V : R d → R が
lim t ↓ 0 sup x ∈ R d E x [ ∫ 0 t ∣ V ( X s ) ∣ d s ] = 0 \lim_{t\downarrow0}\sup_{x\in\mathbb{R}^d}E_x\Bigl[\int_0^t|V(X_s)|\,ds\Bigr]=0 t ↓ 0 lim x ∈ R d sup E x [ ∫ 0 t ∣ V ( X s ) ∣ d s ] = 0 を満たすとき、V V V は Kato クラス に属するという。
E x ∫ 0 t ∣ V ( X s ) ∣ d s = ∫ ∣ V ( y ) ∣ ∫ 0 t p s ( x − y ) d s d y E_x\int_0^t|V(X_s)|ds=\int|V(y)|\int_0^tp_s(x-y)\,ds\,dy E x ∫ 0 t ∣ V ( X s ) ∣ d s = ∫ ∣ V ( y ) ∣ ∫ 0 t p s ( x − y ) d s d y (Tonelli)なので、これは ∣ V ∣ |V| ∣ V ∣ と「短時間の Green 関数」の畳み込みが一様に小さいという条件です。d ≥ 3 d\ge3 d ≥ 3 では ∫ 0 ∞ p s ( z ) d s \int_0^\infty p_s(z)ds ∫ 0 ∞ p s ( z ) d s が c d ∣ z ∣ 2 − d c_d|z|^{2-d} c d ∣ z ∣ 2 − d になることから、この条件は Schrödinger 作用素 の 定義 4.1[シュレーディンガー作用素の自己共役性と束縛状態] と同値になります(Aizenman–Simon の定理。本稿では同値性を証明せず、確率論的な定義のほうを使います)。式に出てくる ∫ 0 t p s ( x − y ) d s \int_0^tp_s(x-y)\,ds ∫ 0 t p s ( x − y ) d s は、粒子が時間 t t t までに y y y の近くで過ごす時間の、単位体積あたりの平均を表します。
例 4.4 (Coulomb ポテンシャルは Kato クラス )
d = 3 d=3 d = 3 、V ( x ) = − 1 / ∣ x ∣ V(x)=-1/|x| V ( x ) = − 1/∣ x ∣ とします。Z ∼ N ( 0 , I 3 ) Z\sim N(0,I_3) Z ∼ N ( 0 , I 3 ) について E ∣ Z ∣ − 1 = ∫ 0 ∞ r − 1 ( 2 π ) − 3 / 2 e − r 2 / 2 4 π r 2 d r = 2 / π E|Z|^{-1}=\int_0^\infty r^{-1}(2\pi)^{-3/2}e^{-r^2/2}4\pi r^2dr=\sqrt{2/\pi} E ∣ Z ∣ − 1 = ∫ 0 ∞ r − 1 ( 2 π ) − 3/2 e − r 2 /2 4 π r 2 d r = 2/ π なので、E 0 ∣ X s ∣ − 1 = 2 / ( π s ) E_0|X_s|^{-1}=\sqrt{2/(\pi s)} E 0 ∣ X s ∣ − 1 = 2/ ( π s ) です。x ≠ 0 x\ne0 x = 0 では、∣ y ∣ − 1 |y|^{-1} ∣ y ∣ − 1 の球面平均が 1 / max ( ∣ x ∣ , r ) ≤ 1 / r 1/\max(|x|,r)\le1/r 1/ max ( ∣ x ∣ , r ) ≤ 1/ r であること(Newton の定理。Lieb–Loss の第 9 章)から E x ∣ X s ∣ − 1 ≤ E 0 ∣ X s ∣ − 1 E_x|X_s|^{-1}\le E_0|X_s|^{-1} E x ∣ X s ∣ − 1 ≤ E 0 ∣ X s ∣ − 1 。したがって
sup x E x ∫ 0 t d s ∣ X s ∣ ≤ ∫ 0 t 2 π s d s = 2 2 t π → 0. \sup_xE_x\int_0^t\frac{ds}{|X_s|}\le\int_0^t\sqrt{\frac{2}{\pi s}}\,ds=2\sqrt{\frac{2t}{\pi}}\to0 . x sup E x ∫ 0 t ∣ X s ∣ d s ≤ ∫ 0 t π s 2 d s = 2 π 2 t → 0. 一方 V = − 1 / ∣ x ∣ 2 V=-1/|x|^2 V = − 1/∣ x ∣ 2 では E 0 ∣ X s ∣ − 2 = 1 / s E_0|X_s|^{-2}=1/s E 0 ∣ X s ∣ − 2 = 1/ s で、∫ 0 t d s / s = ∞ \int_0^t ds/s=\infty ∫ 0 t d s / s = ∞ となり Kato クラスに入りません。
Coulomb ポテンシャルは原点付近で非有界なので、道ごとに ∫ ∣ V ∣ \int|V| ∫ ∣ V ∣ を押さえることはできず、e ∫ ∣ V ∣ e^{\int|V|} e ∫ ∣ V ∣ の平均が有限かどうかは自明ではありません。次の補題は、平均 E x ∫ W E_x\int W E x ∫ W が 1 より小さく押さえられれば、指数の平均も押さえられる、と述べます。証明は指数を ∑ n I n / n ! \sum_nI^n/n! ∑ n I n / n ! に展開し、n n n 重積分を帰納法で α n \alpha^n α n 以下と評価します。
補題 4.5 (Khasminskii の補題 ) ✓
W ≥ 0 W\ge0 W ≥ 0 を可測関数とし、t > 0 t>0 t > 0 について α = sup x E x ∫ 0 t W ( X s ) d s < 1 \alpha=\sup_xE_x\int_0^tW(X_s)\,ds<1 α = sup x E x ∫ 0 t W ( X s ) d s < 1 とする。このとき
sup x E x [ e ∫ 0 t W ( X s ) d s ] ≤ 1 1 − α . \sup_xE_x\Bigl[e^{\int_0^tW(X_s)ds}\Bigr]\le\frac1{1-\alpha}. x sup E x [ e ∫ 0 t W ( X s ) d s ] ≤ 1 − α 1 .
証明(補題 4.5) I = ∫ 0 t W ( X s ) d s I=\int_0^tW(X_s)ds I = ∫ 0 t W ( X s ) d s とします。被積分関数の対称性から
I n n ! = ∫ 0 < s 1 < ⋯ < s n < t W ( X s 1 ) ⋯ W ( X s n ) d s 1 ⋯ d s n . \frac{I^n}{n!}=\int_{0<s_1<\cdots<s_n<t}W(X_{s_1})\cdots W(X_{s_n})\,ds_1\cdots ds_n . n ! I n = ∫ 0 < s 1 < ⋯ < s n < t W ( X s 1 ) ⋯ W ( X s n ) d s 1 ⋯ d s n . g n ( r , y ) = E y ∫ 0 < s 1 < ⋯ < s n < r ∏ i W ( X s i ) d s g_n(r,y)=E_y\int_{0<s_1<\cdots<s_n<r}\prod_iW(X_{s_i})\,ds g n ( r , y ) = E y ∫ 0 < s 1 < ⋯ < s n < r ∏ i W ( X s i ) d s とおき、sup y g n ( t , y ) ≤ α n \sup_yg_n(t,y)\le\alpha^n sup y g n ( t , y ) ≤ α n を帰納法で示します。n = 1 n=1 n = 1 は仮定です。g n ( r , y ) g_n(r,y) g n ( r , y ) は r r r について非減少です。s 1 s_1 s 1 を外に出して内側の積分を Φ t − s 1 ∘ θ s 1 \Phi_{t-s_1}\circ\theta_{s_1} Φ t − s 1 ∘ θ s 1 (Φ r ( w ) = ∫ 0 < u 2 < ⋯ < u n < r ∏ W ( w ( u i ) ) d u \Phi_r(w)=\int_{0<u_2<\cdots<u_n<r}\prod W(w(u_i))du Φ r ( w ) = ∫ 0 < u 2 < ⋯ < u n < r ∏ W ( w ( u i )) d u )と書くと、Tonelli の定理と 定理 2.2 (非負の汎関数)から
g n ( t , x ) = ∫ 0 t E x [ W ( X s 1 ) g n − 1 ( t − s 1 , X s 1 ) ] d s 1 ≤ sup y g n − 1 ( t , y ) ⋅ E x ∫ 0 t W ( X s 1 ) d s 1 ≤ α n − 1 ⋅ α . g_n(t,x)=\int_0^tE_x\bigl[W(X_{s_1})\,g_{n-1}(t-s_1,X_{s_1})\bigr]ds_1\le\sup_yg_{n-1}(t,y)\cdot E_x\int_0^tW(X_{s_1})ds_1\le\alpha^{n-1}\cdot\alpha . g n ( t , x ) = ∫ 0 t E x [ W ( X s 1 ) g n − 1 ( t − s 1 , X s 1 ) ] d s 1 ≤ y sup g n − 1 ( t , y ) ⋅ E x ∫ 0 t W ( X s 1 ) d s 1 ≤ α n − 1 ⋅ α . 単調収束定理で和をとると E x e I = ∑ n E x [ I n / n ! ] ≤ ∑ n α n = 1 / ( 1 − α ) E_xe^I=\sum_n E_x[I^n/n!]\le\sum_n\alpha^n=1/(1-\alpha) E x e I = ∑ n E x [ I n / n !] ≤ ∑ n α n = 1/ ( 1 − α ) 。
∎
定理 4.6 (Kato クラスのポテンシャルの Feynman–Kac 半群 ) ✓
V V V を Kato クラスの関数とする。
定数 C , c C,c C , c があって sup x E x [ e ∫ 0 t ∣ V ( X s ) ∣ d s ] ≤ C e c t \sup_xE_x\bigl[e^{\int_0^t|V(X_s)|ds}\bigr]\le Ce^{ct} sup x E x [ e ∫ 0 t ∣ V ( X s ) ∣ d s ] ≤ C e c t (すべての t ≥ 0 t\ge0 t ≥ 0 )。特に P x P_x P x -ほとんど確実に ∫ 0 t ∣ V ( X s ) ∣ d s < ∞ \int_0^t|V(X_s)|ds<\infty ∫ 0 t ∣ V ( X s ) ∣ d s < ∞ 。
有界可測な f f f について T t V f T^V_tf T t V f は well-defined で ∥ T t V f ∥ ∞ ≤ C e c t ∥ f ∥ ∞ \|T^V_tf\|_\infty\le Ce^{ct}\|f\|_\infty ∥ T t V f ∥ ∞ ≤ C e c t ∥ f ∥ ∞ 、かつ T t + s V = T t V T s V T^V_{t+s}=T^V_tT^V_s T t + s V = T t V T s V 。
V n V_n V n が可測で ∣ V n ∣ ≤ ∣ V ∣ |V_n|\le|V| ∣ V n ∣ ≤ ∣ V ∣ 、ルベーグ測度についてほとんど至るところ V n → V V_n\to V V n → V ならば、すべての t > 0 t>0 t > 0 、x x x と有界可測な f f f について T t V n f ( x ) → T t V f ( x ) T^{V_n}_tf(x)\to T^V_tf(x) T t V n f ( x ) → T t V f ( x ) 。特に V V V の切断 V n = max ( − n , min ( V , n ) ) V_n=\max(-n,\min(V,n)) V n = max ( − n , min ( V , n )) について成り立つ。
証明(定理 4.6) (1) t 0 > 0 t_0>0 t 0 > 0 を sup x E x ∫ 0 t 0 ∣ V ( X s ) ∣ d s ≤ 1 2 \sup_xE_x\int_0^{t_0}|V(X_s)|ds\le\frac12 sup x E x ∫ 0 t 0 ∣ V ( X s ) ∣ d s ≤ 2 1 となるように取ると、補題 4.5 から sup x E x e ∫ 0 t 0 ∣ V ∣ ≤ 2 \sup_xE_xe^{\int_0^{t_0}|V|}\le2 sup x E x e ∫ 0 t 0 ∣ V ∣ ≤ 2 。補題 4.1 の証明と同じく Markov 性で区間 [ 0 , t 0 ] [0,t_0] [ 0 , t 0 ] 、[ t 0 , 2 t 0 ] [t_0,2t_0] [ t 0 , 2 t 0 ] 、… に分けると sup x E x e ∫ 0 k t 0 ∣ V ∣ ≤ 2 k \sup_xE_xe^{\int_0^{kt_0}|V|}\le2^k sup x E x e ∫ 0 k t 0 ∣ V ∣ ≤ 2 k 。k = ⌈ t / t 0 ⌉ k=\lceil t/t_0\rceil k = ⌈ t / t 0 ⌉ として C = 2 C=2 C = 2 、c = ( log 2 ) / t 0 c=(\log2)/t_0 c = ( log 2 ) / t 0 が取れます。
(2) (1) から c V ( t ) ≤ C e c t < ∞ c_V(t)\le Ce^{ct}<\infty c V ( t ) ≤ C e c t < ∞ なので 補題 4.1 です。
(3) 収束しない点の集合 N N N はルベーグ零集合です。s > 0 s>0 s > 0 では X s X_s X s の P x P_x P x の下の分布は密度 p s ( ⋅ − x ) p_s(\cdot-x) p s ( ⋅ − x ) をもつので、E x ∫ 0 t 1 N ( X s ) d s = 0 E_x\int_0^t\mathbf 1_N(X_s)ds=0 E x ∫ 0 t 1 N ( X s ) d s = 0 。よって P x P_x P x -ほとんど確実に、ほとんどすべての s ∈ [ 0 , t ] s\in[0,t] s ∈ [ 0 , t ] で X s ∉ N X_s\notin N X s ∈ / N で、そこでは V n ( X s ) → V ( X s ) V_n(X_s)\to V(X_s) V n ( X s ) → V ( X s ) 。(1) から ∫ 0 t ∣ V ( X s ) ∣ d s < ∞ \int_0^t|V(X_s)|ds<\infty ∫ 0 t ∣ V ( X s ) ∣ d s < ∞ も成り立つので、s s s についての優収束定理から ∫ 0 t V n ( X s ) d s → ∫ 0 t V ( X s ) d s \int_0^tV_n(X_s)ds\to\int_0^tV(X_s)ds ∫ 0 t V n ( X s ) d s → ∫ 0 t V ( X s ) d s (a.s.)。∣ e − ∫ V n f ( X t ) ∣ ≤ e ∫ ∣ V ∣ ∥ f ∥ |e^{-\int V_n}f(X_t)|\le e^{\int|V|}\|f\| ∣ e − ∫ V n f ( X t ) ∣ ≤ e ∫ ∣ V ∣ ∥ f ∥ は (1) により可積分なので、ω \omega ω についての優収束定理から結論を得ます。
∎
Feynman–Kac 半群は T t V 1 ≠ 1 T^V_t1\ne1 T t V 1 = 1 なので確率を保存しません。正の固有関数で「ゲージ変換」すると、確率を保存する半群に戻せます。
定義 5.1 (正値核の半群 ) ↗
可測空間 ( E , E ) (E,\mathcal{E}) ( E , E ) 上の核の族 ( Q t ) t ≥ 0 (Q_t)_{t\ge0} ( Q t ) t ≥ 0 (各 Q t ( x , ⋅ ) Q_t(x,\cdot) Q t ( x , ⋅ ) は有限測度、各 A A A で x ↦ Q t ( x , A ) x\mapsto Q_t(x,A) x ↦ Q t ( x , A ) は可測、sup x Q t ( x , E ) < ∞ \sup_xQ_t(x,E)<\infty sup x Q t ( x , E ) < ∞ )が Q 0 ( x , ⋅ ) = δ x Q_0(x,\cdot)=\delta_x Q 0 ( x , ⋅ ) = δ x と Chapman–Kolmogorov 方程式 Q t + s ( x , A ) = ∫ Q t ( x , d y ) Q s ( y , A ) Q_{t+s}(x,A)=\int Q_t(x,dy)Q_s(y,A) Q t + s ( x , A ) = ∫ Q t ( x , d y ) Q s ( y , A ) を満たすとき、正値核の半群 という。Q t f ( x ) = ∫ f ( y ) Q t ( x , d y ) Q_tf(x)=\int f(y)Q_t(x,dy) Q t f ( x ) = ∫ f ( y ) Q t ( x , d y ) と書く。Q t ( x , E ) = 1 Q_t(x,E)=1 Q t ( x , E ) = 1 がすべての t , x t,x t , x で成り立つとき Markov 半群 という。
V V V が Kato クラスなら T t V ( x , A ) = E x [ e − ∫ 0 t V ( X s ) d s ; X t ∈ A ] T^V_t(x,A)=E_x[e^{-\int_0^tV(X_s)ds};X_t\in A] T t V ( x , A ) = E x [ e − ∫ 0 t V ( X s ) d s ; X t ∈ A ] は正値核の半群です(定理 4.6 (2))。
h ( y ) / h ( x ) h(y)/h(x) h ( y ) / h ( x ) を掛けるのは、Q t h = e λ t h Q_th=e^{\lambda t}h Q t h = e λ t h のとき、e − λ t e^{-\lambda t} e − λ t と合わせて核の全質量がちょうど 1 1 1 になるようにするためです。確率を保存する半群に戻す操作と考えてください。
証明(定理 5.2) Q t h ( x , E ) = e − λ t Q t h ( x ) / h ( x ) = 1 Q^h_t(x,E)=e^{-\lambda t}Q_th(x)/h(x)=1 Q t h ( x , E ) = e − λ t Q t h ( x ) / h ( x ) = 1 。可測性は Q t ( x , ⋅ ) Q_t(x,\cdot) Q t ( x , ⋅ ) の性質から従います。Chapman–Kolmogorov 方程式は
∫ Q t h ( x , d y ) Q s h ( y , A ) = ∫ e − λ t h ( y ) h ( x ) Q t ( x , d y ) e − λ s Q s ( h 1 A ) ( y ) h ( y ) = e − λ ( t + s ) h ( x ) Q t + s ( h 1 A ) ( x ) = Q t + s h ( x , A ) . \int Q^h_t(x,dy)Q^h_s(y,A)=\int e^{-\lambda t}\frac{h(y)}{h(x)}Q_t(x,dy)\,e^{-\lambda s}\frac{Q_s(h\mathbf 1_A)(y)}{h(y)}=\frac{e^{-\lambda(t+s)}}{h(x)}Q_{t+s}(h\mathbf 1_A)(x)=Q^h_{t+s}(x,A). ∫ Q t h ( x , d y ) Q s h ( y , A ) = ∫ e − λ t h ( x ) h ( y ) Q t ( x , d y ) e − λ s h ( y ) Q s ( h 1 A ) ( y ) = h ( x ) e − λ ( t + s ) Q t + s ( h 1 A ) ( x ) = Q t + s h ( x , A ) . 生成子の式は Q t h f − f = e − λ t h − 1 ( Q t ( h f ) − h f ) + ( e − λ t − 1 ) f Q^h_tf-f=e^{-\lambda t}h^{-1}\bigl(Q_t(hf)-hf\bigr)+(e^{-\lambda t}-1)f Q t h f − f = e − λ t h − 1 ( Q t ( h f ) − h f ) + ( e − λ t − 1 ) f を t t t で割れば得られます。
∎
例 5.3 (ドリフトをもつ Brown 運動 )
Q t = P t Q_t=P_t Q t = P t (Brown 運動)、θ ∈ R d \theta\in\mathbb{R}^d θ ∈ R d 、h ( x ) = e ⟨ θ , x ⟩ h(x)=e^{\langle\theta,x\rangle} h ( x ) = e ⟨ θ , x ⟩ とします。正規分布の積率母関数から P t h = e ∣ θ ∣ 2 t / 2 h P_th=e^{|\theta|^2t/2}h P t h = e ∣ θ ∣ 2 t /2 h です。平方完成すると
P t h f ( x ) = e − ∣ θ ∣ 2 t / 2 ∫ f ( y ) e ⟨ θ , y − x ⟩ p t ( y − x ) d y = ∫ f ( y ) p t ( y − x − θ t ) d y = E 0 f ( x + θ t + X t ) , P^h_tf(x)=e^{-|\theta|^2t/2}\int f(y)e^{\langle\theta,y-x\rangle}p_t(y-x)\,dy=\int f(y)p_t(y-x-\theta t)\,dy=E_0f(x+\theta t+X_t), P t h f ( x ) = e − ∣ θ ∣ 2 t /2 ∫ f ( y ) e ⟨ θ , y − x ⟩ p t ( y − x ) d y = ∫ f ( y ) p t ( y − x − θ t ) d y = E 0 f ( x + θ t + X t ) , つまり h 変換はドリフト θ \theta θ の Brown 運動 x + θ t + B t x+\theta t+B_t x + θ t + B t で、生成子は 1 2 Δ + θ ⋅ ∇ \frac12\Delta+\theta\cdot\nabla 2 1 Δ + θ ⋅ ∇ です。道の空間では、F t \mathcal{F}_t F t 上で P x P_x P x に密度 e ⟨ θ , X t − x ⟩ − ∣ θ ∣ 2 t / 2 e^{\langle\theta,X_t-x\rangle-|\theta|^2t/2} e ⟨ θ , X t − x ⟩ − ∣ θ ∣ 2 t /2 を掛けた測度がこの過程の法則で、密度は指数マルチンゲール(定理 6.1[マルチンゲールとブラウン運動] (c))です。
Schrödinger 作用素 H = − 1 2 Δ + V H=-\frac12\Delta+V H = − 2 1 Δ + V の基底状態 ψ > 0 \psi>0 ψ > 0 、H ψ = E ψ H\psi=E\psi H ψ = E ψ は、Feynman–Kac 半群の固有関数 T t V ψ = e − E t ψ T^V_t\psi=e^{-Et}\psi T t V ψ = e − E t ψ になるはずです。これを確かめてから h 変換を適用します。
定理 5.4 (基底状態変換 ) ✓
V V V を R d \mathbb{R}^d R d 上の有界かつ一様連続な関数とし、ψ ∈ C 2 ( R d ) \psi\in C^2(\mathbb{R}^d) ψ ∈ C 2 ( R d ) は ψ > 0 \psi>0 ψ > 0 で、ψ \psi ψ と 2 階までの偏導関数が有界であり、ある E ∈ R E\in\mathbb{R} E ∈ R について
− 1 2 Δ ψ + V ψ = E ψ -\tfrac12\Delta\psi+V\psi=E\psi − 2 1 Δ ψ + V ψ = E ψ を満たすとする。このとき
すべての t ≥ 0 t\ge0 t ≥ 0 で T t V ψ = e − E t ψ T^V_t\psi=e^{-Et}\psi T t V ψ = e − E t ψ 。
P t ψ f = e E t ψ − 1 T t V ( ψ f ) P^\psi_tf=e^{Et}\psi^{-1}T^V_t(\psi f) P t ψ f = e E t ψ − 1 T t V ( ψ f ) は Markov 半群である。
f ∈ C c 2 ( R d ) f\in C^2_c(\mathbb{R}^d) f ∈ C c 2 ( R d ) なら、各点 x x x で lim t ↓ 0 P t ψ f ( x ) − f ( x ) t = 1 2 Δ f ( x ) + ∇ log ψ ( x ) ⋅ ∇ f ( x ) \displaystyle\lim_{t\downarrow0}\frac{P^\psi_tf(x)-f(x)}{t}=\frac12\Delta f(x)+\nabla\log\psi(x)\cdot\nabla f(x) t ↓ 0 lim t P t ψ f ( x ) − f ( x ) = 2 1 Δ f ( x ) + ∇ log ψ ( x ) ⋅ ∇ f ( x ) 。
さらに ψ ∈ L 2 ( R d ) \psi\in L^2(\mathbb{R}^d) ψ ∈ L 2 ( R d ) なら、測度 ψ ( x ) 2 d x \psi(x)^2dx ψ ( x ) 2 d x は ( P t ψ ) (P^\psi_t) ( P t ψ ) の不変測度である:有界可測な f f f について ∫ P t ψ f ψ 2 d x = ∫ f ψ 2 d x \int P^\psi_tf\,\psi^2dx=\int f\psi^2dx ∫ P t ψ f ψ 2 d x = ∫ f ψ 2 d x 。
証明(定理 5.4) (1) x x x を固定し g ( s ) = T s V ψ ( x ) g(s)=T^V_s\psi(x) g ( s ) = T s V ψ ( x ) とします。道の連続性と優収束定理(被積分関数は e s ∥ V ∥ ∥ ψ ∥ e^{s\|V\|}\|\psi\| e s ∥ V ∥ ∥ ψ ∥ 以下)から g g g は連続です。定理 4.6 (2) の半群性(有界な V V V は Kato クラス)から g ( s + h ) − g ( s ) = T s V ( T h V ψ − ψ ) ( x ) g(s+h)-g(s)=T^V_s(T^V_h\psi-\psi)(x) g ( s + h ) − g ( s ) = T s V ( T h V ψ − ψ ) ( x ) です。各点 y y y で
T h V ψ ( y ) − ψ ( y ) h = P h ψ ( y ) − ψ ( y ) h + E y [ e − ∫ 0 h V ( X s ) d s − 1 h ψ ( X h ) ] → 1 2 Δ ψ ( y ) − V ( y ) ψ ( y ) = − E ψ ( y ) . \frac{T^V_h\psi(y)-\psi(y)}{h}=\frac{P_h\psi(y)-\psi(y)}{h}+E_y\Bigl[\frac{e^{-\int_0^hV(X_s)ds}-1}{h}\psi(X_h)\Bigr]\to\tfrac12\Delta\psi(y)-V(y)\psi(y)=-E\psi(y). h T h V ψ ( y ) − ψ ( y ) = h P h ψ ( y ) − ψ ( y ) + E y [ h e − ∫ 0 h V ( X s ) d s − 1 ψ ( X h ) ] → 2 1 Δ ψ ( y ) − V ( y ) ψ ( y ) = − E ψ ( y ) . 第 1 項の収束は 定理 3.5 の生成子の計算を点 y y y で行えば(Hesse 行列の連続性だけで)得られ、Taylor の定理から ∣ P h ψ − ψ ∣ / h ≤ d 2 2 ∥ ∇ 2 ψ ∥ max |P_h\psi-\psi|/h\le\frac{d^2}2\|\nabla^2\psi\|_{\max} ∣ P h ψ − ψ ∣/ h ≤ 2 d 2 ∥ ∇ 2 ψ ∥ m a x と一様に有界です。第 2 項の収束は V V V 、ψ \psi ψ の連続性と優収束定理により、有界性は ∥ V ∥ e h ∥ V ∥ ∥ ψ ∥ \|V\|e^{h\|V\|}\|\psi\| ∥ V ∥ e h ∥ V ∥ ∥ ψ ∥ です。T s V T^V_s T s V は有界な重み e − ∫ 0 s V e^{-\int_0^sV} e − ∫ 0 s V をもつ期待値なので、優収束定理から右微分 D + g ( s ) = T s V ( − E ψ ) ( x ) = − E g ( s ) D^+g(s)=T^V_s(-E\psi)(x)=-Eg(s) D + g ( s ) = T s V ( − E ψ ) ( x ) = − E g ( s ) 。F ( s ) = e E s g ( s ) F(s)=e^{Es}g(s) F ( s ) = e E s g ( s ) は連続で D + F = 0 D^+F=0 D + F = 0 なので、補題 3.3 により F ( s ) = F ( 0 ) = ψ ( x ) F(s)=F(0)=\psi(x) F ( s ) = F ( 0 ) = ψ ( x ) です。
(2) (1) と 定理 5.2 (Q t = T t V Q_t=T^V_t Q t = T t V 、h = ψ h=\psi h = ψ 、λ = − E \lambda=-E λ = − E )から従います。
(3) ψ f ∈ C c 2 ⊂ D ( A ) \psi f\in C^2_c\subset D(A) ψ f ∈ C c 2 ⊂ D ( A ) なので、定理 4.2 から ( T t V ( ψ f ) − ψ f ) / t → 1 2 Δ ( ψ f ) − V ψ f (T^V_t(\psi f)-\psi f)/t\to\frac12\Delta(\psi f)-V\psi f ( T t V ( ψ f ) − ψ f ) / t → 2 1 Δ ( ψ f ) − V ψ f (一様に)。定理 5.2 の生成子の式に代入し、1 2 Δ ( ψ f ) = 1 2 f Δ ψ + ∇ ψ ⋅ ∇ f + 1 2 ψ Δ f \frac12\Delta(\psi f)=\frac12f\Delta\psi+\nabla\psi\cdot\nabla f+\frac12\psi\Delta f 2 1 Δ ( ψ f ) = 2 1 f Δ ψ + ∇ ψ ⋅ ∇ f + 2 1 ψ Δ f と 1 2 Δ ψ − V ψ + E ψ = 0 \frac12\Delta\psi-V\psi+E\psi=0 2 1 Δ ψ − V ψ + E ψ = 0 を使うと
1 2 Δ ( ψ f ) − V ψ f ψ + E f = 1 2 Δ f + ∇ ψ ψ ⋅ ∇ f + f ψ ( 1 2 Δ ψ − V ψ + E ψ ) = 1 2 Δ f + ∇ log ψ ⋅ ∇ f . \frac{\frac12\Delta(\psi f)-V\psi f}{\psi}+Ef=\frac12\Delta f+\frac{\nabla\psi}{\psi}\cdot\nabla f+\frac{f}{\psi}\Bigl(\frac12\Delta\psi-V\psi+E\psi\Bigr)=\frac12\Delta f+\nabla\log\psi\cdot\nabla f . ψ 2 1 Δ ( ψ f ) − V ψ f + E f = 2 1 Δ f + ψ ∇ ψ ⋅ ∇ f + ψ f ( 2 1 Δ ψ − V ψ + E ψ ) = 2 1 Δ f + ∇ log ψ ⋅ ∇ f . (4) 定理 2.1[Feynman–Kac 表示と場の模型の汎関数積分] (M = 1 M=1 M = 1 )により、L 2 L^2 L 2 の元 ϕ \phi ϕ について e − t H ϕ = T t V ϕ e^{-tH}\phi=T^V_t\phi e − t H ϕ = T t V ϕ (a.e.)で、e − t H e^{-tH} e − t H は自己共役です。ψ \psi ψ と ψ f \psi f ψ f は L 2 L^2 L 2 に属するので、(1) を使って
∫ P t ψ f ψ 2 d x = e E t ⟨ ψ , e − t H ( ψ f ) ⟩ = e E t ⟨ e − t H ψ , ψ f ⟩ = e E t ⟨ e − E t ψ , ψ f ⟩ = ∫ f ψ 2 d x . \int P^\psi_tf\,\psi^2dx=e^{Et}\langle\psi,e^{-tH}(\psi f)\rangle=e^{Et}\langle e^{-tH}\psi,\psi f\rangle=e^{Et}\langle e^{-Et}\psi,\psi f\rangle=\int f\psi^2dx . ∫ P t ψ f ψ 2 d x = e E t ⟨ ψ , e − t H ( ψ f )⟩ = e E t ⟨ e − t H ψ , ψ f ⟩ = e E t ⟨ e − E t ψ , ψ f ⟩ = ∫ f ψ 2 d x . ∎
基底状態変換で、量子力学の基底状態 ψ \psi ψ は「ドリフト ∇ log ψ \nabla\log\psi ∇ log ψ で原点付近に引き戻される拡散」に、確率密度 ψ 2 \psi^2 ψ 2 はその拡散の平衡分布 (時間発展で分布が変わらない測度。定理 5.4 (4))に対応します。Feynman–Kac 半群の「確率が減る(V > 0 V>0 V > 0 )・増える(V < 0 V<0 V < 0 )」という重みが、ドリフトに吸収されたわけです。
例 5.5 (Pöschl–Teller ポテンシャル )
d = 1 d=1 d = 1 、V ( x ) = − sech 2 x V(x)=-\operatorname{sech}^2x V ( x ) = − sech 2 x 、ψ ( x ) = sech x \psi(x)=\operatorname{sech}x ψ ( x ) = sech x とします。ψ ′ = − sech x tanh x \psi'=-\operatorname{sech}x\tanh x ψ ′ = − sech x tanh x 、ψ ′ ′ = sech x ( 1 − 2 sech 2 x ) \psi''=\operatorname{sech}x\,(1-2\operatorname{sech}^2x) ψ ′′ = sech x ( 1 − 2 sech 2 x ) なので
− 1 2 ψ ′ ′ + V ψ = − 1 2 sech x + sech 3 x − sech 3 x = − 1 2 ψ , -\tfrac12\psi''+V\psi=-\tfrac12\operatorname{sech}x+\operatorname{sech}^3x-\operatorname{sech}^3x=-\tfrac12\psi , − 2 1 ψ ′′ + V ψ = − 2 1 sech x + sech 3 x − sech 3 x = − 2 1 ψ , すなわち E = − 1 2 E=-\frac12 E = − 2 1 です。V V V は導関数が有界なので一様連続、ψ \psi ψ と導関数は有界で ψ ∈ L 2 \psi\in L^2 ψ ∈ L 2 なので、定理 5.4 がすべて使えます。変換後の生成子は 1 2 f ′ ′ − tanh ( x ) f ′ \frac12f''-\tanh(x)f' 2 1 f ′′ − tanh ( x ) f ′ で、原点から離れるほど強く(最大で速さ 1 1 1 で)引き戻すドリフトをもち、平衡分布は 1 2 sech 2 x d x \frac12\operatorname{sech}^2x\,dx 2 1 sech 2 x d x です(∫ sech 2 x d x = 2 \int\operatorname{sech}^2x\,dx=2 ∫ sech 2 x d x = 2 )。
Pöschl–Teller ポテンシャルの基底状態変換:平衡密度 sech²(x)/2(曲線)とドリフト −tanh x(矢印)。ドリフトは原点に向かい、|x| が大きいところで速さ 1 に近づく
調和振動子 V ( x ) = 1 2 x 2 V(x)=\frac12x^2 V ( x ) = 2 1 x 2 、ψ = e − x 2 / 2 \psi=e^{-x^2/2} ψ = e − x 2 /2 に同じ計算をすると、生成子 1 2 f ′ ′ − x f ′ \frac12f''-xf' 2 1 f ′′ − x f ′ の Ornstein–Uhlenbeck 過程と平衡分布 N ( 0 , 1 2 ) N(0,\frac12) N ( 0 , 2 1 ) が得られます。ただし V V V が非有界なので 定理 5.4 の仮定を満たさず、この場合の正当化には別の議論(Mehler 核の直接計算など)が要ります。
以下 d = 1 d=1 d = 1 とし、( B t ) (B_t) ( B t ) を標準 Brown 運動、( F t ) (\mathcal{F}_t) ( F t ) をその自然なフィルトレーションとします。[ 0 , T ] [0,T] [ 0 , T ] の分割 Δ : 0 = t 0 < t 1 < ⋯ < t n = T \Delta:0=t_0<t_1<\cdots<t_n=T Δ : 0 = t 0 < t 1 < ⋯ < t n = T の幅を ∣ Δ ∣ = max i ( t i + 1 − t i ) |\Delta|=\max_i(t_{i+1}-t_i) ∣Δ∣ = max i ( t i + 1 − t i ) 、Δ B i = B t i + 1 − B t i \Delta B_i=B_{t_{i+1}}-B_{t_i} Δ B i = B t i + 1 − B t i 、h i = t i + 1 − t i h_i=t_{i+1}-t_i h i = t i + 1 − t i と書きます。
定義 6.1 (Stratonovich 積分 ) ↗
連続関数 f f f について、台形則の Riemann 和
S Δ ( f ) = ∑ i = 0 n − 1 f ( B t i ) + f ( B t i + 1 ) 2 Δ B i S_\Delta(f)=\sum_{i=0}^{n-1}\frac{f(B_{t_i})+f(B_{t_{i+1}})}{2}\,\Delta B_i S Δ ( f ) = i = 0 ∑ n − 1 2 f ( B t i ) + f ( B t i + 1 ) Δ B i が ∣ Δ ∣ → 0 |\Delta|\to0 ∣Δ∣ → 0 のとき確率収束するならば、その極限を ∫ 0 T f ( B s ) ∘ d B s \int_0^Tf(B_s)\circ dB_s ∫ 0 T f ( B s ) ∘ d B s と書き、Stratonovich 積分 という。
左端点 f ( B t i ) f(B_{t_i}) f ( B t i ) を使う和は伊藤積分に収束し(下の証明の第 1 段)、右端点との平均を使う和は別の極限をもちます。Riemann–Stieltjes 積分なら端点の選び方によらないはずのところで差が出るのは、Brown 運動の道が有界変動でないからです(系 3.2[確率微分方程式] )。
定理 6.2 (Stratonovich 積分と伊藤積分の関係 ) ✓
f ∈ C 2 ( R ) f\in C^2(\mathbb{R}) f ∈ C 2 ( R ) で f ′ f' f ′ 、f ′ ′ f'' f ′′ が有界ならば、∫ 0 T f ( B s ) ∘ d B s \int_0^Tf(B_s)\circ dB_s ∫ 0 T f ( B s ) ∘ d B s が存在して
∫ 0 T f ( B s ) ∘ d B s = ∫ 0 T f ( B s ) d B s + 1 2 ∫ 0 T f ′ ( B s ) d s . \int_0^Tf(B_s)\circ dB_s=\int_0^Tf(B_s)\,dB_s+\frac12\int_0^Tf'(B_s)\,ds . ∫ 0 T f ( B s ) ∘ d B s = ∫ 0 T f ( B s ) d B s + 2 1 ∫ 0 T f ′ ( B s ) d s . 収束は L 1 ( P ) L^1(P) L 1 ( P ) で起こる。
証明(定理 6.2) S Δ ( f ) = L Δ + R Δ S_\Delta(f)=L_\Delta+R_\Delta S Δ ( f ) = L Δ + R Δ 、L Δ = ∑ i f ( B t i ) Δ B i L_\Delta=\sum_if(B_{t_i})\Delta B_i L Δ = ∑ i f ( B t i ) Δ B i 、R Δ = 1 2 ∑ i ( f ( B t i + 1 ) − f ( B t i ) ) Δ B i R_\Delta=\frac12\sum_i\bigl(f(B_{t_{i+1}})-f(B_{t_i})\bigr)\Delta B_i R Δ = 2 1 ∑ i ( f ( B t i + 1 ) − f ( B t i ) ) Δ B i と分けます。
(1) L Δ → ∫ 0 T f ( B ) d B L_\Delta\to\int_0^Tf(B)\,dB L Δ → ∫ 0 T f ( B ) d B (L 2 L^2 L 2 )。 ∣ f ( x ) ∣ ≤ ∣ f ( 0 ) ∣ + ∥ f ′ ∥ ∣ x ∣ |f(x)|\le|f(0)|+\|f'\||x| ∣ f ( x ) ∣ ≤ ∣ f ( 0 ) ∣ + ∥ f ′ ∥∣ x ∣ なので f ( B ) f(B) f ( B ) は L 2 [ 0 , T ] \mathcal{L}^2[0,T] L 2 [ 0 , T ] に属し、X t Δ = ∑ i f ( B t i ) 1 ( t i , t i + 1 ] ( t ) X^\Delta_t=\sum_if(B_{t_i})\mathbf 1_{(t_i,t_{i+1}]}(t) X t Δ = ∑ i f ( B t i ) 1 ( t i , t i + 1 ] ( t ) は単純過程で L Δ = ∫ X Δ d B L_\Delta=\int X^\Delta dB L Δ = ∫ X Δ d B です(定義 4.2[確率微分方程式] )。E ∣ f ( B t ) − f ( B t i ) ∣ 2 ≤ ∥ f ′ ∥ 2 ( t − t i ) E|f(B_t)-f(B_{t_i})|^2\le\|f'\|^2(t-t_i) E ∣ f ( B t ) − f ( B t i ) ∣ 2 ≤ ∥ f ′ ∥ 2 ( t − t i ) から E ∫ 0 T ∣ f ( B t ) − X t Δ ∣ 2 d t ≤ ∥ f ′ ∥ 2 T ∣ Δ ∣ → 0 E\int_0^T|f(B_t)-X^\Delta_t|^2dt\le\|f'\|^2T|\Delta|\to0 E ∫ 0 T ∣ f ( B t ) − X t Δ ∣ 2 d t ≤ ∥ f ′ ∥ 2 T ∣Δ∣ → 0 。伊藤等長性(定理 4.7[確率微分方程式] )により E ∣ L Δ − ∫ f ( B ) d B ∣ 2 → 0 E|L_\Delta-\int f(B)dB|^2\to0 E ∣ L Δ − ∫ f ( B ) d B ∣ 2 → 0 。
(2) R Δ → 1 2 ∫ 0 T f ′ ( B s ) d s R_\Delta\to\frac12\int_0^Tf'(B_s)ds R Δ → 2 1 ∫ 0 T f ′ ( B s ) d s (L 1 L^1 L 1 )。 Taylor の定理から f ( B t i + 1 ) − f ( B t i ) = f ′ ( B t i ) Δ B i + r i f(B_{t_{i+1}})-f(B_{t_i})=f'(B_{t_i})\Delta B_i+r_i f ( B t i + 1 ) − f ( B t i ) = f ′ ( B t i ) Δ B i + r i 、∣ r i ∣ ≤ 1 2 ∥ f ′ ′ ∥ ( Δ B i ) 2 |r_i|\le\frac12\|f''\|(\Delta B_i)^2 ∣ r i ∣ ≤ 2 1 ∥ f ′′ ∥ ( Δ B i ) 2 。よって
2 R Δ = ∑ i f ′ ( B t i ) h i + ∑ i f ′ ( B t i ) ( ( Δ B i ) 2 − h i ) + ∑ i r i Δ B i . 2R_\Delta=\sum_if'(B_{t_i})h_i+\sum_if'(B_{t_i})\bigl((\Delta B_i)^2-h_i\bigr)+\sum_ir_i\Delta B_i . 2 R Δ = i ∑ f ′ ( B t i ) h i + i ∑ f ′ ( B t i ) ( ( Δ B i ) 2 − h i ) + i ∑ r i Δ B i . 第 1 項は連続関数 s ↦ f ′ ( B s ) s\mapsto f'(B_s) s ↦ f ′ ( B s ) の Riemann 和なので各道で ∫ 0 T f ′ ( B s ) d s \int_0^Tf'(B_s)ds ∫ 0 T f ′ ( B s ) d s に収束し、∥ f ′ ∥ T \|f'\|T ∥ f ′ ∥ T で押さえられるので L 1 L^1 L 1 でも収束します。第 2 項の各項 ξ i \xi_i ξ i について、i < j i<j i < j なら ξ i f ′ ( B t j ) \xi_if'(B_{t_j}) ξ i f ′ ( B t j ) は F t j \mathcal{F}_{t_j} F t j -可測で ( Δ B j ) 2 − h j (\Delta B_j)^2-h_j ( Δ B j ) 2 − h j は F t j \mathcal{F}_{t_j} F t j と独立で平均 0 0 0 なので、命題 4.5[条件付き期待値] から E [ ξ i ξ j ] = 0 E[\xi_i\xi_j]=0 E [ ξ i ξ j ] = 0 。E [ ξ i 2 ] ≤ ∥ f ′ ∥ 2 V a r ( ( Δ B i ) 2 ) = 2 ∥ f ′ ∥ 2 h i 2 E[\xi_i^2]\le\|f'\|^2\mathrm{Var}((\Delta B_i)^2)=2\|f'\|^2h_i^2 E [ ξ i 2 ] ≤ ∥ f ′ ∥ 2 Var (( Δ B i ) 2 ) = 2∥ f ′ ∥ 2 h i 2 なので、E ∣ ∑ ξ i ∣ 2 ≤ 2 ∥ f ′ ∥ 2 T ∣ Δ ∣ → 0 E|\sum\xi_i|^2\le2\|f'\|^2T|\Delta|\to0 E ∣ ∑ ξ i ∣ 2 ≤ 2∥ f ′ ∥ 2 T ∣Δ∣ → 0 。第 3 項は E ∑ ∣ r i Δ B i ∣ ≤ 1 2 ∥ f ′ ′ ∥ ∑ i E ∣ Δ B i ∣ 3 = 1 2 ∥ f ′ ′ ∥ E ∣ Z ∣ 3 ∑ i h i 3 / 2 ≤ 1 2 ∥ f ′ ′ ∥ E ∣ Z ∣ 3 T ∣ Δ ∣ 1 / 2 → 0 E\sum|r_i\Delta B_i|\le\frac12\|f''\|\sum_iE|\Delta B_i|^3=\frac12\|f''\|\,E|Z|^3\sum_ih_i^{3/2}\le\frac12\|f''\|E|Z|^3\,T|\Delta|^{1/2}\to0 E ∑ ∣ r i Δ B i ∣ ≤ 2 1 ∥ f ′′ ∥ ∑ i E ∣Δ B i ∣ 3 = 2 1 ∥ f ′′ ∥ E ∣ Z ∣ 3 ∑ i h i 3/2 ≤ 2 1 ∥ f ′′ ∥ E ∣ Z ∣ 3 T ∣Δ ∣ 1/2 → 0 (Z ∼ N ( 0 , 1 ) Z\sim N(0,1) Z ∼ N ( 0 , 1 ) )。
∎
系 6.3 (Stratonovich 積分の連鎖律 ) ✓
F ∈ C 3 ( R ) F\in C^3(\mathbb{R}) F ∈ C 3 ( R ) で F ′ ′ F'' F ′′ 、F ′ ′ ′ F''' F ′′′ が有界ならば
F ( B T ) − F ( B 0 ) = ∫ 0 T F ′ ( B s ) ∘ d B s . F(B_T)-F(B_0)=\int_0^TF'(B_s)\circ dB_s . F ( B T ) − F ( B 0 ) = ∫ 0 T F ′ ( B s ) ∘ d B s .
証明(系 6.3) 伊藤の公式(定理 5.2[確率微分方程式] 、f ( t , x ) = F ( x ) f(t,x)=F(x) f ( t , x ) = F ( x ) 、X = B X=B X = B )から F ( B T ) − F ( B 0 ) = ∫ 0 T F ′ ( B ) d B + 1 2 ∫ 0 T F ′ ′ ( B ) d s F(B_T)-F(B_0)=\int_0^TF'(B)dB+\frac12\int_0^TF''(B)ds F ( B T ) − F ( B 0 ) = ∫ 0 T F ′ ( B ) d B + 2 1 ∫ 0 T F ′′ ( B ) d s 。定理 6.2 を f = F ′ f=F' f = F ′ (f ′ = F ′ ′ f'=F'' f ′ = F ′′ 、f ′ ′ = F ′ ′ ′ f''=F''' f ′′ = F ′′′ は有界)に使えば、右辺は ∫ 0 T F ′ ( B ) ∘ d B \int_0^TF'(B)\circ dB ∫ 0 T F ′ ( B ) ∘ d B です。
∎
伊藤積分はマルチンゲールになる(未来を見ない)代わりに連鎖律に補正項 1 2 F ′ ′ \frac12F'' 2 1 F ′′ が付き、Stratonovich 積分は連鎖律がそのまま成り立つ代わりにマルチンゲール性を失います。微分方程式で書けば、Stratonovich 形 d X = σ ( X ) ∘ d B + b ( X ) d t dX=\sigma(X)\circ dB+b(X)dt d X = σ ( X ) ∘ d B + b ( X ) d t は伊藤形 d X = σ ( X ) d B + ( b + 1 2 σ σ ′ ) ( X ) d t dX=\sigma(X)dB+\bigl(b+\frac12\sigma\sigma'\bigr)(X)dt d X = σ ( X ) d B + ( b + 2 1 σ σ ′ ) ( X ) d t と同じ解をもち、生成子は 1 2 ( σ ∂ x ) 2 + b ∂ x \frac12(\sigma\partial_x)^2+b\partial_x 2 1 ( σ ∂ x ) 2 + b ∂ x と「ベクトル場の 2 乗」で書けます(( σ ∂ x ) 2 f = σ 2 f ′ ′ + σ σ ′ f ′ (\sigma\partial_x)^2f=\sigma^2f''+\sigma\sigma'f' ( σ ∂ x ) 2 f = σ 2 f ′′ + σ σ ′ f ′ )。座標変換で形が変わらないので、多様体上の Brown 運動や、磁場中の粒子の Feynman–Kac–伊藤公式は Stratonovich 形で書くのが自然です。滑らかな近似の極限が Stratonovich 解になるという Wong–Zakai の定理(Brown 運動の道を滑らかな道で近似した方程式の解が、近似を細かくしたときに収束する先を述べる定理)も、この形の自然さを示しています(本稿では証明しません)。
演習 7.1 標準
d = 1 d=1 d = 1 とし、g ∈ C 0 ( R ) g\in C_0(\mathbb{R}) g ∈ C 0 ( R ) 、λ > 0 \lambda>0 λ > 0 、κ = 2 λ \kappa=\sqrt{2\lambda} κ = 2 λ とする。u ( x ) = 1 κ ∫ e − κ ∣ x − y ∣ g ( y ) d y u(x)=\frac1\kappa\int e^{-\kappa|x-y|}g(y)\,dy u ( x ) = κ 1 ∫ e − κ ∣ x − y ∣ g ( y ) d y が C 2 C^2 C 2 級で λ u − 1 2 u ′ ′ = g \lambda u-\frac12u''=g λ u − 2 1 u ′′ = g を満たすことを確かめよ。命題 3.2 (2) と合わせると、これは Brown 運動の半群の Laplace 変換 ∫ 0 ∞ e − λ t P t g d t \int_0^\infty e^{-\lambda t}P_tg\,dt ∫ 0 ∞ e − λ t P t g d t の明示式を与える(u ∈ C 0 2 u\in C^2_0 u ∈ C 0 2 のとき)。
解答 u ( x ) = 1 κ ( e − κ x ∫ − ∞ x e κ y g ( y ) d y + e κ x ∫ x ∞ e − κ y g ( y ) d y ) u(x)=\frac1\kappa\bigl(e^{-\kappa x}\int_{-\infty}^xe^{\kappa y}g(y)dy+e^{\kappa x}\int_x^\infty e^{-\kappa y}g(y)dy\bigr) u ( x ) = κ 1 ( e − κ x ∫ − ∞ x e κ y g ( y ) d y + e κ x ∫ x ∞ e − κ y g ( y ) d y ) と書くと、微分して u ′ ( x ) = − e − κ x ∫ − ∞ x e κ y g d y + e κ x ∫ x ∞ e − κ y g d y u'(x)=-e^{-\kappa x}\int_{-\infty}^xe^{\kappa y}g\,dy+e^{\kappa x}\int_x^\infty e^{-\kappa y}g\,dy u ′ ( x ) = − e − κ x ∫ − ∞ x e κ y g d y + e κ x ∫ x ∞ e − κ y g d y (境界項 1 κ g ( x ) − 1 κ g ( x ) \frac1\kappa g(x)-\frac1\kappa g(x) κ 1 g ( x ) − κ 1 g ( x ) は打ち消す)。もう一度微分すると u ′ ′ = κ 2 u − 2 g u''=\kappa^2u-2g u ′′ = κ 2 u − 2 g 。よって 1 2 u ′ ′ = λ u − g \frac12u''=\lambda u-g 2 1 u ′′ = λ u − g 。u ∈ C 0 2 u\in C^2_0 u ∈ C 0 2 のとき、( λ − A ) u = g (\lambda-A)u=g ( λ − A ) u = g と ( λ − A ) (\lambda-A) ( λ − A ) の単射性(命題 3.2 (4))から u = ( λ − A ) − 1 g = ∫ 0 ∞ e − λ t P t g d t u=(\lambda-A)^{-1}g=\int_0^\infty e^{-\lambda t}P_tg\,dt u = ( λ − A ) − 1 g = ∫ 0 ∞ e − λ t P t g d t です。
演習 7.3 易
∫ 0 T B s ∘ d B s = 1 2 B T 2 \int_0^TB_s\circ dB_s=\frac12B_T^2 ∫ 0 T B s ∘ d B s = 2 1 B T 2 を定義から直接示し、伊藤積分 ∫ 0 T B s d B s \int_0^TB_s\,dB_s ∫ 0 T B s d B s を求めよ。
解答 f ( x ) = x f(x)=x f ( x ) = x では S Δ = ∑ i 1 2 ( B t i + B t i + 1 ) ( B t i + 1 − B t i ) = 1 2 ∑ i ( B t i + 1 2 − B t i 2 ) = 1 2 B T 2 S_\Delta=\sum_i\frac12(B_{t_i}+B_{t_{i+1}})(B_{t_{i+1}}-B_{t_i})=\frac12\sum_i(B_{t_{i+1}}^2-B_{t_i}^2)=\frac12B_T^2 S Δ = ∑ i 2 1 ( B t i + B t i + 1 ) ( B t i + 1 − B t i ) = 2 1 ∑ i ( B t i + 1 2 − B t i 2 ) = 2 1 B T 2 と、分割によらず厳密に一致します。定理 6.2 から ∫ 0 T B d B = 1 2 B T 2 − 1 2 ∫ 0 T 1 d s = 1 2 ( B T 2 − T ) \int_0^TB\,dB=\frac12B_T^2-\frac12\int_0^T1\,ds=\frac12(B_T^2-T) ∫ 0 T B d B = 2 1 B T 2 − 2 1 ∫ 0 T 1 d s = 2 1 ( B T 2 − T ) 。
演習 7.4 標準
d = 1 d=1 d = 1 、W = c 1 [ − 1 , 1 ] W=c\mathbf 1_{[-1,1]} W = c 1 [ − 1 , 1 ] (c > 0 c>0 c > 0 )とする。sup x E x ∫ 0 t W ( X s ) d s ≤ c t \sup_xE_x\int_0^tW(X_s)ds\le ct sup x E x ∫ 0 t W ( X s ) d s ≤ c t を示し、c t < 1 ct<1 c t < 1 のとき 補題 4.5 の評価と、自明な評価 E x e ∫ 0 t W ≤ e c t E_xe^{\int_0^tW}\le e^{ct} E x e ∫ 0 t W ≤ e c t を比べよ。この W W W では Khasminskii の補題が不要な理由と、Coulomb ポテンシャル(例 4.4 )では必要な理由を説明せよ。
解答 0 ≤ W ≤ c 0\le W\le c 0 ≤ W ≤ c なので ∫ 0 t W ( X s ) d s ≤ c t \int_0^tW(X_s)ds\le ct ∫ 0 t W ( X s ) d s ≤ c t 、したがって E x e ∫ W ≤ e c t E_xe^{\int W}\le e^{ct} E x e ∫ W ≤ e c t が各道の評価から出ます。Khasminskii の評価 1 / ( 1 − c t ) 1/(1-ct) 1/ ( 1 − c t ) は e c t e^{ct} e c t 以上(e − c t ≥ 1 − c t e^{-ct}\ge1-ct e − c t ≥ 1 − c t )なので、有界な W W W では自明な評価のほうが良いです。Coulomb の場合は ∫ 0 t ∣ X s ∣ − 1 d s \int_0^t|X_s|^{-1}ds ∫ 0 t ∣ X s ∣ − 1 d s が道ごとには有界でなく(原点の近くを通る道ほど大きい)、各道の評価が使えません。期待値 E x ∫ 0 t ∣ X s ∣ − 1 d s E_x\int_0^t|X_s|^{-1}ds E x ∫ 0 t ∣ X s ∣ − 1 d s は一様に小さい(例 4.4 )ので、Khasminskii の補題がその平均の情報を指数モーメントに持ち上げます。
Stewart N. Ethier, Thomas G. Kurtz『Markov Processes: Characterization and Convergence』Wiley、1986 — 第 4 章(Feller 半群、正の最大値原理、生成子による特徴づけ)。
Daniel Revuz, Marc Yor『Continuous Martingales and Brownian Motion』(第 3 版)Springer、1999 — 第 III 章(Markov 過程と Feller 過程)。
Barry Simon, “Schrödinger semigroups”, Bulletin of the American Mathematical Society (N.S.) 7 (1982), 447–526.
M. Aizenman, B. Simon, “Brownian motion and Harnack inequality for Schrödinger operators”, Communications on Pure and Applied Mathematics 35 (1982), 209–273.
J. L. Doob, “Conditional Brownian motion and the boundary limits of harmonic functions”, Bulletin de la Société Mathématique de France 85 (1957), 431–458.
Elliott H. Lieb, Michael Loss『Analysis』(第 2 版)American Mathematical Society、2001 — 第 9 章(Newton の定理)。
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 2 回(Claude Code)・最終 2026-10-12