The R Coronae Borealis (RCB) stars are extremely hydrogen-deficient carbon stars which produce large amounts of dust, causing sudden deep declines in brightness. They are believed to be formed primarily through white dwarf mergers. In this paper, we use MESA to investigate how post-merger objects with a range of He-burning shell temperatures from 2.1 - 5.4 x 10^8 K with solar and subsolar metallicities evolve into RCB stars. The most successful model of these has subsolar metallicity and a temperature near 3 x 10^8 K. We find a strong dependence on He-burning shell temperature for surface abundances of elements involved in the CNO cycle, as well as differences in effective temperature and radius of RCBs. We find around 1 dex diminished surface abundances for the models with 10% subsolar metallicity, which is expected. Models with subsolar metallicities also exhibit longer lifetimes than their solar counterparts. Additionally, we find that convective mixing of the burned material occurs only in the first few years of post-merger evolution, after which the surface abundances are constant during and after the RCB phase, providing evidence for why these stars show a strong enhancement of partial He-burning products.

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 primary's main-sequence phase, while it was expanding toward becoming a red supergiant star and typically with radius ~200 - 300 Rsun, then it's 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.

This paper presents data and analysis of SN 2010kd, a low-redshift ($z = 0.101$) H-deficient superluminous supernova (SLSN), based on ultraviolet/optical photometry and optical spectroscopy spanning between $-$28 and +194 days relative to $\mathit{B}$ band maximum light. The $\mathit{B}$ band light curve comparison of SN 2010kd with a subset of well-studied SLSNe I at comparable redshifts indicates that it is a slow-decaying PTF12dam like SLSN. Analytical light-curve modeling using the $\mathtt{Minim}$ code suggests that the bolometric light curve of SN 2010kd favors circumstellar matter interaction for the powering mechanism. $\mathtt{SYNAPPS}$ modeling of the early-phase spectra does not identify broad H or He lines, whereas the photospheric-phase spectra are dominated by O I, O II, C II, C IV and Si II, particularly, presence of both low and high-velocity components of O II and Si II lines. The nebular-phase spectra of SN 2010kd are dominated by O I and Ca II emission lines similar to those seen in other SLSNe I. The line velocities in SN 2010kd exhibit flatter evolution curves similar to SN 2015bn but with comparatively higher values. SN 2010kd shows a higher single-zone local thermodynamic equilibrium temperature in comparison to PTF12dam and SN 2015bn, and it has an upper O I ejected mass limit of $\sim 10~M_\odot$. The host of SN 2010kd is a dwarf galaxy with a high star-formation rate ($\sim 0.18 \pm 0.04~M_\odot$ yr$^{-1}$) and extreme emission lines.

We present the newly-incorporated gray radiation hydrodynamics capabilities of the FLASH code based on a radiation flux-limiter aware hydrodynamics numerical implementation designed specifically for applications in astrophysical problems. The newly incorporated numerical methods consist of changes in the unsplit hydrodynamics solver and adjustments in the flux-limited radiation diffusion unit. Our method can treat problems in both the strong and weak radiation-matter coupling limits as well as transitions between the two regimes. Appropriate extensions in the "Helmholtz" equation of state are implemented to treat two-temperature astrophysical plasmas involving the interaction between radiation and matter and the addition of a new opacity unit based on the OPAL opacity database, commonly used for astrophysical fluids. A set of radiation-hydrodynamics test problems is presented aiming to showcase the new capabilities of FLASH and to provide direct comparison to other codes like CASTRO. To illustrate the capacity of FLASH to simulate phenomena occurring in stellar explosions, such as shock break-out, radiative precursors and supernova ejecta heating due to the decays of radioactive nickel-56 and cobalt-56, we also present 1D supernova simulations and compare the computed lightcurves to those of the SNEC code. The latest public release of FLASH with these enhanced capabilities is available for download and use by the broader astrophysics community.

