Related papers: Cosmic Rays Propagation with HelMod: Difference be…
The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and…
We study the time-dependent modulation effect and derive the local interstellar spectra (LIS) for the cosmic ray (CR) proton, helium, boron and carbon. A two-dimensional modulation model including the variation of the interplanetary…
Light propagation in turbid media is driven by the equation of radiative transfer. We give a formal probabilistic representation of its solution in the framework of biological tissues and we implement algorithms based on Monte Carlo methods…
The ACE-CRIS collaboration has recently released the measurement of radioactive $^{60}$Fe nuclei abundance in Galactic Cosmic Rays, in the energy range $\sim 195-500$ MeV per nucleon. We model Cosmic Ray propagation and derive from this…
Numerical solutions of the cosmic-ray (CR) magneto-hydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by…
We introduce Preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilises a Normalising Flow (NF) in order to…
An important factor limiting our ability to understand the production and propagation of cosmic rays pertains to the effects of heliospheric forces, commonly known as solar modulation. The solar wind is capable of generating time and…
We present a Monte-Carlo (MC) calculation of the propagation of cosmic ray protons in the Galaxy for energies above 1 PeV. We discuss the relative strengths of competing effects such as parallel/perpendicular diffusion and drifts in toy…
The Diffusion Monte Carlo method is devoted to the computation of electronic ground-state energies of molecules. In this paper, we focus on implementations of this method which consist in exploring the configuration space with a {\bf fixed}…
We introduce a new Monte Carlo method by incorporating a guided distribution function to the conventional Monte Carlo method. In this way, the efficiency of Monte Carlo methods is drastically improved. To further speed up the algorithm, we…
We present a new model of the heliospheric interface - the region of the solar wind interaction with the local interstellar medium. This new model performs a multi-component treatment of charged particles in the heliosphere. All charged…
Radiative transfer in curved spacetimes has become increasingly important to understanding high-energy astrophysical phenomena and testing general relativity in the strong field limit. The equations of radiative transfer are physically…
We use cosmography to present constraints on the kinematics of the Universe, without postulating any underlying theoretical model. To this end, we use a Monte Carlo Markov Chain analysis to perform comparisons to the supernova Ia Union 2…
We discuss the efficiency of Monte Carlo methods in solving continuum radiative transfer problems. The sampling of the radiation field and convergence of dust temperature calculations in the case of optically thick clouds are both studied.…
Particle methods based on evolving the spatial derivatives of the solution were originally introduced to simulate reaction-diffusion processes, inspired by vortex methods for the Navier--Stokes equations. Such methods, referred to as…
We consider an alternative to the Monte Carlo method for dust continuous radiative transfer simulations: the Quasi-Monte Carlo method. We briefly discuss what it is, its history, and possible implementations. We compare the Monte Carlo…
In this paper, we propose and analyze a new stochastic homogenization method for diffusion equations with random and fast oscillatory coefficients. In the proposed method, the homogenized solutions are sought through a two-stage procedure.…
During the past years several variance reduction techniques for Monte Carlo electron transport have been developed in order to reduce the electron computation time transport for absorbed dose distribution. We have implemented the Macro…
We solve numerically exactly a simple toy model to quantum general relativity or more properly to path integral on a curved space. We consider the thermal equilibrium of a quantum many body problem on the sphere, the surface of constant…
Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration.…