Monte Carlo Radiative Transfer Code Research Articles

Lyα emitters in a cosmological volume – I. The impact of radiative transfer

Lyman-{\alpha} emitters (LAEs) are a promising target to probe the large scale structure of the Universe at high redshifts, $z\gtrsim 2$. However, their detection is sensitive to radiative transfer effects that depend on local astrophysical conditions. Thus, modeling the bulk properties of this galaxy population remains challenging for theoretical models. Here we develop a physically-motivated scheme to predict LAEs in cosmological simulations. The escape of Ly{\alpha} photons is computed using a Monte Carlo radiative transfer code which outputs a Ly{\alpha} escape fraction. To speed-up the process of assigning escape fractions to individual galaxies, we employ fitting formulae that approximate the full Monte Carlo results within an accuracy of 10% for a broad range of column densities, gas metallicities and gas bulk velocities. We apply our methodology to the semi-analytical model GALFORM on a large N-body simulation. The Ly{\alpha} photons escape through an outflowing neutral gas medium, implemented assuming different geometries. This results in different predictions for the typical column density and outflow velocities of the LAE population. To understand the impact of radiative transfer on our predictions, we contrast our models against a simple abundance matching assignment. Our full models populate LAEs in less massive haloes than what is obtained with abundance matching. Overall, radiative transfer effects result in better agreement when confronting the properties of LAEs against observational measurements. This suggest that incorporating the effects of Ly{\alpha} radiative transfer in the analysis of this galaxy population, including their clustering, can be important for obtaining an unbiased interpretation of future datasets.

Predicting the broad-lines polarization emitted by supermassive binary black holes

Context. Some Type-1 active galactic nuclei (AGN) show extremely asymmetric Balmer lines with the broad peak redshifted or blueshifted by thousands of km s−1. These AGN may be good candidates for supermassive binary black holes (SMBBHs). The complex line shapes can be due to the complex kinematics of the two broad line regions (BLRs). Therefore other methods should be applied to confirm the SMBBHs. One of them is spectropolarimetry.Aims. We rely on numerical modeling of the polarimetry of binary black holes systems, since polarimetry is highly sensitive to geometry, in order to find the specific influence of supermassive binary black hole (SMBBH) geometry and dynamics on polarized parameters across the broad line profiles. We apply our method to SMBBHs in which both components are assumed to be AGN with distances at the subparsec scale.Methods. We used a Monte Carlo radiative transfer code that simulates the geometry, dynamics, and emission pattern of a binary system where two black holes are getting increasingly close. Each gravitational well is accompanied by its own BLR and the whole system is surrounded by an accretion flow from the distant torus. We examined the emission line deformation and predicted the associated polarization that could be observed.Results. We modeled scattering-induced broad line polarization for various BLR geometries with complex kinematics. We find that the presence of SMBBHs can produce complex polarization angle profilesφand strongly affect the polarized and unpolarized line profiles. Depending on the phase of the SMBBH, the resulting double-peaked emission lines either show red or blue peak dominance, or both the peaks can have the same intensity. In some cases, the whole line profile appears as a single Gaussian line, hiding the true nature of the source.Conclusions. Our results suggest that future observation with the high resolution spectropolarimetry of optical broad emission lines could play an important role in detecting subparsec SMBBHs.

Spectral modeling of type II supernovae

We present substantial extensions to the Monte Carlo radiative transfer code TARDIS to perform spectral synthesis for type II supernovae. By incorporating a non-LTE ionization and excitation treatment for hydrogen, a full account of free–free and bound-free processes, a self-consistent determination of the thermal state and by improving the handling of relativistic effects, the improved code version includes the necessary physics to perform spectral synthesis for type II supernovae to high precision as required for the reliable inference of supernova properties. We demonstrate the capabilities of the extended version of TARDIS by calculating synthetic spectra for the prototypical type II supernova SN1999em and by deriving a new and independent set of dilution factors for the expanding photosphere method. We have investigated in detail the dependence of the dilution factors on photospheric properties and, for the first time, on changes in metallicity. We also compare our results with the previously published sets of dilution factors and discuss the potential sources of the discrepancies between studies.

Photospheric Emission from Variable Engine Gamma-Ray Burst Simulations

