ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0$ mag (100d)$^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $^{56}$Ni mass of $M_{\rm Ni} \sim0.4\rm M_\odot$ and ejecta with high kinetic energy $E_{\rm kin}= (5-9)\times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3\, M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $^{56}$Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $^{56}$Ni yields. The earliest spectrum shows "flash ionization" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4} \, \rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 \rm~km~s^{-1}$ for H$\alpha$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.7-3.1 $\rm M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of 19-24 $\rm M_\odot$.

A crucial ingredient in population synthesis studies involving massive stars is the determination of whether they explode or implode in the end. While the final fate of a massive star is sensitive to its core structure at the onset of collapse, the existing binary population synthesis studies do not reach core-collapse. Instead, they employ simple prescriptions to infer their final fates without knowing the presupernova core structure. We explore a potential solution to this problem by treating the carbon-oxygen (CO) core independently from the rest of the star. Using the implicit hydrodynamics code KEPLER, we have computed an extensive grid of 3496 CO-core models from a diverse range of initial conditions, each evolved from carbon ignition until core-collapse. The final core structure, and thus the explodability, varies non-monotonically and depends sensitively on both the mass and initial composition of the CO-core. Although bare CO-cores are not perfect substitutes for cores embedded in massive stars, our models compare well both with MESA and full hydrogenic and helium star calculations. Our results can be used to infer the presupernova core structures from population synthesis estimates of CO-core properties, thus to determine the final outcomes based on the results of modern neutrino-driven explosion simulations. A sample application is presented for a population of Type-IIb supernova progenitors. All of our models are available at https://doi.org/10.5281/zenodo.3785377.

The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5$\%$ to 30$\%$, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32 - 1.37 $M_\odot$ regardless of metallicity. The median black hole mass for solar metallicity is typically 8 to 9 $M_\odot$ if only initial helium cores below 40 $M_\odot$ (ZAMS mass less than 80 $M_\odot$) are counted, and 9 - 13 $M_\odot$, in most cases, if helium cores with initial masses up to 150 $M_\odot$ (ZAMS mass less than 300 $M_\odot$) contribute. As long as the mass-loss rate as a function of mass exhibits no strong non-linearities, the black hole birth function from 15 to 35 $M_\odot$ has a slope that depends mostly on the initial mass function for main sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational wave radiation. The expected rotation rates of the black holes are briefly discussed.

We present photometric and spectroscopic observations of the unusual Type Ia supernova ASASSN-18tb, including a series of SALT spectra obtained over the course of nearly six months and the first observations of a supernova by the Transiting Exoplanet Survey Satellite (TESS). We confirm a previous observation by Kollmeier et al. (2019) showing that ASASSN-18tb is the first relatively normal Type Ia supernova to exhibit clear broad ($\sim1000$ km s$^{-1}$) H$\alpha$ emission in its nebular phase spectra. We find that this event is best explained as a sub-Chandrasekhar mass explosion with $M_{Ni} \approx 0.3\; \rm{M}_\odot$. Despite the strong H$\alpha$ signature at late times, we find that the early rise of the supernova shows no evidence for deviations from a single-component power-law and is best fit with a moderately shallow power-law of index $1.69\pm0.04$. We find that the H$\alpha$ luminosity remains approximately constant after its initial detection at phase +37 d, and that the H$\alpha$ velocity evolution does not trace that of the Fe~III$~\lambda4660$ emission. These suggest that the H$\alpha$ emission arises from circumstellar medium (CSM) rather than swept up material from a non-degenerate companion. However, ASASSN-18tb is strikingly different from other known CSM-interacting Type Ia supernovae in a number of significant ways. Those objects typically show an H$\alpha$ luminosity two orders of magnitude higher than what is seen in ASASSN-18tb, pushing them away from the empirical light-curve relations that define "normal" Type Ia supernovae. Conversely, ASASSN-18tb exhibits a fairly typical light curve and luminosity for an underluminous or transitional SN Ia, with $M_R \approx -18.1$ mag. Moreover, ASASSN-18tb is the only SN Ia showing H$\alpha$ from CSM interaction to be discovered in an early-type galaxy.

