Skip to content

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

Raw

This content is not available in your language yet.

この年度の問題冊子は数学2題と英語2題を合わせた1冊で、4問すべてが必答、試験時間は9時から11時までの120分です。この記事では数学の2題を扱います。第1問は2次の実対称行列の固有値・固有ベクトル・対角化という定型の流れですが、最後に幾何的な設問が置かれていて、行列を A=IγBA=I-\gamma BBB は鏡映)と見抜けるかどうかで解答の見通しが大きく変わります。第2問は3次元ラプラス方程式の変数分離で、動径方程式がオイラー型になること、球面調和関数の直交性で境界条件を係数に翻訳することという二つの手続きを、球殻の具体的な境界値問題まで通して実行させる構成です。どちらも計算量は軽く、独立性の根拠や係数が一意に決まる根拠をどこまで書けるかが点差になります。

問題分野主題
第1問線形代数鏡映を含む2次実対称行列の対角化と接線方向シアー
第2問微分方程式・特殊関数・フーリエ解析ラプラス方程式の動径方程式と球面調和関数展開

第1問 2次実対称行列の対角化と円の変形

Section titled “第1問 2次実対称行列の対角化と円の変形”

2行2列の実行列

A=(1γcos2ϕγsin2ϕγsin2ϕ1+γcos2ϕ),0<γ<1A=\begin{pmatrix} 1-\gamma\cos 2\phi & -\gamma\sin 2\phi \\ -\gamma\sin 2\phi & 1+\gamma\cos 2\phi \end{pmatrix}, \qquad 0<\gamma<1

について、固有値 λ±\lambda_\pmλ+>λ\lambda_+>\lambda_-)、規格化された固有ベクトル w±\vec w_\pmUAU1=diag(λ,λ+)UAU^{-1}=\mathrm{diag}(\lambda_-,\lambda_+) となる UU を求め、最後に AA が引き起こす図形の変形を調べます。ϕ\phi は実パラメータで、設問(iv) では単位円 CC 上の点 P=(cosϕ,sinϕ)\mathrm P=(\cos\phi,\sin\phi) の偏角として現れます。

以下で繰り返し使うので、分解を先に書いておきます。

A=IγB,B=(cos2ϕsin2ϕsin2ϕcos2ϕ)A=I-\gamma B,\qquad B=\begin{pmatrix}\cos 2\phi & \sin 2\phi \\ \sin 2\phi & -\cos 2\phi\end{pmatrix}

BB は原点を通り xx 軸と角 ϕ\phi をなす直線に関する鏡映の行列です。実際 B2=IB^2=ItrB=0\operatorname{tr}B=0detB=cos22ϕsin22ϕ=1\det B=-\cos^2 2\phi-\sin^2 2\phi=-1 で、BB は実対称かつ直交、固有値は +1+11-1 です。AA は実対称なので固有値は実で、異なる固有値に属する固有ベクトルは直交します。

トレースと行列式から特性方程式を作ります。

trA=2,detA=(1γcos2ϕ)(1+γcos2ϕ)γ2sin22ϕ=1γ2\operatorname{tr}A=2,\qquad \det A=(1-\gamma\cos 2\phi)(1+\gamma\cos 2\phi)-\gamma^2\sin^2 2\phi =1-\gamma^2

ϕ\phi 依存性が detA\det A で完全に消えるところがこの問題の仕掛けです。よって

det(AλI)=λ22λ+(1γ2)=(λ1)2γ2=0\det(A-\lambda I)=\lambda^2-2\lambda+(1-\gamma^2)=(\lambda-1)^2-\gamma^2=0

となり λ=1±γ\lambda=1\pm\gamma を得ます。0<γ<10<\gamma<1 より 1+γ>1γ>01+\gamma>1-\gamma>0 ですから、答えは

λ+=1+γ,λ=1γ\lambda_+=1+\gamma,\qquad \lambda_-=1-\gamma

です。どちらも正なので AA は正定値で、detA=λ+λ=1γ2>0\det A=\lambda_+\lambda_-=1-\gamma^2>0λ++λ=2=trA\lambda_++\lambda_-=2=\operatorname{tr}A と整合します。A=IγBA=I-\gamma BBB の固有値 ±1\pm 1 からも同じ結果が直ちに出ます。

