Error Propogation for Half-Life

Join the discussion
Registration is free. Start your own thread to ask a follow-up.
4 replies · 3K views · last reply
Alejandro Golob
Messages
2
Reaction score
0
In this homework question we are told to calculate the half-life of an isotope based on count-rates before and after a given time interval, the relevant equation is given below.

Half-life = t1/2 = tln(2)/ln (A/A0)

The second part asks to determine the standard deviation of the half-life due to counting statistics. We have an uncertainty for the value A and the value A0.

I have found the following solution but cannot quite follow it, so was hoping someone might be able to explain/walk through it with me.

d/dt (t1/2) = − (ln (2) t/(A/A0))/ln((A/A0)2)

Plugging in t = 24 and R = 2.875 yields 5.248 hr−1

σ2t1/2 = (5.248)2 ((9.14/(41.4)2 + 16.83/(118.3)2) *(2.875)2

Thus, we find that expected standard deviation of half-life is σ(t1/2) = 1.22 h

Thanks in advance for any help.
 
Physics news on Phys.org
d/dt does not make sense at that point. If you calculate the time-derivative of the expression for the half-life, you get a different result.
Alejandro Golob said:
Plugging in t = 24 and R = 2.875
24 what, and what is R?
The time-derivative of a time should be dimensionless.

You can calculate the uncertainty on the ratio first, and then propagate this to the half-life measurement. ##-\frac{\ln( 2) t }{A/A_0 \ln(A/A_0)^2}## can be useful there as some part of a different formula.

Can you give the numbers you use for your calculations?
 
Thank you for the reply.

t=24hr and R=A/A0=(41.4minute-1)/118.3minute-1

I certainly agree with what you are saying about the time derivative of a time, I wasn't sure on this myself. The half-life itself is not a function of time anyways, so not sure how the time derivative makes sense or why it is invoked. It is possible that this is an error. I found this solution here and have been trying to see how it was arrived at. See the image of the solution below. The equation you have arrived at certainly is consistent with the first part of this solution and makes sense to me, however I am still unclear on how that carries over to the second part.
Handwritten-style typeset solution: S2 = 914/10 - 50 = 41.4 min-1, sigma S1 = sqrt(1683)/10 = 16.83, sigma S2 = sqrt(914)/10 = 9.14, R = S2/S1 = 2.875, half-life = t ln2 / ln R, t1/2 = 15.84 h; part (b) derivative of t1/2 with respect to t and the variance expression yielding sigma t1/2 = 1.22 h.

Many thanks for your help.

Best Regards,
 
Attachments
  • Printed lab solution page showing half-life error propagation: t_1/2 = 15.84 h, and a derivative expression for the uncertainty with result sigma = 1.22 h, followed by Problem 23 on Poisson event intervals with result 0.632.
    upload_2015-12-1_15-25-11.png
    14.5 KB · Views: 535
  • Printed lab solution page showing half-life calculation steps with uncertainty propagation, including the derivative of half-life with respect to R and the resulting standard deviation of 1.22 h.
    upload_2015-12-1_15-25-11.png
    14.5 KB · Views: 554
  • Handwritten-style worked solution computing half-life t = t ln2 / ln R = 15.84 h, then its derivative and standard deviation of 1.22 h via error propagation with S2 = 41.4 min-1, R = 2.875.
    Untitled.png
    12.4 KB · Views: 957
Alejandro Golob said:
I certainly agree with what you are saying about the time derivative of a time, I wasn't sure on this myself. The half-life itself is not a function of time anyways, so not sure how the time derivative makes sense or why it is invoked. It is possible that this is an error.
If you look at the right hand side of the d/dt(t1/2) equation you can see that the derivative taken was with respect to R, not t.
 
There are more issues. The ##\sigma S_i## calculations mix standard deviations (left and middle) and variance (right).
##\sigma R## should have been calculated in a clearer way, and so on. It is the left two brackets in the last formula.
5.248 seems to be a bad approximation of the expression above, I get 5.188.