The R Coronae Borealis (RCB) stars are rare hydrogen--deficient, carbon--rich supergiants. They undergo extreme, irregular declines in brightness of many magnitudes due to the formation of thick clouds of carbon dust. It is thought that RCB stars result from the mergers of CO/He white dwarf (WD) binaries. We constructed post--merger spherically asymmetric models computed with the MESA code, and then followed the evolution into the region of the HR diagram where the RCB stars are located. We also investigated nucleosynthesis in the dynamically accreting material of CO/He WD mergers which may provide a suitable environment for significant production of 18O and the very low 16O/18O values observed. We have also discovered that the N abundance depends sensitively on the peak temperature in the He--burning shell. Our MESA modeling consists of engineering the star by adding He--WD material to an initial CO--WD model, and then following the post--merger evolution using a nuclear--reaction network to match the observed RCB abundances as it expands and cools to become an RCB star. These new models are more physical because they include rotation, mixing, mass-loss, and nucleosynthesis within MESA. We follow the later evolution beyond the RCB phase to determine the stars' likely lifetimes. The relative numbers of known RCB and Extreme Helium (EHe) stars correspond well to the lifetimes predicted from the MESA models. In addition, most of computed abundances agree very well with the observed range of abundances for the RCB class.

Superluminous supernova (SLSN) lightcurves exhibit a superior diversity compared to their regular luminosity counterparts in terms of rise and decline timescales, peak luminosities and overall shapes. It remains unclear whether this striking variety arises due to a dominant power input mechanism involving many underlying parameters, or due to contributions by different progenitor channels. In this work, we propose that a systematic quantitative study of SLSN lightcurve timescales and shape properties, such as symmetry around peak luminosity, can be used to characterize these enthralling stellar explosions. We find that applying clustering analysis on the properties of model SLSN lightcurves, powered by either a magnetar spin-down or a supernova ejecta-circumstellar interaction mechanism, can yield a distinction between the two, especially in terms of lightcurve symmetry. We show that most events in the observed SLSN sample with well-constrained lightcurves and early detections strongly associate with clusters dominated by circumstellar interaction models. Magnetar spin-down models also show association at a lower degree but have difficulty in reproducing fast-evolving and fully symmetric lightcurves. We believe this is due to the truncated nature of the circumstellar interaction shock energy input as compared to decreasing but continuous power input sources like magnetar spin-down and radioactive Ni-56 decay. Our study demonstrates the importance of clustering analysis in characterizing SLSNe based on high-cadence photometric observations that will be made available in the near future by surveys like LSST, ZTF and Pan-STARRS.

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.

We present photometry, spectra, and spectropolarimetry of supernova (SN) 2012ab, mostly obtained over the course of $\sim 300$ days after discovery. SN 2012ab was a Type IIn (SN IIn) event discovered near the nucleus of spiral galaxy 2MASXJ12224762+0536247. While its light curve resembles that of SN 1998S, its spectral evolution does not. We see indications of CSM interaction in the strong intermediate-width emission features, the high luminosity (peak at absolute magnitude $M=-19.5$), and the lack of broad absorption features in the spectrum. The H$\alpha$ emission undergoes a peculiar transition. At early times it shows a broad blue emission wing out to $-14{,}000$ km $\mathrm{s^{-1}}$ and a truncated red wing. Then at late times ($>$ 100$\,$days) it shows a truncated blue wing and a very broad red emission wing out to roughly $+20{,}000$ km $\mathrm{s^{-1}}$. This late-time broad red wing probably arises in the reverse shock. Spectra also show an asymmetric intermediate-width H$\alpha$ component with stronger emission on the red side at late times. The evolution of the asymmetric profiles requires a density structure in the distant CSM that is highly aspherical. Our spectropolarimetric data also suggest asphericity with a strong continuum polarization of $\sim 1-3$% and depolarization in the H$\alpha$ line, indicating asphericity in the CSM at a level comparable to that in other SNe IIn. We estimate a mass-loss rate of $\dot{M} = 0.050\, {\rm M}_{\odot}\,\mathrm{yr^{-1}}$ for $v_{\rm pre} = 100$$\,$km$\,$$\mathrm{s^{-1}}$ extending back at least 75$\,$yr prior to the SN. The strong departure from axisymmetry in the CSM of SN 2012ab may suggest that the progenitor was an eccentric binary system undergoing eruptive mass loss.

