平成 16 年度大学院入学試験問題 I (３時間)注意(1) 問題 I-1、I-2、I-3 の解答はそれぞれ別の用紙１枚に記入せよ（裏面を用いてもよい）。
(2) 各解答用紙ごとに、指定の位置に問題番号、受験番号、氏名を記入せよ。
——————————————————————————————————

I-1 (弦の振動) (100 点)
長さ 2l 離れた、水平面内にある固定点間に張られた弦について考える。両固定点を通る方
向を x 軸にとり、固定点の中央を x = 0 とする。鉛直下方を y 軸にとり、時刻 t での弦の形
状を y(x, t) で表す。両端では y(±l, t) = 0 とする。弦の単位長さあたりの質量を m、張力
を T 、重力加速度を g とおく。|y(x, t)| ¿ l の場合には、y(x, t) は方程式2∂ 2 y(x, t)
2 ∂ y(x, t)
−c
=g
∂t2
∂x2

に従う。ここで、c =(A)q

T /m である。

(1) 振動していない平衡状態では、弦は重力のためにわずかに下方にたれている。このとき、弦の平衡形状 y0 (x) はgy0 (x) = 2 (l2 − x2 )(B)cと表されることを示せ。

(2) 問 (1) の平衡状態にある弦の中央をつまんで、微小量 h だけ下に引っ張ったときの弦の形状を y1 (x) とする。u1 (x) ≡ y1 (x) − y0 (x) を求めよ。
この状態で弦を静かに放すと、弦は振動を始める。y 方向の変位を u(x, t) ≡ y(x, t) − y0 (x)として、以下の問に答えよ。

(3) この振動にあらわれる固有振動の振動数を求めよ。
(4) 初期条件を満足する u(x, t) を求めよ。
(5) この振動における、もっとも低い固有振動のモ−ドと、その次に低い固有振動のモ−
ドの振幅の２乗の比を求めよ。

I-2（数学）(100 点)(1) 関数 y(x) についての以下の二つの微分方程式の一般解を求めよ。なお、
dy
d2 y
y0 ≡
、y 00 ≡ 2 の意味である。
dx
dx
(a)

y 00 + 2y 0 − 3y = 9e2x(b)x2 y 00 − 6y = 5x4(2) Bernoulli 数 Bn (n = 0, 1, 2, . . .) を以下の無限級数展開で定義する。
∞
X
x
Bn n
=
x
ex − 1 n=0 n!

このとき、 B0 、B1 、B2 を求めよ。 また、

B3 = B5 = B7 = · · · = 0を証明せよ。

(3) cot x の無限級数展開、
cot x =

∞
1X
An n
x
x n=0 n!

の係数 An を Bernoulli 数を用いて書け。

(4) cos(kx) の x ∈ [−π, π] における Fourier 展開、
½

¾

∞
C0 X
+
Cn cos(nx) + Sn sin(nx)cos(kx) =2n=1において、 係数 Cn と Sn を求めよ。

I-3 (交換相互作用) (100 点)
二つの原子に電子の軌道 1, 2 がそれぞれ局在している。その軌道部分の規格化された波動関数を φ1 (x), φ2 (x) とし、それらは互いに直交するものとする。今 φ1 , φ2にそれぞれ電子が１つ入っている。また 2 つの電子のスピン演算子を、 s1 と s2
(1)(2)
とし、スピン部分の波動関数を ψsz , ψsz とする。

(1) 全スピン演算子を S = s1 + s2 、全スピンを S 、その磁気量子数を Sz としたとき、その規格化された固有関数 ΨS,Sz をすべて求めよ。また電子がフェルミ統計に従うことを考慮し、S = 1, S = 0 のそれぞれについて、許される軌道部分の規格化された波動関数を書け。ただし二つの電子の座標を x1 と
x2 とする。
(2) 演算子 P1 , P0 を
P1 Ψ1,Sz = Ψ1,Sz , P1 Ψ0,0 = 0,
P0 Ψ0,0 = Ψ0,0 , P0 Ψ1,Sz = 0を満たすものとする。P0 , P1 を Pi = ci + di s1 · s2 という形で書け。ただしci , di は実数である。

(3) 電子間にクーロン相互作用が働く。このハミルトニアンを H = e2 /r12 とする。ここで r12 ≡ |x1 − x2 | である。このハミルトニアンの期待値を S = 0,
S = 1 の両方の場合について計算し、それがhHi = a + bhs1 · s2 i|Sの形にまとまり、 b がZ
e2
b = −2 d3 x1 d3 x2 φ∗1 (x2 )φ∗2 (x1 ) φ1 (x1 )φ2 (x2 )
r12で与えられることを示せ。ただし、hs1 · s2 i|S は、 全スピンが S の場合の
s1 · s2 の固有値であり、a, b は共に S によらない量とする。

