コンテンツにスキップ

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

対象読者: 学部の物理(力学・解析力学、電磁気学、量子力学、熱・統計力学)を一通り学び、東大物理学専攻の院試問題を自分で解いた、または解こうとしている人。問題文は載せていないので、出典の問題を手元に置いて読みます。

生 Markdown

4時間で4問という構成の試験です。第1問は Aharonov–Casher のゼロモード数え上げを、演算子代数から正則関数の議論まで誘導つきで辿らせる問題で、最後は磁束量子数の整数部分という位相的な答えに落ちます。第2問と第3問は基礎的な統計力学と相対論ですが、第3問後半の4元運動量保存は微小量の扱いを丁寧にやらないと合いません。選択問題の3問はいずれも実験に近く、桁の評価とグラフの読み取りが得点を分けます。

問題分野主題
第1問量子力学・複素解析磁場中の2次元電子とゼロエネルギー束縛状態
第2問統計力学調和振動子系・1次元連成振動子系の内部エネルギーと比熱
第3問力学・相対論相対論的ロケットの質量減少
第4問電磁気学・波動プラズマ中の電磁波と干渉計による電子密度計測
第5問原子核高エネルギー電子弾性散乱と原子核の荷電半径
第6問電気回路フィルタ回路・整流回路とダイオードの非線形性

第1問から第3問が必答で、第4問から第6問のうち1問を選択する形式ですが、ここでは全6問の解答を載せます。

第1問 磁場中の2次元電子とゼロエネルギー束縛状態

Section titled “第1問 磁場中の2次元電子とゼロエネルギー束縛状態”

xyxy 平面内を動く電子を考えます。zz 方向の磁束密度は半径 RR の円盤内で一定 b>0b>0、外では 00 です。

Bz(x,y)={b>0,(r<R),0,(r≥R),r=x2+y2B_z(x,y)=\begin{cases} b>0, & (r<R),\\ 0, & (r\ge R),\end{cases}\qquad r=\sqrt{x^2+y^2}

波動関数は 2 成分 Ψ=(ψ↑,ψ↓)T\Psi=(\psi_\uparrow,\psi_\downarrow)^{\mathrm T} で、ハミルトニアン(エネルギーを表す演算子)は

H=px 2+py 22m12×2+ℏeBz(x,y)2mσz,px=−iℏ∂∂x+eAx,py=−iℏ∂∂y+eAyH=\frac{p_x^{\,2}+p_y^{\,2}}{2m}\mathbf{1}_{2\times2}+\frac{\hbar e B_z(x,y)}{2m}\sigma_z,\qquad p_x=-i\hbar\frac{\partial}{\partial x}+eA_x,\quad p_y=-i\hbar\frac{\partial}{\partial y}+eA_y

です。−e<0-e<0 が電子の電荷、Bz=∂xAy−∂yAxB_z=\partial_x A_y-\partial_y A_x。目標は E=0E=0 の束縛状態の個数を数えることです。

パウリ行列を直接掛けると

σxσy=(i00−i)=iσz,σyσx=−iσz\sigma_x\sigma_y=\begin{pmatrix}i&0\\0&-i\end{pmatrix}=i\sigma_z,\qquad \sigma_y\sigma_x=-i\sigma_z

なので [σx,σy]=2iσz[\sigma_x,\sigma_y]=2i\sigma_z です。よって

[Sx,Sy]=ℏ24[σx,σy]=ℏ24⋅2iσz=iℏ⋅ℏ2σz=iℏSz.[S_x,S_y]=\frac{\hbar^2}{4}[\sigma_x,\sigma_y]=\frac{\hbar^2}{4}\cdot 2i\sigma_z=i\hbar\cdot\frac{\hbar}{2}\sigma_z=i\hbar S_z .

[H,Sz][H,S_z] については、HH の第1項はスピン空間で単位行列に比例するので σz\sigma_z と可換、第2項は Bz(x,y)B_z(x,y) という座標の関数(スピン空間では数)に σz\sigma_z が掛かったものなので、やはり σz\sigma_z と可換です。px,pyp_x,p_y に含まれる Ax,AyA_x,A_y も座標の関数でスピンには作用しません。したがって

[H,Sz]=0.[H,S_z]=0 .

答えは [Sx,Sy]=iℏSz[S_x,S_y]=i\hbar S_z、[H,Sz]=0[H,S_z]=0 です。SzS_z は保存量で、以下スピン成分ごとに議論できます。

D=pxσx+pyσyD=p_x\sigma_x+p_y\sigma_y とします。px,pyp_x,p_y はスピンに作用しないので σz\sigma_z と可換であり、σxσz=−σzσx\sigma_x\sigma_z=-\sigma_z\sigma_x、σyσz=−σzσy\sigma_y\sigma_z=-\sigma_z\sigma_y より

Dσz=pxσxσz+pyσyσz=−σz(pxσx+pyσy)=−σzD.D\sigma_z=p_x\sigma_x\sigma_z+p_y\sigma_y\sigma_z=-\sigma_z(p_x\sigma_x+p_y\sigma_y)=-\sigma_z D .

次に D2D^2 を計算します。σx2=σy2=1\sigma_x^2=\sigma_y^2=\mathbf{1}、σxσy=iσz\sigma_x\sigma_y=i\sigma_z、σyσx=−iσz\sigma_y\sigma_x=-i\sigma_z を使うと

D2=px 2σx2+py 2σy2+pxpy σxσy+pypx σyσx=(px 2+py 2)12×2+i [px,py] σz.\begin{aligned} D^2&=p_x^{\,2}\sigma_x^2+p_y^{\,2}\sigma_y^2+p_xp_y\,\sigma_x\sigma_y+p_yp_x\,\sigma_y\sigma_x\\ &=(p_x^{\,2}+p_y^{\,2})\mathbf{1}_{2\times2}+i\,[p_x,p_y]\,\sigma_z . \end{aligned}

交換子は、任意の関数への作用を見れば

[px,py]=[−iℏ∂x, eAy]+[eAx, −iℏ∂y]=−iℏe(∂Ay∂x−∂Ax∂y)=−iℏeBz[p_x,p_y]=[-i\hbar\partial_x,\,eA_y]+[eA_x,\,-i\hbar\partial_y] =-i\hbar e\left(\frac{\partial A_y}{\partial x}-\frac{\partial A_x}{\partial y}\right)=-i\hbar e B_z

です。これを代入すると i[px,py]σz=ℏeBzσzi[p_x,p_y]\sigma_z=\hbar e B_z\sigma_z となり

D2=(px 2+py 2)12×2+ℏeBzσz=2mH,D^2=(p_x^{\,2}+p_y^{\,2})\mathbf{1}_{2\times2}+\hbar e B_z\sigma_z=2mH,

すなわち H=D2/(2m)H=D^2/(2m) が成り立ちます。ゼーマン項は運動量の非可換性から自動的に出てくる、というのがこの形の意味です。

∣Φn⟩=D∣Ψn⟩|\Phi_n\rangle=D|\Psi_n\rangle で、DD はエルミートなので

⟨Φn∣Φn⟩=⟨Ψn∣D†D∣Ψn⟩=⟨Ψn∣D2∣Ψn⟩=2m⟨Ψn∣H∣Ψn⟩=2mEn⟨Ψn∣Ψn⟩.\langle\Phi_n|\Phi_n\rangle=\langle\Psi_n|D^\dagger D|\Psi_n\rangle=\langle\Psi_n|D^2|\Psi_n\rangle =2m\langle\Psi_n|H|\Psi_n\rangle=2mE_n\langle\Psi_n|\Psi_n\rangle .

3番目の等号で式 (9)、最後で H∣Ψn⟩=En∣Ψn⟩H|\Psi_n\rangle=E_n|\Psi_n\rangle を使いました。

設問3の等式を EnE_n について解くと

En=⟨Φn∣Φn⟩2m⟨Ψn∣Ψn⟩E_n=\frac{\langle\Phi_n|\Phi_n\rangle}{2m\langle\Psi_n|\Psi_n\rangle}

です。内積は正定値なので分子は ≥0\ge 0、分母は >0>0(∣Ψn⟩|\Psi_n\rangle は規格化可能な非零状態)。したがって En≥0E_n\ge0 であり、仮定 En≠0E_n\ne0 と合わせて En>0E_n>0 が従います。

次に ∣Ψn⟩|\Psi_n\rangle が H∣Ψn⟩=En∣Ψn⟩H|\Psi_n\rangle=E_n|\Psi_n\rangle(En>0E_n>0)かつ σz∣Ψn⟩=+∣Ψn⟩\sigma_z|\Psi_n\rangle=+|\Psi_n\rangle を満たすとします。まず ∣Φn⟩=D∣Ψn⟩|\Phi_n\rangle=D|\Psi_n\rangle は零ベクトルではありません。実際、設問3より ⟨Φn∣Φn⟩=2mEn⟨Ψn∣Ψn⟩>0\langle\Phi_n|\Phi_n\rangle=2mE_n\langle\Psi_n|\Psi_n\rangle>0 です。式 (9) から H=D2/(2m)H=D^2/(2m) は DD と可換なので

H∣Φn⟩=D22mD∣Ψn⟩=D D22m∣Ψn⟩=D H∣Ψn⟩=En∣Φn⟩,H|\Phi_n\rangle=\frac{D^2}{2m}D|\Psi_n\rangle=D\,\frac{D^2}{2m}|\Psi_n\rangle=D\,H|\Psi_n\rangle=E_n|\Phi_n\rangle,

エネルギーは変わりません。一方、式 (8) の反交換関係から

σz∣Φn⟩=σzD∣Ψn⟩=−Dσz∣Ψn⟩=−D∣Ψn⟩=−∣Φn⟩\sigma_z|\Phi_n\rangle=\sigma_z D|\Psi_n\rangle=-D\sigma_z|\Psi_n\rangle=-D|\Psi_n\rangle=-|\Phi_n\rangle

となり、σz\sigma_z の固有値は −1-1 です。よって ∣Φn⟩|\Phi_n\rangle は同じエネルギー EnE_n で σz=−1\sigma_z=-1 の同時固有状態です。En>0E_n>0 の固有値はスピン上向きと下向きが必ず対になって現れ、DD がその間を移す、というのがこの結果の内容です。

E=0E=0 のとき設問3の等式は ⟨Φ∣Φ⟩=0\langle\Phi|\Phi\rangle=0 を与えます。内積が正定値なのでベクトル自身が零、すなわち

D∣Ψ⟩=0.D|\Psi\rangle=0 .

DD を行列で書くと

D=pxσx+pyσy=(0px−ipypx+ipy0)D=p_x\sigma_x+p_y\sigma_y=\begin{pmatrix}0& p_x-ip_y\\ p_x+ip_y&0\end{pmatrix}

なので、DΨ=0D\Psi=0 は 2 成分について分離した 2 本の式

(px+ipy)ψ↑=0,(px−ipy)ψ↓=0(p_x+ip_y)\psi_\uparrow=0,\qquad (p_x-ip_y)\psi_\downarrow=0

になります。ここで Ax=−∂yρA_x=-\partial_y\rho、Ay=+∂xρA_y=+\partial_x\rho とおくと Bz=∂xAy−∂yAx=∇2ρB_z=\partial_xA_y-\partial_yA_x=\nabla^2\rho で、また Ax+iAy=i(∂x+i∂y)ρA_x+iA_y=i(\partial_x+i\partial_y)\rho、Ax−iAy=−i(∂x−i∂y)ρA_x-iA_y=-i(\partial_x-i\partial_y)\rho です。したがって

px+ipy=−iℏ(∂x+i∂y)+e(Ax+iAy)=−iℏ[(∂x+i∂y)−eℏ(∂x+i∂y)ρ],px−ipy=−iℏ(∂x−i∂y)+e(Ax−iAy)=−iℏ[(∂x−i∂y)+eℏ(∂x−i∂y)ρ].\begin{aligned} p_x+ip_y&=-i\hbar(\partial_x+i\partial_y)+e(A_x+iA_y)=-i\hbar\left[(\partial_x+i\partial_y)-\frac{e}{\hbar}(\partial_x+i\partial_y)\rho\right],\\ p_x-ip_y&=-i\hbar(\partial_x-i\partial_y)+e(A_x-iA_y)=-i\hbar\left[(\partial_x-i\partial_y)+\frac{e}{\hbar}(\partial_x-i\partial_y)\rho\right]. \end{aligned}

そこで ψ↑=f↑exp⁡[(e/ℏ)ρ]\psi_\uparrow=f_\uparrow\exp[(e/\hbar)\rho]、ψ↓=f↓exp⁡[−(e/ℏ)ρ]\psi_\downarrow=f_\downarrow\exp[-(e/\hbar)\rho] と書くと、ρ\rho の微分から出る項がちょうど打ち消して

(px+ipy)ψ↑=−iℏ[(∂x+i∂y)f↑]eeℏρ,(px−ipy)ψ↓=−iℏ[(∂x−i∂y)f↓]e−eℏρ(p_x+ip_y)\psi_\uparrow=-i\hbar\left[(\partial_x+i\partial_y)f_\uparrow\right]e^{\frac{e}{\hbar}\rho},\qquad (p_x-ip_y)\psi_\downarrow=-i\hbar\left[(\partial_x-i\partial_y)f_\downarrow\right]e^{-\frac{e}{\hbar}\rho}

となります。よって求める方程式は

(∂∂x+i∂∂y)f↑(x,y)=0,(∂∂x−i∂∂y)f↓(x,y)=0\left(\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}\right)f_\uparrow(x,y)=0,\qquad \left(\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right)f_\downarrow(x,y)=0

