Nebular phase spectra of core-collapse supernovae (SNe) provide critical and unique information on the progenitor massive star and its explosion. We present a set of 1-D steady-state non-local thermodynamic equilibrium radiative transfer calculations of type II SNe at 300d after explosion. Guided by results for a large set of stellar evolution simulations, we craft ejecta models for type II SNe from the explosion of a 12, 15, 20, and 25Msun star. The ejecta density structure and kinetic energy, the 56Ni mass, and the level of chemical mixing are parametrized. Our model spectra are sensitive to the adopted line Doppler width, a phenomenon we associate with the overlap of FeII and OI lines with Lyalpha and Lybeta. Our spectra show a strong sensitivity to 56Ni mixing since it determines where decay power is absorbed. Even at 300d after explosion, the H-rich layers reprocess the radiation from the inner metal rich layers. In a given progenitor model, variations in 56Ni mass and distribution impact the ejecta ionization, which can modulate the strength of all lines. Such ionization shifts can quench CaII line emission. In our set of models, the OI6300 doublet strength is the most robust signature of progenitor mass. However, we emphasize that convective shell merging in the progenitor massive star interior can pollute the O-rich shell with Ca, which will weaken the OI6300 doublet flux in the resulting nebular SN II spectrum. This process may occur in Nature, with a greater occurrence in higher mass progenitors, and may explain in part the preponderance of progenitor masses below 17Msun inferred from nebular spectra.

We present observations of the unusually luminous Type II supernova (SN) 2016gsd. With a peak absolute magnitude of V = $-$19.95 $\pm$ 0.08, this object is one of the brightest Type II SNe, and lies in the gap of magnitudes between the majority of Type II SNe and the superluminous SNe. Its light curve shows little evidence of the expected drop from the optically thick phase to the radioactively powered tail. The velocities derived from the absorption in H$\alpha$ are also unusually high with the blue edge tracing the fastest moving gas initially at 20000 km s$^{-1}$, and then declining approximately linearly to 15000 km s$^{-1}$ over $\sim$100 d. The dwarf host galaxy of the SN indicates a low-metallicity progenitor which may also contribute to the weakness of the metal lines in its spectra. We examine SN 2016gsd with reference to similarly luminous, linear Type II SNe such as SNe 1979C and 1998S, and discuss the interpretation of its observational characteristics. We compare the observations with a model produced by the JEKYLL code and find that a massive star with a depleted and inflated hydrogen envelope struggles to reproduce the high luminosity and extreme linearity of SN 2016gsd. Instead, we suggest that the influence of interaction between the SN ejecta and circumstellar material can explain the majority of the observed properties of the SN. The high velocities and strong H$\alpha$ absorption present throughout the evolution of the SN may imply a circumstellar medium configured in an asymmetric geometry.

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors.

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]

A non-local-thermodynamic-equilibrium (NLTE) level population model of the first and second ionisation stages of iron, nickel and cobalt is used to fit a sample of XShooter optical + near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia). From the ratio of the NIR lines to the optical lines limits can be placed on the temperature and density of the emission region. We find a similar evolution of these parameters across our sample. Using the evolution of the Fe II 12$\,$570$\,\mathring{A}\,$to 7$\,$155$\,\mathring{A}\,$line as a prior in fits of spectra covering only the optical wavelengths we show that the 7200$\,\mathring{A}\,$feature is fully explained by [Fe II] and [Ni II] alone. This approach allows us to determine the abundance of Ni II$\,$/$\,$Fe II for a large sample of 130 optical spectra of 58 SNe Ia with uncertainties small enough to distinguish between Chandrasekhar mass (M$_{\text{Ch}}$) and sub-Chandrasekhar mass (sub-M$_{\text{Ch}}$) explosion models. We conclude that the majority (85$\%$) of normal SNe Ia have a Ni/Fe abundance that is in agreement with predictions of sub-M$_{\text{Ch}}$ explosion simulations of $\sim Z_\odot$ progenitors. Only a small fraction (11$\%$) of objects in the sample have a Ni/Fe abundance in agreement with M$_{\text{Ch}}$ explosion models.

Hydrogen-rich (type II) supernovae (SNe) exhibit considerable photometric and spectroscopic diversity. Extending previous work that focused exclusively on photometry, we simultaneously model the multi-band light curves and optical spectra of Type II SNe using RSG progenitors that are characterized by their H-rich envelope masses or the mass and extent of an enshrouding cocoon at the star's surface. Reducing the H-rich envelope mass yields faster declining light curves, a shorter duration of the photospheric phase, broader line profiles at early times, but only a modest boost in early-time optical brightness. Increasing the mass of the circumstellar material (CSM) is more effective at boosting the early-time brightness and producing a fast-declining light curve while leaving the duration of the photospheric phase intact. It also makes the optical color bluer, delays the onset of recombination, and can severely reduce the speed of the fastest ejecta material. The early ejecta interaction with CSM is conducive to producing featureless spectra at $10-20$ d and a weak or absent H$\alpha$ absorption during the recombination phase. The slow decliners SNe 1999em, 2012aw, and 2004et can be explained with a $1.2 \times 10^{51}$ erg explosion in a compact ($\sim$600 R$_\odot$) RSG star from a 15 M$_\odot$ stellar evolution model. A small amount of CSM ($<0.2$ M$_\odot$) improves the match to the SN photometry at times $<$ 10 d. With more extended RSG progenitors, one predicts lower ejecta kinetic energies, but the SN color stays blue for too long and the spectral line widths are too narrow. The fast decliners SNe 2013ej and 2014G may require $0.5-1.0$ M$_\odot$ of CSM, although this depends on the CSM structure. A larger boost to the luminosity (as for fast decliners SNe 1979C or 1998S) requires interaction with a more spatially extended CSM, which might also be detached from the star.

Core-collapse supernova (SN) ejecta are probably structured on both small and large scales, with greater deviations from spherical symmetry nearer the explosion site. Here, we present 2D and 3D gray radiation-hydrodynamics simulations of type II SN light curves from red (RSG) and blue supergiant (BSG) star explosions to investigate the impact on SN observables of inhomogeneities in density or composition, with a characteristic scale set to a few percent of the local radius. Clumping is found to hasten the release of stored radiation, boosting the early time luminosity and shortening the photospheric phase. Around the photosphere, radiation leaks between the clumps where the photon mean free path is greater. Since radiation is stored uniformly in volume, a greater clumping can increase this leakage by storing more and more mass into smaller and denser clumps containing less and less radiation energy. An inhomogeneous medium in which different regions recombine at different temperatures can also impact the light curve. Clumping can thus be a source of diversity in SN brightness. Clumping may lead to a systematic underestimate of ejecta masses from light curve modeling, although a significant offset seems to require a large density contrast of a few tens between clumps and interclump medium.

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] $\lambda\lambda$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] $\lambda\lambda$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.

We present a new relativistic radiative-transfer code for $\gamma$-rays of energy less than 5 MeV in supernova (SN) ejecta. This code computes the opacities, the prompt emissivity (i.e. decay), and the scattering emissivity, and solves for the intensity in the co-moving frame. Because of the large expansion velocities of SN ejecta, we ignore redistribution effects associated with thermal motions. The energy deposition is calculated from the energy removed from the radiation field by scattering or photoelectric absorption. This new code yields comparable results to an independent Monte Carlo code. However, both yield non-trivial differences with the results from a pure absorption treatment of $\gamma$-ray transport. A synthetic observer's frame spectrum is also produced from the CMF intensity. At early times when the optical depth to $\gamma$-rays is large, the synthetic spectrum show asymmetric line profiles with redshifted absorption as seen in SN 2014J. This new code is integrated within CMFGEN and allows for an accurate and fast computation of the decay energy deposition in SN ejecta.

Much controversy surrounds the inferred progenitor masses of Type II-Plateau (II-P) supernovae (SNe). The debate is nourished by the discrepant results from radiation-hydrodynamics simulations, from pre-explosion imaging, and from studies of host stellar populations. Here, we present a controlled experiment using four solar metallicity models with zero-age main-sequence masses of 12, 15, 20, and 25Msun. Because of the effects of core burning and surface mass loss, these models reach core collapse as red-supergiant (RSG) stars with a similar H-rich envelope mass of 8 to 9Msun but with final masses in the range 11 to 16Msun. We explode the progenitors using a thermal bomb, adjusting the energy deposition to yield an asymptotic ejecta kinetic energy of 1.25 x 10^51 erg and an initial 56Ni mass of 0.04Msun. The resulting SNe produce similar photometric and spectroscopic properties from 10 to 200d. The spectral characteristics are degenerate. The scatter in early-time color results from the range in progenitor radii, while the differences in late-time spectra reflect the larger oxygen yields in more massive progenitors. Because the progenitors have a comparable H-rich envelope mass, the photospheric phase duration is comparable for all models; the difference in He-core mass is invisible. As different main-sequence masses can produce progenitors with a similar H-rich envelope mass, light curve modeling cannot provide a robust and unique solution for the ejecta mass of Type II-P SNe. The numerous uncertainties in massive star evolution and wind mass loss also prevent a robust association with a main-sequence star mass. Light curve modeling can at best propose compatibility.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Since the discovery of SN (supernova) 1987A, the number of Type II-peculiar SNe has grown, revealing a rich diversity in photometric and spectroscopic properties. In this study, using a single 15Msun low-metallicity progenitor that dies as a blue supergiant (BSG), we have generated explosions with a range of energies and 56Ni masses. We then performed the radiative transfer modeling with CMFGEN from 1d until 300d after explosion. Our models yield light curves that rise to optical maximum in ~100d, with a similar brightening rate, and with a peak absolute V-band magnitude spanning from -14 to -16.5mag. All models follow a similar color evolution, entering the recombination phase within a few days of explosion, and reddening further until the nebular phase. Their spectral evolution is analogous, mostly differing in line profile width. With this model set, we study the Type II-pec SNe 1987A, 2000cb, 2006V, 2006au, 2009E, and 2009mw. Their photometric and spectroscopic diversity suggest that there is no prototypical Type II-pec SN. These SNe brighten to maximum faster than our model set, except perhaps SN2009mw. The spectral evolution of SN1987A conflicts with other observations and with model predictions from 20d until maximum: Halpha narrows and weakens while BaII lines strengthen faster than predicted, which we interpret as signatures of clumping. SN2000cb rises to maximum in only 20d and shows weak BaII lines. Its spectral evolution is well matched by an energetic ejecta but the light curve may require asymmetry. The persistent blue color, narrow lines, and weak Halpha absorption, seen in SN2006V conflicts with expectations for a BSG explosion powered by 56Ni and may require an alternative power source. In addition to diversity arising from different BSG progenitors, we surmise that their ejecta are asymmetric, clumped, and, in some cases, not solely powered by 56Ni decay [abridged].