固有空間を、固有値を代入した行列を因数分解する形で取り出します。半角の公式 1cos2ϕ=2sin2ϕ1-\cos 2\phi=2\sin^2\phi1+cos2ϕ=2cos2ϕ1+\cos 2\phi=2\cos^2\phisin2ϕ=2sinϕcosϕ\sin 2\phi=2\sin\phi\cos\phi を使うと

AλI=γ(1cos2ϕsin2ϕsin2ϕ1+cos2ϕ)=2γ(sinϕcosϕ)(sinϕcosϕ),Aλ+I=γ(1+cos2ϕsin2ϕsin2ϕ1cos2ϕ)=2γ(cosϕsinϕ)(cosϕsinϕ)\begin{aligned} A-\lambda_- I&=\gamma\begin{pmatrix}1-\cos 2\phi & -\sin 2\phi\\ -\sin 2\phi & 1+\cos 2\phi\end{pmatrix} =2\gamma\begin{pmatrix}\sin\phi\\ -\cos\phi\end{pmatrix}\begin{pmatrix}\sin\phi & -\cos\phi\end{pmatrix},\\[2pt] A-\lambda_+ I&=-\gamma\begin{pmatrix}1+\cos 2\phi & \sin 2\phi\\ \sin 2\phi & 1-\cos 2\phi\end{pmatrix} =-2\gamma\begin{pmatrix}\cos\phi\\ \sin\phi\end{pmatrix}\begin{pmatrix}\cos\phi & \sin\phi\end{pmatrix} \end{aligned}

と書けます。γ0\gamma\neq 0 なのでどちらも階数1で、核はちょうど1次元です。AλIA-\lambda_-I の核は (sinϕ,cosϕ)(\sin\phi,-\cos\phi) に直交するベクトル、Aλ+IA-\lambda_+I の核は (cosϕ,sinϕ)(\cos\phi,\sin\phi) に直交するベクトルなので、規格化された固有ベクトルは

w+=(sinϕcosϕ),w=(cosϕsinϕ)\vec w_+=\begin{pmatrix}-\sin\phi\\ \cos\phi\end{pmatrix},\qquad \vec w_-=\begin{pmatrix}\cos\phi\\ \sin\phi\end{pmatrix}

です(全体の符号は任意で、w±=1\lvert \vec w_\pm\rvert=1w+w=0\vec w_+\cdot\vec w_-=0)。

念のため直接代入して確かめます。加法定理 cos2ϕcosϕ+sin2ϕsinϕ=cosϕ\cos 2\phi\cos\phi+\sin 2\phi\sin\phi=\cos\phisin2ϕcosϕcos2ϕsinϕ=sinϕ\sin 2\phi\cos\phi-\cos 2\phi\sin\phi=\sin\phi を使うと

Aw=(cosϕγ(cos2ϕcosϕ+sin2ϕsinϕ)sinϕγ(sin2ϕcosϕcos2ϕsinϕ))=(1γ)w,Aw+=(sinϕγ(sin2ϕcosϕcos2ϕsinϕ)cosϕ+γ(cos2ϕcosϕ+sin2ϕsinϕ))=(1+γ)w+\begin{aligned} A\vec w_-&=\begin{pmatrix}\cos\phi-\gamma(\cos 2\phi\cos\phi+\sin 2\phi\sin\phi)\\ \sin\phi-\gamma(\sin 2\phi\cos\phi-\cos 2\phi\sin\phi)\end{pmatrix} =(1-\gamma)\vec w_-,\\[2pt] A\vec w_+&=\begin{pmatrix}-\sin\phi-\gamma(\sin 2\phi\cos\phi-\cos 2\phi\sin\phi)\\ \cos\phi+\gamma(\cos 2\phi\cos\phi+\sin 2\phi\sin\phi)\end{pmatrix} =(1+\gamma)\vec w_+ \end{aligned}

となり、確かに固有値は λ\lambda_\mp です。

幾何的には w\vec w_- は偏角 ϕ\phi の点における動径方向の単位ベクトル e^r\hat e_rw+\vec w_+ はそれに直交する方位角方向の単位ベクトル e^ϕ\hat e_\phi です。動径方向が縮む向き(固有値 1γ<11-\gamma<1)、方位角方向が伸びる向き(固有値 1+γ>11+\gamma>1)である、というこの対応が設問(iv) の答えを決めます。