We present photometric and spectroscopic follow-up observations of the highly luminous Type Ibn supernova ASASSN-14ms, which was discovered on UT 2014-12-26.61 at $m_V \sim 16.5$. With a peak absolute $V$-band magnitude brighter than $-20.5$, a peak bolometric luminosity of $1.7 \times 10^{44}$ ergs s$^{-1}$, and a total radiated energy of $2.1 \times 10^{50}$ ergs, ASASSN-14ms is one of the most luminous Type Ibn supernovae yet discovered. In simple models, the most likely power source for this event is a combination of the radioactive decay of $^{56}$Ni and $^{56}$Co at late times and the interaction of supernova ejecta with the progenitor's circumstellar medium at early times, although we cannot rule out the possibility of a magnetar-powered light curve. The presence of a dense circumstellar medium is indicated by the intermediate-width He I features in the spectra. The faint ($m_g \sim 21.6$) host galaxy SDSS J130408.52+521846.4 has an oxygen abundance below $12+\log(O/H) \lesssim 8.3$, a stellar mass of $M_* \sim 2.6 \times 10^8 M_{\odot}$, and a star formation rate of $\textrm{SFR} \sim 0.02$ $M_{\odot}$ yr$^{-1}$.

Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored.

Recent studies on direct imaging of Type II core-collapse supernova progenitors indicate a possible threshold around $M_{\rm ZAMS}\sim 16-20$ M$_\odot$, where red supergiants with larger birth masses do not appear to result in supernova explosions and instead implode directly into a black hole. In this study we argue that it is not a coincidence that this threshold closely matches the critical transition of central Carbon burning in massive stars from the convective to radiative regime. In lighter stars, Carbon burns convectively in the center and result in compact final presupernova cores that are likely to result in explosions, while in heavier stars after the transition, it burns as a radiative flame and the stellar cores become significantly harder to explode. Using the KEPLER code we demonstrate the sensitivity of this transition to the rate of $^{12}$C$(\alpha,\gamma)^{16}$O reaction and the overshoot mixing efficiency, and we argue that the upper mass limit of exploding red supergiants could be employed to constrain uncertain input physics of massive stellar evolution calculations. The initial mass corresponding to the central Carbon burning transition range from 14 to 26 M$_\odot$ in recently published models from various groups and codes, and only a few are in agreement with the estimates inferred from direct imaging studies.

LOX is a lunar-orbiting astrophysics mission that will probe the cosmos at MeV energies. It is guided by open questions regarding thermonuclear, or Type-Ia, supernovae (SNeIa) and will characterize these inherently radioactive objects by enabling a systematic survey of SNeIa at gamma-ray energies for the first time. Astronomical investigations from lunar orbit afford new opportunities to advance our understanding of the cosmos. The foundation of LOX is an observational approach well suited to the all-sky monitoring demands of supernova investigations and time-domain astronomy. Its inherently wide field-of-view and continuous all-sky monitoring provides an innovative way of addressing decadal survey questions at MeV energies (0.1-10 MeV). The LOX approach achieves high sensitivity with a simple, high-heritage instrument design that eliminates the need for complex, position-sensitive detectors, kinematic event reconstruction, masks, or other insensitive detector mass, while also mitigating technology development, implementation complexity, and their associated costs. LOX can be realized within existing programs, like Explorer.

Supernova remnants (SNRs) that contain pulsar wind nebulae (PWNe) are characterized by distinct evolutionary stages. In very young systems, the PWN drives a shock into the innermost supernova (SN) material, giving rise to low-excitation lines and an infrared (IR) continuum from heated dust grains. These observational signatures make it possible to cleanly measure the properties of the deepest SN ejecta layers that can, in turn, provide constraints on the SN progenitor. We present Herschel Space Observatory far-IR observations of the PWN in the Galactic SNR Kes 75, containing the youngest known pulsar that exhibited magnetar-like activity. We detect highly-broadened oxygen and carbon line emission that arises from the SN ejecta encountered by the PWN. We also detect a small amount (a few thousandths of a solar mass) of shock-heated dust that spatially coincides with the ejecta material and was likely formed in the SN explosion. We use hydrodynamical models to simulate the evolution of Kes 75 and find that the PWN has so far swept up 0.05-0.1 solar masses of SN ejecta. Using explosion and nucleosynthesis models for different progenitor masses in combination with shock models, we compare the predicted far-IR emission with the observed line intensities and find that lower mass and explosion energy SN progenitors with mildly mixed ejecta profiles and comparable abundance fractions of carbon and oxygen are favored over higher mass ones. We conclude that Kes 75 likely resulted from an 8-12 solar-mass progenitor, providing further evidence that lower energy explosions of such progenitors can give rise to magnetars.