Numerous superluminous supernovae (SLSNe) of Type Ic have been discovered and monitored in the last decade. The favored mechanism at their origin is a sustained power injection from a magnetar. This study presents non-local thermodynamic equilibrium time-dependent radiative transfer simulations of various single carbon-rich Wolf-Rayet star explosions influenced by magnetars of diverse properties and covering from a few days to one or two years after explosion. Nonthermal processes are treated; the magnetar-power deposition profile is prescribed; dynamical effects are ignored. In this context, the main influence of the magnetar power is to boost the internal energy of the ejecta on week-long time scales, enhancing the ejecta temperature and ionization, shifting the spectral energy distribution to the near-UV (even for the adopted solar metallicity), creating blue optical colors. Varying the ejecta and magnetar properties introduces various stretches and shifts to the light curve (rise time, peak or nebular luminosity, light curve width). At maximum, all models show the presence of OII and CII lines in the optical, and more rarely OIII and CIII lines. Non-thermal effects are found to be negligible during the high-brightness phase. After maximum, higher energy explosions are hotter and more ionized, and produce spectra that are optically bluer. Clumping is a source of spectral diversity after maximum. Clumping is essential to trigger ejecta recombination and yield the presence of OI, CaII, and FeII lines from a few weeks after maximum until nebular times. The UV and optical spectrum of Gaia16apd at maximum or the nebular spectrum of LSQ14an at +410d are compatible with some models that assume no clumping. However, most observed SLSNe Ic seem to require clumping from early post-maximum to nebular times (e.g., SN2007bi at +46 and +367d; Gaia16apd at +43d).

In the first weeks-to-months of a type II-P supernova (SN), the spectrum formation region is within the hydrogen-rich envelope of the exploding star. Optical spectra taken within a few days of the SN explosion, when the photosphere is hot, show features of ionised carbon, nitrogen and oxygen, as well as hydrogen and helium. Quantitative analysis of this very early phase may therefore constrain the chemical abundances of the stellar envelope at the point of core-collapse. Using existing and new evolutionary calculations for Red Supergiants (RSGs), we show that the predictions for the terminal surface [C/N] ratio is correlated with the initial mass of the progenitor star. Specifically, a star with an initial mass above 20M$_{\odot}$ exploding in the RSG phase should have an unequivocal signal of a low [C/N] abundance. Furthermore, we show that the model predictions are relatively insensitive to uncertainties in the treatment of convective mixing. Although there is a dependence on initial rotation, this can be dealt with in a probabilistic sense by convolving the model predictions with the observed distribution of stellar rotation rates. Using numerical experiments, we present a strategy for using very early-time spectroscopy to determine the upper limit to the progenitor mass distribution for type II-P SNe.

We investigate the effect of mixing of radioactive nickel ($^{56}$Ni) on the early-time color evolution of Type Ib and Ic supernovae (SNe Ib/Ic) using multi-group radiation hydrodynamics simulations. We consider both helium-rich and helium-poor progenitors. Mixing of $^{56}$Ni is parameterized using a Gaussian distribution function. We find that the early-time color evolution with a weak $^{56}$Ni mixing is characterized by three different phases: initial rapid reddening, blueward evolution due to the delayed effect of $^{56}$Ni heating, and redward evolution thereafter until the transition to the nebular phase. With a strong $^{56}$Ni mixing, the second phase disappears. We compare our models with the early-time color evolution of several SNe Ib/Ic (SN1999ex, SN 2008D, SN 2009jf, iPTF13bvn, SN 1994I, SN 2007gr, SN 2013ge, and 2017ein) and find signatures of relatively weak and strong $^{56}$Ni mixing for SNe Ib and SNe Ic, respectively. This suggests that SNe Ib progenitors are distinct from SN Ic progenitors in terms of helium content and that $^{56}$Ni mixing is generally stronger in the carbon-oxygen core and weaker in the helium-rich envelope. We conclude that the early-time color evolution is a powerful probe of $^{56}$Ni mixing in SNe Ib/Ic.

Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.

There is both observational and theoretical evidence that the ejecta of core-collapse supernovae (SNe) are structured. Rather than being smooth and homogeneous, the material is made of over-dense and under-dense regions of distinct composition. Here, we explore the effect of clumping on the SN radiation during the photospheric phase using 1-D non-local thermodynamic equilibrium radiative transfer and an ejecta model arising from a blue-supergiant explosion (yielding a Type II-peculiar SN). Neglecting chemical segregation, we adopt a velocity-dependent volume-filling factor approach that assumes that the clumps are small but does not change the column density along any sightline. We find that clumping boosts the recombination rate in the photospheric layers, leading to a faster recession of the photosphere, an increase in bolometric luminosity, and a reddening of the SN colors through enhanced blanketing. The SN bolometric light curve peaks earlier and transitions faster to the nebular phase. On the rise to maximum, the strongest luminosity contrast between our clumped and smooth models is obtained at the epoch when the photosphere has receded to ejecta layers where the clumping factor is only 0.5 -- this clumping factor may be larger in Nature. Clumping is seen to have a similar influence in a Type II-Plateau SN model. As we neglect both porosity and chemical segregation our models underestimate the true impact of clumping. These results warrant further study of the influence of clumping on the observables of other SN types during the photospheric phase.

The origin of dust in the early Universe has been the subject of considerable debate. Core-collapse supernovae (ccSNe), which occur several million years after their massive progenitors form, could be a major source of that dust, as in the local universe several ccSNe have been observed to be copious dust producers. Here we report nine near-infrared (0.8 - 2.5 micron) spectra of the Type II-P SN 2017eaw in NGC 6946, spanning the time interval 22 - 205 days after discovery. The specta show the onset of CO formation and continuum emission at wavelengths greater than 2.1 micron from newly-formed hot dust, in addition to numerous lines of hydrogen and metals, which reveal the change in ionization as the density of much of the ejecta decreases. The observed CO masses estimated from an LTE model are typically 10^{-4} Msun during days 124 - 205, but could be an order of magnitude larger if non-LTE conditions are present in the emitting region. The timing of the appearance of CO is remarkably consistent with chemically controlled dust models of Sarangi & Cherchneff.

We present an analysis of a new sample of type II core-collapse supernovae (SNe II) occurring within low-luminosity galaxies, comparing these with a sample of events in brighter hosts. Our analysis is performed comparing SN II spectral and photometric parameters and estimating the influence of metallicity (inferred from host luminosity differences) on SN II transient properties. We measure the SN absolute magnitude at maximum, the light-curve plateau duration, the optically thick duration, and the plateau decline rate in the V-band, together with expansion velocities and pseudo-equivalent-widths (pEWs) of several absorption lines in the SN spectra. For the SN host galaxies, we estimate the absolute magnitude and the stellar mass, a proxy for the metallicity of the host galaxy. SNe II exploding in low luminosity galaxies display weaker pEWs of Fe II $\lambda5018$, confirming the theoretical prediction that metal lines in SN II spectra should correlate with metallicity. We also find that SNe II in low-luminosity hosts have generally slower declining light curves and display weaker absorption lines. We find no relationship between the plateau duration or the expansion velocities with SN environment, suggesting that the hydrogen envelope mass and the explosion energy are not correlated with the metallicity of the host galaxy. This result supports recent predictions that mass-loss for red supergiants is independent of metallicity.

Red supergiants have been confirmed as the progenitor stars of the majority of hydrogen-rich type II supernovae. However, while such stars are observed with masses >25M$_\odot$, detections of >18M$_\odot$ progenitors remain elusive. Red supergiants are also expected to form at all metallicities, but discoveries of explosions from low-metallicity progenitors are scarce. Here, we report observations of the type II supernova, SN 2015bs, for which we infer a progenitor metallicity of $\leq$0.1Z$_\odot$ from comparison to photospheric-phase spectral models, and a Zero Age Main-Sequence mass of 17-25M$_\odot$ through comparison to nebular-phase spectral models. SN 2015bs displays a normal 'plateau' light-curve morphology, and typical spectral properties, implying a red supergiant progenitor. This is the first example of such a high mass progenitor for a 'normal' type II supernova, suggesting a link between high mass red supergiant explosions and low-metallicity progenitors.

Magnetar power is believed to be at the origin of numerous super-luminous supernovae (SNe) of Type Ic, arising from compact, hydrogen-deficient, Wolf-Rayet type stars. Here, we investigate the properties that magnetar power would have on standard-energy SNe associated with 15-20Msun blue or red supergiant (BSG/RSG) stars. We use a combination of Eulerian grey radiation-hydrodynamics and non-LTE steady-state radiative transfer to study their dynamical, photometric, and spectroscopic properties. Adopting magnetar fields of 1, 3.5, 7 x 10^14G and rotational energies of 0.4, 1, and 3 x 10^51erg, we produce bolometric light curves with a broad maximum covering 50-150d and a magnitude of 10^43-10^44erg/s. The spectra at maximum light are analogous to those of standard SNe II-P but bluer. Although the magnetar energy is channelled roughly in equal proportion between SN kinetic energy and SN luminosity, the latter may be boosted by a factor 10-100 compared to a standard SN II. This influence breaks the observed relation between brightness and expansion rate of standard Type II SNe. Magnetar energy injection also delays recombination and may even cause re-ionization, with a reversal in photospheric temperature and velocity. Depositing the magnetar energy in a narrow mass shell at the ejecta base leads to the formation of a dense shell at a few 1000km/s, which causes a light-curve bump at the end of the photospheric phase. Depositing this energy over a broad range of mass in the inner ejecta, to mimic the effect of multi-dimensional fluid instabilities, prevents the formation of a dense shell and produces an earlier-rising and smoother light curve. The magnetar influence on the SN radiation is generally not visible prior to 20-30d, during which one may discern a BSG from a RSG progenitor. We propose a magnetar model for the super-luminous Type II SN OGLE-SN14-073.

We present a study of observed Type II supernova (SN~II) colours using optical/near-infrared photometric data from the \textit{Carnegie Supernovae Project-I}. We analyse four colours ($B-V$, $u-g$, $g-r$, and $g-Y$) and find that SN~II colour curves can be described by two linear regimes during the photospheric phase. The first ($s_{\rm 1,colour}$) is steeper and has a median duration of $\sim 40$ days. The second, shallower slope ($s_{\rm 2,colour}$) lasts until the end of the "plateau" ($\sim 80$ days). The two slopes correlate in the sense that steeper initial colour curves also imply steeper colour curves at later phases. As suggested by recent studies, SNe~II form a continuous population of objects from the colour point of view as well. We investigate correlations between the observed colours and a range of photometric and spectroscopic parameters including the absolute magnitude, the $V$-band light-curve slopes, and metal-line strengths. We find that less luminous SNe~II appear redder, a trend that we argue is not driven by uncorrected host-galaxy reddening. While there is significant dispersion, we find evidence that redder SNe~II (mainly at early epochs) display stronger metal-line equivalent widths. Host-galaxy reddening does not appear to be a dominant parameter, neither driving observed trends nor dominating the dispersion in observed colours. Intrinsic SN~II colours are most probably dominated by photospheric temperature differences, with progenitor metallicity possibly playing a minor role. Such temperature differences could be related to differences in progenitor radius, together with the presence or absence of circumstellar material close to the progenitor stars.

