Radioactive Decay Differential Equation Question

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
31 replies · 6K views
TFM
Messages
1,016
Reaction score
0

Homework Statement



Unobtainium is a radioactive (and fictional) element, with its rate of decay being proportional to the amount of Unobtainium, [tex]u' = −k u[/tex]. The half-life of Unobtainium, in other words the time it takes for half the initial amount to decay, is [tex]2.3 × 10^9[/tex]yr.

(a) Find the value of k.

(b) How many decay events will there be per second in a 1 kg block of Unobtainium, given
that it has an atomic mass of 415?

Homework Equations



N/A

The Attempt at a Solution



I think I am doing the right thing, but I can't seem to work out k. I have a feeling I need limits to me integrals. so far I have:

[tex]\frac{du}{dt} = -ku[/tex]

[tex]\frac{1}{u} dt = -k dt[/tex]

[tex]ln u = -kt + C[/tex]

Where C is a constant of integration

taking exponentials of both sides gives:

[tex]u = e^{-kt} + C[/tex]

If we had a [tex]u_{0}[/tex] somewhere, we know the half life, and thus could say that [tex]u = \frac{1}{2} u_o[/tex], stick some values for u and u0 in, and thus could probably work out k.

Any ideas where I am going wrong?

TFM
 
Physics news on Phys.org
That was rather silly of me, thanks.

I still need a [tex]u_0[/tex] in there somewhere though to calculate k, don't I...?

I also need to work out the constant C as well.

TFM
 
TFM said:
That was rather silly of me, thanks.

I still need a [tex]u_0[/tex] in there somewhere though to calculate k, don't I...?

I also need to work out the constant C as well.

Are you tired? :redface:

How about u0 = C?
 
How did you get C to equal [tex]u_0[/tex]

?

TFM
 
I see, so:

[tex]u = Ce^{-k*0}[/tex]

[tex]u = Ce^0[/tex]

and:

[tex]e^0 = 1[/tex]

Thus:

[tex]u = C[/tex] at t = 0

thus C = [tex]u_0[/tex]

Thanks,

So now:

[tex]u = u_0e^{-kt}[/tex]

Insert values for u and u_0 that follow half-life, eg u_0 = 10, u = 5:

[tex]5 = 10e^{-k*2.3 * 10^9}[/tex]

Divide both sides by 10:

[tex]\frac{5}{10} = e^{-k*2.3 * 10^9}[/tex]

Take logs of both sides:

[tex]ln(\frac{1}{2}) = -k*2.3 * 10^9[/tex]

divide by half life:

[tex]\frac{ln\frac{1}{2}}{2.3 * 10^9} = -k[/tex]

thus

[tex]k = -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

Does this look right?

TFM
 
TFM said:
I see, so:

[tex]u = Ce^{-k*0}[/tex]

[tex]u = Ce^0[/tex]

and:

[tex]e^0 = 1[/tex]

Thus:

[tex]u = C[/tex] at t = 0

thus C = [tex]u_0[/tex]

Thanks,

So now:

[tex]u = u_0e^{-kt}[/tex]

Insert values for u and u_0 that follow half-life, eg u_0 = 10, u = 5:

[tex]5 = 10e^{-k*2.3 * 10^9}[/tex]

Divide both sides by 10:

[tex]\frac{5}{10} = e^{-k*2.3 * 10^9}[/tex]

Take logs of both sides:

[tex]ln(\frac{1}{2}) = -k*2.3 * 10^9[/tex]

divide by half life:

[tex]\frac{ln\frac{1}{2}}{2.3 * 10^9} = -k[/tex]

thus

[tex]k = -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

Does this look right?

TFM

Yes (but a bit long-winded :wink:).

Now simplify -ln(1/2). :smile:
 
So:

[tex]k = -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

[tex]-ln\frac{1}{2}[/tex]

Thi simplifies to:

[tex]ln(1) - ln(2)[/tex]

ln(1) = 0

thus:

[tex]k = -\frac{-ln2}{2.3 * 10^9} = \frac{ln2}{2.3 * 10^9}[/tex]

TFM
 
Excellent. Thanks.

So for part b)

How many decay events will there be per second in a 1 kg block of Unobtainium, given
that it has an atomic mass of 415?


We know now that:

[tex]u = u_0e^{\frac{ln\frac{1}{2}}{2.3 * 10^9}t}[/tex]

Do we insert t = 1 second, and u_0 = 1?

TFM
 
