NumRoom02 の履歴(No.3)

L1
表題 (C1-60)。標準出力またはターミナル上にそのまま書き出されます。C61 以降は無視されます。
L2-L4
プログラム中で使用するファイルを指定する部分です。
L5
人間用の見出しです。
L6
ユーザのコメント (C1-50)。C51 以降は無視されます。
L7 Z1, Z2, T を、それぞれ 10桁の F 型実数で指定します。
L8
人間用の見出しです。
L9
計算結果を出力するファイル機番、角度の最小値、同最大値、同刻みを、それぞれ 10桁の F 型実数で指定します。角度に関する数値は、degree 単位で指定します。
L10
桁数を見やすくするための補助行です。

↑

出力の説明†

↑

rscat.outlist†

標準出力ファイルです。

↑

rscat.dat†

左が散乱角 θ (単位: degree) で、右が微分断面積 dσ/dΩ (単位: mb/sr) です。rscat.cnt を用いて計算した結果 rscat.dat が、図2.2 の実線でプロットされているものです。

↑

注意点と補足†

このプログラムの構造は、第1章で使用した hsphere と基本的に同じです。hsphere の説明と共通の部分は適当に省略しています。

0° ではラザフォード散乱の微分断面積は発散します*1。プログラムでは、θ が 10^-10 よりも小さい場合には、θ = 10^-10 (単位: degree) として計算します。これは「0 の代わりに極めて小さな値を用いて発散を防ぐ」という処理ですので、このとき、微分断面積の計算結果にはそれほど意味はありません (小さな値としてどういった数字を入れるかで答えが変わるから)。

↑

発展課題†

入力パラメータを変更し、図2.2 にプロットされている破線と点線の結果を算出してください。
公開しているプログラムでは、Z1 に負の値が代入されるとエラーになります。これと同様の処理を、Z2 と T に対しても行うようにしてください。
- ところで、Z1 あるいは Z2 が負のとき、エラーとする必要はあるのでしょうか?
教科書の式は、座標原点にある標的粒子が入射粒子よりもはるかに重い状況を想定しています。入力データに対して、この仮定が成立しているかどうかを判定する処理を追加しましょう。たとえば Z2 が Z1 の 10倍よりも小さいとき、エラーを返すように書き換えてみてください (粒子の質量を入力していないので、ここでは原子番号を質量の指標として利用します)。
トムソンのレーズンパン模型を模倣する計算を行い、微分断面積を求めてみましょう。具体的には、次のようにします。
- 標的粒子を半径 1Å の球体とし、+Z2 e の電荷 (e は電気素量) が、球全体に一様に分布しているものとする。
- 衝突径数が b のとき、標的粒子の正電荷のうち、半径 b の球体の内側にあるものだけが散乱に関与するものとする。すなわち、b ごとに Z2 が変化するものとし、ラザフォードの公式において Z2 --> Z2(b) と修正する。
- 標的粒子内の電子の影響は、ラザフォード散乱と同様に無視する。

↑

参考文献†

[4] F. T. Baker, A. Scott, R. C. Styles, T. H. Kruse, K. Jones, and R. Suchannek, Nucl. Phys. A 351, 63
　　(1981).

[5] W. Karcz, I. Kluska and Z. Sanok, Acta Phys. Pol. B3, 525 (1972).

*1 0° 散乱における断面積の定義を見直せば発散は避けられますが、ここでは立ち入りません。教科書の付録 A を参照してください。