We explore simple semi-analytic fits to the bolometric light curve of Gaia17biu/SN 2017egm, the most nearby hydrogen-deficient superluminous supernova (SLSN I) yet discovered. SN 2017egm has a quasi-bolometric light curve that is uncharacteristic of other SLSN I by having a nearly linear rise to maximum and decline from peak, with a very sharp transition. Magnetar models have difficulty explaining the sharp peak and may tend to be too bright 20 d after maximum. Light curves powered only by radioactive decay of nickel fail on similar grounds and because they demand greater nickel mass than ejecta mass. Simple models based on circumstellar interaction do have a sharp peak corresponding to the epoch when the forward shock breaks out of the optically-thick circumstellar medium or the reverse shock reaches the inside of the ejecta. We find that models based on circumstellar interaction with a constant-density shell provide an interesting fit to the bolometric light curve from 15 d before to 15 d after peak light of SN 2017egm and that both magnetar and radioactive decay models fail to fit the sharp peak. Future photometric observations should easily discriminate basic CSI models from basic magnetar models. The implications of a CSI model are briefly discussed.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

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.

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

Pair Instability Supernovae have been suggested as candidates for some Super Luminous Supernovae, such as SN 2007bi, and as one of the dominant types of explosion occurring in the early Universe from massive, zero-metallicity Population III stars. The progenitors of such events can be rapidly rotating, therefore exhibiting different evolutionary properties due to the effects of rotationally-induced mixing and mass-loss. Proper identification of such events requires rigorous radiation hydrodynamics and radiative transfer calculations that capture not only the behavior of the light curve but also the spectral evolution of these events. We present radiation hydrodynamics and radiation transport calculations for 90-300 Msun rotating pair-instability supernovae covering both the shock break-out and late light curve phases. We also investigate cases of different initial metallicity and rotation rate to determine the impact of these parameters on the detailed spectral characteristics of these events. In agreement with recent results on non-rotating pair instability supernovae, we find that for a range of progenitor masses and rotation rates these events have intrinsically red colors in contradiction with observations of super-luminous supernovae. The spectroscopic properties of rotating pair instability supernovae are similar to those of non-rotating events with stripped hydrogen and helium envelopes. We find that the progenitor metallicity and rotation rate properties are erased after the explosion and cannot be identified in the resulting model spectra. It is the combined effects of pre-supernova mass-loss and the basic properties of the supernova ejecta such as mass, temperature and velocity that have the most direct impact in the model spectra of pair instability supernovae.

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

The magnetorotational instability (MRI) is key physics in accretion disks and is widely considered to play some role in massive-star core collapse. Models of rotating massive stars naturally develop very strong shear at composition boundaries, a necessary condition for MRI instability, and the MRI is subject to triply-diffusive destabilizing effects in radiative regions. We have used the MESA stellar evolution code to compute magnetic effects due to the Spruit-Taylor mechanism and the MRI, separately and together, in a sample of massive star models. We find that the MRI can be active in the later stages of massive star evolution, leading to mixing effects that are not captured in models that neglect the MRI. The MRI and related magneto-rotational effects can move models of given ZAMS mass across "boundaries" from degenerate CO cores to degenerate O/Ne/Mg cores and from degenerate O/Ne/Mg cores to iron cores, thus affecting the final evolution and the physics of core collapse. The MRI acting alone can slow the rotation of the inner core in general agreement with the observed "initial" rotation rates of pulsars. The MRI analysis suggests that localized fields ~10^{12} G may exist at the boundary of the iron core. With both the ST and MRI mechanisms active in the 20 solar mass model, we find that the helium shell mixes entirely out into the envelope. Enhanced mixing could yield a population of yellow or even blue supergiant supernova progenitors that would not be standard SN IIP.

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.

We present follow-up observations of an optical transient (OT) discovered by ROTSE on Jan. 21, 2009. Photometric monitoring was carried out with ROTSE-IIIb in the optical and Swift in the UV up to +70 days after discovery. The light curve showed a fast rise time of ~10 days followed by a steep decline over the next 60 days, which was much faster than that implied by 56Ni - 56Co radioactive decay. The SDSS DR10 database contains a faint, red object at the position of the OT, which appears slightly extended. This and other lines of evidence suggest that the OT is of extragalactic origin, and this faint object is likely the host galaxy. A sequence of optical spectra obtained with the 9.2-m Hobby-Eberly Telescope (HET) between +8 and +45 days after discovery revealed a hot, blue continuum with no visible spectral features. A few weak features that appeared after +30 days probably originated from the underlying host. Fitting synthetic templates to the observed spectrum of the host galaxy revealed a redshift of z = 0.19. At this redshift the peak magnitude of the OT is close to -22.5, similar to the brightest super-luminous supernovae; however, the lack of identifiable spectral features makes the massive stellar death hypothesis less likely. A more plausible explanation appears to be the tidal disruption of a sun-like star by the central super-massive black hole. We argue that this transient likely belongs to a class of super-Eddington tidal disruption events.

We apply the mathematical formalism of vector spherical harmonics decomposition to convective stellar velocity fields from multi-dimensional hydrodynamics simulations, and show that the resulting power spectra furnish a robust and stable statistical description of stellar convective turbulence. Analysis of the power spectra help identify key physical parameters of the convective process such as the dominant scale of the turbulent motions that influence the structure of massive evolved pre-supernova stars. We introduce the numerical method that can be used to calculate vector spherical harmonics power spectra from 2D and 3D convective shell simulation data. Using this method we study the properties of oxygen shell burning and convection for a 15 Msun star simulated by the hydrodynamics code FLASH in 2D and 3D. We discuss the importance of realistic initial conditions to achieving successful core-collapse supernova explosions in multi-dimensional simulations. We show that the calculated power spectra can be used to generate realizations of the velocity fields of pre-supernova convective shells. We find that the slope of the solenoidal mode power spectrum remains mostly constant throughout the evolution of convection in the oxygen shell in both 2D and 3D simulations. We also find that the characteristic radial scales of the convective elements are smaller in 3D than in 2D while the angular scales are larger in 3D.

We study the effects of rotation on the dynamics, energetics and Ni-56 production of Pair Instability Supernova explosions by performing rotating two-dimensional ("2.5-D") hydrodynamics simulations. We calculate the evolution of eight low metallicity (Z = 10^-3, 10^-4 Zsun) massive (135-245 Msun) PISN progenitors with initial surface rotational velocities 50% that of the critical Keplerian value using the stellar evolution code MESA. We allow for both the inclusion and the omission of the effects of magnetic fields in the angular momentum transport and in chemical mixing, resulting in slowly-rotating and rapidly-rotating final carbon-oxygen cores, respectively. Increased rotation for carbon-oxygen cores of the same mass and chemical stratification leads to less energetic PISN explosions that produce smaller amounts of Ni-56 due to the effect of the angular momentum barrier that develops and slows the dynamical collapse. We find a non-monotonic dependence of Ni-56 production on rotational velocity in situations when smoother composition gradients form at the outer edge of the rotating cores. In these cases, the PISN energetics are determined by the competition of two factors: the extent of chemical mixing in the outer layers of the core due to the effects of rotation in the progenitor evolution and the development of angular momentum support against collapse. Our 2.5-D PISN simulations with rotation are the first presented in the literature. They reveal hydrodynamic instabilities in several regions of the exploding star and increased explosion asymmetries with higher core rotational velocity.

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.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

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.

