Source author record

Luc Dessart

Luc Dessart appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2025arXiv

An optical to infrared study of type II SN2024ggi at nebular times

We present 0.3-21mic observations at ~275d and ~400d for Type II supernova (SN) 2024ggi, combining ground-based optical and near-infrared data from the Keck I/II telescopes and space-based infrared data from the James Webb Space Telescope. Although the optical regions dominate the observed flux, SN2024ggi is bright at infrared wavelengths (65%/35% falls each side of 1mic). SN2024ggi exhibits a plethora of emission lines from H, He, intermediate-mass elements (O, Na, Mg, S, Ar, Ca), and iron-group elements (IGEs; Fe, Co, and Ni) -- all lines have essentially the same width, suggesting efficient macroscopic chemical mixing of the inner ejecta at <~2000km/s and little mixing of 56Ni at larger velocities. Molecular emission in the infrared range is dominated by the CO fundamental, which radiates about 5% of the total SN luminosity. A molecule-free radiative-transfer model based on a standard red-supergiant star explosion (i.e., ~1e51erg, 0.06Msun of 56Ni from a 15.2Msun progenitor) yields a satisfactory match throughout the optical and infrared at both epochs. The SN2024ggi CO luminosity is comparable to the fractional decay-power absorbed in the model C/O-rich shell -- accounting for CO cooling would likely resolve the model overestimate of the [OI]0.632mic flux. The relative weakness of the molecular emission in SN2024ggi and the good overall match obtained with our molecule-free model suggests negligible microscopic mixing -- about 95% of the SN luminosity is radiated by atoms and ions. Lines from IGEs, which form from explosion ashes at such late times, are ideal diagnostics of the magnitude of 56Ni mixing in core-collapse SN ejecta. Stable Ni, clearly identified in SN2024ggi (e.g., [NiII]6.634mic), is probably a common product of massive-star explosions.

preprint2025arXiv

Infrared diagnostics of late-time core-collapse supernova spectra

We present nonlocal thermodynamic equilibrium radiative transfer calculations of red supergiant and He-star explosions, extending previous work to focus on the infrared emission from atoms and ions in the ejecta during the nebular-phase (i.e., ~200 to ~500d) -- molecules and dust are ignored. We cover non-rotating solar-metallicity progenitors spanning an initial mass between 10 and about 40Msun and exploding as Type II or Ibc supernovae (SNe). Both photometrically and spectroscopically, the SN II models evolve distinctly from the SN Ibc models primarily because of the greater ejecta kinetic-energy-to-mass ratio in the latter, which leads to a greater gamma-ray escape together with a lower density and a higher ionization in our H-deficient ejecta. Type II SN models remain optically luminous at all times, whereas Type Ibc models progressively brighten in the infrared (which holds 80% of their luminosity at 500d), causing strong infrared lines such as [NeII]12.810mic and [NiII]6.634mic to evolve essentially at constant luminosity. The strength of [NeII]12.810mic exhibits a complicated dependence with either He-core or preSN mass because of the additional impact of ejecta ionization -- this line radiates alone up to 20% of the SN luminosity after ~300d in our Type Ibc SN models. The numerous infrared Ni lines are found to be good tracers of the material that underwent explosive nucleosynthesis and can thus be used directly to constrain the level of 56Ni mixing in core-collapse SNe. The evolution of the integrated flux in infrared Fe and Co lines shows a great amount of diversity, which compromises their use as a diagnostic of the 56Ni-decay power source in our models. Future spectroscopic observations of core-collapse SNe by JWST will provide unprecedented information on the emission from atoms and ions in their ejecta, delivering critical constraints on massive star explosions.

preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program.

preprint2022arXiv

Modeling the signatures of interaction in Type II supernovae: UV emission, high-velocity features, broad-boxy profiles

Because mass loss is a fundamental phenomenon in massive stars, interaction with circumstellar material (CSM) should be universal in core-collapse supernovae (SNe). Leaving aside the extreme CSM density, extent, or mass typically encountered in Type IIn SNe, we investigate the diverse long-term radiative signatures of interaction between a Type II SN ejecta and CSM corresponding to mass loss rates up to 10$^{-3}$ $M_{\odot}$ yr$^{-1}$. Because these CSM are relatively tenuous and optically-thin to electron-scattering beyond a few stellar radii, radiation hydrodynamics is not essential and one may treat the interaction directly as an additional power source in the non-local thermodynamic equilibrium radiative transfer problem. The CSM accumulated since shock breakout forms a dense shell in the outer ejecta and leads to high-velocity absorption features in spectral lines, even for negligible shock power. Besides Balmer lines, such features may appear in NaID, HeI lines etc. A stronger interaction strengthens the continuum flux (preferentially in the UV), quenches the absorption of P-Cygni profiles, boosts the MgII $λλ$ $2795,2802$ doublet, and fosters the production of a broad boxy H$α$ emission component. The rise in ionization in the outer ejecta may quench some lines (e.g., the CaII near-infrared triplet). The interaction power emerges preferentially in the UV, in particular at later times, shifting the optical color to the blue, but increasing modestly the optical luminosity. Strong thermalization and clumping seem to be required to make an interaction superluminous in the optical. The UV range contains essential signatures that provide critical constraints to infer the mass loss history and inner workings of core-collapse SN progenitors at death.

preprint2022arXiv

Stable nickel production in Type Ia supernovae: A smoking gun for the progenitor mass?

At present, there are strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit (MCh ~ 1.4 Msun) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between MCh and sub-MCh events. In particular, it is thought that the higher-density ejecta of MCh WDs, which favours the synthesis of stable isotopes of nickel, results in prominent [Ni II] lines in late-time spectra. We study the explosive nucleosynthesis of stable nickel in SNe Ia resulting from MCh and sub-MCh progenitors. We explore the potential for lines of [Ni II] at 7378 Å and 1.94 microns in late-time spectra to serve as a diagnostic of the exploding WD mass, using nonlocal thermodynamic equilibrium radiative-transfer simulations with the CMFGEN code. We find that the radiative proton-capture reaction 57Co(p,gamma)58Ni is the dominant production mode for 58Ni in both MCh and sub-MCh models, while the alpha-capture reaction on 54Fe has a negligible impact on the final 58Ni yield. More importantly, we demonstrate that the lack of [Ni II] lines in late-time spectra of sub-MCh events is not always due to an under-abundance of stable Ni; rather, it results from the higher ionization of Ni in the inner ejecta. Conversely, the strong [Ni II] lines predicted in our 1D MCh models are completely suppressed when 56Ni is sufficiently mixed with the innermost layers, which are rich in stable iron-group elements. [Ni II] lines in late-time SN Ia spectra have a complex dependency on the abundance of stable Ni, which limits their use in distinguishing among MCh and sub-MCh progenitors. However, we argue that a low-luminosity SN Ia displaying strong [Ni II] lines would most likely result from a Chandrasekhar-mass progenitor. [Abridged]

