パラメータ設定
同位体選択
崩壊連鎖
プリセット
理論メモ
$N(t) = N_0 \, e^{-\lambda t}$
$\lambda = \dfrac{\ln 2}{T_{1/2}}$
$A = \lambda N \; [\text{Bq}]$
崩壊連鎖 (Bateman):
$\dfrac{dN_B}{dt} = \lambda_A N_A - \lambda_B N_B$
実務応用: 医療用I-131は甲状腺がん治療（T½=8日で短期間に崩壊）、C-14は年代測定（T½=5730年）、Cs-137は食品汚染モニタリング（T½=30年）に関係します。
原子ドット可視化（青=未崩壊、赤=崩壊済み）— 確認時刻でのスナップショット
放射性崩壊・半減期シミュレーターとは
🧑🎓
半減期って、教科書だと「原子数が半分になる時間」って書いてあるけど、実際にどう減っていくのかイメージが湧かないんです。このシミュレーターで何がわかるんですか？
🎓
ざっくり言うと、時間とともに原子が「指数関数的」に減っていく様子が目で見てわかるんだ。例えば、上の「同位体選択」で「C-14」を選んで「初期原子数 N₀」を100万個に設定して「実行」ボタンを押してみて。5730年（半減期）ごとに、グラフ上の原子数がキレイに半分、また半分になっていくのが確認できるよ。
🧑🎓
え、そうなんですか！C-14は半減期が長すぎて実感わかないかも…。もっと身近な例はありますか？
🎓
実務では医療用の放射性医薬品で半減期の短いものが多いね。例えば「I-131」を選んでみて。半減期は約8日だ。シミュレーターの「時間単位」を「日」に変えて、「表示時間幅」を半減期の5倍くらいにすると、1ヶ月ちょっとでほとんど原子が崩壊しきるのがわかる。これが治療に使われる理由で、体に取り込まれた放射性物質が短期間で消えてくれるんだ。
🧑🎓
なるほど！ところで、親が崩壊して子ができる「崩壊系列」って何ですか？グラフの「バテマン方程式を解く」にチェックを入れると、線が2本出てきました。
🎓
良いところに気づいたね。例えば「Co-60」は、崩壊して「Ni-60」になるけど、その過程でガンマ線を出す。この親核種（Co-60）と子核種（Ni-60）の量が時間とともにどう変わるかが、バテマン方程式で計算できるんだ。パラメータの「Bの半減期」を変えると、子核種が溜まっていく様子や、親子で平衡状態になる様子が観察できるよ。触ってみると、教科書の図の意味がすごく理解しやすくなる。
物理モデルと主要な数式
放射性崩壊の基本は、単位時間あたりの崩壊確率が一定（崩壊定数λ）であることです。これから導かれる、時間t後の残存原子数N(t)を表す式が以下です。
$$N(t) = N_0 \, e^{-\lambda t}$$
ここで、$N_0$は初期原子数、$\lambda$は崩壊定数です。崩壊の速さを表す半減期$T_{1/2}$とは、$\lambda = \dfrac{\ln 2}{T_{1/2}}$の関係があります。半減期が短いほどλは大きく、崩壊が速いことを意味します。
放射能（ベクレルBq）は、1秒間に崩壊する原子核の数です。原子数Nと崩壊定数λから以下のように計算されます。
$$A = \lambda N = \frac{\ln 2}{T_{1/2}}N \; [\text{Bq}]$$
原子数が多いほど、また半減期が短い（λが大きい）ほど、放射能Aは高くなります。シミュレーターでは、この値もリアルタイムで計算・表示されています。
ある核種Aが崩壊して別の放射性核種Bになり、さらにBが崩壊するような連鎖（崩壊系列）を記述するのがバテマン方程式です。子核種Bの時間変化は以下の式で表されます。
$$\dfrac{dN_B}{dt} = \lambda_A N_A - \lambda_B N_B$$
右辺第1項$\lambda_A N_A$は親Aからの生成項、第2項$-\lambda_B N_B$は子B自身の崩壊項です。この方程式を解くことで、親と子の原子数が時間とともにどのように推移するかを追跡できます。
実世界での応用
放射性医薬品・核医学治療：甲状腺がんの治療に使われるI-131（半減期8日）は、その半減期の短さが鍵です。体内に取り込まれた放射性ヨウ素が、治療に必要な期間だけ放射線を出した後、比較的速やかに崩壊して無害化します。投与量の計画には半減期と崩壊の計算が不可欠です。
放射性年代測定：考古学や地質学で使われるC-14法は、生物が死んでから大気とのC-14交換が止まり、その含有量が半減期5730年で減っていく原理を利用します。試料中のC-14と安定同位体C-12の比率を測定し、上記の指数減衰式に当てはめることで年代を推定します。
環境モニタリング・除染：原子力事故後に問題となるCs-137（半減期約30年）は、半減期が長いため環境中に長期に残留します。土壌や食品中の汚染レベル（放射能Bq/kg）を評価し、除染計画や摂取制限の指標を定める際に、半減期を考慮した減衰予測が行われます。
放射線源・非破壊検査：工業用の非破壊検査（透過写真）などに使われるCo-60線源（半減期約5年）は、その強度（放射能）が時間とともに減衰します。検査の精度を保つためには、半減期に基づいて線源の強度を補正したり、交換時期を計画したりする必要があります。
よくある誤解と注意点
まず、「半減期が過ぎれば放射能はゼロになる」というのは大きな誤解だ。半減期は「半分になる時間」だから、1回で1/2、2回で1/4、3回で1/8…と減っていく。例えば、初期原子数100万個のCo-60（半減期5.27年）をシミュレートすると、10半減期（約53年）経っても約1000個は残る。実務では、放射性廃棄物の管理期間を考える際、この「事実上無視できるレベル」に至るまでの長さを認識することが極めて重要になる。
次に、「放射能（Bq）と原子数は比例するが、半減期にも依存する」点だ。ツールで試してみてほしい。初期原子数を100万個に固定し、半減期が8日のI-131と5730年のC-14を比べると、初期の放射能はI-131の方が圧倒的に高い。逆に、同じ放射能1MBq（メガベクレル）を出すために必要な原子数は、半減期の長いC-14の方が莫大になる。線源の取り扱いでは、この「放射能」の値が安全基準の直接的な指標となることを忘れずに。
最後に、バテマン方程式の「子核種の半減期」設定の落とし穴。親よりも子の半減期が極端に長い場合（例えば、親の半減期を1日、子を100年に設定）、子核種がほとんど崩壊せずに溜まり続け、長期にわたる管理対象となる。これは、原子炉で生成される長寿命放射性核種の問題そのものだ。シミュレーションパラメータをいじる時は、現実の核種データを当てはめて考える癖をつけよう。
関連する工学分野
このシミュレーターの根幹をなす「指数関数的減衰」と「連立微分方程式の数値解法」は、CAEを支える広範な工学分野に直接応用されている。まず挙げるのは、「化学反応工学」だ。放射性崩壊の一次反応速度式 $ -dC/dt = kC $ は、酵素反応や薬物の体内動態（ファーマコキネティクス）のモデルと数学的に同一だ。シミュレーターでλ（崩壊定数）を変える感覚は、反応速度定数kを変える感覚にそのまま通じる。
次に、「信頼性工学」や「故障物理学」だ。大量の部品やシステムが、時間とともに一定の確率で故障する現象は、原子の崩壊と酷似している。例えば、半減期を「平均故障間隔（MTBF）」に置き換えれば、時間経過による稼働部品数の減り方を予測する信頼度曲線を描ける。バテマン方程式は、直列システムの故障や、在庫管理のモデル化にも使われる連立微分方程式の典型例なんだ。
さらに、「熱流体解析（CFD）」における物質拡散や化学種輸送の計算でも、生成と消滅を記述するソース項の扱いはバテマン方程式と本質的に同じだ。放射性物質の大気拡散シミュレーションは、これらの方程式を流体場に組み合わせた応用例と言える。一見特殊な「放射性崩壊」のモデルが、実は工学の共通言語であることを感じ取ってほしい。
発展的な学習のために
このツールに慣れたら、次のステップとして「微分方程式の数値解法」そのものを学ぶことを強く勧める。このシミュレーターでは、おそらくオイラー法やルンゲ・クッタ法といった手法で $dN/dt = -λN$ をコンピュータに解かせている。ExcelやPython（NumPy/SciPy）で自分で簡単な数値計算コードを書いてみると、「なぜ時間ステップを細かく取る必要があるのか」といった計算精度の課題が体感できる。これがCAEエンジニアの必須スキルだ。
数学的背景としては、「ラプラス変換」の勉強が役立つ。バテマン方程式のような連立線形常微分方程式は、ラプラス変換を使うと代数的に解く（解析解を求める）ことができる。ツールのグラフが描く曲線の形を、数式の導出から理解できるようになれば、現象の本質をより深く捉えられる。まずは1階と2階の方程式で試してみよう。
次のトピックとしては、「放射線の物質との相互作用（遮蔽計算）」に進むのが現実的だ。崩壊で出たガンマ線が、鉄やコンクリートを通過する時にどのように減衰するかは、ここで学んだ指数減衰則（ $I = I_0 e^{-μx}$ ）で記述される。半減期が「時間」に関する減衰なら、こちらは「距離」に関する減衰だ。シミュレーターで身につけた指数関数への直感が、遮蔽設計の基礎で再び活きてくるはずだ。