UAU1=DUAU^{-1}=DD=diag(λ,λ+)D=\mathrm{diag}(\lambda_-,\lambda_+)UA=DUUA=DU と同値で、成分で書けば UU の第 iiuiu_iuiA=λiuiu_iA=\lambda_i u_i を満たす、つまり uiu_i は左固有ベクトルです。AA は対称なので左固有ベクトルは右固有ベクトルの転置に一致します。γ>0\gamma>0 より λ+λ\lambda_+\neq\lambda_- で固有空間はいずれも1次元ですから、第1行は wT\vec w_-^{\mathsf T}、第2行は w+T\vec w_+^{\mathsf T} のスカラー倍に限られます。行ベクトルを単位ベクトルに取れば

U=(cosϕsinϕsinϕcosϕ),U1=UT=(cosϕsinϕsinϕcosϕ)U=\begin{pmatrix}\cos\phi & \sin\phi\\ -\sin\phi & \cos\phi\end{pmatrix}, \qquad U^{-1}=U^{\mathsf T}=\begin{pmatrix}\cos\phi & -\sin\phi\\ \sin\phi & \cos\phi\end{pmatrix}

です。UTU=IU^{\mathsf T}U=IdetU=1\det U=1 を満たします(各行を独立に 00 でない定数倍しても対角化は成り立つので、UU はこの規格化のもとで符号を除いて一意です)。

これが確かに要求された形になることを、鏡映の表示から確認します。BB は角 ϕ\phi の直線に関する鏡映なので、回転行列 R(ϕ)=U1R(\phi)=U^{-1} を使って B=R(ϕ)diag(1,1)R(ϕ)1B=R(\phi)\,\mathrm{diag}(1,-1)\,R(\phi)^{-1} と書けます。したがって

UAU1=IγUBU1=Iγdiag(1,1)=(1γ001+γ)=(λ00λ+)UAU^{-1}=I-\gamma\,UBU^{-1}=I-\gamma\,\mathrm{diag}(1,-1)=\begin{pmatrix}1-\gamma & 0\\ 0 & 1+\gamma\end{pmatrix} =\begin{pmatrix}\lambda_- & 0\\ 0 & \lambda_+\end{pmatrix}

となり、λ\lambda_- が左上、λ+\lambda_+ が右下という指定どおりの並びです。

UU が表す変換について。U=R(ϕ)U=R(-\phi) なので、ベクトルに作用させる立場では原点まわりの角 ϕ-\phi の回転(時計まわりに ϕ\phi だけ回す操作)です。座標系を取り直す立場で読むほうが意味がはっきりします。(Uv)1=ve^r(U\vec v)_1=\vec v\cdot\hat e_r(Uv)2=ve^ϕ(U\vec v)_2=\vec v\cdot\hat e_\phi ですから、UU は標準基底での成分を、点 P\mathrm P における動径方向 e^r\hat e_r と方位角方向 e^ϕ\hat e_\phi を軸とする正規直交系での成分に書き換える操作、すなわち座標軸を +ϕ+\phi だけ回転させて極座標系に合わせる操作です。detU=1\det U=1 かつ直交行列なので、鏡映も伸縮も含まない純粋な回転であり、長さと角度を保ちます。

単位円 CC 上の8個の点

Pk=(cosϕk, sinϕk),ϕk=kπ4(k=0,1,,7)\mathrm P_k=(\cos\phi_k,\ \sin\phi_k),\qquad \phi_k=\frac{k\pi}{4}\quad(k=0,1,\dots,7)

をそれぞれ中心とする半径 ρ\rho の円(図の点線の円、ρ\rhoCC の半径より十分小さい共通の値)が描かれています。各点では ϕ=ϕk\phi=\phi_k を代入した AA が定義されていて、写像は「Pk\mathrm P_k を始点とする変位ベクトル w\vec wAwA\vec w に置き換える」というものです。したがって w=0\vec w=\vec 0 は動かず、円の中心 Pk\mathrm P_k は写像の前後で同じ位置に留まります。

各点で e^r=w\hat e_r=\vec w_-e^ϕ=w+\hat e_\phi=\vec w_+ を基底に取り、w=ξe^r+ηe^ϕ\vec w=\xi\,\hat e_r+\eta\,\hat e_\phi と分解します。点線の円は ξ2+η2=ρ2\xi^2+\eta^2=\rho^2 です。設問(ii) より