preprint2022arXiv

Stripped-Envelope Stars in Different Metallicity Environments I. Evolutionary Phases, Classification and Populations

Massive stars that become stripped of their hydrogen envelope through binary interaction or winds can be observed either as Wolf-Rayet stars, if they have optically thick winds, or as transparent-wind stripped-envelope stars. We approximate their evolution through evolutionary models of single helium stars, and compute detailed model grids in the initial mass range 1.5 to 70 M$_{\odot}$ for metallicities between 0.01 and 0.04, from core helium ignition until core collapse. Throughout their lifetime, some stellar models expose the ashes of helium burning. We propose that models that have nitrogen-rich envelopes are candidate WN stars, while models with a carbon-rich surface are candidate WC stars during core helium burning, and WO stars afterwards. We measure metallicity dependance of the total lifetime of our models and the duration of their evolutionary phases. We propose an analytic estimate of the wind optical depth to distinguish models of Wolf-Rayet stars from transparent-wind stripped-envelope stars, and find that the luminosity ranges at which WN, WC and WO type stars can exist is a strong function of metallicity. We find that all carbon-rich models produced in our grids have optically thick winds and match the luminosity distribution of observed populations. We construct population models and predict the numbers of transparent-wind stripped-envelope stars and Wolf-Rayet stars, and derive their number ratios at different metallicities. We find that as metallicity increases, the number of transparent-wind stripped-envelope stars decreases and the number of Wolf-Rayet stars increases. At high metallicities WC and WO type stars become more common. We apply our population models to nearby galaxies, and find that populations are more sensitive to the transition luminosity between Wolf-Rayet stars and transparent-wind helium stars than to the metallicity dependent mass loss rates.

preprint2021arXiv

Final Moments I: Precursor Emission, Envelope Inflation, and Enhanced Mass loss Preceding the Luminous Type II Supernova 2020tlf

We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discovery, "flash" spectroscopy of SN 2020tlf revealed prominent narrow symmetric emission lines ($v_w < 300$ km s$^{-1}$) that resulted from the photo-ionization of unshocked circumstellar material (CSM) shedded in progenitor mass loss episodes in the final weeks to months before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a RSG progenitor with ZAMS mass of only 10-12 M$_{\odot}$, as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative transfer code CMFGEN and radiation-hydrodynamical (RHD) code HERACLES suggests a dense CSM limited to $r \approx 10^{15}$ cm, and mass loss rate of $10^{-2}$ M$_{\odot}$ yr$^{-1}$. The subsequent luminous light-curve plateau and persistent blue excess indicates an extended progenitor, compatible with a RSG model with $R_{\star} = 1100$ R$_{\odot}$. Inferences from the limits on the shock-powered X-ray and radio luminosity are consistent with these conclusions and suggest a CSM density of $ρ< 2 \times 10^{-16}$ g cm$^{-3}$ for distances of $r \approx 5 \times 10^{15}$ cm, as well as a mass loss rate of $\dot M<1.3 \times 10^{-5}\,\rm{M_{\odot}\,yr^{-1}}$ at larger distances. A promising power source for the observed precursor emission is the ejection of stellar material following energy disposition into the stellar envelope as a result of gravity waves emitted during either neon/oxygen burning or a nuclear flash from silicon combustion.

preprint2021arXiv

Helium stars exploding in circumstellar material and the origin of Type Ibn supernovae

Type Ibn supernovae (SNe) are a mysterious class of transients whose spectra exhibit persistently narrow HeI lines, and whose bolometric light curves are typically fast evolving and overluminous at peak relative to standard Type Ibc SNe. We explore the interaction scenario of such Type Ibn SNe by performing radiation-hydrodynamics and radiative-transfer calculations. We find that standard-energy helium-star explosions within dense wind-like circumstellar material (CSM) can reach on day timescales a peak luminosity of a few 10^44 erg/s, reminiscent of exceptional events like AT2018cow. Similar interactions but with weaker winds can lead to Type Ibc SNe with double-peak light curves and peak luminosities in the range ~10^42.2 to ~10^43 erg/s. In contrast, the narrow spectral lines and modest peak luminosities of most Type Ibn SNe are suggestive of a low-energy explosion in an initially <~5 Msun helium star, most likely arising from interacting binaries, and colliding with a massive helium-rich, probably ejecta-like, CSM at ~10^15 cm. Nonlocal thermodynamic equilibrium radiative-transfer simulations of a slow-moving dense shell born out and powered by the interaction compare favorably to Type Ibn SNe like 2006jc, 2011hw, or 2018bcc at late times and suggest a composition made of about 50% helium, a solar metallicity, and a total ejecta/CSM mass of 1-2Msun. A lower fractional helium abundance leads to weak or absent HeI lines and thus excludes more massive configurations for observed Type Ibn SNe. Further, the dominance of FeII emission below 5500A seen in Type Ibn SNe at late times is not predicted at low metallicity. Hence, despite their promising properties, Type Ibn SNe from pulsational-pair instability in very massive stars, which requires low metallicity, have probably not yet been observed.

preprint2020arXiv

Spectral signatures of H-rich material stripped from a non-degenerate companion by a Type Ia supernova

The single-degenerate scenario for Type Ia supernovae (SNe Ia) should yield metal-rich ejecta that enclose some stripped material from the non-degenerate H-rich companion star. We present a large grid of non-local thermodynamic equilibrium steady-state radiative transfer calculations for such hybrid ejecta and provide analytical fits for the Halpha luminosity and equivalent width. Our set of models covers a range of masses for 56Ni and the ejecta, for the stripped material (Mst), and post-explosion epochs from 100 to 300d. The brightness contrast between stripped material and metal-rich ejecta challenges the detection of HI and HeI lines prior to ~100d. Intrinsic and extrinsic optical depth effects also influence the radiation emanating from the stripped material. This inner denser region is marginally thick in the continuum and optically thick in all Balmer lines. The overlying metal-rich ejecta blanket the inner regions, completely below about 5000A, and more sparsely at longer wavelengths. As a consequence, Hbeta should not be observed for all values of Mst through at least 300 days, while Halpha should be observed after ~100d for all Mst >= 0.01Msun. This contrasts with the case of circumstellar (CSM) interaction, not subject to external blanketing, which should produce Halpha and Hbeta lines with a strength dependent primarily on CSM density. We confirm previous analyses that suggest low values of order 0.001Msun for Mst to explain the observations of the two SNe Ia with nebular-phase Halpha detection, in conflict with the much greater stripped mass predicted by hydrodynamical simulations for the single-degenerate scenario. A more likely solution is the double-degenerate scenario, together with CSM interaction, or enclosed material from a tertiary star in a triple system or from a giant planet. [Abridged]