です。ゲージ場が指数因子に押し込められ、残るのは自由な微分方程式になります。

w=x+iyw=x+iy、w‾=x−iy\overline w=x-iy と書くと

∂∂w‾=12(∂∂x+i∂∂y),∂∂w=12(∂∂x−i∂∂y)\frac{\partial}{\partial \overline w}=\frac12\left(\frac{\partial}{\partial x}+i\frac{\partial}{\partial y}\right),\qquad \frac{\partial}{\partial w}=\frac12\left(\frac{\partial}{\partial x}-i\frac{\partial}{\partial y}\right)

なので、上の 2 式はそれぞれ ∂f↑/∂w‾=0\partial f_\uparrow/\partial\overline w=0、∂f↓/∂w=0\partial f_\downarrow/\partial w=0 です。f=u+ivf=u+iv と実部虚部に分けると (∂x+i∂y)f=0(\partial_x+i\partial_y)f=0 は

∂u∂x=∂v∂y,∂v∂x=−∂u∂y\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y},\qquad \frac{\partial v}{\partial x}=-\frac{\partial u}{\partial y}

すなわち w=x+iyw=x+iy を変数とみたときのコーシー・リーマン関係式です。したがって f↑f_\uparrow は ww の正則関数です。同様に (∂x−i∂y)f↓=0(\partial_x-i\partial_y)f_\downarrow=0 は w‾=x−iy\overline w=x-iy を変数とみたコーシー・リーマン関係式で、f↓f_\downarrow は w‾\overline w の正則関数です。

r≥Rr\ge R では ρ=(bR2/2)(1/2+log⁡(r/R))\rho=(bR^2/2)\left(1/2+\log(r/R)\right) なので、ν\nu を

ν≡ebR22ℏ=eΦ2πℏ=Φh/e,Φ≡πR2b\nu\equiv\frac{ebR^2}{2\hbar}=\frac{e\Phi}{2\pi\hbar}=\frac{\Phi}{h/e},\qquad \Phi\equiv\pi R^2 b

(Φ\Phi は円盤を貫く全磁束、h/eh/e は磁束量子)とおくと

eeℏρ=eν/2(rR)ν,e−eℏρ=e−ν/2(rR)−ν(r≥R)e^{\frac{e}{\hbar}\rho}=e^{\nu/2}\left(\frac{r}{R}\right)^{\nu},\qquad e^{-\frac{e}{\hbar}\rho}=e^{-\nu/2}\left(\frac{r}{R}\right)^{-\nu}\qquad (r\ge R)

です。b>0b>0、e>0e>0 なので ν>0\nu>0、つまりスピン上向きの指数因子は r→∞r\to\infty で rνr^{\nu} のように発散し、下向きは r−νr^{-\nu} のように減衰します。

スピン上向きの場合を考えます。f↑f_\uparrow は ww の正則関数なので f↑=∑k≥0ckwkf_\uparrow=\sum_{k\ge0}c_k w^k とべき級数に展開でき、wk=rkeikθw^k=r^ke^{ik\theta} の角度積分で異なる kk は直交するので、各単項式が別々に規格化可能でなければなりません。すると

∫∣ψ↑∣2 d2r=2π∑k∣ck∣2∫0∞r2ke2eℏρ r dr\int |\psi_\uparrow|^2\,d^2r = 2\pi\sum_k |c_k|^2\int_0^\infty r^{2k}e^{\frac{2e}{\hbar}\rho}\,r\,dr

で、r→∞r\to\infty での被積分関数は r2k+2ν+1r^{2k+2\nu+1} となり、k≥0k\ge0、ν>0\nu>0 ではべきが正なので必ず発散します。規格化可能であるためにはすべての ckc_k が 00、すなわち f↑≡0f_\uparrow\equiv0 でなければなりません。したがってスピンが上向きでエネルギーがゼロの束縛状態は存在しません。上向きスピンは磁場中でゼーマンエネルギーを損する側であり、その分が指数因子の発散として現れています。

スピン下向きでは f↓=∑k≥0ckw‾ kf_\downarrow=\sum_{k\ge0}c_k\overline w^{\,k}、w‾ k=rke−ikθ\overline w^{\,k}=r^k e^{-ik\theta} です。角度積分で異なる kk は直交するので

∫∣ψ↓∣2d2r=2π∑k∣ck∣2∫0∞r2ke−2eℏρ r dr\int|\psi_\downarrow|^2 d^2r=2\pi\sum_k |c_k|^2\int_0^\infty r^{2k}e^{-\frac{2e}{\hbar}\rho}\,r\,dr

となり、各単項式が別々に規格化可能である必要があります。r→∞r\to\infty では被積分関数が r2k−2ν+1r^{2k-2\nu+1} なので、収束条件は 2k−2ν+1<−12k-2\nu+1<-1、すなわち

k<ν−1=ebR22ℏ−1.k<\nu-1=\frac{ebR^2}{2\hbar}-1 .

(2k−2ν+1=−12k-2\nu+1=-1 の境界の場合は ∫dr/r\int dr/r の対数発散なので除かれます。原点近傍は e−ebr2/(2ℏ)e^{-ebr^2/(2\hbar)} と r2k+1r^{2k+1} なので問題ありません。)よって一次独立な波動関数は f↓=1,w‾,w‾ 2,…f_\downarrow=1,\overline w,\overline w^{\,2},\dots のうち k<ν−1k<\nu-1 を満たすものだけで、その個数は

N0=⌈ν⌉−1=⌈ebR22ℏ⌉−1(ν>1),N0=0(ν≤1)N_0=\lceil \nu\rceil-1=\left\lceil \frac{ebR^2}{2\hbar}\right\rceil-1\qquad(\nu>1),\qquad N_0=0\quad(\nu\le1)

です。言い換えると、ν\nu が整数でなければ N0=⌊ν⌋N_0=\lfloor\nu\rfloor、整数なら N0=ν−1N_0=\nu-1 個です。ν=Φ/(h/e)\nu=\Phi/(h/e) は円盤を貫く磁束を磁束量子 h/eh/e で測った数なので、ゼロエネルギー束縛状態の個数は磁束量子数の整数部分(境界の場合はそれより 1 つ少ない数)で決まります。これが Aharonov–Casher の結果です。導出で使ったのは ρ\rho の遠方での対数的な振舞い、すなわち総磁束だけであり、磁場分布の詳細(ここでは一定 bb の円盤)には個数が依存しません。

