Related papers: Acceleration time scale in an ultrarelativistic sh…
The rate at which particles are accelerated by the first-order Fermi mechanism in shocks depends on the angle, \teq{\Tbone}, that the upstream magnetic field makes with the shock normal. The greater the obliquity the greater the rate, and…
Electron acceleration to non-thermal, ultra-relativistic energies (~ 10-100 TeV) is revealed by radio and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock acceleration (DSA) mechanism is usually invoked…
Finite detector resolution and limited acceptance require to apply unfolding methods in high energy physics experiments. Information on the detector resolution is usually given by a set of Monte Carlo events. Based on the experience with a…
We explore an analytic model of the accretion shock in the post bounce phase of a core-collapse supernova explosion. We find growing oscillations of the shock in the l=1 and l=2 modes, in agreement with a variety of existing numerical…
We develop a new Monte Carlo method that solves hyperbolic transport equations with stiff terms, characterized by a (small) scaling parameter. In particular, we focus on systems which lead to a reduced problem of parabolic type in the limit…
Due to a hard dependency between time steps, large-scale simulations of gas using the Direct Simulation Monte Carlo (DSMC) method proceed at the pace of the slowest processor. Scalability is therefore achievable only by ensuring that the…
A Monte Carlo method is presented to evaluate quantum states with many particles moving in the continuum. The scattering state is generated at each time by a Monte Carlo random sampling algorithm. The same calculation are repeated until the…
We describe the large-time moment asymptotics for the parabolic Anderson model where the speed of the diffusion is coupled with time, inducing an acceleration or deceleration. We find a lower critical scale, below which the mass flow gets…
Adaptive Kinetic Monte Carlo combines the simplicity of Kinetic Monte Carlo (KMC) with a Molecular Dynamics (MD) based saddle point search algorithm in order to simulate metastable systems. Key to making Adaptive KMC effective is a stopping…
Particle acceleration at plasma shocks appears to be ubiquitous in the universe, spanning systems in the heliosphere, supernova remnants, and relativistic jets in distant active galaxies and gamma-ray bursts. This review addresses some of…
We analyse a multilevel Monte Carlo method for the approximation of distribution functions of univariate random variables. Since, by assumption, the target distribution is not known explicitly, approximations have to be used. We provide an…
The extremes of a stationary time series typically occur in clusters. A primary measure for this phenomenon is the extremal index, representing the reciprocal of the expected cluster size. Both a disjoint and a sliding blocks estimator for…
When implementing Markov Chain Monte Carlo (MCMC) algorithms, perturbation caused by numerical errors is sometimes inevitable. This paper studies how perturbation of MCMC affects the convergence speed and Monte Carlo estimation accuracy.…
Computing systems interacting with real-world processes must safely and reliably process uncertain data. The Monte Carlo method is a popular approach for computing with such uncertain values. This article introduces a framework for…
In this paper we discuss the possibility of using multilevel Monte Carlo (MLMC) methods for weak approximation schemes. It turns out that by means of a simple coupling between consecutive time discretisation levels, one can achieve the same…
This project investigates the applicability of quasi-Monte Carlo methods to Euclidean lattice systems in order to improve the asymptotic error scaling of observables for such theories. The error of an observable calculated by averaging over…
Interval jitter and spike resampling methods are used to analyze the time scale on which temporal correlations occur. They allow the computation of jitter corrected cross correlograms and the performance of an associated statistically…
Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection…
Generalized method of moments estimators based on higher-order moment conditions derived from independent shocks can be used to identify and estimate the simultaneous interaction in structural vector autoregressions. This study highlights…
We investigate accelerated Unruh-deWitt detectors as a model for particle decay. We find non-trivial decay rates, including a pattern of peaks in decay rate that extends to lower accelerations. Applying our model to the alpha decay of…