Related papers: Testing a hadronic rescattering model for RHIC col…
In this work, we derive particle schemes, based on micro-macro decomposition, for linear kinetic equations in the diffusion limit. Due to the particle approximation of the micro part, a splitting between the transport and the collision part…
For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated…
Reaction rates are a complicated function of molecular interactions, which can be selected from vast chemical design spaces. Seeking the design that optimizes a rate is a particularly challenging problem since the rate calculation for any…
``Glauber'' models are used to calculate geometric quantities in the initial state of heavy ion collisions, such as impact parameter, number of participating nucleons and initial eccentricity. The four RHIC experiments have different…
This paper studies a non-random-walk Markov Chain Monte Carlo method, namely the Hamiltonian Monte Carlo (HMC) method in the context of Subset Simulation used for structural reliability analysis. The HMC method relies on a deterministic…
We review the di-hadron correlation results from RHIC. A consistent physical picture was constructed based on the correlation landscape in $p_T$, $\Delta\phi$, $\Delta\eta$ and particle species. We show that the data are consistent with…
We propose a novel class of Sequential Monte Carlo (SMC) algorithms, appropriate for inference in probabilistic graphical models. This class of algorithms adopts a divide-and-conquer approach based upon an auxiliary tree-structured…
Recent work on leading order multiparton calculations for hadronic collisions is reviewed, with special emphasis on the ALPGEN event generator. Some problems connected with the interface of the partonic events generated via matrix elements…
When a Monte Carlo algorithm is used to evaluate a physical observable A, it is possible to slightly modify the algorithm so that it evaluates simultaneously A and the derivatives $\partial$ $\varsigma$ A of A with respect to each…
We show that particle species dependence of enhanced hadron production at intermediate transverse momentum (p_{T}) for d+Au collisions at RHIC can be understood in terms of the hadronic rescatterings in the final state. A multiphase…
Particle-in-cell methods with stochastic collision models are commonly used to simulate collisional plasma dynamics, with applications ranging from hypersonic flight to semiconductor manufacturing. Code verification of such methods is…
We have developed a Monte Carlo simulation describing the $2 \rightarrow 2$ scatterings of perturbatively produced, non-eikonally propagating high-energy partons with the quarks and gluons of the expanding QCD medium created in…
We propose new way of heavy ion collisions experiment data analysis. We analyze physical parameters of fireball created in RHIC experiment based on Grand Canonical Distribution and different Lattice QCD data available at the moment. Our…
In this article we propose a model based on the Statistical Bootstrap approach to estimate the cross sections of different hadronic reactions up to a few GeV in c.m.s energy. The method is based on the idea, when two particles collide a so…
Despite more than 40 years of research in condensed-matter physics, state-of-the-art approaches for simulating the radial distribution function (RDF) g(r) still rely on binning pair-separations into a histogram. Such methods suffer from…
The effect of hadronic rescattering on $K^{*0}$ resonance yield can be studied by measuring $K^{*0}/K$ ratio as a function of centrality or multiplicity. This study investigates how the size of the system and the chemical composition…
We study Hamiltonian Monte Carlo (HMC) samplers based on splitting the Hamiltonian $H$ as $H_0(\theta,p)+U_1(\theta)$, where $H_0$ is quadratic and $U_1$ small. We show that, in general, such samplers suffer from stepsize stability…
The event topology in relativistic heavy ion collisions is determined by various multi-particle production mechanisms. The simultaneous model treatment of different collective nuclear effects at high energies (such as a hard multi-parton…
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…
This paper investigates a class of algorithms for numerical integration of a function in d dimensions over a compact domain by Monte Carlo methods. We construct a histogram approximation to the function using a partition of the integration…