Abstract By coupling radiation transfer calculations to hydrodynamic simulations, there have been major advancements in understanding the long gamma-ray burst (LGRB) prompt emission. Building upon these achievements, we present an analysis of photospheric emission acquired by using the Monte Carlo Radiation Transfer (MCRaT) code on hydrodynamic simulations with variable jet profiles. MCRaT propagates and Compton-scatters individual photons that have been injected into the collimated outflow in order to produce synthetic light curves and spectra. These light curves and spectra allow us to compare our results to LGRB observational data. We find excellent agreement between our fitted time-resolved β parameters and those that are observed. Additionally, our simulations show that photospheric emission, under certain conditions, is able to create the observationally expected Band α parameter. Finally, we show that the simulations are consistent with the Golenetskii correlation but exhibit some strain with the Amati and Yonetoku correlations.

Environmental effects on the ionisation of brown dwarf atmospheres

Context. Brown dwarfs emit bursts of Hα, white-light flares, and show radio flares and quiescent radio emission. They are suggested to form aurorae, similar to planets in the solar system, but much more energetic. All these processes require a source gas with an appropriate degree of ionisation, which, so far, is mostly postulated to be sufficient. Aims. We aim to demonstrate that the Galactic environment influences atmospheric ionisation, and that it hence amplifies or enables the magnetic coupling of the atmospheres of ultra-cool objects, like brown dwarfs and free-floating planets. Methods. We build on our previous work on thermal ionisation of ultra-cool atmospheres and explore the effect of environmental high-energy radiation on the degree of ionisation in the atmosphere. We consider the effect of photoionisation by Lyman-continuum radiation in three different environments: in the interstellar radiation field (ISRF), O and B stars in star-forming regions, and in white dwarf companions in binary systems. We apply our Monte Carlo radiation transfer code to investigate the effect of Lyman-continuum photoionisation for prescribed atmosphere structures for very low-mass objects. Results. The external radiation environment plays an important role for the atmospheric ionisation of very low-mass, ultra-cool objects. Lyman-continuum irradiation greatly increases the level of ionisation in the uppermost atmospheric regions. Our results suggest that a shell of an almost fully ionised atmospheric gas emerges for brown dwarfs in star-forming regions and brown dwarfs in white dwarf binary systems. As a consequence, brown dwarf atmospheres can be magnetically coupled, which is the presumption for chromospheric heating to occur and for aurorae to emerge. First tests for assumed chromosphere-like temperature values suggest that the resulting free-free X-ray luminosities are comparable with those observed from non-accreting brown dwarfs in star-forming regions.

Radiative transfer of ionizing radiation through gas and dust: the stellar source case

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H II regions and on small scale reionization. We find that H II regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions.

Type Iax supernovae as a few-parameter family

We present direct spectroscopic modeling of five Type Iax supernovae (SNe) with the one dimensional Monte Carlo radiative transfer code TARDIS. The abundance tomography technique is used to map the chemical structure and physical properties of the SN atmosphere. Through via fitting of multiple spectral epochs with self-consistent ejecta models, we can then constrain the location of some elements within the ejecta. The synthetic spectra of the best-fit models are able to reproduce the flux continuum and the main absorption features in the whole sample. We find that the mass fractions of IGEs and IMEs show a decreasing trend toward the outer regions of the atmospheres using density profiles similar to those of deflagration models in the literature. Oxygen is the only element, which could be dominant at higher velocities. The stratified abundance structure contradicts the well-mixed chemical profiles predicted by pure deflagration models. Based on the derived densities and abundances, a template model atmosphere is created for the SN Iax class and compared to the observed spectra. Free parameters are the scaling of the density profile, the velocity shift of the abundance template, and the peak luminosity. The results of this test support the idea that all SNe Iax can be described by a similar internal structure, which argues for a common origin of this class of explosions.

The Impact of Neutral Intergalactic Gas on Lyα Intensity Mapping during Reionization

Abstract We present the first simulations of the high-redshift Lyα intensity field that account for scattering in the intergalactic medium (IGM). Using a 3D Monte Carlo radiative transfer code, we find that Lyα scattering smooths spatial fluctuations in the Lyα intensity on small scales, and that the spatial dependence of this smoothing depends on the mean neutral fraction of the IGM. Our simulations find a strong effect of reionization on k = 0.1–1 Mpc−1, with for and for , in contrast to after reionization. At wavenumbers of k > 1 Mpc−1, we find that the signal is sensitive to the emergent Lyα line profiles from galaxies. We also demonstrate that the cross-correlation between a Lyα intensity map and a future galaxy redshift survey could be detected on large scales by an instrument similar to SPHEREx, and over a wide range of scales by a hypothetical intensity mapping instrument in the vein of the Cosmic Dawn Intensity Mapper (CDIM).