(4) 系の体積を Ω としたとき
1 X
φ∗1 (r)φ2 (r) = √aq eiq·rΩ qと展開できることを用い、 b が負であることを示せ。ただしZZ
iq·r−²r
iq·r
3 e
3 e
≡ lim
dr
dr
²→+0
|r|
|r|
とする。

(5) S = 1, S = 0 のうちエネルギーが低い状態はどちらか。

平成 16 年度大学院入学試験問題 II (3 時間 30 分)
注意(1) 問題 II-1、II-2、II-3 の解答はそれぞれ別の用紙 1 枚に記入せよ。（裏面を用いてもよい。）(2) 各用紙ごとに、指定の位置に問題番号、受験番号、氏名を記入せよ。

II-1 (格子振動) (100 点)
下の図の様に N 個の同じ質量 m を持つ質点が同じ強さ k のばね (N + 1) 本でつながれている場合を考えよう。ただし、両端のばねは固定されている。

¨¨¨
¥
¨
¥
¨
¥¥¥ ~¨¨¨¥¨¥¨¥¥¥ ~¨¨¨¥¨¥¨¥¥¥ p p p p p p p p p p p p
§§§
¦
§
¦
§
¦¦¦
§§§
¦
§
¦
§
¦¦¦
§
§§
¦
§
¦
§
¦¦¦

¨
¨¨
¥
¨
¥
¨
¥¥¥ ~¨¨¨¥¨¥¨¥¥¥ ~¨¨¨
¥
¨
¥
¨
¥
¥¥
§
§§
¦
§
¦
§
¦¦¦
§§
§
¦
§
¦
§
¦¦¦
§§§
¦
§
¦
§
¦
¦¦

(1) 左から n 番目の質点の任意の時刻 t における変位を xn (t) として、系のラグ
ランジアンを書け。ただし、右の方向を正にとる。
(2) xn (t) を第 n 成分とする N 次元ベクトルを x(t) と書く。 系の運動方程式を
d2
x = Lxdt2と表すとき、N 次の対称行列 L を書け。
(p)(3) すべての質点が同じ角振動数で振動する運動 x(p)n (t) = an cos ωp t (n = 1, 2, . . . , N )を基準振動という。独立な基準振動は N 個ある（p = 1, 2, . . . , N ）。基準角
振動数 ωp と振幅 a(p)n はそれぞれωp = Ω sin

θp
,
2

a(p)
n = A sin nθp(A)qと与えられることを示せ。ただし、Ω = 2 k/m、θp = pπ/(N + 1)、そして
A は任意定数である。
）（ヒント：a(p)n = a sin(nθ + α) と置いて、境界条件から θ と α を求めよ。
(p)

(p)

(p)
は次のように正規(4) a(p) = t (a1 , a2 , . . . , a(p)n ) とベクトル表示するとき、aq直交系になっていることを示せ。ただし、式 (A) において A = 2/(N + 1)と選ぶものとする。
N
X

0

(p )
a(p)
= δpp0n an(p, p0 = 1, 2, . . . , N )

n=1
N
X

ここで、δpp0 はクロネッカーのデルタである。また、
である。

n=1

sin2 nθp = (N + 1)/2

(5) 各質点にさらに速度に比例する抵抗 −2mγ ẋn （ただし、γ > 0； 定数）と外力 Fn (t) = mfn cos ωe t が加わっている場合を考える。このときの運動方程式はd2dx − Lx + 2γ x = f cos ωe t2dtdtとなる。ただし、f = t (f1 , f2 , . . . , fN )。基準座標 qp (t) ≡

(B)
N
X
n=1

う運動方程式を f (p) =

N
X

a(p)
n fn を使って表せ。

n=1(6) 十分時間が経過したときの qp (t) を求めよ。

a(p)n xn (t) の従II-2 (静電容量) (100 点)(1) 図１のような同軸円筒コンデンサーがある。内側の円筒極板の半径は a、外側の円筒極板の半径は b、軸方向の長さは L であり、両円筒極板の厚みは無視できるものとする。この同軸コンデンサーの静電容量
C0 を求めよ。ただし、両円筒極板間は真空であり、真空の誘電率を ²0 とする。またコンデンサーは十分長く、両端で電場の乱れる効果は無視してよい。

a
L
b
図１