I show that the characteristic diffusion timescale and the gamma-ray escape timescale, of SN Ia supernova ejecta, are related with each other through the time when the bolometric luminosity, $L_{\rm bol}$, intersects with instantaneous radioactive decay luminosity, $L_\gamma$, for the second time after the light-curve peak. Analytical arguments, numerical radiation-transport calculations, and observational tests show that $L_{\rm bol}$ generally intersects $L_\gamma$ at roughly $1.7$ times the characteristic diffusion timescale of the ejecta. This relation implies that the gamma-ray escape timescale is typically 2.7 times the diffusion timescale, and also implies that the bolometric luminosity 15 days after the peak, $L_{\rm bol}(t_{15})$, must be close to the instantaneous decay luminosity at that time, $L_\gamma (t_{15})$. With the employed calculations and observations, the accuracy of $L_{\rm bol}=L_\gamma$ at $t=t_{15}$ is found to be comparable to the simple version of "Arnett's rule" ($L_{\rm bol}=L_\gamma$ at $t=t_{\rm peak}$). This relation aids the interpretation of SN Ia supernova light curves and may also be applicable to general hydrogen-free explosion scenarios powered by other central engines.

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.

The density structure surrounding the iron core of a massive star when it dies is known to have a major effect on whether or not the star explodes. Here we repeat previous surveys of presupernova evolution with some important corrections to code physics and four to ten times better mass resolution in each star. The number of presupernova masses considered is also much larger. Over 4,000 models are calculated in the range from 12 to 60 Msun with varying mass loss rates. The core structure is not greatly affected by the increased spatial resolution. The qualitative patterns of compactness measures and their extrema are the same, but with the increased number of models, the scatter seen in previous studies is replaced by several localized branches. More physics-based analyses by Ertl et al. (2016) and M$\"u$eller et al. (2016) show these branches with less scatter than the single parameter characterization of O'Connor & Ott (2011). These branches are particularly apparent for stars in the mass ranges 14 - 19 Msun and 22 - 24 Msun. The multi-valued solutions are a consequence of interference between several carbon and oxygen burning shells during the late stages of evolution. For a relevant range of masses, whether a star explodes or not may reflect more the small, almost random differences in its late evolution than its initial mass. The large number of models allow statistically meaningful statements about the radius, luminosity, and effective temperatures of presupernova stars, their core structures, and their remnant mass distribution.

We study the optical emission from heavy element ejecta in the oxygen-rich young supernova remnant (SNR) 1E 0102.2-7219 (1E 0102) in the Small Magellanic Cloud. We have used the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph at the Very Large Telescope (VLT) on Cerro Paranal and the wide field spectrograph (WiFeS) at the ANU 2.3 m telescope at Siding Spring Observatory to obtain deep observations of 1E 0102. Our observations cover the entire extent of the remnant from below 3500{\AA} to 9350{\AA}. Our observations unambiguously reveal the presence of fast-moving ejecta emitting in [S II], [S III], [Ar III], and [Cl II]. The sulfur-rich ejecta appear more asymmetrically distributed compared to oxygen or neon, a product of carbon-burning. In addition to the forbidden line emission from products of oxygen burning (S, Ar, Cl), we have also discovered H{\alpha} and H{\beta} emission from several knots of low surface brightness, fast-moving ejecta. The presence of fast-moving hydrogen points towards a progenitor that had not entirely shed its hydrogen envelope prior to the supernova. The explosion that gave rise to 1E 0102 is therefore commensurate with a Type IIb supernova.