第2問 調和振動子系と1次元連成振動子系の比熱

Section titled “第2問 調和振動子系と1次元連成振動子系の比熱”

自然長 aa、バネ定数 mω2m\omega^2 のバネと質量 mm の質点からなる系を温度 TT の平衡状態で考えます。前半は NN 個の質点が独立に

H=∑j=1N(pj22m+mω22(xj−a)2)H=\sum_{j=1}^{N}\left(\frac{p_j^2}{2m}+\frac{m\omega^2}{2}(x_j-a)^2\right)

に従う場合、後半は NN 個の質点を N+1N+1 個のバネで 1 次元的につなぎ、両端を全体の自然長 (N+1)a(N+1)a で固定した場合です。β=1/(kBT)\beta=1/(k_{\mathrm B}T) と書きます。

1 質点の古典分配関数は、u=x1−au=x_1-a とおいて

z=1h∫dp e−βp2/2m∫du e−βmω2u2/2=1h2πmβ2πβmω2=2πhβωz=\frac1h\int dp\,e^{-\beta p^2/2m}\int du\,e^{-\beta m\omega^2u^2/2} =\frac1h\sqrt{\frac{2\pi m}{\beta}}\sqrt{\frac{2\pi}{\beta m\omega^2}}=\frac{2\pi}{h\beta\omega}

です。z∝β−1z\propto\beta^{-1} なので

U=−∂ln⁡z∂β=1β=kBT,C=∂U∂T=kB.U=-\frac{\partial \ln z}{\partial\beta}=\frac1\beta=k_{\mathrm B}T,\qquad C=\frac{\partial U}{\partial T}=k_{\mathrm B}.

運動エネルギーと位置エネルギーがそれぞれ kBT/2k_{\mathrm B}T/2 を持つ、等分配則そのものです。

位置の平均は、分布 ∝e−βmω2(x1−a)2/2\propto e^{-\beta m\omega^2(x_1-a)^2/2} が x1=ax_1=a について対称なので

⟨x1⟩=a.\langle x_1\rangle=a .

分散は α=βmω2/2\alpha=\beta m\omega^2/2 として与えられた公式を使うと

⟨x12⟩−⟨x1⟩2=⟨(x1−a)2⟩=∫u2e−αu2du∫e−αu2du=12α=kBTmω2.\langle x_1^2\rangle-\langle x_1\rangle^2=\langle (x_1-a)^2\rangle =\frac{\int u^2e^{-\alpha u^2}du}{\int e^{-\alpha u^2}du}=\frac1{2\alpha}=\frac{k_{\mathrm B}T}{m\omega^2}.

これも 12mω2⟨u2⟩=12kBT\frac12 m\omega^2\langle u^2\rangle=\frac12 k_{\mathrm B}T と等分配則から従う値です。

量子論では 1 質点のエネルギー準位が ℏω(n+1/2)\hbar\omega(n+1/2) なので、分配関数は等比級数で

zQ=∑n=0∞e−βℏω(n+1/2)=e−βℏω/21−e−βℏω,z_{\mathrm Q}=\sum_{n=0}^{\infty}e^{-\beta\hbar\omega(n+1/2)}=\frac{e^{-\beta\hbar\omega/2}}{1-e^{-\beta\hbar\omega}}, UQ=−∂ln⁡zQ∂β=ℏω2+ℏωeℏω/kBT−1,U_{\mathrm Q}=-\frac{\partial\ln z_{\mathrm Q}}{\partial\beta} =\frac{\hbar\omega}{2}+\frac{\hbar\omega}{e^{\hbar\omega/k_{\mathrm B}T}-1}, CQ=∂UQ∂T=kB(ℏωkBT)2eℏω/kBT(eℏω/kBT−1)2.C_{\mathrm Q}=\frac{\partial U_{\mathrm Q}}{\partial T} =k_{\mathrm B}\left(\frac{\hbar\omega}{k_{\mathrm B}T}\right)^2 \frac{e^{\hbar\omega/k_{\mathrm B}T}}{\left(e^{\hbar\omega/k_{\mathrm B}T}-1\right)^2}.

位置については X=x1−aX=x_1-a を式 (4) の XX と同一視します。X∣n⟩X|n\rangle は ∣n±1⟩|n\pm1\rangle の重ね合わせで ∣n⟩|n\rangle 成分を持たないので ⟨n∣X∣n⟩=0\langle n|X|n\rangle=0、よって

⟨x1⟩=a\langle x_1\rangle=a

です。分散は式 (4) を 2 回使って ⟨n∣X2∣n⟩=ℏ2mω(2n+1)\langle n|X^2|n\rangle=\dfrac{\hbar}{2m\omega}(2n+1) となるので、⟨n⟩=1/(eβℏω−1)\langle n\rangle=1/(e^{\beta\hbar\omega}-1) を代入して

⟨x12⟩−⟨x1⟩2=ℏ2mω(2⟨n⟩+1)=ℏ2mωcoth⁡ ⁣(ℏω2kBT).\langle x_1^2\rangle-\langle x_1\rangle^2=\frac{\hbar}{2m\omega}\left(2\langle n\rangle+1\right) =\frac{\hbar}{2m\omega}\coth\!\left(\frac{\hbar\omega}{2k_{\mathrm B}T}\right).

比較します。高温極限 kBT≫ℏωk_{\mathrm B}T\gg\hbar\omega では eβℏω−1≃βℏωe^{\beta\hbar\omega}-1\simeq\beta\hbar\omega より

CQ→kB(βℏω)21(βℏω)2=kB=CC_{\mathrm Q}\to k_{\mathrm B}(\beta\hbar\omega)^2\frac{1}{(\beta\hbar\omega)^2}=k_{\mathrm B}=C

で古典結果に一致します(分散も coth⁡x≃1/x\coth x\simeq1/x より kBT/(mω2)k_{\mathrm B}T/(m\omega^2) に一致)。低温極限 kBT≪ℏωk_{\mathrm B}T\ll\hbar\omega では

CQ≃kB(ℏωkBT)2e−ℏω/kBT⟶0C_{\mathrm Q}\simeq k_{\mathrm B}\left(\frac{\hbar\omega}{k_{\mathrm B}T}\right)^2e^{-\hbar\omega/k_{\mathrm B}T}\longrightarrow 0

と指数関数的に消え、C=kBC=k_{\mathrm B} とはまったく違います。準位間隔 ℏω\hbar\omega より熱エネルギーが小さいと励起が凍結するためで、古典論が第三法則を満たさないことの現れです。分散も ℏ/(2mω)\hbar/(2m\omega) という零点ゆらぎに飽和し、T→0T\to0 でも 00 になりません。

平衡位置からのずれ uj≡xj−jau_j\equiv x_j-ja(j=1,…,Nj=1,\dots,N)を導入し、u0≡0u_0\equiv0、uN+1≡0u_{N+1}\equiv0 と約束します。すると x1−a=u1x_1-a=u_1、xj−xj−1−a=uj−uj−1x_j-x_{j-1}-a=u_j-u_{j-1}、(N+1)a−xN−a=−uN(N+1)a-x_N-a=-u_N なので、式 (5) は

H=∑j=1Npj22m+mω22∑j=1N+1(uj−uj−1)2H=\sum_{j=1}^{N}\frac{p_j^2}{2m}+\frac{m\omega^2}{2}\sum_{j=1}^{N+1}(u_j-u_{j-1})^2

と書けます。ポテンシャルは mω22∑i,juiKijuj\dfrac{m\omega^2}{2}\sum_{i,j}u_iK_{ij}u_j の形で、KK は対角成分 22、隣接成分 −1-1 の三重対角行列です。

KK の固有ベクトルを vj(ℓ)=2N+1sin⁡ ⁣(jℓπN+1)v_j^{(\ell)}=\sqrt{\dfrac{2}{N+1}}\sin\!\left(\dfrac{j\ell\pi}{N+1}\right) とすると、θℓ=ℓπ/(N+1)\theta_\ell=\ell\pi/(N+1) として

−vj−1(ℓ)+2vj(ℓ)−vj+1(ℓ)=(2−2cos⁡θℓ) vj(ℓ)-v_{j-1}^{(\ell)}+2v_j^{(\ell)}-v_{j+1}^{(\ell)}=(2-2\cos\theta_\ell)\,v_j^{(\ell)}

が成り立ちます(sin⁡((j±1)θ)\sin((j\pm1)\theta) を加法定理で展開すれば直ちに出ます)。境界条件は v0(ℓ)=0v_0^{(\ell)}=0 が自動的に満たされ、vN+1(ℓ)=sin⁡((N+1)θℓ)=0v_{N+1}^{(\ell)}=\sin((N+1)\theta_\ell)=0 が θℓ=ℓπ/(N+1)\theta_\ell=\ell\pi/(N+1) から満たされます。独立な非自明解は ℓ=1,2,…,N\ell=1,2,\dots,N の NN 個です。さらに

