跳到主要内容

東北大学 理学研究科 天文学専攻 2025年度 物理(4)

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

問1

質量 mm、角振動数 ω\omega の一次元調和振動子を考える。

(a)

波動関数を ϕ(x)\phi(x)、エネルギーを EE として、時間に依存しない Schrödinger 方程式を書け。

(b)

無次元量

ξ=mωx,ϵ=2Eω\xi=\sqrt{\frac{m\omega}{\hbar}}x, \qquad \epsilon=\frac{2E}{\hbar\omega}

を用いて方程式を無次元化せよ。

(c)

ξ|\xi| が大きいとき、規格化可能な波動関数が

ϕ(ξ)exp(ξ22)\phi(\xi)\simeq\exp\left(-\frac{\xi^2}{2}\right)

と漸近することを示せ。

(d)

ϕ(ξ)=H(ξ)eξ2/2\phi(\xi)=H(\xi)e^{-\xi^2/2} とおき、HH の微分方程式を求めよ。

(e)

H(ξ)=n=0anξnH(\xi)=\sum_{n=0}^{\infty}a_n\xi^n として、ana_n の漸化式を求めよ。

(f)

規格化可能性から級数が有限次数で終わることを用いて、

En=ω(n+12)E_n=\hbar\omega\left(n+\frac12\right)

を示せ。

問2

zz 方向の一様な磁場 B=Bez\boldsymbol B=B\boldsymbol e_z 中の電子を考える。スピンは無視し、電子の質量を mm、電荷を e-e とする。

(a)

電磁場が、任意のスカラー関数 λ\lambda に対する

ϕ=ϕ+λt,A=Aλ\phi'=\phi+\frac{\partial\lambda}{\partial t}, \qquad \boldsymbol A'=\boldsymbol A-\nabla\lambda

で不変であることを示し、この性質の名称を答えよ。

(b)

B=Bez\boldsymbol B=B\boldsymbol e_z を与え、Ax=Az=0A_x=A_z=0 を満たす最も簡単な AyA_y を求めよ。

電子の Hamiltonian は

H^=12m[p^+eA(x^)]2\hat H=\frac1{2m} [\hat{\boldsymbol p}+e\boldsymbol A(\hat{\boldsymbol x})]^2

とする。

(c)

zz 方向の運動量 pzp_z が保存することを示せ。

(d)

xxyy 平面内では

H^=12m[p^x2+(p^y+eAy)2]\hat H=\frac1{2m}\left[\hat p_x^2+(\hat p_y+eA_y)^2\right]

となる。[p^y,H^]=0[\hat p_y,\hat H]=0 を示せ。

(e)

p^y\hat p_y の固有値を kyk_y とするとき、xx 方向の運動が調和振動子に対応することを示せ。

(f)

量子化されたエネルギー準位(Landau 準位)を求めよ。

(g)

基底状態と第 nn 励起状態のエネルギー差を

ΔEn=nX(B1 T)eV\Delta E_n=nX\left(\frac{B}{1\ \mathrm{T}}\right)\mathrm{eV}

と書き、XX を有効数字1桁で求めよ。必要なら

=6.6×1016 eVs,m=9.1×1031 kg,e=1.6×1019 C\hbar=6.6\times10^{-16}\ \mathrm{eV\,s},\quad m=9.1\times10^{-31}\ \mathrm{kg},\quad e=1.6\times10^{-19}\ \mathrm C

を用いてよい。

(h)

ある天体の X 線スペクトルに、ほぼ等間隔の吸収構造が観測された。図の矢印の概略位置は次の通りである。これらを電子の Landau 準位によるものとして磁場 BB を推定せよ。

吸収構造エネルギー
128 keV28\ \mathrm{keV}
255 keV55\ \mathrm{keV}
380 keV80\ \mathrm{keV}

Kai

問1

(a)

22md2ϕdx2+12mω2x2ϕ=Eϕ.-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} +\frac12m\omega^2x^2\phi=E\phi.

(b)

d2ϕdξ2+ξ2ϕ=ϵϕ.-\frac{d^2\phi}{d\xi^2}+\xi^2\phi=\epsilon\phi.

(c)

ϕ=eS\phi=e^{S} とおくと

