Related papers: Rare Event Sampling Improves Mercury Instability S…
The simulation of extreme Mach astrophysical flows is traditionally viewed through the lens of deterministic positivity-preserving schemes. However, due to Kelvin--Helmholtz instabilities and shock anomalies, the multi-dimensional Euler…
Contemporary scientific studies often rely on the understanding of complex quantum systems via computer simulation. This paper initiates the statistical study of quantum simulation and proposes a Monte Carlo method for estimating…
We derive a formula for the entropy of two dimensional incompressible inviscid flow, by determining the volume of the space of vorticity distributions with fixed values for the moments Q_k= \int_w(x)^k d^2 x. This space is approximated by a…
We consider a class of queries called durability prediction queries that arise commonly in predictive analytics, where we use a given predictive model to answer questions about possible futures to inform our decisions. Examples of…
In plasma edge simulations, kinetic Monte Carlo (MC) is often used to simulate neutral particles and estimate source terms. For large-sized reactors, like ITER and DEMO, high particle collision rates lead to a substantial computational cost…
The uncertainty of Compton backscattering process is studied by virtue of analytical formulas, and the special effects of variant energy spread and energy drift on the systematic uncertainty estimation are also studied with Monte Carlo…
A revision of Stodolkiewicz's Monte-Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first…
We study the problem of predicting rare critical transition events for a class of slow-fast nonlinear dynamical systems. The state of the system of interest is described by a slow process, whereas a faster process drives its evolution and…
Rare events such as financial crashes, climate extremes, and biological anomalies are notoriously difficult to model due to their scarcity and heavy-tailed distributions. Classical deep generative models often struggle to capture these rare…
We present a Monte Carlo model for the structure of low mass (total mass < 25 earth mass) planetary systems that form by the in situ gravitational assembly of planetary embryos into final planets. Our model includes distributions of mass,…
The problem of precise evaluation of perturbative QCD predictions at moderate energies is addressed. In order to improve stability of the predictions with respect to change of the renormalization scheme it is proposed to replace the…
Recent applications of Quantum Monte Carlo (QMC) technique to Fe-based superconductors opened a way to directly verify the applicability of the itinerant scenario for these systems. Fe-based superconductors undergo various instabilities…
We study the phase-space behaviour of nearby trajectories in integrable potentials. We show that the separation of nearby orbits initially diverges very fast, mimicking a nearly exponential behaviour, while at late times it grows linearly.…
Quasi-Monte Carlo (QMC) methods for estimating integrals are attractive since the resulting estimators typically converge at a faster rate than pseudo-random Monte Carlo. However, they can be difficult to set up on arbitrary posterior…
One of the open challenges in quantum computing is to find meaningful and practical methods to leverage quantum computation to accelerate classical machine learning workflows. A ubiquitous problem in machine learning workflows is sampling…
Finding and sampling rare trajectories in dynamical systems is a difficult computational task underlying numerous problems and applications. In this paper we show how to construct Metropolis- Hastings Monte Carlo methods that can…
The efficient calculation of rare-event kinetics in complex dynamical systems, such as the rate and pathways of ligand dissociation from a protein, is a generally unsolved problem. Markov state models can systematically integrate ensembles…
Quantum Monte Carlo (QMC) techniques are used to calculate the one-body density matrix and excitation energies for the valence electrons of bulk silicon. The one-body density matrix and energies are obtained from a Slater-Jastrow wave…
Intractable generative models are models for which the likelihood is unavailable but sampling is possible. Most approaches to parameter inference in this setting require the computation of some discrepancy between the data and the…
The forecast of the time of arrival of a coronal mass ejection (CME) to Earth is of critical importance for our high-technology society and for any future manned exploration of the Solar System. As critical as the forecast accuracy is the…