4.7 Article

Automatic differentiation for error analysis of Monte Carlo data

Journal

COMPUTER PHYSICS COMMUNICATIONS
Volume 238, Issue -, Pages 19-35

Publisher

ELSEVIER
DOI: 10.1016/j.cpc.2018.12.020

Keywords

Lattice QCD; Monte Carlo; Error analysis

Ask authors/readers for more resources

Automatic Differentiation (AD) allows to determine exactly the Taylor series of any function truncated at any order. Here we propose to use AD techniques for Monte Carlo data analysis. We discuss how to estimate errors of a general function of measured observables in different Monte Carlo simulations. Our proposal combines the Gamma-method with Automatic differentiation, allowing exact error propagation in arbitrary observables, even those defined via iterative algorithms. The case of special interest where we estimate the error in fit parameters is discussed in detail. We also present a freely available Fortran reference implementation of the ideas discussed in this work. (C) 2019 Elsevier B.V. All rights reserved.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Physics, Particles & Fields

Symanzik improvement of the gradient flow in lattice gauge theories

Alberto Ramos, Stefan Sint

EUROPEAN PHYSICAL JOURNAL C (2016)

Article Physics, Multidisciplinary

QCD Coupling from a Nonperturbative Determination of the Three-Flavor Λ Parameter

Mattia Bruno, Mattia Dalla Brida, Patrick Fritzsch, Tomasz Korzec, Alberto Ramos, Stefan Schaefer, Hubert Simma, Stefan Sint, Rainer Sommer

PHYSICAL REVIEW LETTERS (2017)

Article Physics, Particles & Fields

A non-perturbative exploration of the high energy regime in Nf=3 QCD

Mattia Dalla Brida, Patrick Fritzsch, Tomasz Korzec, Alberto Ramos, Stefan Sint, Rainer Sommer

EUROPEAN PHYSICAL JOURNAL C (2018)

Article Physics, Particles & Fields

Non-perturbative quark mass renormalisation and running in Nf=3 QCD

I. Campos, P. Fritzsch, C. Pena, D. Preti, A. Ramos, A. Vladikas

EUROPEAN PHYSICAL JOURNAL C (2018)

Article Physics, Particles & Fields

The gradient flow coupling at high-energy and the scale of SU(3) Yang-Mills theory

Mattia Dalla Brida, Alberto Ramos

EUROPEAN PHYSICAL JOURNAL C (2019)

Article Physics, Particles & Fields

Light quark masses in Nf=2+1 lattice QCD withWilson fermions ALPHA Collaboration

M. Bruno, I. Campos, P. Fritzsch, J. Koponen, C. Pena, D. Preti, A. Ramos, A. Vladikas

EUROPEAN PHYSICAL JOURNAL C (2020)

Review Physics, Particles & Fields

FLAG Review 2019

S. Aoki, Y. Aoki, D. Becirevic, T. Blum, G. Colangelo, S. Collins, M. Della Morte, P. Dimopoulos, S. Duerr, H. Fukaya, M. Golterman, Steven Gottlieb, R. Gupta, S. Hashimoto, U. M. Heller, G. Herdoiza, R. Horsley, A. Juttner, T. Kaneko, C. -J. D. Lin, E. Lunghi, R. Mawhinney, A. Nicholson, T. Onogi, C. Pena, A. Portelli, A. Ramos, S. R. Sharpe, J. N. Simone, S. Simula, R. Sommer, R. Van de Water, A. Vladikas, U. Wenger, H. Wittig

EUROPEAN PHYSICAL JOURNAL C (2020)

Article Physics, Particles & Fields

openQ*D code: a versatile tool for QCD plus QED simulations:collaboration

Isabel Campos, Patrick Fritzsch, Martin Hansen, Marina Krstic Marinkovic, Agostino Patella, Alberto Ramos, Nazario Tantalo

EUROPEAN PHYSICAL JOURNAL C (2020)

Article Astronomy & Astrophysics

Non-perturbative renormalization by decoupling

Mattia Dalla Brida, Roman Hoellwieser, Francesco Knechtli, Tomasz Korzec, Alberto Ramos, Rainer Sommer

PHYSICS LETTERS B (2020)

Article Physics, Particles & Fields

An analysis of systematic effects in finite size scaling studies using the gradient flow

Alessandro Nada, Alberto Ramos

Summary: The proposal introduces a new strategy for determining the step scaling function sigma (u) in finite size scaling studies using the gradient flow, showing better control over continuum extrapolations in pure gauge theory. Additionally, it reexamines the running coupling at high energies and the determination of the Lambda-parameter, focusing on perturbative truncation uncertainties.

EUROPEAN PHYSICAL JOURNAL C (2021)

Review Physics, Multidisciplinary

Lattice determinations of the strong coupling

Luigi Del Debbio, Alberto Ramos