∑j=1Nvj(ℓ)vj(ℓ′)=δℓℓ′\sum_{j=1}^{N}v_j^{(\ell)}v_j^{(\ell')}=\delta_{\ell\ell'}

なので、{v(ℓ)}\{v^{(\ell)}\} は正規直交基底をなします。

そこで直交変換 uj=∑ℓvj(ℓ)Qℓu_j=\sum_\ell v_j^{(\ell)}Q_\ell、pj=∑ℓvj(ℓ)Pℓp_j=\sum_\ell v_j^{(\ell)}P_\ell を行うと、直交性から ∑jpj2=∑ℓPℓ2\sum_j p_j^2=\sum_\ell P_\ell^2 であり、ポテンシャルは固有値で対角化されて

H=∑ℓ=1N(Pℓ22m+mωℓ22Qℓ2),ωℓ=ω2−2cos⁡(ℓN+1π)H=\sum_{\ell=1}^{N}\left(\frac{P_\ell^2}{2m}+\frac{m\omega_\ell^2}{2}Q_\ell^2\right),\qquad \omega_\ell=\omega\sqrt{2-2\cos\left(\frac{\ell}{N+1}\pi\right)}

となります。これで系が NN 個の独立な基準振動の和で表され、固有振動数が式 (6) で与えられることが示せました。なお ωℓ=2ω∣sin⁡ℓπ2(N+1)∣\omega_\ell=2\omega\left|\sin\dfrac{\ell\pi}{2(N+1)}\right| と書けて、0<ωℓ<2ω0<\omega_\ell<2\omega です。

古典論では、系は NN 個の独立な 1 次元調和振動子(振動数はばらばらでよい)の集まりです。各振動子は運動量と座標の 2 つの 2 次形式を持つので等分配則から kBTk_{\mathrm B}T を与え、内部エネルギーは NkBTNk_{\mathrm B}T、質点 1 個当たりの比熱は

C1=kBC_1=k_{\mathrm B}

です。振動数 ωℓ\omega_\ell に依存せず、設問1の CC と同じ値になります。古典比熱は 2 次形式の自由度の数だけで決まり、バネのつなぎ方には依存しません。

質点 1 個当たりの量子比熱は、設問2の 1 振動子の結果を基準振動について足したものです。xℓ=βℏωℓx_\ell=\beta\hbar\omega_\ell と書くと

C1Q=kBN∑ℓ=1Nxℓ2exℓ(exℓ−1)2.C_{1\mathrm Q}=\frac{k_{\mathrm B}}{N}\sum_{\ell=1}^{N}\frac{x_\ell^2 e^{x_\ell}}{(e^{x_\ell}-1)^2}.

被和関数 x2ex/(ex−1)2x^2e^x/(e^x-1)^2 は x≫1x\gg1 で指数関数的に小さくなるので、低温 kBT≪ℏωk_{\mathrm B}T\ll\hbar\omega では ℏωℓ≲kBT\hbar\omega_\ell\lesssim k_{\mathrm B}T を満たす ωℓ≪ω\omega_\ell\ll\omega、すなわち ℓ≪N\ell\ll N の長波長モードだけが寄与します。そこでは

ωℓ=2ωsin⁡ℓπ2(N+1)≃πω ℓN+1\omega_\ell=2\omega\sin\frac{\ell\pi}{2(N+1)}\simeq\frac{\pi\omega\,\ell}{N+1}

と線形分散になります。NN が非常に大きければ xℓx_\ell の間隔 Δ=πβℏω/(N+1)\Delta=\pi\beta\hbar\omega/(N+1) は 11 に比べて微小なので、和を積分に置き換えられます(この近似の条件は πℏω/(N+1)≪kBT≪ℏω\pi\hbar\omega/(N+1)\ll k_{\mathrm B}T\ll\hbar\omega です)。

C1Q≃kBN⋅1Δ∫0∞x2ex(ex−1)2dx=kBN⋅N+1πβℏω⋅π23≃πkB3βℏω.C_{1\mathrm Q}\simeq\frac{k_{\mathrm B}}{N}\cdot\frac{1}{\Delta}\int_0^{\infty}\frac{x^2e^x}{(e^x-1)^2}dx =\frac{k_{\mathrm B}}{N}\cdot\frac{N+1}{\pi\beta\hbar\omega}\cdot\frac{\pi^2}{3} \simeq\frac{\pi k_{\mathrm B}}{3\beta\hbar\omega}.

与えられた公式から ∫0∞x2ex(ex−1)2dx=2⋅π26=π23\displaystyle\int_0^\infty\frac{x^2e^x}{(e^x-1)^2}dx=2\cdot\frac{\pi^2}{6}=\frac{\pi^2}{3} を使いました。1/β=kBT1/\beta=k_{\mathrm B}T を戻すと

C1Q≃πkB23ℏω T,C_{1\mathrm Q}\simeq\frac{\pi k_{\mathrm B}^2}{3\hbar\omega}\,T,

すなわち

b=1,A=πkB23ℏωb=1,\qquad A=\frac{\pi k_{\mathrm B}^2}{3\hbar\omega}

です。1 次元デバイ模型の T1T^1 則で、次元も [kB2T/(ℏω)]=(J/K)2⋅K/J=J/K[k_{\mathrm B}^2T/(\hbar\omega)]=(\mathrm{J/K})^2\cdot\mathrm{K}/\mathrm{J}=\mathrm{J/K} と比熱の次元になっています。設問4の古典値 kBk_{\mathrm B} と比べると、kBT≪ℏωk_{\mathrm B}T\ll\hbar\omega で C1Q≪C1C_{1\mathrm Q}\ll C_1 であり、独立振動子の場合の指数関数的な消え方(設問2)と違ってべき的にしか消えないのは、いくら低温でも ωℓ→0\omega_\ell\to0 の音響モードが必ず励起されているからです。

外力を受けず xx 軸上を運動するロケットが、自分から見て常に一定の速度 −vgas-v_{\mathrm{gas}}(vgas>0v_{\mathrm{gas}}>0、進行方向と逆向き)でガスを噴射します。噴射量は、ロケットが常に一定の加速度 α\alpha を持つように自動調整されます。ある慣性系(この系での速度を対地速度と呼ぶ)で議論します。前半は非相対論、後半は特殊相対論です。γ(v)=1/1−v2/c2\gamma(v)=1/\sqrt{1-v^2/c^2} と書きます。

時刻 tt から t+dtt+dt の間に、質量は m→m+dmm\to m+dm(dm<0dm<0)、速度は v→v+dvv\to v+dv になり、質量 −dm>0-dm>0 のガスが対地速度 v−vgasv-v_{\mathrm{gas}} で放出されます。運動量保存則は

mv=(m+dm)(v+dv)+(−dm)(v−vgas)mv=(m+dm)(v+dv)+(-dm)(v-v_{\mathrm{gas}})

で、微小量の 2 次を落として整理すると

0=m dv+vgas dm.0=m\,dv+v_{\mathrm{gas}}\,dm .

dv/dt=αdv/dt=\alpha を代入して

dmdt=−αvgas m(t)  (<0).\frac{dm}{dt}=-\frac{\alpha}{v_{\mathrm{gas}}}\,m(t)\;(<0).

上の式は dln⁡m/dt=−α/vgasd\ln m/dt=-\alpha/v_{\mathrm{gas}} なので、m(0)=m0m(0)=m_0 のもとで

m(t)=m0exp⁡(−αtvgas).m(t)=m_0\exp\left(-\frac{\alpha t}{v_{\mathrm{gas}}}\right).

v(t)=αtv(t)=\alpha t を使うと v=vgasln⁡(m0/m)v=v_{\mathrm{gas}}\ln(m_0/m)、ツィオルコフスキーの式に一致します。

イ について。ds2=−c2dt2+dx2ds^2=-c^2dt^2+dx^2 より dτ=−ds2/c=dt2−dx2/c2d\tau=\sqrt{-ds^2}/c=\sqrt{dt^2-dx^2/c^2} なので、入るのは

1−v2c2  (=1γ(v)).\sqrt{1-\frac{v^2}{c^2}}\;\left(=\frac{1}{\gamma(v)}\right).

ロ について。uμ⊙uμ=−c2dt2+dx2dτ2=ds2dτ2u^\mu\odot u^\mu=\dfrac{-c^2dt^2+dx^2}{d\tau^2}=\dfrac{ds^2}{d\tau^2} で、dτ2=−ds2/c2d\tau^2=-ds^2/c^2 なので、入るのは −c2-c^2 です。

ハ について。ロの式を τ\tau で微分すると 2 uμ⊙aμ=02\,u^\mu\odot a^\mu=0 なので、入るのは 00 です。

ニ と ホ について。ハより −u0a0+u1a1=0-u^0a^0+u^1a^1=0 すなわち a0=u1a1/u0a^0=u^1a^1/u^0。これを −(a0)2+(a1)2=α2-(a^0)^2+(a^1)^2=\alpha^2 に入れると

(a1)2(u0)2−(u1)2(u0)2=α2.(a^1)^2\frac{(u^0)^2-(u^1)^2}{(u^0)^2}=\alpha^2 .

ロより (u0)2−(u1)2=c2(u^0)^2-(u^1)^2=c^2 なので、a1>0a^1>0 をとって、ニ と ホ に入るのはそれぞれ

a0=α u1c,a1=α u0c.a^0=\frac{\alpha\,u^1}{c},\qquad a^1=\frac{\alpha\,u^0}{c}.

瞬間静止系では uμ=(c,0)u^\mu=(c,0) なので aμ=(0,α)a^\mu=(0,\alpha) となり、α\alpha が固有加速度であることが確かめられます。

設問3の結果は

du0dτ=αcu1,du1dτ=αcu0\frac{du^0}{d\tau}=\frac{\alpha}{c}u^1,\qquad \frac{du^1}{d\tau}=\frac{\alpha}{c}u^0

です。和と差をとると ddτ(u0±u1)=±αc(u0±u1)\dfrac{d}{d\tau}(u^0\pm u^1)=\pm\dfrac{\alpha}{c}(u^0\pm u^1) なので、τ=0\tau=0 で v=0v=0(u0=c, u1=0u^0=c,\ u^1=0)という初期条件のもとで

u0±u1=c e±ατ/c,u^0\pm u^1=c\,e^{\pm\alpha\tau/c},

すなわち

u0=ccosh⁡ατc,u1=csinh⁡ατc.u^0=c\cosh\frac{\alpha\tau}{c},\qquad u^1=c\sinh\frac{\alpha\tau}{c}.

対地速度は v=dxdt=dx/dτdt/dτ=c u1u0v=\dfrac{dx}{dt}=\dfrac{dx/d\tau}{dt/d\tau}=\dfrac{c\,u^1}{u^0} なので

v(τ)=ctanh⁡(ατc).v(\tau)=c\tanh\left(\frac{\alpha\tau}{c}\right).

ατ≪c\alpha\tau\ll c では v≃ατv\simeq\alpha\tau で非相対論の等加速度運動に戻り、τ→∞\tau\to\infty で v→cv\to c に漸近して光速を超えません。

4元運動量は pμ=muμ=(mγc, mγv)p^\mu=m u^\mu=(m\gamma c,\ m\gamma v) です。噴射前のロケットは (m,v)(m,v)、噴射後のロケットは (m+dm, v+dv)(m+dm,\ v+dv)、ガスは静止質量 dm′dm'、対地速度 v′v' なので、保存則は

m uμ(v)=(m+dm) uμ(v+dv)+dm′ uμ(v′)m\,u^\mu(v)=(m+dm)\,u^\mu(v+dv)+dm'\,u^\mu(v')

