コンテンツにスキップ

平成20年度 東大院 物理学専攻 修士 数学 解答

生 Markdown

大問2題で、2問すべてが必答です(試験時間は60分)。第1問は3次元回転の集合 SO(3)SO(3) を題材に、群としての性質の真偽判定から始めて、回転の固有値、無限小回転の生成子、交代行列の指数関数が回転になることまでを一本の線でたどります。第2問は fn(x)=ex2(d/dx)nex2f_n(x)=e^{x^2}(d/dx)^n e^{-x^2} で定義される多項式(符号を除いてエルミート多項式)の問題で、部分積分によるガウス型積分、コーシーの積分表示から母関数を導く部分、母関数が満たす偏微分方程式から各 fnf_n の常微分方程式を取り出す部分が山になります。どちらの問題も「示せ」「説明せよ」が多く、極限操作や場合分けの詰めがそのまま得点になります。

問題分野主題
第1問線形代数・群論・微分方程式3次元回転群、回転の固有値、交代行列と行列指数関数
第2問特殊関数・複素解析・微分方程式エルミート多項式、ガウス型積分、複素積分による母関数

第1問 3次元回転群と行列指数関数

Section titled “第1問 3次元回転群と行列指数関数”

RR は3次元空間の回転を表す実 3×33\times 3 行列全体の集合で、単位行列 II も含みます。ここで回転とは、原点を固定する線形変換のうち、長さ(したがって内積)を保ち、かつ右手系を右手系に写す(行列式が正の)もの、という意味で使います。この解釈のもとで R=SO(3)R=SO(3) です。

Ω(n,θ)\Omega(\boldsymbol{n},\theta) は単位ベクトル n\boldsymbol{n} を軸とする角 θ\theta の回転行列で、回転の正の向きは n\boldsymbol{n} の方向に進む右ねじを回す向きとします。任意の実3次元ベクトル v0\boldsymbol{v}_0 に対して v(θ)=Ω(n,θ)v0\boldsymbol{v}(\theta)=\Omega(\boldsymbol{n},\theta)\boldsymbol{v}_0 と置きます。また A(a)A(\boldsymbol{a}) は、任意の u\boldsymbol{u} に対して A(a)u=a×uA(\boldsymbol{a})\boldsymbol{u}=\boldsymbol{a}\times\boldsymbol{u} を満たす実 3×33\times 3 行列です。標準基底を e1,e2,e3\boldsymbol{e}_1,\boldsymbol{e}_2,\boldsymbol{e}_3 と書きます。

以下で繰り返し使う道具を1つ先に用意します。単位ベクトル n\boldsymbol{n} に対して、(f1,f2,n)(\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) が右手系の正規直交系(つまり f1×f2=n\boldsymbol{f}_1\times\boldsymbol{f}_2=\boldsymbol{n}n×f1=f2\boldsymbol{n}\times\boldsymbol{f}_1=\boldsymbol{f}_2f2×n=f1\boldsymbol{f}_2\times\boldsymbol{n}=\boldsymbol{f}_1)となるように f1,f2\boldsymbol{f}_1,\boldsymbol{f}_2 を取ると、Ω(n,θ)\Omega(\boldsymbol{n},\theta) の作用は

Ω(n,θ)n=n,Ω(n,θ)f1=cosθf1+sinθf2,Ω(n,θ)f2=sinθf1+cosθf2\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

で尽くされます。これは「軸を動かさず、軸に垂直な平面内を右ねじの向きに角 θ\theta だけ回す」という定義をそのまま書き下したものです。

(a) ○ です。X,YX,Y が長さを保つなら合成 XYXY も長さを保ちます。また det(XY)=detXdetY>0\det(XY)=\det X\det Y>0 なので向きも保ちます。よって XYXY も回転で、RR は行列の積について閉じています(実際 RR は群をなします)。

(b) × です。軸の異なる回転は一般に交換しません。X=Ω(e1,π/2)X=\Omega(\boldsymbol{e}_1,\pi/2)Y=Ω(e3,π/2)Y=\Omega(\boldsymbol{e}_3,\pi/2)、すなわち

X=(100001010),Y=(010100001)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}

を取ると

XY=(010001100),YX=(001100010)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}