Transient surveys have recently revealed the existence of H-rich super-luminous supernovae (SLSN; e.g., iPTF14hls, OGLE-SN14-073) characterized by an exceptionally large time-integrated bolometric luminosity, a sustained blue optical color, and Doppler-broadened HI lines at all times. Here, I investigate the effect that a magnetar (initial rotational energy of 4x10^50erg and field strength of 7x10^13G) would have on the properties of a typical Type II SN ejecta (mass of 13.35Msun, kinetic energy of 1.32x10^51erg, 0.077Msun of 56Ni) produced by the terminal explosion of an H-rich blue-supergiant star. I present a non-LTE time-dependent radiative transfer simulation of the resulting photometric and spectroscopic evolution from 1d until 600d after explosion. With magnetar power, the model luminosity and brightness are enhanced, the ejecta is everywhere hotter and more ionised, and the spectrum formation region is much more extended. This magnetar-powered SN ejecta reproduces most of the observed properties of SLSN iPTF14hls, including the sustained brightness of -18mag in the R band, the blue optical color, and the broad HI lines for 600d. The non-extreme magnetar properties, combined with the standard Type II SN ejecta properties offer an interesting alternative to the pair-unstable super-massive star model recently proposed, which involves a highly-energetic and super-massive ejecta. Hence, such Type II SLSNe may differ from standard Type II SNe exclusively through the influence of a magnetar.

Massive stars with a core-halo structure are interesting objects for stellar physics and hydrodynamics. Using simulations for stellar evolution, radiation hydrodynamics, and radiative transfer, we study the explosion of stars with an extended and tenuous envelope (i.e., stars in which 95% of the mass is contained within 10% of the surface radius or less). We consider both H-rich supergiant and He-giant progenitors resulting from close-binary evolution and dying with a final mass of 2.8-5Msun. An extended envelope causes the supernova (SN) shock to brake and a reverse shock to form, sweeping core material into a dense shell. The shock deposited energy, which suffers little degradation from expansion, is trapped in ejecta layers of moderate optical depth, thereby enhancing the SN luminosity at early times. With the delayed 56Ni heating, we find that the resulting optical and near-IR light curves all exhibit a double-peak morphology. We show how an extended progenitor can explain the blue and featureless optical spectra of some Type IIb and Ib SNe. The dense shell formed by the reverse shock leads to line profiles with a smaller and near-constant width. This ejecta property can explain the statistically narrower profiles of Type IIb compared to Type Ib SNe, as well as the peculiar Halpha profile seen in SN1993J. At early times, our He-giant star explosion model shows a high luminosity, a blue color, and featureless spectra reminiscent of the Type Ib SN2008D, suggesting a low-mass progenitor.

While Chandrasekhar-mass ($M_{\rm Ch}$) models with a low $^{56}\rm{Ni}$ yield can match the peak luminosities of fast-declining, 91bg-like Type Ia supernovae (SNe Ia), they systematically fail to reproduce their faster light-curve evolution. Here we illustrate the impact of a low ejecta mass on the radiative display of low-luminosity SNe Ia, by comparing a sub-$M_{\rm Ch}$ model resulting from the pure central detonation of a C-O White Dwarf (WD) to a $M_{\rm Ch}$ delayed-detonation model with the same $^{56}\rm{Ni}$ yield of 0.12 M$_{\odot}$. Our sub-$M_{\rm Ch}$ model from a 0.90 M$_{\odot}$ WD progenitor has a $\sim$5 day shorter rise time in the integrated UV-optical-IR (uvoir) luminosity, as well as in the $B$-band, and a $\sim$20 per cent higher peak uvoir luminosity ($\sim$1 mag brighter peak $M_B$). This sub-$M_{\rm Ch}$ model also displays bluer maximum-light colours due to the larger specific heating rate, and larger post-maximum uvoir and $B$-band decline rates. The luminosity decline at nebular times is also more pronounced, reflecting the enhanced escape of gamma rays resulting from the lower density of the progenitor WD. The deficit of stable nickel in the innermost ejecta leads to a notable absence of forbidden lines of [Ni II] in the nebular spectra. In contrast, the $M_{\rm Ch}$ model displays a strong line due to [Ni II] 1.939 $\mu\rm{m}$, which could in principle serve to distinguish between different progenitor scenarios. Our sub-$M_{\rm Ch}$ model offers an unprecedented agreement with optical and near-infrared observations of the 91bg-like SN 1999by, making a strong case for a WD progenitor significantly below the Chandrasekhar-mass limit for this event and other low-luminosity SNe Ia.

We present one-dimensional non-local thermodynamic equilibrium time-dependent radiative transfer simulations (using CMFGEN) of two sub-Chandrasekhar (sub-$M_{\rm Ch}$), one $M_{\rm Ch}$ and one super-$M_{\rm Ch}$ Type Ia SN ejecta models. Three originate from $M_{\rm Ch}$ delayed detonation models, and the fourth is a sub-$M_{\rm Ch}$ detonation model. Ejecta masses are 1.02, 1.04, 1.40, and 1.70 M$_\odot$, and all models have 0.62 M$_\odot$ of $^{56}{\rm Ni}$. Sub-$M_{\rm Ch}$ model light curves evolve faster, reaching bolometric maximum 2--3 days earlier and having 3--4 days shorter bolometric half light widths. The models vary by $\sim$12 per cent at maximum bolometric luminosity and by 0.17 mag in $B_{\rm max}$. While $\Delta M_{15}(B)$ increases with ejecta mass it only varies by $\sim$5 per cent around 1 mag. Sub-$M_{\rm Ch}$ models are 0.25 mag bluer in $B-R$ at $B_{\rm max}$. Optical spectra share many similarities, but lower mass models exhibit less UV line blanketing during the photospheric phase. At nebular times, significant NIR spectroscopic differences are seen. In particular, emission lines of the Ca II NIR triplet; [S III] $\lambda\lambda$9068,9530; [Ca II] $\lambda\lambda$7291,7324; [Ar III] $\lambda\lambda$7135,7751; and [Ni II] 1.939 $\mu$m are stronger in higher mass models. The [Ni II] 1.939 $\mu$m line is absent in the sub-$M_{\rm Ch}$ detonation model, and provides a valuable potential tool to distinguish sub-$M_{\rm Ch}$ explosions from $M_{\rm Ch}$ explosions. In general, the nebular phase models are too highly ionized. We attribute this to the neglect of clumping and/or the distribution of intermediate mass and iron group elements. The two sub-$M_{\rm Ch}$ models, while exploded by different mechanisms, can be distinguished in the $J$ and $H$ bands at late times (e.g., $+200$ days).

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present our study of OGLE-2014-SN-073, one of the brightest Type II SN ever discovered, with an unusually broad lightcurve combined with high ejecta velocities. From our hydrodynamical modelling we infer a remarkable ejecta mass of $60^{+42}_{-16}$~M$_\odot$, and a relatively high explosion energy of $12.4^{+13.0}_{-5.9} \times10^{51}$~erg. We show that this object belongs, with a very small number of other hydrogen-rich SNe, to an energy regime that is not explained by standard core-collapse (CC) neutrino-driven explosions. We compare the quantities inferred by the hydrodynamical modelling with the expectations of various exploding scenarios, trying to explain the high energy and luminosity released. We find some qualitative similarities with pair-instabilities SNe, although a prompt injection of energy by a magnetar seems also a viable alternative to explain such extreme event.

The progenitors of low-luminosity Type II-Plateau supernovae (SNe II-P) are believed to be red supergiant (RSG) stars, but there is much disparity in the literature concerning their mass at core collapse and therefore on the main sequence. Here, we model the SN radiation arising from the low-energy explosion of RSG stars of 12, 25, and 27 M$_{\odot}$ on the main sequence and formed through single star evolution. Despite the narrow range in ejecta kinetic energy (2.5$-$4.2$\times$10$^{50}$ erg) in our model set, the SN observables from our three models are significantly distinct, reflecting the differences in progenitor structure (e.g., surface radius, H-rich envelope mass, He-core mass). Our higher mass RSG stars give rise to Type II SNe that tend to have bluer colors at early times, a shorter photospheric phase, and a faster declining $V$-band light curve (LC) more typical of Type II-linear SNe, in conflict with the LC plateau observed for low-luminosity SNe II. The complete fallback of the CO core in the low-energy explosions of our high mass RSG stars prevents the ejection of any ${}^{56}$Ni (nor any core O or Si), in contrast to low-luminosity SNe II-P, which eject at least 0.001 M$_{\odot}$ of ${}^{56}$Ni. In contrast to observations, type II SN models from higher mass RSGs tend to show an H$\alpha$ absorption that remains broad at late times (due to a larger velocity at the base of the H-rich envelope). In agreement with the analyses of pre-explosion photometry, we conclude that low-luminosity SNe II-P likely arise from low-mass rather than high-mass RSG stars.

We present an analysis of observed trends and correlations between a large range of spectral and photometric parameters of more than 100 type II supernovae (SNe~II), during the photospheric phase. We define a common epoch for all SNe of 50 days post-explosion where the majority of the sample is likely to be under similar physical conditions. Several correlation matrices are produced to search for interesting trends between more than 30 distinct light-curve and spectral properties that characterize the diversity of SNe~II. Overall, SNe with higher expansion velocities are brighter, have more rapidly declining light-curves, shorter plateau durations, and higher $^{56}$Ni masses. Using a larger sample than previous studies, we argue that `$Pd$' - the plateau duration from the transition of the initial to `plateau' decline rates to the end of the `plateau' - is a better indicator of the hydrogen envelope mass than the traditionally used optically thick phase duration ($OPTd$: explosion epoch to end of plateau). This argument is supported by the fact that $Pd$ also correlates with s$_3$, the light-curve decline rate at late times: lower $Pd$ values correlate with larger s$_3$ decline rates. Large s$_3$ decline rates are likely related to lower envelope masses that enables gamma-ray escape. We also find a significant anticorrelation between $Pd$ and s$_2$ (the plateau decline rate), confirming the long standing hypothesis that faster declining SNe~II (SNe~IIL) are the result of explosions with lower hydrogen envelope masses and therefore have shorter $Pd$ values.

We present 888 visual-wavelength spectra of 122 nearby type II supernovae (SNe II) obtained between 1986 and 2009, and ranging between 3 and 363 days post explosion. In this first paper, we outline our observations and data reduction techniques, together with a characterization based on the spectral diversity of SNe~II. A statistical analysis of the spectral matching technique is discussed as an alternative to non-detection constraints for estimating SN explosion epochs. The time evolution of spectral lines is presented and analysed in terms of how this differs for SNe of different photometric, spectral, and environmental properties: velocities, pseudo-equivalent widths, decline rates, magnitudes, time durations, and environment metallicity. Our sample displays a large range in ejecta expansion velocities, from $\sim9600$ to $\sim1500$ km s$^{-1}$ at 50 days post explosion with a median H$_{\alpha}$ value of 7300 km s$^{-1}$. This is most likely explained through differing explosion energies. Significant diversity is also observed in the absolute strength of spectral lines, characterised through their pseudo-equivalent widths. This implies significant diversity in both temperature evolution (linked to progenitor radius) and progenitor metallicity between different SNe~II. Around 60\% of our sample show an extra absorption component on the blue side of the H$_{\alpha}$ P-Cygni profile ("Cachito" feature) between 7 and 120 days since explosion. Studying the nature of Cachito, we conclude that these features at early times (before $\sim35$ days) are associated with \ion{Si}{2} $\lambda6355$, while past the middle of the plateau phase they are related to high velocity (HV) features of hydrogen lines.

