大問2題で、2問すべてが必答です(試験時間は60分)。第1問は3次元回転の集合 S O ( 3 ) SO(3) S O ( 3 ) を題材に、群としての性質の真偽判定から始めて、回転の固有値、無限小回転の生成子、交代行列の指数関数が回転になることまでを一本の線でたどります。第2問は f n ( x ) = e x 2 ( d / d x ) n e − x 2 f_n(x)=e^{x^2}(d/dx)^n e^{-x^2} f n ( x ) = e x 2 ( d / d x ) n e − x 2 で定義される多項式(符号を除いてエルミート多項式)の問題で、部分積分によるガウス型積分、コーシーの積分表示から母関数を導く部分、母関数が満たす偏微分方程式から各 f n f_n f n の常微分方程式を取り出す部分が山になります。どちらの問題も「示せ」「説明せよ」が多く、極限操作や場合分けの詰めがそのまま得点になります。
問題 分野 主題 第1問 線形代数・群論・微分方程式 3次元回転群、回転の固有値、交代行列と行列指数関数 第2問 特殊関数・複素解析・微分方程式 エルミート多項式、ガウス型積分、複素積分による母関数
R R R は3次元空間の回転を表す実 3 × 3 3\times 3 3 × 3 行列全体の集合で、単位行列 I I I も含みます。ここで回転とは、原点を固定する線形変換のうち、長さ(したがって内積)を保ち、かつ右手系を右手系に写す(行列式が正の)もの、という意味で使います。この解釈のもとで R = S O ( 3 ) R=SO(3) R = S O ( 3 ) です。
Ω ( n , θ ) \Omega(\boldsymbol{n},\theta) Ω ( n , θ ) は単位ベクトル n \boldsymbol{n} n を軸とする角 θ \theta θ の回転行列で、回転の正の向きは n \boldsymbol{n} n の方向に進む右ねじを回す向きとします。任意の実3次元ベクトル v 0 \boldsymbol{v}_0 v 0 に対して v ( θ ) = Ω ( n , θ ) v 0 \boldsymbol{v}(\theta)=\Omega(\boldsymbol{n},\theta)\boldsymbol{v}_0 v ( θ ) = Ω ( n , θ ) v 0 と置きます。また A ( a ) A(\boldsymbol{a}) A ( a ) は、任意の u \boldsymbol{u} u に対して A ( a ) u = a × u A(\boldsymbol{a})\boldsymbol{u}=\boldsymbol{a}\times\boldsymbol{u} A ( a ) u = a × u を満たす実 3 × 3 3\times 3 3 × 3 行列です。標準基底を e 1 , e 2 , e 3 \boldsymbol{e}_1,\boldsymbol{e}_2,\boldsymbol{e}_3 e 1 , e 2 , e 3 と書きます。
以下で繰り返し使う道具を1つ先に用意します。単位ベクトル n \boldsymbol{n} n に対して、( f 1 , f 2 , n ) (\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) ( f 1 , f 2 , n ) が右手系の正規直交系(つまり f 1 × f 2 = n \boldsymbol{f}_1\times\boldsymbol{f}_2=\boldsymbol{n} f 1 × f 2 = n 、n × f 1 = f 2 \boldsymbol{n}\times\boldsymbol{f}_1=\boldsymbol{f}_2 n × f 1 = f 2 、f 2 × n = f 1 \boldsymbol{f}_2\times\boldsymbol{n}=\boldsymbol{f}_1 f 2 × n = f 1 )となるように f 1 , f 2 \boldsymbol{f}_1,\boldsymbol{f}_2 f 1 , f 2 を取ると、Ω ( n , θ ) \Omega(\boldsymbol{n},\theta) Ω ( n , θ ) の作用は
Ω ( n , θ ) n = n , Ω ( n , θ ) f 1 = cos θ f 1 + sin θ f 2 , Ω ( n , θ ) f 2 = − sin θ f 1 + cos θ f 2 \Omega(\boldsymbol{n},\theta)\boldsymbol{n}=\boldsymbol{n},\qquad
\Omega(\boldsymbol{n},\theta)\boldsymbol{f}_1=\cos\theta\,\boldsymbol{f}_1+\sin\theta\,\boldsymbol{f}_2,\qquad
\Omega(\boldsymbol{n},\theta)\boldsymbol{f}_2=-\sin\theta\,\boldsymbol{f}_1+\cos\theta\,\boldsymbol{f}_2 Ω ( n , θ ) n = n , Ω ( n , θ ) f 1 = cos θ f 1 + sin θ f 2 , Ω ( n , θ ) f 2 = − sin θ f 1 + cos θ f 2
で尽くされます。これは「軸を動かさず、軸に垂直な平面内を右ねじの向きに角 θ \theta θ だけ回す」という定義をそのまま書き下したものです。
(a) ○ です。X , Y X,Y X , Y が長さを保つなら合成 X Y XY X Y も長さを保ちます。また det ( X Y ) = det X det Y > 0 \det(XY)=\det X\det Y>0 det ( X Y ) = det X det Y > 0 なので向きも保ちます。よって X Y XY X Y も回転で、R R R は行列の積について閉じています(実際 R R R は群をなします)。
(b) × です。軸の異なる回転は一般に交換しません。X = Ω ( e 1 , π / 2 ) X=\Omega(\boldsymbol{e}_1,\pi/2) X = Ω ( e 1 , π /2 ) 、Y = Ω ( e 3 , π / 2 ) Y=\Omega(\boldsymbol{e}_3,\pi/2) Y = Ω ( e 3 , π /2 ) 、すなわち
X = ( 1 0 0 0 0 − 1 0 1 0 ) , Y = ( 0 − 1 0 1 0 0 0 0 1 ) X=\begin{pmatrix}1&0&0\\0&0&-1\\0&1&0\end{pmatrix},\qquad
Y=\begin{pmatrix}0&-1&0\\1&0&0\\0&0&1\end{pmatrix} X = 1 0 0 0 0 1 0 − 1 0 , Y = 0 1 0 − 1 0 0 0 0 1
を取ると
X Y = ( 0 − 1 0 0 0 − 1 1 0 0 ) , Y X = ( 0 0 1 1 0 0 0 1 0 ) XY=\begin{pmatrix}0&-1&0\\0&0&-1\\1&0&0\end{pmatrix},\qquad
YX=\begin{pmatrix}0&0&1\\1&0&0\\0&1&0\end{pmatrix} X Y = 0 0 1 − 1 0 0 0 − 1 0 , Y X = 0 1 0 0 0 1 1 0 0
で、たとえば X Y e 1 = ( 0 , 0 , 1 ) XY\boldsymbol{e}_1=(0,0,1) X Y e 1 = ( 0 , 0 , 1 ) に対し Y X e 1 = ( 0 , 1 , 0 ) YX\boldsymbol{e}_1=(0,1,0) Y X e 1 = ( 0 , 1 , 0 ) となり一致しません。
(c) ○ です。X ∈ R X\in R X ∈ R は内積を保つので、任意の u , w \boldsymbol{u},\boldsymbol{w} u , w に対して
t u ( t X X ) w = t ( X u ) ( X w ) = t u w {}^{t}\boldsymbol{u}\,({}^{t}XX)\,\boldsymbol{w}={}^{t}(X\boldsymbol{u})(X\boldsymbol{w})={}^{t}\boldsymbol{u}\,\boldsymbol{w} t u ( t X X ) w = t ( X u ) ( X w ) = t u w
が成り立ちます。u = e i \boldsymbol{u}=\boldsymbol{e}_i u = e i 、w = e j \boldsymbol{w}=\boldsymbol{e}_j w = e j と取れば ( t X X ) i j = δ i j ({}^{t}XX)_{ij}=\delta_{ij} ( t X X ) ij = δ ij 、つまり t X X = I {}^{t}XX=I t X X = I です。これで X X X は正則で逆行列が t X {}^{t}X t X と分かるので、X t X = I X{}^{t}X=I X t X = I も従います。
(d) × です。直交性だけでは行列式の符号が決まりません。反例は鏡映
X = ( 1 0 0 0 1 0 0 0 − 1 ) X=\begin{pmatrix}1&0&0\\0&1&0\\0&0&-1\end{pmatrix} X = 1 0 0 0 1 0 0 0 − 1
で、t X X = X t X = I {}^{t}XX=X{}^{t}X=I t X X = X t X = I を満たしますが det X = − 1 \det X=-1 det X = − 1 であり、右手系を左手系に写すので回転ではありません(X = − I X=-I X = − I でも同じです)。直交行列全体 O ( 3 ) O(3) O ( 3 ) のうち det = + 1 \det=+1 det = + 1 の部分だけが R = S O ( 3 ) R=SO(3) R = S O ( 3 ) です。
(e) ○ です。(c) より X X X は実直交行列なので X † = t X = X − 1 X^{\dagger}={}^{t}X=X^{-1} X † = t X = X − 1 であり、X † X = X X † = I X^{\dagger}X=XX^{\dagger}=I X † X = X X † = I が成り立ちます。すなわち X X X はユニタリ行列であり、特に正規行列です。正規行列はユニタリ行列で対角化できる(スペクトル定理)ので、U † X U U^{\dagger}XU U † X U が対角行列となるユニタリ U U U が存在します。U U U の具体形は設問(ii)(2)で作ります。
まとめると (a)○、(b)×、(c)○、(d)×、(e)○ が答えです。
n = ( 0 , 0 , 1 ) \boldsymbol{n}=(0,0,1) n = ( 0 , 0 , 1 ) のときは設定で述べた作用そのままで、f 1 = e 1 \boldsymbol{f}_1=\boldsymbol{e}_1 f 1 = e 1 、f 2 = e 2 \boldsymbol{f}_2=\boldsymbol{e}_2 f 2 = e 2 と取れて
Ω ( n , θ ) = ( cos θ − sin θ 0 sin θ cos θ 0 0 0 1 ) \Omega(\boldsymbol{n},\theta)=\begin{pmatrix}\cos\theta&-\sin\theta&0\\ \sin\theta&\cos\theta&0\\ 0&0&1\end{pmatrix} Ω ( n , θ ) = cos θ sin θ 0 − sin θ cos θ 0 0 0 1
です。固有方程式は
det ( Ω − λ I ) = ( 1 − λ ) { ( cos θ − λ ) 2 + sin 2 θ } = ( 1 − λ ) ( λ 2 − 2 λ cos θ + 1 ) = 0 \det\bigl(\Omega-\lambda I\bigr)=(1-\lambda)\bigl\{(\cos\theta-\lambda)^2+\sin^2\theta\bigr\}=(1-\lambda)\bigl(\lambda^2-2\lambda\cos\theta+1\bigr)=0 det ( Ω − λ I ) = ( 1 − λ ) { ( cos θ − λ ) 2 + sin 2 θ } = ( 1 − λ ) ( λ 2 − 2 λ cos θ + 1 ) = 0
となり、λ 2 − 2 λ cos θ + 1 = 0 \lambda^2-2\lambda\cos\theta+1=0 λ 2 − 2 λ cos θ + 1 = 0 の根は λ = cos θ ± i sin θ \lambda=\cos\theta\pm i\sin\theta λ = cos θ ± i sin θ です。答えは
λ = 1 , e i θ , e − i θ \lambda=1,\quad e^{i\theta},\quad e^{-i\theta} λ = 1 , e i θ , e − i θ
の3つです。検算として Tr Ω = 1 + 2 cos θ \operatorname{Tr}\Omega=1+2\cos\theta Tr Ω = 1 + 2 cos θ が固有値の和 1 + e i θ + e − i θ 1+e^{i\theta}+e^{-i\theta} 1 + e i θ + e − i θ に、det Ω = 1 \det\Omega=1 det Ω = 1 が積 1 ⋅ e i θ ⋅ e − i θ 1\cdot e^{i\theta}\cdot e^{-i\theta} 1 ⋅ e i θ ⋅ e − i θ に一致します。
任意の単位ベクトル n \boldsymbol{n} n についても固有値は同じで、答えは λ = 1 , e i θ , e − i θ \lambda=1,\ e^{i\theta},\ e^{-i\theta} λ = 1 , e i θ , e − i θ です。
理由は基底の取り替えです。設定で用意した右手系正規直交系 ( f 1 , f 2 , n ) (\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) ( f 1 , f 2 , n ) を列に並べた行列を P = ( f 1 f 2 n ) P=(\boldsymbol{f}_1\ \boldsymbol{f}_2\ \boldsymbol{n}) P = ( f 1 f 2 n ) とすると P P P は直交行列で、Ω ( n , θ ) \Omega(\boldsymbol{n},\theta) Ω ( n , θ ) のこの基底での表現行列は設定で書いた作用より設問(ii)(1) と同じ形の行列 M ( θ ) M(\theta) M ( θ ) になります。すなわち
Ω ( n , θ ) = P M ( θ ) t P , M ( θ ) = ( cos θ − sin θ 0 sin θ cos θ 0 0 0 1 ) \Omega(\boldsymbol{n},\theta)=P\,M(\theta)\,{}^{t}P,\qquad
M(\theta)=\begin{pmatrix}\cos\theta&-\sin\theta&0\\ \sin\theta&\cos\theta&0\\ 0&0&1\end{pmatrix} Ω ( n , θ ) = P M ( θ ) t P , M ( θ ) = cos θ sin θ 0 − sin θ cos θ 0 0 0 1
です。相似変換で固有値は変わらないので、設問(ii)(1) の結果がそのまま使えます。
固有ベクトルも書き下せます。設定の作用から
Ω ( n , θ ) ( f 1 ∓ i f 2 ) = e ± i θ ( f 1 ∓ i f 2 ) , Ω ( n , θ ) n = n \Omega(\boldsymbol{n},\theta)(\boldsymbol{f}_1\mp i\boldsymbol{f}_2)=e^{\pm i\theta}(\boldsymbol{f}_1\mp i\boldsymbol{f}_2),\qquad \Omega(\boldsymbol{n},\theta)\boldsymbol{n}=\boldsymbol{n} Ω ( n , θ ) ( f 1 ∓ i f 2 ) = e ± i θ ( f 1 ∓ i f 2 ) , Ω ( n , θ ) n = n
が直接確かめられます。( f 1 − i f 2 ) / 2 (\boldsymbol{f}_1-i\boldsymbol{f}_2)/\sqrt{2} ( f 1 − i f 2 ) / 2 、( f 1 + i f 2 ) / 2 (\boldsymbol{f}_1+i\boldsymbol{f}_2)/\sqrt{2} ( f 1 + i f 2 ) / 2 、n \boldsymbol{n} n はエルミート内積について正規直交なので、これらを列に並べた U U U はユニタリで
U † Ω ( n , θ ) U = ( e i θ 0 0 0 e − i θ 0 0 0 1 ) U^{\dagger}\,\Omega(\boldsymbol{n},\theta)\,U=\begin{pmatrix}e^{i\theta}&0&0\\0&e^{-i\theta}&0\\0&0&1\end{pmatrix} U † Ω ( n , θ ) U = e i θ 0 0 0 e − i θ 0 0 0 1
となります。これが設問(i)(e) の U U U の具体形です。
同じ軸のまわりの回転は角度が足し算になります。実際、設定の基底 ( f 1 , f 2 , n ) (\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) ( f 1 , f 2 , n ) では Ω ( n , θ ) \Omega(\boldsymbol{n},\theta) Ω ( n , θ ) は M ( θ ) M(\theta) M ( θ ) で表され、三角関数の加法定理から M ( Δ θ ) M ( θ ) = M ( θ + Δ θ ) M(\Delta\theta)M(\theta)=M(\theta+\Delta\theta) M ( Δ θ ) M ( θ ) = M ( θ + Δ θ ) が成り立つので
v ( θ + Δ θ ) = Ω ( n , θ + Δ θ ) v 0 = Ω ( n , Δ θ ) v ( θ ) \boldsymbol{v}(\theta+\Delta\theta)=\Omega(\boldsymbol{n},\theta+\Delta\theta)\boldsymbol{v}_0=\Omega(\boldsymbol{n},\Delta\theta)\,\boldsymbol{v}(\theta) v ( θ + Δ θ ) = Ω ( n , θ + Δ θ ) v 0 = Ω ( n , Δ θ ) v ( θ )
です。つまり、θ \theta θ から θ + Δ θ \theta+\Delta\theta θ + Δ θ への変化は、いまの位置 v ( θ ) \boldsymbol{v}(\theta) v ( θ ) を軸 n \boldsymbol{n} n のまわりに Δ θ \Delta\theta Δ θ だけ回すことに等しくなります。
そこで u = v ( θ ) \boldsymbol{u}=\boldsymbol{v}(\theta) u = v ( θ ) を軸に平行な成分と垂直な成分に分けます。
u ∥ = ( n ⋅ u ) n , u ⊥ = u − u ∥ \boldsymbol{u}_{\parallel}=(\boldsymbol{n}\cdot\boldsymbol{u})\boldsymbol{n},\qquad \boldsymbol{u}_{\perp}=\boldsymbol{u}-\boldsymbol{u}_{\parallel} u ∥ = ( n ⋅ u ) n , u ⊥ = u − u ∥
u ⊥ ≠ 0 \boldsymbol{u}_{\perp}\neq\boldsymbol{0} u ⊥ = 0 の場合を先に扱います。ρ = ∣ u ⊥ ∣ > 0 \rho=|\boldsymbol{u}_{\perp}|>0 ρ = ∣ u ⊥ ∣ > 0 と置き、f 1 = u ⊥ / ρ \boldsymbol{f}_1=\boldsymbol{u}_{\perp}/\rho f 1 = u ⊥ / ρ 、f 2 = n × f 1 \boldsymbol{f}_2=\boldsymbol{n}\times\boldsymbol{f}_1 f 2 = n × f 1 と取れば ( f 1 , f 2 , n ) (\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) ( f 1 , f 2 , n ) は右手系の正規直交系です。回転は軸方向成分を動かさず、垂直成分を f 1 f 2 \boldsymbol{f}_1\boldsymbol{f}_2 f 1 f 2 平面内で Δ θ \Delta\theta Δ θ 回すので
Ω ( n , Δ θ ) u = u ∥ + ρ ( cos Δ θ f 1 + sin Δ θ f 2 ) \Omega(\boldsymbol{n},\Delta\theta)\boldsymbol{u}=\boldsymbol{u}_{\parallel}+\rho\bigl(\cos\Delta\theta\,\boldsymbol{f}_1+\sin\Delta\theta\,\boldsymbol{f}_2\bigr) Ω ( n , Δ θ ) u = u ∥ + ρ ( cos Δ θ f 1 + sin Δ θ f 2 )
となります。ここで ρ f 2 = n × u ⊥ = n × u \rho\boldsymbol{f}_2=\boldsymbol{n}\times\boldsymbol{u}_{\perp}=\boldsymbol{n}\times\boldsymbol{u} ρ f 2 = n × u ⊥ = n × u (n × u ∥ = 0 \boldsymbol{n}\times\boldsymbol{u}_{\parallel}=\boldsymbol{0} n × u ∥ = 0 を使いました)なので
v ( θ + Δ θ ) − v ( θ ) = ( cos Δ θ − 1 ) u ⊥ + sin Δ θ ( n × u ) \boldsymbol{v}(\theta+\Delta\theta)-\boldsymbol{v}(\theta)=(\cos\Delta\theta-1)\,\boldsymbol{u}_{\perp}+\sin\Delta\theta\,\bigl(\boldsymbol{n}\times\boldsymbol{u}\bigr) v ( θ + Δ θ ) − v ( θ ) = ( cos Δ θ − 1 ) u ⊥ + sin Δ θ ( n × u )
が厳密に成り立ちます。cos Δ θ − 1 = − Δ θ 2 / 2 + O ( Δ θ 4 ) \cos\Delta\theta-1=-\Delta\theta^2/2+O(\Delta\theta^4) cos Δ θ − 1 = − Δ θ 2 /2 + O ( Δ θ 4 ) 、sin Δ θ = Δ θ + O ( Δ θ 3 ) \sin\Delta\theta=\Delta\theta+O(\Delta\theta^3) sin Δ θ = Δ θ + O ( Δ θ 3 ) より、Δ θ \Delta\theta Δ θ の1次までで
v ( θ + Δ θ ) − v ( θ ) = n × v ( θ ) Δ θ + O ( Δ θ 2 ) \boldsymbol{v}(\theta+\Delta\theta)-\boldsymbol{v}(\theta)=\boldsymbol{n}\times\boldsymbol{v}(\theta)\,\Delta\theta+O(\Delta\theta^2) v ( θ + Δ θ ) − v ( θ ) = n × v ( θ ) Δ θ + O ( Δ θ 2 )
となり、これが示すべき式です。誤差は Δ θ 2 \Delta\theta^2 Δ θ 2 のオーダーで、軸に平行な向きには生じません。
幾何的には次の絵になります。v ( θ ) \boldsymbol{v}(\theta) v ( θ ) の終点は、軸上の点 u ∥ \boldsymbol{u}_{\parallel} u ∥ を中心とし軸に垂直な平面内にある半径 ρ = ∣ n × v ( θ ) ∣ \rho=|\boldsymbol{n}\times\boldsymbol{v}(\theta)| ρ = ∣ n × v ( θ ) ∣ の円周上を動きます。θ \theta θ が Δ θ \Delta\theta Δ θ 増えると終点はこの円周上を弧長 ρ Δ θ \rho\,\Delta\theta ρ Δ θ だけ進みます。進む向きは円の接線方向で、右ねじの正の向きに取れば単位接ベクトルは f 2 = ( n × v ) / ρ \boldsymbol{f}_2=(\boldsymbol{n}\times\boldsymbol{v})/\rho f 2 = ( n × v ) / ρ です。したがって変位は ρ Δ θ ⋅ ( n × v ) / ρ = n × v Δ θ \rho\,\Delta\theta\cdot(\boldsymbol{n}\times\boldsymbol{v})/\rho=\boldsymbol{n}\times\boldsymbol{v}\,\Delta\theta ρ Δ θ ⋅ ( n × v ) / ρ = n × v Δ θ となります。弧と弦の差が Δ θ \Delta\theta Δ θ の2次であることが、上の展開に対応します。
n v(θ) n × v(θ) Δθ ρ O
残る場合として u ⊥ = 0 \boldsymbol{u}_{\perp}=\boldsymbol{0} u ⊥ = 0 、つまり v ( θ ) \boldsymbol{v}(\theta) v ( θ ) が軸に平行なときは、終点が軸上にあって回転しても動かず v ( θ + Δ θ ) − v ( θ ) = 0 \boldsymbol{v}(\theta+\Delta\theta)-\boldsymbol{v}(\theta)=\boldsymbol{0} v ( θ + Δ θ ) − v ( θ ) = 0 です。一方 n × v ( θ ) = 0 \boldsymbol{n}\times\boldsymbol{v}(\theta)=\boldsymbol{0} n × v ( θ ) = 0 なので、示すべき式は両辺ともゼロベクトルとなって成立します。これで場合が尽くされました。
設問(ii)(3) の式を Δ θ \Delta\theta Δ θ で割って Δ θ → 0 \Delta\theta\to 0 Δ θ → 0 とすると、右辺の誤差項は O ( Δ θ ) → 0 O(\Delta\theta)\to 0 O ( Δ θ ) → 0 なので v ( θ ) \boldsymbol{v}(\theta) v ( θ ) は微分可能で
d v d θ = n × v ( θ ) = A ( n ) v ( θ ) \frac{d\boldsymbol{v}}{d\theta}=\boldsymbol{n}\times\boldsymbol{v}(\theta)=A(\boldsymbol{n})\,\boldsymbol{v}(\theta) d θ d v = n × v ( θ ) = A ( n ) v ( θ )
を満たします。これが定数係数の線形常微分方程式で、初期条件は v ( 0 ) = Ω ( n , 0 ) v 0 = v 0 \boldsymbol{v}(0)=\Omega(\boldsymbol{n},0)\boldsymbol{v}_0=\boldsymbol{v}_0 v ( 0 ) = Ω ( n , 0 ) v 0 = v 0 です。
これを解きます。A = A ( n ) A=A(\boldsymbol{n}) A = A ( n ) と略記します。行列指数関数 e θ A = ∑ k ≥ 0 θ k A k / k ! e^{\theta A}=\sum_{k\ge 0}\theta^kA^k/k! e θ A = ∑ k ≥ 0 θ k A k / k ! は任意の θ \theta θ で絶対収束し(∥ θ k A k / k ! ∥ ≤ ( ∣ θ ∣ ∥ A ∥ ) k / k ! \lVert \theta^kA^k/k!\rVert\le(\lvert\theta\rvert\lVert A\rVert)^k/k! ∥ θ k A k / k !∥ ≤ (∣ θ ∣ ∥ A ∥ ) k / k ! が収束級数で押さえられます)、項別微分して d d θ e θ A = A e θ A = e θ A A \frac{d}{d\theta}e^{\theta A}=Ae^{\theta A}=e^{\theta A}A d θ d e θ A = A e θ A = e θ A A が成り立ちます。また e θ A e − θ A = I e^{\theta A}e^{-\theta A}=I e θ A e − θ A = I なので e θ A e^{\theta A} e θ A は正則です。そこで
w ( θ ) = e − θ A v ( θ ) \boldsymbol{w}(\theta)=e^{-\theta A}\boldsymbol{v}(\theta) w ( θ ) = e − θ A v ( θ )
と置くと
d w d θ = − A e − θ A v + e − θ A d v d θ = e − θ A ( − A v + A v ) = 0 \frac{d\boldsymbol{w}}{d\theta}=-Ae^{-\theta A}\boldsymbol{v}+e^{-\theta A}\frac{d\boldsymbol{v}}{d\theta}=e^{-\theta A}\bigl(-A\boldsymbol{v}+A\boldsymbol{v}\bigr)=\boldsymbol{0} d θ d w = − A e − θ A v + e − θ A d θ d v = e − θ A ( − A v + A v ) = 0
となります(A A A と e − θ A e^{-\theta A} e − θ A が交換することを使いました)。よって w \boldsymbol{w} w は定ベクトルで w ( θ ) = w ( 0 ) = v 0 \boldsymbol{w}(\theta)=\boldsymbol{w}(0)=\boldsymbol{v}_0 w ( θ ) = w ( 0 ) = v 0 、つまり
Ω ( n , θ ) v 0 = v ( θ ) = e θ A ( n ) v 0 \Omega(\boldsymbol{n},\theta)\boldsymbol{v}_0=\boldsymbol{v}(\theta)=e^{\theta A(\boldsymbol{n})}\boldsymbol{v}_0 Ω ( n , θ ) v 0 = v ( θ ) = e θ A ( n ) v 0
です。この式は任意の v 0 ∈ R 3 \boldsymbol{v}_0\in\mathbb{R}^3 v 0 ∈ R 3 について成立するので、v 0 = e 1 , e 2 , e 3 \boldsymbol{v}_0=\boldsymbol{e}_1,\boldsymbol{e}_2,\boldsymbol{e}_3 v 0 = e 1 , e 2 , e 3 と取れば両辺の行列の各列が一致します。ゆえに行列として
Ω ( n , θ ) = e θ A ( n ) \Omega(\boldsymbol{n},\theta)=e^{\theta A(\boldsymbol{n})} Ω ( n , θ ) = e θ A ( n )
が成り立ちます。
確認しておきます。n = ( 0 , 0 , 1 ) \boldsymbol{n}=(0,0,1) n = ( 0 , 0 , 1 ) なら、設問(ii)(5) で求める対応から
A ( n ) = ( 0 − 1 0 1 0 0 0 0 0 ) , A ( n ) 2 = ( − 1 0 0 0 − 1 0 0 0 0 ) , A ( n ) 3 = − A ( n ) A(\boldsymbol{n})=\begin{pmatrix}0&-1&0\\1&0&0\\0&0&0\end{pmatrix},\qquad A(\boldsymbol{n})^2=\begin{pmatrix}-1&0&0\\0&-1&0\\0&0&0\end{pmatrix},\qquad A(\boldsymbol{n})^3=-A(\boldsymbol{n}) A ( n ) = 0 1 0 − 1 0 0 0 0 0 , A ( n ) 2 = − 1 0 0 0 − 1 0 0 0 0 , A ( n ) 3 = − A ( n )
なので、指数関数の級数を偶数次と奇数次に分けると
e θ A = I + sin θ A + ( 1 − cos θ ) A 2 = ( cos θ − sin θ 0 sin θ cos θ 0 0 0 1 ) e^{\theta A}=I+\sin\theta\,A+(1-\cos\theta)A^2=\begin{pmatrix}\cos\theta&-\sin\theta&0\\ \sin\theta&\cos\theta&0\\ 0&0&1\end{pmatrix} e θ A = I + sin θ A + ( 1 − cos θ ) A 2 = cos θ sin θ 0 − sin θ cos θ 0 0 0 1
となり、設問(ii)(1) の Ω \Omega Ω と一致します。なお A ( n ) 3 = − A ( n ) A(\boldsymbol{n})^3=-A(\boldsymbol{n}) A ( n ) 3 = − A ( n ) は一般の単位ベクトル n \boldsymbol{n} n でも成り立ちます(n × ( n × u ) = ( n ⋅ u ) n − u \boldsymbol{n}\times(\boldsymbol{n}\times\boldsymbol{u})=(\boldsymbol{n}\cdot\boldsymbol{u})\boldsymbol{n}-\boldsymbol{u} n × ( n × u ) = ( n ⋅ u ) n − u から従います)。したがって上の3項の表式はどの軸でも使えて、これがロドリゲスの公式です。
A ( a ) A(\boldsymbol{a}) A ( a ) を成分で書きます。( a × u ) 1 = a 2 u 3 − a 3 u 2 (\boldsymbol{a}\times\boldsymbol{u})_1=a_2u_3-a_3u_2 ( a × u ) 1 = a 2 u 3 − a 3 u 2 などから
A ( a ) = ( 0 − a 3 a 2 a 3 0 − a 1 − a 2 a 1 0 ) A(\boldsymbol{a})=\begin{pmatrix}0&-a_3&a_2\\ a_3&0&-a_1\\ -a_2&a_1&0\end{pmatrix} A ( a ) = 0 a 3 − a 2 − a 3 0 a 1 a 2 − a 1 0
です。これを問題の X X X と成分ごとに比べます。( 1 , 2 ) (1,2) ( 1 , 2 ) 成分から − a 3 = x 12 -a_3=x_{12} − a 3 = x 12 、( 1 , 3 ) (1,3) ( 1 , 3 ) 成分から a 2 = x 13 a_2=x_{13} a 2 = x 13 、( 2 , 3 ) (2,3) ( 2 , 3 ) 成分から − a 1 = x 23 -a_1=x_{23} − a 1 = x 23 が決まります。交代性 x i j = − x j i x_{ij}=-x_{ji} x ij = − x j i を使えば残りの ( 2 , 1 ) , ( 3 , 1 ) , ( 3 , 2 ) (2,1),(3,1),(3,2) ( 2 , 1 ) , ( 3 , 1 ) , ( 3 , 2 ) 成分の条件も同じ3式に帰着し、矛盾なく解けます。よって答えは
a 0 = ( x 32 , x 13 , x 21 ) = ( − x 23 , − x 31 , − x 12 ) \boldsymbol{a}_0=(x_{32},\,x_{13},\,x_{21})=(-x_{23},\,-x_{31},\,-x_{12}) a 0 = ( x 32 , x 13 , x 21 ) = ( − x 23 , − x 31 , − x 12 )
です。添字を使えば ( a 0 ) i = − 1 2 ε i j k X j k (\boldsymbol{a}_0)_i=-\tfrac{1}{2}\varepsilon_{ijk}X_{jk} ( a 0 ) i = − 2 1 ε ij k X j k とも書けます。対角成分が 0 0 0 で独立成分が3個という交代行列の自由度が、ベクトルの3成分とちょうど対応しています。写像 a ↦ A ( a ) \boldsymbol{a}\mapsto A(\boldsymbol{a}) a ↦ A ( a ) は線形で、いま見たようにベクトル全体から実3次交代行列全体への全単射です。
X X X を実3次交代行列とし、Y = e X Y=e^X Y = e X とします。設問(ii)(5) より X = A ( a 0 ) X=A(\boldsymbol{a}_0) X = A ( a 0 ) となる実ベクトル a 0 \boldsymbol{a}_0 a 0 が取れます。
a 0 = 0 \boldsymbol{a}_0=\boldsymbol{0} a 0 = 0 のときは X = O X=O X = O 、Y = I Y=I Y = I で、I ∈ R I\in R I ∈ R なので主張は成り立ちます。
a 0 ≠ 0 \boldsymbol{a}_0\neq\boldsymbol{0} a 0 = 0 のときは θ = ∣ a 0 ∣ > 0 \theta=|\boldsymbol{a}_0|>0 θ = ∣ a 0 ∣ > 0 、n = a 0 / θ \boldsymbol{n}=\boldsymbol{a}_0/\theta n = a 0 / θ (単位ベクトル)と置きます。A ( a ) u = a × u A(\boldsymbol{a})\boldsymbol{u}=\boldsymbol{a}\times\boldsymbol{u} A ( a ) u = a × u は a \boldsymbol{a} a について線形なので A ( a 0 ) = A ( θ n ) = θ A ( n ) A(\boldsymbol{a}_0)=A(\theta\boldsymbol{n})=\theta A(\boldsymbol{n}) A ( a 0 ) = A ( θ n ) = θ A ( n ) です。よって設問(ii)(4) の結果が使えて
Y = e X = e θ A ( n ) = Ω ( n , θ ) Y=e^{X}=e^{\theta A(\boldsymbol{n})}=\Omega(\boldsymbol{n},\theta) Y = e X = e θ A ( n ) = Ω ( n , θ )
となります。これは軸 n \boldsymbol{n} n のまわりの角 θ \theta θ の回転を表す行列そのものなので、Y ∈ R Y\in R Y ∈ R です。以上で示せました。軸と回転角が n = a 0 / ∣ a 0 ∣ \boldsymbol{n}=\boldsymbol{a}_0/|\boldsymbol{a}_0| n = a 0 /∣ a 0 ∣ 、θ = ∣ a 0 ∣ \theta=|\boldsymbol{a}_0| θ = ∣ a 0 ∣ と具体的に読み取れる点まで込めての結論です。
行列の性質から直接確認することもできます。X X X が実なら級数の各項が実なので Y Y Y は実行列です。転置は級数の各項ごとに取れて
t Y = t ( e X ) = e t X = e − X = ( e X ) − 1 = Y − 1 {}^{t}Y={}^{t}\bigl(e^{X}\bigr)=e^{{}^{t}X}=e^{-X}=\bigl(e^{X}\bigr)^{-1}=Y^{-1} t Y = t ( e X ) = e t X = e − X = ( e X ) − 1 = Y − 1
なので t Y Y = Y t Y = I {}^{t}YY=Y{}^{t}Y=I t Y Y = Y t Y = I 、つまり Y Y Y は実直交行列です。さらに det e X = e Tr X \det e^{X}=e^{\operatorname{Tr}X} det e X = e Tr X と交代行列の対角成分が 0 0 0 であること(Tr X = 0 \operatorname{Tr}X=0 Tr X = 0 )から det Y = e 0 = 1 \det Y=e^{0}=1 det Y = e 0 = 1 です。設問(i)(d) で見たように直交性だけでは足りず、det Y = + 1 \det Y=+1 det Y = + 1 まで言えて初めて回転だと結論できます。ここまでで Y ∈ S O ( 3 ) Y\in SO(3) Y ∈ S O ( 3 ) が分かり、S O ( 3 ) SO(3) S O ( 3 ) の元が必ずある軸のまわりの回転であること(オイラーの定理)を使えば同じ結論になりますが、上の構成的な議論はその軸と角を明示的に与えている点で強い主張です。
n = 0 , 1 , 2 , … n=0,1,2,\dots n = 0 , 1 , 2 , … に対して
f n ( x ) = e x 2 d n d x n ( e − x 2 ) f_n(x)=e^{x^2}\frac{d^n}{dx^n}\bigl(e^{-x^2}\bigr) f n ( x ) = e x 2 d x n d n ( e − x 2 )
と定義します。以下では g ( x ) = e − x 2 g(x)=e^{-x^2} g ( x ) = e − x 2 と書きます。定義は e − x 2 f n ( x ) = g ( n ) ( x ) e^{-x^2}f_n(x)=g^{(n)}(x) e − x 2 f n ( x ) = g ( n ) ( x ) と同値で、積分の設問ではこの形をそのまま使います。設問(iv) 以降では x x x を複素変数 z z z に置き換えますが、e − z 2 e^{-z^2} e − z 2 が整関数なので f n f_n f n は複素平面全体で定義された多項式になります。
漸化式を先に用意します。定義から
f n + 1 ( x ) = e x 2 d d x ( e − x 2 f n ( x ) ) = e x 2 ( e − x 2 f n ′ ( x ) − 2 x e − x 2 f n ( x ) ) = f n ′ ( x ) − 2 x f n ( x ) f_{n+1}(x)=e^{x^2}\frac{d}{dx}\Bigl(e^{-x^2}f_n(x)\Bigr)=e^{x^2}\Bigl(e^{-x^2}f_n'(x)-2xe^{-x^2}f_n(x)\Bigr)=f_n'(x)-2xf_n(x) f n + 1 ( x ) = e x 2 d x d ( e − x 2 f n ( x ) ) = e x 2 ( e − x 2 f n ′ ( x ) − 2 x e − x 2 f n ( x ) ) = f n ′ ( x ) − 2 x f n ( x )
です。f 0 = 1 f_0=1 f 0 = 1 から順に
f 1 = − 2 x , f 2 = 4 x 2 − 2 , f 3 = − 8 x 3 + 12 x , f 4 = 16 x 4 − 48 x 2 + 12 f_1=-2x,\qquad f_2=4x^2-2,\qquad f_3=-8x^3+12x,\qquad f_4=16x^4-48x^2+12 f 1 = − 2 x , f 2 = 4 x 2 − 2 , f 3 = − 8 x 3 + 12 x , f 4 = 16 x 4 − 48 x 2 + 12
が得られます。エルミート多項式 H n H_n H n とは f n ( x ) = ( − 1 ) n H n ( x ) f_n(x)=(-1)^nH_n(x) f n ( x ) = ( − 1 ) n H n ( x ) という関係にあります。
もう1つ、以下で何度も使う減衰の事実を確認しておきます。g ( k ) ( x ) g^{(k)}(x) g ( k ) ( x ) は(x x x の多項式)× e − x 2 \times e^{-x^2} × e − x 2 の形なので、任意の m , k ≥ 0 m,k\ge 0 m , k ≥ 0 に対して x m g ( k ) ( x ) → 0 x^mg^{(k)}(x)\to 0 x m g ( k ) ( x ) → 0 (∣ x ∣ → ∞ |x|\to\infty ∣ x ∣ → ∞ )であり、この形の関数の R \mathbb{R} R 上の積分は絶対収束します。以下の部分積分の境界項はすべてこれで消えます。
n n n 次です。漸化式 f n + 1 = f n ′ − 2 x f n f_{n+1}=f_n'-2xf_n f n + 1 = f n ′ − 2 x f n による帰納法で示します。f 0 = 1 f_0=1 f 0 = 1 は 0 0 0 次で最高次係数は a 0 = 1 a_0=1 a 0 = 1 です。f n f_n f n が n n n 次で最高次係数 a n ≠ 0 a_n\neq 0 a n = 0 だとすると、右辺の f n ′ f_n' f n ′ は n − 1 n-1 n − 1 次以下、− 2 x f n -2xf_n − 2 x f n は n + 1 n+1 n + 1 次で最高次係数 − 2 a n -2a_n − 2 a n です。次数が異なるので打ち消しは起こらず、f n + 1 f_{n+1} f n + 1 は n + 1 n+1 n + 1 次で a n + 1 = − 2 a n ≠ 0 a_{n+1}=-2a_n\neq 0 a n + 1 = − 2 a n = 0 となります。よって deg f n = n \deg f_n=n deg f n = n で、最高次係数は a n = ( − 2 ) n a_n=(-2)^n a n = ( − 2 ) n です。答えは x x x の n n n 次多項式です。
e − x 2 f n ( x ) = g ( n ) ( x ) e^{-x^2}f_n(x)=g^{(n)}(x) e − x 2 f n ( x ) = g ( n ) ( x ) を使うと
∫ − ∞ ∞ x f n ( x ) e − x 2 d x = ∫ − ∞ ∞ x g ( n ) ( x ) d x \int_{-\infty}^{\infty}xf_n(x)e^{-x^2}dx=\int_{-\infty}^{\infty}x\,g^{(n)}(x)\,dx ∫ − ∞ ∞ x f n ( x ) e − x 2 d x = ∫ − ∞ ∞ x g ( n ) ( x ) d x
です。n > 1 n>1 n > 1 すなわち n ≥ 2 n\ge 2 n ≥ 2 のもとで部分積分すると
∫ − ∞ ∞ x g ( n ) ( x ) d x = [ x g ( n − 1 ) ( x ) ] − ∞ ∞ − ∫ − ∞ ∞ g ( n − 1 ) ( x ) d x = − [ g ( n − 2 ) ( x ) ] − ∞ ∞ = 0 \int_{-\infty}^{\infty}x\,g^{(n)}(x)\,dx=\Bigl[x\,g^{(n-1)}(x)\Bigr]_{-\infty}^{\infty}-\int_{-\infty}^{\infty}g^{(n-1)}(x)\,dx=-\Bigl[g^{(n-2)}(x)\Bigr]_{-\infty}^{\infty}=0 ∫ − ∞ ∞ x g ( n ) ( x ) d x = [ x g ( n − 1 ) ( x ) ] − ∞ ∞ − ∫ − ∞ ∞ g ( n − 1 ) ( x ) d x = − [ g ( n − 2 ) ( x ) ] − ∞ ∞ = 0
となります。1つ目の境界項は設定で確認した減衰から 0 0 0 です。2つ目では n − 1 ≥ 1 n-1\ge 1 n − 1 ≥ 1 なので g ( n − 1 ) g^{(n-1)} g ( n − 1 ) の原始関数として g ( n − 2 ) g^{(n-2)} g ( n − 2 ) (n = 2 n=2 n = 2 のときは g g g 自身)が取れ、その境界値も 0 0 0 です。答えは 0 0 0 です。
n = 1 n=1 n = 1 を除外している理由もはっきりします。n = 1 n=1 n = 1 では ∫ g ( 0 ) d x = ∫ e − x 2 d x = π \int g^{(0)}dx=\int e^{-x^2}dx=\sqrt{\pi} ∫ g ( 0 ) d x = ∫ e − x 2 d x = π が残り、積分値は − π ≠ 0 -\sqrt{\pi}\neq 0 − π = 0 です。n ≥ 2 n\ge 2 n ≥ 2 で消えるのは、重み e − x 2 e^{-x^2} e − x 2 についての直交性から f n f_n f n が 1 1 1 次以下の多項式と直交することの言い換えでもあります。
同じ書き換えで
∫ − ∞ ∞ x n f n ( x ) e − x 2 d x = ∫ − ∞ ∞ x n g ( n ) ( x ) d x \int_{-\infty}^{\infty}x^nf_n(x)e^{-x^2}dx=\int_{-\infty}^{\infty}x^n\,g^{(n)}(x)\,dx ∫ − ∞ ∞ x n f n ( x ) e − x 2 d x = ∫ − ∞ ∞ x n g ( n ) ( x ) d x
です。部分積分を n n n 回繰り返します。j j j 回目の境界項は [ d j − 1 d x j − 1 ( x n ) ⋅ g ( n − j ) ( x ) ] − ∞ ∞ \bigl[\frac{d^{j-1}}{dx^{j-1}}(x^n)\cdot g^{(n-j)}(x)\bigr]_{-\infty}^{\infty} [ d x j − 1 d j − 1 ( x n ) ⋅ g ( n − j ) ( x ) ] − ∞ ∞ の形で、設定の減衰からすべて 0 0 0 です。よって
∫ − ∞ ∞ x n g ( n ) ( x ) d x = ( − 1 ) n ∫ − ∞ ∞ d n d x n ( x n ) g ( x ) d x = ( − 1 ) n n ! ∫ − ∞ ∞ e − x 2 d x \int_{-\infty}^{\infty}x^n\,g^{(n)}(x)\,dx=(-1)^n\int_{-\infty}^{\infty}\frac{d^n}{dx^n}\bigl(x^n\bigr)\,g(x)\,dx=(-1)^n\,n!\int_{-\infty}^{\infty}e^{-x^2}dx ∫ − ∞ ∞ x n g ( n ) ( x ) d x = ( − 1 ) n ∫ − ∞ ∞ d x n d n ( x n ) g ( x ) d x = ( − 1 ) n n ! ∫ − ∞ ∞ e − x 2 d x
となり、答えは
∫ − ∞ ∞ x n f n ( x ) e − x 2 d x = ( − 1 ) n n ! π \int_{-\infty}^{\infty}x^nf_n(x)e^{-x^2}dx=(-1)^n\,n!\,\sqrt{\pi} ∫ − ∞ ∞ x n f n ( x ) e − x 2 d x = ( − 1 ) n n ! π
です。
小さい n n n で検算します。n = 0 n=0 n = 0 では ∫ e − x 2 d x = π \int e^{-x^2}dx=\sqrt{\pi} ∫ e − x 2 d x = π で、式は ( − 1 ) 0 0 ! π = π (-1)^0 0!\sqrt{\pi}=\sqrt{\pi} ( − 1 ) 0 0 ! π = π を与えます。n = 1 n=1 n = 1 では f 1 = − 2 x f_1=-2x f 1 = − 2 x から ∫ ( − 2 x 2 ) e − x 2 d x = − 2 ⋅ π 2 = − π \int(-2x^2)e^{-x^2}dx=-2\cdot\frac{\sqrt{\pi}}{2}=-\sqrt{\pi} ∫ ( − 2 x 2 ) e − x 2 d x = − 2 ⋅ 2 π = − π で、式の値 ( − 1 ) 1 1 ! π (-1)^1\,1!\,\sqrt{\pi} ( − 1 ) 1 1 ! π と合います。n = 2 n=2 n = 2 では f 2 = 4 x 2 − 2 f_2=4x^2-2 f 2 = 4 x 2 − 2 から ∫ x 2 ( 4 x 2 − 2 ) e − x 2 d x = 4 ⋅ 3 π 4 − 2 ⋅ π 2 = 2 π \int x^2(4x^2-2)e^{-x^2}dx=4\cdot\frac{3\sqrt{\pi}}{4}-2\cdot\frac{\sqrt{\pi}}{2}=2\sqrt{\pi} ∫ x 2 ( 4 x 2 − 2 ) e − x 2 d x = 4 ⋅ 4 3 π − 2 ⋅ 2 π = 2 π で、2 ! π 2!\sqrt{\pi} 2 ! π と一致します。
e − ω 2 e^{-\omega^2} e − ω 2 は整関数なので、与えられた表示は z z z を中心とする半径 1 1 1 の円周に対するコーシーの積分公式です。まず n n n 階微分の表示に持ち上げます。z z z を固定し C = { ω : ∣ ω − z ∣ = 1 } C=\{\omega:|\omega-z|=1\} C = { ω : ∣ ω − z ∣ = 1 } (反時計回り)と書きます。∣ ζ − z ∣ < 1 / 2 |\zeta-z|<1/2 ∣ ζ − z ∣ < 1/2 を満たす点 ζ \zeta ζ に対して、ζ \zeta ζ を中心とする半径 1 1 1 の円周は、e − ω 2 e^{-\omega^2} e − ω 2 が正則な領域(複素平面全体)の内部で C C C へ連続変形でき、どちらも ζ \zeta ζ を正の向きに1回だけ囲みます。したがってコーシーの定理から経路を C C C に取り替えられて
e − ζ 2 = 1 2 π i ∮ C e − ω 2 ω − ζ d ω ( ∣ ζ − z ∣ < 1 / 2 ) e^{-\zeta^2}=\frac{1}{2\pi i}\oint_{C}\frac{e^{-\omega^2}}{\omega-\zeta}d\omega\qquad\bigl(|\zeta-z|<1/2\bigr) e − ζ 2 = 2 π i 1 ∮ C ω − ζ e − ω 2 d ω ( ∣ ζ − z ∣ < 1/2 )
が成り立ちます。右辺の経路 C C C は ζ \zeta ζ に依らないので、ζ \zeta ζ について微分できます。実際、ω ∈ C \omega\in C ω ∈ C と ∣ ζ − z ∣ < 1 / 2 |\zeta-z|<1/2 ∣ ζ − z ∣ < 1/2 のとき ∣ ω − ζ ∣ ≥ 1 / 2 |\omega-\zeta|\ge 1/2 ∣ ω − ζ ∣ ≥ 1/2 なので分母が 0 0 0 に近づかず、被積分関数は ζ \zeta ζ について何回でも微分可能でその導関数は ω ∈ C \omega\in C ω ∈ C について一様有界です。C C C はコンパクトなので微分と積分は交換できます。ζ \zeta ζ で n n n 回微分して ζ = z \zeta=z ζ = z と置くと
d n d z n e − z 2 = n ! 2 π i ∮ C e − ω 2 ( ω − z ) n + 1 d ω , つまり f n ( z ) n ! = e z 2 2 π i ∮ C e − ω 2 ( ω − z ) n + 1 d ω \frac{d^n}{dz^n}e^{-z^2}=\frac{n!}{2\pi i}\oint_{C}\frac{e^{-\omega^2}}{(\omega-z)^{n+1}}d\omega,
\qquad\text{つまり}\qquad
\frac{f_n(z)}{n!}=\frac{e^{z^2}}{2\pi i}\oint_{C}\frac{e^{-\omega^2}}{(\omega-z)^{n+1}}d\omega d z n d n e − z 2 = 2 π i n ! ∮ C ( ω − z ) n + 1 e − ω 2 d ω , つまり n ! f n ( z ) = 2 π i e z 2 ∮ C ( ω − z ) n + 1 e − ω 2 d ω
です(グルサの公式)。
これを級数に代入します。
∑ n = 0 ∞ t n n ! f n ( z ) = e z 2 1 2 π i ∮ C e − ω 2 ∑ n = 0 ∞ t n ( ω − z ) n + 1 d ω \sum_{n=0}^{\infty}\frac{t^n}{n!}f_n(z)=e^{z^2}\,\frac{1}{2\pi i}\oint_{C}e^{-\omega^2}\sum_{n=0}^{\infty}\frac{t^n}{(\omega-z)^{n+1}}\,d\omega n = 0 ∑ ∞ n ! t n f n ( z ) = e z 2 2 π i 1 ∮ C e − ω 2 n = 0 ∑ ∞ ( ω − z ) n + 1 t n d ω
順序交換は問題の指示で認められていますが、根拠も明確です。C C C 上では ∣ ω − z ∣ = 1 |\omega-z|=1 ∣ ω − z ∣ = 1 なので ∣ t / ( ω − z ) ∣ = ∣ t ∣ < 1 \bigl|t/(\omega-z)\bigr|=|t|<1 t / ( ω − z ) = ∣ t ∣ < 1 が ω \omega ω に依らない定数で押さえられ、級数は C C C 上一様絶対収束します。e − ω 2 e^{-\omega^2} e − ω 2 は C C C 上有界なので、有界な連続関数を掛けても一様収束は保たれ、コンパクトな経路上での項別積分が正当化されます。
等比級数を足すと
∑ n = 0 ∞ t n ( ω − z ) n + 1 = 1 ω − z ⋅ 1 1 − t ω − z = 1 ω − ( z + t ) \sum_{n=0}^{\infty}\frac{t^n}{(\omega-z)^{n+1}}=\frac{1}{\omega-z}\cdot\frac{1}{1-\dfrac{t}{\omega-z}}=\frac{1}{\omega-(z+t)} n = 0 ∑ ∞ ( ω − z ) n + 1 t n = ω − z 1 ⋅ 1 − ω − z t 1 = ω − ( z + t ) 1
なので
∑ n = 0 ∞ t n n ! f n ( z ) = e z 2 1 2 π i ∮ C e − ω 2 ω − ( z + t ) d ω \sum_{n=0}^{\infty}\frac{t^n}{n!}f_n(z)=e^{z^2}\,\frac{1}{2\pi i}\oint_{C}\frac{e^{-\omega^2}}{\omega-(z+t)}\,d\omega n = 0 ∑ ∞ n ! t n f n ( z ) = e z 2 2 π i 1 ∮ C ω − ( z + t ) e − ω 2 d ω
となります。∣ t ∣ < 1 |t|<1 ∣ t ∣ < 1 より点 z + t z+t z + t は円 C C C の内部にあり、e − ω 2 e^{-\omega^2} e − ω 2 は整関数なので、右辺の積分にコーシーの積分公式が使えて値は e − ( z + t ) 2 e^{-(z+t)^2} e − ( z + t ) 2 です。よって
∑ n = 0 ∞ t n n ! f n ( z ) = e z 2 e − ( z + t ) 2 = e − t 2 − 2 t z \sum_{n=0}^{\infty}\frac{t^n}{n!}f_n(z)=e^{z^2}e^{-(z+t)^2}=e^{-t^2-2tz} n = 0 ∑ ∞ n ! t n f n ( z ) = e z 2 e − ( z + t ) 2 = e − t 2 − 2 t z
が示せました。仮定 ∣ t ∣ < 1 |t|<1 ∣ t ∣ < 1 は、等比級数の収束と z + t z+t z + t が C C C の内部に入ることの両方に使われています。
低次で確認しておきます。右辺を t t t で展開すると e − t 2 − 2 t z = 1 − 2 t z + ( 2 z 2 − 1 ) t 2 + ⋯ e^{-t^2-2tz}=1-2tz+(2z^2-1)t^2+\cdots e − t 2 − 2 t z = 1 − 2 t z + ( 2 z 2 − 1 ) t 2 + ⋯ で、左辺の t 0 , t 1 , t 2 t^0,t^1,t^2 t 0 , t 1 , t 2 の係数は f 0 = 1 f_0=1 f 0 = 1 、f 1 = − 2 z f_1=-2z f 1 = − 2 z 、f 2 / 2 ! = ( 4 z 2 − 2 ) / 2 = 2 z 2 − 1 f_2/2!=(4z^2-2)/2=2z^2-1 f 2 /2 ! = ( 4 z 2 − 2 ) /2 = 2 z 2 − 1 です。一致します。
F ( t , z ) = e − t 2 − 2 t z F(t,z)=e^{-t^2-2tz} F ( t , z ) = e − t 2 − 2 t z と置きます。与えられた式は直接確かめられます。∂ F / ∂ z = − 2 t F \partial F/\partial z=-2tF ∂ F / ∂ z = − 2 tF 、∂ 2 F / ∂ z 2 = 4 t 2 F \partial^2F/\partial z^2=4t^2F ∂ 2 F / ∂ z 2 = 4 t 2 F 、∂ F / ∂ t = − ( 2 t + 2 z ) F \partial F/\partial t=-(2t+2z)F ∂ F / ∂ t = − ( 2 t + 2 z ) F なので
∂ 2 F ∂ z 2 − 2 z ∂ F ∂ z = ( 4 t 2 + 4 t z ) F = − 2 t ( − 2 t − 2 z ) F = − 2 t ∂ F ∂ t \frac{\partial^2F}{\partial z^2}-2z\frac{\partial F}{\partial z}=\bigl(4t^2+4tz\bigr)F=-2t\bigl(-2t-2z\bigr)F=-2t\frac{\partial F}{\partial t} ∂ z 2 ∂ 2 F − 2 z ∂ z ∂ F = ( 4 t 2 + 4 t z ) F = − 2 t ( − 2 t − 2 z ) F = − 2 t ∂ t ∂ F
です。
設問(iv) より ∣ t ∣ < 1 |t|<1 ∣ t ∣ < 1 で F ( t , z ) = ∑ m ≥ 0 t m m ! f m ( z ) F(t,z)=\sum_{m\ge 0}\frac{t^m}{m!}f_m(z) F ( t , z ) = ∑ m ≥ 0 m ! t m f m ( z ) であり、収束冪級数の係数は一意なので
f n ( z ) = ∂ n F ∂ t n ∣ t = 0 f_n(z)=\frac{\partial^nF}{\partial t^n}\bigg|_{t=0} f n ( z ) = ∂ t n ∂ n F t = 0
が成り立ちます。そこで上の式の両辺を t t t で n n n 回微分して t = 0 t=0 t = 0 と置きます。F F F は ( t , z ) (t,z) ( t , z ) の整関数なので偏微分の順序はすべて交換できます。
左辺は z z z 微分と t t t 微分を入れ替えて
∂ n ∂ t n ( ∂ 2 F ∂ z 2 − 2 z ∂ F ∂ z ) ∣ t = 0 = d 2 d z 2 f n ( z ) − 2 z d d z f n ( z ) \frac{\partial^n}{\partial t^n}\Bigl(\frac{\partial^2F}{\partial z^2}-2z\frac{\partial F}{\partial z}\Bigr)\bigg|_{t=0}
=\frac{d^2}{dz^2}f_n(z)-2z\frac{d}{dz}f_n(z) ∂ t n ∂ n ( ∂ z 2 ∂ 2 F − 2 z ∂ z ∂ F ) t = 0 = d z 2 d 2 f n ( z ) − 2 z d z d f n ( z )
となります。右辺はライプニッツ則から
∂ n ∂ t n ( t ∂ F ∂ t ) = t ∂ n + 1 F ∂ t n + 1 + n ∂ n F ∂ t n \frac{\partial^n}{\partial t^n}\Bigl(t\frac{\partial F}{\partial t}\Bigr)=t\frac{\partial^{n+1}F}{\partial t^{n+1}}+n\frac{\partial^{n}F}{\partial t^{n}} ∂ t n ∂ n ( t ∂ t ∂ F ) = t ∂ t n + 1 ∂ n + 1 F + n ∂ t n ∂ n F
なので、t = 0 t=0 t = 0 で第1項が落ちて n f n ( z ) n f_n(z) n f n ( z ) が残り、全体で − 2 n f n ( z ) -2nf_n(z) − 2 n f n ( z ) です。両者を等号で結べば
d 2 d z 2 f n ( z ) − 2 z d d z f n ( z ) + 2 n f n ( z ) = 0 \frac{d^2}{dz^2}f_n(z)-2z\frac{d}{dz}f_n(z)+2n\,f_n(z)=0 d z 2 d 2 f n ( z ) − 2 z d z d f n ( z ) + 2 n f n ( z ) = 0
が得られます。これが示すべき微分方程式で、答えは λ = 2 n \lambda=2n λ = 2 n です。
λ \lambda λ の値は最高次の項だけを見ても確認できます。f n = a n z n + ⋯ f_n=a_nz^n+\cdots f n = a n z n + ⋯ (a n = ( − 2 ) n ≠ 0 a_n=(-2)^n\neq 0 a n = ( − 2 ) n = 0 )を代入すると、z n z^n z n の係数は f n ′ ′ f_n'' f n ′′ から寄与がなく、− 2 z f n ′ -2zf_n' − 2 z f n ′ から − 2 n a n -2na_n − 2 n a n 、λ f n \lambda f_n λ f n から λ a n \lambda a_n λ a n なので、( λ − 2 n ) a n = 0 (\lambda-2n)a_n=0 ( λ − 2 n ) a n = 0 すなわち λ = 2 n \lambda=2n λ = 2 n が必要です。多項式全体で成り立つことは上の議論が保証します。具体例でも確かめられます。f 2 = 4 z 2 − 2 f_2=4z^2-2 f 2 = 4 z 2 − 2 では f 2 ′ ′ = 8 f_2''=8 f 2 ′′ = 8 、− 2 z f 2 ′ = − 16 z 2 -2zf_2'=-16z^2 − 2 z f 2 ′ = − 16 z 2 、4 f 2 = 16 z 2 − 8 4f_2=16z^2-8 4 f 2 = 16 z 2 − 8 の和が 0 0 0 です。f 3 = − 8 z 3 + 12 z f_3=-8z^3+12z f 3 = − 8 z 3 + 12 z では f 3 ′ ′ = − 48 z f_3''=-48z f 3 ′′ = − 48 z 、− 2 z f 3 ′ = 48 z 3 − 24 z -2zf_3'=48z^3-24z − 2 z f 3 ′ = 48 z 3 − 24 z 、6 f 3 = − 48 z 3 + 72 z 6f_3=-48z^3+72z 6 f 3 = − 48 z 3 + 72 z の和が 0 0 0 です。λ = 2 n \lambda=2n λ = 2 n はエルミートの微分方程式の固有値で、設定で述べた f n = ( − 1 ) n H n f_n=(-1)^nH_n f n = ( − 1 ) n H n という対応と整合しています。
出典: 東京大学大学院理学系研究科 物理学専攻 平成20年度 修士課程 入学試験問題 数学。問題文は要約して引用しています。