で、たとえば XYe1=(0,0,1)XY\boldsymbol{e}_1=(0,0,1) に対し YXe1=(0,1,0)YX\boldsymbol{e}_1=(0,1,0) となり一致しません。

(c) ○ です。XRX\in R は内積を保つので、任意の u,w\boldsymbol{u},\boldsymbol{w} に対して

tu(tXX)w=t(Xu)(Xw)=tuw{}^{t}\boldsymbol{u}\,({}^{t}XX)\,\boldsymbol{w}={}^{t}(X\boldsymbol{u})(X\boldsymbol{w})={}^{t}\boldsymbol{u}\,\boldsymbol{w}

が成り立ちます。u=ei\boldsymbol{u}=\boldsymbol{e}_iw=ej\boldsymbol{w}=\boldsymbol{e}_j と取れば (tXX)ij=δij({}^{t}XX)_{ij}=\delta_{ij}、つまり tXX=I{}^{t}XX=I です。これで XX は正則で逆行列が tX{}^{t}X と分かるので、XtX=IX{}^{t}X=I も従います。

(d) × です。直交性だけでは行列式の符号が決まりません。反例は鏡映

X=(100010001)X=\begin{pmatrix}1&0&0\\0&1&0\\0&0&-1\end{pmatrix}

で、tXX=XtX=I{}^{t}XX=X{}^{t}X=I を満たしますが detX=1\det X=-1 であり、右手系を左手系に写すので回転ではありません(X=IX=-I でも同じです)。直交行列全体 O(3)O(3) のうち det=+1\det=+1 の部分だけが R=SO(3)R=SO(3) です。

(e) ○ です。(c) より XX は実直交行列なので X=tX=X1X^{\dagger}={}^{t}X=X^{-1} であり、XX=XX=IX^{\dagger}X=XX^{\dagger}=I が成り立ちます。すなわち XX はユニタリ行列であり、特に正規行列です。正規行列はユニタリ行列で対角化できる(スペクトル定理)ので、UXUU^{\dagger}XU が対角行列となるユニタリ UU が存在します。UU の具体形は設問(ii)(2)で作ります。

まとめると (a)○、(b)×、(c)○、(d)×、(e)○ が答えです。

n=(0,0,1)\boldsymbol{n}=(0,0,1) のときは設定で述べた作用そのままで、f1=e1\boldsymbol{f}_1=\boldsymbol{e}_1f2=e2\boldsymbol{f}_2=\boldsymbol{e}_2 と取れて

Ω(n,θ)=(cosθsinθ0sinθcosθ0001)\Omega(\boldsymbol{n},\theta)=\begin{pmatrix}\cos\theta&-\sin\theta&0\\ \sin\theta&\cos\theta&0\\ 0&0&1\end{pmatrix}

です。固有方程式は

det(ΩλI)=(1λ){(cosθλ)2+sin2θ}=(1λ)(λ22λ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

となり、λ22λcosθ+1=0\lambda^2-2\lambda\cos\theta+1=0 の根は λ=cosθ±isinθ\lambda=\cos\theta\pm i\sin\theta です。答えは

λ=1,eiθ,eiθ\lambda=1,\quad e^{i\theta},\quad e^{-i\theta}

の3つです。検算として TrΩ=1+2cosθ\operatorname{Tr}\Omega=1+2\cos\theta が固有値の和 1+eiθ+eiθ1+e^{i\theta}+e^{-i\theta} に、detΩ=1\det\Omega=1 が積 1eiθeiθ1\cdot e^{i\theta}\cdot e^{-i\theta} に一致します。

任意の単位ベクトル n\boldsymbol{n} についても固有値は同じで、答えは λ=1, eiθ, eiθ\lambda=1,\ e^{i\theta},\ e^{-i\theta} です。

理由は基底の取り替えです。設定で用意した右手系正規直交系 (f1,f2,n)(\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) を列に並べた行列を P=(f1 f2 n)P=(\boldsymbol{f}_1\ \boldsymbol{f}_2\ \boldsymbol{n}) とすると PP は直交行列で、Ω(n,θ)\Omega(\boldsymbol{n},\theta) のこの基底での表現行列は設定で書いた作用より設問(ii)(1) と同じ形の行列 M(θ)M(\theta) になります。すなわち

Ω(n,θ)=PM(θ)tP,M(θ)=(cosθsinθ0sinθcosθ0001)\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}