Aw=(1γ)ξe^r+(1+γ)ηe^ϕA\vec w=(1-\gamma)\,\xi\,\hat e_r+(1+\gamma)\,\eta\,\hat e_\phi

なので、像の成分を Aw=ξe^r+ηe^ϕA\vec w=\xi'\hat e_r+\eta'\hat e_\phi と書けば ξ=(1γ)ξ\xi'=(1-\gamma)\xiη=(1+γ)η\eta'=(1+\gamma)\eta です。0<γ<10<\gamma<1 よりどちらの係数も 00 でないので逆に解けて、ξ=ξ/(1γ)\xi=\xi'/(1-\gamma)η=η/(1+γ)\eta=\eta'/(1+\gamma) を円の式に代入すると

(ξ(1γ)ρ)2+(η(1+γ)ρ)2=1\left(\frac{\xi'}{(1-\gamma)\rho}\right)^2+\left(\frac{\eta'}{(1+\gamma)\rho}\right)^2=1

を得ます。これは中心 Pk\mathrm P_k、動径方向の半軸 (1γ)ρ(1-\gamma)\rho、方位角方向の半軸 (1+γ)ρ(1+\gamma)\rho の楕円です。detA=1γ20\det A=1-\gamma^2\neq 0 なので像はつぶれず、確かに楕円になります。

答えをまとめます。8つの円はいずれも中心を動かさないまま楕円に変形し、長軸は方位角方向(OPk\mathrm{OP}_k に垂直、すなわち CC の接線方向)を向いて長さ 2(1+γ)ρ2(1+\gamma)\rho、短軸は動径方向(OPk\mathrm{OP}_k 方向)を向いて長さ 2(1γ)ρ2(1-\gamma)\rho になります。短軸を長軸で割った軸比は (1γ)/(1+γ)<1(1-\gamma)/(1+\gamma)<1kk に依らず共通、面積は πρ2detA=πρ2(1γ2)\pi\rho^2\det A=\pi\rho^2(1-\gamma^2) と元の円よりわずかに小さくなります。伸びの向きは場所ごとに回転するので、たとえば ϕ=0\phi=0(右端)と ϕ=π\phi=\pi(左端)の円は縦長に、ϕ=π/2\phi=\pi/2(上端)と ϕ=3π/2\phi=3\pi/2(下端)の円は横長に、ϕ=π/4\phi=\pi/4 の円は (1,1)/2(-1,1)/\sqrt 2 方向に長い楕円になります。全体としては、円 CC に沿ってぐるりと接線方向に引き伸ばされた8つの楕円が並ぶ、8回対称の模様です。

CO

図は γ=1/2\gamma=1/2ρ=0.3\rho=0.3 の場合で、点線が元の円、実線が像の楕円です。極限の様子も整合します。γ0\gamma\to 0 では AIA\to I で楕円は円に戻り、γ1\gamma\to 1 では動径方向の半軸が 00 に近づいて楕円は接線方向の線分に潰れていきます。

なお AA は、収束 κ=0\kappa=0、シアー (γ1,γ2)=(γcos2ϕ, γsin2ϕ)(\gamma_1,\gamma_2)=(\gamma\cos 2\phi,\ \gamma\sin 2\phi) の弱重力レンズの変形行列そのもので、この設問の模様は点状質量のまわりに現れる接線方向のシアーパターンに対応します。

第2問 ラプラス方程式の動径方程式と球面調和関数展開

Section titled “第2問 ラプラス方程式の動径方程式と球面調和関数展開”

3次元ラプラス方程式を極座標 (r,θ,φ)(r,\theta,\varphi) で書いた形

(1r2r2rL^2r2)U(r,θ,φ)=0,L^2(1sinθθsinθθ+1sin2θ2φ2)\left(\frac{1}{r}\frac{\partial^2}{\partial r^2}r-\frac{\hat L^2}{r^2}\right)U(r,\theta,\varphi)=0, \qquad \hat L^2\equiv-\left(\frac{1}{\sin\theta}\frac{\partial}{\partial\theta}\sin\theta\frac{\partial}{\partial\theta} +\frac{1}{\sin^2\theta}\frac{\partial^2}{\partial\varphi^2}\right)

