Related papers: Mocassin: A fully three-dimensional Monte Carlo ph…
This paper studies a new and highly efficient Markov chain Monte Carlo (MCMC) methodology to perform Bayesian inference in low-photon imaging problems, with particular attention to situations involving observation noise processes that…
A novel hybrid Monte Carlo transport scheme is demonstrated in a scene with solar illumination, scattering and absorbing 2D atmosphere, a textured reflecting mountain, and a small detector located in the sky (mounted on a satellite or a…
In order to meet the challenges of new accelerator and space projects and further improve modelling of radiation effects in microscopic objects, heavy-ion interaction and transport physics have been recently incorporated into the MARS15…
We develop a Monte Carlo scheme for sampling series of Feynman diagrams for the proper self-energy which are self-consistently expressed in terms of renormalized particle propagators. This approach is used to solve the problem of a single…
We present the analysis of photospheric emission for a set of hydrodynamic simulations of long duration gamma-ray burst jets from massive compact stars. The results are obtained by using the Monte Carlo Radiation Transfer code (MCRaT) to…
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…
In the field of corona discharges, the complex chemical mechanisms inside the ionization region have prompted the development of simplified models to replicate the macroscopic effects of ion generation, thereby reducing the computational…
Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission.…
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…
Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations…
We present a new implementation of the tracer particles algorithm based on a Monte Carlo approach for the Eulerian adaptive mesh refinement code Ramses. The purpose of tracer particles is to keep track of where fluid elements originate in…
We revisit the case of fast Monte-Carlo simulations of galaxy positions for a non-gaussian field. More precisely we address the question of generating a 3D field with a given one-point function (as a log-normal one, but not only) and some…
We present the Monte Carlo integration code MoCaNLO, which computes cross sections and distributions for processes at high-energy colliders like the LHC at leading and next-to-leading order (NLO) in the strong and electroweak couplings. It…
Starting from the radiative transfer equation and its usual boundary conditions, the objective of this work is to design Monte Carlo algorithms estimating the specific intensity spatial and angular derivatives as well as its geometric…
We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust…
The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close…
Markov Chain Monte Carlo based Bayesian data analysis has now become the method of choice for analyzing and interpreting data in almost all disciplines of science. In astronomy, over the last decade, we have also seen a steady increase in…
One of the most demanding calculations is to generate random samples from a specified probability distribution (usually with an unknown normalizing prefactor) in a high-dimensional configuration space. One often has to resort to using a…
We develop a numerical hydrodynamics code using a pseudo-Newtonian formulation that uses the weak field approximation for the geometry, and a generalized source term for the Poisson equation that takes into account relativistic effects. The…
We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all…