です。相似変換で固有値は変わらないので、設問(ii)(1) の結果がそのまま使えます。

固有ベクトルも書き下せます。設定の作用から

Ω(n,θ)(f1if2)=e±iθ(f1if2),Ω(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}

が直接確かめられます。(f1if2)/2(\boldsymbol{f}_1-i\boldsymbol{f}_2)/\sqrt{2}(f1+if2)/2(\boldsymbol{f}_1+i\boldsymbol{f}_2)/\sqrt{2}n\boldsymbol{n} はエルミート内積について正規直交なので、これらを列に並べた UU はユニタリで

UΩ(n,θ)U=(eiθ000eiθ0001)U^{\dagger}\,\Omega(\boldsymbol{n},\theta)\,U=\begin{pmatrix}e^{i\theta}&0&0\\0&e^{-i\theta}&0\\0&0&1\end{pmatrix}

となります。これが設問(i)(e) の UU の具体形です。

同じ軸のまわりの回転は角度が足し算になります。実際、設定の基底 (f1,f2,n)(\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) では Ω(n,θ)\Omega(\boldsymbol{n},\theta)M(θ)M(\theta) で表され、三角関数の加法定理から M(Δθ)M(θ)=M(θ+Δθ)M(\Delta\theta)M(\theta)=M(\theta+\Delta\theta) が成り立つので

v(θ+Δθ)=Ω(n,θ+Δθ)v0=Ω(n,Δθ)v(θ)\boldsymbol{v}(\theta+\Delta\theta)=\Omega(\boldsymbol{n},\theta+\Delta\theta)\boldsymbol{v}_0=\Omega(\boldsymbol{n},\Delta\theta)\,\boldsymbol{v}(\theta)

です。つまり、θ\theta から θ+Δθ\theta+\Delta\theta への変化は、いまの位置 v(θ)\boldsymbol{v}(\theta) を軸 n\boldsymbol{n} のまわりに Δθ\Delta\theta だけ回すことに等しくなります。

そこで u=v(θ)\boldsymbol{u}=\boldsymbol{v}(\theta) を軸に平行な成分と垂直な成分に分けます。

u=(nu)n,u=uu\boldsymbol{u}_{\parallel}=(\boldsymbol{n}\cdot\boldsymbol{u})\boldsymbol{n},\qquad \boldsymbol{u}_{\perp}=\boldsymbol{u}-\boldsymbol{u}_{\parallel}

u0\boldsymbol{u}_{\perp}\neq\boldsymbol{0} の場合を先に扱います。ρ=u>0\rho=|\boldsymbol{u}_{\perp}|>0 と置き、f1=u/ρ\boldsymbol{f}_1=\boldsymbol{u}_{\perp}/\rhof2=n×f1\boldsymbol{f}_2=\boldsymbol{n}\times\boldsymbol{f}_1 と取れば (f1,f2,n)(\boldsymbol{f}_1,\boldsymbol{f}_2,\boldsymbol{n}) は右手系の正規直交系です。回転は軸方向成分を動かさず、垂直成分を f1f2\boldsymbol{f}_1\boldsymbol{f}_2 平面内で Δθ\Delta\theta 回すので

Ω(n,Δθ)u=u+ρ(cosΔθf1+sinΔθf2)\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)

となります。ここで ρf2=n×u=n×u\rho\boldsymbol{f}_2=\boldsymbol{n}\times\boldsymbol{u}_{\perp}=\boldsymbol{n}\times\boldsymbol{u}n×u=0\boldsymbol{n}\times\boldsymbol{u}_{\parallel}=\boldsymbol{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)

が厳密に成り立ちます。cosΔθ1=Δθ2/2+O(Δθ4)\cos\Delta\theta-1=-\Delta\theta^2/2+O(\Delta\theta^4)sinΔθ=Δθ+O(Δθ3)\sin\Delta\theta=\Delta\theta+O(\Delta\theta^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)

となり、これが示すべき式です。誤差は Δθ2\Delta\theta^2 のオーダーで、軸に平行な向きには生じません。

