Monte-Carlo methods for NLTE spectral synthesis of supernovae

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.Comment: Accepted for publication by Astronomy & Astrophysic

Late-time spectral line formation in Type IIb supernovae, with application to SN 1993J, SN 2008ax, and SN 2011dh

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modeling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date - SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modeling of the [O I] 6300, 6364 lines constrains the progenitors of these three SNe to the M_ZAMS=12-16 M_sun range (ejected oxygen masses 0.3-0.9 M_sun), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from M_ZAMS >= 17 M_sun progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low/moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of of 0.02-0.14 M_sun is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO-burning gives strong [N II] 6548, 6583 emission lines that dominate over H-alpha emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable H-alpha emission or absorption after ~150 days, and nebular phase emission seen around 6550 A is in many cases likely caused by [N II] 6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated...(abridged)Comment: Published versio

A metallicity study of 1987A-like supernova host galaxies

The origin of the blue supergiant (BSG) progenitor of Supernova (SN) 1987A has long been debated, along with the role that its sub-solar metallicity played. We now have a sample of 1987A-like SNe that arise from the core collapse (CC) of BSGs. The metallicity of the explosion sites of the known BSG SNe is investigated, as well as their association to star-forming regions. Both indirect and direct metallicity measurements of 13 BSG SN host galaxies are presented, and compared to those of other CC SN types. Indirect measurements are based on the known luminosity-metallicity relation and on published metallicity gradients of spiral galaxies. To provide direct estimates based on strong line diagnostics, we obtained spectra of each BSG SN host both at the SN explosion site and at the positions of other HII regions. Continuum-subtracted Ha images allowed us to quantify the association between BSG SNe and star-forming regions. BSG SNe explode either in low-luminosity galaxies or at large distances from the nuclei of luminous hosts. Therefore, their indirectly measured metallicities are typically lower than those of SNe IIP and Ibc. This is confirmed by the direct estimates, which show slightly sub-solar values (12+log(O/H)=8.3-8.4 dex), similar to that of the Large Magellanic Cloud (LMC), where SN 1987A exploded. However, two SNe (1998A and 2004em) were found at near solar metallicity. SNe IIb have a metallicity distribution similar to that of BSG SNe. Finally, the association to star-forming regions is similar among BSG SNe, SNe IIP and IIn. Our results suggest that LMC metal abundances play a role in the formation of some 1987A-like SNe. This would naturally fit in a single star scenario for the progenitors. However, the existence of two events at nearly solar metallicity suggests that also other channels, e.g. binarity, contribute to produce BSG SNe.Comment: 28 pages, 17 figures; accepted for publication (Astronomy and Astrophysics); abstract abridged for arXiv submissio

Whatever happened to the progenitors of supernovae 2008cn, 2009kr and 2009md?

We present new late-time, high-resolution observations of the sites of supernovae (SNe) 2008cn, 2009kr and 2009md, acquired with the Hubble Space Telescope. In all instances, significant flux from the SNe is still recovered at late times. We show that the previous identification of the progenitor of SN 2008cn was actually a blend of two sources, whose locations are resolved in these new observations. We suggest that the progenitor of SN 2008cn was actually a red supergiant with Minit < 16 M⊙. In the late-time observations of SN 2009kr, we find that the pre-explosion source (previously thought to be a yellow supergiant) is most probably a small compact cluster with mass ∼6000 M⊙. In late-time F814W observations of the site of SN 2009md, we find a single point source with identical brightness to the pre-explosion source, suggesting some caution in assuming that the pre-explosion source was the progenitor

Metallicity at the explosion sites of interacting transients

Context. Some circumstellar-interacting (CSI) supernovae (SNe) are produced by the explosions of massive stars that have lost mass shortly before the SN explosion. There is evidence that the precursors of some SNe IIn were luminous blue variable (LBV) stars. For a small number of CSI SNe, outbursts have been observed before the SN explosion. Eruptive events of massive stars are named as SN impostors (SN IMs) and whether they herald a forthcoming SN or not is still unclear. The large variety of observational properties of CSI SNe suggests the existence of other progenitors, such as red supergiant (RSG) stars with superwinds. Furthermore, the role of metallicity in the mass loss of CSI SN progenitors is still largely unexplored. Aims. Our goal is to gain insight on the nature of the progenitor stars of CSI SNe by studying their environments, in particular the metallicity at their locations. Methods. We obtain metallicity measurements at the location of 60 transients (including SNe IIn, SNe Ibn, and SN IMs), via emission-line diagnostic on optical spectra obtained at the Nordic Optical Telescope and through public archives. Metallicity values from the literature complement our sample. We compare the metallicity distributions among the different CSI SN subtypes and to those of other core-collapse SN types. We also search for possible correlations between metallicity and CSI SN observational properties. Results. We find that SN IMs tend to occur in environments with lower metallicity than those of SNe IIn. Among SNe IIn, SN IIn-L(1998S-like) SNe show higher metallicities, similar to those of SNe IIL/P, whereas long-lasting SNe IIn (1988Z-like) show lower metallicities, similar to those of SN IMs. The metallicity distribution of SNe IIn can be reproduced by combining the metallicity distributions of SN IMs (that may be produced by major outbursts of massive stars like LBVs) and SNe IIP (produced by RSGs). The same applies to the distributions of the Normalized Cumulative Rank (NCR) values, which quantifies the SN association to H II regions. For SNe IIn, we find larger mass-loss rates and higher CSM velocities at higher metallicities. The luminosity increment in the optical bands during SN IM outbursts tend to be larger at higher metallicity, whereas the SN IM quiescent optical luminosities tend to be lower. Conclusions. The difference in metallicity between SNe IIn and SN IMs suggests that LBVs are only one of the progenitor channels for SNe IIn, with 1988Z-like and 1998S-like SNe possibly arising from LBVs and RSGs, respectively. Finally, even though linedriven winds likely do not primarily drive the late mass-loss of CSI SN progenitors, metallicity has some impact on the observational properties of these transients. Key words. supernovae: general - stars: evolution - galaxies: abundancesComment: Submitted to Astronomy and Astrophysics on 28/02/2015; submitted to arXiv after the 1st referee repor

