There are now hundreds of publicly available supernova spectral time series. Radiative transfer modeling of this data gives insights into the physical properties of these explosions such as the composition, the density structure, or the intrinsic luminosity---this is invaluable for understanding the supernova progenitors, the explosion mechanism, or for constraining the supernova distance. However, a detailed parameter study of the available data has been out of reach due to the high dimensionality of the problem coupled with the still significant computational expense. We tackle this issue through the use of machine-learning emulators, which are algorithms for high-dimensional interpolation. These use a pre-calculated training dataset to mimic the output of a complex code but with run times orders of magnitude shorter. We present the application of such an emulator to synthetic type II supernova spectra generated with the TARDIS radiative transfer code. The results show that with a relatively small training set of 780 spectra we can generate emulated spectra with interpolation uncertainties of less than one percent. We demonstrate the utility of this method by automatic spectral fitting of two well-known type IIP supernovae; as an exemplary application, we determine the supernova distances from the spectral fits using the tailored-expanding-photosphere method. We compare our results to previous studies and find good agreement. This suggests that emulation of TARDIS spectra can likely be used to perform automatic and detailed analysis of many transient classes putting the analysis of large data repositories within reach.

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion.

The theory and numerical modelling of radiation processes and radiative transfer play a key role in astrophysics: they provide the link between the physical properties of an object and the radiation it emits. In the modern era of increasingly high-quality observational data and sophisticated physical theories, development and exploitation of a variety of approaches to the modelling of radiative transfer is needed. In this article, we focus on one remarkably versatile approach: Monte Carlo Radiative Transfer (MCRT). We describe the principles behind this approach, and highlight the relative ease with which they can (and have) been implemented for application to a range of astrophysical problems. All MCRT methods have in common a need to consider the adverse consequences of Monte Carlo noise in simulation results. We overview a range of methods used to suppress this noise and comment on their relative merits for a variety of applications. We conclude with a brief review of specific applications for which MCRT methods are currently popular and comment on the prospects for future developments.

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 two previously published sets of dilution factors by Eastman et al. (1996) and by Dessart & Hillier (2005), and discuss the potential sources of the discrepancies between studies.

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.

Upcoming high-cadence transient survey programmes will produce a wealth of observational data for Type Ia supernovae. These data sets will contain numerous events detected very early in their evolution, shortly after explosion. Here, we present synthetic light curves, calculated with the radiation hydrodynamical approach Stella for a number of different explosion models, specifically focusing on these first few days after explosion. We show that overall the early light curve evolution is similar for most of the investigated models. Characteristic imprints are induced by radioactive material located close to the surface. However, these are very similar to the signatures expected from ejecta-CSM or ejecta-companion interaction. Apart from the pure deflagration explosion models, none of our synthetic light curves exhibit the commonly assumed power-law rise. We demonstrate that this can lead to substantial errors in the determination of the time of explosion. In summary, we illustrate with our calculations that even with very early data an identification of specific explosion scenarios is challenging, if only photometric observations are available.

We investigate the consequences of fairly normal Type Ia supernovae being embedded in compact and dense envelopes of carbon and oxygen rich circumstellar material by means of detailed radiation hydrodynamic simulations. Our main focus rests on exploring the effects of the interaction between ejecta and circumstellar material on the ejecta evolution and the broad-band light curve. In our calculations, we find that a strong reverse shock efficiently decelerates and compresses the ejecta material. This leads to a significant broadening of the optical light curve, a longer rise to maximum and a slower decline in the tail phase. During the interaction, substantial radiative energy is generated, which mostly emerges in the extreme ultraviolet and X-ray regime. Only if reprocessing due to radiation-matter interactions is very efficient, a significant boost in the optical light curve is observed. We discuss these findings in particular in the context of the super-luminous event SN 2009dc. As our calculations are able to reproduce a number of its peculiar properties, we conclude that the flavour of the interaction scenario investigated in this work constitutes a promising candidate to explain such 'Super-Chandrasekhar' supernovae.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