Many population synthesis and stellar evolution studies have addressed the evolution of close binary systems in which the primary is a compact remnant and the secondary is filling its Roche lobe, thus triggering mass transfer. Although tidal locking is expected in such systems, most studies have neglected the rotationally-induced mixing that may occur. Here we study the possible effects of mixing in the mass-losing stars for a range in secondary star masses and metallicities. We find that tidal locking can induce rotational mixing prior to contact and thus affect the evolution of the secondary star if the effects of the Spruit-Tayler dynamo are included both for angular momentum and chemical transport. Once contact is made, the effect of mass transfer tends to be more rapid than the evolutionary time scale, so the effects of mixing are no longer directly important, but the mass transfer strips matter to inner layers that may have been affected by the mixing. These effects are enhanced for secondaries of 1-1.2 Msun and for lower metallicities. We discuss the possible implications for the paucity of carbon in the secondaries of the cataclysmic variable SS Cyg and the black hole candidate XTE J1118+480 and for the progenitor evolution of Type Ia supernovae. We also address the issue of the origin of blue straggler stars in globular and open clusters. We find that for models that include rotation consistent with that observed for some blue straggler stars, evolution is chemically homogeneous. This leads to tracks in the HR diagram that are brighter and bluer than the non-rotating main-sequence turn-off point. Rotational mixing could thus be one of the factors that contribute to the formation of blue stragglers.

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 issue of which stars may reach the conditions of electron/positron pair formation instability is of importance to understand the final evolution both of the first stars and of contemporary stars. The criterion to enter the pair instability regime in density and temperature is basically controlled by the mass of the oxygen core. The main sequence masses that produce a given oxygen core mass are, in turn, dependent on metallicity, mass loss, and convective and rotationally-induced mixing. We examine the evolution of massive stars to determine the minimum main sequence mass that can encounter pair-instability effects, either a pulsational pair instability (PPISN) or a full-fledged pair-instability supernova (PISN). We concentrate on zero-metallicity stars with no mass loss subject to the Schwarzschild criterion for convective instability, but also explore solar metallicity and mass loss and the Ledoux criterion. As expected, for sufficiently strong rotationally-induced mixing, the minimum main sequence mass is encountered for conditions that induce effectively homogeneous evolution such that the original mass is converted almost entirely to helium and then to oxygen. For this case, we find that the minimum main sequence mass is ~40 Msun to encounter PPISN and ~65 Msun to encounter a PISN. When mass-loss is taken into account those mass limits become ~50 Msun for PPISN and ~80 Msun for PISN progenitors. The implications of these results for the first stars and for contemporary supernovae is discussed.

We present generalized supernova (SN) light curve (LC) models for a variety of power inputs. We provide an expression for the power input that is produced by self-similar forward and reverse shocks in SN ejecta - circumstellar matter (CSM) interaction. We find that this ejecta-CSM interaction luminosity is in agreement with results from multi-dimensional radiation hydrodynamics simulations in the optically-thin case. We develop a model for the case of an optically-thick CSM by invoking an approximation for the effects of radiative diffusion. In the context of this model, we provide predictions for the time of forward shock break-out from the optically-thick part of the CSM envelope. We also introduce a hybrid LC model that incorporates ejecta-CSM interaction plus Ni-56 and Co-56 radioactive decay input. We fit this hybrid model to the LC of the Super-Luminous Supernova (SLSN) 2006gy. We find that this model provides a better fit to the LC of this event than previously presented models. We also address the relation between Type IIL and Type IIn SN with ejecta-CSM interaction models. Forward and reverse shock power input due to CSM interaction can produce the LCs of Type IIn SNe in terms of duration, shape and decline rate. This model can also produce LCs that are symmetric in shape around peak luminosity. We conclude that the observed LC variety of SNe Type IIn and of the SLSNe is likely to be a byproduct of the large range of conditions relevant to significant ejecta-CSM interaction as a power source.

