Related papers: Evaluating radiation transport errors in merger si…
Monte Carlo radiative transfer (MCRT) simulations are a powerful tool for determining the appearance of astrophysical objects, analyzing the prevalent physical conditions within them, and inferring their properties on the basis of real…
The detection of GW230529_181500 suggested the existence of more symmetric black hole-neutron star mergers where the black hole mass can be as low as 2.6 times that of the neutron star. Black hole-neutron star binaries with even more…
The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents…
A Monte Carlo model of radiative transfer in multi-phase dusty media is applied to the situation of stars and clumpy dust in a sphere or a disk. The distribution of escaping and absorbed photons are shown for various filling factors and…
In the standard approaches to neutrino transport in the simulation of core-collapse supernovae, one will often start from the classical Boltzmann equation for the neutrino's spatial, temporal, and spectral evolution. For each neutrino…
Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, "replica" ensemble of walkers, whose population evolves in…
Monte Carlo criticality simulations are widely used in nuclear safety demonstrations, as they offer an arbitrarily precise estimation of global and local tallies while making very few assumptions. However, since the inception of such…
The extraction efficiency of ultracold neutrons from cryogenic moderators depends critically on the neutron transparency of the moderator material. The Monte Carlo simulation of the probability of the UCN going out from non-ideal…
Neutrinos play an important role in compact star astrophysics: neutrino-heating is one of the main ingredients in core-collapse supernovae, neutrino-matter interactions determine the composition of matter in binary neutron star mergers and…
Inspired by the Boltzmann kinetics, we propose a collision-based dynamics with a Monte Carlo solution algorithm that approximates the solution of the multi-marginal optimal transport problem via randomized pairwise swapping of sample…
General relativistic multi-group and multi-flavor Boltzmann neutrino transport in spherical symmetry adds a new level of detail to the numerical bridge between microscopic nuclear and weak interaction physics and the macroscopic evolution…
Bearing in mind the application to the collapsar models of gamma-ray bursts (GRBs), we develop a numerical scheme and code for estimating the deposition of energy and momentum due to the neutrino pair annihilation ($\nu + {\bar \nu}…
We present new two-dimensional (2D) axisymmetric neutrino radiation/hydrodynamic models of core-collapse supernova (CCSN) cores. We use the CASTRO code, which incorporates truly multi-dimensional, multi-group, flux-limited diffusion (MGFLD)…
Monte Carlo methods are widely used for neutron transport simulations at least partly because of the accuracy they bring to the modeling of these problems. However, the computational burden associated with the slow convergence rate of Monte…
Extending the use of Monte Carlo (MC) event generators to jets in nuclear collisions requires a probabilistic implementation of the non-abelian LPM effect. We demonstrate that a local, probabilistic MC implementation based on the concept of…
We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated…
In this work we investigate replacing standard quadrature techniques used in deterministic linear solvers with a fixed-seed Quasi-Monte Carlo calculation to obtain more accurate and efficient solutions to the neutron transport equation…
The flow inside remnants of binary neutron star (NS) mergers is expected to be turbulent, because of magnetohydrodynamics instability activated at scales too small to be resolved in simulations. To study the large-scale impact of these…
Neutrino mass constraints are a primary focus of current and future large-scale structure (LSS) surveys. Non-linear LSS models rely heavily on cosmological simulations -- the impact of massive neutrinos should therefore be included in these…
Pairing correlations in nuclei play a decisive role in determining nuclear drip-lines, binding energies, and many collective properties. In this work a new Configuration-Space Monte-Carlo (CSMC) method for treating nuclear pairing…