A systematic survey of grain growth in discs around post-AGB binaries with PACS and SPIRE photometry

AbstractPost-AGB stars are the final stage of evolution of low-intermediate mass stars (M < 8 Mʘ). Those in binary systems have stable circumbinary discs. Using data from Herschel (PACS/SPIRE), we extend the SEDs of 50 galactic post-AGB binary systems to sub-millimetre wavelengths and use the slope of the SED as a diagnostic tool to probe the presence of large grains. Using a Monte Carlo radiative transfer code (MCMax), we create a large grid of models to quantify the observed spectral indices, and use the presence of large grains in the disc as a proxy for evolution.

Spinning Dust Emission from Circumstellar Disks and Its Role In Excess Microwave Emission

Abstract Electric dipole emission from rapidly spinning polycyclic aromatic hydrocarbons (PAHs) is widely believed to be an origin of anomalous microwave emission (AME), but recently it has encountered a setback owing to the noncorrelation of AME with PAH abundance seen in a full-sky analysis. Microwave observations for specific regions with well-constrained PAH features would be crucial to test the spinning dust hypothesis. In this paper, we present physical modeling of microwave emission from spinning PAHs from protoplanetary disks (PPDs) around Herbig Ae/Be stars and T Tauri stars where PAH features are well observed. Guided by the presence of 10 μm silicate features in some PPDs, we also model microwave emission from spinning nanosilicates. Thermal emission from big dust grains is computed using the Monte Carlo radiative transfer code (radmc-3d). Our numerical results demonstrate that microwave emission from either spinning PAHs or spinning nanosilicates dominates over thermal dust at frequencies ν < 60 GHz, even in the presence of significant grain growth. Finally, we attempt to fit millimeter–centimeter observational data with both thermal dust and spinning dust for several disks around Herbig Ae/Be stars that exhibit PAH features and find that spinning dust can successfully reproduce the observed excess microwave emission (EME). Future radio observations with ngVLA, SKA, and ALMA Band 1 would be valuable for elucidating the origin of EME and potentially open a new window for probing nanoparticles in circumstellar disks.

The role of external far-ultraviolet irradiation in the survival of astrophysical ices in Elias 29

The survival of astrophysical ices in the star-forming regions depends on the suitability of temperature, density and radiation conditions. In this paper, the role of Interstellar Radiation Field (ISRF) on ices in Elias 29 is addressed. This object is the most luminous protostar in $\rho$ Oph E molecular cloud, and is surrounded by many Young Stellar Objects distant only a few arcminutes. In addition, other two brightness BV stars (S1 and HD 147889) enhance the external irradiation in Elias 29. This study was carried out by using the Monte Carlo radiative transfer code RADMC-3D assuming an internal and external irradiation. As result, we found that HD 147889 dominates the ISRF, instead of the closest protostars, and contributes to enhance the external irradiation in 44 times the standard value. Furthermore, remarkable effects are observed in FIR spectrum, as well as in near-IR image. Additionally, the snowline positions of volatile compounds, such as CO, O$_2$, N$_2$ and CH$_4$ are redefined to a toroidal-shaped morphology in the envelope, with low FUV flux (10$^{-7}$ ergs cm$^{-2}$ s$^{-1}$). In such scenario, the formation of complex molecules as result of hydrogenation or oxygenation of volatile species are expected to be severely affected.

Dust modeling of the combined ALMA and SPHERE datasets of HD 163296

