Commun. Comput. Phys., 15 (2014), pp. 1352-1367. |
The Lognormal Distribution and Quantum Monte Carlo Data Mervlyn Moodley ^{1*} 1 School of Chemistry and Physics, Quantum Research Group, University of KwaZulu-Natal, Westville Campus, Private Bag X54001, Durban, 4000, South Africa.Received 19 March 2013; Accepted (in revised version) 17 October 2013 Available online 26 February 2014 doi:10.4208/cicp.190313.171013a Abstract Quantum Monte Carlo data are often afflicted with distributions that resemble lognormal probability distributions and consequently their statistical analysis cannot be based on simple Gaussian assumptions. To this extent a method is introduced to estimate these distributions and thus give better estimates to errors associated with them. This method entails reconstructing the probability distribution of a set of data, with given mean and variance, that has been assumed to be lognormal prior to undergoing a blocking or renormalization transformation. In doing so, we perform a numerical evaluation of the renormalized sum of lognormal random variables. This technique is applied to a simple quantum model utilizing the single-thread Monte Carlo algorithm to estimate the ground state energy or dominant eigenvalue of a Hamiltonian matrix. Notice: Undefined variable: ams in /var/www/html/readabs.php on line 164 PACS: 02.70.Ss, 02.50.Cw, 02.50.Ng, 02.50.-r Key words: Lognormal distribution, renormalization transformation, quantum Monte Carlo. *Corresponding author. Email: moodleym2@ukzn.ac.za (M. Moodley) |