Three-component modelling of O-rich AGB star winds I. Effects of drift using forsterite

Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic model code T-800 that includes frequency-dependent radiative transfer, dust extinction based on Mie scattering, grain growth and ablation, gas-to-dust drift using one mean grain size, a piston that simulates stellar pulsations, and an accurate high spatial resolution numerical scheme. To enable this study, we calculated new gas opacities based on the ExoMol database, and we extended the model code to handle the formation of minerals that may form in M-type stars. We determine the effects of drift by comparing drift models to our new and extant non-drift models. Three out of four new drift models show high drift velocities, 87-310 km/s. Our new drift model mass-loss rates are 1.7-13 per cent of the corresponding values of our non-drift models, but compared to the results of two extant non-drift models that use the same stellar parameters, these same values are 0.33-1.5 per cent. Meanwhile, a comparison of other properties such as the expansion velocity and grain size show similar values. Our results, which are based on single-component forsterite particles, show that the inclusion of gas-to-drift is of fundamental importance in stellar wind models driven by such transparent grains. Assuming that the drift velocity is insignificant, properties such as the mass-loss rate may be off from more realistic values by a factor of 50 or more.Comment: 15 pages, 6 figures, accepted and in pres

The Unusually Luminous Extragalactic Nova SN 2010U

We present observations of the unusual optical transient SN 2010U, including spectra taken 1.03 days to 15.3 days after maximum light that identify it as a fast and luminous Fe II type nova. Our multi-band light curve traces the fast decline (t_2 = 3.5 days) from maximum light (M_V = -10.2 mag), placing SN 2010U in the top 0.5% of the most luminous novae ever observed. We find typical ejecta velocities of approximately 1100 km/s and that SN 2010U shares many spectral and photometric characteristics with two other fast and luminous Fe II type novae, including Nova LMC 1991 and M31N-2007-11d. For the extreme luminosity of this nova, the maximum magnitude vs. rate of decline relationship indicates a massive white dwarf progenitor with a low pre-outburst accretion rate. However, this prediction is in conflict with emerging theories of nova populations, which predict that luminous novae from massive white dwarfs should preferentially exhibit an alternate spectral type (He/N) near maximum light.Comment: 16 pages, 16 figures. Submitted to the Astrophysical Journa

The Type IIn Supernova SN 2010bt: The Explosion of a Star in Outburst

Indexación: Scopus.It is well known that massive stars (M > 8 M ) evolve up to the collapse of the stellar core, resulting in most cases in a supernova (SN) explosion. Their heterogeneity is related mainly to different configurations of the progenitor star at the moment of the explosion and to their immediate environments. We present photometry and spectroscopy of SN 2010bt, which was classified as a Type IIn SN from a spectrum obtained soon after discovery and was observed extensively for about 2 months. After the seasonal interruption owing to its proximity to the Sun, the SN was below the detection threshold, indicative of a rapid luminosity decline. We can identify the likely progenitor with a very luminous star (log L/L ≈ 7) through comparison of Hubble Space Telescope images of the host galaxy prior to explosion with those of the SN obtained after maximum light. Such a luminosity is not expected for a quiescent star, but rather for a massive star in an active phase. This progenitor candidate was later confirmed via images taken in 2015 (∼5 yr post-discovery), in which no bright point source was detected at the SN position. Given these results and the SN behavior, we conclude that SN 2010bt was likely a Type IIn SN and that its progenitor was a massive star that experienced an outburst shortly before the final explosion, leading to a dense H-rich circumstellar environment around the SN progenitor. © 2018. The American Astronomical Society. All rights reserved.https://iopscience.iop.org/article/10.3847/1538-4357/aac51

Late-time spectral line formation in Type IIb supernovae, with application to SN 1993J, SN 2008ax, and SN 2011dh

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modelling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date − SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modelling of the [O I] λλ6300, 6364 lines constrains the progenitors of these three SNe to the MZAMS = 12−16 M⊙ range (ejected oxygen masses 0.3−0.9 M⊙), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from MZAMS ≳ 17 M⊙ progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low-to-moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of 0.02−0.14 M⊙ is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO burning gives strong [N II] λλ6548, 6583 emission lines that dominate over Hα emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable Hα emission or absorption after ∼150 days, and nebular phase emission seen around 6550 Å is in many cases likely caused by [N II] λλ6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated. Significant line blocking in the metal core remains for several hundred days, which affects the emergent spectrum. These radiative transfer effects lead to early-time blueshifts of the emission line peaks, which gradually disappear as the optical depths decrease with time. The modelled evolution of this effect matches the observed evolution in SN 2011dh.Facultad de Ciencias Astronómicas y GeofísicasInstituto de Astrofísica de La Plat