幾何的には次の絵になります。v(θ)\boldsymbol{v}(\theta) の終点は、軸上の点 u\boldsymbol{u}_{\parallel} を中心とし軸に垂直な平面内にある半径 ρ=n×v(θ)\rho=|\boldsymbol{n}\times\boldsymbol{v}(\theta)| の円周上を動きます。θ\thetaΔθ\Delta\theta 増えると終点はこの円周上を弧長 ρΔθ\rho\,\Delta\theta だけ進みます。進む向きは円の接線方向で、右ねじの正の向きに取れば単位接ベクトルは f2=(n×v)/ρ\boldsymbol{f}_2=(\boldsymbol{n}\times\boldsymbol{v})/\rho です。したがって変位は ρΔθ(n×v)/ρ=n×vΔθ\rho\,\Delta\theta\cdot(\boldsymbol{n}\times\boldsymbol{v})/\rho=\boldsymbol{n}\times\boldsymbol{v}\,\Delta\theta となります。弧と弦の差が Δθ\Delta\theta の2次であることが、上の展開に対応します。

nv(θ)n × v(θ) ΔθρO

残る場合として u=0\boldsymbol{u}_{\perp}=\boldsymbol{0}、つまり v(θ)\boldsymbol{v}(\theta) が軸に平行なときは、終点が軸上にあって回転しても動かず v(θ+Δθ)v(θ)=0\boldsymbol{v}(\theta+\Delta\theta)-\boldsymbol{v}(\theta)=\boldsymbol{0} です。一方 n×v(θ)=0\boldsymbol{n}\times\boldsymbol{v}(\theta)=\boldsymbol{0} なので、示すべき式は両辺ともゼロベクトルとなって成立します。これで場合が尽くされました。

設問(ii)(3) の式を Δθ\Delta\theta で割って Δθ0\Delta\theta\to 0 とすると、右辺の誤差項は O(Δθ)0O(\Delta\theta)\to 0 なので v(θ)\boldsymbol{v}(\theta) は微分可能で

dvdθ=n×v(θ)=A(n)v(θ)\frac{d\boldsymbol{v}}{d\theta}=\boldsymbol{n}\times\boldsymbol{v}(\theta)=A(\boldsymbol{n})\,\boldsymbol{v}(\theta)

を満たします。これが定数係数の線形常微分方程式で、初期条件は v(0)=Ω(n,0)v0=v0\boldsymbol{v}(0)=\Omega(\boldsymbol{n},0)\boldsymbol{v}_0=\boldsymbol{v}_0 です。

これを解きます。A=A(n)A=A(\boldsymbol{n}) と略記します。行列指数関数 eθA=k0θkAk/k!e^{\theta A}=\sum_{k\ge 0}\theta^kA^k/k! は任意の θ\theta で絶対収束し(θkAk/k!(θA)k/k!\lVert \theta^kA^k/k!\rVert\le(\lvert\theta\rvert\lVert A\rVert)^k/k! が収束級数で押さえられます)、項別微分して ddθeθA=AeθA=eθAA\frac{d}{d\theta}e^{\theta A}=Ae^{\theta A}=e^{\theta A}A が成り立ちます。また eθAeθA=Ie^{\theta A}e^{-\theta A}=I なので eθAe^{\theta A} は正則です。そこで

w(θ)=eθAv(θ)\boldsymbol{w}(\theta)=e^{-\theta A}\boldsymbol{v}(\theta)

と置くと

dwdθ=AeθAv+eθAdvdθ=eθA(Av+Av)=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}

となります(AAeθAe^{-\theta A} が交換することを使いました)。よって w\boldsymbol{w} は定ベクトルで w(θ)=w(0)=v0\boldsymbol{w}(\theta)=\boldsymbol{w}(0)=\boldsymbol{v}_0、つまり

Ω(n,θ)v0=v(θ)=eθA(n)v0\Omega(\boldsymbol{n},\theta)\boldsymbol{v}_0=\boldsymbol{v}(\theta)=e^{\theta A(\boldsymbol{n})}\boldsymbol{v}_0

です。この式は任意の v0R3\boldsymbol{v}_0\in\mathbb{R}^3 について成立するので、v0=e1,e2,e3\boldsymbol{v}_0=\boldsymbol{e}_1,\boldsymbol{e}_2,\boldsymbol{e}_3 と取れば両辺の行列の各列が一致します。ゆえに行列として

Ω(n,θ)=eθA(n)\Omega(\boldsymbol{n},\theta)=e^{\theta A(\boldsymbol{n})}

が成り立ちます。