The mass distribution of compact objects provides a fossil record that can be studied to uncover information on the late stages of massive star evolution, the supernova explosion mechanism, and the dense matter equation of state. Observations of neutron star masses indicate a bimodal Gaussian distribution, while the observed black hole mass distribution decays exponentially for stellar-mass black holes. We use these observed distributions to directly confront the predictions of stellar evolution models and the neutrino-driven supernova simulations of Sukhbold et al. (2016). We find excellent agreement between the black hole and low-mass neutron star distributions created by these simulations and the observations. We show that a large fraction of the stellar envelope must be ejected, either during the formation of stellar-mass black holes or prior to the implosion through tidal stripping due to a binary companion, in order to reproduce the observed black hole mass distribution. We also determine the origins of the bimodal peaks of the neutron star mass distribution, finding that the low-mass peak (centered at ~1.4 M_sun) originates from progenitors with M_zams ~ 9-18 M_sun. The simulations fail to reproduce the observed peak of high-mass neutron stars (centered at ~1.8 M_sun) and we explore several possible explanations. We argue that the close agreement between the observed and predicted black hole and low-mass neutron star mass distributions provides new promising evidence that these stellar evolution and explosion models are accurately capturing the relevant stellar, nuclear, and explosion physics involved in the formation of compact objects.

A large fraction of core-collapse supernovae (CCSNe), 30-50%, are expected to originate from the low-mass end of progenitors with $M_{\rm ZAMS}~= 8-12~M_\odot$. However, degeneracy effects make stellar evolution modelling of such stars challenging, and few predictions for their supernova light curves and spectra have been presented. Here we calculate synthetic nebular spectra of a 9 $M_\odot$ Fe CCSN model exploded with the neutrino mechanism. The model predicts emission lines with FWHM$\sim$1000 km/s, including signatures from each deep layer in the metal core. We compare this model to observations of the three subluminous IIP SNe with published nebular spectra; SN 1997D, SN 2005cs, and SN 2008bk. The prediction of both line profiles and luminosities are in good agreement with SN 1997D and SN 2008bk. The close fit of a model with no tuning parameters provides strong evidence for an association of these objects with low-mass Fe CCSNe. For SN 2005cs, the interpretation is less clear, as the observational coverage ended before key diagnostic lines from the core had emerged. We perform a parameterised study of the amount of explosively made stable nickel, and find that none of these three SNe show the high $^{58}$Ni/$^{56}$Ni ratio predicted by current models of electron capture SNe (ECSNe) and ECSN-like explosions. Combined with clear detection of lines from O and He shell material, these SNe rather originate from Fe core progenitors. We argue that the outcome of self-consistent explosion simulations of low-mass stars, which gives fits to many key observables, strongly suggests that the class of subluminous Type IIP SNe is the observational counterpart of the lowest mass CCSNe.

Nascent neutron stars with millisecond periods and magnetic fields in excess of $10^{16}$ Gauss can drive highly energetic and asymmetric explosions known as magnetar-powered supernovae. These exotic explosions are one theoretical interpretation for supernovae Ic-BL which are sometimes associated with long gamma-ray bursts. Twisted magnetic field lines extract the rotational energy of the neutron star and release it as a disk wind or a jet with energies greater than 10$^{52}$ erg over $\sim 20$ sec. What fractions of the energy of the central engine go into the wind and the jet remain unclear. We have performed two-dimensional hydrodynamical simulations of magnetar-powered supernovae (SNe) driven by disk winds and jets with the CASTRO code to investigate the effect of the central engine on nucleosynthetic yields, mixing, and light curves. We find that these explosions synthesize less than 0.05 Msun of Ni and that this mass is not very sensitive to central engine type. The morphology of the explosion can provide a powerful diagnostic of the properties of the central engine. In the absence of a circumstellar medium these events are not very luminous, with peak bolometric magnitudes $M_b \sim -16.5 $ due to low Ni production.

We have observed the oxygen-rich SNR 1E 0102.2-7219 with the integral field spectrograph WiFeS at Siding Spring Observatory and discovered sulfur-rich ejecta for the first time. Follow-up deep DDT observations with MUSE on the VLT (8100 s on source) have led to the additional discovery of fast- moving hydrogen as well as argon-rich and chlorine-rich material. The detection of fast-moving hydrogen knots challenges the interpretation that the progenitor of 1E 0102 was a compact core of a Wolf-Rayet star that had shed its entire envelope. In addition to the detection of hydrogen and the products of oxygen-burning, this unprecedented sharp (0.2" spaxel size at ~0.7" seeing) and deep MUSE view of an oxygen-rich SNR in the Magellanic Clouds reveals further exciting discoveries, including [Fe xiv]{\lambda}5303 and [Fe xi]{\lambda}7892 emission, which we associate with the forward shock. We present this exciting data set and discuss some of its implications for the explosion mechanism and nucleosynthesis of the associated supernova.