です。dγdv=γ3vc2\dfrac{d\gamma}{dv}=\gamma^3\dfrac{v}{c^2} と、γ3v2c2+γ=γ3\gamma^3\dfrac{v^2}{c^2}+\gamma=\gamma^3 を使って微小量の 1 次まで展開すると、時間成分(エネルギー)と空間成分(運動量)がそれぞれ

γ dm+mγ3vc2 dv+γ′ dm′=0,\gamma\,dm+m\gamma^3\frac{v}{c^2}\,dv+\gamma'\,dm'=0, γv dm+mγ3 dv+γ′v′ dm′=0\gamma v\,dm+m\gamma^3\,dv+\gamma' v'\,dm'=0

となります。ここで γ=γ(v)\gamma=\gamma(v)、γ′=γ(v′)\gamma'=\gamma(v') です。これが求める保存則です。

上の 2 式から dm′dm' を消去します。第1式に v′v' を掛けて第2式を引くと

γ(v′−v) dm+mγ3(vv′c2−1)dv=0  ⟹  dmm=γ2 1−vv′/c2v′−v dv.\gamma(v'-v)\,dm+m\gamma^3\left(\frac{vv'}{c^2}-1\right)dv=0 \;\Longrightarrow\; \frac{dm}{m}=\gamma^2\,\frac{1-vv'/c^2}{v'-v}\,dv .

速度変換則 −vgas=v′−v1−vv′/c2-v_{\mathrm{gas}}=\dfrac{v'-v}{1-vv'/c^2} より 1−vv′/c2v′−v=−1vgas\dfrac{1-vv'/c^2}{v'-v}=-\dfrac{1}{v_{\mathrm{gas}}} なので

dmm=−γ2 dvvgas.\frac{dm}{m}=-\frac{\gamma^2\,dv}{v_{\mathrm{gas}}}.

ここに設問4の v=ctanh⁡(ατ/c)v=c\tanh(\alpha\tau/c) を入れます。χ≡ατ/c\chi\equiv\alpha\tau/c とすると γ=cosh⁡χ\gamma=\cosh\chi、dv=c sech2χ dχdv=c\,\mathrm{sech}^2\chi\,d\chi なので γ2dv=c dχ=α dτ\gamma^2dv=c\,d\chi=\alpha\,d\tau となり

dmm=−αvgas dτ  ⟹  m(τ)=m0exp⁡(−ατvgas).\frac{dm}{m}=-\frac{\alpha}{v_{\mathrm{gas}}}\,d\tau \;\Longrightarrow\; m(\tau)=m_0\exp\left(-\frac{\alpha\tau}{v_{\mathrm{gas}}}\right).

設問2の答え m=m0e−αt/vgasm=m_0e^{-\alpha t/v_{\mathrm{gas}}} と比べると、形は完全に同じで、慣性系の時刻 tt が固有時間 τ\tau に置き換わっただけです。これは自然な結果です。α\alpha も vgasv_{\mathrm{gas}} もロケットの瞬間静止系で定義された量であり、その系では微小時間 dτd\tau の間の運動はニュートン力学そのものなので、設問1の関係式がそのまま dτd\tau について成り立つからです。相対論的補正はすべて t→τt\to\tau に吸収されます。

物理的な帰結を 2 つ挙げます。第一に、ατ/c=χ=artanh(v/c)\alpha\tau/c=\chi=\mathrm{artanh}(v/c) を使うと必要な質量比は

m0m=exp⁡(cvgas artanhvc)=(1+v/c1−v/c)c/(2vgas)\frac{m_0}{m}=\exp\left(\frac{c}{v_{\mathrm{gas}}}\,\mathrm{artanh}\frac{v}{c}\right) =\left(\frac{1+v/c}{1-v/c}\right)^{c/(2v_{\mathrm{gas}})}

で、非相対論の ev/vgase^{v/v_{\mathrm{gas}}} より大きく、v→cv\to c で発散します。有限の燃料で光速には到達できません。第二に、時間の遅れにより τ<t\tau<t なので、地上時刻 tt で比べると相対論的ロケットの方が燃料の消費は少ないものの、速度も αt\alpha t に届きません。なお設問5の 2 式から dm′dm' を求めると dm′=−dm1−vgas2/c2<−dmdm'=-dm\sqrt{1-v_{\mathrm{gas}}^2/c^2}<-dm となり、静止質量は保存しません。差分はガスの運動エネルギーに変わっています。

第4問 プラズマ中の電磁波と干渉計による電子密度計測

Section titled “第4問 プラズマ中の電磁波と干渉計による電子密度計測”

磁場のないプラズマ中で、電子密度 nen_{\mathrm e}、カットオフ密度 nc=meε0ω2/e2n_{\mathrm c}=m_{\mathrm e}\varepsilon_0\omega^2/e^2 として屈折率が α=1−ne/nc\alpha=\sqrt{1-n_{\mathrm e}/n_{\mathrm c}}、波数が k=αω/ck=\alpha\omega/c で与えられます。電場の偏光成分は E(x⃗,t)=E0exp⁡(ik⃗⋅x⃗−iωt)E(\vec x,t)=E_0\exp(i\vec k\cdot\vec x-i\omega t) です。ne/nc=ωpe2/ω2n_{\mathrm e}/n_{\mathrm c}=\omega_{\mathrm{pe}}^2/\omega^2 であることに注意します。ここで ωpe=nee2/(ε0me)\omega_{\mathrm{pe}}=\sqrt{n_{\mathrm e}e^2/(\varepsilon_0m_{\mathrm e})} はプラズマ振動数です。数値は me=9.1×10−31 kgm_{\mathrm e}=9.1\times10^{-31}\,\mathrm{kg}、e=1.6×10−19 Ce=1.6\times10^{-19}\,\mathrm{C}、ε0=8.9×10−12 F/m\varepsilon_0=8.9\times10^{-12}\,\mathrm{F/m}、c=3.0×108 m/sc=3.0\times10^{8}\,\mathrm{m/s} を使います。

z≤0z\le0 が真空、z>0z>0 が一様密度 nen_{\mathrm e} のプラズマで、電磁波は zz 軸に沿って入射します。z>0z>0 での電場は E=E0eikz−iωtE=E_0e^{ikz-i\omega t}、k=(ω/c)1−ne/nck=(\omega/c)\sqrt{1-n_{\mathrm e}/n_{\mathrm c}} です。

ne<ncn_{\mathrm e}<n_{\mathrm c} の場合、α\alpha は 0<α<10<\alpha<1 の実数なので kk も実数で、電場は振幅一定のまま zz 方向に伝搬します。空間変化の特徴的な長さは波長

λ=2πk=2πcω1−ne/nc=λvac1−ne/nc\lambda=\frac{2\pi}{k}=\frac{2\pi c}{\omega\sqrt{1-n_{\mathrm e}/n_{\mathrm c}}}=\frac{\lambda_{\mathrm{vac}}}{\sqrt{1-n_{\mathrm e}/n_{\mathrm c}}}

で、真空中の波長 λvac=2πc/ω\lambda_{\mathrm{vac}}=2\pi c/\omega より長くなります。位相速度 ω/k=c/α\omega/k=c/\alpha は cc より大きく、群速度 αc\alpha c は cc より小さいです。

ne>ncn_{\mathrm e}>n_{\mathrm c} の場合、α\alpha は純虚数になります。減衰する解を選んで α=ine/nc−1\alpha=i\sqrt{n_{\mathrm e}/n_{\mathrm c}-1} とすれば

E=E0exp⁡(−zδ)e−iωt,δ=cωne/nc−1=cωpe2−ω2E=E_0\exp\left(-\frac{z}{\delta}\right)e^{-i\omega t},\qquad \delta=\frac{c}{\omega\sqrt{n_{\mathrm e}/n_{\mathrm c}-1}}=\frac{c}{\sqrt{\omega_{\mathrm{pe}}^2-\omega^2}}

となり、電場は位相を進めずに表面から長さ δ\delta で指数関数的に減衰します(エバネッセント波)。伝搬できないので入射波はほぼ全反射され、これがカットオフです。ω≪ωpe\omega\ll\omega_{\mathrm{pe}} では δ→c/ωpe\delta\to c/\omega_{\mathrm{pe}}(プラズマの表皮厚さ)になります。

半径 aa の円柱プラズマの断面の中心を原点とし、直径方向に zz 軸をとります。波長とビーム半径が aa に比べて十分小さいので幾何光学が使えて、経路上の点は円柱軸から距離 r=∣z∣r=|z| にあります。位相変化は

ϕp=∫−aak dz=ωc∫−aa1−ne(r)nc dz.\phi_{\mathrm p}=\int_{-a}^{a}k\,dz=\frac{\omega}{c}\int_{-a}^{a}\sqrt{1-\frac{n_{\mathrm e}(r)}{n_{\mathrm c}}}\,dz .

ne(r)≪ncn_{\mathrm e}(r)\ll n_{\mathrm c} として平方根を 1 次まで展開すると