S(S)2+ξ2=ϵ.-S''-(S')^2+\xi^2=\epsilon.

ξ|\xi|\to\infty では (S)2ξ2(S')^2\simeq\xi^2 だから S±ξ2/2S\simeq\pm\xi^2/2。規格化可能な方を選び、

ϕ(ξ)eξ2/2.\phi(\xi)\simeq e^{-\xi^2/2}.

(d)

代入すると

H2ξH+(ϵ1)H=0.H''-2\xi H'+(\epsilon-1)H=0.

(e)

係数比較により

(n+2)(n+1)an+2+(ϵ12n)an=0,(n+2)(n+1)a_{n+2}+(\epsilon-1-2n)a_n=0,

すなわち

an+2=2n+1ϵ(n+2)(n+1)an.a_{n+2}=\frac{2n+1-\epsilon}{(n+2)(n+1)}a_n.

(f)

級数が n=Nn=N で終わるには

2N+1ϵ=02N+1-\epsilon=0

が必要である。したがって

ϵN=2N+1,EN=ω2ϵN=ω(N+12).\epsilon_N=2N+1, \qquad E_N=\frac{\hbar\omega}{2}\epsilon_N =\hbar\omega\left(N+\frac12\right).

問2

(a)

E=ϕAt=ϕAt=E,B=×A=×A=B.\begin{aligned} \boldsymbol E' &=-\nabla\phi'-\frac{\partial\boldsymbol A'}{\partial t} =-\nabla\phi-\frac{\partial\boldsymbol A}{\partial t} =\boldsymbol E,\\ \boldsymbol B' &=\nabla\times\boldsymbol A' =\nabla\times\boldsymbol A =\boldsymbol B. \end{aligned}

したがって電磁場は不変である。この性質を gauge invariance(ゲージ不変性)という。

(b)

A=(0,Bx,0)\boxed{\boldsymbol A=(0,Bx,0)}

とすれば (×A)z=Ay/x=B(\nabla\times\boldsymbol A)_z=\partial A_y/\partial x=B となる。

(c)

A\boldsymbol AH^\hat Hzz に依存しないから

dp^zdt=1i[p^z,H^]=0.\frac{d\hat p_z}{dt} =\frac1{i\hbar}[\hat p_z,\hat H] =0.

よって pzp_z は保存する。

(d)

Ay=BxA_y=Bxyy に依存せず、[p^y,p^x]=[p^y,x^]=0[\hat p_y,\hat p_x]=[\hat p_y,\hat x]=0 である。したがって

[p^y,H^]=0.\boxed{[\hat p_y,\hat H]=0}.

(e)

p^y\hat p_y の固有値を kyk_y とすると

H^x=p^x22m+(ky+eBx)22m=p^x22m+12mωc2(x+x0)2,\begin{aligned} \hat H_x &=\frac{\hat p_x^2}{2m}+\frac{(k_y+eBx)^2}{2m}\\ &=\frac{\hat p_x^2}{2m} +\frac12m\omega_c^2(x+x_0)^2, \end{aligned}

ただし

ωc=eBm,x0=kyeB.\omega_c=\frac{eB}{m},\qquad x_0=\frac{k_y}{eB}.

よって中心が x0-x_0 の調和振動子である。

(f)

問1の結果より

En=ωc(n+12)=eBm(n+12)(n=0,1,2,).\boxed{E_n=\hbar\omega_c\left(n+\frac12\right) =\frac{\hbar eB}{m}\left(n+\frac12\right)} \qquad(n=0,1,2,\ldots).

(g)

ΔEn=neBm=n(em1T)(B1T).\Delta E_n=n\frac{\hbar eB}{m} =n\left(\frac{\hbar e}{m}\,1\mathrm T\right) \left(\frac{B}{1\mathrm T}\right).

したがって

X=(6.6×1016)(1.6×1019)9.1×1031=1.16×104.X=\frac{(6.6\times10^{-16})(1.6\times10^{-19})} {9.1\times10^{-31}} =1.16\times10^{-4}.

有効数字1桁で

X=1×104.\boxed{X=1\times10^{-4}}.

(h)

吸収構造の間隔は約 28 keV28\ \mathrm{keV} なので

B2.8×104 eV1.16×104 eV/T=2.4×108 T.B\simeq \frac{2.8\times10^4\ \mathrm{eV}} {1.16\times10^{-4}\ \mathrm{eV/T}} =2.4\times10^8\ \mathrm T.

したがって

B2.4×108 T=2.4×1012 G.\boxed{B\simeq2.4\times10^8\ \mathrm T =2.4\times10^{12}\ \mathrm G}.

Reference