問題文
速度選択器と質量分析器における軌道収束とビームプロファイルの広がり 問題文
質量 m m m 、電荷 q q q (q > 0 q > 0 q > 0 ) のイオンを分析する装置を考える。空間には直交座標系 ( x , y , z ) (x, y, z) ( x , y , z ) が設定されている。
装置は以下の2つの領域からなる。
【領域1(速度選択器)】 0 ≤ x ≤ L 0 \le x \le L 0 ≤ x ≤ L
この領域には、y y y 軸正の向きに一様な電場 E E E 、z z z 軸正の向きに一様な磁場 B 1 B_1 B 1 が存在する。
x = L x=L x = L の平面には、y = 0 y=0 y = 0 を中心として y y y 方向に微小な幅 w w w (− w 2 ≤ y ≤ w 2 -\frac{w}{2} \le y \le \frac{w}{2} − 2 w ≤ y ≤ 2 w )を持つスリットが開いている。これ以外の場所に入射した粒子は遮蔽される。
【領域2(質量分析器)】 x > L x > L x > L かつ y > 0 y > 0 y > 0
この領域には電場はなく、z z z 軸負の向きに一様な磁場 B 2 B_2 B 2 が存在する。
また、x = L x=L x = L かつ y > 0 y>0 y > 0 の平面上にはイオンの到達位置を記録する検出器が配置されている。
原点 ( 0 , 0 , 0 ) (0,0,0) ( 0 , 0 , 0 ) から x x x 軸正の向きに、様々な初速度 v v v を持つイオンのビームが入射する。
以下の手順で、検出器上に到達するイオンの分布幅を求めよ。ただし、領域1を通過する時間は十分に短く、粒子が磁場から受ける影響は小さいため、領域1内での y y y 方向の速度成分 v y v_y v y は v y ≪ v v_y \ll v v y ≪ v を満たし、x x x 方向の速度成分は入射時の速度 v v v で一定であると近似してよい。
初速度 v = v c v = v_c v = v c のイオンは、領域1で軌道が曲がることなく直進してスリットの中心 ( L , 0 , 0 ) (L, 0, 0) ( L , 0 , 0 ) を通過した。v c v_c v c を E , B 1 E, B_1 E , B 1 を用いて表せ。
初速度 v = v c + Δ v v = v_c + \Delta v v = v c + Δ v (ただし ∣ Δ v ∣ ≪ v c |\Delta v| \ll v_c ∣Δ v ∣ ≪ v c )のイオンが領域1に入射した。このイオンが x = L x=L x = L に到達したときの y y y 座標 y ( L ) y(L) y ( L ) を、m , q , B 1 , v c , L , Δ v m, q, B_1, v_c, L, \Delta v m , q , B 1 , v c , L , Δ v を用いて表せ。
イオンがスリットを通過できる初速度のずれの最大値を Δ v m a x > 0 \Delta v_{max} > 0 Δ v ma x > 0 とする(すなわち − Δ v m a x ≤ Δ v ≤ Δ v m a x -\Delta v_{max} \le \Delta v \le \Delta v_{max} − Δ v ma x ≤ Δ v ≤ Δ v ma x の粒子がスリットを抜ける)。Δ v m a x \Delta v_{max} Δ v ma x を m , q , B 1 , v c , L , w m, q, B_1, v_c, L, w m , q , B 1 , v c , L , w を用いて表せ。
初速度 v = v c + Δ v v = v_c + \Delta v v = v c + Δ v でスリットを通過したイオンは、領域2に入り半円軌道を描いて検出器に到達する。領域2に入射する際の速度の向きは x x x 軸から微小な角度ずれているが、この微小角の2次以上の影響を無視すると、検出器上での到達位置の y y y 座標 y f i n a l y_{final} y f ina l は、スリット通過時の位置 y ( L ) y(L) y ( L ) と初速度成分 v v v を用いて y f i n a l ≈ y ( L ) + 2 m v q B 2 y_{final} \approx y(L) + \frac{2mv}{qB_2} y f ina l ≈ y ( L ) + q B 2 2 m v と近似できる。
この関係を用いて、y f i n a l y_{final} y f ina l を m , q , B 1 , B 2 , v c , L , Δ v m, q, B_1, B_2, v_c, L, \Delta v m , q , B 1 , B 2 , v c , L , Δ v のみを用いて表せ。
スリットを通過したイオン群が検出器上で広がる y y y 座標の幅(最大値と最小値の差)を δ y \delta y δ y とする。後述の制約条件の数値を代入し、δ y \delta y δ y の値を求めよ。
制約
各変数は以下の現実的なSI単位系の値を持つ。
イオンの質量: m = 1.6 × 10 − 25 kgm = 1.6 \times 10^{-25} \, \text{kg} m = 1.6 × 1 0 − 25 kg (約 100 u100 \, \text{u} 100 u の重イオン)
イオンの電荷: q = 1.6 × 10 − 19 Cq = 1.6 \times 10^{-19} \, \text{C} q = 1.6 × 1 0 − 19 C
領域1の電場: E = 3.8 × 10 2 V/mE = 3.8 \times 10^2 \, \text{V/m} E = 3.8 × 1 0 2 V/m
領域1の磁場: B 1 = 1.0 × 10 − 2 TB_1 = 1.0 \times 10^{-2} \, \text{T} B 1 = 1.0 × 1 0 − 2 T
領域2の磁場: B 2 = 1.0 TB_2 = 1.0 \, \text{T} B 2 = 1.0 T
領域1の長さ: L = 4.0 × 10 − 1 mL = 4.0 \times 10^{-1} \, \text{m} L = 4.0 × 1 0 − 1 m
スリットの幅: w = 3.0 × 10 − 4 mw = 3.0 \times 10^{-4} \, \text{m} w = 3.0 × 1 0 − 4 m
入力形式
δ y \delta y δ y の値を μ m\mu\text{m} μ m (マイクロメートル、10 − 6 m10^{-6} \, \text{m} 1 0 − 6 m )単位で計算し、その数値を自然数で入力せよ。
解説
1. 速度選択器の直進条件
領域1において、速度 v v v で x x x 軸正の向きに進む電荷 q q q の粒子は、y y y 軸正の向きに電場による力 F E = q E F_E = qE F E = q E を受け、フレミングの左手の法則より y y y 軸負の向きにローレンツ力 F B = q v B 1 F_B = qvB_1 F B = q v B 1 を受ける。
直進するためにはこれらの力が釣り合う必要があるため、
q E = q v c B 1 ⟹ v c = E B 1 qE = qv_c B_1 \implies v_c = \frac{E}{B_1} q E = q v c B 1 ⟹ v c = B 1 E
2. 領域1における軌道のずれの導出
速度 v = v c + Δ v v = v_c + \Delta v v = v c + Δ v で入射した粒子に対する y y y 方向の運動方程式を立てる。
x x x 方向の速度成分は v x ≈ v v_x \approx v v x ≈ v で一定と近似できるため、粒子が受ける y y y 方向の力は
F y = q E − q v B 1 = q E − q ( v c + Δ v ) B 1 = q E − q v c B 1 − q Δ v B 1 F_y = qE - qvB_1 = qE - q(v_c + \Delta v)B_1 = qE - qv_cB_1 - q\Delta v B_1 F y = q E − q v B 1 = q E − q ( v c + Δ v ) B 1 = q E − q v c B 1 − q Δ v B 1
q E = q v c B 1 qE = qv_cB_1 q E = q v c B 1 であるから、F y = − q B 1 Δ v F_y = -qB_1 \Delta v F y = − q B 1 Δ v となる。
したがって、y y y 方向の加速度 a y a_y a y は一定であり、
a y = − q B 1 m Δ v a_y = -\frac{qB_1}{m}\Delta v a y = − m q B 1 Δ v
初期条件 y ( 0 ) = 0 , y ˙ ( 0 ) = 0 y(0)=0, \dot{y}(0)=0 y ( 0 ) = 0 , y ˙ ( 0 ) = 0 より、t t t 秒後の y y y 座標は等加速度運動の公式より
y ( t ) = 1 2 a y t 2 = − q B 1 2 m Δ v ⋅ t 2 y(t) = \frac{1}{2} a_y t^2 = -\frac{qB_1}{2m}\Delta v \cdot t^2 y ( t ) = 2 1 a y t 2 = − 2 m q B 1 Δ v ⋅ t 2
粒子が x = L x=L x = L に到達するまでの時間 t L t_L t L は、t L ≈ L v ≈ L v c t_L \approx \frac{L}{v} \approx \frac{L}{v_c} t L ≈ v L ≈ v c L と近似できる。これを代入すると、
y ( L ) = − q B 1 2 m Δ v ( L v c ) 2 = − q B 1 L 2 2 m v c 2 Δ v y(L) = -\frac{qB_1}{2m}\Delta v \left(\frac{L}{v_c}\right)^2 = -\frac{qB_1 L^2}{2mv_c^2} \Delta v y ( L ) = − 2 m q B 1 Δ v ( v c L ) 2 = − 2 m v c 2 q B 1 L 2 Δ v
3. スリットを通過できる速度幅
粒子がスリットを通過する条件は、− w 2 ≤ y ( L ) ≤ w 2 -\frac{w}{2} \le y(L) \le \frac{w}{2} − 2 w ≤ y ( L ) ≤ 2 w である。
これを解くと、
− w 2 ≤ − q B 1 L 2 2 m v c 2 Δ v ≤ w 2 -\frac{w}{2} \le -\frac{qB_1 L^2}{2mv_c^2} \Delta v \le \frac{w}{2} − 2 w ≤ − 2 m v c 2 q B 1 L 2 Δ v ≤ 2 w
∣ Δ v ∣ ≤ m v c 2 w q B 1 L 2 |\Delta v| \le \frac{mv_c^2 w}{qB_1 L^2} ∣Δ v ∣ ≤ q B 1 L 2 m v c 2 w
したがって、通過できる初速度のずれの最大値 Δ v m a x \Delta v_{max} Δ v ma x は
Δ v m a x = m v c 2 w q B 1 L 2 \Delta v_{max} = \frac{mv_c^2 w}{qB_1 L^2} Δ v ma x = q B 1 L 2 m v c 2 w
4. 領域2における到達位置 y f i n a l y_{final} y f ina l
問題文で与えられた近似式 y f i n a l ≈ y ( L ) + 2 m v q B 2 y_{final} \approx y(L) + \frac{2mv}{qB_2} y f ina l ≈ y ( L ) + q B 2 2 m v に、これまでに求めた y ( L ) y(L) y ( L ) と v = v c + Δ v v = v_c + \Delta v v = v c + Δ v を代入する。
y f i n a l = − q B 1 L 2 2 m v c 2 Δ v + 2 m ( v c + Δ v ) q B 2 y_{final} = -\frac{qB_1 L^2}{2mv_c^2} \Delta v + \frac{2m(v_c + \Delta v)}{qB_2} y f ina l = − 2 m v c 2 q B 1 L 2 Δ v + q B 2 2 m ( v c + Δ v )
y f i n a l = 2 m v c q B 2 + ( 2 m q B 2 − q B 1 L 2 2 m v c 2 ) Δ v y_{final} = \frac{2mv_c}{qB_2} + \left( \frac{2m}{qB_2} - \frac{qB_1 L^2}{2mv_c^2} \right) \Delta v y f ina l = q B 2 2 m v c + ( q B 2 2 m − 2 m v c 2 q B 1 L 2 ) Δ v
第一項は速度 v c v_c v c の粒子の理想的な到達位置であり、第二項が初速度のずれ Δ v \Delta v Δ v による到達位置の変動を表す。
5. ビームの広がり幅 δ y \delta y δ y と数値計算
到達位置の広がり幅 δ y \delta y δ y は、Δ v \Delta v Δ v が − Δ v m a x -\Delta v_{max} − Δ v ma x から Δ v m a x \Delta v_{max} Δ v ma x まで変化したときの y f i n a l y_{final} y f ina l の変動幅である。
δ y = ∣ 2 m q B 2 − q B 1 L 2 2 m v c 2 ∣ × 2 Δ v m a x = ∣ 2 m q B 2 − q B 1 L 2 2 m v c 2 ∣ 2 m v c 2 w q B 1 L 2 \delta y = \left| \frac{2m}{qB_2} - \frac{qB_1 L^2}{2mv_c^2} \right| \times 2\Delta v_{max} = \left| \frac{2m}{qB_2} - \frac{qB_1 L^2}{2mv_c^2} \right| \frac{2mv_c^2 w}{qB_1 L^2} δ y = q B 2 2 m − 2 m v c 2 q B 1 L 2 × 2Δ v ma x = q B 2 2 m − 2 m v c 2 q B 1 L 2 q B 1 L 2 2 m v c 2 w
展開して整理すると、
δ y = ∣ 4 m 2 v c 2 q 2 B 1 B 2 L 2 − 1 ∣ w \delta y = \left| \frac{4m^2 v_c^2}{q^2 B_1 B_2 L^2} - 1 \right| w δ y = q 2 B 1 B 2 L 2 4 m 2 v c 2 − 1 w
ここで、v c = E B 1 v_c = \frac{E}{B_1} v c = B 1 E を用いると、
δ y = ∣ 4 m 2 E 2 q 2 B 1 3 B 2 L 2 − 1 ∣ w \delta y = \left| \frac{4m^2 E^2}{q^2 B_1^3 B_2 L^2} - 1 \right| w δ y = q 2 B 1 3 B 2 L 2 4 m 2 E 2 − 1 w
各変数に制約の数値を代入する。まず、係数となる無次元量 K = 4 m 2 E 2 q 2 B 1 3 B 2 L 2 K = \frac{4m^2 E^2}{q^2 B_1^3 B_2 L^2} K = q 2 B 1 3 B 2 L 2 4 m 2 E 2 を計算する。
K = 4 × ( 1.6 × 10 − 25 ) 2 × ( 3.8 × 10 2 ) 2 ( 1.6 × 10 − 19 ) 2 × ( 1.0 × 10 − 2 ) 3 × 1.0 × ( 4.0 × 10 − 1 ) 2 K = \frac{4 \times (1.6 \times 10^{-25})^2 \times (3.8 \times 10^2)^2}{(1.6 \times 10^{-19})^2 \times (1.0 \times 10^{-2})^3 \times 1.0 \times (4.0 \times 10^{-1})^2} K = ( 1.6 × 1 0 − 19 ) 2 × ( 1.0 × 1 0 − 2 ) 3 × 1.0 × ( 4.0 × 1 0 − 1 ) 2 4 × ( 1.6 × 1 0 − 25 ) 2 × ( 3.8 × 1 0 2 ) 2
電場 E E E 以外の項をまとめると 4 ( 4.0 × 10 2 ) 2 \frac{4}{(4.0 \times 10^2)^2} ( 4.0 × 1 0 2 ) 2 4 の形になるため、計算は以下のように簡略化できる。
K = 4 × ( 3.8 × 10 2 4.0 × 10 2 ) 2 = 4 × ( 3.8 4.0 ) 2 = 4 × ( 0.95 ) 2 = 4 × 0.9025 = 3.61 K = 4 \times \left( \frac{3.8 \times 10^2}{4.0 \times 10^2} \right)^2 = 4 \times \left( \frac{3.8}{4.0} \right)^2 = 4 \times (0.95)^2 = 4 \times 0.9025 = 3.61 K = 4 × ( 4.0 × 1 0 2 3.8 × 1 0 2 ) 2 = 4 × ( 4.0 3.8 ) 2 = 4 × ( 0.95 ) 2 = 4 × 0.9025 = 3.61
よって、広がり幅 δ y \delta y δ y は次のように求まる。
δ y = ∣ 3.61 − 1 ∣ w = 2.61 w \delta y = |3.61 - 1| w = 2.61w δ y = ∣3.61 − 1∣ w = 2.61 w
w = 3.0 × 10 − 4 mw = 3.0 \times 10^{-4} \, \text{m} w = 3.0 × 1 0 − 4 m を代入すると、
δ y = 2.61 × 3.0 × 10 − 4 m= 7.83 × 10 − 4 m\delta y = 2.61 \times 3.0 \times 10^{-4} \, \text{m} = 7.83 \times 10^{-4} \, \text{m} δ y = 2.61 × 3.0 × 1 0 − 4 m = 7.83 × 1 0 − 4 m
これを μ m\mu\text{m} μ m 単位に変換する(1 m= 10 6 μ m1 \, \text{m} = 10^6 \, \mu\text{m} 1 m = 1 0 6 μ m )。
δ y = 7.83 × 10 − 4 × 10 6 μ m= 783 μ m\delta y = 7.83 \times 10^{-4} \times 10^6 \, \mu\text{m} = 783 \, \mu\text{m} δ y = 7.83 × 1 0 − 4 × 1 0 6 μ m = 783 μ m
最終的な自然数の答えは 783 である。