Related papers: A simple procedure for generating a Kappa distribu…
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…
In this work, we study a protein synthesis degradation process by defining a general mathematical model. Using generating function technique we present a method that allows exact calculation of joint probability distribution of protein…
We give a algorithm for exact sampling from the Bingham distribution $p(x)\propto \exp(x^\top A x)$ on the sphere $\mathcal S^{d-1}$ with expected runtime of $\operatorname{poly}(d, \lambda_{\max}(A)-\lambda_{\min}(A))$. The algorithm is…
We consider random instances of non-convex perceptron problems in the high-dimensional limit of a large number of examples $M$ and weights $N$, with finite load $\alpha = M/N$. We develop a formalism based on replica theory to predict the…
Simulation has become the evaluation method of choice for many areas of distributing computing research. However, most existing simulation packages have several limitations on the size and complexity of the system being modeled. Fine…
A new approach to simulating warm and hot dense matter that combines density functional theory based calculations of the electronic structure to classical molecular dynamics simulations with pair interaction potentials is presented. The new…
The development and the implementation of a Particle-in-Cell code written in the Unified Parallel C (UPC) language for plasma simulations with application to astrophysics and fusion nuclear energy machines are presented. A simple one…
We study the effects of the velocity distribution functions of the plasma particles on the equilibrium charge of dust grains, acquired through inelastic collisions of the particles with the grains. This paper is the second in a series of…
The paper presents a 1D2V electrostatic PIC model with a drift-kinetic description of all particle types aiming at simulating classical longitudinal plasma transport in axially symmetric open traps. The model generalizes the semi-implicit…
We review the derivation of the Kac master equation model for random collisions of particles, its relationship to the Poisson process, and existing algorithms for simulating values from the marginal distribution of velocity for a single…
We present and compare distributed parallelization strategies for the particle-in-Fourier (PIF) schemes used in kinetic plasma simulations. The different strategies are i) domain decomposition, where both the particles and Fourier modes are…
The Earth's magnetosphere represents a natural plasma laboratory that allows us to study the behavior of particle distribution functions in the absence of Coulomb collisions, typically described by the Kappa distributions. We have…
Rejection sampling is a common tool for low dimensional problems ($d \leq 2$), often touted as an "easy" way to obtain valid samples from a distribution $f(\cdot)$ of interest. In practice it is non-trivial to apply, often requiring…
The paper deals with planar segment processes given by a density with respect to the Poisson process. Parametric models involve reference distributions of directions and/or lengths of segments. These distributions generally do not coincide…
We present SPARC: Simulation Package for Ab-initio Real-space Calculations. SPARC can perform Kohn-Sham density functional theory calculations for isolated systems such as molecules as well as extended systems such as crystals and surfaces,…
We present the first fundamental simulation method for the determination of crater morphology due to femtosecond-pulse laser damage. To this end we have adapted the particle-in-cell (PIC) method commonly used in plasma physics for use in…
We analyze here in some detail, the derivation of the Particle and Monte Carlo methods of plasma simulation, such as Particle in Cell (PIC), Monte Carlo (MC) and Particle in Cell / Monte Carlo (PIC/MC) from formal manipulation of transport…
We propose a homotopy sampling procedure, loosely based on importance sampling. Starting from a known probability distribution, the homotopy procedure generates the unknown normalization of a target distribution. In the context of…
We present an efficient method for evaluating random phase errors in phase shifters within photonic integrated circuits, avoiding the computational cost of traditional Monte Carlo simulations. By modeling spatially correlated manufacturing…
The plasma emission process, which is the mechanism for solar type II and type III radio bursts phenomena, is studied by means of particle-in-cell and weak turbulence simulation methods. By plasma emission, it is meant as a loose…