確認しておきます。n=(0,0,1)\boldsymbol{n}=(0,0,1) なら、設問(ii)(5) で求める対応から

A(n)=(010100000),A(n)2=(100010000),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})

なので、指数関数の級数を偶数次と奇数次に分けると

eθA=I+sinθA+(1cosθ)A2=(cosθsinθ0sinθcosθ0001)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}

となり、設問(ii)(1) の Ω\Omega と一致します。なお A(n)3=A(n)A(\boldsymbol{n})^3=-A(\boldsymbol{n}) は一般の単位ベクトル n\boldsymbol{n} でも成り立ちます(n×(n×u)=(nu)nu\boldsymbol{n}\times(\boldsymbol{n}\times\boldsymbol{u})=(\boldsymbol{n}\cdot\boldsymbol{u})\boldsymbol{n}-\boldsymbol{u} から従います)。したがって上の3項の表式はどの軸でも使えて、これがロドリゲスの公式です。

A(a)A(\boldsymbol{a}) を成分で書きます。(a×u)1=a2u3a3u2(\boldsymbol{a}\times\boldsymbol{u})_1=a_2u_3-a_3u_2 などから

A(a)=(0a3a2a30a1a2a10)A(\boldsymbol{a})=\begin{pmatrix}0&-a_3&a_2\\ a_3&0&-a_1\\ -a_2&a_1&0\end{pmatrix}

です。これを問題の XX と成分ごとに比べます。(1,2)(1,2) 成分から a3=x12-a_3=x_{12}(1,3)(1,3) 成分から a2=x13a_2=x_{13}(2,3)(2,3) 成分から a1=x23-a_1=x_{23} が決まります。交代性 xij=xjix_{ij}=-x_{ji} を使えば残りの (2,1),(3,1),(3,2)(2,1),(3,1),(3,2) 成分の条件も同じ3式に帰着し、矛盾なく解けます。よって答えは

a0=(x32,x13,x21)=(x23,x31,x12)\boldsymbol{a}_0=(x_{32},\,x_{13},\,x_{21})=(-x_{23},\,-x_{31},\,-x_{12})

です。添字を使えば (a0)i=12εijkXjk(\boldsymbol{a}_0)_i=-\tfrac{1}{2}\varepsilon_{ijk}X_{jk} とも書けます。対角成分が 00 で独立成分が3個という交代行列の自由度が、ベクトルの3成分とちょうど対応しています。写像 aA(a)\boldsymbol{a}\mapsto A(\boldsymbol{a}) は線形で、いま見たようにベクトル全体から実3次交代行列全体への全単射です。

XX を実3次交代行列とし、Y=eXY=e^X とします。設問(ii)(5) より X=A(a0)X=A(\boldsymbol{a}_0) となる実ベクトル a0\boldsymbol{a}_0 が取れます。

a0=0\boldsymbol{a}_0=\boldsymbol{0} のときは X=OX=OY=IY=I で、IRI\in R なので主張は成り立ちます。

a00\boldsymbol{a}_0\neq\boldsymbol{0} のときは θ=a0>0\theta=|\boldsymbol{a}_0|>0n=a0/θ\boldsymbol{n}=\boldsymbol{a}_0/\theta(単位ベクトル)と置きます。A(a)u=a×uA(\boldsymbol{a})\boldsymbol{u}=\boldsymbol{a}\times\boldsymbol{u}a\boldsymbol{a} について線形なので A(a0)=A(θn)=θA(n)A(\boldsymbol{a}_0)=A(\theta\boldsymbol{n})=\theta A(\boldsymbol{n}) です。よって設問(ii)(4) の結果が使えて

Y=eX=eθA(n)=Ω(n,θ)Y=e^{X}=e^{\theta A(\boldsymbol{n})}=\Omega(\boldsymbol{n},\theta)

となります。これは軸 n\boldsymbol{n} のまわりの角 θ\theta の回転を表す行列そのものなので、YRY\in R です。以上で示せました。軸と回転角が n=a0/a0\boldsymbol{n}=\boldsymbol{a}_0/|\boldsymbol{a}_0|θ=a0\theta=|\boldsymbol{a}_0| と具体的に読み取れる点まで込めての結論です。

行列の性質から直接確認することもできます。XX が実なら級数の各項が実なので YY は実行列です。転置は級数の各項ごとに取れて