The faster light-curve evolution of low-luminosity Type Ia supernovae (SNe Ia) suggests that they could result from the explosion of white dwarf (WD) progenitors below the Chandrasekhar mass ($M_{\rm Ch}$). Here we present 1D non-LTE time-dependent radiative transfer simulations of pure central detonations of carbon-oxygen WDs with a mass ($M_\rm{tot}$) between 0.88 M$_{\odot}$ and 1.15 M$_{\odot}$, and a $^{56}\rm{Ni}$ yield between 0.08 M$_{\odot}$ and 0.84 M$_{\odot}$. Their lower ejecta density compared to $M_{\rm Ch}$ models results in a more rapid increase of the luminosity at early times and an enhanced $\gamma$-ray escape fraction past maximum light. Consequently, their bolometric light curves display shorter rise times and larger post-maximum decline rates. Moreover, the higher $M(^{56}\rm{Ni})/M_\rm{tot}$ ratio at a given $^{56}\rm{Ni}$ mass enhances the temperature and ionization level in the spectrum-formation region for the less luminous models, giving rise to bluer colours at maximum light and a faster post-maximum evolution of the $B-V$ colour. For sub-$M_{\rm Ch}$ models fainter than $M_B\approx -18.5$ mag at peak, the greater bolometric decline and faster colour evolution lead to a larger $B$-band post-maximum decline rate, $\Delta M_{15}(B)$. In particular, all of our previously-published $M_{\rm Ch}$ models (standard and pulsational delayed detonations) are confined to $\Delta M_{15}(B) < 1.4$ mag, while the sub-$M_{\rm Ch}$ models with $M_\rm{tot}\lesssim 1$ M$_{\odot}$ extend beyond this limit to $\Delta M_{15}(B)\approx 1.65$ mag for a peak $M_B\approx -17$ mag, in better agreement with the observed width-luminosity relation (WLR). Regardless of the precise ignition mechanism, these simulations suggest that fast-declining SNe Ia at the faint end of the WLR could result from the explosion of WDs whose mass is significantly below the Chandrasekhar limit.

Early-time observations of the Type II supernovae (SNe) 2013cu and 2013fs have revealed an interaction of ejecta with material near the star surface. Unlike the Type IIn SN2010jl, which interacts with a dense wind for ~1yr, the interaction ebbs after 2-3d, suggesting a dense and compact circumstellar envelope. Here, we use multi-group radiation-hydrodynamics and non-local-thermodynamic-equilibrium radiative transfer to explore the properties of red supergiant (RSG) star explosions embedded in a variety of dense envelopes. We consider the cases of an extended static atmosphere or a steady-state wind, adopting a range of mass loss rates. The shock-breakout signal, the SN radiation up to 10d, and the ejecta dynamics are strongly influenced by the properties of this nearby environment. This compromises the use of early-time observations to constrain Rstar. The presence of narrow lines for 2-3d in 2013fs and 2013cu require a cocoon of material of ~0.01Msun out to 5-10Rstar. Spectral lines evolve from electron-scattering to Doppler broadened, with a growing blueshift of their emission peaks. Recent studies propose a super-wind phase with a mass loss rate from 0.001 up to 1Msun/yr in the last months/years of the RSG life, although there is no observational constraint that this external material is a steady-state outflow. Alternatively, observations may be explained by the explosion of a RSG star inside its complex atmosphere. Indeed, spatially resolved observations reveal that RSG stars have extended atmospheres, with the presence of downflows and upflows out to several Rstar, even in a standard RSG like Betelgeuse. Mass loading in the region intermediate between star and wind can accommodate the 0.01Msun needed to explain the observations of 2013fs. A puzzling super-wind phase prior to core-collapse may be therefore superfluous [Abridged].

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.

Colliding stellar winds in massive binary systems have been studied through their radio, optical lines and strong X-ray emission for decades. More recently, near-infrared spectrointerferometric observations have become available in a few systems, but isolating the contribution from the individual stars and the wind collision region still remains a challenge. In this paper, we study the colliding wind binary gamma2 Velorum and aim at identifying the wind collision zone from infrared interferometric data, which provide unique spatial information to determine the wind properties. Our analysis is based on multi-epoch VLTI/AMBER data that allows us to separate the spectral components of both stars. First, we determine the astrometric solution of the binary and confirm previous distance measurements. We then analyse the spectra of the individual stars, showing that the O star spectrum is peculiar within its class. Then, we perform three-dimensional hydrodynamic simulations of the system from which we extract model images, visibility curves and closure phases which can be directly compared with the observed data. The hydrodynamic simulations reveal the 3D spiral structure of the wind collision region, which results in phase-dependent emission maps. Our model visibility curves and closure phases provide a good match when the wind collision region accounts for 3 to 10 per cent of the total flux in the near infrared. The dialogue between hydrodynamic simulations, radiative transfer models and observations allows us to fully exploit the observations. Similar efforts will be crucial to study circumstellar environments with the new generation of VLTI instruments like GRAVITY and MATISSE.

We explore properties of Type Ib and IIb SN progenitors that are produced by stable mass transfer in binary systems, using a new grid of stellar evolution models from an initial primary mass in the range of 10 - 18 $\mathrm{M_\odot}$ at solar and Large Magellanic Cloud metallicities. We find that blue and yellow supergiant SN IIb progenitors (e.g., of SN 2008ax, SN 2011dh, SN 2016gkg) have a hydrogen envelope mass less than about 0.15 $\mathrm{M_\odot}$, mostly resulting from early Case B mass transfer with relatively low initial masses and/or low metallicity. Red supergiant (RSG) SN IIb progenitors (e.g., of SN 1993J, SN 2013df) are produced via late Case B mass transfer and have a more massive hydrogen envelope ($M_\mathrm{H,env} > 0.15~\mathrm{M_\odot}$). SN Ib progenitors are predominantly produced by early Case B mass transfer. Our models predict that SN IIb progenitors are systematically more luminous in the optical ($-8.0 \lesssim M_\mathrm{V} \lesssim -5.0$) than the majority of SN Ib progenitors ($M_\mathrm{V} \gtrsim -5.0$) for our considered initial mass range. However, the optically bright progenitor of SN Ib iPTF13bvn (i.e., $M_\mathrm{V} \simeq -6.5$) can be well explained by a relatively low-mass progenitor with a final mass of $\sim 3.0~\mathrm{M_\odot}$. The event rate of blue and yellow SN IIb progenitors would increase as metallicity decreases, while the event rate of SN Ib progenitors would decrease instead. By contrast, the population of RSG SN IIb progenitors would not be significantly affected by metallicity.

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development.

Supernova (SN) 2008bk is a well observed low-luminosity Type II event visually associated with a low-mass red-supergiant progenitor. To model SN 2008bk, we evolve a 12Msun star from the main sequence until core collapse, when it has a total mass of 9.88Msun, a He-core mass of 3.22Msun, and a radius of 502Rsun. We then artificially trigger an explosion that produces 8.29Msun of ejecta with a total energy of 2.5x10^50erg and ~0.009Msun of 56Ni. We model the subsequent evolution of the ejecta with non-Local-Thermodynamic-Equilibrium time-dependent radiative transfer. Although somewhat too luminous and energetic, this model reproduces satisfactorily the multi-band light curves and multi-epoch spectra of SN 2008bk, confirming the suitability of a low-mass massive star progenitor. As in other low-luminosity SNe II, the structured Halpha profile at the end of the plateau phase is probably caused by BaII 6496.9A rather than asphericity. We discuss the sensitivity of our results to changes in progenitor radius and mass, as well as chemical mixing. A 15% increase in progenitor radius causes a 15% increase in luminosity and a 0.2mag V-band brightening of the plateau but leaves its length unaffected. An increase in ejecta mass by 10% lengthens the plateau by ~10d. Chemical mixing introduces slight changes to the bolometric light curve, limited to the end of the plateau, but has a large impact on colours and spectra at nebular times.

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.

We present optical imaging and spectroscopy of supernova (SN) LSQ13fn, a type II supernova with several hitherto-unseen properties. Although it initially showed strong symmetric spectral emission features attributable to \ion{He}{ii}, \ion{N}{iii}, and \ion{C}{iii}, reminiscent of some interacting SNe, it transitioned into an object that would fall more naturally under a type II-Plateau (IIP) classification. However, its spectral evolution revealed several unusual properties: metal lines appeared later than expected, were weak, and some species were conspicuous by their absence. Furthermore, the line velocities were found to be lower than expected given the plateau brightness, breaking the SNe~IIP standardised candle method for distance estimates. We found that, in combination with a short phase of early-time ejecta-circumstellar material interaction, metal-poor ejecta, and a large progenitor radius could reasonably account for the observed behaviour. Comparisons with synthetic model spectra of SNe~IIP of a given progenitor mass would imply a progenitor star metallicity as low as 0.1\,Z$_{\odot}$. LSQ13fn highlights the diversity of SNe~II and the many competing physical effects that come into play towards the final stages of massive star evolution immediately preceding core-collapse.

Spectral modelling of SNII atmospheres indicates a clear dependence of metal line strengths on progenitor metallicity. This motivates further work to evaluate the accuracy with which these SNe can be used as metallicity indicators. To assess this accuracy we present a sample of SNII HII-region spectroscopy, from which environment abundances are derived. These environment abundances are compared to the observed strength of metal lines in SN spectra. Combining our sample with measurements from the literature, we present oxygen abundances of 119 host HII regions, by extracting emission line fluxes and using abundance diagnostics. Then, following Dessart et al., these abundances are compared to equivalent widths of Fe 5018 A at various time and colour epochs. Our distribution of inferred SNII host HII-region abundances has a range of ~0.6 dex. We confirm the dearth of SNeII exploding at metallicities lower than those found (on average) in the Large Magellanic Cloud. The equivalent width of Fe 5018 A at 50 days post explosion shows a statistically significant correlation with host HII-region oxygen abundance. The strength of this correlation increases if one excludes abundance measurements derived far from SN explosion sites. The correlation significance also increases if we only analyse a 'gold' IIP sample, and if a colour epoch is used in place of time. In addition, no evidence is found of correlation between progenitor metallicity and SN light-curve or spectral properties - except for that stated above with respect to Fe 5018 A equivalent width - suggesting progenitor metallicity is not a driving factor in producing the diversity observed in our sample. This study provides observational evidence of the usefulness of SNII as metallicity indicators. We finish with a discussion of the methodology needed to use SN spectra as independent metallicity diagnostics throughout the Universe.

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\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ 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$\mu$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.

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 $\tau$ 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 $\tau<$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 $\tau>$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.

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.

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.

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".

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.

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.

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.

We present a spatially resolved, high-spectral resolution (R=12000) K-band temporal monitoring of Rigel using AMBER at the VLTI. Rigel was observed in the Bracket Gamma line and its nearby continuum in 2006-2007, and 2009-2010. These unprecedented observations were complemented by contemporaneous optical high-resolution spectroscopy. We analyse the near-IR spectra and visibilities with the 1D non-LTE radiative-transfer code CMFGEN. The differential and closure phase signal exhibit asymmetries that are interpreted as perturbations of the wind. A systematic visibility decrease is observed across the Bracket Gamma. During the 2006-2007 period the Bracket Gamma and likely the continuum forming regions were larger than in the 2009-2010 epoch. Using CMFGEN, we infer a mass-loss rate change of about 20% between the two epochs. We further find time variations in the differential visibilities and phases. The 2006-2007 period is characterized by noticeable variations of the differential visibilities in Doppler position and width and by weak variations in differential and closure phase. The 2009-2010 period is much more quiet with virtually no detectable variations in the dispersed visibilities but a strong S-shape signal is observed in differential phase coinciding with a strong ejection event discernible in the optical spectra. The differential phase signal that is sometimes detected is reminiscent of the signal computed from hydrodynamical models of corotating interaction regions. For some epochs the temporal evolution of the signal suggests the rotation of the circumstellar structures.

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].

