4 時間で 4 問という分量の試験です。前半 3 問（量子力学・統計力学・電磁気学）は素直な設定から始まって、最後に一段深いところへ連れて行かれる形になっています。後半 3 問は選択で、第4問と第5問はほとんど実験屋の日常業務（熱容量測定、光検出器の雑音見積り）、第6問は加速器実験の入口です。手が止まりやすいのは第1問設問8以降の確率の勘定、第2問設問8の熱力学的な組み立て、第4問設問7の表からの読み取りでしょう。
問題 分野 主題 第1問 量子力学 二状態系の測定とベルの不等式（Mermin 版） 第2問 統計力学 古典理想気体と自由フェルミ粒子系 第3問 電磁気学 同軸ケーブルの TEM モードと群速度 第4問 熱力学・物性物理 緩和法と交流法による熱容量測定 第5問 電気回路・統計力学 抵抗の熱雑音と並列共振型光検出器 第6問 相対論・素粒子 反陽子生成のしきい値と飛行時間測定
第1問から第3問が必答で、第4問・第5問・第6問から 1 問を選択する形式ですが、ここでは全 6 問の解答を載せます。
観測量は
σ z = ( 1 0 0 − 1 ) , σ x = ( 0 1 1 0 ) , σ ( θ ) = ( cos θ ) σ z + ( sin θ ) σ x \sigma_z=\begin{pmatrix}1&0\\0&-1\end{pmatrix},\qquad
\sigma_x=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad
\sigma(\theta)=(\cos\theta)\,\sigma_z+(\sin\theta)\,\sigma_x σ z = ( 1 0 0 − 1 ) , σ x = ( 0 1 1 0 ) , σ ( θ ) = ( cos θ ) σ z + ( sin θ ) σ x で与えられ、状態ベクトルは ∣ ↑ ⟩ = ( 1 , 0 ) T |\uparrow\rangle=(1,0)^{T} ∣ ↑ ⟩ = ( 1 , 0 ) T ∣ ↓ ⟩ = ( 0 , 1 ) T |\downarrow\rangle=(0,1)^{T} ∣ ↓ ⟩ = ( 0 , 1 ) T σ ( θ ) \sigma(\theta) σ ( θ ) z z z x x x z z z θ \theta θ
まず固有値を確定させます。σ z σ x + σ x σ z = 0 \sigma_z\sigma_x+\sigma_x\sigma_z=0 σ z σ x + σ x σ z = 0 σ z 2 = σ x 2 = I \sigma_z^2=\sigma_x^2=I σ z 2 = σ x 2 = I
σ ( θ ) 2 = cos 2 θ σ z 2 + sin 2 θ σ x 2 + cos θ sin θ ( σ z σ x + σ x σ z ) = I \sigma(\theta)^2=\cos^2\theta\,\sigma_z^2+\sin^2\theta\,\sigma_x^2+\cos\theta\sin\theta\,(\sigma_z\sigma_x+\sigma_x\sigma_z)=I σ ( θ ) 2 = cos 2 θ σ z 2 + sin 2 θ σ x 2 + cos θ sin θ ( σ z σ x + σ x σ z ) = I なので、σ ( θ ) \sigma(\theta) σ ( θ ) ± 1 \pm1 ± 1
∣ θ , + ⟩ = cos θ 2 ∣ ↑ ⟩ + sin θ 2 ∣ ↓ ⟩ , ∣ θ , − ⟩ = − sin θ 2 ∣ ↑ ⟩ + cos θ 2 ∣ ↓ ⟩ |\theta,+\rangle=\cos\frac{\theta}{2}|\uparrow\rangle+\sin\frac{\theta}{2}|\downarrow\rangle,\qquad
|\theta,-\rangle=-\sin\frac{\theta}{2}|\uparrow\rangle+\cos\frac{\theta}{2}|\downarrow\rangle ∣ θ , + ⟩ = cos 2 θ ∣ ↑ ⟩ + sin 2 θ ∣ ↓ ⟩ , ∣ θ , − ⟩ = − sin 2 θ ∣ ↑ ⟩ + cos 2 θ ∣ ↓ ⟩ が固有値 + 1 , − 1 +1,-1 + 1 , − 1 θ = 0 \theta=0 θ = 0 σ z \sigma_z σ z θ = π / 2 \theta=\pi/2 θ = π /2 σ x \sigma_x σ x
後半では二状態系 A , B A,B A , B
∣ Ψ ⟩ = 1 2 ( ∣ ↑ ⟩ A ∣ ↓ ⟩ B − ∣ ↓ ⟩ A ∣ ↑ ⟩ B ) |\Psi\rangle=\frac{1}{\sqrt2}\Big(|\uparrow\rangle_A|\downarrow\rangle_B-|\downarrow\rangle_A|\uparrow\rangle_B\Big) ∣Ψ ⟩ = 2 1 ( ∣ ↑ ⟩ A ∣ ↓ ⟩ B − ∣ ↓ ⟩ A ∣ ↑ ⟩ B ) です。
∣ ↑ ⟩ |\uparrow\rangle ∣ ↑ ⟩ σ z \sigma_z σ z + 1 +1 + 1 σ z \sigma_z σ z ± 1 \pm1 ± 1 − 1 -1 − 1 ∣ ⟨ ↓ ∣ ↑ ⟩ ∣ 2 = 0 |\langle\downarrow|\uparrow\rangle|^2=0 ∣ ⟨ ↓ ∣ ↑ ⟩ ∣ 2 = 0
s z = + 1 s_z=+1 s z = + 1 だけをとります。期待値は ⟨ σ z ⟩ = ⟨ ↑ ∣ σ z ∣ ↑ ⟩ = + 1 \langle\sigma_z\rangle=\langle\uparrow|\sigma_z|\uparrow\rangle=+1 ⟨ σ z ⟩ = ⟨ ↑ ∣ σ z ∣ ↑ ⟩ = + 1
σ x \sigma_x σ x ∣ x , ± ⟩ = ( ∣ ↑ ⟩ ± ∣ ↓ ⟩ ) / 2 |x,\pm\rangle=(|\uparrow\rangle\pm|\downarrow\rangle)/\sqrt2 ∣ x , ± ⟩ = ( ∣ ↑ ⟩ ± ∣ ↓ ⟩) / 2 ∣ ↑ ⟩ = ( ∣ x , + ⟩ + ∣ x , − ⟩ ) / 2 |\uparrow\rangle=(|x,+\rangle+|x,-\rangle)/\sqrt2 ∣ ↑ ⟩ = ( ∣ x , + ⟩ + ∣ x , − ⟩) / 2
s x = + 1 ( 確率 1 / 2 ) , s x = − 1 ( 確率 1 / 2 ) s_x=+1\ (\text{確率}\ 1/2),\qquad s_x=-1\ (\text{確率}\ 1/2) s x = + 1 ( 確率 1/2 ) , s x = − 1 ( 確率 1/2 ) の 2 通りです。期待値は
⟨ σ x ⟩ = ⟨ ↑ ∣ σ x ∣ ↑ ⟩ = ( 1 0 ) ( 0 1 1 0 ) ( 1 0 ) = 0 \langle\sigma_x\rangle=\langle\uparrow|\sigma_x|\uparrow\rangle=(1\ \ 0)\begin{pmatrix}0&1\\1&0\end{pmatrix}\begin{pmatrix}1\\0\end{pmatrix}=0 ⟨ σ x ⟩ = ⟨ ↑ ∣ σ x ∣ ↑ ⟩ = ( 1 0 ) ( 0 1 1 0 ) ( 1 0 ) = 0 です。答えは s x = ± 1 s_x=\pm1 s x = ± 1 0 0 0
設定で求めた固有ベクトルを使うと ⟨ θ , + ∣ ↑ ⟩ = cos ( θ / 2 ) \langle\theta,+|\uparrow\rangle=\cos(\theta/2) ⟨ θ , + ∣ ↑ ⟩ = cos ( θ /2 ) ⟨ θ , − ∣ ↑ ⟩ = − sin ( θ / 2 ) \langle\theta,-|\uparrow\rangle=-\sin(\theta/2) ⟨ θ , − ∣ ↑ ⟩ = − sin ( θ /2 )
s θ = + 1 ( 確率 cos 2 θ 2 ) , s θ = − 1 ( 確率 sin 2 θ 2 ) s_\theta=+1\ \left(\text{確率}\ \cos^2\frac{\theta}{2}\right),\qquad
s_\theta=-1\ \left(\text{確率}\ \sin^2\frac{\theta}{2}\right) s θ = + 1 ( 確率 cos 2 2 θ ) , s θ = − 1 ( 確率 sin 2 2 θ ) の 2 通りです。期待値は
⟨ σ ( θ ) ⟩ = cos 2 θ 2 − sin 2 θ 2 = cos θ \langle\sigma(\theta)\rangle=\cos^2\frac{\theta}{2}-\sin^2\frac{\theta}{2}=\cos\theta ⟨ σ ( θ )⟩ = cos 2 2 θ − sin 2 2 θ = cos θ で、これは ⟨ ↑ ∣ σ ( θ ) ∣ ↑ ⟩ = cos θ \langle\uparrow|\sigma(\theta)|\uparrow\rangle=\cos\theta ⟨ ↑ ∣ σ ( θ ) ∣ ↑ ⟩ = cos θ θ = 0 \theta=0 θ = 0 1 1 1 θ = π / 2 \theta=\pi/2 θ = π /2 0 0 0
∣ Ψ ⟩ |\Psi\rangle ∣Ψ ⟩ ∣ ↑ ⟩ A ∣ ↑ ⟩ B |\uparrow\rangle_A|\uparrow\rangle_B ∣ ↑ ⟩ A ∣ ↑ ⟩ B ∣ ↓ ⟩ A ∣ ↓ ⟩ B |\downarrow\rangle_A|\downarrow\rangle_B ∣ ↓ ⟩ A ∣ ↓ ⟩ B ∣ ↑ ⟩ A ∣ ↓ ⟩ B |\uparrow\rangle_A|\downarrow\rangle_B ∣ ↑ ⟩ A ∣ ↓ ⟩ B ∣ ↓ ⟩ A ∣ ↑ ⟩ B |\downarrow\rangle_A|\uparrow\rangle_B ∣ ↓ ⟩ A ∣ ↑ ⟩ B σ z A σ z B \sigma_z^A\sigma_z^B σ z A σ z B
( s z A , s z B ) = ( + 1 , − 1 ) または ( − 1 , + 1 ) (s_z^A,s_z^B)=(+1,-1)\ \ \text{または}\ \ (-1,+1) ( s z A , s z B ) = ( + 1 , − 1 ) または ( − 1 , + 1 ) の 2 通りで、それぞれ確率 1 / 2 1/2 1/2 − 1 -1 − 1
⟨ s z A s z B ⟩ = − 1 \langle s_z^A s_z^B\rangle=-1 ⟨ s z A s z B ⟩ = − 1 です。演算子で書けば σ z A σ z B ∣ Ψ ⟩ = − ∣ Ψ ⟩ \sigma_z^A\sigma_z^B|\Psi\rangle=-|\Psi\rangle σ z A σ z B ∣Ψ ⟩ = − ∣Ψ ⟩ ∣ Ψ ⟩ |\Psi\rangle ∣Ψ ⟩ − 1 -1 − 1
σ x A σ x B \sigma_x^A\sigma_x^B σ x A σ x B ∣ ↑ ⟩ A ∣ ↓ ⟩ B |\uparrow\rangle_A|\downarrow\rangle_B ∣ ↑ ⟩ A ∣ ↓ ⟩ B ∣ ↓ ⟩ A ∣ ↑ ⟩ B |\downarrow\rangle_A|\uparrow\rangle_B ∣ ↓ ⟩ A ∣ ↑ ⟩ B
σ x A σ x B ∣ Ψ ⟩ = 1 2 ( ∣ ↓ ⟩ A ∣ ↑ ⟩ B − ∣ ↑ ⟩ A ∣ ↓ ⟩ B ) = − ∣ Ψ ⟩ \sigma_x^A\sigma_x^B|\Psi\rangle=\frac{1}{\sqrt2}\Big(|\downarrow\rangle_A|\uparrow\rangle_B-|\uparrow\rangle_A|\downarrow\rangle_B\Big)=-|\Psi\rangle σ x A σ x B ∣Ψ ⟩ = 2 1 ( ∣ ↓ ⟩ A ∣ ↑ ⟩ B − ∣ ↑ ⟩ A ∣ ↓ ⟩ B ) = − ∣Ψ ⟩ となり、∣ Ψ ⟩ |\Psi\rangle ∣Ψ ⟩ σ x A σ x B \sigma_x^A\sigma_x^B σ x A σ x B − 1 -1 − 1 ∣ x , ± ⟩ |x,\pm\rangle ∣ x , ± ⟩
∣ Ψ ⟩ = − 1 2 ( ∣ x , + ⟩ A ∣ x , − ⟩ B − ∣ x , − ⟩ A ∣ x , + ⟩ B ) |\Psi\rangle=-\frac{1}{\sqrt2}\Big(|x,+\rangle_A|x,-\rangle_B-|x,-\rangle_A|x,+\rangle_B\Big) ∣Ψ ⟩ = − 2 1 ( ∣ x , + ⟩ A ∣ x , − ⟩ B − ∣ x , − ⟩ A ∣ x , + ⟩ B ) なので、取りうる組は
( s x A , s x B ) = ( + 1 , − 1 ) または ( − 1 , + 1 ) (s_x^A,s_x^B)=(+1,-1)\ \ \text{または}\ \ (-1,+1) ( s x A , s x B ) = ( + 1 , − 1 ) または ( − 1 , + 1 ) の 2 通り（各確率 1 / 2 1/2 1/2 ⟨ s x A s x B ⟩ = − 1 \langle s_x^A s_x^B\rangle=-1 ⟨ s x A s x B ⟩ = − 1 z z z x x x
c = cos ( θ / 2 ) c=\cos(\theta/2) c = cos ( θ /2 ) s = sin ( θ / 2 ) s=\sin(\theta/2) s = sin ( θ /2 ) ∣ θ , ± ⟩ |\theta,\pm\rangle ∣ θ , ± ⟩ ∣ ↑ ⟩ = c ∣ θ , + ⟩ − s ∣ θ , − ⟩ |\uparrow\rangle=c|\theta,+\rangle-s|\theta,-\rangle ∣ ↑ ⟩ = c ∣ θ , + ⟩ − s ∣ θ , − ⟩ ∣ ↓ ⟩ = s ∣ θ , + ⟩ + c ∣ θ , − ⟩ |\downarrow\rangle=s|\theta,+\rangle+c|\theta,-\rangle ∣ ↓ ⟩ = s ∣ θ , + ⟩ + c ∣ θ , − ⟩ ∣ Ψ ⟩ |\Psi\rangle ∣Ψ ⟩
∣ ↑ ⟩ A ∣ ↓ ⟩ B = c s ∣ + + ⟩ + c 2 ∣ + − ⟩ − s 2 ∣ − + ⟩ − s c ∣ − − ⟩ , ∣ ↓ ⟩ A ∣ ↑ ⟩ B = s c ∣ + + ⟩ − s 2 ∣ + − ⟩ + c 2 ∣ − + ⟩ − c s ∣ − − ⟩ \begin{aligned}
|\uparrow\rangle_A|\downarrow\rangle_B&=cs\,|{+}{+}\rangle+c^2|{+}{-}\rangle-s^2|{-}{+}\rangle-sc\,|{-}{-}\rangle,\\
|\downarrow\rangle_A|\uparrow\rangle_B&=sc\,|{+}{+}\rangle-s^2|{+}{-}\rangle+c^2|{-}{+}\rangle-cs\,|{-}{-}\rangle
\end{aligned} ∣ ↑ ⟩ A ∣ ↓ ⟩ B ∣ ↓ ⟩ A ∣ ↑ ⟩ B = cs ∣ + + ⟩ + c 2 ∣ + − ⟩ − s 2 ∣ − + ⟩ − sc ∣ − − ⟩ , = sc ∣ + + ⟩ − s 2 ∣ + − ⟩ + c 2 ∣ − + ⟩ − cs ∣ − − ⟩ （∣ ± ± ⟩ |{\pm}{\pm}\rangle ∣ ± ± ⟩ ∣ θ , ± ⟩ A ∣ θ , ± ⟩ B |\theta,\pm\rangle_A|\theta,\pm\rangle_B ∣ θ , ± ⟩ A ∣ θ , ± ⟩ B ∣ + + ⟩ |{+}{+}\rangle ∣ + + ⟩ ∣ − − ⟩ |{-}{-}\rangle ∣ − − ⟩ c 2 + s 2 = 1 c^2+s^2=1 c 2 + s 2 = 1
∣ Ψ ⟩ = 1 2 ( ∣ θ , + ⟩ A ∣ θ , − ⟩ B − ∣ θ , − ⟩ A ∣ θ , + ⟩ B ) |\Psi\rangle=\frac{1}{\sqrt2}\Big(|\theta,+\rangle_A|\theta,-\rangle_B-|\theta,-\rangle_A|\theta,+\rangle_B\Big) ∣Ψ ⟩ = 2 1 ( ∣ θ , + ⟩ A ∣ θ , − ⟩ B − ∣ θ , − ⟩ A ∣ θ , + ⟩ B ) です。θ \theta θ ∣ θ , + ⟩ A ∣ θ , + ⟩ B |\theta,+\rangle_A|\theta,+\rangle_B ∣ θ , + ⟩ A ∣ θ , + ⟩ B ∣ θ , − ⟩ A ∣ θ , − ⟩ B |\theta,-\rangle_A|\theta,-\rangle_B ∣ θ , − ⟩ A ∣ θ , − ⟩ B 0 0 0 ( s θ A , s θ B ) (s_\theta^A,s_\theta^B) ( s θ A , s θ B ) ( + 1 , − 1 ) (+1,-1) ( + 1 , − 1 ) ( − 1 , + 1 ) (-1,+1) ( − 1 , + 1 ) s θ A = − s θ B = ± 1 s_\theta^A=-s_\theta^B=\pm1 s θ A = − s θ B = ± 1
まず期待値を求めます。σ A ( θ ) σ B ( φ ) \sigma^A(\theta)\sigma^B(\varphi) σ A ( θ ) σ B ( φ ) σ z A σ x B ∣ Ψ ⟩ = ( ∣ ↑ ⟩ A ∣ ↑ ⟩ B + ∣ ↓ ⟩ A ∣ ↓ ⟩ B ) / 2 \sigma_z^A\sigma_x^B|\Psi\rangle=(|\uparrow\rangle_A|\uparrow\rangle_B+|\downarrow\rangle_A|\downarrow\rangle_B)/\sqrt2 σ z A σ x B ∣Ψ ⟩ = ( ∣ ↑ ⟩ A ∣ ↑ ⟩ B + ∣ ↓ ⟩ A ∣ ↓ ⟩ B ) / 2 ∣ Ψ ⟩ |\Psi\rangle ∣Ψ ⟩ ⟨ σ z A σ x B ⟩ = 0 \langle\sigma_z^A\sigma_x^B\rangle=0 ⟨ σ z A σ x B ⟩ = 0 ⟨ σ x A σ z B ⟩ = 0 \langle\sigma_x^A\sigma_z^B\rangle=0 ⟨ σ x A σ z B ⟩ = 0
⟨ σ A ( θ ) σ B ( φ ) ⟩ = − cos θ cos φ − sin θ sin φ = − cos ( θ − φ ) \langle\sigma^A(\theta)\sigma^B(\varphi)\rangle=-\cos\theta\cos\varphi-\sin\theta\sin\varphi=-\cos(\theta-\varphi) ⟨ σ A ( θ ) σ B ( φ )⟩ = − cos θ cos φ − sin θ sin φ = − cos ( θ − φ ) となります。つまり ⟨ s θ A s φ B ⟩ = − cos ( θ − φ ) \langle s_\theta^A s_\varphi^B\rangle=-\cos(\theta-\varphi) ⟨ s θ A s φ B ⟩ = − cos ( θ − φ )
取りうる組を出すために、設問6 の表式で B B B φ \varphi φ α ≡ ( θ − φ ) / 2 \alpha\equiv(\theta-\varphi)/2 α ≡ ( θ − φ ) /2 ∣ θ , + ⟩ = cos α ∣ φ , + ⟩ + sin α ∣ φ , − ⟩ |\theta,+\rangle=\cos\alpha\,|\varphi,+\rangle+\sin\alpha\,|\varphi,-\rangle ∣ θ , + ⟩ = cos α ∣ φ , + ⟩ + sin α ∣ φ , − ⟩ ∣ θ , − ⟩ = − sin α ∣ φ , + ⟩ + cos α ∣ φ , − ⟩ |\theta,-\rangle=-\sin\alpha\,|\varphi,+\rangle+\cos\alpha\,|\varphi,-\rangle ∣ θ , − ⟩ = − sin α ∣ φ , + ⟩ + cos α ∣ φ , − ⟩
∣ Ψ ⟩ = 1 2 [ − sin α ∣ θ , + ⟩ A ∣ φ , + ⟩ B + cos α ∣ θ , + ⟩ A ∣ φ , − ⟩ B − cos α ∣ θ , − ⟩ A ∣ φ , + ⟩ B − sin α ∣ θ , − ⟩ A ∣ φ , − ⟩ B ] |\Psi\rangle=\frac{1}{\sqrt2}\Big[-\sin\alpha\,|\theta,+\rangle_A|\varphi,+\rangle_B+\cos\alpha\,|\theta,+\rangle_A|\varphi,-\rangle_B-\cos\alpha\,|\theta,-\rangle_A|\varphi,+\rangle_B-\sin\alpha\,|\theta,-\rangle_A|\varphi,-\rangle_B\Big] ∣Ψ ⟩ = 2 1 [ − sin α ∣ θ , + ⟩ A ∣ φ , + ⟩ B + cos α ∣ θ , + ⟩ A ∣ φ , − ⟩ B − cos α ∣ θ , − ⟩ A ∣ φ , + ⟩ B − sin α ∣ θ , − ⟩ A ∣ φ , − ⟩ B ] です。よって θ ≠ φ \theta\neq\varphi θ = φ m o d 2 π \bmod\ 2\pi mod 2 π
( s θ A , s φ B ) = ( + 1 , + 1 ) , ( + 1 , − 1 ) , ( − 1 , + 1 ) , ( − 1 , − 1 ) (s_\theta^A,s_\varphi^B)=(+1,+1),\ (+1,-1),\ (-1,+1),\ (-1,-1) ( s θ A , s φ B ) = ( + 1 , + 1 ) , ( + 1 , − 1 ) , ( − 1 , + 1 ) , ( − 1 , − 1 ) がすべて起こり、確率は同符号が各 1 2 sin 2 α \tfrac12\sin^2\alpha 2 1 sin 2 α 1 2 cos 2 α \tfrac12\cos^2\alpha 2 1 cos 2 α
sin 2 α − cos 2 α = − cos 2 α = − cos ( θ − φ ) \sin^2\alpha-\cos^2\alpha=-\cos2\alpha=-\cos(\theta-\varphi) sin 2 α − cos 2 α = − cos 2 α = − cos ( θ − φ ) で、上の演算子計算と一致します。θ = φ \theta=\varphi θ = φ 0 0 0
θ , φ \theta,\varphi θ , φ { 0 ∘ , 120 ∘ , 240 ∘ } \{0^\circ,120^\circ,240^\circ\} { 0 ∘ , 12 0 ∘ , 24 0 ∘ } 1 / 3 1/3 1/3 θ = φ \theta=\varphi θ = φ θ ≠ φ \theta\neq\varphi θ = φ
( s A , s B ) = ( + 1 , + 1 ) , ( + 1 , − 1 ) , ( − 1 , + 1 ) , ( − 1 , − 1 ) (s^A,s^B)=(+1,+1),\ (+1,-1),\ (-1,+1),\ (-1,-1) ( s A , s B ) = ( + 1 , + 1 ) , ( + 1 , − 1 ) , ( − 1 , + 1 ) , ( − 1 , − 1 ) の 4 通りです。
期待値は設問7 の結果を 9 通りについて平均すればよく、θ − φ \theta-\varphi θ − φ 0 ∘ , 120 ∘ , 240 ∘ 0^\circ,120^\circ,240^\circ 0 ∘ , 12 0 ∘ , 24 0 ∘
⟨ s A s B ⟩ = − 1 9 ∑ θ ∑ φ cos ( θ − φ ) = − 1 9 [ 3 cos 0 ∘ + 3 cos 120 ∘ + 3 cos 240 ∘ ] = − 1 9 [ 3 − 3 2 − 3 2 ] = 0 \langle s^As^B\rangle=-\frac{1}{9}\sum_{\theta}\sum_{\varphi}\cos(\theta-\varphi)
=-\frac{1}{9}\Big[3\cos0^\circ+3\cos120^\circ+3\cos240^\circ\Big]
=-\frac{1}{9}\left[3-\frac32-\frac32\right]=0 ⟨ s A s B ⟩ = − 9 1 θ ∑ φ ∑ cos ( θ − φ ) = − 9 1 [ 3 cos 0 ∘ + 3 cos 12 0 ∘ + 3 cos 24 0 ∘ ] = − 9 1 [ 3 − 2 3 − 2 3 ] = 0 です。見通しよく言えば、cos ( θ − φ ) = n ( θ ) ⋅ n ( φ ) \cos(\theta-\varphi)=\boldsymbol{n}(\theta)\cdot\boldsymbol{n}(\varphi) cos ( θ − φ ) = n ( θ ) ⋅ n ( φ ) n ( θ ) = ( sin θ , cos θ ) \boldsymbol{n}(\theta)=(\sin\theta,\cos\theta) n ( θ ) = ( sin θ , cos θ )
∑ θ ∑ φ n ( θ ) ⋅ n ( φ ) = ( ∑ θ n ( θ ) ) ⋅ ( ∑ φ n ( φ ) ) = 0 \sum_{\theta}\sum_{\varphi}\boldsymbol{n}(\theta)\cdot\boldsymbol{n}(\varphi)=\Big(\sum_\theta\boldsymbol{n}(\theta)\Big)\cdot\Big(\sum_\varphi\boldsymbol{n}(\varphi)\Big)=0 θ ∑ φ ∑ n ( θ ) ⋅ n ( φ ) = ( θ ∑ n ( θ ) ) ⋅ ( φ ∑ n ( φ ) ) = 0 となります。120 ∘ 120^\circ 12 0 ∘ ⟨ s A s B ⟩ = 0 \langle s^As^B\rangle=0 ⟨ s A s B ⟩ = 0
仮説のもとでは、各回の実験の前から ( s 0 ∘ A , s 120 ∘ A , s 240 ∘ A ) (s_{0^\circ}^A,s_{120^\circ}^A,s_{240^\circ}^A) ( s 0 ∘ A , s 12 0 ∘ A , s 24 0 ∘ A ) B B B s φ B = − s φ A s_\varphi^B=-s_\varphi^A s φ B = − s φ A A , B A,B A , B Σ ≡ s 0 ∘ A + s 120 ∘ A + s 240 ∘ A \Sigma\equiv s_{0^\circ}^A+s_{120^\circ}^A+s_{240^\circ}^A Σ ≡ s 0 ∘ A + s 12 0 ∘ A + s 24 0 ∘ A
⟨ s A s B ⟩ = 1 9 ∑ θ ∑ φ s θ A ( − s φ A ) = − 1 9 ( ∑ θ s θ A ) ( ∑ φ s φ A ) = − Σ 2 9 \langle s^As^B\rangle=\frac{1}{9}\sum_{\theta}\sum_{\varphi}s_\theta^A\left(-s_\varphi^A\right)
=-\frac{1}{9}\Big(\sum_\theta s_\theta^A\Big)\Big(\sum_\varphi s_\varphi^A\Big)=-\frac{\Sigma^2}{9} ⟨ s A s B ⟩ = 9 1 θ ∑ φ ∑ s θ A ( − s φ A ) = − 9 1 ( θ ∑ s θ A ) ( φ ∑ s φ A ) = − 9 Σ 2 が成り立ちます。以下これを使います。
(i) ( + 1 , + 1 , + 1 ) (+1,+1,+1) ( + 1 , + 1 , + 1 ) Σ = 3 \Sigma=3 Σ = 3 ⟨ s A s B ⟩ = − 9 / 9 = − 1 \langle s^As^B\rangle=-9/9=-1 ⟨ s A s B ⟩ = − 9/9 = − 1 A A A s A = + 1 s^A=+1 s A = + 1 B B B s B = − 1 s^B=-1 s B = − 1 − 1 -1 − 1
(ii) ( + 1 , + 1 , − 1 ) (+1,+1,-1) ( + 1 , + 1 , − 1 ) Σ = 1 \Sigma=1 Σ = 1 ⟨ s A s B ⟩ = − 1 / 9 \langle s^As^B\rangle=-1/9 ⟨ s A s B ⟩ = − 1/9 ⟨ s A ⟩ = 1 / 3 \langle s^A\rangle=1/3 ⟨ s A ⟩ = 1/3 ⟨ s B ⟩ = − 1 / 3 \langle s^B\rangle=-1/3 ⟨ s B ⟩ = − 1/3 − 1 / 9 -1/9 − 1/9
(iii) Σ \Sigma Σ ± 1 \pm1 ± 1 Σ ≠ 0 \Sigma\neq0 Σ = 0 Σ 2 ≥ 1 \Sigma^2\ge1 Σ 2 ≥ 1
⟨ s A s B ⟩ = − Σ 2 9 ≤ − 1 9 < 0 \langle s^As^B\rangle=-\frac{\Sigma^2}{9}\le-\frac{1}{9}<0 ⟨ s A s B ⟩ = − 9 Σ 2 ≤ − 9 1 < 0 です。8 通りの配位 c c c p c p_c p c ∑ c p c = 1 \sum_c p_c=1 ∑ c p c = 1
⟨ s A s B ⟩ = − 1 9 ∑ c p c Σ c 2 ≤ − 1 9 < 0 \langle s^As^B\rangle=-\frac{1}{9}\sum_c p_c\,\Sigma_c^{\,2}\le-\frac{1}{9}<0 ⟨ s A s B ⟩ = − 9 1 c ∑ p c Σ c 2 ≤ − 9 1 < 0 となり、負であることが示されました。値の範囲は − 1 ≤ ⟨ s A s B ⟩ ≤ − 1 / 9 -1\le\langle s^As^B\rangle\le-1/9 − 1 ≤ ⟨ s A s B ⟩ ≤ − 1/9 0 0 0
質量 m m m N N N V V V T T T
Z = 1 h 3 N N ! ∫ ⋯ ∫ d p 1 x ⋯ d p N z d x 1 ⋯ d z N e − β H , H = ∑ i = 1 N p i 2 2 m Z=\frac{1}{h^{3N}N!}\int\!\cdots\!\int dp_{1x}\cdots dp_{Nz}\,dx_1\cdots dz_N\ e^{-\beta H},
\qquad H=\sum_{i=1}^{N}\frac{p_i^2}{2m} Z = h 3 N N ! 1 ∫ ⋯ ∫ d p 1 x ⋯ d p N z d x 1 ⋯ d z N e − β H , H = i = 1 ∑ N 2 m p i 2 で与えられます。β = 1 / ( k B T ) \beta=1/(k_{\mathrm B}T) β = 1/ ( k B T ) h = 2 π ℏ h=2\pi\hbar h = 2 π ℏ μ \mu μ L L L V = L 3 V=L^3 V = L 3 ε k = ℏ 2 k 2 / ( 2 m ) \varepsilon_k=\hbar^2k^2/(2m) ε k = ℏ 2 k 2 / ( 2 m ) d U = d ′ Q − P d V + μ d N dU=d'Q-PdV+\mu dN d U = d ′ Q − P d V + μ d N d ′ Q = T d S d'Q=TdS d ′ Q = T d S F = U − T S F=U-TS F = U − T S
以下、熱的ド・ブロイ波長
λ ≡ h 2 π m k B T \lambda\equiv\frac{h}{\sqrt{2\pi mk_{\mathrm B}T}} λ ≡ 2 π m k B T h を導入すると式が短くなるので、これを使います。
位置積分は各粒子について V V V a = β / ( 2 m ) a=\beta/(2m) a = β / ( 2 m )
∫ − ∞ ∞ e − β p 2 / ( 2 m ) d p = 2 π m β = 2 π m k B T \int_{-\infty}^{\infty}e^{-\beta p^2/(2m)}dp=\sqrt{\frac{2\pi m}{\beta}}=\sqrt{2\pi mk_{\mathrm B}T} ∫ − ∞ ∞ e − β p 2 / ( 2 m ) d p = β 2 π m = 2 π m k B T となります。運動量成分は 3 N 3N 3 N
Z ( T , V , N ) = 1 h 3 N N ! V N ( 2 π m k B T ) 3 N / 2 = 1 N ! ( V λ 3 ) N Z(T,V,N)=\frac{1}{h^{3N}N!}V^N\left(2\pi mk_{\mathrm B}T\right)^{3N/2}
=\frac{1}{N!}\left(\frac{V}{\lambda^3}\right)^{N} Z ( T , V , N ) = h 3 N N ! 1 V N ( 2 π m k B T ) 3 N /2 = N ! 1 ( λ 3 V ) N です。ヘルムホルツの自由エネルギーは ln N ! ≈ N ln N − N \ln N!\approx N\ln N-N ln N ! ≈ N ln N − N
F = − k B T ln Z = − N k B T [ ln V N λ 3 + 1 ] F=-k_{\mathrm B}T\ln Z=-Nk_{\mathrm B}T\left[\ln\frac{V}{N\lambda^3}+1\right] F = − k B T ln Z = − N k B T [ ln N λ 3 V + 1 ] となり、λ \lambda λ V V V
P = − ( ∂ F ∂ V ) T , N = N k B T V P=-\left(\frac{\partial F}{\partial V}\right)_{T,N}=\frac{Nk_{\mathrm B}T}{V} P = − ( ∂ V ∂ F ) T , N = V N k B T です。答えは Z = ( V / λ 3 ) N / N ! Z=(V/\lambda^3)^N/N! Z = ( V / λ 3 ) N / N ! P = N k B T / V P=Nk_{\mathrm B}T/V P = N k B T / V
満たされなくなるのは示量性（相加性）です。1 / N ! 1/N! 1/ N !
F 誤 = − N k B T ln V λ 3 F_{\text{誤}}=-Nk_{\mathrm B}T\ln\frac{V}{\lambda^3} F 誤 = − N k B T ln λ 3 V となります。自由エネルギーは示量変数について 1 次の同次関数、すなわち F ( T , α V , α N ) = α F ( T , V , N ) F(T,\alpha V,\alpha N)=\alpha F(T,V,N) F ( T , α V , α N ) = α F ( T , V , N )
F 誤 ( T , α V , α N ) = α F 誤 ( T , V , N ) − α N k B T ln α F_{\text{誤}}(T,\alpha V,\alpha N)=\alpha F_{\text{誤}}(T,V,N)-\alpha Nk_{\mathrm B}T\ln\alpha F 誤 ( T , α V , α N ) = α F 誤 ( T , V , N ) − α N k B T ln α となって余分な項が残ります。原因は N ln V N\ln V N ln V ln ( V / N ) \ln(V/N) ln ( V / N ) 1 / N ! 1/N! 1/ N ! ln ( V / N λ 3 ) \ln(V/N\lambda^3) ln ( V / N λ 3 )
設問1 の F F F T T T ln ( V / N λ 3 ) = ln ( V / N ) + 3 2 ln ( 2 π m k B T / h 2 ) \ln(V/N\lambda^3)=\ln(V/N)+\tfrac32\ln(2\pi mk_{\mathrm B}T/h^2) ln ( V / N λ 3 ) = ln ( V / N ) + 2 3 ln ( 2 π m k B T / h 2 )
S = − ( ∂ F ∂ T ) V , N = N k B [ ln V N λ 3 + 1 ] + N k B T ⋅ 3 2 T S=-\left(\frac{\partial F}{\partial T}\right)_{V,N}
=Nk_{\mathrm B}\left[\ln\frac{V}{N\lambda^3}+1\right]+Nk_{\mathrm B}T\cdot\frac{3}{2T} S = − ( ∂ T ∂ F ) V , N = N k B [ ln N λ 3 V + 1 ] + N k B T ⋅ 2 T 3 となり、
S ( T , V , N ) = N k B [ ln V N λ 3 + 5 2 ] = N k B [ ln V N + 3 2 ln 2 π m k B T h 2 + 5 2 ] S(T,V,N)=Nk_{\mathrm B}\left[\ln\frac{V}{N\lambda^3}+\frac52\right]
=Nk_{\mathrm B}\left[\ln\frac{V}{N}+\frac32\ln\frac{2\pi mk_{\mathrm B}T}{h^2}+\frac52\right] S ( T , V , N ) = N k B [ ln N λ 3 V + 2 5 ] = N k B [ ln N V + 2 3 ln h 2 2 π m k B T + 2 5 ] です。これはザッカー・テトロードの式です。
T → 0 T\to0 T → 0 λ ∝ T − 1 / 2 → ∞ \lambda\propto T^{-1/2}\to\infty λ ∝ T − 1/2 → ∞ V / ( N λ 3 ) → 0 V/(N\lambda^3)\to0 V / ( N λ 3 ) → 0 ln \ln ln S → − ∞ S\to-\infty S → − ∞ S → 0 S\to0 S → 0 λ \lambda λ ( V / N ) 1 / 3 (V/N)^{1/3} ( V / N ) 1/3
内部エネルギーは
U = F + T S = − N k B T [ ln V N λ 3 + 1 ] + N k B T [ ln V N λ 3 + 5 2 ] = 3 2 N k B T U=F+TS=-Nk_{\mathrm B}T\left[\ln\frac{V}{N\lambda^3}+1\right]+Nk_{\mathrm B}T\left[\ln\frac{V}{N\lambda^3}+\frac52\right]=\frac32Nk_{\mathrm B}T U = F + T S = − N k B T [ ln N λ 3 V + 1 ] + N k B T [ ln N λ 3 V + 2 5 ] = 2 3 N k B T なので
C V ( T , V , N ) = ( ∂ U ∂ T ) V , N = 3 2 N k B C_V(T,V,N)=\left(\frac{\partial U}{\partial T}\right)_{V,N}=\frac32Nk_{\mathrm B} C V ( T , V , N ) = ( ∂ T ∂ U ) V , N = 2 3 N k B です。C V = T ( ∂ S / ∂ T ) V = T ⋅ 3 2 N k B / T C_V=T(\partial S/\partial T)_V=T\cdot\tfrac32Nk_{\mathrm B}/T C V = T ( ∂ S / ∂ T ) V = T ⋅ 2 3 N k B / T 1 2 k B \tfrac12k_{\mathrm B} 2 1 k B
一辺 L L L e i k ⋅ r e^{i\boldsymbol{k}\cdot\boldsymbol{r}} e i k ⋅ r e i k x L = 1 e^{ik_xL}=1 e i k x L = 1
k x = 2 π n x L , k y = 2 π n y L , k z = 2 π n z L , n x , n y , n z ∈ Z k_x=\frac{2\pi n_x}{L},\quad k_y=\frac{2\pi n_y}{L},\quad k_z=\frac{2\pi n_z}{L},
\qquad n_x,n_y,n_z\in\mathbb{Z} k x = L 2 π n x , k y = L 2 π n y , k z = L 2 π n z , n x , n y , n z ∈ Z です。n n n 0 0 0 k k k ( 2 π / L ) 3 = ( 2 π ) 3 / V (2\pi/L)^3=(2\pi)^3/V ( 2 π / L ) 3 = ( 2 π ) 3 / V
グランドカノニカル分布では Ξ = ∑ 状態 e − β ( E − μ N ) \Xi=\sum_{\text{状態}}e^{-\beta(E-\mu N)} Ξ = ∑ 状態 e − β ( E − μ N ) μ \mu μ
N ‾ = 1 Ξ ∑ 状態 N e − β ( E − μ N ) = 1 β ( ∂ ln Ξ ∂ μ ) T , V \overline{N}=\frac{1}{\Xi}\sum_{\text{状態}}N\,e^{-\beta(E-\mu N)}
=\frac{1}{\beta}\left(\frac{\partial\ln\Xi}{\partial\mu}\right)_{T,V} N = Ξ 1 状態 ∑ N e − β ( E − μ N ) = β 1 ( ∂ μ ∂ ln Ξ ) T , V 与えられた大分配関数の対数は ln Ξ = ∑ k ln ( 1 + e − β ( ε k − μ ) ) \ln\Xi=\sum_k\ln\left(1+e^{-\beta(\varepsilon_k-\mu)}\right) ln Ξ = ∑ k ln ( 1 + e − β ( ε k − μ ) )
N ‾ = 1 β ∑ k β e − β ( ε k − μ ) 1 + e − β ( ε k − μ ) = ∑ k 1 e β ( ε k − μ ) + 1 = ∑ k f ( ε k ) \overline{N}=\frac{1}{\beta}\sum_k\frac{\beta e^{-\beta(\varepsilon_k-\mu)}}{1+e^{-\beta(\varepsilon_k-\mu)}}
=\sum_k\frac{1}{e^{\beta(\varepsilon_k-\mu)}+1}=\sum_kf(\varepsilon_k) N = β 1 k ∑ 1 + e − β ( ε k − μ ) β e − β ( ε k − μ ) = k ∑ e β ( ε k − μ ) + 1 1 = k ∑ f ( ε k ) となります。示されました。各項は波数 k k k ⟨ n k ⟩ \langle n_k\rangle ⟨ n k ⟩ 0 ≤ f ≤ 1 0\le f\le1 0 ≤ f ≤ 1
L L L k k k ( 2 π ) 3 / V (2\pi)^3/V ( 2 π ) 3 / V
∑ k ⟶ V ( 2 π ) 3 ∫ d 3 k \sum_k\ \longrightarrow\ \frac{V}{(2\pi)^3}\int d^3k k ∑ ⟶ ( 2 π ) 3 V ∫ d 3 k です。高温近似 f ≈ e − β ( ε k − μ ) f\approx e^{-\beta(\varepsilon_k-\mu)} f ≈ e − β ( ε k − μ )
N = e β μ V ( 2 π ) 3 ∫ d 3 k e − β ℏ 2 k 2 / ( 2 m ) = e β μ V ( 2 π ) 3 ( 2 π m k B T ℏ 2 ) 3 / 2 N=e^{\beta\mu}\frac{V}{(2\pi)^3}\int d^3k\ e^{-\beta\hbar^2k^2/(2m)}
=e^{\beta\mu}\frac{V}{(2\pi)^3}\left(\frac{2\pi m k_{\mathrm B}T}{\hbar^2}\right)^{3/2} N = e β μ ( 2 π ) 3 V ∫ d 3 k e − β ℏ 2 k 2 / ( 2 m ) = e β μ ( 2 π ) 3 V ( ℏ 2 2 π m k B T ) 3/2 となります。ガウス積分は各成分について ∫ d k x e − β ℏ 2 k x 2 / ( 2 m ) = 2 π m k B T / ℏ 2 \int dk_x\,e^{-\beta\hbar^2k_x^2/(2m)}=\sqrt{2\pi mk_{\mathrm B}T/\hbar^2} ∫ d k x e − β ℏ 2 k x 2 / ( 2 m ) = 2 π m k B T / ℏ 2 ( 2 π ) 3 (2\pi)^3 ( 2 π ) 3
N = V ( m k B T 2 π ℏ 2 ) 3 / 2 e β μ = V λ 3 e β μ N=V\left(\frac{mk_{\mathrm B}T}{2\pi\hbar^2}\right)^{3/2}e^{\beta\mu}=\frac{V}{\lambda^3}e^{\beta\mu} N = V ( 2 π ℏ 2 m k B T ) 3/2 e β μ = λ 3 V e β μ です。したがって
μ = k B T ln N λ 3 V = k B T ln [ N V ( 2 π ℏ 2 m k B T ) 3 / 2 ] \mu=k_{\mathrm B}T\ln\frac{N\lambda^3}{V}
=k_{\mathrm B}T\ln\left[\frac{N}{V}\left(\frac{2\pi\hbar^2}{mk_{\mathrm B}T}\right)^{3/2}\right] μ = k B T ln V N λ 3 = k B T ln [ V N ( m k B T 2 π ℏ 2 ) 3/2 ] が答えです。高温では N λ 3 / V ≪ 1 N\lambda^3/V\ll1 N λ 3 / V ≪ 1 μ < 0 \mu<0 μ < 0 ∣ μ ∣ ≫ k B T |\mu|\gg k_{\mathrm B}T ∣ μ ∣ ≫ k B T
同じ温度領域では e − β ( ε k − μ ) ≤ e β μ ≪ 1 e^{-\beta(\varepsilon_k-\mu)}\le e^{\beta\mu}\ll1 e − β ( ε k − μ ) ≤ e β μ ≪ 1 ln ( 1 + x ) ≈ x \ln(1+x)\approx x ln ( 1 + x ) ≈ x
ln Ξ = ∑ k ln ( 1 + e − β ( ε k − μ ) ) ≈ ∑ k e − β ( ε k − μ ) = N \ln\Xi=\sum_k\ln\left(1+e^{-\beta(\varepsilon_k-\mu)}\right)\approx\sum_ke^{-\beta(\varepsilon_k-\mu)}=N ln Ξ = k ∑ ln ( 1 + e − β ( ε k − μ ) ) ≈ k ∑ e − β ( ε k − μ ) = N です。グランドポテンシャルは J = − k B T ln Ξ = − N k B T J=-k_{\mathrm B}T\ln\Xi=-Nk_{\mathrm B}T J = − k B T ln Ξ = − N k B T P V = N k B T PV=Nk_{\mathrm B}T P V = N k B T
内部エネルギーは同じ近似で
U = ∑ k ε k f ( ε k ) ≈ e β μ V ( 2 π ) 3 ∫ d 3 k ℏ 2 k 2 2 m e − β ℏ 2 k 2 / ( 2 m ) = 3 2 N k B T U=\sum_k\varepsilon_kf(\varepsilon_k)\approx e^{\beta\mu}\frac{V}{(2\pi)^3}\int d^3k\ \frac{\hbar^2k^2}{2m}e^{-\beta\hbar^2k^2/(2m)}=\frac32Nk_{\mathrm B}T U = k ∑ ε k f ( ε k ) ≈ e β μ ( 2 π ) 3 V ∫ d 3 k 2 m ℏ 2 k 2 e − β ℏ 2 k 2 / ( 2 m ) = 2 3 N k B T となります（ガウス積分の比 ∫ k 2 e − a k 2 d 3 k / ∫ e − a k 2 d 3 k = 3 / ( 2 a ) \int k^2e^{-ak^2}d^3k\big/\int e^{-ak^2}d^3k=3/(2a) ∫ k 2 e − a k 2 d 3 k / ∫ e − a k 2 d 3 k = 3/ ( 2 a ) J = U − T S − μ N J=U-TS-\mu N J = U − T S − μ N S S S
S = U − μ N − J T = 3 2 N k B T − μ N + N k B T T = − μ T N + 5 2 N k B S=\frac{U-\mu N-J}{T}=\frac{\tfrac32Nk_{\mathrm B}T-\mu N+Nk_{\mathrm B}T}{T}
=-\frac{\mu}{T}N+\frac52Nk_{\mathrm B} S = T U − μ N − J = T 2 3 N k B T − μ N + N k B T = − T μ N + 2 5 N k B です。ここで − μ / T = ∣ μ ∣ / T ≫ k B -\mu/T=|\mu|/T\gg k_{\mathrm B} − μ / T = ∣ μ ∣/ T ≫ k B ∣ μ ∣ ≫ k B T |\mu|\gg k_{\mathrm B}T ∣ μ ∣ ≫ k B T
S ≈ − μ T N S\approx-\frac{\mu}{T}N S ≈ − T μ N が示されました。設問3・設問7 と突き合わせると − μ N / T = N k B ln ( V / N λ 3 ) -\mu N/T=Nk_{\mathrm B}\ln\left(V/N\lambda^3\right) − μ N / T = N k B ln ( V / N λ 3 ) 5 2 N k B \tfrac52Nk_{\mathrm B} 2 5 N k B ∣ μ ∣ / ( k B T ) ≫ 1 |\mu|/(k_{\mathrm B}T)\gg1 ∣ μ ∣/ ( k B T ) ≫ 1
低温側は C V = γ T C_V=\gamma T C V = γ T
S ( T ) = ∫ 0 T C V ( T ′ ) T ′ d T ′ = ∫ 0 T γ d T ′ = γ T S(T)=\int_0^T\frac{C_V(T')}{T'}dT'=\int_0^T\gamma\,dT'=\gamma T S ( T ) = ∫ 0 T T ′ C V ( T ′ ) d T ′ = ∫ 0 T γ d T ′ = γ T です。積分定数は第三法則 S ( 0 ) = 0 S(0)=0 S ( 0 ) = 0
S ≈ − μ T N = N k B [ ln V N + 3 2 ln 2 π m k B T h 2 ] S\approx-\frac{\mu}{T}N=Nk_{\mathrm B}\left[\ln\frac{V}{N}+\frac32\ln\frac{2\pi mk_{\mathrm B}T}{h^2}\right] S ≈ − T μ N = N k B [ ln N V + 2 3 ln h 2 2 π m k B T ] で、3 2 N k B ln T \tfrac32Nk_{\mathrm B}\ln T 2 3 N k B ln T T F T_{\mathrm F} T F λ ( T F ) ∼ ( V / N ) 1 / 3 \lambda(T_{\mathrm F})\sim(V/N)^{1/3} λ ( T F ) ∼ ( V / N ) 1/3 γ T F ∼ N k B \gamma T_{\mathrm F}\sim Nk_{\mathrm B} γ T F ∼ N k B
グラフの特徴を書き出すと次のようになります。横軸を T T T S S S ( 0 , 0 ) (0,0) ( 0 , 0 ) T ≪ T F T\ll T_{\mathrm F} T ≪ T F γ \gamma γ S = γ T S=\gamma T S = γ T T ∼ T F T\sim T_{\mathrm F} T ∼ T F S ∼ N k B S\sim Nk_{\mathrm B} S ∼ N k B T ≫ T F T\gg T_{\mathrm F} T ≫ T F d 2 S / d T 2 < 0 d^2S/dT^2<0 d 2 S / d T 2 < 0 T → ∞ T\to\infty T → ∞ 0 0 0
設問3 の古典的なエントロピー S 古典 = N k B [ ln ( V / N λ 3 ) + 5 2 ] S_{\text{古典}}=Nk_{\mathrm B}\left[\ln(V/N\lambda^3)+\tfrac52\right] S 古典 = N k B [ ln ( V / N λ 3 ) + 2 5 ] 5 2 N k B \tfrac52Nk_{\mathrm B} 2 5 N k B S 古典 S_{\text{古典}} S 古典 V / ( N λ 3 ) = e − 5 / 2 V/(N\lambda^3)=e^{-5/2} V / ( N λ 3 ) = e − 5/2 T F T_{\mathrm F} T F S = 0 S=0 S = 0 T → 0 T\to0 T → 0 − ∞ -\infty − ∞ 0 0 0
S T 0 TF フェルミ粒子系 古典（設問3） 半径 a a a b b b ϵ \epsilon ϵ μ \mu μ z z z ( r , θ , z ) (r,\theta,z) ( r , θ , z )
考えるのは TEM モード、すなわち
E z = B z = 0 , ∂ ∂ θ ( すべての場 ) = 0 E_z=B_z=0,\qquad \frac{\partial}{\partial\theta}(\text{すべての場})=0 E z = B z = 0 , ∂ θ ∂ ( すべての場 ) = 0 を満たす伝搬モードです。絶縁体中では
∇ × E = − ∂ B ∂ t , ∇ × B = μ ϵ ∂ E ∂ t , ∇ ⋅ B = 0 , ∇ ⋅ E = 0 \nabla\times\boldsymbol{E}=-\frac{\partial\boldsymbol{B}}{\partial t},\quad
\nabla\times\boldsymbol{B}=\mu\epsilon\frac{\partial\boldsymbol{E}}{\partial t},\quad
\nabla\cdot\boldsymbol{B}=0,\quad
\nabla\cdot\boldsymbol{E}=0 ∇ × E = − ∂ t ∂ B , ∇ × B = μ ϵ ∂ t ∂ E , ∇ ⋅ B = 0 , ∇ ⋅ E = 0 が成り立ち、円筒座標での発散と回転は問題文に与えられた表式を使います。最後の 2 問ではパルスの伝搬を考え、ϵ , μ \epsilon,\mu ϵ , μ
まず E θ E_\theta E θ z z z E z E_z E z θ \theta θ B z = 0 B_z=0 B z = 0
1 r [ ∂ ∂ r ( r E θ ) − ∂ E r ∂ θ ] = − ∂ B z ∂ t = 0 ⟹ ∂ ∂ r ( r E θ ) = 0 \frac{1}{r}\left[\frac{\partial}{\partial r}(rE_\theta)-\frac{\partial E_r}{\partial\theta}\right]=-\frac{\partial B_z}{\partial t}=0
\ \Longrightarrow\ \frac{\partial}{\partial r}(rE_\theta)=0 r 1 [ ∂ r ∂ ( r E θ ) − ∂ θ ∂ E r ] = − ∂ t ∂ B z = 0 ⟹ ∂ r ∂ ( r E θ ) = 0 となり、r E θ rE_\theta r E θ r r r r E θ = h ( z , t ) rE_\theta=h(z,t) r E θ = h ( z , t ) r = a r=a r = a E θ ( a , z , t ) = 0 E_\theta(a,z,t)=0 E θ ( a , z , t ) = 0 h ( z , t ) = 0 h(z,t)=0 h ( z , t ) = 0 a ≤ r ≤ b a\le r\le b a ≤ r ≤ b
E θ = 0 E_\theta=0 E θ = 0 が言えました。
次に B r B_r B r ∇ ⋅ B = 0 \nabla\cdot\boldsymbol{B}=0 ∇ ⋅ B = 0 B z = 0 B_z=0 B z = 0 ∂ θ = 0 \partial_\theta=0 ∂ θ = 0
1 r ∂ ∂ r ( r B r ) = 0 ⟹ r B r = g ( z , t ) \frac{1}{r}\frac{\partial}{\partial r}(rB_r)=0
\ \Longrightarrow\ rB_r=g(z,t) r 1 ∂ r ∂ ( r B r ) = 0 ⟹ r B r = g ( z , t ) で、r B r rB_r r B r r r r r r r θ \theta θ
− ∂ E θ ∂ z = − ∂ B r ∂ t , ∂ B r ∂ z = μ ϵ ∂ E θ ∂ t -\frac{\partial E_\theta}{\partial z}=-\frac{\partial B_r}{\partial t},\qquad
\frac{\partial B_r}{\partial z}=\mu\epsilon\frac{\partial E_\theta}{\partial t} − ∂ z ∂ E θ = − ∂ t ∂ B r , ∂ z ∂ B r = μ ϵ ∂ t ∂ E θ を使うと、E θ = 0 E_\theta=0 E θ = 0 ∂ g / ∂ t = 0 \partial g/\partial t=0 ∂ g / ∂ t = 0 ∂ g / ∂ z = 0 \partial g/\partial z=0 ∂ g / ∂ z = 0 g g g B r = g / r B_r=g/r B r = g / r B r ( a ) = 0 B_r(a)=0 B r ( a ) = 0 g = 0 g=0 g = 0
B r = 0 B_r=0 B r = 0 です。g ≠ 0 g\neq0 g = 0 B r = g / r B_r=g/r B r = g / r r r r ℓ \ell ℓ 2 π ℓ g 2\pi\ell g 2 π ℓ g ∇ ⋅ B = 0 \nabla\cdot\boldsymbol{B}=0 ∇ ⋅ B = 0 B r = E θ = 0 B_r=E_\theta=0 B r = E θ = 0 E r E_r E r B θ B_\theta B θ
Faraday の法則の θ \theta θ r r r E z = B z = 0 E_z=B_z=0 E z = B z = 0 ∂ θ = 0 \partial_\theta=0 ∂ θ = 0
∂ E r ∂ z = − ∂ B θ ∂ t , − ∂ B θ ∂ z = μ ϵ ∂ E r ∂ t \frac{\partial E_r}{\partial z}=-\frac{\partial B_\theta}{\partial t},\qquad
-\frac{\partial B_\theta}{\partial z}=\mu\epsilon\frac{\partial E_r}{\partial t} ∂ z ∂ E r = − ∂ t ∂ B θ , − ∂ z ∂ B θ = μ ϵ ∂ t ∂ E r です。前者を z z z B θ B_\theta B θ
∂ 2 E r ∂ z 2 = − ∂ ∂ t ∂ B θ ∂ z = − ∂ ∂ t ( − μ ϵ ∂ E r ∂ t ) = μ ϵ ∂ 2 E r ∂ t 2 \frac{\partial^2E_r}{\partial z^2}=-\frac{\partial}{\partial t}\frac{\partial B_\theta}{\partial z}
=-\frac{\partial}{\partial t}\left(-\mu\epsilon\frac{\partial E_r}{\partial t}\right)
=\mu\epsilon\frac{\partial^2E_r}{\partial t^2} ∂ z 2 ∂ 2 E r = − ∂ t ∂ ∂ z ∂ B θ = − ∂ t ∂ ( − μ ϵ ∂ t ∂ E r ) = μ ϵ ∂ t 2 ∂ 2 E r したがって求める波動方程式は
∂ 2 E r ∂ z 2 − μ ϵ ∂ 2 E r ∂ t 2 = 0 \frac{\partial^2E_r}{\partial z^2}-\mu\epsilon\frac{\partial^2E_r}{\partial t^2}=0 ∂ z 2 ∂ 2 E r − μ ϵ ∂ t 2 ∂ 2 E r = 0 です。r r r ∇ ⋅ E = 0 \nabla\cdot\boldsymbol{E}=0 ∇ ⋅ E = 0 ∂ ( r E r ) / ∂ r = 0 \partial(rE_r)/\partial r=0 ∂ ( r E r ) / ∂ r = 0 E r = A ( z , t ) / r E_r=A(z,t)/r E r = A ( z , t ) / r r r r r r r
∂ 2 E r ∂ r 2 + 1 r ∂ E r ∂ r − E r r 2 = A ( 2 r 3 − 1 r 3 − 1 r 3 ) = 0 \frac{\partial^2E_r}{\partial r^2}+\frac1r\frac{\partial E_r}{\partial r}-\frac{E_r}{r^2}
=A\left(\frac{2}{r^3}-\frac{1}{r^3}-\frac{1}{r^3}\right)=0 ∂ r 2 ∂ 2 E r + r 1 ∂ r ∂ E r − r 2 E r = A ( r 3 2 − r 3 1 − r 3 1 ) = 0 と恒等的に消え、上の 1 次元波動方程式と矛盾しません。
∇ ⋅ E = 0 \nabla\cdot\boldsymbol{E}=0 ∇ ⋅ E = 0
1 r d d r ( r E ( r ) ) = 0 ⟹ r E ( r ) = const \frac{1}{r}\frac{d}{dr}\left(r\mathcal{E}(r)\right)=0
\ \Longrightarrow\ r\mathcal{E}(r)=\text{const} r 1 d r d ( r E ( r ) ) = 0 ⟹ r E ( r ) = const です。E ( a ) = E 0 \mathcal{E}(a)=E_0 E ( a ) = E 0
E ( r ) = E 0 a r ( a ≤ r ≤ b ) \mathcal{E}(r)=E_0\frac{a}{r}\qquad(a\le r\le b) E ( r ) = E 0 r a ( a ≤ r ≤ b ) が答えです。芯線を単位長さあたり Λ \Lambda Λ Λ / ( 2 π ϵ r ) \Lambda/(2\pi\epsilon r) Λ/ ( 2 π ϵr )
E r = E ( r ) e i k z − i ω t E_r=\mathcal{E}(r)e^{ikz-i\omega t} E r = E ( r ) e ik z − iω t − k 2 + μ ϵ ω 2 = 0 -k^2+\mu\epsilon\omega^2=0 − k 2 + μ ϵ ω 2 = 0 ω = k / μ ϵ \omega=k/\sqrt{\mu\epsilon} ω = k / μ ϵ
v p = ω k = 1 μ ϵ v_{\mathrm p}=\frac{\omega}{k}=\frac{1}{\sqrt{\mu\epsilon}} v p = k ω = μ ϵ 1 です。これは絶縁体中の光速そのもので、k k k ϵ , μ \epsilon,\mu ϵ , μ 1 / μ 0 ϵ 0 = c 1/\sqrt{\mu_0\epsilon_0}=c 1/ μ 0 ϵ 0 = c
まず B θ B_\theta B θ E r E_r E r i k E = i ω B ik\mathcal{E}=i\omega\mathcal{B} ik E = iω B
B ( r ) = k ω E ( r ) = μ ϵ E 0 a r \mathcal{B}(r)=\frac{k}{\omega}\mathcal{E}(r)=\sqrt{\mu\epsilon}\,E_0\frac{a}{r} B ( r ) = ω k E ( r ) = μ ϵ E 0 r a です。実場は E 0 E_0 E 0 E r = ( E 0 a / r ) cos ( k z − ω t ) E_r=(E_0a/r)\cos(kz-\omega t) E r = ( E 0 a / r ) cos ( k z − ω t ) B θ = μ ϵ ( E 0 a / r ) cos ( k z − ω t ) B_\theta=\sqrt{\mu\epsilon}\,(E_0a/r)\cos(kz-\omega t) B θ = μ ϵ ( E 0 a / r ) cos ( k z − ω t ) e r × e θ = e z \boldsymbol{e}_r\times\boldsymbol{e}_\theta=\boldsymbol{e}_z e r × e θ = e z
S z = E r B θ μ = ϵ μ ( E 0 a r ) 2 cos 2 ( k z − ω t ) S_z=\frac{E_rB_\theta}{\mu}=\sqrt{\frac{\epsilon}{\mu}}\left(\frac{E_0a}{r}\right)^2\cos^2(kz-\omega t) S z = μ E r B θ = μ ϵ ( r E 0 a ) 2 cos 2 ( k z − ω t ) で、断面（a ≤ r ≤ b a\le r\le b a ≤ r ≤ b
P ( z , t ) = ∫ a b S z 2 π r d r = 2 π ϵ μ E 0 2 a 2 cos 2 ( k z − ω t ) ∫ a b d r r = 2 π ϵ μ E 0 2 a 2 ln b a cos 2 ( k z − ω t ) P(z,t)=\int_a^bS_z\,2\pi r\,dr
=2\pi\sqrt{\frac{\epsilon}{\mu}}\,E_0^2a^2\cos^2(kz-\omega t)\int_a^b\frac{dr}{r}
=2\pi\sqrt{\frac{\epsilon}{\mu}}\,E_0^2a^2\ln\frac{b}{a}\ \cos^2(kz-\omega t) P ( z , t ) = ∫ a b S z 2 π r d r = 2 π μ ϵ E 0 2 a 2 cos 2 ( k z − ω t ) ∫ a b r d r = 2 π μ ϵ E 0 2 a 2 ln a b cos 2 ( k z − ω t ) となります。時間平均をとれば
P ‾ = π ϵ μ E 0 2 a 2 ln b a \overline{P}=\pi\sqrt{\frac{\epsilon}{\mu}}\,E_0^2a^2\ln\frac{b}{a} P = π μ ϵ E 0 2 a 2 ln a b です。μ / ϵ \sqrt{\mu/\epsilon} μ / ϵ Ω \Omega Ω a 2 E 0 2 ϵ / μ a^2E_0^2\sqrt{\epsilon/\mu} a 2 E 0 2 ϵ / μ m 2 × ( V / m ) 2 / Ω = W \mathrm{m^2}\times(\mathrm{V/m})^2/\Omega=\mathrm{W} m 2 × ( V/m ) 2 /Ω = W
TEM モードでは E θ = 0 E_\theta=0 E θ = 0 ∂ θ = 0 \partial_\theta=0 ∂ θ = 0 ( ∇ × E ) z = 0 (\nabla\times\boldsymbol{E})_z=0 ( ∇ × E ) z = 0
V ( z , t ) = ∫ a b E r d r = E 0 a e i k z − i ω t ∫ a b d r r = E 0 a ln b a e i k z − i ω t V(z,t)=\int_a^bE_r\,dr=E_0a\,e^{ikz-i\omega t}\int_a^b\frac{dr}{r}=E_0a\ln\frac{b}{a}\ e^{ikz-i\omega t} V ( z , t ) = ∫ a b E r d r = E 0 a e ik z − iω t ∫ a b r d r = E 0 a ln a b e ik z − iω t 実場で書けば
V ( z , t ) = E 0 a ln b a cos ( k z − ω t ) V(z,t)=E_0a\ln\frac{b}{a}\ \cos(kz-\omega t) V ( z , t ) = E 0 a ln a b cos ( k z − ω t ) です。E r > 0 E_r>0 E r > 0
芯線を流れる電流は Ampère–Maxwell の法則の積分形から求めます。半径 r r r a < r < b a<r<b a < r < b E z = 0 E_z=0 E z = 0
2 π r B θ = μ I ( z , t ) ⟹ I ( z , t ) = 2 π r μ μ ϵ E 0 a r e i k z − i ω t = 2 π a E 0 ϵ μ e i k z − i ω t 2\pi rB_\theta=\mu I(z,t)
\ \Longrightarrow\ I(z,t)=\frac{2\pi r}{\mu}\sqrt{\mu\epsilon}\,\frac{E_0a}{r}e^{ikz-i\omega t}
=2\pi aE_0\sqrt{\frac{\epsilon}{\mu}}\ e^{ikz-i\omega t} 2 π r B θ = μ I ( z , t ) ⟹ I ( z , t ) = μ 2 π r μ ϵ r E 0 a e ik z − iω t = 2 π a E 0 μ ϵ e ik z − iω t です。r r r
Z 0 = V ( z , t ) I ( z , t ) = E 0 a ln ( b / a ) 2 π a E 0 ϵ / μ = 1 2 π μ ϵ ln b a Z_0=\frac{V(z,t)}{I(z,t)}=\frac{E_0a\ln(b/a)}{2\pi aE_0\sqrt{\epsilon/\mu}}
=\frac{1}{2\pi}\sqrt{\frac{\mu}{\epsilon}}\ln\frac{b}{a} Z 0 = I ( z , t ) V ( z , t ) = 2 π a E 0 ϵ / μ E 0 a ln ( b / a ) = 2 π 1 ϵ μ ln a b です。z , t z,t z , t V I VI V I
V I = 2 π a 2 E 0 2 ϵ μ ln b a cos 2 ( k z − ω t ) VI=2\pi a^2E_0^2\sqrt{\frac{\epsilon}{\mu}}\ln\frac{b}{a}\ \cos^2(kz-\omega t) V I = 2 π a 2 E 0 2 μ ϵ ln a b cos 2 ( k z − ω t ) で、設問5 で断面積分したポインティング流と完全に一致します。数値でも、μ 0 / ϵ 0 ≈ 377 Ω \sqrt{\mu_0/\epsilon_0}\approx377\,\Omega μ 0 / ϵ 0 ≈ 377 Ω Z 0 ≈ 60 ln ( b / a ) Ω / ϵ / ϵ 0 Z_0\approx60\ln(b/a)\,\Omega/\sqrt{\epsilon/\epsilon_0} Z 0 ≈ 60 ln ( b / a ) Ω/ ϵ / ϵ 0 2.25 2.25 2.25 b / a = 3.5 b/a=3.5 b / a = 3.5 50 Ω 50\,\Omega 50 Ω
ω 0 ≡ ω ( k 0 ) \omega_0\equiv\omega(k_0) ω 0 ≡ ω ( k 0 )
i k z − i ω ( k ) t = i k z − i ω 0 t − i A ( k − k 0 ) t = i k ( z − A t ) − i ( ω 0 − A k 0 ) t ikz-i\omega(k)t=ikz-i\omega_0t-iA(k-k_0)t=ik(z-At)-i(\omega_0-Ak_0)t ik z − iω ( k ) t = ik z − i ω 0 t − i A ( k − k 0 ) t = ik ( z − A t ) − i ( ω 0 − A k 0 ) t 第 2 項は k k k
V ( z , t ) = e − i ( ω 0 − A k 0 ) t ⋅ 1 2 π ∫ − ∞ ∞ v ( k ) e i k ( z − A t ) d k V(z,t)=e^{-i(\omega_0-Ak_0)t}\cdot\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty}v(k)e^{ik(z-At)}dk V ( z , t ) = e − i ( ω 0 − A k 0 ) t ⋅ 2 π 1 ∫ − ∞ ∞ v ( k ) e ik ( z − A t ) d k ここで F ( ζ ) ≡ 1 2 π ∫ − ∞ ∞ v ( k ) e i k ζ d k F(\zeta)\equiv\dfrac{1}{\sqrt{2\pi}}\displaystyle\int_{-\infty}^{\infty}v(k)e^{ik\zeta}dk F ( ζ ) ≡ 2 π 1 ∫ − ∞ ∞ v ( k ) e ik ζ d k F ( ζ ) = V ( ζ , 0 ) F(\zeta)=V(\zeta,0) F ( ζ ) = V ( ζ , 0 )
V ( z , t ) = e − i ( ω 0 − A k 0 ) t F ( z − A t ) V(z,t)=e^{-i(\omega_0-Ak_0)t}F(z-At) V ( z , t ) = e − i ( ω 0 − A k 0 ) t F ( z − A t ) と書けます。前の因子は絶対値 1 の位相因子で、しかも z z z
∣ V ( z , t ) ∣ = ∣ F ( z − A t ) ∣ = ∣ V ( z − A t , 0 ) ∣ |V(z,t)|=|F(z-At)|=|V(z-At,\,0)| ∣ V ( z , t ) ∣ = ∣ F ( z − A t ) ∣ = ∣ V ( z − A t , 0 ) ∣ となり、振幅の波形は t = 0 t=0 t = 0 A A A A = d ω / d k ∣ k 0 A=d\omega/dk|_{k_0} A = d ω / d k ∣ k 0 ω 0 / k 0 \omega_0/k_0 ω 0 / k 0 ∣ V ∣ |V| ∣ V ∣ ω ( k ) \omega(k) ω ( k ) v ( k ) v(k) v ( k ) k 0 k_0 k 0
熱容量 C C C R k R_k R k T b T_{\mathrm b} T b R R R V V V P = V 2 / R P=V^2/R P = V 2 / R T s T_{\mathrm s} T s
前半（設問1〜3）は矩形パルス状の電圧を加える緩和法、後半（設問4〜7）は交流電圧を長時間加える交流法です。
試料に流れ込む熱は、ヒーターからの P ( t ) P(t) P ( t ) ( T s − T b ) / R k (T_{\mathrm s}-T_{\mathrm b})/R_k ( T s − T b ) / R k C d T s / d t C\,dT_{\mathrm s}/dt C d T s / d t
T = T b , A = C \boxed{\ \mathcal{T}\ }=T_{\mathrm b},\qquad \boxed{\ \mathcal{A}\ }=C T = T b , A = C です。つまり
P ( t ) − 1 R k { T s ( t ) − T b } = C d T s ( t ) d t P(t)-\frac{1}{R_k}\left\{T_{\mathrm s}(t)-T_{\mathrm b}\right\}=C\frac{dT_{\mathrm s}(t)}{dt} P ( t ) − R k 1 { T s ( t ) − T b } = C d t d T s ( t ) が満たすべき方程式です。以下 Δ T s ≡ T s − T b \Delta T_{\mathrm s}\equiv T_{\mathrm s}-T_{\mathrm b} Δ T s ≡ T s − T b τ ≡ R k C \tau\equiv R_kC τ ≡ R k C
0 < t < t 1 0<t<t_1 0 < t < t 1 V 0 V_0 V 0 P = P 0 ≡ V 0 2 / R P=P_0\equiv V_0^2/R P = P 0 ≡ V 0 2 / R
C d Δ T s d t = P 0 − Δ T s R k C\frac{d\Delta T_{\mathrm s}}{dt}=P_0-\frac{\Delta T_{\mathrm s}}{R_k} C d t d Δ T s = P 0 − R k Δ T s となり、初期条件 Δ T s ( 0 ) = 0 \Delta T_{\mathrm s}(0)=0 Δ T s ( 0 ) = 0 t < 0 t<0 t < 0
T s ( t ) = T b + V 0 2 R k R ( 1 − e − t / τ ) ( 0 < t < t 1 ) T_{\mathrm s}(t)=T_{\mathrm b}+\frac{V_0^2R_k}{R}\left(1-e^{-t/\tau}\right)\qquad(0<t<t_1) T s ( t ) = T b + R V 0 2 R k ( 1 − e − t / τ ) ( 0 < t < t 1 ) です。t ≫ τ t\gg\tau t ≫ τ T s → T b + P 0 R k T_{\mathrm s}\to T_{\mathrm b}+P_0R_k T s → T b + P 0 R k
t > t 1 t>t_1 t > t 1 P = 0 P=0 P = 0 C d Δ T s / d t = − Δ T s / R k C\,d\Delta T_{\mathrm s}/dt=-\Delta T_{\mathrm s}/R_k C d Δ T s / d t = − Δ T s / R k t = t 1 t=t_1 t = t 1
T s ( t ) = T b + V 0 2 R k R ( 1 − e − t 1 / τ ) e − ( t − t 1 ) / τ ( t > t 1 ) T_{\mathrm s}(t)=T_{\mathrm b}+\frac{V_0^2R_k}{R}\left(1-e^{-t_1/\tau}\right)e^{-(t-t_1)/\tau}\qquad(t>t_1) T s ( t ) = T b + R V 0 2 R k ( 1 − e − t 1 / τ ) e − ( t − t 1 ) / τ ( t > t 1 ) です。立ち上がりも立ち下がりも同じ時定数 τ = R k C \tau=R_kC τ = R k C
図3 から次の 3 つを読み取ります。熱浴温度は T b = 100 K T_{\mathrm b}=100\,\mathrm{K} T b = 100 K t 1 = 5 s t_1=5\,\mathrm{s} t 1 = 5 s Δ T s = 2.0 K \Delta T_{\mathrm s}=2.0\,\mathrm{K} Δ T s = 2.0 K t = 4 t=4 t = 4 5 5 5 t > t 1 t>t_1 t > t 1 102 K 102\,\mathrm{K} 102 K 101 K 101\,\mathrm{K} 101 K 1.0 s 1.0\,\mathrm{s} 1.0 s 101 K 101\,\mathrm{K} 101 K 100.5 K 100.5\,\mathrm{K} 100.5 K 1.0 s 1.0\,\mathrm{s} 1.0 s t 1 / 2 = 1.0 s t_{1/2}=1.0\,\mathrm{s} t 1/2 = 1.0 s
熱流は
P 0 = V 0 2 R = ( 1.0 V ) 2 720 Ω = 1.39 × 10 − 3 W P_0=\frac{V_0^2}{R}=\frac{(1.0\,\mathrm{V})^2}{720\,\Omega}=1.39\times10^{-3}\,\mathrm{W} P 0 = R V 0 2 = 720 Ω ( 1.0 V ) 2 = 1.39 × 1 0 − 3 W なので、飽和値の関係 Δ T s = P 0 R k \Delta T_{\mathrm s}=P_0R_k Δ T s = P 0 R k
R k = Δ T s P 0 = 2.0 K × 720 Ω ( 1.0 V ) 2 = 1.44 × 10 3 K / W ≈ 1 × 10 3 K / W R_k=\frac{\Delta T_{\mathrm s}}{P_0}=\frac{2.0\,\mathrm{K}\times720\,\Omega}{(1.0\,\mathrm{V})^2}=1.44\times10^{3}\,\mathrm{K/W}
\ \approx\ 1\times10^{3}\,\mathrm{K/W} R k = P 0 Δ T s = ( 1.0 V ) 2 2.0 K × 720 Ω = 1.44 × 1 0 3 K/W ≈ 1 × 1 0 3 K/W です。時定数は τ = t 1 / 2 / ln 2 = 1.0 / 0.7 = 1.4 s \tau=t_{1/2}/\ln2=1.0/0.7=1.4\,\mathrm{s} τ = t 1/2 / ln 2 = 1.0/0.7 = 1.4 s
C = τ R k = 1.4 s 1.44 × 10 3 K / W = 9.9 × 10 − 4 J / K ≈ 1 × 10 − 3 J / K C=\frac{\tau}{R_k}=\frac{1.4\,\mathrm{s}}{1.44\times10^{3}\,\mathrm{K/W}}=9.9\times10^{-4}\,\mathrm{J/K}
\ \approx\ 1\times10^{-3}\,\mathrm{J/K} C = R k τ = 1.44 × 1 0 3 K/W 1.4 s = 9.9 × 1 0 − 4 J/K ≈ 1 × 1 0 − 3 J/K です。1 桁の有効数字で答えると R k ≈ 1 × 10 3 K / W R_k\approx1\times10^{3}\,\mathrm{K/W} R k ≈ 1 × 1 0 3 K/W C ≈ 1 × 10 − 3 J / K C\approx1\times10^{-3}\,\mathrm{J/K} C ≈ 1 × 1 0 − 3 J/K P 0 R k = 2.06 K P_0R_k=2.06\,\mathrm{K} P 0 R k = 2.06 K τ = 1.44 s \tau=1.44\,\mathrm{s} τ = 1.44 s t = 1 , 2 , 3 , 4 , 5 s t=1,2,3,4,5\,\mathrm{s} t = 1 , 2 , 3 , 4 , 5 s 1.03 , 1.55 , 1.80 , 1.93 , 2.00 K 1.03,\ 1.55,\ 1.80,\ 1.93,\ 2.00\,\mathrm{K} 1.03 , 1.55 , 1.80 , 1.93 , 2.00 K
V ( t ) = 1 2 V 0 e i ω t + ( 複素共役 ) = V 0 cos ω t V(t)=\tfrac12V_0e^{i\omega t}+(\text{複素共役})=V_0\cos\omega t V ( t ) = 2 1 V 0 e iω t + ( 複素共役 ) = V 0 cos ω t
P ( t ) = V ( t ) 2 R = V 0 2 R cos 2 ω t = V 0 2 2 R ( 1 + cos 2 ω t ) P(t)=\frac{V(t)^2}{R}=\frac{V_0^2}{R}\cos^2\omega t=\frac{V_0^2}{2R}\left(1+\cos2\omega t\right) P ( t ) = R V ( t ) 2 = R V 0 2 cos 2 ω t = 2 R V 0 2 ( 1 + cos 2 ω t ) です。非振動成分は V 0 2 / ( 2 R ) V_0^2/(2R) V 0 2 / ( 2 R ) 2 ω 2\omega 2 ω V 0 2 / ( 2 R ) V_0^2/(2R) V 0 2 / ( 2 R )
十分に長い時間が経った定常状態では、Δ T s \Delta T_{\mathrm s} Δ T s T 0 T_0 T 0 T 0 / R k = V 0 2 / ( 2 R ) T_0/R_k=V_0^2/(2R) T 0 / R k = V 0 2 / ( 2 R )
T 0 = R k V 0 2 2 R ⟺ R k = 2 R T 0 V 0 2 T_0=\frac{R_kV_0^2}{2R}\qquad\Longleftrightarrow\qquad R_k=\frac{2RT_0}{V_0^2} T 0 = 2 R R k V 0 2 ⟺ R k = V 0 2 2 R T 0 です。答えは R k = 2 R T 0 / V 0 2 R_k=2RT_0/V_0^2 R k = 2 R T 0 / V 0 2
遅れます。試料は熱容量をもつため、熱が入り始めてもすぐには温度が上がらず、応答が時間的に後ろへずれます。式で言えば、次の設問で見るように振動成分の応答は 1 / ( 1 + 2 i ω R k C ) 1/(1+2i\omega R_kC) 1/ ( 1 + 2 iω R k C ) − ϕ -\phi − ϕ ϕ > 0 \phi>0 ϕ > 0 ω → 0 \omega\to0 ω → 0 ω \omega ω π / 2 \pi/2 π /2
Δ T s ( t ) = T 0 + R e [ T ~ e 2 i ω t ] \Delta T_{\mathrm s}(t)=T_0+\mathrm{Re}\!\left[\tilde T e^{2i\omega t}\right] Δ T s ( t ) = T 0 + Re [ T ~ e 2 iω t ]
( 2 i ω C + 1 R k ) T ~ = V 0 2 2 R ⟹ T ~ = V 0 2 R k 2 R ⋅ 1 1 + 2 i ω R k C \left(2i\omega C+\frac{1}{R_k}\right)\tilde T=\frac{V_0^2}{2R}
\ \Longrightarrow\
\tilde T=\frac{V_0^2R_k}{2R}\cdot\frac{1}{1+2i\omega R_kC} ( 2 iω C + R k 1 ) T ~ = 2 R V 0 2 ⟹ T ~ = 2 R V 0 2 R k ⋅ 1 + 2 iω R k C 1 したがって振幅と位相は
T s a c = ∣ T ~ ∣ = V 0 2 R k 2 R 1 + ( 2 ω R k C ) 2 , tan ϕ = 2 ω R k C T_{\mathrm s}^{\mathrm{ac}}=|\tilde T|=\frac{V_0^2R_k}{2R\sqrt{1+\left(2\omega R_kC\right)^2}},
\qquad
\tan\phi=2\omega R_kC T s ac = ∣ T ~ ∣ = 2 R 1 + ( 2 ω R k C ) 2 V 0 2 R k , tan ϕ = 2 ω R k C です。T ~ = ∣ T ~ ∣ e − i ϕ \tilde T=|\tilde T|e^{-i\phi} T ~ = ∣ T ~ ∣ e − i ϕ ϕ > 0 \phi>0 ϕ > 0 2 ω R k C ≪ 1 2\omega R_kC\ll1 2 ω R k C ≪ 1 T s a c → R k V 0 2 / ( 2 R ) = T 0 T_{\mathrm s}^{\mathrm{ac}}\to R_kV_0^2/(2R)=T_0 T s ac → R k V 0 2 / ( 2 R ) = T 0 ϕ → 0 \phi\to0 ϕ → 0 2 ω R k C ≫ 1 2\omega R_kC\gg1 2 ω R k C ≫ 1 T s a c → V 0 2 / ( 4 ω R C ) T_{\mathrm s}^{\mathrm{ac}}\to V_0^2/(4\omega RC) T s ac → V 0 2 / ( 4 ω R C ) ϕ → π / 2 \phi\to\pi/2 ϕ → π /2
表1 の f f f ω = 2 π f \omega=2\pi f ω = 2 π f 2 ω = 4 π f 2\omega=4\pi f 2 ω = 4 π f
tan ϕ = 4 π f R k C , T s a c 1 + tan 2 ϕ = V 0 2 R k 2 R \tan\phi=4\pi fR_kC,\qquad
T_{\mathrm s}^{\mathrm{ac}}\sqrt{1+\tan^2\phi}=\frac{V_0^2R_k}{2R} tan ϕ = 4 π f R k C , T s ac 1 + tan 2 ϕ = 2 R V 0 2 R k と書けます。第 2 式は振幅と位相を組み合わせると周波数に依らない量になるところが便利です。f = 2.0 H z f=2.0\,\mathrm{Hz} f = 2.0 Hz 1 + 1.49 2 = 1.794 \sqrt{1+1.49^2}=1.794 1 + 1.4 9 2 = 1.794
V 0 2 R k 2 R = 0.418 × 1.794 = 0.750 K ⟹ R k = 2 R × 0.750 K V 0 2 = 2 × 1000 × 0.750 1.0 = 1.5 × 10 3 K / W \frac{V_0^2R_k}{2R}=0.418\times1.794=0.750\,\mathrm{K}
\ \Longrightarrow\
R_k=\frac{2R\times0.750\,\mathrm{K}}{V_0^2}=\frac{2\times1000\times0.750}{1.0}=1.5\times10^{3}\,\mathrm{K/W} 2 R V 0 2 R k = 0.418 × 1.794 = 0.750 K ⟹ R k = V 0 2 2 R × 0.750 K = 1.0 2 × 1000 × 0.750 = 1.5 × 1 0 3 K/W です。f = 10.0 H z f=10.0\,\mathrm{Hz} f = 10.0 Hz 0.100 × 1 + 7.45 2 = 0.752 K 0.100\times\sqrt{1+7.45^2}=0.752\,\mathrm{K} 0.100 × 1 + 7.4 5 2 = 0.752 K
熱容量には高周波側の行が使えます。4 π f R k C ≫ 1 4\pi fR_kC\gg1 4 π f R k C ≫ 1
T s a c ≈ V 0 2 2 R ⋅ 1 4 π f C = V 0 2 8 π f R C T_{\mathrm s}^{\mathrm{ac}}\approx\frac{V_0^2}{2R}\cdot\frac{1}{4\pi fC}=\frac{V_0^2}{8\pi fRC} T s ac ≈ 2 R V 0 2 ⋅ 4 π f C 1 = 8 π f R C V 0 2 なので f T s a c fT_{\mathrm s}^{\mathrm{ac}} f T s ac f = 10 , 20 , 100 H z f=10,20,100\,\mathrm{Hz} f = 10 , 20 , 100 Hz 1.00 , 1.00 , 1.01 K H z 1.00,\ 1.00,\ 1.01\,\mathrm{K\,Hz} 1.00 , 1.00 , 1.01 K Hz f T s a c = 1.005 K H z fT_{\mathrm s}^{\mathrm{ac}}=1.005\,\mathrm{K\,Hz} f T s ac = 1.005 K Hz
C = V 0 2 8 π R ( f T s a c ) = 1.0 8 π × 1000 × 1.005 = 3.96 × 10 − 5 J / K C=\frac{V_0^2}{8\pi R\left(fT_{\mathrm s}^{\mathrm{ac}}\right)}=\frac{1.0}{8\pi\times1000\times1.005}=3.96\times10^{-5}\,\mathrm{J/K} C = 8 π R ( f T s ac ) V 0 2 = 8 π × 1000 × 1.005 1.0 = 3.96 × 1 0 − 5 J/K です。2 桁の有効数字で C ≈ 4.0 × 10 − 5 J / K C\approx4.0\times10^{-5}\,\mathrm{J/K} C ≈ 4.0 × 1 0 − 5 J/K R k ≈ 1.5 × 10 3 K / W R_k\approx1.5\times10^{3}\,\mathrm{K/W} R k ≈ 1.5 × 1 0 3 K/W
検算します。この C , R k C,R_k C , R k τ = R k C = 5.94 × 10 − 2 s \tau=R_kC=5.94\times10^{-2}\,\mathrm{s} τ = R k C = 5.94 × 1 0 − 2 s tan ϕ = 4 π f τ \tan\phi=4\pi f\tau tan ϕ = 4 π f τ f = 2.0 , 10.0 H z f=2.0,10.0\,\mathrm{Hz} f = 2.0 , 10.0 Hz 1.49 , 7.46 1.49,\ 7.46 1.49 , 7.46 1.49 , 7.45 1.49,\ 7.45 1.49 , 7.45 T s a c T_{\mathrm s}^{\mathrm{ac}} T s ac f = 2.0 , 10.0 , 20.0 , 100.0 H z f=2.0,10.0,20.0,100.0\,\mathrm{Hz} f = 2.0 , 10.0 , 20.0 , 100.0 Hz 0.418 , 0.0996 , 0.0501 , 0.0100 K 0.418,\ 0.0996,\ 0.0501,\ 0.0100\,\mathrm{K} 0.418 , 0.0996 , 0.0501 , 0.0100 K 0.418 , 0.100 , 0.0502 , 0.0101 0.418,\ 0.100,\ 0.0502,\ 0.0101 0.418 , 0.100 , 0.0502 , 0.0101 f ≥ 20 H z f\ge20\,\mathrm{Hz} f ≥ 20 Hz tan ϕ \tan\phi tan ϕ tan ϕ = 15 , 75 \tan\phi=15,\ 75 tan ϕ = 15 , 75 R k R_k R k 1.4 × 10 3 K / W 1.4\times10^{3}\,\mathrm{K/W} 1.4 × 1 0 3 K/W C C C 1.4 s 1.4\,\mathrm{s} 1.4 s 5.9 × 10 − 2 s 5.9\times10^{-2}\,\mathrm{s} 5.9 × 1 0 − 2 s 2 2 2 100 H z 100\,\mathrm{Hz} 100 Hz
前半では抵抗を、断面積 A A A l l l m m m q q q n n n i i i
m d v i d t + γ v i = f i m\frac{dv_i}{dt}+\gamma v_i=f_i m d t d v i + γ v i = f i に従います。γ \gamma γ V V V
後半では図2 の光検出器回路を扱います。フォトダイオード PD にバイアス電圧をかけ、角周波数 ω \omega ω
I ( t ) = I 0 + I 1 cos ( ω t ) I(t)=I_0+I_1\cos(\omega t) I ( t ) = I 0 + I 1 cos ( ω t ) が流れます。この電流は、コンデンサ C C C L L L R R R R a R_{\mathrm a} R a R b R_{\mathrm b} R b V O U T V_{\mathrm{OUT}} V OUT
長さ l l l V V V E = V / l E=V/l E = V / l q q q
f i = q E = q V l f_i=qE=\frac{qV}{l} f i = q E = l q V です。定常状態では d v i / d t = 0 dv_i/dt=0 d v i / d t = 0 γ v d = f i \gamma v_{\mathrm d}=f_i γ v d = f i
v d = f i γ = q V γ l v_{\mathrm d}=\frac{f_i}{\gamma}=\frac{qV}{\gamma l} v d = γ f i = γ l q V がドリフト速度です。
減衰力は − γ v i -\gamma v_i − γ v i v i = v d v_i=v_{\mathrm d} v i = v d
( − γ v d ) ⋅ v d = − γ v d 2 = − q 2 V 2 γ l 2 (-\gamma v_{\mathrm d})\cdot v_{\mathrm d}=-\gamma v_{\mathrm d}^{\,2}=-\frac{q^2V^2}{\gamma l^2} ( − γ v d ) ⋅ v d = − γ v d 2 = − γ l 2 q 2 V 2 です。負であることは、この力が粒子から仕事を奪って媒質（格子）へ渡していることを意味します。したがって 1 粒子あたり単位時間に散逸するエネルギーは
γ v d 2 = q 2 V 2 γ l 2 \gamma v_{\mathrm d}^{\,2}=\frac{q^2V^2}{\gamma l^2} γ v d 2 = γ l 2 q 2 V 2 です。定常状態では、電場が単位時間に粒子へ与える仕事 f i v d = q V v d / l f_iv_{\mathrm d}=qVv_{\mathrm d}/l f i v d = q V v d / l
導体中の荷電粒子の総数は n A l nAl n A l
W = n A l ⋅ γ v d 2 = n A l ⋅ q 2 V 2 γ l 2 = n A q 2 V 2 γ l W=nAl\cdot\gamma v_{\mathrm d}^{\,2}=nAl\cdot\frac{q^2V^2}{\gamma l^2}=\frac{nAq^2V^2}{\gamma l} W = n A l ⋅ γ v d 2 = n A l ⋅ γ l 2 q 2 V 2 = γ l n A q 2 V 2 です。一方、電流は電荷密度とドリフト速度と断面積の積で
I = n q v d A = n q A ⋅ q V γ l = n q 2 A V γ l I=nqv_{\mathrm d}A=nqA\cdot\frac{qV}{\gamma l}=\frac{nq^2AV}{\gamma l} I = n q v d A = n q A ⋅ γ l q V = γ l n q 2 A V なので、ジュール熱は
I V = n q 2 A V 2 γ l = W IV=\frac{nq^2AV^2}{\gamma l}=W I V = γ l n q 2 A V 2 = W となり、一致します。示されました。ついでに R = V / I = γ l / ( n q 2 A ) R=V/I=\gamma l/(nq^2A) R = V / I = γ l / ( n q 2 A ) σ = n q 2 / γ \sigma=nq^2/\gamma σ = n q 2 / γ
理想オペアンプなので、反転入力と非反転入力の電圧は等しく（仮想短絡）、入力に電流は流れ込みません。非反転入力は点 P につながっているので V − = V + = V 1 V_-=V_+=V_1 V − = V + = V 1 R a R_{\mathrm a} R a R b R_{\mathrm b} R b V O U T V_{\mathrm{OUT}} V OUT
V 1 = V − = V O U T R a R a + R b V_1=V_-=V_{\mathrm{OUT}}\frac{R_{\mathrm a}}{R_{\mathrm a}+R_{\mathrm b}} V 1 = V − = V OUT R a + R b R a です。したがって
G = V O U T V 1 = R a + R b R a = 1 + R b R a G=\frac{V_{\mathrm{OUT}}}{V_1}=\frac{R_{\mathrm a}+R_{\mathrm b}}{R_{\mathrm a}}=1+\frac{R_{\mathrm b}}{R_{\mathrm a}} G = V 1 V OUT = R a R a + R b = 1 + R a R b です。非反転増幅回路なので G ≥ 1 G\ge1 G ≥ 1
3 素子は並列なのでアドミタンスを足します。
Y ( ω ) = 1 R + i ω C + 1 i ω L = 1 R + i ( ω C − 1 ω L ) Y(\omega)=\frac{1}{R}+i\omega C+\frac{1}{i\omega L}=\frac{1}{R}+i\left(\omega C-\frac{1}{\omega L}\right) Y ( ω ) = R 1 + iω C + iω L 1 = R 1 + i ( ω C − ω L 1 ) よってインピーダンスは
Z ( ω ) = 1 Y ( ω ) = 1 1 R + i ( ω C − 1 ω L ) = R 1 + i R ( ω C − 1 ω L ) Z(\omega)=\frac{1}{Y(\omega)}=\frac{1}{\dfrac{1}{R}+i\left(\omega C-\dfrac{1}{\omega L}\right)}
=\frac{R}{1+iR\left(\omega C-\dfrac{1}{\omega L}\right)} Z ( ω ) = Y ( ω ) 1 = R 1 + i ( ω C − ω L 1 ) 1 = 1 + i R ( ω C − ω L 1 ) R です。虚部が消えて ∣ Z ∣ |Z| ∣ Z ∣ ω C = 1 / ( ω L ) \omega C=1/(\omega L) ω C = 1/ ( ω L )
ω 0 = 1 L C \omega_0=\frac{1}{\sqrt{LC}} ω 0 = L C 1 が共振角周波数で、そこでのインピーダンスは
Z ( ω 0 ) = R Z(\omega_0)=R Z ( ω 0 ) = R です。共振点では L C LC L C
∣ Z ∣ = R / 1 + R 2 ( ω C − 1 / ( ω L ) ) 2 |Z|=R/\sqrt{1+R^2\left(\omega C-1/(\omega L)\right)^2} ∣ Z ∣ = R / 1 + R 2 ( ω C − 1/ ( ω L ) ) 2 ∣ Z ∣ = R / 2 |Z|=R/\sqrt2 ∣ Z ∣ = R / 2
∣ ω C − 1 ω L ∣ = 1 R \left|\omega C-\frac{1}{\omega L}\right|=\frac{1}{R} ω C − ω L 1 = R 1 です。左辺を整理すると
ω C − 1 ω L = C ( ω 2 − ω 0 2 ) ω = C ( ω − ω 0 ) ( ω + ω 0 ) ω \omega C-\frac{1}{\omega L}=\frac{C\left(\omega^2-\omega_0^2\right)}{\omega}=\frac{C(\omega-\omega_0)(\omega+\omega_0)}{\omega} ω C − ω L 1 = ω C ( ω 2 − ω 0 2 ) = ω C ( ω − ω 0 ) ( ω + ω 0 ) で、ω = ω 0 + Δ ω \omega=\omega_0+\Delta\omega ω = ω 0 + Δ ω Δ ω / ω 0 ≪ 1 \Delta\omega/\omega_0\ll1 Δ ω / ω 0 ≪ 1 ( ω + ω 0 ) / ω ≈ 2 (\omega+\omega_0)/\omega\approx2 ( ω + ω 0 ) / ω ≈ 2 2 C Δ ω ≈ ± 1 / R 2C\Delta\omega\approx\pm1/R 2 C Δ ω ≈ ± 1/ R ∣ Δ ω ∣ = 1 / ( 2 R C ) |\Delta\omega|=1/(2RC) ∣Δ ω ∣ = 1/ ( 2 R C )
2 Δ ω = 1 R C 2\Delta\omega=\frac{1}{RC} 2Δ ω = R C 1 です。近似を使わずに 2 次方程式 ω 2 ∓ ω / ( R C ) − ω 0 2 = 0 \omega^2\mp\omega/(RC)-\omega_0^2=0 ω 2 ∓ ω / ( R C ) − ω 0 2 = 0 1 / ( R C ) 1/(RC) 1/ ( R C ) Q Q Q ω 0 / ( 2 Δ ω ) = ω 0 R C = R C / L \omega_0/(2\Delta\omega)=\omega_0RC=R\sqrt{C/L} ω 0 / ( 2Δ ω ) = ω 0 R C = R C / L R C / L ≫ 1 R\sqrt{C/L}\gg1 R C / L ≫ 1
光電流の変動成分は共振しているので、並列共振回路のインピーダンスは Z ( ω 0 ) = R Z(\omega_0)=R Z ( ω 0 ) = R
V 1 信号 ( t ) = I 1 R cos ( ω 0 t ) ⟹ ⟨ ( V 1 信号 ) 2 ⟩ = I 1 2 R 2 2 V_1^{\text{信号}}(t)=I_1R\cos(\omega_0t)
\ \Longrightarrow\
\left\langle\left(V_1^{\text{信号}}\right)^2\right\rangle=\frac{I_1^2R^2}{2} V 1 信号 ( t ) = I 1 R cos ( ω 0 t ) ⟹ ⟨ ( V 1 信号 ) 2 ⟩ = 2 I 1 2 R 2 です（正弦波の二乗平均は振幅の二乗の半分）。雑音側は、測定周波数幅を半値全幅にとるので 2 π Δ f = 2 Δ ω = 1 / ( R C ) 2\pi\Delta f=2\Delta\omega=1/(RC) 2 π Δ f = 2Δ ω = 1/ ( R C ) Δ f = 1 / ( 2 π R C ) \Delta f=1/(2\pi RC) Δ f = 1/ ( 2 π R C )
⟨ δ V 1 2 ⟩ = 4 k B T R ⋅ Δ f = 4 k B T R 2 π R C = 2 k B T π C \left\langle\delta V_1^{\ 2}\right\rangle=4k_{\mathrm B}TR\cdot\Delta f=\frac{4k_{\mathrm B}TR}{2\pi RC}=\frac{2k_{\mathrm B}T}{\pi C} ⟨ δ V 1 2 ⟩ = 4 k B T R ⋅ Δ f = 2 π R C 4 k B T R = π C 2 k B T となります。V O U T V_{\mathrm{OUT}} V OUT G G G
I 1 2 R 2 2 = 2 k B T π C ⟹ I 1 = 2 R k B T π C \frac{I_1^2R^2}{2}=\frac{2k_{\mathrm B}T}{\pi C}
\ \Longrightarrow\
I_1=\frac{2}{R}\sqrt{\frac{k_{\mathrm B}T}{\pi C}} 2 I 1 2 R 2 = π C 2 k B T ⟹ I 1 = R 2 π C k B T が答えです。L L L L L L R R R C C C k B T / C k_{\mathrm B}T/C k B T / C J / F = V 2 \mathrm{J/F}=\mathrm{V^2} J/F = V 2 ⋅ \sqrt{\cdot} ⋅ R R R R = 1 M Ω R=1\,\mathrm{M\Omega} R = 1 MΩ C = 1 n F C=1\,\mathrm{nF} C = 1 nF T = 300 K T=300\,\mathrm{K} T = 300 K I 1 ≈ 2 p A I_1\approx2\,\mathrm{pA} I 1 ≈ 2 pA R R R C C C
陽子同士を衝突させて反陽子をつくる反応
p + p → p + p + p ˉ + p p+p\to p+p+\bar p+p p + p → p + p + p ˉ + p を考えます。終状態は陽子 3 個と反陽子 1 個で、静止質量の総和は 4 m p 4m_p 4 m p c c c m p m_p m p L L L
ローレンツ変換は β = v / c \beta=v/c β = v / c γ = 1 / 1 − β 2 \gamma=1/\sqrt{1-\beta^2} γ = 1/ 1 − β 2
c t ′ = γ ( c t − β x ) , x ′ = γ ( x − β c t ) ct'=\gamma\left(ct-\beta x\right),\qquad x'=\gamma\left(x-\beta\,ct\right) c t ′ = γ ( c t − β x ) , x ′ = γ ( x − β c t ) です。これを代入して計算します。
( c t ′ ) 2 − x ′ 2 = γ 2 [ ( c t − β x ) 2 − ( x − β c t ) 2 ] = γ 2 [ ( c t ) 2 − 2 β x c t + β 2 x 2 − x 2 + 2 β c t x − β 2 ( c t ) 2 ] = γ 2 ( 1 − β 2 ) [ ( c t ) 2 − x 2 ] = ( c t ) 2 − x 2 \begin{aligned}
(ct')^2-x'^2
&=\gamma^2\left[(ct-\beta x)^2-(x-\beta ct)^2\right]\\
&=\gamma^2\left[(ct)^2-2\beta x\,ct+\beta^2x^2-x^2+2\beta ct\,x-\beta^2(ct)^2\right]\\
&=\gamma^2\left(1-\beta^2\right)\left[(ct)^2-x^2\right]
=(ct)^2-x^2
\end{aligned} ( c t ′ ) 2 − x ′2 = γ 2 [ ( c t − β x ) 2 − ( x − β c t ) 2 ] = γ 2 [ ( c t ) 2 − 2 β x c t + β 2 x 2 − x 2 + 2 β c t x − β 2 ( c t ) 2 ] = γ 2 ( 1 − β 2 ) [ ( c t ) 2 − x 2 ] = ( c t ) 2 − x 2 最後で γ 2 ( 1 − β 2 ) = 1 \gamma^2(1-\beta^2)=1 γ 2 ( 1 − β 2 ) = 1 ( c t ) 2 − x 2 (ct)^2-x^2 ( c t ) 2 − x 2 S ′ S' S ′
不変量 s ≡ ( ∑ E ) 2 − ( ∑ p c ) 2 s\equiv(\sum E)^2-\left(\sum \boldsymbol{p}c\right)^2 s ≡ ( ∑ E ) 2 − ( ∑ p c ) 2 s = ( 4 m p c 2 ) 2 s=\left(4m_pc^2\right)^2 s = ( 4 m p c 2 ) 2
(i) 正面衝突では 2 つの陽子の運動量が打ち消すので、∑ E = 2 E p \sum E=2E_p ∑ E = 2 E p ∑ p = 0 \sum\boldsymbol{p}=0 ∑ p = 0 s = ( 2 E p ) 2 s=(2E_p)^2 s = ( 2 E p ) 2
( 2 E p ) 2 ≥ ( 4 m p c 2 ) 2 ⟹ E p ≥ 2 m p c 2 (2E_p)^2\ge\left(4m_pc^2\right)^2\ \Longrightarrow\ E_p\ge2m_pc^2 ( 2 E p ) 2 ≥ ( 4 m p c 2 ) 2 ⟹ E p ≥ 2 m p c 2 なので最小値は E p = 2 m p c 2 E_p=2m_pc^2 E p = 2 m p c 2 m p c 2 m_pc^2 m p c 2
(ii) 標的が静止している場合、入射陽子のエネルギー E p E_p E p p p p E p 2 = ( p c ) 2 + ( m p c 2 ) 2 E_p^2=(pc)^2+\left(m_pc^2\right)^2 E p 2 = ( p c ) 2 + ( m p c 2 ) 2
s = ( E p + m p c 2 ) 2 − ( p c ) 2 = E p 2 + 2 E p m p c 2 + ( m p c 2 ) 2 − [ E p 2 − ( m p c 2 ) 2 ] = 2 m p c 2 ( E p + m p c 2 ) s=\left(E_p+m_pc^2\right)^2-(pc)^2
=E_p^2+2E_pm_pc^2+\left(m_pc^2\right)^2-\left[E_p^2-\left(m_pc^2\right)^2\right]
=2m_pc^2\left(E_p+m_pc^2\right) s = ( E p + m p c 2 ) 2 − ( p c ) 2 = E p 2 + 2 E p m p c 2 + ( m p c 2 ) 2 − [ E p 2 − ( m p c 2 ) 2 ] = 2 m p c 2 ( E p + m p c 2 ) これが ( 4 m p c 2 ) 2 = 16 ( m p c 2 ) 2 \left(4m_pc^2\right)^2=16\left(m_pc^2\right)^2 ( 4 m p c 2 ) 2 = 16 ( m p c 2 ) 2
2 m p c 2 E p + 2 ( m p c 2 ) 2 ≥ 16 ( m p c 2 ) 2 ⟹ E p ≥ 7 m p c 2 2m_pc^2E_p+2\left(m_pc^2\right)^2\ge16\left(m_pc^2\right)^2
\ \Longrightarrow\
E_p\ge7m_pc^2 2 m p c 2 E p + 2 ( m p c 2 ) 2 ≥ 16 ( m p c 2 ) 2 ⟹ E p ≥ 7 m p c 2 となり、最小値は E p = 7 m p c 2 E_p=7m_pc^2 E p = 7 m p c 2 6 m p c 2 ≈ 5.6 G e V 6m_pc^2\approx5.6\,\mathrm{GeV} 6 m p c 2 ≈ 5.6 GeV m p c 2 m_pc^2 m p c 2 s ∝ E p \sqrt{s}\propto\sqrt{E_p} s ∝ E p 6.2 G e V 6.2\,\mathrm{GeV} 6.2 GeV
a) 電荷が保存しません。始状態の電荷は 2 e 2e 2 e ( + 1 ) + ( + 1 ) + ( − 1 ) + 0 = + 1 (+1)+(+1)+(-1)+0=+1 ( + 1 ) + ( + 1 ) + ( − 1 ) + 0 = + 1 e e e 1 + 1 = 2 1+1=2 1 + 1 = 2 1 + 1 − 1 + 1 = 2 1+1-1+1=2 1 + 1 − 1 + 1 = 2
b) バリオン数が保存しません。始状態のバリオン数は 2 2 2 p , p , p ˉ , π + p,p,\bar p,\pi^+ p , p , p ˉ , π + + 1 , + 1 , − 1 , 0 +1,+1,-1,0 + 1 , + 1 , − 1 , 0 1 1 1 π \pi π − 1 -1 − 1 2 e 2e 2 e ( + 1 ) + ( + 1 ) + ( − 1 ) + ( + 1 ) = + 2 (+1)+(+1)+(-1)+(+1)=+2 ( + 1 ) + ( + 1 ) + ( − 1 ) + ( + 1 ) = + 2
一様磁場中で運動量に垂直な磁場を受ける荷電粒子は円運動し、円運動の式は相対論的にも ∣ d p d t ∣ = ∣ q ∣ v B \left|\dfrac{d\boldsymbol{p}}{dt}\right|=|q|vB d t d p = ∣ q ∣ v B ∣ d p d t ∣ = p v r \left|\dfrac{d\boldsymbol{p}}{dt}\right|=\dfrac{pv}{r} d t d p = r p v
p = ∣ q ∣ B r p=|q|Br p = ∣ q ∣ B r となります（p = γ m v p=\gamma mv p = γ m v e e e
p c = e B r c = ( 1.602 × 10 − 19 C ) × ( 1.0 T ) × ( 0.50 m ) × ( 3.0 × 10 8 m / s ) = 2.40 × 10 − 11 J pc=eBrc=(1.602\times10^{-19}\,\mathrm{C})\times(1.0\,\mathrm{T})\times(0.50\,\mathrm{m})\times(3.0\times10^{8}\,\mathrm{m/s})=2.40\times10^{-11}\,\mathrm{J} p c = e B r c = ( 1.602 × 1 0 − 19 C ) × ( 1.0 T ) × ( 0.50 m ) × ( 3.0 × 1 0 8 m/s ) = 2.40 × 1 0 − 11 J です。e V \mathrm{eV} eV e e e e e e
p c = e B r c e e V = B r c e V = 1.0 × 0.50 × 3.0 × 10 8 = 1.5 × 10 8 e V pc=\frac{eBrc}{e}\,\mathrm{eV}=Brc\,\mathrm{eV}=1.0\times0.50\times3.0\times10^{8}=1.5\times10^{8}\,\mathrm{eV} p c = e e B r c eV = B r c eV = 1.0 × 0.50 × 3.0 × 1 0 8 = 1.5 × 1 0 8 eV となります。したがって
p = 0.15 G e V / c p=0.15\ \mathrm{GeV}/c p = 0.15 GeV / c です。有効数字 2 桁で 0.15 G e V / c 0.15\,\mathrm{GeV}/c 0.15 GeV / c p [ G e V / c ] = 0.3 B [ T ] r [ m ] p[\mathrm{GeV}/c]=0.3\,B[\mathrm{T}]\,r[\mathrm{m}] p [ GeV / c ] = 0.3 B [ T ] r [ m ] B = 1.0 B=1.0 B = 1.0 r = 0.50 r=0.50 r = 0.50
光電子増倍管では、光電面で放出された光電子が第1ダイノードへ加速され、そこで数個の二次電子を叩き出し、それが第2ダイノードへ加速されてさらに増える、という増倍が段ごとに繰り返されます。したがって回路に要求されるのは、光電面から各ダイノードを経て陽極に向かって電位が単調に高くなっていること、そして各段の電位差がほぼ等しいことです。
高電圧源が − 1 k V -1\,\mathrm{kV} − 1 kV − 1 k V -1\,\mathrm{kV} − 1 kV 0 V 0\,\mathrm{V} 0 V 1000 / 12 ≈ 83 V 1000/12\approx83\,\mathrm{V} 1000/12 ≈ 83 V − 1000 V -1000\,\mathrm{V} − 1000 V − 917 V -917\,\mathrm{V} − 917 V 0 V 0\,\mathrm{V} 0 V
信号は最も電位の高い陽極から取り出します。陽極を直流的に接地電位に置いたまま、結合コンデンサを介して信号を取り出すか、陽極と接地の間に負荷抵抗を入れてその両端の電圧を読みます。分圧器に流す定常電流は、信号によるパルス電流よりも十分大きく（電極電位がパルス中に動かないよう）とる必要があり、後段のダイノードには電荷を供給するためのバイパスコンデンサを並列に入れるのが普通です。
p c = 0.3 G e V pc=0.3\,\mathrm{GeV} p c = 0.3 GeV E = ( p c ) 2 + ( m c 2 ) 2 E=\sqrt{(pc)^2+\left(mc^2\right)^2} E = ( p c ) 2 + ( m c 2 ) 2 β = v / c = p c / E \beta=v/c=pc/E β = v / c = p c / E
反陽子（m c 2 = 0.9 G e V mc^2=0.9\,\mathrm{GeV} m c 2 = 0.9 GeV
E = 0.3 2 + 0.9 2 = 0.90 = 0.949 G e V , β p ˉ = 0.3 0.949 = 0.316 ≈ 0.3 E=\sqrt{0.3^2+0.9^2}=\sqrt{0.90}=0.949\,\mathrm{GeV},\qquad
\beta_{\bar p}=\frac{0.3}{0.949}=0.316\approx0.3 E = 0. 3 2 + 0. 9 2 = 0.90 = 0.949 GeV , β p ˉ = 0.949 0.3 = 0.316 ≈ 0.3 π − \pi^- π − m c 2 = 0.15 G e V mc^2=0.15\,\mathrm{GeV} m c 2 = 0.15 GeV
E = 0.3 2 + 0.15 2 = 0.1125 = 0.335 G e V , β π = 0.3 0.335 = 0.894 ≈ 0.9 E=\sqrt{0.3^2+0.15^2}=\sqrt{0.1125}=0.335\,\mathrm{GeV},\qquad
\beta_{\pi}=\frac{0.3}{0.335}=0.894\approx0.9 E = 0. 3 2 + 0.1 5 2 = 0.1125 = 0.335 GeV , β π = 0.335 0.3 = 0.894 ≈ 0.9 です。答えは反陽子が光速の約 0.3 0.3 0.3 π − \pi^- π − 0.9 0.9 0.9
飛行時間は t = L / ( β c ) t=L/(\beta c) t = L / ( β c ) 1 / β = E / ( p c ) 1/\beta=E/(pc) 1/ β = E / ( p c )
Δ t = L c ( 1 β p ˉ − 1 β π ) = 0.3 m 3.0 × 10 8 m / s ( 3.162 − 1.118 ) = 1.0 × 10 − 9 × 2.044 s \Delta t=\frac{L}{c}\left(\frac{1}{\beta_{\bar p}}-\frac{1}{\beta_\pi}\right)
=\frac{0.3\,\mathrm{m}}{3.0\times10^{8}\,\mathrm{m/s}}\left(3.162-1.118\right)
=1.0\times10^{-9}\times2.044\ \mathrm{s} Δ t = c L ( β p ˉ 1 − β π 1 ) = 3.0 × 1 0 8 m/s 0.3 m ( 3.162 − 1.118 ) = 1.0 × 1 0 − 9 × 2.044 s すなわち
Δ t ≈ 2 × 10 − 9 s = 2 n s \Delta t\approx2\times10^{-9}\,\mathrm{s}=2\,\mathrm{ns} Δ t ≈ 2 × 1 0 − 9 s = 2 ns です。個別には反陽子が 3.2 n s 3.2\,\mathrm{ns} 3.2 ns π − \pi^- π − 1.1 n s 1.1\,\mathrm{ns} 1.1 ns β = 0.3 , 0.9 \beta=0.3,\ 0.9 β = 0.3 , 0.9 3.3 − 1.1 = 2.2 n s 3.3-1.1=2.2\,\mathrm{ns} 3.3 − 1.1 = 2.2 ns 2 n s 2\,\mathrm{ns} 2 ns
b) の検出器が適切です。図3 の横軸のスケールが判断材料になります。a) は 200 n s 200\,\mathrm{ns} 200 ns n s \mathrm{ns} ns 200 n s 200\,\mathrm{ns} 200 ns 10 n s 10\,\mathrm{ns} 10 ns 20 20 20 30 n s 30\,\mathrm{ns} 30 ns n s \mathrm{ns} ns n s \mathrm{ns} ns
設問6 で見たように識別すべき飛行時間差は 2 n s 2\,\mathrm{ns} 2 ns n s \mathrm{ns} ns 2 n s 2\,\mathrm{ns} 2 ns 0.1 V 0.1\,\mathrm{V} 0.1 V 0.5 V 0.5\,\mathrm{V} 0.5 V 1 n s 1\,\mathrm{ns} 1 ns
出典: 東京大学大学院理学系研究科 物理学専攻 令和2年度 修士課程 入学試験問題 物理学。問題文は要約して引用しています。
© 2026 夢現技研合同会社 ・Feeding the text to an LLM is welcome. Code samples are MIT licensed.