tY=t(eX)=etX=eX=(eX)1=Y1{}^{t}Y={}^{t}\bigl(e^{X}\bigr)=e^{{}^{t}X}=e^{-X}=\bigl(e^{X}\bigr)^{-1}=Y^{-1}

なので tYY=YtY=I{}^{t}YY=Y{}^{t}Y=I、つまり YY は実直交行列です。さらに deteX=eTrX\det e^{X}=e^{\operatorname{Tr}X} と交代行列の対角成分が 00 であること(TrX=0\operatorname{Tr}X=0)から detY=e0=1\det Y=e^{0}=1 です。設問(i)(d) で見たように直交性だけでは足りず、detY=+1\det Y=+1 まで言えて初めて回転だと結論できます。ここまでで YSO(3)Y\in SO(3) が分かり、SO(3)SO(3) の元が必ずある軸のまわりの回転であること(オイラーの定理)を使えば同じ結論になりますが、上の構成的な議論はその軸と角を明示的に与えている点で強い主張です。

第2問 エルミート多項式の母関数と微分方程式

Section titled “第2問 エルミート多項式の母関数と微分方程式”

n=0,1,2,n=0,1,2,\dots に対して

fn(x)=ex2dndxn(ex2)f_n(x)=e^{x^2}\frac{d^n}{dx^n}\bigl(e^{-x^2}\bigr)

と定義します。以下では g(x)=ex2g(x)=e^{-x^2} と書きます。定義は ex2fn(x)=g(n)(x)e^{-x^2}f_n(x)=g^{(n)}(x) と同値で、積分の設問ではこの形をそのまま使います。設問(iv) 以降では xx を複素変数 zz に置き換えますが、ez2e^{-z^2} が整関数なので fnf_n は複素平面全体で定義された多項式になります。

漸化式を先に用意します。定義から

fn+1(x)=ex2ddx(ex2fn(x))=ex2(ex2fn(x)2xex2fn(x))=fn(x)2xfn(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)

です。f0=1f_0=1 から順に

f1=2x,f2=4x22,f3=8x3+12x,f4=16x448x2+12f_1=-2x,\qquad f_2=4x^2-2,\qquad f_3=-8x^3+12x,\qquad f_4=16x^4-48x^2+12

が得られます。エルミート多項式 HnH_n とは fn(x)=(1)nHn(x)f_n(x)=(-1)^nH_n(x) という関係にあります。

もう1つ、以下で何度も使う減衰の事実を確認しておきます。g(k)(x)g^{(k)}(x) は(xx の多項式)×ex2\times e^{-x^2} の形なので、任意の m,k0m,k\ge 0 に対して xmg(k)(x)0x^mg^{(k)}(x)\to 0x|x|\to\infty)であり、この形の関数の R\mathbb{R} 上の積分は絶対収束します。以下の部分積分の境界項はすべてこれで消えます。

nn 次です。漸化式 fn+1=fn2xfnf_{n+1}=f_n'-2xf_n による帰納法で示します。f0=1f_0=100 次で最高次係数は a0=1a_0=1 です。fnf_nnn 次で最高次係数 an0a_n\neq 0 だとすると、右辺の fnf_n'n1n-1 次以下、2xfn-2xf_nn+1n+1 次で最高次係数 2an-2a_n です。次数が異なるので打ち消しは起こらず、fn+1f_{n+1}n+1n+1 次で an+1=2an0a_{n+1}=-2a_n\neq 0 となります。よって degfn=n\deg f_n=n で、最高次係数は an=(2)na_n=(-2)^n です。答えは xxnn 次多項式です。

ex2fn(x)=g(n)(x)e^{-x^2}f_n(x)=g^{(n)}(x) を使うと

xfn(x)ex2dx=xg(n)(x)dx\int_{-\infty}^{\infty}xf_n(x)e^{-x^2}dx=\int_{-\infty}^{\infty}x\,g^{(n)}(x)\,dx

です。n>1n>1 すなわち n2n\ge 2 のもとで部分積分すると

