Skip to main content
SpringerLink
Log in

Adaptive Monte Carlo algorithm for Wigner kernel evaluation

  • Original Article
  • Published:
Neural Computing and Applications Aims and scope Submit manuscript

Abstract

In this paper, we study numerically various approaches, namely an adaptive Monte Carlo algorithm, a particular rank-1 lattice algorithm based on generalized Fibonacci numbers and a Monte Carlo algorithm based on Latin hypercube sampling for computing multidimensional integrals. We compare the performance of the algorithms over three case studies—multidimensional integrals from Bayesian statistics, the so-called Genz test functions and the Wigner kernel—an important issue in quantum mechanics represented by multidimensional integrals. A comprehensive study and an analysis of the computational complexity of the algorithms under consideration has been presented. Adaptive strategy is well-established as an efficient and reliable tool for multidimensional integration of integrands functions with computational peculiarities like peaks. The presented adaptive Monte Carlo algorithm gives reliable results in computing the Wigner kernel by a stochastic approach that has significantly lower computational complexity than the existing deterministic approaches.

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

Access this article

Log in via an institution

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

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2

Similar content being viewed by others

References

  1. Antonov I, Saleev V (1979) An economic method of computing \(LP_{\tau }\)-sequences. USSR Comput Math Phys 19:252–256

    MATH  Google Scholar 

  2. Baraniuk RG, Jones DL (1993) A signal-dependent time-frequency representation: optimal kernel design. IEEE Trans Signal Process 41(4):1589–1602

    Article  MATH  Google Scholar 

  3. Berntsen J, Espelid TO, Genz A (1991) An adaptive algorithm for the approximate calculation of multiple integrals, ACM Trans. Math Softw 17:437–451

    Article  MATH  Google Scholar 

  4. Bratley P, Fox B (1988) Algorithm 659: implementing Sobol’s Quasirandom sequence generator. ACM Trans Math Softw 14(1):88–100

    Article  MathSciNet  MATH  Google Scholar 

  5. Briol F-X, Oates CJ, Girolami M, Osborne MA, Sejdinovic D (2019) Probabilistic integration: a role in statistical computation? Stat Sci 34(1):1–22

    Article  MathSciNet  MATH  Google Scholar 

  6. Cull P, Holloway JL (1989) Computing Fibonacci numbers quickly. Inf Process Lett 32(3):143–149

    Article  MathSciNet  Google Scholar 

  7. Davis PJ, Rabinowitz P (1984) Methods of numerical integration, 2nd edn. Academic Press, London

    MATH  Google Scholar 

  8. Dimov I (2008) Monte Carlo methods for applied scientists. World Scientific, New Jersey, p 291

    MATH  Google Scholar 

  9. Dimov I, Karaivanova A, Georgieva R, Ivanovska S, (2003) Parallel Importance Separation and Adaptive Monte Carlo Algorithms for Multiple Integrals, 5th Int. conf. on NMA, (August 2002) Borovets, Bulgaria, Springer Lecture Notes in Computer Science, 2542. Springer-Verlag, Berlin, Heidelberg, New York, pp 99–107

  10. Dimov I, Georgieva R (2010) Monte Carlo algorithms for evaluating Sobol’ sensitivity indices. Math Comput Simul 81(3):506–514

    Article  MathSciNet  MATH  Google Scholar 

  11. Eglajs V, Audze P (1977) New approach to the design of multifactor experiments: problems of dynamics and strengths. 35 (in Russian). Riga: Zinatne Publishing House pp 104–107

  12. Ermakov SM (1985) Monte Carlo methods and mixed problems. Nauka, Moscow

    Google Scholar 

  13. Feynman RP (1948) Space-time approach to non-relativistic quantum mechanics. Rev Mod Phys: 20

  14. Fox B (1986) Algorithm 647: implementation and relative efficiency of Quasirandom sequence generators. ACM Trans Math Softw 12(4):362–376

    Article  MATH  Google Scholar 

  15. Genz A (1984) Testing multidimensional integration routines. Tools, Methods and Languages for Scientific and Engineering Computation, pp 81–94

  16. Hua LK, Wang Y (1981) Applications of number theory to numerical analysis. Springer, Berlin

    MATH  Google Scholar 

  17. Jarosz W (2008) Efficient Monte Carlo methods for light transport in scattering media, PhD dissertation, UCSD

  18. Joe S, Kuo F (2003) Remark on algorithm 659: implementing Sobol’s Quasirandom sequence generator. ACM Trans Math Softw 29(1):49–57

    Article  MathSciNet  MATH  Google Scholar 

  19. Karaivanova A, Dimov I, Ivanovska S (2001) A quasi-monte carlo method for integration with improved convergence. LNCS 2179:158–165

    MATH  Google Scholar 

  20. Larkin FM (1972) Gaussian measure in Hilbert space and applications in numerical analysis. Rocky Mt J Math 2(3):379–421

    Article  MathSciNet  MATH  Google Scholar 

  21. Larkin FM (1974) Probabilistic error estimates in spline interpolation and quadrature. In: IFIP Congress, pp 605–609

  22. Li Wei, Lingyun Lu, Xie Xiaotian, Yang Ming (2017) A novel extension algorithm for optimized Latin hypercube sampling. J Stat Comput Simul 87(13):2549–2559. https://doi.org/10.1080/00949655.2017.1340475

    Article  MathSciNet  MATH  Google Scholar 

  23. Lin S (2011) Algebraic methods for evaluating integrals in Bayesian statistics, Ph.D. dissertation, UC Berkeley

  24. Lin S, Sturmfels B, Xu Z (2009) Marginal likelihood integrals for mixtures of independence models. J Mach Learn Res 10:1611–1631

    MathSciNet  MATH  Google Scholar 

  25. McKay MD, Beckman RJ, Conover WJ (1979) A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 21(2):239–245. https://doi.org/10.2307/1268522

    Article  MathSciNet  MATH  Google Scholar 

  26. Metropolis N, Ulam S (1949) The Monte Carlo method. J Am Stat Assoc 44(247):335–341

    Article  MATH  Google Scholar 

  27. Minasny BB, McBratney AB (2006) A conditioned Latin hypercube method for sampling in the presence of ancillary information. J Comput Geosci Arch 32(9):1378–1388

    Article  Google Scholar 

  28. Minasny B, McBratney AB (2010) Conditioned Latin hypercube sampling for calibrating soil sensor data to soil properties, Chapter: Proximal Soil Sensing, Progress in Soil Science pp 111–119

  29. O’Hagan A (1991) Bayes–Hermite quadrature. J Stat Plan Inference 29(3):245–260

    Article  MathSciNet  MATH  Google Scholar 

  30. Schurer R (2001) Parallel high-dimensional integration: Quasi-monte carlo versus adaptive cubature rules, Computational Science, ICCS 2001. Springer, Berlin

    MATH  Google Scholar 

  31. Sellier JM (2015) A signed particle formulation of non-relativistic quantum mechanics. J Comput Phys 297:254–265

    Article  MathSciNet  MATH  Google Scholar 

  32. Sellier JM, Dimov I (2016) On a full Monte Carlo approach to quantum mechanics. Phys A Stat Mech Appl 463:45–62

    Article  MathSciNet  MATH  Google Scholar 

  33. Sellier JM, Dimov I (2014) The many-body Wigner Monte Carlo method for time-dependent ab-initio quantum simulations. J Comput Phys 273:589–597

    Article  MathSciNet  MATH  Google Scholar 

  34. Sellier JM, Nedjalkov M, Dimov I (2015) An introduction to applied quantum mechanics in the Wigner Monte Carlo formalism. Phys Rep 577:1–34

    Article  MathSciNet  MATH  Google Scholar 

  35. Shao S, Lu T, Cai W (2011) Adaptive conservative cell average spectral element methods for transient Wigner equation in quantum transport. Commun Comput Phys 9:711–739

    Article  MathSciNet  MATH  Google Scholar 

  36. Shao S, Sellier JM (2015) Comparison of deterministic and stochastic methods for time-dependent Wigner simulations. J Comput Phys 300:167–185

    Article  MathSciNet  MATH  Google Scholar 

  37. Sloan IH, Kachoyan PJ (1987) Lattice methods for multiple integration: theory, error analysis and examples, SIAM. J Numer Anal 24:116–128

    Article  MathSciNet  MATH  Google Scholar 

  38. Sloan IH, Joe S (1994) Lattice methods for multiple integration. Oxford University Press, Oxford

    MATH  Google Scholar 

  39. Sobol I (1973) Numerical methods Monte Carlo. Nauka, Moscow

    MATH  Google Scholar 

  40. Sobol IM (1967) “Distribution of points in a cube and approximate evaluation of integrals”. Zh. Vych. Mat. Mat. Fiz. 7: 784–802 (in Russian); U.S.S.R Comput. Maths. Math. Phys. 7:86–112

  41. Sobol IM (1989) Quasi-Monte Carlo methods. In: Dimov IT, Sendov BI (eds) International youth workshop on Monte Carlo methods and parallel algorithms. World Scientific, Singapore, pp 75–81

    Google Scholar 

  42. Xiong Y, Chen Z, Shao S (2016) An advective-spectral-mixed method for time-dependent many-body Wigner simulations. SIAM J Sci Comput, to appear, [arXiv:1602.08853]

  43. Wang Y, FJ Hickernell (2002) An historical overview of lattice point sets

  44. Wigner E (1932) On the quantum correction for thermodynamic equilibrium. Phys Rev 40:749

    Article  MATH  Google Scholar 