preprint2019arXiv

Understanding nebular spectra of Type Ia supernovae

In this study, we present one-dimensional, non-local-thermodynamic-equilibrium, radiative transfer simulations (using CMFGEN) in which we introduce micro-clumping at nebular times into two Type Ia supernova ejecta models. We use one sub-Chandrasekhar (sub-M$_{\rm Ch}$) ejecta with 1.02 M$_\odot$ and one M$_{\rm Ch}$ ejecta model with 1.40 M$_\odot$. We introduce clumping factors $f=$0.33,0.25, and 0.10 which are constant throughout the ejecta and compared to the unclumped $f=1.0$ case. We find that clumping is a natural mechanism to reduce the ionization of the ejecta, reducing emission from [Fe III], [Ar III], and [S III] by a factor of a few. For decreasing values of the clumping factor $f$, the [Ca II] $λλ$7291,7324 doublet became a dominant cooling line for our M$_{\rm Ch}$ model but still weak in our sub-M$_{\rm Ch}$ model. Strong [Ca II] $λλ$7291,7324 indicates non-thermal heating in that region and may constrain explosion modelling. Due to the low abundance of stable nickel, our sub-M$_{\rm Ch}$ model never showed the [Ni II] 1.939 micron diagnostic feature for all clumping values.

preprint2016arXiv

Inferring supernova IIb/Ib/Ic ejecta properties from light curves and spectra: Correlations from radiative-transfer models

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875Å at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772Å line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$μ$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

preprint2016arXiv

Models of interacting supernovae and their spectral diversity

Using radiation-hydrodynamics and radiative-transfer simulations, we explore the origin of the spectral diversity of interacting supernovae (SNe) of type IIn. We revisit SN1994W and investigate the dynamical configurations that can give rise to spectra with narrow lines at all times. We find that a standard ~10Msun 10^51erg SN ejecta ramming into a 0.4Msun dense CSM is inadequate for SN1994W, as it leads to the appearance of broad lines at late times. This structure, however, generates spectra that exhibit the key morphological changes seen in SN1998S. For SN1994W, we consider a completely different configuration, which involves the interaction at a large radius of a low mass inner shell with a high mass outer shell. Such a structure may arise in an 8-12Msun star from a nuclear flash (e.g., of Ne) followed within a few years by core collapse. Our simulations show that the large mass of the outer shell leads to the complete braking of the inner shell material, the formation of a slow dense shell, and the powering of a luminous SN IIn, even for a low inner shell energy. Early on, our model line profiles are typical of SNe IIn, exhibiting narrow cores and broad electron-scattering wings. As observed in SN1994W, they also remain narrow at late times. Our SN1994W model invokes two low energy ejections, both atypical of observed massive stars, and illustrates the diversity of configurations leading to SNe IIn. These results also highlight the importance of spectra to constrain the dynamical properties and understand the origin of SNe IIn.

preprint2016arXiv

On Variations Of Pre-Supernova Model Properties

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

preprint2016arXiv

Two-dimensional radiation-hydrodynamics simulations of super-luminous interacting supernovae of type IIn

Some interacting supernovae (SNe) of type IIn show a sizeable continuum polarisation suggestive of a large scale asymmetry in the circumstellar medium (CSM) and/or the SN ejecta. Here, we extend the recent work of Dessart et al. on super-luminous SNe IIn and perform axially-symmetric (i.e., 2D) multi-group radiation hydrodynamics simulations to explore the impact of an imposed large scale density asymmetry. When the CSM is asymmetric, the latitudinal variation of the radial optical depth $τ$ introduces a strong flux redistribution from the higher-density CSM regions, where the shock luminosity is larger, towards the lower-density CSM regions where photons escape more freely --- this redistribution ceases when $τ<$1. Along directions where the CSM density is larger, the shock deceleration is stronger and its progression slower, producing a non-spherical cold-dense shell (CDS). For an oblate CSM density distribution, the photosphere (CDS) has an oblate (prolate) morphology when $τ>$1. When the CSM is symmetric and the ejecta asymmetric, the flux redistribution within the CSM now tends to damp the latitudinal variation of the luminosity at the shock. It then requires a larger ejecta asymmetry to produce a sizeable latitudinal variation in the emergent flux. When the interaction is between a SN ejecta and a relic disk, the luminosity boost at early times scales with the disk opening angle -- forming a super-luminous SN IIn this way requires an unrealistically thick disk. In contrast, interaction with a disk of modest thickness/mass can yield a power that rivals radioactive decay of a standard SN II at nebular times.

preprint2015arXiv

A correlation between SN II metal line strength and host HII region oxygen abundance

Dessart et al., demonstrated that type II supernova (SN II) model spectra present increasing metal line strength with increasing progenitor metallicity. To confront these models with observations, we obtained a large sample of SN II host HII region emission line spectroscopy. We show that inferred SN II host HII region metallicities have a statistically significant correlation with the strength of SN II metal lines, specifically FeII 5018A.

preprint2015arXiv

A one-dimensional Chandrasekhar-mass delayed-detonation model for the broad-lined Type Ia supernova 2002bo

We present 1D non-local thermodynamic equilibrium (non-LTE) time-dependent radiative-transfer simulations of a Chandrasekhar-mass delayed-detonation model which synthesizes 0.51 Msun of 56Ni, and confront our results to the Type Ia supernova (SN Ia) 2002bo over the first 100 days of its evolution. Assuming only homologous expansion, this same model reproduces the bolometric and multi-band light curves, the secondary near-infrared (NIR) maxima, and the optical and NIR spectra. The chemical stratification of our model qualitatively agrees with previous inferences by Stehle et al., but reveals significant quantitative differences for both iron-group and intermediate-mass elements. We show that +/-0.1 Msun (i.e., +/-20 per cent) variations in 56Ni mass have a modest impact on the bolometric and colour evolution of our model. One notable exception is the U-band, where a larger abundance of iron-group elements results in less opaque ejecta through ionization effects, our model with more 56Ni displaying a higher near-UV flux level. In the NIR range, such variations in 56Ni mass affect the timing of the secondary maxima but not their magnitude, in agreement with observational results. Moreover, the variation in the I, J, and K_s magnitudes is less than 0.1 mag within ~10 days from bolometric maximum, confirming the potential of NIR photometry of SNe Ia for cosmology. Overall, the delayed-detonation mechanism in single Chandrasekhar-mass white dwarf progenitors seems well suited for SN 2002bo and similar SNe Ia displaying a broad Si II 6355 A line. Whatever multidimensional processes are at play during the explosion leading to these events, they must conspire to produce an ejecta comparable to our spherically-symmetric model.

preprint2015arXiv