が与えられます。左の第1項は rUrU2/r2\partial^2/\partial r^2 を作用させてから rr で割る、という意味です。球面調和関数 Ym(θ,φ)Y_{\ell m}(\theta,\varphi)L^2Ym=(+1)Ym\hat L^2Y_{\ell m}=\ell(\ell+1)Y_{\ell m}=0,1,2,\ell=0,1,2,\dotsm=,,m=-\ell,\dots,\ell)を満たし、球面上で正規直交系をなします。

設問(iii)(iv) では、b>a>0b>a>0 として球殻 arba\le r\le b

U(a,θ,φ)=0,U(b,θ,φ)=sinθcosφU(a,\theta,\varphi)=0,\qquad U(b,\theta,\varphi)=\sin\theta\cos\varphi

というディリクレ境界条件を課します。L^2\hat L^2rr を含まず、動径部分の演算子は角度変数を含まないので、変数分離がそのまま通ります。

U=R(r)Ym(θ,φ)U=R_\ell(r)Y_{\ell m}(\theta,\varphi) を代入します。動径演算子は YmY_{\ell m} を素通りし、L^2\hat L^2 は固有値 (+1)\ell(\ell+1) を出すので

[1rd2dr2(rR(r))(+1)r2R(r)]Ym(θ,φ)=0\left[\frac{1}{r}\frac{d^2}{dr^2}\bigl(rR_\ell(r)\bigr)-\frac{\ell(\ell+1)}{r^2}R_\ell(r)\right]Y_{\ell m}(\theta,\varphi)=0

となります。YmY_{\ell m} は恒等的に 00 ではなく、Ym0Y_{\ell m}\neq 0 となる (θ,φ)(\theta,\varphi) は球面上で稠密ですから、角括弧の中が r>0r>000 でなければなりません。よって求める微分方程式は

1rd2dr2(rR)(+1)r2R=0(r>0)\frac{1}{r}\frac{d^2}{dr^2}\bigl(rR_\ell\bigr)-\frac{\ell(\ell+1)}{r^2}R_\ell=0 \qquad(r>0)

です。(rR)=rR+2R(rR_\ell)''=rR_\ell''+2R_\ell' を使って展開すれば

R+2rR(+1)r2R=0,すなわちr2R+2rR(+1)R=0R_\ell''+\frac{2}{r}R_\ell'-\frac{\ell(\ell+1)}{r^2}R_\ell=0, \qquad\text{すなわち}\qquad r^2R_\ell''+2rR_\ell'-\ell(\ell+1)R_\ell=0

と書けます。1r2ddr(r2dRdr)=R+2rR\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dR_\ell}{dr}\right)=R_\ell''+\frac{2}{r}R_\ell' なので、これは見慣れた動径方程式に一致します。方程式に mm が現れないこと、つまり同じ \ell に属する 2+12\ell+1 個の mm が同じ動径方程式を共有することが、解を RR_\ell と書いてよい理由です。

得られた方程式は r2R+2rR(+1)R=0r^2R''+2rR'-\ell(\ell+1)R=0 という同次のオイラー型方程式なので、R=rsR=r^s を試します。

r2s(s1)rs2+2rsrs1(+1)rs=[s2+s(+1)]rsr^2s(s-1)r^{s-2}+2rsr^{s-1}-\ell(\ell+1)r^s=\bigl[s^2+s-\ell(\ell+1)\bigr]r^s

より指数方程式は s2+s(+1)=0s^2+s-\ell(\ell+1)=0 です。判別式は 1+4(+1)=(2+1)21+4\ell(\ell+1)=(2\ell+1)^2 と完全平方になり

s=1±(2+1)2=  または  (+1)s=\frac{-1\pm(2\ell+1)}{2}=\ell\ \ \text{または}\ \ -(\ell+1)

を得ます。二つの根が一致するのは 2+1=02\ell+1=0 のときだけで、=0,1,2,\ell=0,1,2,\dots ではあり得ません。したがって重根の場合(対数を含む解)は生じず、すべての \ell に対して

A(r)=r,B(r)=r(+1)A_\ell(r)=r^{\ell},\qquad B_\ell(r)=r^{-(\ell+1)}

が答えです。=0\ell=0 なら A0=1A_0=1B0=1/rB_0=1/r です。

独立性を確認します。ロンスキアンは