(2) このコンデンサーの軸を鉛直方向に保持しながら、下端を大きな容器にためた非導電性の液体表面に接 α'Cするように設置する。このとき極板間の間隔が十分
狭ければ、液面の極板との接触線 C の接線方向に働αθ0く表面張力 α （図２）（注）の反作用 α によって、液面は高さ h1 だけ上昇する（毛管現象）（図３）。こ
図２の h1 を上昇した液体に働く重力 (重力加速度を g とせよ) と α0 との釣り合いを考えることによって求めよ。ただし、液面形状の水平面からの変形は高さ h1と比べて無視できるものとし、接触線 C での接触角は常に θ であるものとする。また液体の密度は ρ で与えられ、液体より上にある気体の密度は無視する。
（注）二つの媒質が接する界面（表面）を変形させる時に必要な仕事は、体積変化に伴う仕事 (＝体積変化×圧力差) と表面積を変化させるのに伴う仕事(＝表面積変化×表面エネルギー α) に分類される。
この表面エネルギーは表面張力とも呼ばれ（エネルギー÷面積）の次元を持つ量である。

図３

V(3) このときのコンデンサーの静電容量 C1 を測定した。
液の高さ h1 を C1 より求めよ。ただし液体の誘電率は ²、気体の誘電率は ²0 と書けるものとする。

h2h1(4) 極板間に電位差 V を与えると液面はさらに h2 だけ上昇した（図４）。この状態の静電容量 C2 を用いて
h2 および液体の密度 ρ を求めよ。
(5) コンデンサーを液体中に完全に沈めたときの静電容量 C3 を測定すれば、液体の誘電率 ² を知ることができる。このことを利用して a, b, C0 , C1 , C2 , C3 , V ,²0 , θ を用いて表面張力 α を表せ。

h1

図４

II-3 (Ising モデル) (100 点)(1) 体積 V の結晶の全原子数を N とする。各原子が磁気モーメント µ を持ち、
磁場 H に対して平行 (+) か反平行 (−) のいずれかの状態のみを取るものとする。各格子点の状態を表す変数を σi (i = 1, ..., N ) とし、磁場 H に対して平行な時 +1、反平行な時 −1 とする。また隣り合うスピンの間のみに相互作用が働くものとする。この時の系全体の Hamiltonian は
H = −JX
hiji

σi σj − µH

X

σi

(A)

i

と書ける。ただし、hiji は隣り合うスピンに対する和である。ここで i 番目の格子点にあるスピンの Hamiltonian Hi を σi の平均値 hσi を用いて次の形で書けると仮定する。

Hi = −Jzhσiσi − µHσi

(B)

ここで z は最隣接原子の数である。また Jzhσi = µH 0 と置けば、一つのスピンには有効磁場 Heff = H + H 0 が働いているとみなすことができる。

(a) 一つのスピンに対する分配関数 Z1 を Heff を用いて書け。
(b) 各スピンが独立であるとして、系全体についての Helmholtz の自由エネ
ルギー F を求めよ。更にエントロピー S と内部エネルギー U の表式を導け。
(c) この系の磁化の強さ M (= Nµhσi) が次の式を満たすことを示せ。
V
M=

³ µH
Nµ
JzV M ´
tanh
+
V
kT
kT µN

(C)

(d) 式 (C) が H = 0 で M = 0 以外の解を持つ条件を求めよ。ここで、Curie
温度 Tc = Jz/k （k は Boltzmann 定数）を用いよ。
(2) 同じことを統計力学的に考える。原子の総数を N とし、磁場 H に対して平行 (+)、反平行 (−) にあるスピンの数をそれぞれ N+ 、N− とする。またパラ
メータ Φ を N+ /N = (1 + Φ)/2、N− /N = (1 − Φ)/2 により定義する。隣接原子対（総数 zN/2）を +, − で区別し、隣接原子の両方が磁場に平行である
対の数を N++ 、以下同様に N−− 、N+− と置くと、この時の相互作用エネル
ギー U はU = −J(N++ + N−− − N+− )(D)となる。ここで一つの格子点が + または − により占められる確率 p+ =N+ /N 、p− = N− /N を導入し、U の式の中の N++ 等をその平均値

1z
N++ = zN+ p+ = N (1 + Φ)228z
N+− = zN+ p− = N (1 − Φ2 )41z
N−− = zN− p− = N (1 − Φ)2
2
8

(E)
(F)
(G)

で置き換える近似を行う。