In classical P-Cygni profiles, theory predicts emission to peak at zero rest velocity. However, supernova spectra exhibit emission that is generally blue shifted. While this characteristic has been reported in many supernovae, it is rarely discussed in any detail. Here we present an analysis of H-alpha emission-peaks using a dataset of 95 type II supernovae, quantifying their strength and time evolution. Using a post-explosion time of 30d, we observe a systematic blueshift of H-alpha emission, with a mean value of -2000 kms-1. This offset is greatest at early times but vanishes as supernovae become nebular. Simulations of Dessart et al. (2013) match the observed behaviour, reproducing both its strength and evolution in time. Such blueshifts are a fundamental feature of supernova spectra as they are intimately tied to the density distribution of ejecta, which falls more rapidly than in stellar winds. This steeper density structure causes line emission/absorption to be much more confined; it also exacerbates the occultation of the receding part of the ejecta, biasing line emission to the blue for a distant observer. We conclude that blue-shifted emission-peak offsets of several thousand kms-1 are a generic property of observations, confirmed by models, of photospheric-phase type II 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.

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.

We explore a method for metallicity determinations based on quantitative spectroscopy of type II-Plateau (II-P) supernovae (SNe). For consistency, we first evolve a set of 15Msun main sequence stars at 0.1, 0.4, 1, and 2 x the solar metallicity. At the onset of core collapse, we trigger a piston-driven explosion and model the resulting ejecta and radiation. Our theoretical models of such red-supergiant-star explosions at different metallicity show that synthetic spectra of SNe II-P possess optical signatures during the recombination phase that are sensitive to metallicity variations. This sensitivity can be quantified and the metallicity inferred from the strength of metal-line absorptions. Furthermore, these signatures are not limited to O, but also include Na, Ca, Sc, Ti, or Fe. When compared to a sample of SNe II-P from the Carnegie SN Project and previous SN followup programs, we find that most events lie at a metallicity between 0.4 and 2 x solar, with a marked scarcity of SN II-P events at SMC metallicity. This most likely reflects the paucity of low metallicity star forming regions in the local Universe. SNe II-P have high plateau luminosities that make them observable spectroscopically at large distances. Because they exhibit signatures of diverse metal species, in the future they may offer a means to constrain the evolution of the composition (e.g., the O/Fe ratio) in the Universe out to a redshift of one and beyond.

We initiated long-term optical interferometry monitoring of the diameters of unstable yellow hypergiants (YHG) with the goal of detecting both the long-term evolution of their radius and shorter term formation related to large mass-loss events. We observed HR5171 A with AMBER/VLTI. We also examined archival photometric data in the visual and near-IR spanning more than 60 years, as well as sparse spectroscopic data. HR5171A exhibits a complex appearance. Our AMBER data reveal a surprisingly large star for a YHG R*=1315+/-260Rsun\ (~6.1AU) at the distance of 3.6+/-0.5kpc. The source is surrounded by an extended nebulosity, and these data also show a large level of asymmetry in the brightness distribution of the system, which we attribute to a newly discovered companion star located in front of the primary star. The companion's signature is also detected in the visual photometry, which indicates an orbital period of Porb=1304+/-6d. Modeling the light curve with the NIGHTFALL program provides clear evidence that the system is a contact or possibly over-contact eclipsing binary. A total current system mass of 39^{+40}_{-22} solar mass and a high mass ratio q>10 is inferred for the system. The low-mass companion of HR5171 A is very close to the primary star that is embedded within its dense wind. Tight constraints on the inclination and vsini of the primary are lacking, which prevents us from determining its influence precisely on the mass-loss phenomenon, but the system is probably experiencing a wind Roche-Lobe overflow. Depending on the amount of angular momentum that can be transferred to the stellar envelope, HR5171 A may become a fast-rotating B[e]/Luminous Blue Variable (LBV)/Wolf-Rayet star. In any case, HR5171 A highlights the possible importance of binaries for interpreting the unstable YHGs and for massive star evolution in general.

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.

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.

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.

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.

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.

We incorporate non-thermal excitation and ionization processes arising from non-thermal electrons that result from \gamma-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 \mu 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\alpha. 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.

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.

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

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.

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.

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.

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$\alpha$ 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.

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.

BA-type supergiants are amongst the most optically-bright stars. They are observable in extragalactic environments, hence potential accurate distance indicators. Emission activity in the Halpha line of the BA supergiants Rigel (B8Ia) and Deneb (A2Ia) is indicative of presence of localized time-dependent mass ejections. Here, we employ optical interferometry to study the Halpha line-formation region in these stellar environments. High spatial- (0.001 arcsec) and spectral- (R=30 000) resolution observations of Halpha were obtained with the visible recombiner VEGA installed on the CHARA interferometer, using the S1S2 array-baseline (34m). Six independent observations were done on Deneb over the years 2008 and 2009, and two on Rigel in 2009. We analyze this dataset with the 1D non-LTE radiative-transfer code CMFGEN, and assess the impact of the wind on the visible and near-IR interferometric signatures, using both Balmer-line and continuum photons. We observe a visibility decrease in Halpha for both Rigel and Deneb, suggesting that the line-formation region is extended (1.5-1.75 R*). We observe a significant visibility decrease for Deneb in the SiII6371 line. We witness time variations in the differential phase for Deneb, implying an inhomogeneous and unsteady circumstellar environment, while no such variability is seen in differential visibilities. Radiative-transfer modeling of Deneb, with allowance for stellar-wind mass loss, accounts fairly well for the observed decrease in the Halpha visibility. Based on the observed differential visibilities, we estimate that the mass-loss rate of Deneb has changed by less than 5%.

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.

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.

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.

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

- Constraining the cosmological parameters and understanding Dark Energy have tremendous implications for the nature of the Universe and its physical laws. - The pervasive limit of systematic uncertainties reached by cosmography based on Cepheids and Type Ia supernovae (SNe Ia) warrants a search for complementary approaches. - Type II SNe have been shown to offer such a path. Their distances can be well constrained by luminosity-based or geometric methods. Competing, complementary, and concerted efforts are underway, to explore and exploit those objects that are extremely well matched to next generation facilities. Spectroscopic follow-up will be enabled by space- based and 20-40 meter class telescopes. - Some systematic uncertainties of Type II SNe, such as reddening by dust and metallicity effects, are bound to be different from those of SNe Ia. Their stellar progenitors are known, promising better leverage on cosmic evolution. In addition, their rate - which closely tracks the ongoing star formation rate - is expected to rise significantly with look- back time, ensuring an adequate supply of distant examples. - These data will competitively constrain the dark energy equation of state, allow the determination of the Hubble constant to 5%, and promote our understanding of the processes involved in the last dramatic phases of massive stellar evolution.

We present a multi-epoch quantitative spectroscopic analysis of the Type IIn SN 1994W, an event interpreted by Chugai et al. as stemming from the interaction between the ejecta of a SN and a 0.4Msun circumstellar shell ejected 1.5yr before core collapse. During the brightening phase, our models suggest that the source of optical radiation is not unique, perhaps associated with an inner optically-thick Cold Dense Shell (CDS) and outer optically-thin shocked material. During the fading phase, our models support a single source of radiation, an hydrogen-rich optically-thick layer with a near-constant temperature of ~7000K that recedes from a radius of 4.3x10^15 at peak to 2.3x10^15cm 40 days later. We reproduce the hybrid narrow-core broad-wing line profile shapes of SN 1994W at all times, invoking an optically-thick photosphere exclusively (i.e., without any external optically-thick shell). In SN 1994W, slow expansion makes scattering with thermal electrons a key escape mechanism for photons trapped in optically-thick line cores, and allows the resulting broad incoherent electron-scattering wings to be seen around narrow line cores. In SNe with larger expansion velocities, the thermal broadening due to incoherent scattering is masked by the broad profile and the dominant frequency redshift occasioned by bulk motions. Given the absence of broad lines at all times and the very low 56Ni yields, we speculate whether SN 1994W could have resulted from an interaction between two ejected shells without core collapse. The high conversion efficiency of kinetic to thermal energy may not require a SN-like energy budget for SN1994W.

We report the serendipitous detection by GALEX of fast (<1 day) rising (>1 mag) UV emission from two Type II plateau (II-P) supernovae (SNe) at z=0.185 and 0.324 discovered by the Supernova Legacy Survey. Optical photometry and VLT spectroscopy 2 weeks after the GALEX detections link the onset of UV emission to the time of shock breakout. Using radiation hydrodynamics and non-LTE radiative transfer simulations, and starting from a standard red supergiant (RSG; Type II-P SN progenitor) star evolved self-consistently from the main sequence to iron core collapse, we model the shock breakout phase and the 55 hr that follow. The small scale height of our RSG atmosphere model suggests that the breakout signature is a thermal soft X-ray burst (lambda_peak ~ 90\AA) with a duration of <~ 2000 s. Longer durations are possible but require either an extended and tenuous non-standard envelope, or an unusually dense RSG wind with \dot{M} ~ 10^(-3) Msun yr^(-1). The GALEX observations miss the peak of the luminous (M_FUV ~ -20) UV burst but unambiguously capture the rise of the emission and a subsequent 2 day long plateau. The postbreakout, UV-bright plateau is a prediction of our model in which the shift of the peak of the spectral energy distribution (SED) from ~100 to ~1000\AA and the ejecta expansion both counteract the decrease in bolometric luminosity from ~10^11 to ~10^9 L_sun over that period. Based on the observed detection efficiency of our study we make predictions for the breakout detection rate of the GALEX Time Domain Survey.

Massive stars deeply influence their surroundings by their luminosity and the injection of kinetic energy. So far, they have mostly been studied with spatially unresolved observations, although evidence of geometrical complexity of their wind are numerous. Interferometry can provide spatially resolved observations of massive stars and their immediate vicinity. Specific geometries (disks, jets, latitude-dependent winds) can be probed by this technique. The first observation of a Wolf-Rayet (WR) star (\gamma^2 Vel) with the AMBER/VLTI instrument yielded to a re-evaluation of its distance and an improved characterization of the stellar components, from a very limited data-set. This motivated our team to increase the number of WR targets observed with AMBER. We present here new preliminary results that encompass several spectral types, ranging from early WN to evolved dusty WC. We present unpublished data on WR79a, a massive star probably at the boundary between the O and Wolf- Rayet type, evidencing some Wolf-Rayet broad emission lines from an optically thin wind. We also present new data obtained on \gamma^2 Vel that can be compared to the up-to-date interferometry-based orbital parameters from North et al. (2007). We discuss the presence of a wind-wind collision zone in the system and provide preliminary analysis suggesting the presence of such a structure in the data. Then, we present data obtained on 2 dusty Wolf-Rayet stars: WR48a-b and WR118, the latter exhibiting some clues of a pinwheel-like structure from the visibility variations.