Numerical simulations of super-luminous supernovae of type IIn

We present numerical simulations that include 1-D Eulerian multi-group radiation-hydrodynamics, 1-D non-LTE radiative transfer, and 2-D polarised radiative transfer for super-luminous interacting supernovae (SNe). Our reference model is a ~10Msun inner shell with 10^51erg ramming into a ~3Msun cold outer shell (the circumstellar-medium, or CSM) that extends from 10^15cm to 2x10^16cm and moves at 100km/s. We discuss the light curve evolution, which cannot be captured adequately with a grey approach. In these interactions, the shock-crossing time through the optically-thick CSM is much longer than the photon diffusion time. Radiation is thus continuously leaking from the shock through the CSM, in disagreement with the shell-shocked model that is often invoked. Our spectra redden with time, with a peak distribution in the near-UV during the first month gradually shifting to the optical range over the following year. Initially Balmer lines exhibit a narrow line core and the broad line wings that are characteristic of electron scattering in the SNe IIn atmospheres (CSM). At later times they also exhibit a broad blue shifted component which arises from the cold dense shell. Our model results are broadly consistent with the bolometric light curve and spectral evolution observed for SN2010jl. Invoking a prolate pole-to-equator density ratio in the CSM, we can also reproduce the ~2% continuum polarisation, and line depolarisation, observed in SN2010jl. By varying the inner shell kinetic energy and the mass and extent of the outer shell, a large range of peak luminosities and durations, broadly compatible with super-luminous SNe IIn like 2010jl or 2006gy, can be produced.

preprint2015arXiv

Radiative-transfer models for supernovae IIb/Ib/Ic from binary-star progenitors

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

preprint2014arXiv

[CoIII] versus NaID in type Ia supernova spectra

The high metal content and fast expansion of supernova (SN) Ia ejecta lead to considerable line overlap in their optical spectra. Uncertainties in composition and ionization further complicate the process of line identification. In this paper, we focus on the 5900A emission feature seen in SN Ia spectra after bolometric maximum, a line which in the last two decades has been associated with [CoIII]5888A or NaID. Using non-LTE time-dependent radiative-transfer calculations based on Chandrasekhar-mass delayed-detonation models, we find that NaID line emission is extremely weak at all post-maximum epochs. Instead, we predict the presence of [CoIII]5888A after maximum in all our SN Ia models, which cover a range from 0.12 to 0.87Msun of 56Ni. We also find that the [CoIII]5888A forbidden line is present within days of bolometric maximum, and strengthens steadily for weeks thereafter. Both predictions are confirmed by observations. Rather than trivial taxonomy, these findings confirm that it is necessary to include forbidden-line transitions in radiative-transfer simulations of SNe Ia, both to obtain the correct ejecta cooling rate and to match observed optical spectra.

preprint2014arXiv

Constraints on the explosion mechanism and progenitors of type Ia supernovae

Observations of SN 2011fe at early times reveal an evolution analogous to a fireball model of constant color. In contrast, our unmixed delayed detonations of Chandrasekhar-mass white dwarfs (DDC series) exhibit a faster brightening concomitant with a shift in color to the blue. In this paper, we study the origin of these discrepancies. We find that strong chemical mixing largely resolves the photometric mismatch at early times, but it leads to an enhanced line broadening that contrasts, for example, with the markedly narrow SiII6355A line of SN 2011fe. We also explore an alternative configuration with pulsational-delayed detonations (PDDEL model series). Because of the pulsation, PDDEL models retain more unburnt carbon, have little mass at high velocity, and have a much hotter outer ejecta after the explosion. The pulsation does not influence the inner ejecta, so PDDEL and DDC models exhibit similar radiative properties beyond maximum. However, at early times, PDDEL models show bluer optical colors and a higher luminosity, even for weak mixing. Their early-time radiation is derived primarily from the initial shock-deposited energy in the outer ejecta rather than radioactive decay heating. Furthermore, PDDEL models show short-lived CII lines, reminiscent of SN 2013dy. They typically exhibit lines that are weaker, narrower, and of near-constant width, reminiscent of SN 2011fe. In addition to multi-dimensional effects, varying configurations for such ``pulsations" offer a source of spectral diversity amongst SNe Ia. PDDEL and DDC models also provide one explanation for low- and high-velocity gradient SNe Ia.

preprint2014arXiv

Critical ingredients of supernova Ia radiative-transfer modeling

We explore the physics of SN Ia light curves and spectra using the 1-D non-LTE time-dependent radiative-transfer code CMFGEN. Rather than adjusting ejecta properties to match observations, we select as input one "standard" 1-D Chandrasekhar-mass delayed-detonation hydrodynamical model, and then explore the sensitivity of radiation and gas properties on radiative-transfer modeling assumptions. The correct computation of SN Ia radiation is not exclusively a solution to an "opacity problem", characterized by the treatment of a large number of lines. It is also key to treat important atomic processes consistently. Besides handling line blanketing in non-LTE, we show that including forbidden line transitions of metals is increasingly important for the temperature and ionization of the gas beyond maximum light. Non-thermal ionization and excitation are also critical since they affect the color evolution and the Delta-M15 of our model. While impacting little the bolometric luminosity, a more complete treatment of decay routes leads to enhanced line blanketing, e.g., associated with 48Ti in the U and B bands. Overall, we find that SN Ia radiation properties are influenced in a complicated way by the atomic data we employ, so that obtaining converged results is a challenge. We nonetheless obtain a good match to the golden standard type Ia SN 2005cf in the optical and near-IR, from 5 to 60d after explosion, suggesting that assuming spherical symmetry is not detrimental to SN Ia radiative-transfer modeling at these times. Multi-D effects no doubt matter, but they are perhaps less important than accurately treating non-LTE processes [abridged].

preprint2014arXiv

One-dimensional non-LTE time-dependent radiative transfer of an He-detonation model and the connection to faint and fast-decaying supernovae

We present non-LTE time-dependent radiative transfer simulations for ejecta produced by the detonation of an helium shell at the surface of a low-mass carbon/oxygen white dwarf (WD). This mechanism is one possible origin for supernovae (SNe) with faint and fast-decaying light curves, such as .Ia SNe and Ca-rich transients. Our initial ejecta conditions at 1d are given by the 0.18B explosion model COp45HEp2 of Waldman et al.. The 0.2Msun ejecta initially contains 0.11Msun of He, 0.03Msun of Ca, and 0.03Msun of Ti. We obtain a ~5d rise to a bolometric maximum of 3.59x10^41 erg/s, primarily powered by 48V decay. Multi-band light curves show distinct morphologies, with a rise to maximum magnitude (-14.3 to -16.7mag) that varies between 3 to 9d from the U to the K bands. Near-IR light curves show no secondary maximum. Because of the presence of both HeI and SiII lines at early times we obtain a hybrid Type Ia/Ib classification. During the photospheric phase line blanketing is caused primarily by TiII. At nebular times, the spectra show strong CaII lines in the optical (but no [OI] 6300--6364A emission), and TiII in the near-IR. Overall, these results match qualitatively the very disparate properties of .Ia SNe and Ca-rich transients. Although the strong TiII blanketing and red colors that we predict are rarely observed, they are seen, for example, in OGLE-2013- SN-079. Furthermore, we obtain a faster light-curve evolution than, for example, PTF10iuv, indicating an ejecta mass >0.2Msun. An alternate scenario may be the merger of two WDs, one or both composed of He.