We present observations and interpretation of the Type IIn supernova SN 2008am discovered by the ROTSE Supernova Verification Project (RSVP). SN 2008am peaked at approximately -22.3 mag at a redshift of z=0.2338, giving it a peak luminosity of 3 x 10^{44}erg/s and making it one of the most luminous supernovae ever observed. The total radiated energy is ~ 2 x 10^{51} erg. Photometric observations in the ultraviolet, optical and infrared bands (J,H,Ks) constrain the SED evolution. We obtained six optical spectra of the supernova, five on the early decline from maximum light and a sixth nearly a year later plus a very late-time spectrum (~2 yr) of the host galaxy. The spectra of SN 2008am show strong Balmer-line and He I lambda 5876A emission with intermediate widths (~25A) in the first ~40 days after optical maximum. We examine a variety of models for the line wings and conclude that multiple scattering is most likely, implying that our spectra contain no specific information on the bulk flow velocity. We examine a variety of models for the ROTSE light curve subject to the rise time and the nature of the spectra, including radioactive decay, shocks in optically-thick and optically-thin circumstellar media (CSM) and a magnetar. The most successful model is one for which the CSM is optically-thick and in which diffusion of forward shock-deposited luminosity gives rise to the observed light curve. Diffusion of the shock-deposited energy from the forward shock is found to be important to account for the rising part of the light curve. Although there are differences in detail, SN 2008am appears to be closely related to other super-luminous Type IIn supernovae, SN 2006gy, SN 2006tf and perhaps SN 2008iy, that may represent the deaths of very massive LBV-type progenitors and for which the luminosity is powered by the interaction of the ejecta with a dense circumstellar medium.

SN 2007if was the third over-luminous SN Ia detected after 2003fg and 2006gz. We present the photometric and spectroscopic observations of the supernova and its host by ROTSE-III, HET and Keck. From the H_alpha line identified in the host spectra, we determine a redshift of 0.0736. At this distance, the supernova reached an absolute magnitude of -20.4, brighter than any other SNe Ia ever observed. If the source of luminosity is radioactive decay, a large amount of radioactive nickel (~1.5 solar masses) is required to power the peak luminosity, more than can be produced realistically in a Chandrasekhar mass progenitor. Low expansion velocity, similar to that of 2003fg, is also measured around the maximum light. The observations may suggest that SN 2007if was from a massive white dwarf progenitor, plausibly exploding with mass well beyond 1.4 solar masses. Alternatively, we investigate circumstellar interaction that may contribute to the excess luminosity.

We consider simple models based on core collapse or pair-formation supernovae to account for the light curve of the transient SCP06F6. A radioactive decay diffusion model provides estimates of the mass of the required radioactive nickel and the ejecta as functions of the unknown redshift. An opacity change such as by dust formation or a recombination front may account for the rapid decline from maximum. We particularly investigate two specific redshifts: $z=0.143$, for which Gaensicke et al. (2008) have proposed that the unidentified broad absorption features in the spectrum of SCP06F6 are C$_{2}$ Swan bands, and $z=0.57$ based on a crude agreement with the Ca H&K and UV iron-peak absorption features that are characteristic of supernovae of various types. The ejected masses and kinetic energies are smaller for a more tightly constrained model invoking envelope recombination. We also discuss the possibilities of circumstellar matter (CSM) shell diffusion and shock interaction models. In general, optically-thick CSM diffusion models can fit the data with the underlying energy coming from an energetic buried supernova. Models in which the CSM is of lower density so that the shock energy is both rapidly thermalized and radiated tend not to be self-consistent. We suggest that a model of SCP06F6 worth futher exploration is one in which the redshift is $\sim$ 0.57, the spectral features are Ca and iron peak elements, and the light curve is powered by the diffusive release of a substantial amount of energy from nickel decay or from an energetic supernova buried in the ejecta of an LBV-like event.

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.

We present a study of the mid-infrared properties and dust content of a sample of 27 HII ``blobs'', a rare class of compact HII regions in the Magellanic Clouds. A unique feature of this sample is that even though these HII regions are of high and low excitation they have nearly the same physical sizes ~1.5-3 pc. We base our analysis on archival 3-8 microns infrared imagery obtained with the Infrared Array Camera (IRAC) on board the Spitzer Space Telescope. We find that despite their youth, sub-solar metallicity and varied degrees of excitation, the mid-infrared colors of these regions are similar to those of typical HII regions. Higher excitation ``blobs'' (HEBs) display stronger 8 micron emission and redder colors than their low-excitation counterparts (LEBs).