Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
On the sampling of step length in Monte Carlo simulation of trajectories with very small mean free path
Stockholm University, Faculty of Science, Department of Physics.
2012 (English)In: Radiation Physics and Chemistry, ISSN 0969-806X, E-ISSN 1879-0895, Vol. 81, no 11, p. 1703-1709Article in journal (Refereed) Published
Abstract [en]

In an event-by-event simulation of the trajectory of a particle moving in matter it is usually assumed that the probability for the particle to travel a distance s without interaction is exp(-s/lambda), where lambda = (n . sigma)(-1) is the total mean free path, n the number of scatterers per unit volume and sigma the total cross section per scatterer. The step length s between scattering events is then generated by means of a sampling formula s = -lambda ln(1-R), where R a random number in the interval 0 < R < 1. It is here argued that this conventional sampling method, which basically assumes that the scattering medium may be regarded as a homogeneous continuum, may be erroneous unless lambda is much larger than the average distance d(nn) between nearest neighbour scatterers, estimated by d(nn) = n(-1/3). An alternative sampling method (M sampling) is proposed with a fixed step length D = d(nn) and a finite probability I = 1-exp(-D/lambda) of a single elastic or inelastic scattering event at the end of each step. According to this method, conventional sampling may exaggerate the number of events per unit path length; the corrected mean free path between events is found to be lambda(c) = D/(1-exp(-D/lambda)). The correction is substantial when lambda is comparable to or smaller than D, in practice for very low energy particles in liquids and solids. Consequently, quantities like stopping power may then be overestimated, while transport mean free path may be underestimated. In the opposite limit lambda >> D, conventional and M sampling produce the same result. Present results further indicate that conventional sampling using the corrected total mean free path lambda(c) is a good approximation to M sampling.

Place, publisher, year, edition, pages
2012. Vol. 81, no 11, p. 1703-1709
Keywords [en]
Track simulation, Mean free path, Step length
National Category
Physical Chemistry
Identifiers
URN: urn:nbn:se:su:diva-81280DOI: 10.1016/j.radphyschem.2012.06.058ISI: 000308118100006Scopus ID: 2-s2.0-84864629022OAI: oai:DiVA.org:su-81280DiVA, id: diva2:560515
Note

AuthorCount:1;

Available from: 2012-10-15 Created: 2012-10-15 Last updated: 2022-10-03Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

Liljequist, David

Search in DiVA

By author/editor
Liljequist, David
By organisation
Department of Physics
In the same journal
Radiation Physics and Chemistry
Physical Chemistry

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 42 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf