Monte Carlo errors with less errors
read the original abstract
We explain in detail how to estimate mean values and assess statistical errors for arbitrary functions of elementary observables in Monte Carlo simulations. The method is to estimate and sum the relevant autocorrelation functions, which is argued to produce more certain error estimates than binning techniques and hence to help toward a better exploitation of expensive simulations. An effective integrated autocorrelation time is computed which is suitable to benchmark efficiencies of simulation algorithms with regard to specific observables of interest. A Matlab code is offered for download that implements the method. It can also combine independent runs (replica) allowing to judge their consistency.
This paper has not been read by Pith yet.
Forward citations
Cited by 4 Pith papers
-
Scalable Generative Sampling and Multilevel Estimation for Lattice Field Theories Near Criticality
A hierarchical generative model for critical lattice scalar field theories achieves orders-of-magnitude lower autocorrelation times than HMC while enabling exact multilevel Monte Carlo.
-
Bounding statistical errors in lattice field theory simulations
Introduces bounds-based stopping criteria and automatic windowing for autocorrelation integrals to estimate statistical errors in lattice field theory Monte Carlo simulations.
-
SU(2) gauge theory with one and two adjoint fermions towards the continuum limit
Extended lattice simulations yield continuum-limit anomalous dimensions γ* = 0.170(6) for Nf=1 and γ* = 0.291(9) for Nf=2 adjoint SU(2), with chiral perturbation theory ruling out spontaneous chiral symmetry breaking.
-
$f_K/f_{\pi}$ in iso-symmetric QCD and the CKM matrix unitarity
Lattice QCD result for f_K/f_π in isoQCD gives |V_us|/|V_ud| consistent with CKM unitarity once isospin and QED effects are included.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.