相关论文: The Monte Carlo photoionization and moving-mesh ra…
We analyze the mixed frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field, and derive the minimal set of evolution equations that includes all terms that are of leading…
We present a Monte Carlo (MC) code we wrote to simulate the photospheric process and to study the photospheric spectrum above the peak energy. Our simulations were performed with a photon to electron ratio $N_{\gamma}/N_{e} = 10^{5}$, as…
Heavy-ion collisions can be simulated by means of comprehensive approaches, to include the many different reaction mechanisms which may contribute. QMD models and their relativistic extensions are examples of these approaches based on Monte…
Complex organic molecules (COMs) have been widely observed in molecular clouds and protostellar environments. One of the formation mechanisms of COMs is radical reactions on the icy grain surface driven by UV irradiation. While many…
The three-dimensional Monte Carlo photoionization code Mocassin has been applied to construct a realistic model of the planetary nebula NGC 3918. Three different geometric models were tried. The effects of the interaction of the diffuse…
While boundary plasmas in present-day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo…
We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE)…
Ongoing MeV telescopes such as INTEGRAL/SPI and Fermi/GBM, and proposed telescopes including the recently accepted COSI and the e-ASTROGAM and AMEGO missions, provide another window in understanding transients. Their signals contain…
Understanding the spatio-temporal evolution of radiolytic species created by high-energy electrons in water underpins key applications from radiotherapy and nuclear safety to environmental processing and electron microscopy. Here, using the…
By depleting the ionizing photon budget available to expand cosmic HII regions, recombining systems (or Lyman limit systems) can have a large impact during (and following) cosmic reionization. Unfortunately, directly resolving such…
We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass…
Emission and absorption lines from elements heavier than helium (metals) represent one of our strongest probes of galaxy formation physics across nearly all redshifts accessible to observations. The vast majority of simulations that model…
An ionization front (I-front) that propagates through an inhomogeneous medium is slowed down by self-shielding and recombinations. We perform cosmological radiation hydrodynamics simulations of the I-front propagation during the epoch of…
Three non-equilibrium photoionisation algorithms for hydrodynamical grid-based simulation codes are compared in terms of accuracy, timestepping criteria, and parallel scaling. Explicit methods with first-order time accuracy for photon…
We describe a major update to the public GIZMO code. GIZMO has been used in simulations of cosmology; galaxy and star formation and evolution; black hole accretion and feedback; proto-stellar disk dynamics and planet formation; fluid…
We present a radiation-hydrodynamics (RHD) scheme that enables 3D simulations resolving both protostellar interiors and their surrounding accretion flows within a single framework, to clarify how a protostar evolves while interacting with…
Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have…
We have developed a modular semi-numerical code that computes the time and spatially dependent ionization of neutral hydrogen (HI), neutral (HeI) and singly ionized helium (HeII) in the intergalactic medium (IGM). The model accounts for…
We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive…
We compare numerical methods for solving the radiative transfer equation in the context of the photoionization of intergalactic gaseous hydrogen and helium by a central radiating source. Direct integration of the radiative transfer equation…