過去ログ

uchinyan
(若干,追加修正しました。)
次に,#33629,半減期の問題,に関してです。
この問題を解くには,若干の予備知識が必要です。
なお,今回も現れる対数は自然対数です。

まず,ある物質が崩壊する物質であったとし,各時点で N 個あったとすると,
その減少する様子は N に比例することが知られています。
つまり,λ > 0 を比例定数,崩壊定数,壊変定数などといいます,t を時間とすると,
dN/dt = - λN
です。t = 0 で N = N0 としてこの微分方程式を解くと,
1/N * dN/dt = - λ
∫(1/N * dN/dt)dt = ∫(-λ)dt
∫(1/N)dN = ∫(-λ)dt
log|N| = - λ * t + 積分定数
N = 定数 * e^(-λt)
N = N0 * e^(-λt)
です。ここで,t -> t + T となったとき N -> N/2 となる T を半減期といいます。
N = N0 * e^(-λt)
N/2 = N0 * e^(-λ(t + T))
なので,
2 = e^(λT)
λ = log2/T,T = log2/λ
になります。そこで,
N = N0 * e^(- log2/T * t) = N0 * e^(- log2 * t/T) = N0 * (e^log(1/2))^(t/T)
log の定義から e^(log2) = 2 なので,
N = N0 * (1/2)^(t/T)
とも書けます。

さて,以上の準備の下に解いていきましょう。
物質 A,B,C の t 年後の重さを A,B,C とします。単位は kg です。
また,A,B の崩壊定数をλa,λb とし,A,B の半減期を Ta,Tb とします。
ここで,注意というか仮定ですが,先ほどの予備知識では個数でしたが,今回は重さです。
これは明らかに違うものなので,厳密には,その間をつなぐもの,原子 1 個あたりの重さなど,
が必要でしょうが,提示された問題では特に言及がないので,
物質の崩壊後もこれは変わらない,A,B,C の原子 1 個あたりの重さはほとんど同じ,と仮定しました。
すると,全体に定数が掛かるだけで,先ほどの崩壊の話が重さにもそのまま成立します。
以下では,この仮定の下に解きます。

まず,A ですが,これは明らかに,
dA/dt = - λa * A
がいえます。これを t = 0,A = 1 で解くと,
A = e^(- λa * t) = (1/2)^(t/Ta)
です。

次に,B です。B は,A からの崩壊分 λa * A が増え,C への崩壊分 λb * B が減るので,
dB/dt = λa * A - λb * B
dB/dt = λa * e^(- λa * t) - λb * B
になります。これを,t = 0,B = 0 で解けばいいです。
dB/dt + λb * B = λa * e^(- λa * t)
e^(λb * t) * dB/dt + λb * e^(λb * t) * B = λa * e^((λb - λa) * t)
d(e^(λb * t) * B)/dt = λa * e^((λb - λa) * t)
e^(λb * t) * B = λa/(λb - λa) * e^((λb - λa) * t) + 積分定数
t = 0,B = 0 より,
e^(λb * t) * B = λa/(λb - λa) * e^((λb - λa) * t) - λa/(λb - λa)
B = λa/(λb - λa) * (e^(- λa * t) - e^(- λb * t))
B = (1/Ta)/(1/Tb - 1/Ta) * ((1/2)^(t/Ta) - (1/2)^(t/Tb))
B = Tb/(Tb - Ta) * ((1/2)^(t/Tb) - (1/2)^(t/Ta))

最後に,C です。C は B からの崩壊分 λb * B が増えるだけなので,
dC/dt = λb * B
dC/dt = λaλb/(λb - λa) * (e^(- λa * t) - e^(- λb * t))
になります。これを,t = 0,C = 0 で解けばいいです。
C = λaλb/(λb - λa) * (- 1/λa * (e^(- λa * t) - 1) + 1/λb * (e^(- λb * t) - 1))
C = 1 - λb/(λb - λa) * e^(- λa * t) + λa/(λb - λa) * e^(- λb * t)
C = 1 - (1/Tb)/(1/Tb - 1/Ta) * (1/2)^(t/Ta) + (1/Ta)/(1/Tb - 1/Ta) * (1/2)^(t/Tb)
C = 1 - Tb/(Tb - Ta) * (1/2)^(t/Tb) + Ta/(Tb - Ta) * (1/2)^(t/Ta)

以上で,A,B,C の一般形が求まりました。

さて,Ta = 100,Tb = 300 なので,
A = (1/2)^(t/100)
B = 3/2 * (1/2)^(t/300) - 3/2 * (1/2)^(t/100)
C = 1 - 3/2 * (1/2)^(t/300) + 1/2 * (1/2)^(t/100)
t = 500 を求めると,
A = (1/2)^5
B = 3/2 * (1/2)^(5/3) - 3/2 * (1/2)^5 = 3 * (1/2)^(8/3) - 3 * (1/2)^6
C = 1 - 3/2 * (1/2)^(5/3) + 1/2 * (1/2)^5 = 1 - 3 * (1/2)^(8/3) + (1/2)^6
となります。

なお,式の立て方から明らかですが,質量欠損は考慮しなかったので,
A + B + C = 1
が成立しています。

これぐらいも,予備知識つきで,大学入試の数学に出てもおかしくはないかなぁ...
それとも,物理? 微分方程式が絡むので,物理では,ちょっと無理かな。
ネコの住む家   1月14日(水) 18:34:39   MAIL:uchi@sco.bekkoame.ne.jp   33639

過去問題一覧へ戻る