Related papers: Unbiased Monte Carlo continuum radiative transfer …
Aims: We present POLARIS (POLArized RadIation Simulator), a newly developed three-dimensional Monte-Carlo radiative transfer code. POLARIS was designed to calculate dust temperature, polarization maps, and spectral energy distributions. It…
An algorithm for creating synthetic telescope images of Smoothed Particle Hydrodynamics (SPH) density fields is presented, which utilises the adaptive nature of the SPH formalism in full. The imaging process uses Monte Carlo Radiative…
Sequential Monte Carlo (SMC), or particle filtering, is widely used in nonlinear state-space systems, but its performance often suffers from poorly approximated proposal and state-transition distributions. This work introduces a…
We simulate convection near the solar surface, where the continuum optical depth is of order unity. Hence, to determine the radiative heating and cooling in the energy conservation equation, we must solve the radiative transfer equation…
Hamiltonian Monte Carlo is a prominent Markov Chain Monte Carlo algorithm, which employs symplectic integrators to sample from high dimensional target distributions in many applications, such as statistical mechanics, Bayesian statistics…
Over recent years, deep learning based image registration has achieved impressive accuracy in many domains, including medical imaging and, specifically, human neuroimaging with magnetic resonance imaging (MRI). However, the uncertainty…
Monte Carlo simulation was carried out to understand the influence of morphological inhomogeneity on carrier diffusion in organic thin films. The morphological inhomogeneity was considered in the simulation by incorporating the regions of…
Aims. Monte Carlo Radiative Transfer (MCRT) simulations are a powerful tool for understanding the role of dust in astrophysical systems and its influence on observations. However, due to the strong coupling of the radiation field and medium…
We develop a method of analyzing radio frequency spectral line observations to derive data on the temperature, density, velocity, and molecular abundance of the emitting gas. The method incorporates a radiative transfer code with a new…
We describe and analyze some Monte Carlo methods for manifolds in Euclidean space defined by equality and inequality constraints. First, we give an MCMC sampler for probability distributions defined by un-normalized densities on such…
Artifacts arise in the simulations of electrolytes using periodic boundary conditions (PBC). We show the origin of these artifacts are the periodic image charges and the constraint of charge neutrality inside the simulation box, both of…
Sub-millimetre dust emission is often used to derive the column density N of dense interstellar clouds. The observations consist of data at several wavelengths but of variable resolution. We examine two procedures that been proposed for the…
Propagation of charged cosmic-rays in the Galaxy depends on the transport parameters, whose number can be large depending on the propagation model under scrutiny. A standard approach for determining these parameters is a manual scan,…
We develop exact Markov chain Monte Carlo methods for discretely-sampled, directly and indirectly observed diffusions. The qualification "exact" refers to the fact that the invariant and limiting distribution of the Markov chains is the…
In this paper, we present a sparse grid-based Monte Carlo method for solving high-dimensional semi-linear nonlocal diffusion equations with volume constraints. The nonlocal model is governed by a class of semi-linear partial…
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…
We derive analytical expressions to describe how light is frequency-shifted when interacting with ultrasound within scattering media, due to the modulation of the refractive index induced by the ultrasound pressure waves. The model is…
Preliminary low-thrust spacecraft mission design is a global search problem characterized by a complex solution landscape, multiple objectives, and numerous local minima. During this phase, mission parameters are often not yet fully…
We present an ``equation-free'' multiscale approach to the simulation of unsteady diffusion in a random medium. The diffusivity of the medium is modeled as a random field with short correlation length, and the governing equations are cast…
Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we…