Skip to main content
Log in

Stochastic Simulation Algorithms for Solving Transient Anisotropic Diffusion-recombination Equations and Application to Cathodoluminescence Imaging

  • Published:
Methodology and Computing in Applied Probability Aims and scope Submit manuscript

Abstract

A meshless Random Walk on arbitrary parallelepipeds simulation algorithm is developed and implemented for solving transient anisotropic diffusion-reaction equations. In contrast to the conventional Feynman-Kac based algorithm the suggested method does not use small time step simulations of the relevant diffusion processes. Instead, exact simulation of large random jumps over a set of appropriately constructed parallelepipeds in the domain is carried out. This decreases the cost of simulations considerably especially for domains with complicated boundary shape. Application to the problem of time-resolved cathodoluminescence intensity calculations for semiconductor materials with a set of threading dislocations is given. Important issues are the construction of an efficient sampling method from the first passage time density and the position distribution on the surface of an arbitrary parallelepiped. We combine a rejection algorithm and a probability density tabulation approach to construct optimal sampling methods from different densities including the random time a particle spends in a parallelepiped before it is absorbed inside it. We present in the last section results of computer simulation for the evaluation of the exciton flux to dislocations and a plane substrate, the cathodoluminescence intensity for threading dislocations imaging, and the concentration of the survived excitons. In addition, to validate the developed algorithms we have compared the computer simulations with the exact results, and obtained a perfect agreement.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9

Similar content being viewed by others

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

  • Courant R, Friedrichs K, Lewy H (1928) Überber die partiellen Differenzengleichungen der mathematischen Physik. Math Annalen 100:32–74

    Article  MathSciNet  MATH  Google Scholar 

  • Deaconu M, Lejay A (2006) A random walk on rectangles algorithm. Methodol Comput Appl Probab 8(1):135

    Article  MathSciNet  MATH  Google Scholar 

  • Devroye L (1981) The series method for random variate generation and its application to the Kolmogorov-Smirnov distribution. Am J Math Manag Sci 1(4):359–379

    Article  MathSciNet  MATH  Google Scholar 

  • Dynkin EB (2006) Theory of Markov Processes. Dover, Mineola, NY, USA

    MATH  Google Scholar 

  • Elepov BS, Kronberg AA, Mikhailov GA, Sabelfeld KK (1980) Solution of boundary value problems by the Monte Carlo methods. Nauka, Novosibirsk

    Google Scholar 

  • Ermakov SM, Nekrutkin VV, Sipin AS (1989) Random Processes for Classical Equations of Mathematical Physics. Kluwer Academic Publishers, Dordrecht, Springer Science & Business Media

    Book  MATH  Google Scholar 

  • Ermakov SM, Sipin AS (2009) The random walk on semi-spheres and applications for solving boundary value problems, Vest. St.Petersbourg University, ser. 1: N3, 9–18 (in russian)

  • Haji-Sheikh A, Sparrow EM (1966) The floating random walk and its application to Monte Carlo solutions of heat equations. SIAM J Appl Math 14(2):570–589

    Article  MathSciNet  MATH  Google Scholar 

  • Ito K, Mckean P (1965) Diffusion processes and their sample paths. Springer, Berlin

    MATH  Google Scholar 

  • Kac M (1951) On some connections between probability theory and differential and integral equations. 2-d Berk. Symp. 189–215, University of California Press

  • Kloeden P, Platen E, Schurz H (2012) Numerical Solution of Stochastic Differential Equations. Springer, Heidelberg-Berlin

    MATH  Google Scholar 

  • Liu W, Carlin JF, Grandjean N, Deveaud B, Jacopin G (2016) Exciton dynamics at a single dislocation in GaN probed by picosecond time-resolved cathodoluminescence. Appl Phys Lett 109(4):042101

    Article  Google Scholar 

  • Milewski S (2018) Combination of the meshless finite difference approach with the Monte Carlo random walk technique for solution of elliptic problems. Comp Math Appl 76(4): 854–876. https://doi.org/10.1016/j.camwa.2018.05.025

  • Muller ME (1956) Some continuous Monte Carlo methods for the Dirichlet problem. Ann Math Statist 27(3):569–589

    Article  MathSciNet  MATH  Google Scholar 

  • Polyanin AD (2002) Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC, 2002, N.Y. 800 pp

  • Sabelfeld KK (1991) Monte Carlo Methods in Boundary Value Problems. Springer, Berlin

    Google Scholar 

  • Sabelfeld KK (2016) Random walk on spheres method for solving drift-diffusion problems. Monte Carlo Methods Appl 22:N(4):265–281

  • Sabelfeld KK (2017) Random walk on spheres method for solving transient drift-diffusion-reaction problems, Monte Carlo Methods Appl 23:N(3):189–212

  • Sabelfeld KK (2017) A mesh free floating random walk method for solving diffusion imaging problems. Statist. Probab. Lett. 121:6–11

    Article  MathSciNet  MATH  Google Scholar 

  • Sabelfeld KK, Kaganer VM, Pfüller C, Brandt O (2017) Dislocation contrast in cathodoluminescence and electron-beam induced current maps on GaN0001. J Phys D: Appl Phys 50(N40):405101 (11pp)

  • Sabelfeld K, Kireeva A (2017) Probability distribution of the life time of a drift-diffusion-reaction process inside a sphere with applications to transient cathodoluminescence imaging. Monte Carlo Methods Appl 24(2):79–92

    Article  MathSciNet  MATH  Google Scholar 

  • Sabelfeld KK, Simonov NA (2016) Stochastic Methods for Boundary Value Problems. Numerics for High-dimensional PDEs and Applications, De Gruyter, Berlin

    Book  MATH  Google Scholar 

Download references

Acknowledgements

Support of the Russian Science Foundation under Grant 19-11-00019 is gratefully acknowledged.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Karl K. Sabelfeld.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sabelfeld, K.K., Kireeva, A.E. Stochastic Simulation Algorithms for Solving Transient Anisotropic Diffusion-recombination Equations and Application to Cathodoluminescence Imaging. Methodol Comput Appl Probab 24, 3029–3048 (2022). https://doi.org/10.1007/s11009-022-09968-9

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11009-022-09968-9

Keywords

PACS

Mathematics Subject Classification (2000)

Navigation