Context. Multiwavelength observations are indispensable in studying disk geometry and dust evolution processes in protoplanetary disks. Aims. We aim to construct a three-dimensional model of HD 163296 that is capable of reproducing simultaneously new observations of the disk surface in scattered light with the SPHERE instrument and thermal emission continuum observations of the disk midplane with ALMA. We want to determine why the spectral energy distribution of HD 163296 is intermediary between the otherwise well-separated group I and group II Herbig stars. Methods. The disk was modeled using the Monte Carlo radiative transfer code MCMax3D. The radial dust surface density profile was modeled after the ALMA observations, while the polarized scattered light observations were used to constrain the inclination of the inner disk component and turbulence and grain growth in the outer disk. Results. While three rings are observed in the disk midplane in millimeter thermal emission at ~80, 124, and 200 AU, only the innermost of these is observed in polarized scattered light, indicating a lack of small dust grains on the surface of the outer disk. We provide two models that are capable of explaining this difference. The first model uses increased settling in the outer disk as a mechanism to bring the small dust grains on the surface of the disk closer to the midplane and into the shadow cast by the first ring. The second model uses depletion of the smallest dust grains in the outer disk as a mechanism for decreasing the optical depth at optical and near-infrared wavelengths. In the region outside the fragmentation-dominated regime, such depletion is expected from state-of-the-art dust evolution models. We studied the effect of creating an artificial inner cavity in our models, and conclude that HD 163296 might be a precursor to typical group I sources.

OGLE14-073 – a promising pair-instability supernova candidate

The recently discovered bright type II supernova OGLE14-073 evolved very slowly. The light curve rose to maximum for 90 days from discovery and then declined at a rate compatible with the radioactive decay of 56Co. In this study, we show that a pair-instability supernova is a plausible mechanism for this event. We calculate explosion models and light curves with the radiation hydrodynamics code STELLA starting from two M(ZAMS) = 150 solar masses, Z=0.001 progenitors. We obtain satisfactory fits to OGLE14-073 broadband light curves by including additional 56Ni in the centre of the models and mixing hydrogen down into the inner layers of the ejecta to a radial mass coordinate of 10 solar masses. The extra 56Ni required points to a slightly more massive progenitor star. The mixing of hydrogen could be due to large scale mixing during the explosion. We also present synthetic spectra for our models simulated with the Monte Carlo radiative transfer code ARTIS. The synthetic spectra reproduce the main features of the observed spectra of OGLE14-073. We conclude that OGLE14-073 is one of the most promising candidates for a pair-instability explosion.

Differences in the Gas and Dust Distribution in the Transitional Disk of a Sun-like Young Star, PDS 70

Abstract We present ALMA 0.87 mm continuum, HCO+ J = 4–3 emission line, and CO J = 3–2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two-component transitional disk system with a radial dust gap of 0.″42 ± 0.″05, an azimuthal gap in the HCO+ J = 4–3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analog in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D using a two-disk component. We find that there is a radial gap that extends from 15 to 60 au in all grain sizes, which differs from previous work.

Unlocking the Full Potential of Extragalactic Lyα through Its Polarization Properties

Abstract Lyα is a powerful astrophysical probe. Not only is it ubiquitous at high redshifts, it is also a resonant line, making Lyα photons scatter. This scattering process depends on the physical conditions of the gas through which Lyα propagates, and these conditions are imprinted on observables such as the Lyα spectrum and its surface brightness profile. In this work, we focus on a less-used observable capable of probing any scattering process: polarization. We implement the density matrix formalism of polarization into the Monte Carlo radiative transfer code tlac. This allows us to treat it as a quantum mechanical process where single photons develop and lose polarization from scatterings in arbitrary gas geometries. We explore static and expanding ellipsoids, biconical outflows, and clumpy multiphase media. We find that photons become increasingly polarized as they scatter and diffuse into the wings of the line profiles, making scattered Lyα polarized in general. The degree and orientation of Lyα polarization depends on the kinematics and distribution of the scattering H i gas. We find that it generally probes spatial or velocity space asymmetries and aligns itself tangentially to the emission source. We show that the mentioned observables, when studied separately, can leave similar signatures for different source models. We conclude by revealing how a joint analysis of the Lyα spectra, surface brightness profiles, and polarization can break these degeneracies and help us extract unique physical information on galaxies and their environments from their strongest, most prominent emission line.

Variability of Disk Emission in Pre-main-sequence and Related Stars. IV. Investigating the Structural Changes in the Inner Disk Region of MWC 480