xg(n)(x)dx=[xg(n1)(x)]g(n1)(x)dx=[g(n2)(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

となります。1つ目の境界項は設定で確認した減衰から 00 です。2つ目では n11n-1\ge 1 なので g(n1)g^{(n-1)} の原始関数として g(n2)g^{(n-2)}n=2n=2 のときは gg 自身)が取れ、その境界値も 00 です。答えは 00 です。

n=1n=1 を除外している理由もはっきりします。n=1n=1 では g(0)dx=ex2dx=π\int g^{(0)}dx=\int e^{-x^2}dx=\sqrt{\pi} が残り、積分値は π0-\sqrt{\pi}\neq 0 です。n2n\ge 2 で消えるのは、重み ex2e^{-x^2} についての直交性から fnf_n11 次以下の多項式と直交することの言い換えでもあります。

同じ書き換えで

xnfn(x)ex2dx=xng(n)(x)dx\int_{-\infty}^{\infty}x^nf_n(x)e^{-x^2}dx=\int_{-\infty}^{\infty}x^n\,g^{(n)}(x)\,dx

です。部分積分を nn 回繰り返します。jj 回目の境界項は [dj1dxj1(xn)g(nj)(x)]\bigl[\frac{d^{j-1}}{dx^{j-1}}(x^n)\cdot g^{(n-j)}(x)\bigr]_{-\infty}^{\infty} の形で、設定の減衰からすべて 00 です。よって

xng(n)(x)dx=(1)ndndxn(xn)g(x)dx=(1)nn!ex2dx\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

となり、答えは

xnfn(x)ex2dx=(1)nn!π\int_{-\infty}^{\infty}x^nf_n(x)e^{-x^2}dx=(-1)^n\,n!\,\sqrt{\pi}

です。

小さい nn で検算します。n=0n=0 では ex2dx=π\int e^{-x^2}dx=\sqrt{\pi} で、式は (1)00!π=π(-1)^0 0!\sqrt{\pi}=\sqrt{\pi} を与えます。n=1n=1 では f1=2xf_1=-2x から (2x2)ex2dx=2π2=π\int(-2x^2)e^{-x^2}dx=-2\cdot\frac{\sqrt{\pi}}{2}=-\sqrt{\pi} で、式の値 (1)11!π(-1)^1\,1!\,\sqrt{\pi} と合います。n=2n=2 では f2=4x22f_2=4x^2-2 から x2(4x22)ex2dx=43π42π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} で、2!π2!\sqrt{\pi} と一致します。

eω2e^{-\omega^2} は整関数なので、与えられた表示は zz を中心とする半径 11 の円周に対するコーシーの積分公式です。まず nn 階微分の表示に持ち上げます。zz を固定し C={ω:ωz=1}C=\{\omega:|\omega-z|=1\}(反時計回り)と書きます。ζz<1/2|\zeta-z|<1/2 を満たす点 ζ\zeta に対して、ζ\zeta を中心とする半径 11 の円周は、eω2e^{-\omega^2} が正則な領域(複素平面全体)の内部で CC へ連続変形でき、どちらも ζ\zeta を正の向きに1回だけ囲みます。したがってコーシーの定理から経路を CC に取り替えられて

eζ2=12πiCeω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)

が成り立ちます。右辺の経路 CCζ\zeta に依らないので、ζ\zeta について微分できます。実際、ωC\omega\in Cζz<1/2|\zeta-z|<1/2 のとき ωζ1/2|\omega-\zeta|\ge 1/2 なので分母が 00 に近づかず、被積分関数は ζ\zeta について何回でも微分可能でその導関数は ωC\omega\in C について一様有界です。CC はコンパクトなので微分と積分は交換できます。ζ\zetann 回微分して ζ=z\zeta=z と置くと

dndznez2=n!2πiCeω2(ωz)n+1dω,つまりfn(z)n!=ez22πiCeω2(ωz)n+1dω\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

です(グルサの公式)。

これを級数に代入します。

n=0tnn!fn(z)=ez212πiCeω2n=0tn(ωz)n+1dω\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

順序交換は問題の指示で認められていますが、根拠も明確です。CC 上では ωz=1|\omega-z|=1 なので t/(ωz)=t<1\bigl|t/(\omega-z)\bigr|=|t|<1ω\omega に依らない定数で押さえられ、級数は CC 上一様絶対収束します。eω2e^{-\omega^2}CC 上有界なので、有界な連続関数を掛けても一様収束は保たれ、コンパクトな経路上での項別積分が正当化されます。

等比級数を足すと