preprint2013arXiv

The Arduous Journey to Black-Hole Formation in Potential Gamma-Ray Burst Progenitors

We present a quantitative study on the properties at death of fast-rotating massive stars evolved at low-metallicity, objects that are proposed as likely progenitors of long-duration gamma-ray bursts (LGRBs). We perform 1D+rotation stellar-collapse simulations on the progenitor models of Woosley & Heger (2006) and critically assess their potential for the formation of a black hole and a Keplerian disk (namely a collapsar) or a proto-magnetar. We note that theoretical uncertainties in the treatment of magnetic fields and the approximate handling of rotation compromises the accuracy of stellar-evolution models. We find that only the fastest rotating progenitors achieve sufficient compactness for black-hole formation while the bulk of models possess a core density structure typical of garden-variety core-collapse supernova (SN) progenitors evolved without rotation and at solar metallicity. Of the models that do have sufficient compactness for black-hole formation, most of them also retain a large amount of angular momentum in the core, making them prone to a magneto-rotational explosion, therefore preferentially leaving behind a proto-magnetar. A large progenitor angular-momentum budget is often the sole criterion invoked in the community today to assess the suitability for producing a collapsar. This simplification ignores equally important considerations such as the core compactness, which conditions black-hole formation, the core angular momentum, which may foster a magneto-rotational explosion preventing black-hole formation, or the metallicity and the residual envelope mass which must be compatible with inferences from observed LGRB/SNe. Our study suggests that black-hole formation is non trivial, that there is room for accommodating both collapsars and proto-magnetars as LGRB progenitors, although proto-magnetars seem much more easily produced by current stellar-evolutionary models.

preprint2013arXiv

Type II-Plateau supernova radiation: dependencies on progenitor and explosion properties

We explore the properties of Type II-Plateau (II-P) supernovae (SNe) together with their red-supergiant (RSG) star progenitors. Using MESA STAR, we modulate the parameters (e.g., mixing length, overshoot, rotation, metallicity) that control the evolution of a 15Msun main-sequence star to produce a variety of physical pre-SN models and SN II-P ejecta. We extend previous modeling of SN II-P radiation to include photospheric and nebular phases, as well as multi-band light curves and spectra. Our treatment does not assume local thermodynamic equilibrium, is time dependent, treats explicitly the effects of line blanketing, and incorporates non-thermal processes. We find that the color properties of SNe II-P require large model atoms for FeI and FeII, much larger than adopted in Dessart & Hillier (2011). The color properties also imply RSG progenitors of limited extent (~500Rsun) --- larger progenitor stars produce a SN II-P radiation that remains too blue for too long. This finding calls for a reduction of RSG radii, perhaps through a strengthening of convective energy transport in RSG envelopes. Increased overshoot and rotation reduce the ratio of ejecta to helium-core mass, similarly to an increase in main-sequence mass, and thus complicate the inference ofprogenitor masses. In contrast to the great sensitivity on progenitor radius, SN II-P color evolution appears insensitive to variations in explosion energy. Finally, we document the numerous SN II-P signatures that vary with progenitor metallicity, revealing their potential for metallicity determinations in the nearby and distant Universe.

preprint2012arXiv

Non-thermal excitation and ionization in supernovae

We incorporate non-thermal excitation and ionization processes arising from non-thermal electrons that result from γ-ray energy deposition, into our radiative transfer code CMFGEN. The non-thermal electron distribution is obtained by solving the Spencer-Fano equation using the procedure of Kozma & Fransson (1992). We applied the non-thermal calculations to the blue supergiant explosion model whose early evolution was studied in Dessart & Hillier (2010). Non-thermal processes generally increase excitation and ionization and decrease the temperature of the ejecta. We confirm that non-thermal processes are crucial for modeling the nebular spectra. Both optical HI and HeI lines are significantly strengthened. While optical HeI lines are not easily discerned in observational spectra due to severe blending with other lines, HeI 2.058 μm provides an excellent opportunity to infer the influence of non-thermal processes. We also discuss the processes controlling the formation of the HeI lines during the nebular epoch. Most lines of other species are only slightly affected. We also show that the inclusion of FeI has substantial line-blanketing effects on the optical spectra. Our model spectra and synthetic light curves are compared to the observations of SN 1987A. The spectral evolution shows broad agreement with the observations, especially Hα. The uncertainties of the non-thermal solver are studied, and are expected to be small. With this new addition of non-thermal effects in CMFGEN, we now treat all known important processes controlling the radiative transfer of a supernova ejecta, whatever the type and the epoch.

preprint2012arXiv

On the nature of supernovae Ib and Ic

Using non-LTE time-dependent radiative-transfer calculations, we study the impact of mixing and non-thermal processes associated with radioactive decay on SN IIb/Ib/Ic light curves (LCs) and spectra. Starting with short-period binary models of \leq5Msun He-rich stars (18-25Msun on the main-sequence), we produce 1.2B ejecta which we artificially mix to alter the chemical stratification. While the total 56Ni mass influences the LC peak, the spatial distribution of 56Ni, controlled by mixing processes, impacts both the multi-band LCs and spectra. With enhanced mixing, our synthetic LCs start their post-breakout re-brightening phase earlier, follow a more gradual rise to peak, appear redder, and fade faster after peak due to enhanced gamma-ray escape. Non-thermal electrons, crucial for the production of HeI lines, deposit a dominant fraction of their energy as heat. Because energy deposition is generally local well after the LC peak, the broad HeI lines characteristic of maximum-light SN IIb/Ib spectra require mixing that places 56Ni and helium nuclei to within a gamma-ray mean-free-path. This requirement indicates that SNe IIb and Ib most likely arise from the explosion of stripped-envelope massive stars (main-sequence masses \leq25Msun) that have evolved through mass-transfer in a binary system, rather than from more massive single WR stars. In contrast, the lack of HeI lines in SNe Ic may result from a variety of causes: A genuine helium deficiency; strongly-asymmetric mixing; weak mixing; or a more massive, perhaps single, progenitor characterized by a larger oxygen-rich core. Our models, subject to different mixing magnitudes, can produce a variety of SN types, including IIb, IIc, Ib, and Ic. As it is poorly constrained by explosion models, mixing challenges our ability to infer the progenitor and explosion properties of SNe IIb/Ib/Ic.