Download references

Author information

Authors and Affiliations

  1. Information Modeling Department, Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. Georgi Bonchev Str., Block 8, 1113, Sofia, Bulgaria

    Venelin Todorov

  2. Department of Parallel Algorithms, Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, Acad. Georgi Bonchev Str., Block 25 A, 1113, Sofia, Bulgaria

    Venelin Todorov, Ivan Dimov & Rayna Georgieva

  3. Department of Mathematics, Statistics and Computer Science, University of Illinois at Chicago, 851 S Morgan St, Chicago, IL, 60607, USA

    Stoyan Dimitrov

Authors
  1. Venelin Todorov
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ivan Dimov
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Rayna Georgieva
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Stoyan Dimitrov
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Venelin Todorov.

Ethics declarations

Conflict of interest

The authors declare that there is no conflict of interest.

Additional information

Publisher's Note

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

Venelin Todorov is supported by the Bulgarian National Science Fund under Project KP-06-PM32/4—2019 “Advanced Stochastic and Deterministic Approaches for Large-Scale Problems of Computational Mathematics” and by the National Scientific Program “Information and Communication Technologies for a Single Digital Market in Science, Education and Security (ICT in SES)”, contract No DO1-205/23.11.2018, financed by the Ministry of Education and Science in Bulgaria. Stoyan Dimitrov is supported by the Project KP-06-H22/3 financed by the National Scientific Fund of Bulgaria at the Ministry of Education and Science. The work is also partially supported by the Bulgarian National Science Fund under Project DN 12/5-2017. The authors are grateful to Dr. Jean Michel Sellier for the valuable discussions concerning quantum mechanics and especially Wigner function.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Todorov, V., Dimov, I., Georgieva, R. et al. Adaptive Monte Carlo algorithm for Wigner kernel evaluation. Neural Comput & Applic 32, 9953–9964 (2020). https://doi.org/10.1007/s00521-019-04519-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00521-019-04519-9

Keywords

Mathematics Subject Classification

Search

Navigation