We present VULCAN/2D multi-group flux-limited-diffusion radiation hydrodynamics simulations of binary neutron star (BNS) mergers, using the Shen equation of state, covering ~100 ms, and starting from azimuthal-averaged 2D slices obtained from 3D SPH simulations of Rosswog & Price for 1.4 Msun (baryonic) neutron stars with no initial spins, co-rotating spins, and counter-rotating spins. Snapshots are post-processed at 10 ms intervals with a multi-angle neutrino-transport solver. We find polar-enhanced neutrino luminosities, dominated by $\bar{\nu}_e$ and ``$\nu_\mu$'' neutrinos at peak, although $\nu_e$ emission may be stronger at late times. We obtain typical peak neutrino energies for $\nu_e$, $\bar{\nu}_e$, and ``$\nu_\mu$'' of ~12, ~16, and ~22 MeV. The super-massive neutron star (SMNS) formed from the merger has a cooling timescale of ~1 s. Charge-current neutrino reactions lead to the formation of a thermally-driven bipolar wind with <$\dot{M}$> ~10$^{-3}$ Msun/s, baryon-loading the polar regions, and preventing any production of a GRB prior to black-hole formation. The large budget of rotational free energy suggests magneto-rotational effects could produce a much greater polar mass loss. We estimate that ~10$^{-4}$ Msun of material with electron fraction in the range 0.1-0.2 become unbound during this SMNS phase as a result of neutrino heating. We present a new formalism to compute the $\nu_i\bar{\nu}_i$ annihilation rate based on moments of the neutrino specific intensity computed with our multi-angle solver. Cumulative annihilation rates, which decay as $t^{-1.8}$, decrease over our 100 ms window from a few 10$^{50}$ to ~10$^{49}$ erg/s, equivalent to a few 10$^{54}$ to ~10$^{53}$ $e^-e^+$ pairs per second.

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.

The peculiar Type Ib supernova (SN) 2006jc has been observed with the UV/Optical Telescope (UVOT) and X-Ray Telescope (XRT) on board the Swift observatory over a period of 19 to 183 days after the explosion. Signatures of interaction of the outgoing SN shock with dense circumstellar material (CSM) are detected, such as strong X-ray emission (L_{0.2-10} > E39 erg/s) and the presence of MgII 2800A line emission visible in the UV spectra. In combination with a Chandra observation obtained on day 40 after the explosion, the X-ray light curve is constructed, which shows a unique rise of the X-ray emission by a factor of ~5 over a period of ~4 months, followed by a rapid decline. We interpret the unique X-ray and UV properties as a result of the SN shock interacting with a shell of material that was deposited by an outburst of the SN progenitor two years prior to the explosion. Our results are consistent with the explosion of a Wolf-Rayet star that underwent an episodic mass ejection qualitatively similar to those of luminous blue variable stars prior to its explosion. This led to the formation of a dense (>E7 cm**-3) shell at a distance of ~E16 cm from the site of the explosion, which expands with the WR wind at a velocity of (1300+/-300) km/s.

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars yield iron cores at collapse endowed with high angular momentum. It is thought that high angular momentum and black hole formation are critical ingredients of the collapsar model of long-soft gamma-ray bursts (GRBs). Here, we present 2D multi-group, flux-limited-diffusion MHD simulations of the collapse, bounce, and immediate post-bounce phases of a 35-Msun collapsar-candidate model of Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the differentially-rotating surface layers of the millisecond-period neutron star, a magnetically-driven explosion ensues during the proto-neutron star phase, in the form of a baryon-loaded non-relativistic jet, and that a black hole, central to the collapsar model, does not form. Paradoxically, and although much uncertainty surrounds stellar mass loss, angular momentum transport, magnetic fields, and the MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, explosion in the immediate post-bounce phase. We surmise that fast rotation in the iron core may inhibit, rather than enable, collapsar formation, which requires a large angular momentum not in the core but above it. Variations in the angular momentum distribution of massive stars at core collapse might explain both the diversity of Type Ic supernovae/hypernovae and their possible association with a GRB. A corollary might be that, rather than the progenitor mass, the angular momentum distribution, through its effect on magnetic field amplification, distinguishes these outcomes.

We analyze the Type II Plateau supernovae (SN II-P) 2005cs and 2006bp with the non-LTE model atmosphere code CMFGEN. We fit 13 spectra in the first month for SN 2005cs and 18 for SN 2006bp. {\sl Swift} ultraviolet photometry and ground-based optical photometry calibrate each spectrum. Our analysis shows both objects were discovered less than 3 days after they exploded, making these the earliest SN II-P spectra ever studied. They reveal broad and very weak lines from highly-ionized fast ejecta with an extremely steep density profile. We identify He{\sc ii} 4686\AA emission in the SN 2006bp ejecta. Days later, the spectra resemble the prototypical Type II-P SN 1999em, which had a supergiant-like photospheric composition. Despite the association of SN 2005cs with possible X-ray emission, the emergent UV and optical light comes from the photosphere, not from circumstellar emission. We surmise that the very steep density fall-off we infer at early times may be a fossil of the combined actions of the shock wave passage and radiation driving at shock breakout. Based on tailored CMFGEN models, the direct-fitting technique and the Expanding Photosphere Method both yield distances and explosion times that agree within a few percent. We derive a distance to NGC 5194, the host of SN 2005cs, of 8.9$\pm$0.5 Mpc and 17.5$\pm$0.8 Mpc for SN 2006bp in NGC 3953. The luminosity of SN 2006bp is 1.5 times that of SN 1999em, and 6 times that of SN 2005cs. Reliable distances to Type II-P supernovae that do not depend on a small range in luminosity provide an independent route to the Hubble Constant and improved constraints on other cosmological parameters.

Spectroscopic modeling of Type II supernovae (SNe) generally assumes steady-state. Following the recent suggestion of Utrobin & Chugai, but using the 1D non-LTE line-blanketed model atmosphere code CMFGEN, we investigate the effects of including time-dependent terms that appear in the statistical and radiative equilibrium equations. We base our discussion on the ejecta properties and the spectroscopic signatures obtained from time-dependent simulations, investigating different ejecta configurations, and covering their evolution from one day to six weeks after shock breakout. Compared to equivalent steady-state models, our time-dependent models produce SN ejecta that are systematically over-ionized, affecting helium at one week after explosion, but ultimately affecting all ions after a few weeks. While the continuum remains essentially unchanged, time-dependence effects on observed spectral lines are large. At the recombination epoch, HI lines and NaID are considerably stronger and broader than in equivalent steady-state models, while CaII8500A is weakened. If time dependence is allowed for, the HeI lines at 5875A and 10830A appear about 3 times stronger at one week, and HeI10830A persists as a blue-shifted absorption feature even at 6 weeks after explosion. Time dependence operates through the energy gain from changes in ionization and excitation, and, perhaps more universally across SN types, from the competition between recombination and expansion, which in-turn, can be affected by optical-depth effects. Our time-dependent models compare well with observations of the low-luminosity low-velocity SN 1999br and the more standard SN 1999em, reproducing the Halpha line strength at the recombination epoch, and without the need for setting unphysical requirements on the magnitude of nickel mixing.

We present 2D multi-group flux-limited diffusion magnetohydrodynamics (MHD) simulations of the Accretion-Induced Collapse (AIC) of a rapidly-rotating white dwarf. We focus on the dynamical role of MHD processes after the formation of a millisecond-period protoneutron star. We find that including magnetic fields and stresses can lead to a powerful explosion with an energy of a few Bethe, rather than a weak one of at most 0.1 Bethe, with an associated ejecta mass of ~0.1Msun, instead of a few 0.001Msun. The core is spun down by ~30% within 500ms after bounce, and the rotational energy extracted from the core is channeled into magnetic energy that generates a strong magnetically-driven wind, rather than a weak neutrino-driven wind. Baryon loading of the ejecta, while this wind prevails, precludes it from becoming relativistic. This suggests that a GRB is not expected to emerge from such AICs during the early protoneutron star phase, except in the unlikely event that the massive white dwarf has sufficient mass to lead to black hole formation. In addition, we predict both negligible 56Ni-production (that should result in an optically-dark, adiabatically-cooled explosion) and the ejection of 0.1Msun of material with an electron fraction of 0.1-0.2. Such pollution by neutron-rich nuclei puts strong constraints on the possible rate of such AICs. Moreover, being free from ``fallback,'' such highly-magnetized millisecond-period protoneutron stars may later become magnetars, and the magnetically-driven winds may later transition to Poynting-flux-dominated, relativistic winds, eventually detectable as GRBs at cosmological distances. However, the low expected event rate of AICs will constrain them to be, at best, a small subset of GRB and/or magnetar progenitors.

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.

We present results on the X-ray and optical/UV emission from the type IIP SN 2006bp and the interaction of the SN shock with its environment, obtained with the X-Ray Telescope (XRT) and UV/Optical Telescope (UVOT) on-board the Swift observatory. SN 2006bp is detected in X-rays at a 4.5 sigma level of significance in the merged XRT data from days 1 to 12 after the explosion. If the X-ray luminosity of (1.8+/-0.4)E39 ergs/s is caused by interaction of the SN shock with circumstellar material (CSM), deposited by a stellar wind from the progenitor's companion star, a mass-loss rate of ~E-05 M_sun/yr is inferred. The mass-loss rate is consistent with the non-detection in the radio with the VLA on days 2, 9, and 11 after the explosion and characteristic of a red supergiant progenitor with a mass around 12-15 M_sun prior to the explosion. In combination with a follow-up XMM-Newton observation obtained on day 21 after the explosion, an X-ray rate of decline with index 1.2+/-0.6 is inferred. Since no other SN has been detected in X-rays prior to the optical peak and since type IIP SNe have an extended 'plateau' phase in the optical, we discuss the scenario that the X-rays might be due to inverse Compton scattering of photospheric optical photons off relativistic electrons produced in circumstellar shocks. However, due to the high required value of the Lorentz factor (~10-100) we conclude that Inverse Compton scattering is an unlikely explanation for the observed X-ray emission. The fast evolution of the optical/ultraviolet spectral energy distribution and the spectral changes observed with Swift reveal the onset of metal line-blanketing and cooling of the expanding photosphere during the first few weeks after the outburst.

We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of $div(\bB)$ by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant.

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.

In this paper, we bring together various of our published and unpublished findings from our recent 2D multi-group, flux-limited radiation hydrodynamic simulations of the collapse and explosion of the cores of massive stars. Aided by 2D and 3D graphical renditions, we motivate the acoustic mechanism of core-collapse supernova explosions and explain, as best we currently can, the phases and phenomena that attend this mechanism. Two major foci of our presentation are the outer shock instability and the inner core g-mode oscillations. The former sets the stage for the latter, which damp by the generation of sound. This sound propagates outward to energize the explosion and is relevant only if the core has not exploded earlier by some other means. Hence, it is a more delayed mechanism than the traditional neutrino mechanism that has been studied for the last twenty years since it was championed by Bethe and Wilson. We discuss protoneutron star convection, accretion-induced-collapse, gravitational wave emissions, pulsar kicks, the angular anisotropy of the neutrino emissions, a subset of numerical issues, and a new code we are designing that should supercede our current supernova code VULCAN/2D. Whatever ideas last from this current generation of numerical results, and whatever the eventual mechanism(s), we conclude that the breaking of spherical symmetry will survive as one of the crucial keys to the supernova puzzle.

