Related papers: General Relativistic Implicit Monte Carlo Radiatio…
We present quasi-exact ab initio path integral Monte Carlo (PIMC) results for the partial static density responses and local field factors of hydrogen in the warm dense matter regime, from solid density conditions to the strongly compressed…
We study the implicit Langevin Monte Carlo (iLMC) method, which simulates the overdamped Langevin equation via an implicit iteration rule. In many applications, iLMC is favored over other explicit schemes such as the (explicit) Langevin…
Recently, iterative Quasi-Monte Carlo (iQMC) was introduced as a new method of neutron transport which combines deterministic iterative methods and quasi-Monte Carlo simulation for more efficient solutions to the neutron transport equation.…
We investigate the ionization of the diffuse interstellar medium by cosmic rays by modeling their propagation along the wandering magnetic fields using a Monte Carlo method. We study how low-energy cosmic rays propagate in turbulent,…
Lattice simulations are an important class of problems in crystalline solids, surface science, alloys, adsorption, absorption, separation, catalysis, to name a few. We describe a fast computational method for performing lattice…
The Iterative Quasi-Monte Carlo method, or iQMC, replaces standard quadrature techniques used in deterministic linear solvers with Quasi-Monte Carlo simulation for more accurate and efficient solutions to the neutron transport equation.…
We present PEACOCK, a three-dimensional Monte Carlo radiative transfer (RT) framework designed to self-consistently model rest-frame ultraviolet emission and absorption lines arising from multiphase, clumpy galactic winds. Applied to deep…
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…
Radiative transfer and radiation hydrodynamics use the relativistic Boltzmann equation to describe the kinetics of photons. It is difficult to solve the six-dimensional time-dependent transfer equation unless the problem is highly symmetric…
Radiative processes such as synchrotron radiation and Compton scattering play an important role in astrophysics. Radiative processes are fundamentally stochastic in nature, and the best tools currently used for resolving these processes…
We propose a second-order accurate semi-implicit and well-balanced finite volume scheme for the equations of ideal magnetohydrodynamics (MHD) including gravitational source terms. The scheme treats all terms associated with the acoustic…
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…
We present a formulation and numerical algorithm to extend the scheme for grey radiation magneto-hydrodynamics (MHD) developed by Jiang (2021) to include the frequency dependence via the multi-group approach. The entire frequency space can…
We present the implementation of an implicit-explicit (IMEX) Runge-Kutta numerical scheme for general relativistic hydrodynamics coupled to an optically thick radiation field in two existing GR-hydrodynamics codes. We argue that the…
We present a numerical implementation for the solution of the relativistic radiation hydrodynamics and magnetohydrodynamics equations, designed as an independent module within the freely available code PLUTO. The radiation transfer…
We have developed a new massively-parallel radiation-hydrodynamics code (Cosmos) for Newtonian and relativistic astrophysical problems that also includes radiative cooling, self-gravity, and non-equilibrium, multi-species chemistry. Several…
For a long time, people have been focusing on how to extract more information, such as off-diagonal observables, from the quantum Monte Carlo (QMC) simulation of the partition function, but there have been numerous difficulties, and many of…
Modeling correctly the transport of neutrinos is crucial in some astrophysical scenarios such as core-collapse supernovae and binary neutron star mergers. In this paper, we focus on the truncated-moment formalism, considering only the first…
Radiative transfer describes the propagation of electromagnetic radiation through an interacting medium. This process is often simulated by the use of the Monte Carlo method, which involves the probabilistic determination and tracking of…
We present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation…