Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Numerical integration of Bethe formula in Excel: help!

  1. Feb 21, 2013 #1
    Hello,

    I am a PhD student in geology in need of help from a physicist! Can somebody spot the mistake in my spreadsheet?

    I am using an Excel spreadsheet (attached .xls) to integrate the Bethe formula in order to estimate the range (in cm) of alpha particles in quartz. I am using the version of the formula given in Groom and Klein's "Passage of Particles through Matter": http://pdg.lbl.gov/2009/reviews/rpp2009-rev-passage-particles-matter.pdf (but ignoring the density effect correction because I am not interested in high energies).

    In Excel, I can plot the Bragg peak (dE/dX against x) or the energy of the alpha particle (E against x). Something is very wrong: the Bragg curve just goes up and up and never comes down again (don't be misled by Excel plotting #NUM! as 0). I would expect dE/dX to start out at something like 500 and then smoothly increase up to the peak (like the Bragg peaks in Google image search).

    I must be doing something right, though, because it looks as if the stopping distance is going to be close to what I would expect (about 0.04 cm).

    Can anybody find the problem in the attached spreadsheet? Reward: credit in my thesis acknowledgements!

    Many thanks!

    Moho
     

    Attached Files:

  2. jcsd
  3. Feb 21, 2013 #2

    mfb

    User Avatar
    2016 Award

    Staff: Mentor

    For low energy (βγ<0.03), the Bethe-Bloch formula gives wrong results, as shown in figure 27.1 in the pdf.
    Apart from that: After the energy is 0, the energy loss should be 0 as well, giving a nice Bragg peak.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Numerical integration of Bethe formula in Excel: help!
  1. Bethe formula question (Replies: 3)

  2. Bethe Bloch formula (Replies: 1)

Loading...