In this work, we present the first AMBER observations, of the Wolf-Rayet and O (WR+O) star binary system y2 Velorum. The AMBER instrument was used with the telescopes UT2, UT3, and UT4 on baselines ranging from 46m to 85m. It delivered spectrally dispersed visibilities, as well as differential and closure phases, with a resolution R = 1500 in the spectral band 1.95-2.17 micron. We interpret these data in the context of a binary system with unresolved components, neglecting in a first approximation the wind-wind collision zone flux contribution. We show that the AMBER observables result primarily from the contribution of the individual components of the WR+O binary system. We discuss several interpretations of the residuals, and speculate on the detection of an additional continuum component, originating from the free-free emission associated with the wind-wind collision zone (WWCZ), and contributing at most to the observed K-band flux at the 5% level. The expected absolute separation and position angle at the time of observations were 5.1&plusmn;0.9mas and 66&plusmn;15&deg; respectively. However, we infer a separation of 3.62+0.11-0.30 mas and a position angle of 73+9-11&deg;. Our analysis thus implies that the binary system lies at a distance of 368+38-13 pc, in agreement with recent spectrophotometric estimates, but significantly larger than the Hipparcos value of 258+41-31 pc.

In the context of 2D, axisymmetric, multi-group, radiation/hydrodynamic simulations of core-collapse supernovae over the full 180$^{\circ}$ domain, we present an exploration of the progenitor dependence of the acoustic mechanism of explosion. All progenitor models we have tested with our Newtonian code explode. We investigate the roles of the Standing-Accretion-Shock-Instability (SASI), the excitation of core g-modes, the generation of core acoustic power, the ejection of matter with r-process potential, the wind-like character of the explosion, and the fundamental anisotropy of the blasts. We find that the breaking of spherical symmetry is central to the supernova phenomenon and the blasts, when top-bottom asymmetric, are self-collimating. We see indications that the initial explosion energies are larger for the more massive progenitors, and smaller for the less massive progenitors, and that the neutrino contribution to the explosion energy may be an increasing function of progenitor mass. The degree of explosion asymmetry we obtain is completely consistent with that inferred from the polarization measurements of Type Ic supernovae. Furthermore, we calculate for the first time the magnitude and sign of the net impulse on the core due to anisotropic neutrino emission and suggest that hydrodynamic and neutrino recoils in the context of our asymmetric explosions afford a natural mechanism for observed pulsar proper motions. [abridged]

We have incorporated time-dependent terms into the statistical and radiative equilibrium calculations of the non-LTE line-blanketed radiative transfer code CMFGEN. To illustrate the significant improvements in spectral fitting achieved for photospheric phase Type II SN, and to document the effects associated with time dependence, we model the outer 6.1Msun of ejecta of a BSG/RSG progenitor star. Hopping by 3-day increments, we compute the UV to near-IR spectral evolution for both continuum and lines, from the fully ionized conditions at one week to the partially recombined conditions at 6 weeks after the explosion. We confirm the importance of allowing for time-dependence in the modeling of Type-II SN, as recently discussed by Utrobin & Chugai for SN1987A. However, we allow in our approach for the full interaction between the radiation field and level populations, and study the effects on the full spectrum. At the recombination epoch, HI lines and NaD are stronger and broader than in equivalent steady-state models, while CaII is weakened. Former successes of steady-state CMFGEN models are unaffected, while former discrepancies are cured. Time dependence affects all lines, while the continuum, from the UV to the optical, changes only moderately. We identify two key effects: First, time dependence together with the energy gain through changes in ionization and excitation lead to an over-ionization in the vicinity of the photosphere, dramatically affecting line optical depths and profiles. Second, the ionization is frozen-in at large radii/velocities. This stems solely from the time-scale contrast between recombination and expansion and will occur, modulo non-thermal excitation effects, in all SN types. [abridged]

We present a new theory for the gravitational wave signatures of core-collapse supernovae. Previous studies identified axisymmetric rotating core collapse, core bounce, postbounce convection, and anisotropic neutrino emission as the primary processes and phases for the radiation of gravitational waves. Our results, which are based on axisymmetric, Newtonian radiation-hydrodynamics supernova simulations (Burrows et al. 2006), indicate that the dominant emission process of gravitational waves in core-collapse supernovae may be the oscillations of the protoneutron star core. The oscillations are predominantly of g-mode character, are excited hundreds of milliseconds after bounce, and typically last for several hundred milliseconds. Our results suggest that even nonrotating core-collapse supernovae should be visible to current LIGO-class detectors throughout the Galaxy, and depending on progenitor structure, possibly out to Megaparsec distances.

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-Msun and a 1.92-Msun model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50-100ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a neutrino-driven wind that develops within the white dwarf, in a cone of ~40deg opening angle about the poles, with a mass loss rate of 5-8 x 10^{-3} Msun/yr. The ejecta have an entropy on the order of 20-50 k_B/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, at >600ms after bounce, the explosion energy has reached 3-4 x 10^{49}erg and the total ejecta mass has reached a few times 0.001Msun. We estimate the asymptotic explosion energies to be lower than 10^{50}erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion. We document in detail the numerous effects of the fast rotation of the progenitors: The neutron stars are aspherical; the ``nu_mu'' and anti-nu_e neutrino luminosities are reduced compared to the nu_e neutrino luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction depend strongly upon latitude (a la von Zeipel); and a quasi-Keplerian 0.1-0.5-Msun accretion disk is formed.

Based on multi-dimensional, multi-group, flux-limited-diffusion hydrodynamic simulations of core-collapse supernovae with the VULCAN/2D code, we study the physical conditions within and in the vicinity of the nascent protoneutron star (PNS). Our numerical study follows the evolution of the collapsing envelope of the 11-Msun model of Woosley & Weaver, from ~200 ms before bounce to ~300ms after bounce, on a spatial grid that switches from Cartesian at the PNS center to spherical above a 10km radius. As has been shown previously, we do not see any large-scale overturn of the inner PNS material. Convection, directly connected to the PNS, is found to occur in two distinct regions: between 10 and 20km, coincident with the region of negative lepton gradient, and exterior to the PNS above 50km. Separating these two regions, an interface, with no sizable outward or inward motion, is the site of gravity waves, emerging at ~200-300ms after core bounce, excited by the convection in the outer convective zone. In the PNS, convection is always confined within the neutrinospheric radii for all neutrino energies above just a few MeV. We find that such convective motions do not appreciably enhance the nu_e neutrino luminosity, and that they enhance the anti-nu_e and "nu_mu" luminosities modestly, by ~15% and ~30%, respectively, during the first post-bounce 100--200ms. Moreover, we see no evidence of doubly-diffusive instabilities in the PNS, expected to operate on diffusion timescales of at least a second, much longer than the millisecond timescale associated with PNS convection. PNS convection is thus found to be a secondary feature of the core-collapse phenomenon, rather than a decisive ingredient for a successful explosion.

We present a new mechanism for core-collapse supernova explosions that relies upon acoustic power generated in the inner core as the driver. In our simulation using an 11-solar-mass progenitor, a strong advective-acoustic oscillation a la Foglizzo with a period of ~25-30 milliseconds (ms) arises ~200 ms after bounce. Its growth saturates due to the generation of secondary shocks, and kinks in the resulting shock structure funnel and regulate subsequent accretion onto the inner core. However, this instability is not the primary agent of explosion. Rather, it is the acoustic power generated in the inner turbulent region and most importantly by the excitation and sonic damping of core g-mode oscillations. An l=1 mode with a period of ~3 ms grows to be prominent around ~500 ms after bounce. The accreting protoneutron star is a self-excited oscillator. The associated acoustic power seen in our 11-solar-mass simulation is sufficient to drive the explosion. The angular distribution of the emitted sound is fundamentally aspherical. The sound pulses radiated from the core steepen into shock waves that merge as they propagate into the outer mantle and deposit their energy and momentum with high efficiency. The core oscillation acts like a transducer to convert accretion energy into sound. An advantage of the acoustic mechanism is that acoustic power does not abate until accretion subsides, so that it is available as long as it may be needed to explode the star. [abridged]

This work presents a detailed quantitative spectroscopic analysis of, and the determination of the distance to, the type II supernovae (SN) SN1999em with CMFGEN (Dessart & Hillier 2005a), based on spectrophotometric observations at eight dates up to 40 days after discovery. We use the same iron-group metal content for the ejecta, the same power-law density distribution (with exponent n~10), and a Hubble-velocity law at all times. We adopt a H/He/C/N/O abundance pattern compatible with CNO-cycle equilibrium values for a RSG/BSG progenitor, with C/O enhanced and N depleted at later times. Based on our synthetic fits to spectrophotometric observations of SN1999em, we obtain a distance of 11.5Mpc, similar to that of Baron et al. (2004) and the Cepheid distance to the galaxy host of 11.7Mpc (Leonard et al. 2003). Similarly, based on such models, the Expanding Photosphere Method (EPM) delivers a distance of 11.6Mpc, with negligible scatter between photometric bandpass sets; there is thus nothing wrong with the EPM as such. Previous determinations using the tabulated correction factors of Eastman et al. (1996) all led to 30-50% underestimates: we find that this is caused by 1) an underestimate of the correction factors compared to the only other study of the kind by Dessart & Hillier (2005b), 2) a neglect of the intrinsic >20% scatter of correction factors, and 3) the use of the EPM at late times when severe line blanketing makes the method inaccurate. The need of detailed model computations for reliable EPM distance estimates thus defeats the appeal and simplicity of the method. However, detailed fits to SN optical spectra, based on tailored models for individual SN observations, offers a promising approach to obtaining distances with 10-20% accuracy, either through the EPM or a la Baron et al. (2004).

Using archival data of low-redshift (z < 0.01) Type Ia supernovae (SN Ia) and recent observations of high-redshift (0.16 < z <0.64; Matheson et al. 2005) SN Ia, we study the "uniformity'' of the spectroscopic properties of nearby and distant SN Ia. We find no difference in the measures we describe here. In this paper, we base our analysis solely on line-profile morphology, focusing on measurements of the velocity location of maximum absorption (vabs) and peak emission (vpeak). We find that the evolution of vabs and vpeak for our sample lines (Ca II 3945, Si II 6355, and S II 5454, 5640) is similar for both the low- and high-redshift samples. We find that vabs for the weak S II 5454, 5640 lines, and vpeak for S II 5454, can be used to identify fast-declining [dm15 > 1.7] SN Ia, which are also subluminous. In addition, we give the first direct evidence in two high-z SN Ia spectra of a double-absorption feature in Ca II 3945, an event also observed, though infrequently, in low-redshift SN Ia spectra (6/22 SN Ia in our local sample). We report for the first time the unambiguous and systematic intrinsic blueshift of peak emission of optical P-Cygni line profiles in Type Ia spectra, by as much as 8000 km/s. All the high-z SN Ia analyzed in this paper were discovered and followed up by the ESSENCE collaboration, and are now publicly available.

