Monte Carlo Radiative Transfer Research Articles

On the evolution of the observed mass-to-length relationship for star-forming filaments

ABSTRACT The interstellar medium is threaded by a hierarchy of filaments from large scales (∼100 pc) to small scales (∼0.1 pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the ‘Cloud Factory’ galactic-scale ISM suite in combination with 3D Monte Carlo radiative transfer code polaris to investigate how filamentary structure evolves over time. We produce synthetic dust continuum observations in three regions with a series of snapshots and use the filfinder algorithm to identify filaments in the dust derived column density maps. When the synthetic filaments mass and length are plotted on an mass–length (M–L) plot, we see a scaling relation of L ∝ M0.45 similar to that seen in observations, and find that the filaments are thermally supercritical. Projection effects systematically affect the masses and lengths measured for the filaments, and are particularly severe in crowded regions. In the filament M–L diagram we identify three main evolutionary mechanisms: accretion, segmentation, and dispersal. In particular we find that the filaments typically evolve from smaller to larger masses in the observational M–L plane, indicating the dominant role of accretion in filament evolution. Moreover, we find a potential correlation between line mass and filament growth rate. Once filaments are actively star forming they then segment into smaller sections, or are dispersed by internal or external forces.

Improving Monte Carlo radiative transfer in the regime of high optical depths: The minimum scattering order

Radiative transfer simulations are a powerful tool that enables the calculation of synthetic images of a wide range of astrophysical objects. These simulations are often based on the Monte Carlo method, as it provides the needed versatility that allows the consideration of the diverse and often complex conditions found in those objects. However, this method faces fundamental problems in the regime of high optical depths which may result in overly noisy images and severely underestimated flux values. In this study, we propose an advanced Monte Carlo radiative transfer method, namely, an enforced minimum scattering order that is aimed at providing a minimum quality of determined flux estimates. For that purpose, we extended our investigations of the scattering order problem and derived an analytic expression for the minimum number of interactions that depends on the albedo and optical depth of the system, which needs to be considered during a simulation to achieve a certain coverage of the scattering order distribution. The method is based on the utilization of this estimated minimum scattering order and enforces the consideration of a sufficient number of interactions during a Monte Carlo radiative transfer simulation. Moreover, we identified two notably distinct cases that shape the kind of complexity that arises in such simulations: the albedo-dominated case and the optical depth-dominated case. Based on that, we analyzed the implications related to the best usage of a stretching method as a means to alleviate the scattering order problem. We find that its most suitable application requires taking into account the value of the albedo as well as the optical depth. Furthermore, we argue that the derived minimum scattering order can be used to assess the performance of a stretching method with regard to the scattering orders its usage promotes. Finally, we stress the need for developing advanced pathfinding techniques to fully solve the problem of Monte Carlo radiative transfer simulations in the regime of high optical depths.

Changing streetlighting schemes and the ecological availability of darkness.

Artificial light at night (ALAN), from streetlights and other sources, has a wide variety of impacts on the natural environment. A significant challenge remains, however, to predict at intermediate spatial extents (e.g. across a city) the ALAN that organisms experience under different lighting regimes. Here we use Monte Carlo radiative Transfer to model the three-dimensional lighting environment at, and just above, ground level, on the spatial scales at which animals and humans experience it. We show how this technique can be used to model a suite of both real and hypothetical lighting environments, mimicking the transition of public infrastructure between different lighting technologies. We then demonstrate how the behaviour of animals experiencing these simulated lighting environments can be emulated to probe the availability of darkness, and dark corridors, within them. Our simulations show that no single lighting technology provides an unmitigated alleviation of negative impacts within urban environments, and that holistic treatments of entire lighting environments should be employed when understanding how animals use and traverse them.

Numerical Simulation of Photospheric Emission in Long Gamma-Ray Bursts: Prompt Correlations, Spectral Shapes, and Polarizations

We explore the properties of photospheric emission in the context of long gamma-ray bursts (GRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the observed Yonetoku, Amati, and Golenetskii correlations. It is also shown that the spectral peak energy and luminosity correlate with the bulk Lorentz factor, as indicated in the literature. On the other hand, synthetic spectral shapes tend to be narrower than those of the observations. This result indicates that an additional physical process that can provide nonthermal broadening is needed to reproduce the spectral features. Furthermore, the polarization analysis finds that, while the degree of polarization is low for the emission from the jet core (Π < 4%), it tends to increase with viewing angle outside of the core and can be as high as Π ∼ 20%–40% in an extreme case. This suggests that the typical GRBs show systematically low polarization compared to softer, dimmer counterparts (X-ray-rich GRBs and X-ray flashes). Interestingly, our simulations indicate that photospheric emission exhibits large temporal variation in the polarization position angle (Δψ ∼ 90°), which may be compatible with those inferred in observations. A notable energy dependence of the polarization property is another characteristic feature found in the current study. Particularly, the difference in the position angle among different energy bands can be as large as ∼90°.

A Morphokinematic Study of the Enigmatic Emission Nebula NGC 6164/5 Surrounding the Magnetic O-type Star HD 148937

HD 148937 is a peculiar massive star (Of?p) with a strong magnetic field (1 kG). The hourglass-shaped emission nebula NGC 6164/5 surrounds this star. This nebula is presumed to originate from episodic mass-loss events of the central O-type star, but the detailed formation mechanism is not yet well understood. Grasping its three-dimensional structure is essential to uncovering the origin of this nebula. Here we report the high-resolution multiobject spectroscopic observations of NGC 6164/5 using the GIRAFFE on the 8.2 m Very Large Telescope. Integrated intensity maps constructed from several spectral lines delineate well the overall shape of this nebula, such as the two bright lobes and the inner gas region. The position–velocity diagrams show that the two bright lobes are found to be redshifted and blueshifted, respectively, while the inner region has multiple layers. We consider a geometric model composed of a bilateral outflow harboring nitrogen-enriched knots and expanding inner shells. Its spectral features are then simulated by using a Monte Carlo radiative transfer technique for different sets of velocities. Some position–velocity diagrams from simulations are very similar to the observed ones. According to the model that best reproduces the observational data, the two bright lobes and the nitrogen-enriched knots are moving away from HD 148937 at about 120 km s−1. Their minimum kinematic age is estimated to be about 7500 yr. We discuss possible formation mechanisms of this nebula in the context of binary interaction.

Diagnosis of Circumstellar Matter Structure in Interaction-powered Supernovae with Hydrogen Line Features

Some supernovae (SNe) are powered by the collision of the SN ejecta with dense circumstellar matter (CSM). Their emission spectra show characteristic line shapes of combined broad emission and narrow P Cygni lines, which should closely relate to the CSM structure and the mass-loss mechanism that creates the dense CSM. We quantitatively investigate the relationship between the line shape and the CSM structure by Monte Carlo radiative transfer simulations, considering two representative cases of dense CSM formed by steady and eruptive mass loss. Comparing the Hα emission between the two cases, we find that a narrow P Cygni line appears in the eruptive case but does not appear in the steady case due to the difference in the velocity gradient in the dense CSM. We also reproduce the blueshifted photon excess observed in some Type IIn SNe, which is formed by photon transport across the shock wave, and find the relationship between the velocity of the shocked matter and the amount of blueshift of the photon excess. We conclude that the presence or absence of narrow P Cygni lines can distinguish the mass-loss mechanism and suggest high-resolution spectroscopic observations with λ/Δλ ≳ 104 after the light-curve peak for applying this diagnostic method.

Improving Monte Carlo radiative transfer simulations: A shift of framework

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 observations. Consequently, a broad variety of codes has been implemented and optimized with the goal of solving this task efficiently. To that end, two distinct frameworks have emerged, namely, the extinction and the scattering framework, which form the basis of the path determination procedures of those codes. These procedures affect the step length of simulated photon packages and are used for determining flux estimates. Despite the fact that these simulations play an important role at present and thus require significant computational resources, little attention has been paid to the benefits and the drawbacks of both frameworks so far. In this study, we investigate their differences and assess their performance with regard to the quality of thereby obtained flux estimates, with a particular focus on the required computational demand. To that end, we use a testbed composed of an infinite plane-parallel slab, illuminated from one side, and we determine transmitted intensity using MCRT simulations for both frameworks. We find that there are vast differences between the frameworks with regard to their convergence speed. The scattering framework outperforms the extinction framework across all considered optical depths and albedos when solving this task, particularly in the regime of high optical depths. Its implementation can therefore greatly benefit all modern MCRT codes as it has the potential to significantly reduce required computation times. Thus, we highly recommend its consideration for various tasks that require MCRT simulations.

Shaping the CO snowline in protoplanetary disks

Context. Characterizing the dust thermal structure in protoplanetary disks is a fundamental task because the dust surface temperature can affect both the planetary formation and the chemical evolution. Because the temperature depends on many parameters, including the grain size, it can be challenging to properly model the grain temperature structure. Many chemistry disk models usually employ a sophisticated single dust structure designed to reproduce the effect of a realistic population presumably composed of a large diversity of sizes. This generally represents a good approximation in most cases. Nonetheless, it dilutes the effects of the complex radiative interactions between the different grain populations on the resulting dust temperature, and thus, the chemistry. Aims. We seek to show that the radiative interactions between dust grains of different sizes can induce a nontrivial dust temperature structure that cannot be reproduced by a single dust population and that can significantly affect the chemical outcome. Methods. The disk thermal structures were computed using the Monte Carlo radiative transfer code RADMC-3D. The thermal structures were postprocessed using the gas-grain code NAUTILUS to calculate the evolution of the chemical abundance. Results. We find that simultaneously using at least two independent dust grain populations in disk models produces a complex temperature structure due to the starlight that is intercepted by the upper layers of the disk. In particular, we find that micron-sized dust grains are warmer than larger grains and can even show a radial temperature bump in some conditions. This dust temperature spread between the grain populations results in the segregation of the CO snowline and in an unexpected CO gas hole in the midplane. We compare the results with observed close to edge-on class I/II disks. Conclusions. Our study shows that the size dependence of the dust temperature significantly impacts the chemistry, and that two dust populations at least are required to account for this property of the thermal structure in protoplanetary disk models over a wide range of disk masses and dust properties.

Disc Tearing in a Be Star: Predicted 3D Observations

Abstract We build on our previous work involving smoothed particle hydrodynamic simulations of Be stars, by using the model that exhibited disc tearing as input into the three-dimensional Monte Carlo radiative transfer code hdust to predict observables from a variety of viewing angles throughout the disc tearing process. We run one simulation at the start of each orbital period from 20 to 72 orbital periods, which covers two complete disc tearing events. The resulting trends in observables are found to be dependent on the relative position of the observer and the tearing disc. The H$\rm \alpha$ equivalent width, V magnitude, and polarization can all increase or decrease in any combination depending on the viewpoint of the observer. The H$\rm \alpha$ line profile also displays changes in strength and peak separation throughout the tearing process. We show how the outer disc of the torn system can have a large effect on the H$\rm \alpha$ line profile, and also contributes to a wavelength-dependent polarization position angle, resulting in a similar sawtooth shape to the polarization percentage. Finally, we compare our predictions to Pleione (28 Tau) where evidence has suggested that a disc tearing event has occurred in the past. We find that our tearing disc model can broadly match the trends seen in Pleione’s observables, as well as produce the two-component H$\rm \alpha$ lines observed in Pleione. This is the strongest evidence, thus far, of Pleione’s disc having indeed experienced a tearing event.

Monte Carlo radiative transfer peel off mechanism for spatially extended detectors

We present an extension to the well-known peel off optimization for Monte Carlo radiative transfer simulations. The classical method is only applicable when the distance between the detector and the peel off event is much bigger than the size of the detector. We use two alternatives to the classical method and calculate the peel off intensity via a subdivision and an integration method. We compare their performance for a realistic scenario and derive guidelines for a general treatment. This allows for precise peel off calculations at any distance to the detector surface.

Synthetic light curves and spectra from a self-consistent 2D simulation of an ultra-strippped supernova

ABSTRACT Spectroscopy is an important tool for providing insights into the structure of core-collapse supernova explosions. We use the Monte Carlo radiative transfer code artis to compute synthetic spectra and light curves based on a two-dimensional explosion model of an ultra-stripped supernova. These calculations are designed both to identify observable fingerprints of ultra-stripped supernovae and as a proof of principle for using synthetic spectroscopy to constrain the nature of stripped-envelope supernovae more broadly. We predict characteristic spectral and photometric features for our ultra-stripped explosion model, and find that these do not match observed ultra-stripped supernova candidates like SN 2005ek. With a peak bolometric luminosity of $6.8\times 10^{41}\, \mathrm{erg}\, \mathrm{s}^{-1}$, a peak magnitude of $-15.9\, \mathrm{mag}$ in R band, and Δm15,R = 3.50, the model is even fainter and evolves even faster than SN 2005ek as the closest possible analogue in photometric properties. The predicted spectra are extremely unusual. The most prominent features are Mg ii lines at $2 {,}800\, {\mathring{\rm A}}$ and $4 {,}500\, {\mathring{\rm A}}$ and the infrared Ca triplet at late times. The Mg lines are sensitive to the multidimensional structure of the model and are viewing-angle dependent. They disappear due to line blanketing by iron group elements in a spherically averaged model with additional microscopic mixing. In future studies, multi-D radiative transfer calculations need to be applied to a broader range of models to elucidate the nature of observed Type Ib/c supernovae.

A disc wind model for blueshifts in quasar broad emission lines

ABSTRACT Blueshifts – or, more accurately, blue asymmetries – in broad emission lines such as C iv λ1550 are common in luminous quasars and correlate with fundamental properties such as Eddington ratio and broad absorption line (BAL) characteristics. However, the formation of these blueshifts is still not understood, and neither is their physical connection to the BAL phenomenon or accretion disc. In this work, we present Monte Carlo radiative transfer and photoionization simulations using parametrized biconical disc-wind models. We take advantage of the azimuthal symmetry of a quasar and show that we can reproduce C iv blueshifts provided that (i) the disc-mid-plane is optically thick out to radii beyond the line formation region, so that the receding wind bicone is obscured; and (ii) the system is viewed from relatively low (that is, more face-on) inclinations (≲40°). We show that C iv emission-line blueshifts and BALs can form in the same wind structure. The velocity profile of the wind has a significant impact on the location of the line formation region and the resulting line profile, suggesting that the shape of the emission lines can be used as a probe of wind-driving physics. While we are successful at producing blueshifts/blue asymmetries in outflows, we struggle to match the detailed shape or skew of the observed emission-line profiles. In addition, our models produce redshifted emission-line asymmetries for certain viewing angles. We discuss our work in the context of the C iv λ1550 emission blueshift versus equivalent-width space and explore the implications for quasar disc wind physics.

Modelling variable linear polarization produced by Co-rotating Interaction Regions (CIRs) across optical recombination lines of Wolf–Rayet stars

ABSTRACT Massive star winds are structured both stochastically (‘clumps’) and often coherently (Co-rotation Interaction Regions, or CIRs). Evidence for CIRs threading the winds of Wolf–Rayet (WR) stars arises from multiple diagnostics including linear polarimetry. Some observations indicate changes in polarization position angle across optical recombination emission lines from a WR star wind but limited to blueshifted Doppler velocities. We explore a model involving a spherical wind with a single conical CIR stemming from a rotating star as qualitative proof-of-concept. To obtain a realistic distribution of limb polarization and limb darkening across the pseudo-photosphere formed in the optically thick wind of a WR star, we used Monte Carlo radiative transfer (MCRT). Results are shown for a parameter study. For line properties similar to WR 6 (EZ CMa; HD 50896), combining the MCRT results, a simple model for the CIR, and the Sobolev approximation for the line formation, we were able to reproduce variations in both polarization amplitude and position angle commensurate with observations. Characterizing CIRs in WR winds has added importance for providing stellar rotation periods since the vsin i values are unobtainable because the pseudo-photosphere forms in the wind itself.

Non-LTE Monte Carlo radiative transfer – III. The thermal properties of tilted and warped Be star discs

ABSTRACT We use the three-dimensional Monte Carlo radiative transfer code hdust to model Be stars where the disc is tilted from the equatorial plane of the star. We compute 128 models across four spectral types, B0, B2, B5, and B8, tilting the disc by 0, 10○, 20○, and 40○, respectively, while varying disc density according to spectral type. We also compute every model for an average and high stellar rotation rate. We first discuss non-tilted disc temperatures and show its nonlinear dependence on stellar and disc parameters. We find that tilting the disc minimally affects the density-weighted average disc temperature, but tilting does create a temperature asymmetry in disc cross-sections, which is more pronounced for a faster rotation rate. We also investigate the effect tilting has on V-band magnitude, polarization, and the H$\rm \alpha$ line. Tilting the disc does affect these observables, but the changes are entirely dependent on the position of the observer relative to the direction of tilt. We find the observables that distinguish tilting from a change in density or geometry are the H$\rm \alpha$ line shapes, where it can transition between single-peaked and double-peaked, and the polarization position angle, whose value is dependent on the projected major elongation axis of the disc on the sky. We also present one early- and one late-type model with warped discs. We find their temperature structure varies a small amount from the uniformly tilted models, and the different observables correspond to different tilt angles, consistent with their expected volume of origin within the disc.

Emulating radiative transfer with artificial neural networks

ABSTRACT Forward-modeling observables from galaxy simulations enables direct comparisons between theory and observations. To generate synthetic spectral energy distributions (SEDs) that include dust absorption, re-emission, and scattering, Monte Carlo radiative transfer is often used in post-processing on a galaxy-by-galaxy basis. However, this is computationally expensive, especially if one wants to make predictions for suites of many cosmological simulations. To alleviate this computational burden, we have developed a radiative transfer emulator using an artificial neural network (ANN), ANNgelina, that can reliably predict SEDs of simulated galaxies using a small number of integrated properties of the simulated galaxies: star formation rate, stellar and dust masses, and mass-weighted metallicities of all star particles and of only star particles with age &amp;lt;10 Myr. Here, we present the methodology and quantify the accuracy of the predictions. We train the ANN on SEDs computed for galaxies from the IllustrisTNG project’s TNG50 cosmological magnetohydrodynamical simulation. ANNgelina is able to predict the SEDs of TNG50 galaxies in the ultraviolet (UV) to millimetre regime with a typical median absolute error of ∼7 per cent. The prediction error is the greatest in the UV, possibly due to the viewing-angle dependence being greatest in this wavelength regime. Our results demonstrate that our ANN-based emulator is a promising computationally inexpensive alternative for forward-modeling galaxy SEDs from cosmological simulations.

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydrodynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts the electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affects the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high-energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Impact of hot exozodiacal dust on the polarimetric analysis of close-in exoplanets

Context. Hot exozodiacal dust (HEZD) found around main-sequence stars through interferometric observations in the photometric bands H to L is located close to the dust sublimation radius, potentially at orbital radii comparable to those of close-in exoplanets. Consequently, HEZD has a potential influence on the analysis of the scattered-light polarization of close-in exoplanets and vice versa. Aims. We analyze the impact of HEZD on the polarimetric characterization of close-in exoplanets. This study is motivated in particular by the recently proven feasibility of exoplanet polarimetry. Methods. Applying the 3D Monte Carlo radiative transfer code POLARIS in an extended and optimized version for radiative transfer in exoplanetary atmospheres and an analytical tool for modeling the HEZD, we simulated and compared the polarization characteristics of the wavelength-dependent scattered-light polarization of HEZD and close-in exoplanets. As a starting point for our analysis, we defined a reference model consisting of a close-in exoplanet with a scattered-light polarization consistent with the upper limit determined for WASP-18b, and a HEZD consistent with the near-infrared excess detected for HD 22484 (10 Tau). Results. The varied parameters are the planetary phase angle (0°–180°), the dust grain radius (0.02 µm−10 µm), the HEZD mass (10−10 M⊙−10−8 M⊙), the orbital inclination (0°−90°), the composition of the planetary atmosphere (Mie and Rayleigh scattering atmosphere), the orbital radius of HEZD (0.02 au−0.4 au), and the planetary orbital radius (0.01 au−0.05 au). The dust grain radius has the strongest influence on the polarimetric analysis due to its significant impact on the wavelength-dependent polarization characteristics and the total order of magnitude of the scattered-light polarization. In certain scenarios, the scattered-light polarization of the HEZD even exceeds that of the close-in exoplanet, for example for a dust grain radius of 0.1 µm, a HEZD mass of 8 × 10−10 M⊙, an orbital radius of HEZD of 0.04 au and an orbital inclination of 90°. Conclusions. The presence of HEZD potentially has a significant impact on the polarimetric investigations of close-in exoplanets. Furthermore, interferometric observations are required to better constrain the parameter space for HEZD and thus the possible resulting scattered-light polarization.

Systematic Investigation of Very-early-phase Spectra of Type Ia Supernovae

It has been widely accepted that Type Ia supernovae (SNe Ia) are thermonuclear explosions of a CO white dwarf. However, the natures of the progenitor system(s) and explosion mechanism(s) are still unclarified. Thanks to the recent development of transient observations, they are now frequently discovered shortly after the explosion, followed by rapid spectroscopic observations. In this study, by modeling very-early-phase spectra of SNe Ia, we try to constrain the explosion models of SNe Ia. By using the Monte Carlo radiation transfer code, TARDIS, we estimate the properties of their outermost ejecta. We find that the photospheric velocity of normal-velocity supernovae (NV SNe) in the first week is ∼15,000 km s−1. The outer velocity, to which the carbon burning extends, spans the range between ∼20,000 and 25,000 km s−1. The ejecta density of NV SNe also shows a large diversity. For high-velocity supernovae (HV SNe) and 1999aa-like SNe, the photospheric velocity is higher, ∼20,000 km s−1. They have different photospheric densities, with HV SNe having higher densities than 1999aa-like SNe. For all these types, we show that the outermost composition is closely related to the outermost ejecta density; the carbon-burning layer and the unburnt carbon layer are found in the higher-density and lower-density objects, respectively. This finding suggests that there might be two sequences, the high-density and carbon-poor group (HV SNe and some NV SNe) and the low-density and carbon-rich group (1999aa-like and other NV SNe), which may be associated with different progenitor channels.

Exploring the polarization of axially symmetric supernovae with unsupervised deep learning

ABSTRACT The measurement of non-zero polarization can be used to infer the presence of departures from spherical symmetry in supernovae (SNe). The origin of the majority of the intrinsic polarization observed in SNe is in electron scattering, which induces a wavelength-independent continuum polarization that is generally observed to be low ($\lesssim\!\! 1{{\ \rm per\, cent}}$) for all SN types. The key indicator of asymmetry in SNe is the polarization observed across spectral lines, in particular the characteristic ‘inverse P Cygni’ profile. The results of a suite of 900 Monte Carlo radiative transfer simulations are presented here. These simulations cover a range of possible axisymmetric structures (including unipolar, bipolar, and equatorial enhancements) for the line-forming region of the Ca ii infrared triplet. Using a variational auto-encoder, seven key latent parameters are learned that describe the relationship between Stokes I and q, under the assumption of an axially symmetric line-forming region and resonant scattering. Likelihood-free inference techniques are used to invert the Stokes I and q line profiles, in the latent space, to derive the underlying geometries. For axially symmetric structures that yield an observable ‘dominant axis’ on the Stokes q–u plane, we propose the existence of a geometry ‘conjugate’ (which is indistinguishable under a rotation of π/2). Using this machine learning infrastructure, we attempt to identify possible geometries associated with spectropolarimetric observations of the Type Ib SN 2017gax.

A Hydrodynamic Study of the Escape of Metal Species and Excited Hydrogen from the Atmosphere of the Hot Jupiter WASP-121b

In the near-UV and optical transmission spectrum of the hot Jupiter WASP-121b, recent observations have detected strong absorption features of Mg, Fe, Ca, and Hα, extending outside of the planet’s Roche lobe. Studying these atomic signatures can directly trace the escaping atmosphere and constrain the energy balance of the upper atmosphere. To understand these features, we introduce a detailed forward model by expanding the capability of a one-dimensional model of the upper atmosphere and hydrodynamic escape to include important processes of atomic metal species. The hydrodynamic model is coupled to a Lyα Monte Carlo radiative transfer calculation to simulate the excited hydrogen population and associated heating/ionization effects. Using this model, we interpret the detected atomic features in the transmission spectrum of WASP-121b and explore the impact of metals and excited hydrogen on its upper atmosphere. We demonstrate the use of multiple absorption lines to impose stronger constraints on the properties of the upper atmosphere than the analysis of a single transmission feature can provide. In addition, the model shows that line broadening due to atmospheric outflow driven by Roche lobe overflow is necessary to explain the observed line widths and highlights the importance of the high mass-loss rate caused by Roche lobe overflow, which requires careful consideration of the structure of the lower and middle atmosphere. We also show that metal species and excited-state hydrogen can play an important role in the thermal and ionization balance of ultrahot Jupiter thermospheres.