Related papers: Mocassin: A fully three-dimensional Monte Carlo ph…
We present the first results of our 3D Monte Carlo maser radiative transfer code, used to model the 1612 MHz OH maser shell and the amplification of emission from the stellar radio-photosphere.
Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly…
On-going measurements of the cosmic radiation (nuclear, electronic, and gamma-ray) are shedding new light on cosmic-ray physics. A comprehensive picture of these data relies on an accurate determination of the transport and source…
Computational codes based on the Diffusion Monte Carlo method can be used to determine the quantum state of two-electron systems confined by external potentials of various nature and geometry. In this work, we show how the application of…
A new parallelized simulation code is presented, which uses a Monte Carlo method to determine particle spectra in the KATRIN source. Reaction chains are generated from the decay of tritium within the source. The code includes all relevant…
The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes…
We combine a Monte Carlo radiative transfer code with an SPH code, so that -- assuming thermal equilibrium -- we can calculate dust-temperature fields, spectral energy distributions, and isophotal maps, for the individual time-frames…
We use a 3D Monte Carlo radiative transfer code to study the polarimetric and photometric variability from stationary corotating interaction regions (CIR) in the wind of massive stars. Our CIRs are approximated by Archimedean spirals of…
Comptonization is the process in which photon spectrum changes due to multiple Compton scatterings in the electronic plasma. It plays an important role in the spectral formation of astrophysical X-ray and gamma-ray sources. There are…
We present a set of new numerical methods that are relevant to calculating radiation pressure terms in hydrodynamics calculations, with a particular focus on massive star formation. The radiation force is determined from a Monte Carlo…
The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative…
Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated…
We present a review of the TORUS radiation transfer and hydrodynamics code. TORUS uses a 1-D, 2-D or3-D adaptive mesh refinement scheme to store and manipulate the state variables, and solves the equation of radiative transfer using Monte…
A numerical scheme is proposed for the solution of the three-dimensional radiative transfer equation with variable optical depth. We show that time-dependent ray tracing is an attractive choice for simulations of astrophysical ionization…
HYPERION is a new three-dimensional dust continuum Monte-Carlo radiative transfer code that is designed to be as generic as possible, allowing radiative transfer to be computed through a variety of three-dimensional grids. The main part of…
Monte Carlo PDE solvers have become increasingly popular for solving heat-related partial differential equations in geometry processing and computer graphics due to their robustness in handling complex geometries. While existing methods can…
Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the…
Photoionization models obtained with numerical codes are widely used to study the physics of the interstellar medium (Planetary Nebulae, H II regions, etc). Grid of models are performed to understand what are the effects of the different…
We develop and implement a Gaussian approach to calculate partial cross-sections and asymmetry parameters for molecular photoionization. Optimal sets of complex Gaussian-type orbitals (cGTOs) are first obtained by non-linear optimization,…
We present the theoretical framework and numerical methods we have implemented to solve the problem of the generation and transfer of polarized radiation in spectral lines out of local thermodynamical equilibrium, while accounting for…