Multiple observational lines of evidence support a connection between hydrogen-poor superluminous supernovae (SLSNe) and long duration gamma-ray bursts (GRBs). Both events require a powerful central energy source, usually attributed to a millisecond magnetar or an accreting black hole. The GRB-SLSN link raises several theoretical questions: What distinguishes the engines responsible for these different phenomena? Can a single engine power both a GRB and a luminous SN in the same event? We propose a new unifying model for magnetar thermalization and jet formation: misalignment between the rotation (${\bf \Omega}$) and magnetic dipole (${\bf \mu}$) axes thermalizes a fraction of the spindown power by reconnection in the striped equatorial wind, providing a guaranteed source of "thermal" emission to power the supernova. The remaining un-thermalized power energizes a relativistic jet. In this picture, the GRB-SLSN dichotomy is directly linked to ${\bf \Omega \cdot \mu}$. We extend earlier work to show that even weak relativistic jets of luminosity $\sim10^{46}$ erg s$^{-1}$ can escape the expanding SN ejecta hours after the explosion, implying that escaping relativistic jets may accompany many SLSNe. We calculate the observational signature of these jets. We show that they may produce transient UV cocoon emission lasting a few hours when the jet breaks out of the ejecta surface. A longer-lived optical/UV signal may originate from a mildly-relativistic wind driven from the interface between the jet and the ejecta walls. This provides a new explanation for the secondary early-time maximum observed in some SLSNe light curves, such as LSQ14bdq. This scenario also predicts a population of GRB from on-axis jets with extremely long durations, potentially similar to the population of "jetted tidal disruption events", in coincidence with a small subset of SLSNe.

We investigate the properties of Type IIP supernovae that are dominantly powered by the rotational kinetic energy of the newly born neutron star. While the spin-down of a magnetar has previously been proposed as a viable energy source in the context of super-luminous supernovae, we show that a similar mechanism could produce both normal and peculiar Type IIP supernova light curves from red supergiant progenitors for a range of initial spin periods and equivalent dipole magnetic field strengths. Although the formation channel for such magnetars in a typical red supergiant progenitor is unknown, it is tantalizing that this proof of concept model is capable of producing ordinary Type IIP lightcurve properties, perhaps implying that rotation rate and magnetic field strength may play important roles in some ordinary looking Type IIP supernova explosions.

Previous studies have shown that the radiation emitted by a rapidly rotating magnetar embedded in a young supernova can greatly amplify its luminosity. These one-dimensional studies have also revealed the existence of an instability arising from the piling up of radiatively accelerated matter in a thin dense shell deep inside the supernova. Here we examine the problem in two dimensions and find that, while instabilities cause mixing and fracture this shell into filamentary structures that reduce the density contrast, the concentration of matter in a hollow shell persists. The extent of the mixing depends upon the relative energy input by the magnetar and the kinetic energy of the inner ejecta. The light curve and spectrum of the resulting supernova will be appreciably altered, as will the appearance of the supernova remnant, which will be shellular and filamentary. A similar pile up and mixing might characterize other events where energy is input over an extended period by a centrally concentrated source, e.g. a pulsar, radioactive decay, a neutrino-powered wind, or colliding shells. The relevance of our models to the recent luminous transient ASASSN-15lh is briefly discussed.

