Is Betelgeuse the Outcome of a Past Merger?

We explore the possibility that the star alpha Orionis (Betelgeuse) is the outcome of a merger that occurred in a low mass ratio (q = M2/M1 = 0.07 - 0.25) binary system some time in the past hundreds of thousands of years. To that goal, we present a simple analytical model to approximate the perturbed internal structure of a post-merger object following the coalescence of a secondary in the mass range 1-4 Msun into the envelope of a 15-17 Msun primary. We then compute the long-term evolution of post-merger objects for a grid of initial conditions and make predictions about their surface properties for evolutionary stages that are consistent with the observed location of Betelgeuse in the Hertzsprung-Russell diagram. We find that if a merger occurred after the end of the primarys main-sequence phase, while it was expanding toward becoming a red supergiant star and typically with radius ~200 - 300 Rsun, then its envelope is spun-up to values which remain in a range consistent with the Betelgeuse observations for thousands of years of evolution. We argue that the best scenario that can explain both the fast rotation of Betelgeuse and its observed large space velocity is one where a binary was dynamically ejected by its parent cluster a few million years ago and then subsequently merged. An alternative scenario in which the progenitor of Betelgeuse was spun up by accretion in a binary and released by the supernova explosion of the companion requires a finely tuned set of conditions but cannot be ruled out.

Synthetic Spectra of Pair-Instability Supernovae in 3D

Pair-Instability Supernovae (PISNe) may signal the deaths of extremely massive stars in the local Universe or massive primordial stars after the end of the Cosmic Dark Ages. Hydrodynamic simulations of these explosions, performed in 1D, 2D, and 3D geometry, have revealed the strong dependence of mixing in the PISN ejecta on dimensionality. This chemical rearrangement is mainly driven by Rayleigh-Taylor instabilities that start to grow shortly after the collapse of the carbon-oxygen core. We investigate the effects of such mixing on the spectroscopic evolution of PISNe by post-processing explosion profiles with the radiation diffusion-equilibrium code SNEC and the implicit Monte Carlo-discrete diffusion Monte Carlo (IMC-DDMC) radiation transport code SuperNu. The first 3D radiation transport calculation of a PISN explosion is presented yielding viewing angle-dependent synthetic spectra and lightcurves. We find that while 2D and 3D mixing does not significantly affect the lightcurves of PISNe, their spectroscopic and color evolution is impacted. Strong features of intermediate mass elements dominated by silicon, magnesium and oxygen appear at different phases and reach different intensities depending on the extent of mixing in the silicon/oxygen interface of the PISN ejecta. On the other hand, we do not find a significant dependence of PISN lightcurves and spectra on viewing angle. Our results showcase the capabilities of SuperNu to handle 3D radiation transport and highlight the importance of modeling time-series of spectra in identifying PISNe with future missions.

Extreme Supernova Models for the Superluminous Transient ASASSN-15lh

The recent discovery of the unprecedentedly superluminous transient ASASSN-15lh (or SN 2015L) with its UV-bright secondary peak challenges all the power-input models that have been proposed for superluminous supernovae. Here we examine some of the few viable interpretations of ASASSN-15lh in the context of a stellar explosion, involving combinations of one or more power inputs. We model the lightcurve of ASASSN-15lh with a hybrid model that includes contributions from magnetar spin-down energy and hydrogen-poor circumstellar interaction. We also investigate models of pure circumstellar interaction with a massive hydrogen-deficient shell and discuss the lack of interaction features in the observed spectra. We find that, as a supernova ASASSN-15lh can be best modeled by the energetic core-collapse of a ~40 Msun star interacting with a hydrogen-poor shell of ~20 Msun. The circumstellar shell and progenitor mass are consistent with a rapidly rotating pulsational pair-instability supernova progenitor as required for strong interaction following the final supernova explosion. Additional energy injection by a magnetar with initial period of 1-2 ms and magnetic field of 0.1-1 x 10^14 G may supply the excess luminosity required to overcome the deficit in single-component models, but this requires more fine-tuning and extreme parameters for the magnetar, as well as the assumption of efficient conversion of magnetar energy into radiation. We thus favor a single-input model where the reverse shock formed in a strong SN ejecta-CSM interaction following a very powerful core-collapse SN explosion can supply the luminosity needed to reproduce the late-time UV-bright plateau.

Finding the First Cosmic Explosions. IV. 90 - 140 M$_{odot}$ Pair-Instability Supernovae

Population III stars that die as pair-instability supernovae are usually thought to fall in the mass range of 140 - 260 M$_{odot}$. But several lines of work have now shown that rotation can build up the He cores needed to encounter the pair instability at stellar masses as low as 90 $_{odot}$. Depending on the slope of the initial mass function of Population III stars, there could be 4 - 5 times as many stars from 90 - 140 $_{odot}$ in the primordial universe than in the usually accepted range. We present numerical simulations of the pair-instability explosions of such stars performed with the MESA, FLASH and RAGE codes. We find that they will be visible to supernova factories such as Pan-STARRS and LSST in the optical out to z $sim$ 1 - 2 and to JWST and the 30 m-class telescopes in the NIR out to $z sim$ 7 - 10. Such explosions will thus probe the stellar populations of the first galaxies and cosmic star formation rates in the era of cosmological reionization. These supernovae are also easily distinguished from more massive pair-instability explosions, underscoring the fact that there is far greater variety to the light curves of these events than previously understood.