ϕp≃ωc∫−aa(1−ne(r)2nc)dz=ωc⋅2a−ω2cnc∫−aane(r) dz.\phi_{\mathrm p}\simeq\frac{\omega}{c}\int_{-a}^{a}\left(1-\frac{n_{\mathrm e}(r)}{2n_{\mathrm c}}\right)dz =\frac{\omega}{c}\cdot2a-\frac{\omega}{2cn_{\mathrm c}}\int_{-a}^{a}n_{\mathrm e}(r)\,dz .

真空の場合は ϕv=(ω/c)⋅2a\phi_{\mathrm v}=(\omega/c)\cdot2a なので

Δϕ=ϕv−ϕp=ω2cnc∫−aane(r) dz=e22ε0mec ω∫−aane(r) dz.\Delta\phi=\phi_{\mathrm v}-\phi_{\mathrm p}=\frac{\omega}{2cn_{\mathrm c}}\int_{-a}^{a}n_{\mathrm e}(r)\,dz =\frac{e^2}{2\varepsilon_0 m_{\mathrm e}c\,\omega}\int_{-a}^{a}n_{\mathrm e}(r)\,dz .

nc=meε0ω2/e2n_{\mathrm c}=m_{\mathrm e}\varepsilon_0\omega^2/e^2 を代入しました。Δϕ>0\Delta\phi>0、つまりプラズマがあると位相の進みが減ります。この 1 次近似では Δϕ\Delta\phi が電子密度の線積分だけで決まり、密度分布の形によらないというのが計測原理です。Δϕ∝1/ω∝λ\Delta\phi\propto1/\omega\propto\lambda なので、波長が長いほど感度が上がります。

ne(r)=ne0[1−(r/a)2]n_{\mathrm e}(r)=n_{\mathrm e0}[1-(r/a)^2]、a=0.1 ma=0.1\,\mathrm{m}、ne0=1×1019 m−3n_{\mathrm e0}=1\times10^{19}\,\mathrm{m^{-3}} とします。

周波数の下限を決めるのはカットオフです。電磁波が経路上のどこでも伝搬できなければ透過測定にならないので、最も密度が高い中心 r=0r=0 で ne0<nc(ω)n_{\mathrm e0}<n_{\mathrm c}(\omega)、すなわち

ω>ωpe(ne0)=ne0e2ε0me\omega>\omega_{\mathrm{pe}}(n_{\mathrm e0})=\sqrt{\frac{n_{\mathrm e0}e^2}{\varepsilon_0 m_{\mathrm e}}}

が必要です。これより低い周波数では中心部で α\alpha が虚数になり、波は反射されて反対側に抜けません。さらに設問2の 1 次近似(Δϕ\Delta\phi が線積分に比例すること)を成り立たせるには ne0≪ncn_{\mathrm e0}\ll n_{\mathrm c} が望ましく、この意味でも下限はカットオフ周波数で与えられます。数値は

ωpe=1019×(1.6×10−19)28.9×10−12×9.1×10−31=1.8×1011 rad/s,fmin⁡=ωpe2π≃3×1010 Hz\omega_{\mathrm{pe}}=\sqrt{\frac{10^{19}\times(1.6\times10^{-19})^2}{8.9\times10^{-12}\times9.1\times10^{-31}}} =1.8\times10^{11}\,\mathrm{rad/s},\qquad f_{\min}=\frac{\omega_{\mathrm{pe}}}{2\pi}\simeq3\times10^{10}\,\mathrm{Hz}

で、オーダーは 1010 Hz10^{10}\,\mathrm{Hz}(数十 GHz、ミリ波)です。

この周波数 f≃3×1010 Hzf\simeq3\times10^{10}\,\mathrm{Hz}(ω=1.9×1011 rad/s\omega=1.9\times10^{11}\,\mathrm{rad/s}、nc=1.1×1019 m−3n_{\mathrm c}=1.1\times10^{19}\,\mathrm{m^{-3}})で Δϕ\Delta\phi を求めます。経路上で r=∣z∣r=|z| なので

∫−aane dz=ne0∫−aa(1−z2a2)dz=43ne0a=43×1019×0.1=1.3×1018 m−2\int_{-a}^{a}n_{\mathrm e}\,dz=n_{\mathrm e0}\int_{-a}^{a}\left(1-\frac{z^2}{a^2}\right)dz=\frac43 n_{\mathrm e0}a =\frac43\times10^{19}\times0.1=1.3\times10^{18}\,\mathrm{m^{-2}}

であり

Δϕ=2e2ne0a3ε0mec ω=2×(1.6×10−19)2×1019×0.13×8.9×10−12×9.1×10−31×3.0×108×1.9×1011≃4×101 rad.\Delta\phi=\frac{2e^2n_{\mathrm e0}a}{3\varepsilon_0m_{\mathrm e}c\,\omega} =\frac{2\times(1.6\times10^{-19})^2\times10^{19}\times0.1} {3\times8.9\times10^{-12}\times9.1\times10^{-31}\times3.0\times10^{8}\times1.9\times10^{11}} \simeq4\times10^{1}\,\mathrm{rad}.

答えは Δϕ≈4×101 rad\Delta\phi\approx4\times10^1\,\mathrm{rad}(約 40 ラジアン)です。2π2\pi よりはるかに大きいので、実際には位相を連続追跡(アンラップ)して縞の本数を数える必要があります。

参照波とプローブ波の重ね合わせは Et=E1e−iωt+E2e−iωt+iΔϕE_{\mathrm t}=E_1e^{-i\omega t}+E_2e^{-i\omega t+i\Delta\phi} なので

∣Et∣2=∣E1+E2eiΔϕ∣2=E12+E22+2E1E2cos⁡Δϕ.|E_{\mathrm t}|^2=\left|E_1+E_2e^{i\Delta\phi}\right|^2=E_1^2+E_2^2+2E_1E_2\cos\Delta\phi .

ここで E1,E2E_1,E_2 は実の振幅としました。

参照波を ω1\omega_1、プローブ波を ω2\omega_2 にずらすと Et=E1e−iω1t+E2e−iω2t+iΔϕ(t)E_{\mathrm t}=E_1e^{-i\omega_1t}+E_2e^{-i\omega_2t+i\Delta\phi(t)} で、検出器出力は

∣Et∣2=E12+E22+2E1E2cos⁡[Δω t−Δϕ(t)],Δω=ω2−ω1|E_{\mathrm t}|^2=E_1^2+E_2^2+2E_1E_2\cos\left[\Delta\omega\,t-\Delta\phi(t)\right],\qquad \Delta\omega=\omega_2-\omega_1

となります。検出器は Δω\Delta\omega に追随できるので、この角周波数 Δω\Delta\omega のビート(うなり)が信号として観測されます。

(a) が解決できる理由。位相情報はビート成分の位相にのみ現れ、E1,E2E_1,E_2 のゆらぎや迷光などの位相に無関係な成分は直流成分とビートの振幅に現れます。したがって Δω\Delta\omega 近傍の帯域だけを取り出し、別に用意した Δω\Delta\omega の基準信号との位相差(例えばゼロ交差の時刻)を測れば、振幅の時間変動とは独立に Δϕ(t)\Delta\phi(t) が決まります。設問4の Δω=0\Delta\omega=0 の場合は出力の大きさそのものが Δϕ\Delta\phi を担うので、振幅変動と区別できませんでした。

(b) が解決できる理由。Δω=0\Delta\omega=0 では出力が cos⁡Δϕ\cos\Delta\phi で Δϕ→−Δϕ\Delta\phi\to-\Delta\phi に対して不変(偶関数)なので符号が決まりません。ビートがあると位相は Δωt−Δϕ(t)\Delta\omega t-\Delta\phi(t) と時間の関数になり、瞬時周波数は Δω−d(Δϕ)/dt\Delta\omega-d(\Delta\phi)/dt です。Δϕ\Delta\phi が増えるときはビートが遅れ(周波数が下がり)、減るときは進む(周波数が上がる)ので、その向きから符号が判別できます。位相を連続的に追跡できるので 2π2\pi を超える変化も測れます。

必要な条件は

Δω≫∣d(Δϕ)dt∣\Delta\omega\gg\left|\frac{d(\Delta\phi)}{dt}\right|

です。これが破れると瞬時周波数 Δω−d(Δϕ)/dt\Delta\omega-d(\Delta\phi)/dt が符号を変えたり、ビートの 1 周期の間に Δϕ\Delta\phi が大きく変わってしまい、位相を定義して読み取ることができません。言い換えれば、測りたい密度変化の帯域より十分高いビート周波数を選ぶ必要があります。

第5問 高エネルギー電子弾性散乱と原子核の荷電半径

Section titled “第5問 高エネルギー電子弾性散乱と原子核の荷電半径”

核子や原子核の荷電半径を高エネルギー電子の弾性散乱で測る問題です。c=3×108 m/sc=3\times10^{8}\,\mathrm{m/s}、ℏc=0.2 GeV⋅fm\hbar c=0.2\,\mathrm{GeV\cdot fm}、陽子質量 Mc2=1 GeVMc^2=1\,\mathrm{GeV} を使い、概算は有効数字 1 桁とします。

(a) 静電加速器は適しません。電極間の静電圧で得られるエネルギーは端子電圧そのままなので、数百 MeV には数百 MV の直流電圧が必要です。実際には絶縁破壊(放電)のため端子電圧は十数 MV 程度が限界で、タンデム型でもその 2 倍程度にしかなりません。

(b) サイクロトロンは適しません。電子の静止エネルギーは 0.511 MeV0.511\,\mathrm{MeV} しかないので、数 MeV でもすでに γ≫1\gamma\gg1 です。回転周期 T=2πγme/(eB)T=2\pi\gamma m_{\mathrm e}/(eB) は γ\gamma に比例して伸びるため、一定磁場・一定周波数では加速電圧との同期(等時性)がすぐに崩れます。数百 MeV では γ∼600\gamma\sim600 にもなり、まったく同期できません。

