Monte Carlo Model of Charge Recombination in Feldspars

Feldspars (aluminosilicates of potassium, sodium and calcium) are very common geological materials that emit light when heated or after exposure to visible or infrared radiation. These processes can occur in nature. Two different modes of light emission are thermoluminescence (TL) and optically stimulated luminescence (OSL). When feldspars are irradiated, separation of charge can occur, creating closely spaced electron-hole pairs. These pairs can recombine to emit light by quantum tunneling.
The study of how these electron-hole pairs recombine over geological times can be applied to luminescence dating and luminescence dosimetry.
This Demonstration uses Monte Carlo methods to simulate the recombination process, taking account of the physical parameters that characterize feldspar crystals. We model 10 electron-hole pairs located at random positions in a cube. The plot shows how these pairs recombine in the natural crystal over millions of years. The tunneling process is very slow and takes place over geological timescales. The horizontal error bars are the standard deviation of Monte Carlo computations, where . The horizontal axis represents the elapsed geological time.
  • Contributed by: Christopher W. Kulp and Vasilis Pagonis
  • (Department of Astronomy & Physics, Lycoming College, Williamsport, PA, USA and Department of Physics, McDaniel College, Westminster, MD, USA)



  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Use the buttons to set values for the potential barrier constant , the tunneling frequency , the size of the cubic cell and the number of cells.
As the tunneling frequency increases or as the length of the side of the cube decreases, the charge recombination in the crystal increases. The process depends very strongly on the potential barrier penetration constant . The solid line in the graph represents an analytical solution [1]. For more details of the simulation, see [2].
Function Details
Pairing produces a list of times (in seconds) that it takes for each electron-hole pair to recombine. The arguments for Pairing include the list of initial locations of the electrons (epos0) and holes (hpos0), the tunneling frequency in Hz (stun) and the potential barrier penetration constant in (α).
ChargeLossSimulation performs the Monte Carlo simulation of the recombination of charge due to ground-state tunneling. The output of ChargeLossSimulation is the number of charge pairs in the cell as a function of time. In addition, ChargeLossSimulation computes the analytical equation for the number of pairs remaining as a function of time. The arguments for ChargeLossSimulation include the number of cells used in the simulation (nsolids), the tunneling frequency in Hz (stun), the potential barrier penetration constant in (α) and the length of the cell side in nanometers (sidelengthnm).
[1] G. Kitis and V. Pagonis, "Analytical Solutions for Stimulated Luminescence Emission from Tunneling Recombination in Random Distributions of Defects," Journal of Luminescence, 137, 2013 pp. 109–115. doi:10.1016/j.jlumin.2012.12.042.
[2] V. Pagonis and C. Kulp, "Monte Carlo Simulations of Tunneling Phenomena and Nearest Neighbor Hopping Mechanism in Feldspars," Journal of Luminesence, 181, 2017 pp. 114–120. doi:10.1016/j.jlumin.2016.09.014.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+