Due to their high intrinsic brightness, type II supernovae (SN) can be used as lighthouses to constrain distances in the Universe using variants of the Baade-Wesselink method. Based on a large set of CMFGEN models (Hillier & Miller 1998) covering the photospheric phase of type II SN, we study the various concepts entering one such technique, the Expanding Photosphere Method (EPM). We compute correction factors $\xi$ needed to approximate the synthetic Spectral Energy Distribution (SED) with that of a blackbody at temperature $T$. Our $\xi$, although similar, are systematically greater, by $\sim 0.1$, than the values obtained by Eastman et al. (1996). This translate into a systematic enhancement of 10-20% in EPM-distances. We find that line emission and absorption, not directly linked to color temperature variations, can considerably alter the synthetic magnitude, and cause above-unity $\xi$-values. Optically-thick lines, present in the optical at late-times, can also introduce a strong wavelength-dependence of the photospheric radius, invalidating the use of the Baade method. Both the impact of line-blanketing on the SED and the photospheric radius at low $T$ suggest that the EPM is best used at early times, when the outflow is fully ionized and line-opacity mostly confined to the UV range. We also investigate how reliably one can measure the photospheric velocity from P-Cygni line profiles. Contrary to the usually held belief, the velocity at maximum absorption in the P-Cygni trough of optically-thick lines can both {\it overestimate or underestimate} the photospheric velocity, with a magnitude that depends on the SN outflow density gradient and the optical thickness of the line.

We present first results on the quantitative spectroscopic analysis of the photospheric-phase of Type II supernovae (SN). The analyses are based on the non-LTE model atmosphere code, CMFGEN, of Hillier & Miller (1998). A key asset of CMFGEN is its thorough treatment of line-blanketing due to metal species. From its applicability to hot star environments, the main modifications to the source code were to allow a linear velocity law, a power-law density distribution, an adaptive grid to handle the steep H recombination/ionization front occurring in cooler SN models, and a routine to compute the gray temperature structure in the presence of large velocities. In this first paper we demonstrate the ability of CMFGEN to reproduce, with a high level of accuracy, the UV and optical observations of a sample of well observed Type II SN, i.e. SN1987A and SN1999em, at representative stages of their photospheric evolution. Two principal stages of SN are modeled - that where hydrogen is fully ionized, and that in which H is only partially ionized. We provide detailed and accurate line identifications for our sample of epochs. We discuss the impact on the SED of varying the density exponent, the spatial scale of the SN, the value of Rmax and the metallicity; we illustrate the formation sites for selected lines. Both Hydrogen Balmer and Helium lines (when present) are very well reproduced, even for our adopted CNO-equilibrium abundance values adequate for a massive supergiant progenitor star. NII lines are also predicted in the blue-wing of both Hbeta and HeI 5875A, and thus not associated with peculiar absorption/emission far above the photosphere. We finally present a pedagogical discussion of P-Cygni profile formation and explain, e.g., the noticeable blue-shift of Halpha in early time spectra (much abridged abstract).

We present initial attempts to include the multi-dimensional nature of radiation transport in hydrodynamical simulations of the small-scale structure that arises from the line-driven instability in hot-star winds. Compared to previous 1D or 2D models that assume a purely radial radiation force, we seek additionally to treat the lateral momentum and transport of diffuse line-radiation, initially here within a 2D context. A key incentive is to study the damping effect of the associated diffuse line-drag on the dynamical properties of the flow, focusing particularly on whether this might prevent lateral break-up of shell structures at scales near the lateral Sobolev angle of ca. $1^{\rm o}$. We first explore nonlinear simulations that cast the lateral diffuse force in the simple, local form of a parallel viscosity. Second, to account for the lateral mixing of radiation associated with the radial driving, we next explore models in which the radial force is azimuthally smoothed over a chosen scale. Third, to account for both the lateral line-drag and the lateral mixing in a more self-consistent way, we explore further a method first proposed by Owocki (1999), which uses a restricted 3-ray approach that combines a radial ray with two oblique rays set to have an impact parameter $p < R_{\ast}$ within the stellar core. From numerical simulations, we find that, compared to equivalent 1-ray simulations, the high-resolution 3-ray models show systematically a much higher lateral coherence.... (Full abstract in paper)

The migration of profile sub-peaks identified in time-monitored optical emission lines of Wolf-Rayet star spectra provides a direct diagnostic of the dynamics of their stellar winds via a measured line-of-sight velocity change per unit time. Inferring the associated wind acceleration scale from such an apparent acceleration then relies on the adopted intrinsic velocity of the wind material at the origin of this variable pattern. We develop radiative transfer tools to characterise the Line Emission Region (LER) of the program lines, including turbulence and line-optical depth effects. We find that monitored-lines can be fitted well with a pure optically thin formation mechanism, line-broadening resulting from the finite velocity extent of the LER rather than turbulence.Our new estimates of LER velocity centroids are systematically shifted outwards closer to terminal velocity compared to previous determinations, now suggesting WR-wind acceleration length scales of ca. 10-20Rsun, a factor of a few smaller than previously inferred. Based on radiation-hydrodynamics simulations of the line-driven-instability mechanism, we compute synthetic line profile variability for CIII5696A, for WR111 The results match well the measured observed migration of 20-30 m/s2, However, our model stellar radius of 19 Rsun, typical of an O-type supergiant, is a factor 2--10 larger than generally expected for WR core radii.Such small radii leave inferred acceleration scales to be more extended than expected from dynamical models of line driving, but the severity of the discrepancy is substantially reduced compared to previous analyses. We conclude with a discussion of how using lines formed deeper in the wind would provide a stronger constraint on the key wind dynamics as well as information on the onset of wind clumping.

We use ultraviolet space-based (FUSE, HST) and optical/IR ground-based (2.3m MSSSO, NTT) spectroscopy to determine the physical parameters of six WC4-type Wolf-Rayet stars in the Large Magellanic Cloud. Stellar parameters are revised significantly relative to Grafener et al. (1998) based on improved observations and more sophisticated model atmosphere codes, which account for line blanketing and clumping. We find that stellar luminosities are revised upwards by up to 0.4 dex, with surface abundances spanning a lower range of 0.1<C/He<0.35 (20-45% carbon by mass) and O/He<0.06 (<10% oxygen by mass). Relative to Galactic WC5-8 stars at known distance, and analysed in a similar manner, LMC WC4 stars possess systematically higher stellar luminosities, 0.2dex lower wind densities, yet a similar range of surface chemistries. We illustrate how the classification CIII 5696 line is extremely sensitive to wind density, such that this is the principal difference between the subtype distribution of LMC and Galactic early-type WC stars. Temperature differences do play a role, but carbon abundance does not affect WC spectral types. We illustrate the effect of varying temperature and mass-loss rate on the WC spectral type for HD32257 (WC4, LMC) and HD156385 (WC7, Galaxy) which possess similar abundances and luminosities. Using the latest evolutionary models, pre-supernova stellar masses in the range 11-19 Mo are anticipated for LMC WC4 stars, with 7-14 Mo for Galactic WC stars with known distances. These values are consistent with pre-cursors of bright Type-Ic supernovae such as SN1998bw (alias GRB 980425) for which a minimum total mass of C and O of 14Mo has been independently derived.

We present Far-Ultraviolet Spectroscopic Explorer (FUSE) spectroscopy of Sand 2, a LMC WO-type Wolf-Rayet star, revealing the OVI resonance P Cygni doublet at 1032-38A. These data are combined with HST/FOS ultraviolet and Mt Stromlo 2.3m optical spectroscopy, and analysed using a spherical, non-LTE, line-blanketed code. Our study reveals exceptional stellar parameters: T*=150,000K, v_inf=4100 km/s, log (L/Lo)=5.3, and Mdot=10^-5 Mo/yr if we adopt a volume filling factor of 10%. Elemental abundances of C/He=0.7+-0.2 and O/He=0.15(-0.05+0.10) by number qualitatively support previous recombination line studies. We confirm that Sand 2 is more chemically enriched in carbon than LMC WC stars, and is expected to undergo a supernova explosion within the next 50,000 yr.

In this paper we derive stellar parameters for the Wolf-Rayet star in the gamma Velorum binary system (WR11), from a detailed non-LTE model of its optical and infrared spectra. Compared to the study of Schaerer et al., the parameters of the WC8 star are revised to a hotter effective temperature (Teff~57kK),a higher luminosity (log(L/Lsun = 5.00), and a lower mass loss rate (log(Mdot / Msun/yr) = -5.0 - using a 10% clumping filling factor). These changes lead to a significant decrease in wind efficiency number, from 144 to 7, so that the driving mechanism of the wind of this WR star may be simply radiation pressure on lines. The derived spectroscopic luminosity is found to be 40% lower than that derived by De Marco & Schmutz through the mass-luminosity relationship for WR stars (log(L/Lsun = 5.2). The paper furthermore presents a comparison of the independently-developed modelling programs, CMFGEN and ISA-Wind . Overall, there seems to be very reasonable agreement between the derived parameters for WR11, except for the carbon content, which is 2 times higher for CMFGEN (C/He=0.15 vs. 0.06, by number).The comparison also confirms a disparity in the predicted flux at lambda<$400A, found by Crowther et al., which will have effects on several nebular line strengths. The paper also presents the first independent check of the photon loss mechanism proposed by Schmutz. We conclude that, not only is it important to include very many lines to realistically model line blanketing, but in particular those ones that critically interact with strong resonance lines (e.g. HeII lambda303.78). The inclusion of these latter lines may significantly alter the wind ionization structure.

Neon abundances are derived in four Galactic WC stars -- gamma Vel (WR11, WC8+O7.5III), HD156385 (WR90, WC7), HD192103 (WR135, WC8), and WR146 (WC5+O8) - using mid-infrared fine structure lines obtained with ISO/SWS. Stellar parameters for each star are derived using a non-LTE model atmospheric code (Hillier & Miller 1998) together with ultraviolet (IUE), optical (INT, AAT) and infrared (UKIRT, ISO) spectroscopy. In the case of gamma Vel, we adopt results from De Marco et al. (2000), who followed an identical approach. ISO/SWS datasets reveal the [NeIII] 15.5um line in each of our targets, while [NeII] 12.8um, [SIV] 10.5um and [SIII] 18.7um are observed solely in gamma Vel. Using a method updated from Barlow et al. (1988) to account for clumped winds, we derive Ne/He=3-4x10^-3 by number, plus S/He=6x10^-5 for gamma Vel. Neon is highly enriched, such that Ne/S in gamma Vel is eight times higher than cosmic values. However, observed Ne/He ratios are a factor of two times lower than predictions of current evolutionary models of massive stars. An imprecise mass-loss and distance were responsible for the much greater discrepancy in neon content identified by Barlow et al. Our sample of WC5--8 stars span a narrow range in T* (=55--71kK), with no trend towards higher temperature at earlier spectral type, supporting earlier results for a larger sample by Koesterke & Hamann (1995). Stellar luminosities range from 100,000 to 500,000 Lo, while 10^-5.1 < Mdot/(Mo/yr) < 10^-4.5, adopting clumped winds, in which volume filling factors are 10%. In all cases, wind performance numbers are less than 10, significantly lower than recent estimates. Carbon abundances span 0.08 < C/He < 0.25 by number, while oxygen abundances remain poorly constrained.