(c) 線形加速器は適します。電子は 1 MeV 程度で既に v≃cv\simeq c になるので、高周波の進行波の位相速度を cc に合わせておけば全区間で同じ位相に乗せ続けられ、電極を継ぎ足すだけでエネルギーを積み上げられます。円軌道を回らないのでシンクロトロン放射による損失もありません。実際に数百 MeV から数十 GeV の電子加速には線形加速器が使われます。

入射電子数は Ne=It/eN_{\mathrm e}=It/e です。標的中の原子核の数密度は ρNA/A\rho N_{\mathrm A}/A なので、厚さ ℓ\ell の板の単位面積当たりの原子核数(面密度)は ρℓNA/A\rho\ell N_{\mathrm A}/A です。薄い標的で 1 回散乱のみを考えれば、立体角 ω\omega の検出器に入る電子数は

n=Ne×ρℓNAA×dσdΩ(θ)×ωn=N_{\mathrm e}\times\frac{\rho\ell N_{\mathrm A}}{A}\times\frac{d\sigma}{d\Omega}(\theta)\times\omega

です。これを解いて

dσdΩ=n e AI t ρ ℓ NA ω.\frac{d\sigma}{d\Omega}=\frac{n\,e\,A}{I\,t\,\rho\,\ell\,N_{\mathrm A}\,\omega}.

θ\theta は測定器を置いた方向を指定するだけで式には現れず、ZZ と cc も不要です。次元は、ρ\rho を g/cm3\mathrm{g/cm^3}、AA を g/mol\mathrm{g/mol}、NAN_{\mathrm A} を mol−1\mathrm{mol^{-1}} で測れば cm2/sr\mathrm{cm^2/sr} となり断面積になります。

一様磁場中の円運動では p=eBrp=eBr なので B=p/(er)B=p/(er) です。p=300 MeV/cp=300\,\mathrm{MeV}/c を SI に直すと

p=300×106×1.6×10−193×108=1.6×10−19 kg⋅m/s,p=\frac{300\times10^{6}\times1.6\times10^{-19}}{3\times10^{8}}=1.6\times10^{-19}\,\mathrm{kg\cdot m/s}, B=1.6×10−191.6×10−19×1.0≃1 T.B=\frac{1.6\times10^{-19}}{1.6\times10^{-19}\times1.0}\simeq1\,\mathrm{T}.

答えは B≈1 TB\approx1\,\mathrm{T} です。B[T]=p[GeV/c]/(0.3 r[m])B[\mathrm T]=p[\mathrm{GeV}/c]/(0.3\,r[\mathrm m]) という覚え方でも同じ値になります。

電子質量を無視すると電子の 4 元運動量は光子と同じく p2=0p^2=0 を満たし、∣p⃗∣=E/c|\vec p|=E/c です。静止した陽子(4 元運動量 P=(Mc,0⃗)P=(Mc,\vec 0))との弾性散乱で、終状態の陽子の 4 元運動量が P′=P+pi−pfP'=P+p_{\mathrm i}-p_{\mathrm f}、P′2=P2P'^2=P^2 であることを使うと

pi⋅pf=P⋅(pi−pf)  ⟹  EiEfc2(1−cos⁡θ)=M(Ei−Ef)p_{\mathrm i}\cdot p_{\mathrm f}=P\cdot(p_{\mathrm i}-p_{\mathrm f}) \;\Longrightarrow\; \frac{E_{\mathrm i}E_{\mathrm f}}{c^2}(1-\cos\theta)=M(E_{\mathrm i}-E_{\mathrm f})

が得られます。両辺を MEiEfME_{\mathrm i}E_{\mathrm f} で割ると、コンプトン散乱と同じ形の関係

1Ef−1Ei=1−cos⁡θMc2⟺Ef=Ei1+EiMc2(1−cos⁡θ)\frac{1}{E_{\mathrm f}}-\frac{1}{E_{\mathrm i}}=\frac{1-\cos\theta}{Mc^2} \qquad\Longleftrightarrow\qquad E_{\mathrm f}=\frac{E_{\mathrm i}}{1+\dfrac{E_{\mathrm i}}{Mc^2}(1-\cos\theta)}

になります。これに Ei=300 MeVE_{\mathrm i}=300\,\mathrm{MeV}、Mc2=1 GeVMc^2=1\,\mathrm{GeV}、θ=90∘\theta=90^\circ(cos⁡θ=0\cos\theta=0)を入れると

Ef=3001+0.3≃2.3×102 MeV.E_{\mathrm f}=\frac{300}{1+0.3}\simeq2.3\times10^{2}\,\mathrm{MeV}.

答えは Ef≈230 MeVE_{\mathrm f}\approx230\,\mathrm{MeV} で、70 MeV70\,\mathrm{MeV} 程度が陽子の反跳に移ります。

ρ(r)\rho(r) は球対称で ∫ρ d3r=1\int\rho\,d^3r=1 に規格化されています。qq が小さい領域で式 (1) の指数関数を展開します。

F(q2)=∫ρ(r)[1+iq⃗⋅r⃗ℏ−(q⃗⋅r⃗)22ℏ2+⋯ ]d3r.F(q^2)=\int\rho(r)\left[1+i\frac{\vec q\cdot\vec r}{\hbar}-\frac{(\vec q\cdot\vec r)^2}{2\hbar^2}+\cdots\right]d^3r .

第1項は規格化条件から 11。第2項は ρ\rho が球対称なので r⃗\vec r の 1 次の角度積分が消えて 00。第3項は q⃗\vec q の向きを極軸にとると

∫ρ (q⃗⋅r⃗)2d3r=q2∫ρ r2cos⁡2θ d3r=q23∫r2ρ(r) d3r=q2⟨r2⟩3\int\rho\,(\vec q\cdot\vec r)^2d^3r=q^2\int\rho\,r^2\cos^2\theta\,d^3r=\frac{q^2}{3}\int r^2\rho(r)\,d^3r=\frac{q^2\langle r^2\rangle}{3}

です(⟨cos⁡2θ⟩=1/3\langle\cos^2\theta\rangle=1/3)。よって

F(q2)≃1−q2⟨r2⟩6ℏ2,∣F(q2)∣2≃1−q2⟨r2⟩3ℏ2F(q^2)\simeq1-\frac{q^2\langle r^2\rangle}{6\hbar^2},\qquad \left|F(q^2)\right|^2\simeq1-\frac{q^2\langle r^2\rangle}{3\hbar^2}

となり、∣F∣2|F|^2 は q2q^2 の 1 次で直線的に減少します。その傾きを

s≡d∣F(q2)∣2d(q2)=−⟨r2⟩3ℏ2s\equiv\frac{d|F(q^2)|^2}{d(q^2)}=-\frac{\langle r^2\rangle}{3\hbar^2}

と読み取れば

⟨r2⟩=ℏ3 ∣s∣\sqrt{\langle r^2\rangle}=\hbar\sqrt{3\,\lvert s\rvert}

として荷電半径が求まります。q2→0q^2\to0 で ∣F∣2→1|F|^2\to1(点電荷の断面積に一致)という切片も、この展開の帰結として図で確認できます。

図1 の q2q^2 が小さい側の点は、およそ (q2,∣F∣2)=(0,1.00),(0.003,0.87),(0.008,0.75),(0.012,0.62)(q^2,|F|^2)=(0,1.00),(0.003,0.87),(0.008,0.75),(0.012,0.62) です(q2q^2 の単位は (GeV/c)2(\mathrm{GeV}/c)^2)。この直線部分の傾きは

∣s∣≃1.00−0.620.012≃3×101 (GeV/c)−2\lvert s\rvert\simeq\frac{1.00-0.62}{0.012}\simeq3\times10^{1}\,(\mathrm{GeV}/c)^{-2}

です。ℏ2\hbar^2 を (GeV/c)−2(\mathrm{GeV}/c)^{-2} と組ませると ℏ2/(GeV/c)2=(ℏc)2/GeV2=(0.2 GeV⋅fm)2/GeV2=0.04 fm2\hbar^2/(\mathrm{GeV}/c)^2=(\hbar c)^2/\mathrm{GeV}^2=(0.2\,\mathrm{GeV\cdot fm})^2/\mathrm{GeV}^2=0.04\,\mathrm{fm^2} なので

⟨r2⟩=3ℏ2∣s∣≃3×30×0.04=3.6 fm2,⟨r2⟩≃1.9 fm≈2 fm.\langle r^2\rangle=3\hbar^2\lvert s\rvert\simeq3\times30\times0.04=3.6\,\mathrm{fm^2}, \qquad \sqrt{\langle r^2\rangle}\simeq1.9\,\mathrm{fm}\approx2\,\mathrm{fm}.

答えは重陽子の荷電半径 ⟨r2⟩≈2 fm\sqrt{\langle r^2\rangle}\approx2\,\mathrm{fm} です。核子 1 個の半径 0.8 fm0.8\,\mathrm{fm} 程度より大きく、重陽子が結合が弱くて広がった 2 核子系であることと整合します。

第6問 交流回路とダイオード回路

Section titled “第6問 交流回路とダイオード回路”

交流回路の設問です。α\alpha は入力電圧 V1V_1 の振幅に対する出力電圧 V2V_2 の振幅の比、ω=2πf\omega=2\pi f とします。数値は有効数字 1 桁の概算で求めます。

図1 は RR と CC の直列回路で、CC の両端が出力です。分圧の式から

α=∣1/(iωC)R+1/(iωC)∣=∣11+iωRC∣=11+(ωRC)2.\alpha=\left|\frac{1/(i\omega C)}{R+1/(i\omega C)}\right|=\left|\frac{1}{1+i\omega RC}\right| =\frac{1}{\sqrt{1+(\omega RC)^2}} .

漸近線は ωRC≪1\omega RC\ll1 で α=1\alpha=1、ωRC≫1\omega RC\gg1 で α=1/(ωRC)\alpha=1/(\omega RC) です。両者が交わるのは ωRC=1\omega RC=1、つまり