Analytical Light Curve Models of Super-Luminous Supernovae: chi^2-Minimizations of Parameter Fits

We present fits of generalized semi-analytic supernova (SN) light curve (LC) models for a variety of power inputs including Ni-56 and Co-56 radioactive decay, magnetar spin-down, and forward and reverse shock heating due to supernova ejecta-circumstellar matter (CSM) interaction. We apply our models to the observed LCs of the H-rich Super Luminous Supernovae (SLSN-II) SN 2006gy, SN 2006tf, SN 2008am, SN 2008es, CSS100217, the H-poor SLSN-I SN 2005ap, SCP06F6, SN 2007bi, SN 2010gx and SN 2010kd as well as to the interacting SN 2008iy and PTF09uj. Our goal is to determine the dominant mechanism that powers the LCs of these extraordinary events and the physical conditions involved in each case. We also present a comparison of our semi-analytical results with recent results from numerical radiation hydrodynamics calculations in the particular case of SN 2006gy in order to explore the strengths and weaknesses of our models. We find that CS shock heating produced by ejecta-CSM interaction provides a better fit to the LCs of most of the events we examine. We discuss the possibility that collision of supernova ejecta with hydrogen-deficient CSM accounts for some of the hydrogen-deficient SLSNe (SLSN-I) and may be a plausible explanation for the explosion mechanism of SN 2007bi, the pair-instability supernova (PISN) candidate. We characterize and discuss issues of parameter degeneracy.

Hydrogen-Poor Circumstellar Shells from Pulsational Pair-Instability Supernovae with Rapidly Rotating Progenitors

In certain mass ranges, massive stars can undergo a violent pulsation triggered by the electron/positron pair instability that ejects matter, but does not totally disrupt the star. After one or more of these pulsations, such stars are expected to undergo core-collapse to trigger a supernova explosion. The mass range susceptible to this pulsational phenomena may be as low as 50-70 Msun if the progenitor is of very low metallicity and rotating sufficiently rapidly to undergo nearly homogeneous evolution. The mass, dynamics, and composition of the matter ejected in the pulsation are important aspects to determine the subsequent observational characteristics of the explosion. We examine the dynamics of a sample of stellar models and rotation rates and discuss the implications for the first stars, for LBV-like phenomena, and for superluminous supernovae. We find that the shells ejected by pulsational pair-instability events with rapidly rotating progenitors (>30% the critical value) are hydrogen-poor and helium and oxygen-rich.

SN 2006oz: rise of a super-luminous supernova observed by the SDSS-II SN Survey

We study SN 2006oz, a newly-recognized member of the class of H-poor, super-luminous supernovae. We present multi-color light curves from the SDSS-II SN Survey, that cover the rise time, as well as an optical spectrum that shows that the explosion occurred at z~0.376. We fitted black body functions to estimate the temperature and radius evolution of the photosphere and used the parametrized code SYNOW to model the spectrum. We constructed a bolometric light curve and compared it with explosion models. The very early light curves show a dip in the g- and r-bands and a possible initial cooling phase in the u-band before rising to maximum light. The bolometric light curve shows a precursor plateau with a duration of 6-10 days in the rest-frame. A lower limit of M_u < -21.5 can be placed on the absolute peak luminosity of the SN, while the rise time is constrained to be at least 29 days. During our observations, the emitting sphere doubled its radius to 2x10^15 cm, while the temperature remained hot at 15000 K. As for other similar SNe, the spectrum is best modeled with elements including O II and Mg II, while we tentatively suggest that Fe III might be present. We suggest that the precursor plateau might be related to a recombination wave in a circumstellar medium (CSM) and discuss whether this is a common property of all similar explosions. The subsequent rise can be equally well described by input from a magnetar or by ejecta-CSM interaction, but the models are not well constrained owing to the lack of post-maximum observations, and CSM interaction has difficulties accounting for the precursor plateau self-consistently. Radioactive decay is less likely to be the mechanism that powers the luminosity. The host galaxy, detected in deep imaging with the 10 m GTC, is a moderately young and star-forming, but not a starburst, galaxy. It has an absolute magnitude of M_g = -16.9.

The long-term X-ray variability properties of AGN in the Lockman Hole region

We present the results from a detailed X-ray variability analysis of 66 AGN in the Lockman Hole, which have optical spectroscopic identifications. We compare, quantitatively, their variability properties with the properties of local AGN, and we study the variability-luminosity relation as a function of redshift, and the variability-redshift relation in two luminosity bins. We use archival data from the last 10 XMM observations of the Lockman Hole field to extract light curves in the rest frame, 2-10 keV band. We use the normalized excess variance to quantify the variability amplitude. Using the latest results regarding the AGN power spectral shape and its dependence on black hole mass and accretion rate, we are able to compute model variability-luminosity curves, which we compare with the relations we observe. When we consider all the sources in our sample, we find that their variability amplitude decreases with increasing redshift and luminosity. These global anti-correlations are less pronounced when we split the objects in various luminosity and redshift bins. We do not find a significant correlation between variability amplitude and spectral slope. The variability-luminosity relation that we detect has a larger amplitude when compared to that of local AGN. We also find that, at a given luminosity, the variability amplitude increases with redshift up to z~1, and then stays roughly constant. Our results imply that the AGN X-ray mechanism operates in the same way at all redshifts. Among objects with the same luminosity in our sample, the black hole mass decreases and the accretion rate increases with larger redshift.