1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Poisson distribution & exponential decay

  1. Sep 16, 2013 #1
    1. The problem statement, all variables and given/known data
    t(s) = 1 15 30 45 60 75 90 105 120 135
    N(counts) = 106 80 98 75 74 73 49 38 37 22
    Consider a decaying radioactive source whose activity is measured at intervals of 15 seconds. the total counts during each period are given. What is the lifetime τ of the source? What are the uncertainties on N_0 and τ?

    2. Relevant equations
    N(t) = N_0*exp(-t/τ)


    3. The attempt at a solution
    ln[N(t)]=ln[N_0]+(-t/τ) gives an affine function, so if I plot it I can get N_0 and τ from the graph. However this doesn't feel satisfactory to me, I'd prefer to do it mathematically, but I can't figure it out. I don't know how to do a linear fit for a Poisson distributed set of data. I think the error on the N(t) is given by (σ/√n) where n is the number of data points and σ=√avg[N(t)]. Any help?
     
  2. jcsd
  3. Sep 17, 2013 #2

    mfb

    User Avatar
    2016 Award

    Staff: Mentor

    There are formulas how to get the uncertainties for a linear fit, if you know the uncertainties of the individual values. Alternatively, use a computer, fit programs now those formulas.
    The individual values follow (approximately) a Poisson distribution.

    Strictly speaking, this direct approach will lead to slightly wrong values. This can be seen with datapoints that are zero (not present here, but they are just the most extreme case) - the calculated uncertainty would be zero, which is obviously wrong. There are better methods to treat this properly, but I think that is beyond the scope of this task.
     
  4. Sep 17, 2013 #3

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    I'm a bit confused. There are ten timestamps, starting at (effectively) 0, so that should be nine intervals, but there are ten counts.
     
  5. Sep 17, 2013 #4

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    Glossing over the counts versus intervals problem, I tried attacking it from first principles using MLE. It gets nasty.
    Suppose there are N potential decay events initially, all with parameter lambda. The prob that a given decay occurs in (ti, ti+1) = ##e^{-\lambda t_i}-e^{-\lambda t_{i+1}}##.
    The prob that ki occur in each such interval is ##\left(\stackrel{N}{k_1, k_2, ... k_r}\right)\prod \left(e^{-\lambda t_{i}}-e^{-\lambda t_{i+1}}\right)^{k_i}##
    Taking logs and differentiating wrt lambda:
    ##\Sigma k_i \frac{t_{i+1}e^{-\lambda t_{i+1}}-t_{i}e^{-\lambda t_{i}}}{e^{-\lambda t_{i}}-e^{-\lambda t_{i+1}}}=0##
    Solving that for lambda will not be fun.
     
  6. Sep 17, 2013 #5

    mfb

    User Avatar
    2016 Award

    Staff: Mentor

    Why do you want to solve it for lambda? I would use a numerical approximation.
     
  7. Sep 17, 2013 #6

    TSny

    User Avatar
    Homework Helper
    Gold Member

    Interesting calculation.

    Since in this problem, ##t_{i+1} = t_i + \Delta t## where ##\Delta t## is fixed at 15 s, it doesn't seem too hard to solve your equation for ##\lambda##. [Unless I'm making a stupid mistake.]
     
  8. Sep 18, 2013 #7

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

     
  9. Sep 18, 2013 #8

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    You're right - I looked at that too briefly. If the time intervals are all T, looks like it reduces to ##e^{-\lambda T} = \frac{x}{x+1}## where ##x = \frac{\Sigma i k_i}{\Sigma k_ i}##. Is that what you get?
    Not sure how to go about determining the uncertainty though.
     
  10. Sep 18, 2013 #9

    mfb

    User Avatar
    2016 Award

    Staff: Mentor

    I don't think this has to be done analytically, and it is a standard fit problem. Find lambda where the likelihood is maximal, find lambda+, lambda- where the likelihood gets worse by a factor of [I would have to check that].
     
  11. Sep 18, 2013 #10

    TSny

    User Avatar
    Homework Helper
    Gold Member

    Yes, that's what I get.
     
  12. Sep 18, 2013 #11

    haruspex

    User Avatar
    Science Advisor
    Homework Helper
    Gold Member
    2016 Award

    I don't understand what you mean by a standard fit problem in this context. The decays are counts in an interval, not instantaneous decay rates. Doesn't that complicate things a little? Moreover, the counts in one interval are not independent of those in the next.
     
  13. Sep 18, 2013 #12

    mfb

    User Avatar
    2016 Award

    Staff: Mentor

    The number of decays in an interval is a binomial distribution where N is the number of atoms at the start and p follows from lambda (actually, we can estimate p and find lambda afterwards - p is all we need).

    Good point. Well, it looks easy to integrate that in the likelihood. Just calculate it interval by interval. The total number of atoms is a second free parameter, of course.
    Once you have the likelihood function, use a standard fit algorithm to find the maximum and the uncertainty on the parameters.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: Poisson distribution & exponential decay
  1. Exponential decay (Replies: 2)

  2. Exponential decay (Replies: 3)

Loading...