preprint2012arXiv

One-dimensional delayed-detonation models of Type Ia supernovae: Confrontation to observations at bolometric maximum

The delayed-detonation explosion mechanism applied to a Chandrasekhar-mass white dwarf offers a very attractive model to explain the inferred characteristics of Type Ia supernovae (SNe Ia). The resulting ejecta are chemically stratified, have the same mass and roughly the same asymptotic kinetic energy, but exhibit a range in 56Ni mass. We investigate the contemporaneous photometric and spectroscopic properties of a sequence of delayed-detonation models, characterized by 56Ni masses between 0.18 and 0.81 Msun. Starting at 1d after explosion, we perform the full non-LTE, time-dependent radiative transfer with the code CMFGEN, with an accurate treatment of line blanketing, and compare our results to SNe Ia at bolometric maximum. Despite the 1D treatment, our approach delivers an excellent agreement to observations. We recover the range of SN Ia luminosities, colours, and spectral characteristics from the near-UV to 1 micron, for standard as well as low-luminosity 91bg-like SNe Ia. Our models predict an increase in rise time to peak with increasing 56Ni mass, from ~15 to ~21d, yield peak bolometric luminosities that match Arnett's rule to within 10%, and reproduce the much smaller scatter in near-IR magnitudes compared to the optical. We reproduce the morphology of individual spectral features, the stiff dependence of the R(Si) spectroscopic ratio on 56Ni mass, and the onset of blanketing from TiII/ScII in low-luminosity SNe Ia with a 56Ni mass <0.3 Msun. We find that ionization effects, which often dominate over abundance variations, can produce high-velocity features in CaII lines, even in 1D. Distinguishing between different SN Ia explosion mechanisms is a considerable challenge but the results presented here provide additional support to the viability of the delayed-detonation model.

preprint2012arXiv

Radiative Properties of Pair-instability Supernova Explosions

We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.

preprint2012arXiv

Super-luminous supernovae: 56Ni power versus magnetar radiation

Much uncertainty surrounds the origin of super-luminous supernovae (SNe). Motivated by the discovery of the Type Ic SN2007bi, we study its proposed association with a pair-instability SN (PISN). We compute stellar-evolution models for primordial ~200Msun stars, simulating the implosion/explosion due to the pair-production instability, and use them as inputs for detailed non-LTE time-dependent radiative-transfer simulations that include non-local energy deposition and non-thermal processes. We retrieve the basic morphology of PISN light curves from red-supergiant, blue-supergiant, and Wolf-Rayet (WR) star progenitors. Although we confirm that a progenitor 100Msun helium core (PISN model He100) fits well the SN2007bi light curve, the low ratios of its kinetic energy and 56Ni mass to the ejecta mass, similar to standard core-collapse SNe, conspire to produce cool photospheres, red spectra subject to strong line blanketing, and narrow line profiles, all conflicting with SN2007bi observations. He-core models of increasing 56Ni-to-ejecta mass ratio have bluer spectra, but still too red to match SN2007bi, even for model He125 -- the effect of 56Ni heating is offset by the associated increase in blanketing. In contrast, the delayed injection of energy by a magnetar represents a more attractive alternative to reproduce the blue, weakly-blanketed, and broad-lined spectra of super-luminous SNe. The extra heat source is free of blanketing and is not explicitly tied to the ejecta. Experimenting with a ~9Msun WR-star progenitor, initially exploded to yield a ~1.6B SN Ib/c ejecta but later influenced by tunable magnetar-like radiation, we produce a diversity of blue spectral morphologies reminiscent of SN2007bi, the peculiar Type Ib SN2005bf, and super-luminous SN2005ap-like events.

preprint2012arXiv

Time Dependent Radiative Transfer Calculations for Supernovae

In previous papers we discussed results from fully time-dependent radiative transfer models for core-collapse supernova (SN) ejecta, including the Type II-peculiar SN 1987A, the more "generic" SN II-Plateau, and more recently Type IIb/Ib/Ic SNe. Here we describe the modifications to our radiative modeling code, CMFGEN, which allowed those studies to be undertaken. The changes allow for time-dependent radiative transfer of SN ejecta in homologous expansion. In the modeling we treat the entire SN ejecta, from the innermost layer that does not fall back on the compact remnant out to the progenitor surface layers. From our non-LTE time-dependent line-blanketed synthetic spectra, we compute the bolometric and multi-band light curves: light curves and spectra are thus calculated simultaneously using the same physical processes and numerics. These upgrades, in conjunction with our previous modifications which allow the solution of the time dependent rate equations, will improve the modeling of SN spectra and light curves, and hence facilitate new insights into SN ejecta properties, the SN progenitors and the explosion mechanism(s). CMFGEN can now be applied to the modeling of all SN types

preprint2011arXiv

A Spectropolarimetric Comparison of the Type II-Plateau Supernovae SN 2008bk and SN 2004dj

The Type II-Plateau supernova (SN II-P) SN 2004dj was the first SN II-P for which spectropolarimetry data were obtained with fine temporal sampling before, during, and after the fall off of the photometric plateau -- the point that marks the transition from the photospheric to the nebular phase in SNe II-P. Unpolarized during the plateau, SN 2004dj showed a dramatic spike in polarization during the descent off of the plateau, and then exhibited a smooth polarization decline over the next two hundred days. This behavior was interpreted by Leonard et al. (2006) as evidence for a strongly non-spherical explosion mechanism that had imprinted asphericity only in the innermost ejecta. In this brief report, we compare nine similarly well-sampled epochs of spectropolarimetry of the Type II-P SN 2008bk to those of SN 2004dj. In contrast to SN 2004dj, SN 2008bk became polarized well before the end of the plateau and also retained a nearly constant level of polarization through the early nebular phase. Curiously, although the onset and persistence of polarization differ between the two objects, the detailed spectropolarimetric characteristics at the epochs of recorded maximum polarization for the two objects are extremely similar, feature by feature. We briefly interpret the data in light of non-Local-Thermodynamic Equilibrium, time-dependent radiative-transfer simulations specifically crafted for SN II-P ejecta.

preprint2011arXiv

Core-collapse explosions of Wolf-Rayet stars and the connection to type IIb/Ib/Ic supernovae

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

preprint2011arXiv

Synthetic line and continuum linear-polarisation signatures of axisymmetric type II supernova ejecta