(a) N 個のスピンのうち、N+ 個が平行スピン、N− 個が反平行スピンである
ときの系全体のエントロピーを Φ を用いて書き表せ。ただし、Stirlingの公式 (log x! ' x log x − x) を用いること。
(b) 磁場 H をゼロとしたときの Helmholtz の自由エネルギー F を Φ を用いて書け。
(c) Φ の平衡値を求める式が、式 (C) で H = 0 としたときと同じ形になることを示せ。またパラメータ Φ の物理的意味を説明せよ。
(d) 定積比熱 Cv の Tc 近傍での振舞いがどうなるか、グラフを用いて説明せよ。

平成 16 年度大学院入学試験問題 III（3 時間）
注意(1) 問題は III-1 から III-8 まで 8 問ある。これから 3 問選択せよ。
(2) 選択した問題の回答はそれぞれ別の用紙一枚に記入せよ（裏面を用いてもよい）。
(3) 各用紙ごとに、左上に問題番号、右上に受験番号と氏名を記入せよ。

III-1 (選択)(振り子の実験) (100 点)
剛体を水平な固定軸で支えた実体振り子で重力加速度 g を求める実験を行うものとする（Borda の振り子）。この振り子は図１のように細い針金で吊られた半径r、質量 M の金属球からなり、支持体のナイフエッジ K を支点として振動する。
針金の長さを l、振り子が最大振幅になるときの角度を α、振動の周期を T とすると、重力加速度 g は次の式で求めることができる。
"

4π 2 (l + r)
2r2
g=
1
+
T2
5(l + r)2

#"

α2
1+
8

#

(A)

K

l

L

(1) 金属球の直径をノギス（キャリパー）を用いて測定したところ、図 2 のような表示だった。直径の値を誤差を含めて書け。

(2) 振り子の周期 T を測定するため、振り子を 190 回振動させて 10 回ごとのラップタイムを測った。表はその結果である。1 列目は 0 から 90 回目、2 列目は100 から 190 回目のラップタイムで、3 列目は 1 列目と 2 列目の差にあたる
100 周期分の時間である。
回数
0
10
20
30
40
50
60
70
80
90

ラップタイム1 分 40 秒 5
2’01.22’22.02’43.1
3’03.8
3’24.83’45.84’06.64’27.6
4’48.6

回数
100
110
120
130
140
150
160
170
180
190

ラップタイム
5 分 09 秒 7
5’30.7
5’53.7
6’14.56’35.5
6’56.5
7’17.3
7’38.3
7’59.5
8’20.5

100 回ごとの時間
3 分 29 秒 2=209.2 秒
3’29.5=209.6
3’31.7=211.7
3’29.4=209.4
3’29.7=209.7
3’29.7=209.7
3’29.5=209.5
3’29.7=209.7
3’29.9=209.9
3’29.9=209.9

このデータを用いて、平均周期 T とその誤差 ∆T を求めよ。測定データの取り扱いには十分注意すること。また、∆T は平均値の標準偏差より求め、
q
q
4.32
0.43
= 0.219089,= 0.077281 を用いて良い。
90
72
長尺ノギスを用いた測定により、支点 K から金属球の下端までの長さ L は 110.96±
0.01 cm であった。また周期の測定前後において求めた振り子の振幅を求めそれらの幾何平均をとることにより、 α2 = 1.54 × 10−3 を得た。

(3) 針金の長さ l はいくらか。誤差を含めて示せ。
(4) g の計算式を書き下せ。ただし、数値計算は行わなくて良い。
(5) 問 (4) の計算を電卓を用いて行ったところ、979.60184 と言う値が出た。g の値はいくらか。
(6) l, r, T の誤差を ∆l、∆r、∆T としたとき、g の相対誤差が次の形で近似できることを示せ。
¯¯¯¯¯¯¯ ∆g ¯¯ ∆L + ∆r ¯
¯ ∆T ¯
¯
¯
¯
¯
¯
¯
¯
¯=¯
+ 2¯
¯ g ¯
L−r ¯
T ¯

また、有効数字を考えて g の誤差 ∆g を求めよ。

(B)III-2 (選択)(回路) (100 点)
超伝導電磁石の「永久モード」と「永久モードスイッチ」に関する以下の記述を読んで、下記の問いに答えよ。
一般に常伝導線材でできた電磁石に大電流を流すと、線材の電気抵抗に伴う自己発熱のため、発生磁界の限界が小さくなる。超伝導転移点以下に冷却された超伝導電磁石ではこのような問題はなく、小型の電磁石で強磁界を発生できる。また、超伝導物質のみで作られた閉回路に流れる定常電流は減衰しない。この特質を利用したものが超伝導電磁石の永久モードである。超伝導電磁石の巻き線の両端を常伝導物質を介さずに接続しておけば、外部電源から電流を供給し続けなくても、電流および磁界を高精度に維持することができる。しかしながら、当然、超伝導体の閉回路中の電磁石に外部から電流を流し込むことはできない。このままでは超伝導体となった時に流れていた電流が維持されるだけである。永久モード
スイッチがこの問題点を打開する。超伝導電磁石の両端を結ぶ超伝導線材の一部に、局所的に熱を加えて超伝導転移温度以上に加熱保持できる仕組みが取り付けられてある。この部分が超伝導状態の時は永久モード、常伝導の時には外部から電流を流し込める通常モードとなり、あたかも永久モードのスイッチのように働く。

(1) 超伝導線を巻いた直径 10 cm、長さ 50 cm の巻き線の厚みの無視できる空芯
ソレノイド形の電磁石がある。この磁石に 100 A の電流を流した時、中心部
で 10 T の磁束密度が得られた。この磁石のインダクタンス L はいくらか。
ただし、このソレノイドは十分細長く、無限長近似を適用してよいものとする。ソレノイド内部の空間の透磁率としては µ0 = 4π×10−7 kg·m·C−2 を用いよ。
(2) 永久モードスイッチの常伝導部分の線材の電気抵抗率が 1×10−8 Ω·m、線材の半径 0.1 mm、長さ 1 m とした時の電気抵抗を R とする。今、超伝導電磁石に流れる電流が 0 の状態の時に永久モードスイッチのヒーターに熱を入れ通常モードにした。時刻 t=0 にこの電磁石に接続した電源より電流 I0（一定値）を流しはじめたとする。時刻 t において磁石中心部における磁束密度 B
を I0 , L, R を用いて求め、時間変化の様子を L, R の数値を用いてグラフに示せ。
(3) 問 (2) の状態で十分時間が経って定常状態になったところで永久モードにすれば、その時ソレノイドに流れる電流を永久に保持することになる。しかし、現実的な超伝導磁石では、超伝導の完全な閉回路を作ることは難しい。技術的な理由で超伝導線は多数本の超伝導体の細線を束ねて銅線の中に埋め込んだ構造をとっており、超伝導線の閉回路を作ろうとしても接続部分において超伝導体同士の間にこの銅が挟まれてしまうからである。しかしながら、ソレノイドの両端からの二本の超伝導線を十分な長さだけ平行にハンダ付けす
ることにより接続部分の電気抵抗 r を小さくし、擬似的に永久モード動作をさせることができる。今、永久モードにしてから電源を切り離した状態で、磁場が１時間あたり 1 ppm(=1×10−6 ) の精度で安定していることを要求するには接続部分の電気抵抗 r の値の上限はいくらか。

(4) 永久モードスイッチのついた超伝導電磁石では問 (2) のように外部から流した電流が即電磁石に流れる電流にはならないことに注意しなければならない。今、超伝導電磁石に流れる電流が 0 の状態で永久モードスイッチのヒー
ターに熱を入れ、通常モードにしてあるものとする。時刻 t=0 にこの電磁石に接続した電源より時間に対して線形に変化する電流 I0 =βt を流し始めたとする。時刻 t において磁石中心部の磁束密度 B を I0 , L, R, β を用いて求め、
時間変化の様子を β=0.1 A/sec の場合について L, R の数値を用いてグラフに示せ。
(5) 100A の電流を流して永久モードにしてからスイッチ (S1) で電源を切り離した状態の超伝導電磁石がある (図 1)。この超伝導電磁石の永久モードスイッ
チ (PSW) を誤操作により通常モードにしてしまった。このとき、電磁石に貯えられていた磁気エネルギーは常伝導になった部分でジュール熱として発散される。この発熱で超伝導転移温度を超える領域が周辺の線材部分に拡大し、磁気エネルギーはすべて熱となってしまう。このとき局所的に著しく発熱することとなるため、最悪の事態では線材の融点を超えて電磁石巻き線を焼損することになる。このような事態を防ぐため、通常は保護回路 (D1,D2,R1
よりなる) が取り付けられている。この保護回路の動作を説明せよ。ただし
D1,D2 の電圧電流特性は図 2 に示されているものを使用するとし、R1 の抵
抗値は 0.02 Ω とする。また、この保護回路を付加することによって超伝導電磁石の使用条件に制約を受けることになるが、その条件を説明せよ。
I

⿥
વ
ዉ
㔚
⏛
⍹

60A

S1
2
5
9

D2

R1
D1

㔚
Ḯ
0.6V

図１

1.2V

図２

V

III-3 (選択)(調和振動子) (100 点)
一次元調和振動子のハミルトニアンは

H=

´
1³ 2
p + ω 2 x22である。ここで、x, p は座標および運動量の演算子を表わし、ω は角振動数である。また
質量は簡単のため 1 にとった。系が量子力学に従うとして

a = (2ω)−1/2 (p − iωx)a† = (2ω)−1/2 (p + iωx)と表わされる消滅演算子 a、生成演算子 a† を導入する。ただし、h̄ = 1 とした。

(1) a と a† の交換関係を求めよ。また H はµ1
H =ω a a+2¶†のように書き換えられることを示せ。

(2) 系の基底状態 |0i は式
a|0i = 0を満たす。この式から基底状態の波動関数 ϕ0 (x) = hx|0i についての微分方程式を導き、その規格化された解を求めよ。なお、必要ならば次のガウス積分を用いてよい。
Z ∞
−∞

r
−ax2

e

dx =

π
a

(3) 第一励起状態 |1i = a† |0i の規格化した波動関数 ϕ1 (x) = hx|1i を求めよ。

次に、二次元の調和振動子を考える。そのハミルトニアン H0 は

H0 =

´
´
1 ³
1³ 2
px + p2y + ω 2 x2 + y 222と与えられる。この系に·¸1
H = −λ (~p · ~e) − (p2x + p2y )202の摂動が加わったとする。ただし、p~ = (px , py ) であり、~e は単位ベクトル (cos θ, sin θ) を表わしている。

(4) θ = π/4 のとき、基底状態のエネルギ−を H 0 の２次の摂動項まで求めよ。

摂動が加わらないときには第一励起状態には２重の縮退がある。その縮退は H 0 によりとれる。

(5) θ = 0 のとき、エネルギ−の低い方の状態を表わす波動関数とそのエネルギ−を一次の摂動により求めよ。λ > 0 と λ < 0 の２つの場合に分けて考えよ。なお、波動関数の答えは ϕ0 (x), ϕ0 (y), ϕ1 (x), ϕ1 (y) を用いて書くこと。
(6) 一般の θ の場合にエネルギ−の低い方の状態を表わす波動関数とそのエネルギ−を一次の摂動により求めよ。
(7) 問 (5) の結果を考慮して、問 (6) の結果がなぜそのようになったか説明せよ。

III-4 (選択) (電磁誘導) (100 点)
透磁率 µ0 の真空中に、z 軸に沿った無限に長いまっすぐな導線と、それから距離 x(> a) を隔てて半径 a の円形一巻きコイルとが x − z 平面内にある (右図)。以下の問に答
えよ。






(1) 導線に電流 I1 、コイルに電流 I2 を図の向きに流す。
このとき、電流 I1 によって図の P 点の位置に生じる磁界の向きと大きさ H を答えよ。
(2) 問 (1) の磁界によって、P 点の線要素 ds = adθ を流れる電流に作用する力のx 軸方向成分 dFx と z 軸方向成分 dFz を求めよ。
(3) コイル全体に作用する力の方向と大きさ F を求めよ。必要ならば、
Z π
0

π
dθ
=√ 2
u − t cos θu − t2(u > t)を用いてよい。

(4) この後、一旦両方の電流を切り、導線に I1 = I sin ωt と時間変化する電流を流した。導線とコイルの相互インダクタンス M を求め、このときコイルに発生する起電力 e (時計周りを正にとる) を求めよ。
(5) 次に導線に I1 = I の定常電流を流し、コイルを x 軸の正の方向に一定速度 v0で動かした。このとき、コイルが場所 x(> a) を通過する時点でのコイルに流れる電流 I2 の向きを答え、大きさを求めよ。コイルの抵抗を R とし、また自己インダクタンスは無視できるものとする。
(6) 問 (5) の場合、コイルを x = d1 (> a) から x = d2 まで動かしたときの仕事 Wを求めよ。またこの仕事が、発生するジュール熱 Q と等しいことを示せ。

III-5 (選択) (量子統計) (100 点)
ボゾン、フェルミオンといった量子力学的粒子は、状態の占有に関して古典的なMaxwell-Boltzmann 粒子とは大きく異なったルールに従っている。このため、さま
ざまな熱力学的性質にも統計性の効果が顔を出してくる。以下では、状態方程式を
例にとり、
1 2
p
(A)
ε(p) =
2mのような分散を持って３次元空間中を運動する質量 m の Maxwell-Boltzmann 粒子からなる気体（古典気体）と量子力学的粒子の気体とが、どのように異なっているかを考察してみよう。以下の問題では、粒子間の相互作用は無視するものとし、プランク定数を h（あるいは ~）、ボルツマン定数を kB などとするものとする。

(1) まず、古典気体を考えよう。体積 V の箱に閉じこめられ、温度 T で熱平衡にある一個の古典粒子の分配関数 Z1 を計算せよ。ただし、必要があれば Gauss
積分公式:
Z ∞
√
2
dx e−x = π−∞を用いてよい。

(2) 上の問題と同じ設定で、同種粒子が N 個あるときの分配関数 ZN も計算せよ。
(ただし、気体分子の位置は空間内に固定されていないので、エネルギー、エントロピーなどの物理量が正しく示量的になるように、分配関数をある因子で割る必要があることに注意せよ。)(3) この系の粒子数 N が変化できるようにする。化学ポテンシャルを µ として、大分配関数 Ξ を計算せよ。
（Ω = −pV = −kB T ln Ξ はグランドポテンシャル）(4) 熱力学関係式 N = − ∂Ω∂µにより、古典気体の化学ポテンシャル µ = µcl を粒子数 N 、温度 T 、体積 Vの関数として定めると、·µ ¶ ¸Nµcl = kB T ln
λ3TVとなることを示し、状態方程式を導け。ただし、熱的ドブロイ波長 λT は、
s
h2
λT =
2πmkB Tで与えられる。

(5) 次に量子力学的粒子の場合を考えよう。古典粒子と同じ質量 m、式 (A) で与えられる同一の分散 ε(p) を持つボゾン、フェルミオンについて、そのグラン
ドポテンシャルがZ ∞ze−βε2 V m3/2
dε ε3/2Ω=− √(B)
3 2π 2 ~3 0
1 − gze−βεで与えられることを示せ。ただし、β = kB1T で、z = eβµ（µ は化学ポテンシャ
ル）はフガシティ、統計性を表す定数 g は、ボゾンのとき+1、フェルミオン
のとき −1 をとるものとする。

(6) 状態方程式を求めるには、式 (B) の右辺に現れる z を、N 、V 、T などで表す必要がある。以下、フガシティ z = eβµ ¿ 1 の状況を考えよう。式 (B) の右
辺を z について２次まで展開し、N = − ∂Ωを用いて、z と N とを結びつける∂µ
方程式を作れ。ただし、最終的な結果は、積分を全て実行した上で、λT 、V 、
z 、N 、g だけを用いてできるだけ簡潔な形で表せ。ただし、必要ならば Γ 関
数の公式:
Z ∞
√
Γ(z) =
dx xz−1 e−x , Γ(z + 1) = zΓ(z) , Γ(1/2) = π0を用いてよい。

(7) 今の方程式で条件: z = eβµ → 0 を仮定して、z について最低次の項のみを残すと、問 (2) で計算した古典気体のフガシティz = exp(βµcl ) が得られることを示し、量子気体が古典気体で近似できるためには粒子数密度に対してどのような条件を課す必要があるか述べよ。
(8) 問 (6)、(7) の結果を考慮すると、フガシティ z を、微小パラメータ: x ≡ N
λ3
V Tについての巾級数の形に求めると良いことが示唆される。問 (6) で得られた方程式の解 z を x の巾級数として求めよ（x の２次まででよい）。
(9) また、物理量: NpVを古典近似の次のオーダーまで計算し、ボゾン、フェルkB T
ミオンの両方について、圧力 p が量子効果により古典気体からどのようにずれるかを調べよ。さらに、その結果を粒子の統計性と関係付けながら物理的に解釈せよ。

III-7 (選択)(天文学:恒星) (100 点)
恒星からの定常的な球対称質量放出（恒星風）について考える。恒星風は理想気体からなり、温度一定としよう。恒星風中では重力と圧力勾配による力のみが働くものとし、また、
質量放出量は小さく、恒星本体の質量は一定とみなせるものとする。恒星本体の質量を M 、
恒星中心からの距離が r である点での恒星風の圧力、密度、速度をおのおの p(r), ρ(r), v(r)として以下の問に答えよ。

(1) この恒星風における質量保存の式と運動方程式を書け。
(2) 恒星風中の速度 v は
2

2a
− GM
1 dv
r
r2
= 2
v dr
v − a2で記述されることを示せ。ここで、a ≡q(A)p/ρ は音速、G は重力定数である。

式 (A) の解のうち、右辺の分母がゼロとなる点を通る解は、その点で分子もゼロにならなければならない。右辺の分子、分母ともにゼロとなる点を臨界点 (critical point) と呼ぶ。
右図に、距離を臨界点の位置 rcrit 、速度を音速 a を単位にして、式 (A) のすべての解を示してある。臨界点を通る解には、距離 r が大きくなるにつれて、亜音速から超音速になる
解と超音速から亜音速になる解があるが、実際の恒星風では前者が実現している。

(3) v と r に対して分離型の式となっていることを利用して、式 (A) を積分せよ。
(4) 問 (3) の積分形から、臨界点を通る解がどのように得られるか説明せよ。
(5) 太陽風は高温のコロナから流れ出す恒星風である。太陽風を等温とみなした場合、臨界点の位置は太陽半径の何倍か、また臨界点での流出速度はいくらか、以下の数値を用いて答えよ。
重力定数 6.7 × 10−8 dyn cm2 g−2
ガス定数 8.3 × 107 erg deg−1 mol−1
平均分子量 0.60太陽質量 2.0 × 1033 g
太陽半径 7.0 × 1010 cm
恒星風の温度 1.0 × 106 K

III-8 (選択) (天文学：太陽) (100 点)(注意) この問題の解答は、問題冊子の最終ページに綴じてある解答用紙を用いること。
皆既日食で見られる太陽コロナの温度が百万度を越える高温であることは、1940 年代に入るまで誰にも気付かれなかった。コロナの分光スペクトルに見られる緑色輝線や赤色輝線などは、実験室では非常に弱いために長い間同定できなかったので、未知の元素の線と考えられ、コロナに因んでコロ
ニュームの線と呼ばれていた。これらの輝線が高階電離の Fe XIV イオンや Fe X イオ
ンから出る輝線であることが判るきっかけとなったのは、実験室における極紫外線領域のデータであった。Edlen によって同定された極紫外線領域の輝線 (図 1) に着目し 図１: Fe X のエネルギーレベルと輝線波長て、Grotorian はこの Fe X イオンの基底２重項のレベル間の遷移に対応する波長が、ちょうどコロニュームの赤色輝線の波長と合致することを 1939 年に初めて示したのである。
(1) 図１の基底２重項レベル間の遷移による赤色輝線の波長 λR を計算せよ。
(2) これら Fe X や Fe XIV イオンの電離ポテンシャルの平均値が約 300eV であることを用いて、コロナの電離平衡温度の近似値を有効数字 2 桁で計算せよ。但し、
1eV=1.60×10−12 erg、ボルツマン定数 k=1.38×10−16 erg·K−1 を用いよ。
皆既日食の時以外でも、コロナグラフを 3000m 級の高山に設置すれば、地上でコロナを観測できる。図２はコロナグラフの原理を示している (光軸方向に圧縮して作図している)。
L1 は対物レンズで、L2 はリレーレンズ、L3 はコロナ像を最終焦点面 F 上に再結像するレンズである。この光学系では、L1 によって出来た明るい太陽本体の像を遮光円錐 C で除去して、コロナからの光だけが通過するようにしている。破線は太陽中心からの光線が遮光円錐 C の中心に結像し、実線はコロナ中のある 1 点からの光線が点 A に結像していることを示している。3 本ずつ描かれている破線と実線のうち、両端の２本はそれぞれ対物レンズ L1 の縁より少し内側を通る光線を示しており、真中の線はレンズ L1 の中心を通る光線を示している。また、L2 によって、対物レンズ L1 の像を絞り D (真中を円形にくり抜いた板) の位置に作り、対物レンズの縁やレンズ枠から生ずる散乱光を、この絞りによって除去している。

(3) 図２で光線は途中の L2 の位置までしか書かれていないが、これを延長して、コロナ像上の 1 点 A が焦点面 F 上の A’ 点に結像するまでの光路を解答用紙の図に実線で示せ。この際、L1 の縁より少し内側から出た光は、絞り D の開口の端を通過することに留意せよ。
コロナグラフを用いてコロナの分光スペクトルを観測すると、例えば Fe XIV の 5303Å 緑
色輝線は、図３のように波長方向に広がった輪郭を示しており、この広がりからもコロナの温度を求めることが出来る。

図2
温度 T の熱平衡にある Fe XIV イオンプラズマの速度分布は等方的な Maxwell 分布をしていると考えられるが、今視線方向のみの運動を考えると、FeXIV イオン粒子が v と v + dv との間の速度を持つ確率は、µ

m
f (v)dv =
2πkT

¶1/2

2

e−mv /(2kT ) dv

で与えられる。但し、ここで m は鉄イオンの質量、k はボル
ツマン定数である。

(4) この分布における Fe XIV イオン粒子の平均運動エネル
ギーは 12 kT で与えられることを示せ。
図３の Fe XIV 5303Å 輝線の輪郭の幅を広げる主たる要因は、望遠鏡の視線方向に沿って存在する Fe XIV イオン粒子の速度分布であると考えることが出来る。すなわち、視線方向に 図 3: Fe XIV 5303Å 輝線
沿って運動している Fe XIV イオン粒子から放射される光は、の輪郭
ドップラー効果によって、中心波長 λ0 =5303Å から運動速度に比例した量だけずれた波長で観測されるので、図３のようなガウス分布の線輪郭となる。

(5) 輝線の強度は図３のように輝線中心から外側に向けて弱くなるが、中心強度の 1/eにまで強度が落ちた波長での線輪郭幅 4λD をドップラー幅と呼んでいる。このドップラー幅が温度のみによって決まっているとした時は、それらの関係は次の式で与えられることを示せ。但し c は光速である。
λ0
4λD =
c

Ã

2kT
m

!1/2

(6) 核融合によって生成された太陽中心部の熱エネルギーは、放射と対流によって運ばれて、太陽表面 (光球) では約 6 千度の熱力学的平衡状態が保たれている。同じ物理過程で、光球の外側に広がっているコロナを百万度まで加熱することができないことは明らかであるにもかかわらず、コロナが百万度に達しているのはなぜか。光球からコロナに、どの様な物理過程によってエネルギーが輸送されてコロナを加熱しているかを答えよ。