Related papers: Loading a relativistic kappa distribution in parti…
Monte-Carlo simulation of physical processes is an important tool for detector development as it allows to predict signal pulse amplitude and timing, time resolution, efficiency ... Yet despite the fact they are very common, full…
We extend the recently-developed explicit, energy-conserving particle-in-cell (PIC) scheme of [1] to the relativistic Vlasov-Maxwell system. As in the non-relativistic case, the method is built on an optimization problem that is…
Recently popularized randomized methods for principal component analysis (PCA) efficiently and reliably produce nearly optimal accuracy --- even on parallel processors --- unlike the classical (deterministic) alternatives. We adapt one of…
This paper studies sufficient conditions for deriving the kappa distribution in polytropic plasmas by an improved method compared with the previous work [R. Guo, Phys. Plasmas \textbf{27}, 122104 (2020)]. We find that the polytropic…
In this paper we develop an encounter-based model of a run-and-tumble particle (RTP) confined to a finite interval $[0,L]$ with partially absorbing, sticky boundaries at both ends. We assume that the particle switches between two constant…
We investigate the properties of a sequential Monte Carlo method where the particle weight that appears in the algorithm is estimated by a positive, unbiased estimator. We present broadly-applicable convergence results, including a central…
The mutual derivation between arbitrary distribution forms of momenta and momentum components of particles produced in an isotropic emission source are systematically studied in terms of probability theory and mathematical statistics. The…
Increasing shares of fluctuating renewable energy sources induce higher and higher power flow variability at the transmission level. The question arises as to what extent existing networks can absorb additional fluctuating power injection…
A computational procedure is developed for determining the conversion probability for reaction-diffusion systems in which a first-order catalytic reaction is performed over active particles. We apply this general method to systems on metric…
We present an efficient Monte Carlo method to simulate reaction-diffusion processes with spatially varying particle annihilation or transformation rates as it occurs for instance in the context of motor-driven intracellular transport. Like…
We consider the problem of sampling a high dimensional multimodal target probability measure. We assume that a good proposal kernel to move only a subset of the degrees of freedoms (also known as collective variables) is known a priori.…
We present a new Monte Carlo scheme for the efficient simulation of multi-polymer systems. The method permits chains to be inserted into the system using a biased growth technique. The growth proceeds via the use of a retractable feeler,…
The $\Delta$-scaling method has been applied to the total multiplicity distribution of the relativistic ion collisions of p+p, C+C and Pb+Pb which were simulated by a Monte Carlo package, LUCIAE 3.0. It is found that the $\Delta$-scaling…
Levy walk at the finite velocity is considered. To analyze the spatial and temporal characteristics of this process, the method of moments has been used. The asymptotic distributions of the moments (at $t\to\infty$) have been obtained for…
We study semiparametric factor models in high-dimensional panels where the factor loadings consist of a nonparametric component explained by observed covariates and an idiosyncratic component capturing unobserved heterogeneity. A key…
In this paper, we suggest a novel sampling method for Monte Carlo molecular simulations. In order to perform efficient sampling of molecular systems, it is advantageous to avoid extremely high energy configurations while also retaining the…
In the study of phase transitions a very few models are accessible to exact solution. In the most cases analytical simplifications have to be done or some numerical technique has to be used to get insight about their critical properties.…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
We derive an equation for energy transfer from relativistic charged particles to a cold background plasma appropriate for finite-size particles that are used in particle-in-cell simulation codes. Expressions for one-, two-, and…
A number of problems in a variety of fields are characterised by target distributions with a multimodal structure in which the presence of several isolated local maxima dramatically reduces the efficiency of Markov Chain Monte Carlo…