Summary: Lattice QCD has established itself as a mature field, providing precise descriptions of the standard model and determining essential quantities such as the strong coupling constant. In addition, lattice calculations will be crucial in future phenomenological studies.

PHYSICS REPORTS-REVIEW SECTION OF PHYSICS LETTERS (2021)

Proceedings Paper Physics, Multidisciplinary

Controlling quark mass determinations non-perturbatively in three-flavour QCD

Isabel Campos, Patrick Fritzsch, Carlos Pena, David Preti, Alberto Ramos, Anastassios Vladikas

XIITH QUARK CONFINEMENT AND THE HADRON SPECTRUM (2017)

Article Astronomy & Astrophysics

Slow running of the gradient flow coupling from 200 MeV to 4 GeV in Nf=3 QCD

Mattia Dalla Brida, Patrick Fritzsch, Tomasz Korzec, Alberto Ramos, Stefan Sint, Rainer Sommer

PHYSICAL REVIEW D (2017)

Proceedings Paper Physics, Particles & Fields

The determination of αs by the ALPHA collaboration

Mattia Bruno, Mattia Dalla Brida, Patrick Fritzsch, Tomasz Korzec, Alberto Ramos, Stefan Schaefer, Hubert Simma, Stefan Sint, Rainer Sommer

NUCLEAR AND PARTICLE PHYSICS PROCEEDINGS (2017)

Article Computer Science, Interdisciplinary Applications

Modeling and meshing for tokamak edge plasma simulations

Usman Riaz, E. Seegyoung Seol, Robert Hager, Mark S. Shephard

Summary: The accurate representation and effective discretization of a problem domain into a mesh are crucial for achieving high-quality simulation results and computational efficiency. This work presents recent developments in extending an automated tokamak modeling and meshing infrastructure to better support the near flux field following meshing requirements of the XGC Gyro-kinetic Code.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Electron-phonon coupling from GW perturbation theory: Practical workflow combining BerkeleyGW, ABINIT, and EPW

Zhenglu Li, Gabriel Antonius, Yang-Hao Chan, Steven G. Louie

Summary: This article presents a workflow for practical calculations of electron-phonon coupling and includes the effect of many-electron correlations using GW perturbation theory. The workflow combines different software packages to enable accurate calculations at the level of quasiparticle band structures.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

MASAP: A package for atomic scattering amplitude in solids

Akihiro Koide, Sara Rabouli, Pierre Le Meur, Sylvain Tricot, Philippe Schieffer, Didier Sebilleau, Calogero R. Natoli

Summary: We present the MsSpec Atomic Scattering Amplitude Package (MASAP), which includes a computation program and a graphical interface for generating atomic scattering amplitude (ASA). The study investigates the applicability of plane wave (PW) and curved spherical wave (SW) scattering in describing electron propagation. The results show that the imaginary part of the optical potential enhances the elastic scattering in the forward direction but causes damping effects in other directions.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

A Bi-directional method for evaluating integrals involving higher transcendental functions. HyperRAF: A Julia package for new hyper-radial functions

A. Bagci, Gustavo A. Aucar

Summary: The electron repulsion integrals over Slater-type orbitals with non-integer principal quantum numbers are investigated in this study. These integrals are important in calculations of many-electron systems. New relationships free from hyper-geometric functions are derived to simplify the calculations. With the use of auxiliary functions and straightforward recurrence relationships, these integrals can be efficiently computed, providing initial conditions for the evaluation of expectation values and potentials.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

First-principles calculations of specular reflection of high-energy electrons during the two-dimensional crystal growth

Andrzej Daniluk

Summary: RHEED_DIFF_2D is an open-source software for qualitative numerical simulations of RHEED oscillation intensity changes with layer deposition, used for interpreting heteroepitaxial structures under different scattering crystal potential models.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

An incremental singular value decomposition approach for large-scale spatially parallel & distributed but temporally serial data - applied to technical flows ☆

Niklas Kuehl, Hendrik Fischer, Michael Hinze, Thomas Rung

Summary: The article presents a strategy and algorithm for simulation-accompanying, incremental Singular Value Decomposition (SVD) for time-evolving, spatially parallel discrete data sets. The proposed method improves computational efficiency by introducing a bunch matrix, resulting in higher accuracy and practical applicability.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

TRAVOLTA: GPU acceleration and algorithmic improvements for constructing quantum optimal control fields in photo-excited systems

Jose M. Rodriguez-Borbon, Xian Wang, Adrian P. Dieguez, Khaled Z. Ibrahim, Bryan M. Wong

Summary: This paper presents an open-source software package called TRAVOLTA for massively parallelized quantum optimal control calculations on GPUs. The TRAVOLTA package is an improvement on the previous NIC-CAGE algorithm and incorporates algorithmic improvements for faster convergence. Three different variants of GPU parallelization are examined to evaluate their performance in constructing optimal control fields in various quantum systems. The benchmarks show that the GPU-enhanced TRAVOLTA code produces the same results as previous CPU-based algorithms but with a speedup of more than ten times. The GPU enhancements and algorithmic improvements allow large quantum optimal control calculations to be efficiently executed on modern multi-core computational hardware.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