W=ABAB=r((+1)r2)r1r1=2+1r2W=A_\ell B_\ell'-A_\ell'B_\ell =r^\ell\bigl(-(\ell+1)r^{-\ell-2}\bigr)-\ell r^{\ell-1}r^{-\ell-1} =-\frac{2\ell+1}{r^2}

で、r>0r>0 では決して 00 になりません。方程式を正規形 R+(2/r)R(+1)r2R=0R''+(2/r)R'-\ell(\ell+1)r^{-2}R=0 に直すと係数は r>0r>0 で連続なので、解空間は2次元です。ロンスキアンが消えない2解はその基底をなし、{A,B}\{A_\ell,B_\ell\} が基本解系になります。

sinθcosφ=12sinθ(eiφ+eiφ)\sin\theta\cos\varphi=\dfrac{1}{2}\sin\theta\left(e^{i\varphi}+e^{-i\varphi}\right) と書き、与えられた

Y1±1(θ,φ)=38πsinθe±iφY_{1\,\pm1}(\theta,\varphi)=\mp\sqrt{\frac{3}{8\pi}}\sin\theta\,e^{\pm i\varphi}

を逆に解きます。sinθeiφ=8π/3Y11\sin\theta\,e^{i\varphi}=-\sqrt{8\pi/3}\,Y_{11}sinθeiφ=8π/3Y11\sin\theta\,e^{-i\varphi}=\sqrt{8\pi/3}\,Y_{1\,-1} ですから

sinθcosφ=128π3(Y11Y11)=2π3(Y11(θ,φ)Y11(θ,φ))\sin\theta\cos\varphi=\frac{1}{2}\sqrt{\frac{8\pi}{3}}\bigl(Y_{1\,-1}-Y_{11}\bigr) =\sqrt{\frac{2\pi}{3}}\,\bigl(Y_{1\,-1}(\theta,\varphi)-Y_{11}(\theta,\varphi)\bigr)

が答えです。=1\ell=1 の成分だけで書けるのは、sinθcosφ=x/r\sin\theta\cos\varphi=x/r が単位球面に制限した1次の調和多項式だからです。

二通りに検算します。第一に実数性です。上の規約では Y11=Y11Y_{1\,-1}^{*}=-Y_{11}Y11=Y11Y_{11}^{*}=-Y_{1\,-1} なので (Y11Y11)=Y11Y11(Y_{1\,-1}-Y_{11})^{*}=Y_{1\,-1}-Y_{11} となり、右辺は実です。第二にノルムです。左辺を直接積分すると

dΩsin2θcos2φ=0πsin3θdθ02πcos2φdφ=43π=4π3\int d\Omega\,\sin^2\theta\cos^2\varphi =\int_0^\pi\sin^3\theta\,d\theta\int_0^{2\pi}\cos^2\varphi\,d\varphi =\frac{4}{3}\cdot\pi=\frac{4\pi}{3}

で、右辺は直交規格性から 2π3(1+1)=4π3\frac{2\pi}{3}(1+1)=\frac{4\pi}{3} です。一致します。

球殻 arba\le r\le b では rr^\ellr1r^{-\ell-1} も正則なので、両方を残した

U(r,θ,φ)==0m=[αmr+βmr1]Ym(θ,φ)U(r,\theta,\varphi)=\sum_{\ell=0}^{\infty}\sum_{m=-\ell}^{\ell} \Bigl[\alpha_{\ell m}r^{\ell}+\beta_{\ell m}r^{-\ell-1}\Bigr]Y_{\ell m}(\theta,\varphi)

を使います。境界条件の各式に YmY_{\ell'm'}^{*} を掛けて球面全体で積分し、直交条件 dΩYmYm=δδmm\int d\Omega\,Y_{\ell'm'}^{*}Y_{\ell m}=\delta_{\ell'\ell}\delta_{m'm} を使うと、各 (,m)(\ell,m) ごとに独立な2元連立1次方程式

αma+βma1=0,αmb+βmb1=cmdΩYm(θ,φ)sinθcosφ\begin{aligned} \alpha_{\ell m}a^{\ell}+\beta_{\ell m}a^{-\ell-1}&=0,\\ \alpha_{\ell m}b^{\ell}+\beta_{\ell m}b^{-\ell-1}&=c_{\ell m}\equiv\int d\Omega\,Y_{\ell m}^{*}(\theta,\varphi)\sin\theta\cos\varphi \end{aligned}

が得られます。設問(iii) の表示と直交条件から右辺は