We present synthetic single-line and continuum linear-polarisation signatures due to electron scattering in axially-symmetric Type II supernovae (SNe) which we calculate using a Monte Carlo and a long-characteristic radiative-transfer code. Aspherical ejecta are produced by prescribing a latitudinal scaling or stretching of SN ejecta inputs obtained from 1-D non-LTE time-dependent calculations. We study polarisation signatures as a function of inclination, shape factor, wavelength, line identity, post-explosion time. At early times, cancellation and optical-depth effects make the polarisation intrinsically low, causing complicated sign reversals with inclination or continuum wavelength, and across line profiles. While the line polarisation is positive (negative) for an oblate (prolate) morphology at the peak and in the red wing, the continuum polarisation may be of any sign. These complex polarisation variations are produced not just by the asymmetric distribution of scatterers but also of the flux. Our early-time signatures are in contradiction with predictions for a centrally illuminated aspherical nebula, although this becomes a better approximation at nebular times. For a fixed asymmetry, our synthetic continuum polarisation is generally low, may evolve non-monotonically during the plateau phase, but it systematically rises as the ejecta become optically thin. Changes in polarization over time do not necessarily imply a change in the asymmetry of the ejecta. The SN structure (e.g., density/ionization) critically influences the level of polarisation. Importantly, a low polarisation (<0.5%) at early times does not necessarily imply a low degree of asymmetry as usually assumed. Asphericity influences line-profile morphology and the luminosity, which may compromise the accuracy of SN characteristics inferred from these.

preprint2010arXiv

Determining the main-sequence mass of Type II supernova progenitors

We present radiation-hydrodynamics simulations of core-collapse supernova (SN) explosions, artificially generated by driving a piston at the base of the envelope of a rotating or non-rotating red-supergiant progenitor star. We search for trends in ejecta kinematics in the resulting Type II-Plateau (II-P) SN, exploring dependencies with explosion energy and pre-SN stellar-evolution model. We recover the trivial result that larger explosion energies yield larger ejecta velocities in a given progenitor. However, we emphasise that for a given explosion energy, the increasing helium-core mass with main-sequence mass of such Type II-P SN progenitors leads to ejection of core-embedded oxygen-rich material at larger velocities. We find that the photospheric velocity at 15d after shock breakout is a good and simple indicator of the explosion energy in our selected set of pre-SN models. This measurement, combined with the width of the nebular-phase OI6303-6363A line, can be used to place an upper-limit on the progenitor main-sequence mass. Using the results from our simulations, we find that the current, but remarkably scant, late-time spectra of Type II-P SNe support progenitor main-sequence masses inferior to ~20Msun and thus, corroborate the inferences based on the direct, but difficult, progenitor identification in pre-explosion images. The narrow width of OI6303-6363A in Type II-P SNe with nebular spectra does not support high-mass progenitors in the range 25-30Msun. Combined with quantitative spectroscopic modelling, such diagnostics offer a means to constrain the main-sequence mass of the progenitor, the mass fraction of the core ejected, and thus, the mass of the compact remnant formed.

preprint2010arXiv

High spatial resolution monitoring of the activity of BA supergiant winds

There are currently two optical interferometry recombiners that can provide spectral resolutions better than 10000, AMBER/VLTI operating in the H-K bands, and VEGA/CHARA, recently commissioned, operating in the visible. These instruments are well suited to study the wind activity of the brightest AB supergiants in our vicinity, in lines such as H$α$ or BrGamma. We present here the first observations of this kind, performed on Rigel (B8Ia) and Deneb (A2Ia). Rigel was monitored by AMBER in two campaigns, in 2006-2007 and 2009-2010, and observed in 2009 by VEGA; whereas Deneb was monitored in 2008-2009 by VEGA. The extension of the Halpha and BrGamma line forming regions were accurately measured and compared with CMFGEN models of both stars. Moreover, clear signs of activity were observed in the differential visibility and phases. These pioneer observations are still limited, but show the path for a better understanding of the spatial structure and temporal evolution of localized ejections using optical interferometry.

preprint2010arXiv

Non-LTE time-dependent spectroscopic modelling of type II-plateau supernovae from the photospheric to the nebular phase: case study for 15 and 25Msun progenitor stars

We present the first non-LTE time-dependent radiative-transfer simulations of supernovae (SNe) II-Plateau (II-P) covering both the photospheric and nebular phases, from ~10 to >~1000d after the explosion, and based on 1.2B piston-driven ejecta produced from a 15Msun and a 25Msun non-rotating solar-metallicity star. The radial expansion of the gradually cooling photosphere gives rise to a near-constant luminosity up to >~100d after explosion. The photosphere remains in the outer 0.5Msun of the ejecta for up to ~50d after explosion. As the photosphere reaches the edge of the helium core, the SN luminosity drops by an amount mitigated by the progenitor radius and the 56Ni mass. Synthetic light-curves exhibit a bell-shape morphology, evolving faster for more compact progenitors, and with an earlier peak and narrower width in bluer filters. UV and U-band fluxes are very sensitive to line-blanketing, the metallicity, and the adopted model atoms. During the recombination epoch synthetic spectra are dominated by HI and metal lines, and are largely insensitive to the differing H/He/C/N/O composition of our two progenitor stars. In contrast, synthetic nebular-phase spectra reveal a broader/stronger OI doublet line in the higher-mass progenitor model, reflecting the larger masses of oxygen and nickel that are ejected. Our simulations overestimate the typical luminosity and the visual rise time of standard SNe II-P, likely a consequence of our progenitor stars being too big and/or too hydrogen rich. Comparison of our simulations with photospheric-phase observations of SN1999em of the same color are satisfactory. Our neglect of non-thermal processes leads to a fast disappearance of continuum radiation and Balmer-line emission at the end of the plateau phase. With the exception of HI lines, our nebular spectra show a striking similarity to contemporaneous observations of SN1999em.

preprint2010arXiv

Shock-heating of stellar envelopes: A possible common mechanism at the origin of explosions and eruptions in massive stars

Observations of transient phenomena in the Universe reveal a spectrum of mass-ejection properties associated with massive stars, covering from Type II/Ib/Ic core-collapse supernovae (SNe) to giant eruptions of Luminous Blue Variables (LBV) and optical transients. Here, we hypothesize that a fraction of these phenomena may have an explosive origin, the distinguishing ingredient being the ratio of the prompt energy release E_dep to the envelope binding energy E_binding. Using one-dimensional one-group radiation hydrodynamics and a set of 10-25Msun, massive-star models, we explore the dynamical response of a stellar envelope subject to a strong, sudden, and deeply-rooted energy release. Following energy deposition, a shock systematically forms, crosses the progenitor envelope on a day timescale, and breaks-out with a signal of hour-to-days duration and a 10^5-10^11 Lsun luminosity. For E_dep > E_binding, full envelope ejection results with a SN-like bolometric luminosity and kinetic energy, modulations being commensurate to the energy deposited and echoing the diversity of Type II-Plateau SNe. For E_dep ~ E_binding, partial envelope ejection results with a small expansion speed, and a more modest but year-long luminosity plateau, reminiscent of LBV eruptions or so-called SN impostors. For E_dep < E_binding, we obtain a "puffed-up" star, secularly relaxing back to thermal equilibrium. In parallel with gravitational collapse and Type II SNe, we argue that the thermonuclear combustion of merely a few 0.01Msun of C/O could power a wide range of explosions/eruptions in loosely-bound massive stars, as those in the 8-12Msun range, or in more massive ones owing to their proximity to the Eddington limit and/or critical rotation.

