Related papers: Error estimate of the u-series method for molecula…
The effective on-site Coulomb interaction (Hubbard $U$) between localized \textit{d} electrons in 3\textit{d}, 4\textit{d}, and 5\textit{d} transition metals is calculated employing a new parameter-free realization of the constrained…
Large molecular dynamics simulations (millions of atoms, tens of microseconds, thousands of processors) hit the strong scalability wall: simulation on twice as many processors does not take half the time. Inspired by large N-body space…
In this study, we introduce a refined method for ascertaining error estimations in numerical simulations of dynamical systems via an innovative application of composition techniques. Our approach involves a dual application of a basic…
The characterization of errors in a quantum system is a fundamental step for two important goals. First, learning about specific sources of error is essential for optimizing experimental design and error correction methods. Second,…
Simulations of asteroid binaries commonly use mutual gravitational potentials approximated by series expansions, leading to truncation errors, and also preventing correct computations of the forces and torques when the bodies are close. We…
We review a family of local algorithms that permit the simulation of charged particles with purely local dynamics. Molecular dynamics formulations lead to discretizations similar to those of ``particle in cell'' methods in plasma physics.…
Sources of energy errors resulting from the replacement of the physical Coulomb interaction by its long-range $\mathrm{erfc}(\mu r)/r$ approximation are explored. It is demonstrated that the results can be dramatically improved and the…
Machine-learning models of atomic-scale interactions achieve the accuracy of the quantum mechanical calculations on which they are trained, but at a dramatically lower computational cost. Their predictions can be made trustworthy by…
Parametric stochastic simulators are ubiquitous in science, often featuring high-dimensional input parameters and/or an intractable likelihood. Performing Bayesian parameter inference in this context can be challenging. We present a neural…
Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterise uncertainty in model inputs and how…
This work presents a technique for statistically modeling errors introduced by reduced-order models. The method employs Gaussian-process regression to construct a mapping from a small number of computationally inexpensive `error indicators'…
Representing a strongly interacting multi-particle wave function in a finite product basis leads to errors. Simple rescaling of the contact interaction can preserve the low-lying energy spectrum and long-wavelength structure of wave…
We develop a random batch Ewald (RBE) method for molecular dynamics simulations of particle systems with long-range Coulomb interactions, which achieves an $O(N)$ complexity in each step of simulating the $N$-body systems. The RBE method is…
The method of intermediate-energy Coulomb excitation has been widely used to determine absolute B(E2; 0+ -> 2+) quadrupole excitation strengths in exotic nuclei with even numbers of protons and neutrons. Transition rates measured with…
We develop a quantum algorithm for estimating the free energy as well as the total Gibbs state of interacting quantum Coulomb gases and molecular systems in dimensions $d \in \{2,3\}$ at finite temperature. These systems lie beyond the…
We deduce a procedure to apply balanced truncation to parameter-dependent differential-algebraic systems. For that we solve multiple projected Lyapunov equations for different parameter values to compute the Gramians that are required for…
We study the relationship between universal effects and the Coulomb interaction. Here, we present our approach and a first illustrative example. Understanding the relationship of universality and the Coulomb interaction is important for…
We consider dynamically constrained Monte-Carlo dynamics and show that this leads to the generation of long ranged effective interactions. This allows us to construct a local algorithm for the simulation of charged systems without ever…
Gibbs measures, such as Coulomb gases, are popular in modelling systems of interacting particles. Recently, we proposed to use Gibbs measures as randomized numerical integration algorithms with respect to a target measure $\pi$ on $\mathbb…
A number of problems arise when long-range forces, such as those governed by Bessel functions, are used in particle-particle simulations. If a simple cut-off for the interaction is used, the system may find an equilibrium configuration at…