We identify a pre-explosion counterpart to the nearby Type IIP supernova ASASSN-16fq (SN 2016cok) in archival Hubble Space Telescope (HST) data. The source appears to be a blend of several stars that prevents obtaining accurate photometry. However, with reasonable assumptions about the stellar temperature and extinction, the progenitor almost certainly had an initial mass M<17Msun, and was most likely in the mass range 8-12Msun. Observations once ASASSN-16fq has faded will have no difficulty accurately determining the properties of the progenitor. In 8 years of Large Binocular Telescope (LBT) data, no significant progenitor variability is detected to RMS limits of roughly 0.03 mag. Of the six nearby SN with constraints on low level variability, SN 1987A, SN 1993J, SN 2008cn, SN 2011dh, SN 2013ej and ASASSN-16fq, only the slowly fading progenitor of SN 2011dh showed clear evidence of variability. Excluding SN 1987A, the 90% confidence limit implied by these sources on the number of outbursts over the last decade before the SN that last longer than 0.1 years (FWHM) and are brighter than M_R<-8 mag is approximately N<3. Our continuing LBT monitoring program will steadily improve constraints on pre-SN progenitor variability at amplitudes far lower than achievable by SN surveys.

Recent observations have revealed an amazing diversity of extremely luminous supernovae, seemingly increasing in radiant energy without bound. We consider here the physical limits of what existing models can provide for the peak luminosity and total radiated energy for non-relativistic, isotropic stellar explosions. The brightest possible supernova is a Type I explosion powered by a sub-millisecond magnetar. Such models can reach a peak luminosity of $\rm 2\times10^{46}\ erg\ s^{-1}$ and radiate a total energy of $\rm 4 \times10^{52}\ erg$. Other less luminous models are also explored, including prompt hyper-energetic explosions in red supergiants, pulsational-pair instability supernovae, and pair-instability supernovae. Approximate analytic expressions and limits are given for each case. Excluding magnetars, the peak luminosity is near $\rm 1\times10^{44}\ erg\ s^{-1}$ for the brightest models. The corresponding limits on total radiated power are $\rm3 \times 10^{51}\ erg$ (Type I) and $\rm1 \times 10^{51}\ erg$ (Type II). A magnetar-based model for the recent transient event, ASASSN-15lh is presented that strains, but does not exceed the limits of what the model can provide.

Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.

The success or failure of the neutrino-transport mechanism for producing a supernova in an evolved massive star is known to be sensitive not only to the mass of the iron core that collapses, but also to the density gradient in the silicon and oxygen shells surrounding that core. Here we study the systematics of a presupernova core's "compactness" (O'Connor & Ott 2011) as a function of the mass of the star and the physics used in its calculation. Fine-meshed surveys of presupernova evolution are calculated for stars from 15 to 65 Msun. The metallicity and the efficiency of semiconvection and overshoot mixing are both varied and bare carbon-oxygen cores are explored as well as full hydrogenic stars. Two different codes, KEPLER and MESA, are used for the study. A complex interplay of carbon and oxygen burning, especially in shells, can cause rapid variations in the compactness for stars of very nearly the same mass. On larger scales, the distribution of compactness with main sequence mass is found to be robustly non-monotonic, implying islands of "explodability", particularly around 8 to 20 Msun and 25 to 30 Msun. The carbon-oxygen (CO) core mass of a presupernova star is a better, though still ambiguous discriminant of its core structure than the main sequence mass.

We present the results of a search for periodic variables within 4078 time-series light curves and an analysis of the period-color plane for stars in the field of the open cluster NGC2301. One hundred thirty-eight periodic variables were discovered, of which five are eclipsing binary candidates with unequal minima. The remaining 133 periodic variables appear to consist mainly of late-type stars whose variation is due to rotation modulated by star spot activity. The determined periods range from less than a day to over 14 days and have nearly unreddened B-R colors in the range of 0.8 to 2.8. The Barnes (2003) interpretation of the period-color plane of late type stars is tested with our data. Our data did not show distinct I and C sequences, likely due to nonmember field stars contaminating in the background, as we estimate the total contamination to be 43%. Using different assumptions, the gyrochronological age of the cluster is calculated to be 210+/-25 Myr, which falls in the range of age values (164-250 Myr) determined by previous studies. Finally, we present evidence which nullifies the earlier suggestion that two of the variable stars in NGC2301 might be white dwarfs.

Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision observations of stellar parameters to constrain theoretical modeling is also described.

Recently, Dufour et al. (2007) reported the unexpected discovery that a few white dwarfs found in the Sloan Digital Sky Survey had an atmosphere dominated by carbon with little or no trace of hydrogen and helium. Here we present a progress report on these new objects based on new high signal-to-noise follow-up spectroscopic observations obtained at the 6.5m MMT telescope on Mount Hopkins, Arizona.