Abstract We present five epochs of near-IR observations of the protoplanetary disk around MWC 480 (HD 31648) obtained with the SpeX spectrograph on NASA’s Infrared Telescope Facility between 2007 and 2013, inclusive. Using the measured line fluxes in the Pa β and Br γ lines, we found the mass accretion rates to be (1.26–2.30) × 10−7 M ⊙ yr−1 and (1.4–2.01) × 10−7 M ⊙ yr−1, respectively, but which varied by more than 50% from epoch to epoch. The spectral energy distribution reveals a variability of about 30% between 1.5 and 10 μm during this same period of time. We investigated the variability using of the continuum emission of the disk in using the Monte-Carlo Radiative Transfer Code HOCHUNK3D. We find that varying the height of the inner rim successfully produces a change in the NIR flux but lowers the far-IR emission to levels below all measured fluxes. Because the star exhibits bipolar flows, we utilized a structure that simulates an inner disk wind to model the variability in the near-IR, without producing flux levels in the far-IR that are inconsistent with existing data. For this object, variable near-IR emission due to such an outflow is more consistent with the data than changing the scale height of the inner rim of the disk.

Polarization simulations of stellar wind bow-shock nebulae – I. The case of electron scattering

Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow shocks are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionisation state of the material within the structure. In this paper we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow shock-shaped region of enhanced density surrounding a stellar source. We present results assuming electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes we consider resolved and unresolved cases. We discuss the implication of these results as well as their possible use along with observational data to constrain the properties of observed bow shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the scattering region, and the relative luminosities of the star and shock.

Radiative transfer calculations of the diffuse ionized gas in disc galaxies with cosmic ray feedback

The large vertical scale heights of the diffuse ionised gas (DIG) in disc galaxies are challenging to model, as hydrodynamical models including only thermal feedback seem to be unable to support gas at these heights. In this paper, we use a three dimensional Monte Carlo radiation transfer code to post-process disc simulations of the Simulating the Life-Cycle of Molecular Clouds (SILCC) project that include feedback by cosmic rays. We show that the more extended discs in simulations including cosmic ray feedback naturally lead to larger scale heights for the DIG which are more in line with observed scale heights. We also show that including a fiducial cosmic ray heating term in our model can help to increase the temperature as a function of disc scale height, but fails to reproduce observed DIG nitrogen and sulphur forbidden line intensities. We show that, to reproduce these line emissions, we require a heating mechanism that affects gas over a larger density range than is achieved by cosmic ray heating, which can be achieved by fine tuning the total luminosity of ionising sources to get an appropriate ionising spectrum as a function of scale height. This result sheds a new light on the relation between forbidden line emissions and temperature profiles for realistic DIG gas distributions.

A Monte Carlo Radiation Transfer Study of Photospheric Emission in Gamma-Ray Bursts

Abstract We present the analysis of photospheric emission for a set of hydrodynamic simulations of long duration gamma-ray burst jets from massive compact stars. The results are obtained by using the Monte Carlo Radiation Transfer code (MCRaT) to simulate thermal photons scattering through the collimated outflows. MCRaT allows us to study explicitly the time evolution of the photosphere within the photospheric region, as well as the gradual decoupling of the photon and matter counterparts of the jet. The results of the radiation transfer simulations are also used to construct light curves and time-resolved spectra at various viewing angles, which are then used to make comparisons with observed data and outline the agreement and strain points between the photospheric model and long duration gamma-ray burst observations. We find that our fitted time-resolved spectral Band β parameters are in agreement with observations, even though we do not consider the effects of nonthermal particles. Finally, the results are found to be consistent with the Yonetoku correlation, but bear some strain with the Amati correlation.

Synthetic observations of protostellar multiple systems

Observations of protostars are often compared with synthetic observations of models in order to infer the underlying physical properties of the protostars. The majority of these models have a single protostar, attended by a disc and an envelope. However, observational and numerical evidence suggests that a large fraction of protostars form as multiple systems. This means that fitting models of single protostars to observations may be inappropriate. We produce synthetic observations of protostellar multiple systems undergoing realistic, non-continuous accretion. These systems consist of multiple protostars with episodic luminosities, embedded self-consistently in discs and envelopes. We model the gas dynamics of these systems using smoothed particle hydrodynamics and we generate synthetic observations by post-processing the snapshots using the spamcart Monte Carlo radiative transfer code. We present simulation results of three model protostellar multiple systems. For each of these, we generate 4 × 104 synthetic spectra at different points in time and from different viewing angles. We propose a Bayesian method, using similar calculations to those presented here, but in greater numbers, to infer the physical properties of protostellar multiple systems from observations.