c11=2π3,c11=2π3,cm=0  ((,m)(1,±1))c_{1\,-1}=\sqrt{\frac{2\pi}{3}},\qquad c_{11}=-\sqrt{\frac{2\pi}{3}},\qquad c_{\ell m}=0\ \ \bigl((\ell,m)\neq(1,\pm1)\bigr)

です。係数行列の行列式は

aa1bb1=ab1a1b=(ab)1(a2+1b2+1)\begin{vmatrix} a^{\ell} & a^{-\ell-1}\\ b^{\ell} & b^{-\ell-1}\end{vmatrix} =a^{\ell}b^{-\ell-1}-a^{-\ell-1}b^{\ell} =(ab)^{-\ell-1}\bigl(a^{2\ell+1}-b^{2\ell+1}\bigr)

で、0<a<b0<a<b かつ 2+112\ell+1\ge 1 より a2+1b2+1a^{2\ell+1}\neq b^{2\ell+1} なので、どの (,m)(\ell,m) でもこれは 00 になりません。したがって解は一意で、cm=0c_{\ell m}=0 となる (,m)(\ell,m) については αm=βm=0\alpha_{\ell m}=\beta_{\ell m}=0 が唯一の解です。=1\ell=1m=±1m=\pm1 だけが残ります。

=1\ell=1 では第1式から β=αa3\beta=-\alpha a^{3}、これを第2式に入れて

α(ba3b2)=cα=b2b3a3c,β=a3b2b3a3c\alpha\left(b-\frac{a^3}{b^2}\right)=c \quad\Longrightarrow\quad \alpha=\frac{b^{2}}{b^{3}-a^{3}}\,c,\qquad \beta=-\frac{a^{3}b^{2}}{b^{3}-a^{3}}\,c

です。c1±1c_{1\,\pm1} を入れると、00 でない係数は次の4つだけになります。

α11=2π3b2b3a3,β11=+2π3a3b2b3a3,α11=+2π3b2b3a3,β11=2π3a3b2b3a3\begin{aligned} \alpha_{11}&=-\sqrt{\frac{2\pi}{3}}\,\frac{b^{2}}{b^{3}-a^{3}}, & \beta_{11}&=+\sqrt{\frac{2\pi}{3}}\,\frac{a^{3}b^{2}}{b^{3}-a^{3}},\\[2pt] \alpha_{1\,-1}&=+\sqrt{\frac{2\pi}{3}}\,\frac{b^{2}}{b^{3}-a^{3}}, & \beta_{1\,-1}&=-\sqrt{\frac{2\pi}{3}}\,\frac{a^{3}b^{2}}{b^{3}-a^{3}} \end{aligned}

他はすべて αm=βm=0\alpha_{\ell m}=\beta_{\ell m}=0 です。b>ab>a なので分母 b3a3b^3-a^3 は正で、係数は有限です。

これらを展開に戻し、設問(iii) の関係を逆向きに使うと解は

U(r,θ,φ)=b2b3a3(ra3r2)sinθcosφU(r,\theta,\varphi)=\frac{b^{2}}{b^{3}-a^{3}}\left(r-\frac{a^{3}}{r^{2}}\right)\sin\theta\cos\varphi

とまとまります。検算します。r=ar=a では括弧が aa3/a2=0a-a^3/a^2=0 なので U=0U=0r=br=b では b2b3a3b3a3b2=1\frac{b^2}{b^3-a^3}\cdot\frac{b^3-a^3}{b^2}=1 より U=sinθcosφU=\sin\theta\cos\varphi で、境界条件を両方満たします。また rsinθcosφ=xr\sin\theta\cos\varphi=xr2sinθcosφ=x/r3r^{-2}\sin\theta\cos\varphi=x/r^{3} はどちらも r>0r>0 で調和関数なので、UU は球殻内でラプラス方程式を満たします。有界領域のディリクレ問題の解は一意なので、これが求める解です。a0a\to 0 の極限では Ursinθcosφ=xU\to r\sin\theta\cos\varphi=x となり、内側の球面が消えて外側の境界値をそのまま延ばした調和関数になる、という素直な振る舞いを示します。

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

Report an error in this article ・Operated by: Mugen Giken LLCPricingTermsLegal notice

© 2026 夢現技研合同会社 ・Feeding the text to an LLM is welcome. Code samples are MIT licensed.