The question states that it has an atomic mass of 415. This means, if I remember correctly, it has 415 atoms per kilogram, thus 1 kg = 415...?

TFM
 
TFM said:
The question states that it has an atomic mass of 415. This means, if I remember correctly, it has 415 atoms per kilogram, thus 1 kg = 415...?

TFM

oooh, that's chemistry … I've no idea! :redface:

though that does see a rather heavy atom …

(and, come to think of it :rolleyes:, doesn't hydrogen have an atomic mass of 1?)
 
Well on wikipedia:

The atomic mass of an 56Fe isotope is 55.935 u and one mole of 56Fe will in theory weigh 55.935g,

and gives an equation:

[tex]m_{\rm{u}}={m_{\rm{grams}} \over N_{A}}[/tex]

Is this useful?

TFM
 
Hey,

It should be [tex]415 g*mol^{-1}[/tex]. To find the number of mols divide the total mass by the molar mass. To find the number of atoms multiply the number of mols by Avogadro's number. Then put the number of atoms as your value of u to find the number of atoms decaying per second.
 
Thanks, so its

[tex]415 grams m^{-1}[/tex]

we have 1000 grams

so:

[tex]\frac{1000}{415} = 2.4 mols[/tex]

So now we times this number by Avagaro's Constant, [tex]6.022 * 10^{23}[/tex]

This gives the number of particles to be:

[tex]1.45 * 10^{24}[/tex]

So the equation is:

[tex]u = u_0e^{\frac{ln\frac{1}{2}}{2.3 * 10^9}t}[/tex]

t = 1
[tex]u_0 = 1.45 * 10^{24}[/tex]

[tex]u = (1.45 * 10^{24})e^{\frac{ln\frac{1}{2}}{2.3 * 10^9}1}[/tex]

This gives the answer to be

U = 1.45e^24

:mad:

But that's because the answer is so small.

we need to find the number of decay events, so I assume we take away u from u_o,

so this gives:

[tex]4.27*10^{14}[/tex] events per second. Seem right?

TFM
 
Your amount for [tex]u_{0}[/tex] is correct.

The rest is incorrect. [tex]u \&\ u_{0}[/tex] are simply the number of atoms in your sample and as such have no units. You want the number of atoms decaying per second in a 1 kg sample. You have the total number of atoms in 1kg, [tex]u_{0}[/tex]. You have already found the decay constant, in this case k. You need only to multiply these two values to find your answer. Can you see why? (Consider the units of k.)

Also don't forget to change k from per year to per second.
 
Well the decay constant is 1/mean lifetime

So the k at the moment is:

[tex]k = -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

Converting into seconds gives a value of [tex]-9.4923*10^{-3}[/tex]

So multiplying u0 with k(secs) gives: [tex]1.376*10^22[/tex]

?

TFM
 
To convert k to per seconds you should times your k value by 1 year and divide by the number of seconds in one year, 60*60*24*365=3.1536*10^7.
 
this is probably going to sound rather silly but...How do you times it by 1 year? Wouldn't it just be times by 1?

?

TFM :redface:
 
K is the decay constant in this question, so the mean liftime is 1/k.
 
Yeah you just times it by one. I find it helps if you always consider the units you are using, so that per year times year removes the unit per year from your equation.
 
Okay, that does actually make a lot of sense.

So:

[tex]k = -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

[tex]1 * -\frac{ln\frac{1}{2}}{2.3 * 10^9}[/tex]

[tex]\frac{-\frac{ln\frac{1}{2}}{2.3 * 10^9}}{3.1536*10^7}[/tex]

?

TFM
 
Looks good to me. Then just times that by the number of atoms and that should give you the number of atoms decaying per second.
 
So:

[tex]\frac{-\frac{ln\frac{1}{2}}{2.3 * 10^9}}{3.1536*10^7}[/tex]

[tex]\frac{-\frac{ln\frac{1}{2}}{2.3 * 10^9}}{3.1536*10^7}* 1.45 * 10^{24}[/tex]

I'm sure that this can be canceled down?

TFM
 
Err, just plug the numbers into a calculator to get the answer. That figure you have there is somewhat intimidating, :eek:
 
Now why didn't I think of that...?

I think its because usually they say to keep the logs as logs

So I get a value of [tex]1.39 * 10^7[/tex]

Look about right?

TFM
 
I suppose you could collect termsto get,

[tex]\frac{\ln(2)}{5.002*10^{8}}[/tex]