n=0tn(ωz)n+1=1ωz11tω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=0tnn!fn(z)=ez212πiCeω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

となります。t<1|t|<1 より点 z+tz+t は円 CC の内部にあり、eω2e^{-\omega^2} は整関数なので、右辺の積分にコーシーの積分公式が使えて値は e(z+t)2e^{-(z+t)^2} です。よって

n=0tnn!fn(z)=ez2e(z+t)2=et22tz\sum_{n=0}^{\infty}\frac{t^n}{n!}f_n(z)=e^{z^2}e^{-(z+t)^2}=e^{-t^2-2tz}

が示せました。仮定 t<1|t|<1 は、等比級数の収束と z+tz+tCC の内部に入ることの両方に使われています。

低次で確認しておきます。右辺を tt で展開すると et22tz=12tz+(2z21)t2+e^{-t^2-2tz}=1-2tz+(2z^2-1)t^2+\cdots で、左辺の t0,t1,t2t^0,t^1,t^2 の係数は f0=1f_0=1f1=2zf_1=-2zf2/2!=(4z22)/2=2z21f_2/2!=(4z^2-2)/2=2z^2-1 です。一致します。

F(t,z)=et22tzF(t,z)=e^{-t^2-2tz} と置きます。与えられた式は直接確かめられます。F/z=2tF\partial F/\partial z=-2tF2F/z2=4t2F\partial^2F/\partial z^2=4t^2FF/t=(2t+2z)F\partial F/\partial t=-(2t+2z)F なので

2Fz22zFz=(4t2+4tz)F=2t(2t2z)F=2tFt\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}

です。

設問(iv) より t<1|t|<1F(t,z)=m0tmm!fm(z)F(t,z)=\sum_{m\ge 0}\frac{t^m}{m!}f_m(z) であり、収束冪級数の係数は一意なので

fn(z)=nFtnt=0f_n(z)=\frac{\partial^nF}{\partial t^n}\bigg|_{t=0}

が成り立ちます。そこで上の式の両辺を ttnn 回微分して t=0t=0 と置きます。FF(t,z)(t,z) の整関数なので偏微分の順序はすべて交換できます。

左辺は zz 微分と tt 微分を入れ替えて

ntn(2Fz22zFz)t=0=d2dz2fn(z)2zddzfn(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)

となります。右辺はライプニッツ則から

ntn(tFt)=tn+1Ftn+1+nnFtn\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=0t=0 で第1項が落ちて nfn(z)n f_n(z) が残り、全体で 2nfn(z)-2nf_n(z) です。両者を等号で結べば

d2dz2fn(z)2zddzfn(z)+2nfn(z)=0\frac{d^2}{dz^2}f_n(z)-2z\frac{d}{dz}f_n(z)+2n\,f_n(z)=0

が得られます。これが示すべき微分方程式で、答えは λ=2n\lambda=2n です。

λ\lambda の値は最高次の項だけを見ても確認できます。fn=anzn+f_n=a_nz^n+\cdotsan=(2)n0a_n=(-2)^n\neq 0)を代入すると、znz^n の係数は fnf_n'' から寄与がなく、2zfn-2zf_n' から 2nan-2na_nλfn\lambda f_n から λan\lambda a_n なので、(λ2n)an=0(\lambda-2n)a_n=0 すなわち λ=2n\lambda=2n が必要です。多項式全体で成り立つことは上の議論が保証します。具体例でも確かめられます。f2=4z22f_2=4z^2-2 では f2=8f_2''=82zf2=16z2-2zf_2'=-16z^24f2=16z284f_2=16z^2-8 の和が 00 です。f3=8z3+12zf_3=-8z^3+12z では f3=48zf_3''=-48z2zf3=48z324z-2zf_3'=48z^3-24z6f3=48z3+72z6f_3=-48z^3+72z の和が 00 です。λ=2n\lambda=2n はエルミートの微分方程式の固有値で、設定で述べた fn=(1)nHnf_n=(-1)^nH_n という対応と整合しています。

出典: 東京大学大学院理学系研究科 物理学専攻 平成20年度 修士課程 入学試験問題 数学。問題文は要約して引用しています。

この記事の誤りを報告する ・運営: 夢現技研合同会社料金プラン利用条件特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・本文の LLM への入力は自由です。コード例は MIT ライセンスです。