fc=12πRCf_{\mathrm c}=\frac{1}{2\pi RC}

で、そこでの実測値は α=1/2≃0.7\alpha=1/\sqrt2\simeq0.7 になります。図2 の 2 本の破線の交点は f≃2×103 Hzf\simeq2\times10^{3}\,\mathrm{Hz}(そこで実線が約 0.7 を通ることも確認できます)なので fc≃2×103 Hzf_{\mathrm c}\simeq2\times10^{3}\,\mathrm{Hz}、したがって R=1 kΩR=1\,\mathrm{k}\Omega に対して

C=12πfcR=12π×2×103×103≃8×10−8 FC=\frac{1}{2\pi f_{\mathrm c}R}=\frac{1}{2\pi\times2\times10^{3}\times10^{3}}\simeq8\times10^{-8}\,\mathrm{F}

です。答えは C≈8×10−8 FC\approx8\times10^{-8}\,\mathrm{F}(約 0.08 μF0.08\,\mu\mathrm F)です。

抵抗 1 個とコイル 1 個で同じ低周波透過特性を作るには、直列側にコイル、並列(出力)側に抵抗を置きます。すなわち、交流源の一方の端子からコイル LL を直列に通して出力の上側端子へ、そこから抵抗 RR を通して交流源のもう一方の端子(出力の下側端子)へ戻す回路で、RR の両端を V2V_2 とします。図1 の RR をコイルに、CC を抵抗に置き換えた形です。このとき

α=∣RR+iωL∣=11+(ωL/R)2\alpha=\left|\frac{R}{R+i\omega L}\right|=\frac{1}{\sqrt{1+(\omega L/R)^2}}

となり、ωRC\omega RC が ωL/R\omega L/R に置き換わっただけで設問1と同一の関数形です。折れ点は ωL/R=1\omega L/R=1、すなわち fc=R/(2πL)f_{\mathrm c}=R/(2\pi L) なので、図2 の fc≃2×103 Hzf_{\mathrm c}\simeq2\times10^{3}\,\mathrm{Hz} と R=100 ΩR=100\,\Omega から

L=R2πfc=1002π×2×103≃8×10−3 HL=\frac{R}{2\pi f_{\mathrm c}}=\frac{100}{2\pi\times2\times10^{3}}\simeq8\times10^{-3}\,\mathrm{H}

です。答えは L≈8×10−3 HL\approx8\times10^{-3}\,\mathrm{H}(約 8 mH)です。

図3 は、RR を直列に通した先に LL と CC が並列に入り、その並列部分の電圧が V2V_2 です。並列部のインピーダンスは

Z=(1iωL+iωC)−1=iωL1−ω2LCZ=\left(\frac{1}{i\omega L}+i\omega C\right)^{-1}=\frac{i\omega L}{1-\omega^2LC}

なので

α=∣ZR+Z∣=ωLR2(1−ω2LC)2+ω2L2.\alpha=\left|\frac{Z}{R+Z}\right| =\frac{\omega L}{\sqrt{R^2\left(1-\omega^2LC\right)^2+\omega^2L^2}} .

ω0=1/LC\omega_0=1/\sqrt{LC} で分母の第1項が消えて α=1\alpha=1 の鋭い共鳴ピークになります。漸近形は ω≪ω0\omega\ll\omega_0 で α≃ωL/R\alpha\simeq\omega L/R(傾き +1+1)、ω≫ω0\omega\gg\omega_0 で α≃1/(ωRC)\alpha\simeq1/(\omega RC)(傾き −1-1)で、両漸近線の交点は ω=ω0\omega=\omega_0 かつ

αasym=ω0LR=1RLC\alpha_{\mathrm{asym}}=\frac{\omega_0L}{R}=\frac{1}{R}\sqrt{\frac{L}{C}}

です。図4 からピーク位置 f0≃1×105 Hzf_0\simeq1\times10^{5}\,\mathrm{Hz}(実線のピーク値は 1)、破線の交点の高さ αasym≃0.2\alpha_{\mathrm{asym}}\simeq0.2 と読めます。C=10 nF=1×10−8 FC=10\,\mathrm{nF}=1\times10^{-8}\,\mathrm{F} として

L=1(2πf0)2C=1(6.28×105)2×10−8≃2.5×10−4 H,L=\frac{1}{(2\pi f_0)^2C}=\frac{1}{(6.28\times10^{5})^2\times10^{-8}}\simeq2.5\times10^{-4}\,\mathrm{H}, R=1αasymLC=10.22.5×10−410−8≃8×102 Ω.R=\frac{1}{\alpha_{\mathrm{asym}}}\sqrt{\frac{L}{C}}=\frac{1}{0.2}\sqrt{\frac{2.5\times10^{-4}}{10^{-8}}}\simeq8\times10^{2}\,\Omega .

答えは L≈2.5×10−4 HL\approx2.5\times10^{-4}\,\mathrm{H}(約 0.25 mH)、R≈8×102 ΩR\approx8\times10^{2}\,\Omega です。左端 f=104 Hzf=10^{4}\,\mathrm{Hz} で α≃0.02=0.2/10\alpha\simeq0.02=0.2/10 となっているのも、α∝ω\alpha\propto\omega の漸近形と整合しています。

図5 では、交流源の 2 つの端子が B と D に出ており、出力側は A が V2V_2 の高電位側、C が基準です(V1=VB−VDV_1=V_{\mathrm B}-V_{\mathrm D}、V2=VA−VCV_2=V_{\mathrm A}-V_{\mathrm C})。V2=∣V1∣V_2=|V_1| とするには、A を 2 個のダイオードのカソード側、C を 2 個のダイオードのアノード側にした全波整流ブリッジを組みます。すなわち、順方向を B から A へ、D から A へ、C から B へ、C から D へ向けて、ダイオードを 1 個ずつ入れます(記号の三角形の向きがそれぞれ A、A、B、D を向く)。

動作を確認します。V1>0V_1>0(B が高電位)のときは、B から A へ入り RR を通って C へ抜け、C から D へ戻る経路だけが順方向で導通します。このとき VA=VBV_{\mathrm A}=V_{\mathrm B}、VC=VDV_{\mathrm C}=V_{\mathrm D} なので V2=VB−VD=V1>0V_2=V_{\mathrm B}-V_{\mathrm D}=V_1>0 です。V1<0V_1<0(D が高電位)のときは、D から A へ入り RR を通って C へ抜け、C から B へ戻る経路が導通し、V2=VD−VB=−V1>0V_2=V_{\mathrm D}-V_{\mathrm B}=-V_1>0 です。どちらの半周期でも A が高電位側になるので V2=∣V1∣V_2=|V_1| が得られます。

A–C 間にコンデンサを入れると、V2V_2 はピーク値 Vmax⁡=12 VV_{\max}=12\,\mathrm{V}(理想ダイオードなので順方向降下なし)まで充電され、その後は RR を通じて放電します。全波整流なので充電される山は 1 周期に 2 回、放電時間は

T2=12f=12×50=1×10−2 s\frac{T}{2}=\frac{1}{2f}=\frac{1}{2\times50}=1\times10^{-2}\,\mathrm{s}

です。時定数は RC=500×3×10−3=1.5 sRC=500\times3\times10^{-3}=1.5\,\mathrm{s} で放電時間よりはるかに長いので、指数関数を 1 次まで展開して

ΔV=Vmax⁡(1−e−T/2RC)≃Vmax⁡T2RC=12×1×10−21.5≃8×10−2 V.\Delta V=V_{\max}\left(1-e^{-\frac{T/2}{RC}}\right)\simeq V_{\max}\frac{T}{2RC} =12\times\frac{1\times10^{-2}}{1.5}\simeq8\times10^{-2}\,\mathrm{V}.

答えは ΔV≈0.08 V\Delta V\approx0.08\,\mathrm{V} です。12 V12\,\mathrm V に対して 1 % 弱のリプルに収まっています。

図8 は交流源、ダイオード、コンデンサの直列閉回路で、コンデンサの両端が V2V_2 です。コンデンサの容量が十分大きく V2V_2 が直流とみなせるなら、定常状態でコンデンサに流れ込む電荷の 1 周期平均は 00、つまりダイオードを流れる電流の時間平均が 00 でなければなりません。

ダイオードに加わる電圧は V=V1−V2V=V_1-V_2 で、V1=V0sin⁡ωtV_1=V_0\sin\omega t(V0=10 mVV_0=10\,\mathrm{mV})です。与えられた特性 I=gV(1+V2VT)I=gV\left(1+\dfrac{V}{2V_{\mathrm T}}\right) の時間平均をとると、⟨V⟩=−V2\langle V\rangle=-V_2、⟨V2⟩=V022+V22\langle V^2\rangle=\dfrac{V_0^2}{2}+V_2^2 なので

⟨I⟩=g[⟨V⟩+⟨V2⟩2VT]=g[−V2+V02/2+V222VT]=0.\langle I\rangle=g\left[\langle V\rangle+\frac{\langle V^2\rangle}{2V_{\mathrm T}}\right] =g\left[-V_2+\frac{V_0^2/2+V_2^2}{2V_{\mathrm T}}\right]=0 .

V2V_2 は V0V_0 の 2 次の微小量なので V22V_2^2 を落とすと

V2≃V024VT=(10 mV)24×25 mV=1 mV.V_2\simeq\frac{V_0^2}{4V_{\mathrm T}}=\frac{(10\,\mathrm{mV})^2}{4\times25\,\mathrm{mV}}=1\,\mathrm{mV}.

答えは V2≈1 mVV_2\approx1\,\mathrm{mV} です(V22V_2^2 を残して解くと 1.02 mV1.02\,\mathrm{mV} で、近似は妥当です)。線形項の係数 gg は打ち消えて答えに入りません。直流電圧を生むのは特性の 2 次項だけであり、その大きさが入力振幅の 2 乗に比例する(二乗検波になる)ことがこの回路の特徴です。

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

執筆: Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

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

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。