preprint2010arXiv

Supernova radiative-transfer modeling: A new approach using non-LTE and full time dependence

We discuss a new one-dimensional non-LTE time-dependent radiative-transfer technique for the simulation of supernova (SN) spectra and light curves. Starting from a hydrodynamical input characterizing the homologously-expanding ejecta at a chosen post-explosion time, we model the evolution of the entire ejecta, including gas and radiation. Non-LTE, which holds in all regions at and above the photosphere, is accounted for. The effects of line blanketing on the radiation field are explicitly included, using complex model atoms and solving for all ion level populations appearing in the statistical-equilibrium equations. Here, we present results for SN1987A, evolving the model "lm18a7Ad" of Woosley from 0.27 to 20.8d. The fastest evolution occurs prior to day 1, with a spectral energy distribution peaking in the range 300-2000A, subject to line blanketing from highly ionized metal and CNO species. After day 1, our synthetic multi-band light curve and spectra reproduce the observations to within 10-20% in flux in the optical, with a greater mismatch for the faint UV flux. We do not encounter any of the former discrepancies associated with the HeI and HI lines in the optical, which can be fitted well with a standard Blue-supergiant-star surface composition and no contribution from radioactive decay. The effects of time dependence on the ionization structure, discussed in Dessart & Hillier, are recovered, and thus nicely integrated in this new scheme. Despite the 1D nature of our approach, its high physical consistency and accuracy will allow reliable inferences to be made on explosion properties and pre-SN star evolution.

preprint2008arXiv

2D Multi-Angle, Multi-Group Neutrino Radiation-Hydrodynamic Simulations of Postbounce Supernova Cores

We perform axisymmetric (2D) multi-angle, multi-group neutrino radiation-hydrodynamic calculations of the postbounce phase of core-collapse supernovae using a genuinely 2D discrete-ordinate (S_n) method. We follow the long-term postbounce evolution of the cores of one nonrotating and one rapidly-rotating 20-solar-mass stellar model for ~400 milliseconds from 160 ms to ~550 ms after bounce. We present a multi-D analysis of the multi-angle neutrino radiation fields and compare in detail with counterpart simulations carried out in the 2D multi-group flux-limited diffusion (MGFLD) approximation to neutrino transport. We find that 2D multi-angle transport is superior in capturing the global and local radiation-field variations associated with rotation-induced and SASI-induced aspherical hydrodynamic configurations. In the rotating model, multi-angle transport predicts much larger asymptotic neutrino flux asymmetries with pole to equator ratios of up to ~2.5, while MGFLD tends to sphericize the radiation fields already in the optically semi-transparent postshock regions. Along the poles, the multi-angle calculation predicts a dramatic enhancement of the neutrino heating by up to a factor of 3, which alters the postbounce evolution and results in greater polar shock radii and an earlier onset of the initially rotationally weakened SASI. In the nonrotating model, differences between multi-angle and MGFLD calculations remain small at early times when the postshock region does not depart significantly from spherical symmetry. At later times, however, the growing SASI leads to large-scale asymmetries and the multi-angle calculation predicts up to 30% higher average integral neutrino energy deposition rates than MGFLD.

preprint2007arXiv

Simulations of Magnetically-Driven Supernova and Hypernova Explosions in the Context of Rapid Rotation

We present here the first 2D rotating, multi-group, radiation magnetohydrodynamics (RMHD) simulations of supernova core collapse, bounce, and explosion. In the context of rapid rotation, we focus on the dynamical effects of magnetic stresses and the creation and propagation of MHD jets. We find that a quasi-steady state can be quickly established after bounce, during which a well-collimated MHD jet is maintained by continuous pumping of power from the differentially rotating core. If the initial spin period of the progenitor core is $\sles$ 2 seconds, the free energy reservoir in the secularly evolving protoneutron star is adequate to power a supernova explosion, and may be enough for a hypernova. The jets are well collimated by the infalling material and magnetic hoop stresses, and maintain a small opening angle. We see evidence of sausage instabilities in the emerging jet stream. Neutrino heating is sub-dominant in the rapidly rotating models we explore, but can contribute 10$-$25% to the final explosion energy. Our simulations suggest that even in the case of modest or slow rotation, a supernova explosion might be followed by a secondary, weak MHD jet explosion, which, because of its weakness may to date have gone unnoticed in supernova debris. Furthermore, we suggest that the generation of a non-relativistic MHD precursor jet during the early protoneutron star/supernova phase is implicit in both the collapsar and "millisecond magnetar" models of GRBs. The multi-D, multi-group, rapidly rotating RMHD simulations we describe here are a start along the path towards more realistic simulations of the possible role of magnetic fields in some of Nature's most dramatic events.

preprint2006arXiv

Early Ultraviolet, Optical and X-Ray Observations of the Type IIP SN 2005cs in M51 with Swift

We report early photospheric-phase observations of the Type IIP Supernova (SN) 2005cs obtained by Swift's Ultraviolet-Optical and X-Ray Telescopes. Observations started within two days of discovery and continued regularly for three weeks. During this time the V-band magnitude remained essentially constant, while the UV was initially bright but steadily faded until below the brightness of an underlying UV-bright HII region. This UV decay is similar to SNe II observed by the International Ultraviolet Explorer. UV grism spectra show the P-Cygni absorption of MgII 2798A, indicating a photospheric origin of the UV flux. Based on non-LTE model atmosphere calculations with the CMFGEN code, we associate the rapid evolution of the UV flux with the cooling of the ejecta, the peak of the spectral energy distribution (SED) shifting from ~700A on June 30th to ~1200A on July 5th. Furthermore, the corresponding recombination of the ejecta, e.g., the transition from FeIII to FeII, induces considerable strengthening of metal line-blanketing at and above the photosphere, blocking more effectively this fading UV flux. SN2005cs was not detected in X-rays, and the upper limit to the X-ray luminosity yields a limit to the mass loss rate of the progenitor of about 10^-5 solar masses per year. Overall, Swift presents a unique opportunity to capture the early and fast evolution of Type II SNe in the UV, providing additional constraints on the reddening, the SED shortward of 4000A, and the ionization state and temperature of the photon-decoupling regions.