MCNOX: A code for computing and interpreting ultrafast nonlinear X-ray spectra of molecules at the multiconfigurational level

Weijie Hua

Summary: This work introduces a program called MCNOX for computing and analyzing ultrafast nonlinear X-ray spectra. It is designed for cutting-edge applications in photochemistry/photophysics enabled by X-ray free-electron lasers and high harmonic generation light sources. The program can calculate steady-state X-ray absorption spectroscopy and three types of ultrafast nonlinear X-ray spectra, and it is capable of identifying major electronic transitions and providing physical and chemical insights from complex signals.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

PLQ-sim: A computational tool for simulating photoluminescence quenching dynamics in organic donor/acceptor blends

Leandro Benatto, Omar Mesquita, Lucimara S. Roman, Rodrigo B. Capaz, Graziani Candiotto, Marlus Koehler

Summary: Photoluminescence Quenching Simulator (PLQ-Sim) is a user-friendly software for studying the dynamics of photoexcited states at the interface between organic semiconductors. It provides important information on organic photovoltaic and photothermal devices and calculates transfer rates and quenching efficiency.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

A method of calculating bandstructure in real-space with application to all-electron and full potential

Dongming Li, James Kestyn, Eric Polizzi

Summary: This study introduces a practical and efficient approach to calculate the all-electron full potential band structure in real space using a finite element basis. Instead of the k-space method, this method solves the Kohn-Sham equation self-consistently within a larger finite system enclosing the unit-cell. Non-self-consistent calculations are then performed in the Brillouin zone to obtain the band structure results, which are found to be in excellent agreement with the pseudopotential k-space method. Furthermore, the study successfully observes the band bending of core electrons.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

EUTERPE: A global gyrokinetic code for stellarator geometry

R. Kleiber, M. Borchardt, R. Hatzky, A. Koenies, H. Leyh, A. Mishchenko, J. Riemann, C. Slaby, J. M. Garcia-Regana, E. Sanchez, M. Cole

Summary: This paper describes the current state of the EUTERPE code, focusing on the implemented models and their numerical implementation. The code is capable of solving the multi-species electromagnetic gyrokinetic equations in a three-dimensional domain. It utilizes noise reduction techniques and grid resolution transformation for efficient computation. Additionally, various hybrid models are implemented for comparison and the study of plasma-particle interactions. The code is parallelized for high scalability on multiple CPUs.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

SMIwiz: An integrated toolbox for multidimensional seismic modelling and imaging

Pengliang Yang

Summary: This paper presents an open source software called SMIwiz, which combines seismic modelling, reverse time migration, and full waveform inversion into a unified computer implementation. SMIwiz supports both 2D and 3D simulations and provides various computational recipes for efficient calculation. Its independent processing and batchwise job scheduling ensure scalability, and its viability is demonstrated through applications on benchmark models.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Generating and grading 34 optimised norm-conserving Vanderbilt pseudopotentials for actinides and super-heavy elements in the PseudoDojo

Christian Tantardini, Miroslav Ilias, Matteo Giantomassi, Alexander G. Kvashnin, Valeria Pershina, Xavier Gonze

Summary: Material discovery has been an active research field, and this study focuses on developing pseudopotentials for actinides and super-heavy elements. These pseudopotentials are crucial for accurate first-principles calculations and simulations.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Generalisation of splitting methods based on modified potentials to nonlinear evolution equations of parabolic and Schrödinger type

S. Blanes, F. Casas, C. Gonzalez, M. Thalhammer

Summary: This paper explores the extension of modified potential operator splitting methods to specific classes of nonlinear evolution equations. Numerical experiments confirm the advantages of the proposed fourth-order modified operator splitting method over traditional splitting methods in dealing with Gross-Pitaevskii systems.

COMPUTER PHYSICS COMMUNICATIONS (2024)

Article Computer Science, Interdisciplinary Applications

Pole-fitting for complex functions: Enhancing standard techniques by artificial-neural-network classifiers and regressors *

Siegfried Kaidisch, Thomas U. Hilger, Andreas Krassnigg, Wolfgang Lucha

Summary: Motivated by a use case in theoretical hadron physics, this paper revisits an application of a pole-sum fit to dressing functions of a confined quark propagator. Specifically, it investigates approaches to determine the number and positions of singularities closest to the origin for a function known numerically on a specific grid on the positive real axis. Comparing the efficiency of standard techniques to a pure artificial-neural-network approach and a combination of both, it finds that the combined approach is more efficient. This approach can be applied to similar situations where the positions of poles need to be estimated quickly and reliably from real-axis information alone.

COMPUTER PHYSICS COMMUNICATIONS (2024)