We report the results of our follow-up campaign of the supernova impostor PSN J09132750+7627410, based on optical data covering $\sim250\,\rm{d}$. From the beginning, the transient shows prominent narrow Balmer lines with P-Cygni profiles, with a blue-shifted absorption component becoming more prominent with time. Along the $\sim3\,\rm{months}$ of the spectroscopic monitoring, broad components are never detected in the hydrogen lines, suggesting that these features are produced in slowly expanding material. The transient reaches an absolute magnitude $M_r=-13.60\pm0.19\,\rm{mag}$ at maximum, a typical luminosity for supernova impostors. Amateur astronomers provided $\sim4\,\rm{years}$ of archival observations of the host galaxy, NGC 2748. The detection of the quiescent progenitor star in archival images obtained with the Hubble Space Telescope suggests it to be an $18-20$\msun white-yellow supergiant.

We present spectroscopic and photometric observations for the Type Ibn supernova (SN) dubbed PSN J07285387+3349106. Using data provided by amateur astronomers, we monitored the photometric rise of the SN to maximum light, occurred on 2015 February 18.8 UT (JD(max,V) = 2457072.0 +- 0.8). PSN J07285387+3349106 exploded in the inner region of an infrared luminous galaxy, and is the most reddened SN Ibn discovered so far. We apply multiple methods to derive the total reddening to the SN, and determine a total colour excess E(B-V)(tot) = 0.99 +- 0.48 mag. Accounting for the reddening correction, which is affected by a large uncertainty, we estimate a peak absolute magnitude of M(V) = -20.30 +- 1.50. The spectra are dominated by continuum emission at early phases, and He I lines with narrow P-Cygni profiles are detected. We also identify weak Fe III and N II features. All these lines show an absorption component which is blue-shifted by about 900-1000 km/s. The spectra also show relatively broad He I line wings with low contrast, which extend to above 3000 km/s. From about 2 weeks past maximum, broad lines of O I, Mg II and the Ca II near-infrared triplet are identified. The composition and the expansion velocity of the circumstellar material, and the presence of He I and alpha-elements in the SN ejecta indicate that PSN J07285387+3349106 was produced by the core-collapse of a stripped-envelope star. We suggest that the precursor was WNE-type Wolf-Rayet star in its dense, He-rich circumstellar cocoon.

We present the results of the spectroscopic and photometric monitoring campaign of ASASSN-15ed. The transient was discovered quite young by the All Sky Automated Survey for SuperNovae (ASAS-SN). Amateur astronomers allowed us to sample the photometric SN evolution around maximum light, which we estimate to have occurred on JD = 2457087.4 +- 0.6 in the r-band. Its apparent r-band magnitude at maximum was r = 16.91 +- 0.10, providing an absolute magnitude M(r) ~ -20.04 +- 0.20, which is slightly more luminous than the typical magnitudes estimated for Type Ibn SNe. The post-peak evolution was well monitored, and the decline rate (being in most bands around 0.1 mag/d during the first 25 d after maximum) is marginally slower than the average decline rates of SNe Ibn during the same time interval. The object was initially classified as a Type Ibn SN because early-time spectra were characterized by a blue continuum with superimposed narrow P-Cygni lines of He I, suggesting the presence of a slowly moving (1200-1500 km/s), He-rich circumstellar medium. Later on, broad P-Cygni He I lines became prominent. The inferred velocities, as measured from the minimum of the broad absorption components, were between 6000 and 7000 km/s. As we attribute these broad features to the SN ejecta, this is the first time we have observed the transition of a Type Ibn SN to a Type Ib SN.

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite-volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multi-line effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach is demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally we discuss some of the current limitations of our method, some possible extensions and potential future applications.

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions.

The UV spectra of nova-like variables are dominated by emission from the accretion disk, modified by scattering in a wind emanating from the disk. Here we model the spectra of RW Tri and UX UMa, the only two eclipsing nova-likes which have been observed with the Hubble Space Telescope in the far-ultraviolet, in an attempt to constrain the geometry and the ionization structure of their winds. Using our Monte Carlo radiative transfer code we computed spectra for simply-parameterized axisymmetric biconical outflow models and were able to find plausible models for both systems. These reproduce the primary UV resonance lines - N V, Si IV, and C IV - in the observed spectra in and out of eclipse. The distribution of these ions in the wind models is similar in both cases as is the extent of the primary scattering regions in which these lines are formed. The inferred mass loss rates are 6% to 8% of the mass accretion rates for the systems. We discuss the implication of our point models for our understanding of accretion disk winds in cataclysmic variables.