Interacting Supernovae and Supernova Impostors: SN 2009ip, is this the End?

We report the results of a three-year-long dedicated monitoring campaign of a restless luminous blue variable (LBV) in NGC 7259. The object, named SN 2009ip, was observed photometrically and spectroscopically in the optical and near-infrared domains. We monitored a number of erupting episodes in the past few years, and increased the density of our observations during eruptive episodes. In this paper, we present the full historical data set from 2009 to 2012 with multi-wavelength dense coverage of the two high-luminosity events between 2012 August and September. We construct bolometric light curves and measure the total luminosities of these eruptive or explosive events. We label them the 2012a event (lasting ~50 days) with a peak of 3 × 1041 erg s-1, and the 2012b event (14 day rise time, still ongoing) with a peak of 8 × 1042 erg s-1. The latter event reached an absolute R-band magnitude of about -18, comparable to that of a core-collapse supernova (SN). Our historical monitoring has detected high-velocity spectral features (~13,000 km s-1) in 2011 September, one year before the current SN-like event. This implies that the detection of such high-velocity outflows cannot, conclusively, point to a core-collapse SN origin. We suggest that the initial peak in the 2012a event was unlikely to be due to a faint core-collapse SN. We propose that the high intrinsic luminosity of the latest peak, the variability history of SN 2009ip, and the detection of broad spectral lines indicative of high-velocity ejecta are consistent with a pulsational pair-instability event, and that the star may have survived the last outburst. The question of the survival of the LBV progenitor star and its future fate remain open issues, only to be answered with future monitoring of this historically unique explosion.

Ultraluminous Supernovae as a New Probe of the Interstellar Medium in Distant Galaxies

We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultraluminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z = 1.566 with a peak brightness of M UV ≈ -22.3 mag. PS1-11bam is one of the highest redshift spectroscopically confirmed SNe known to date. The spectrum exhibits broad absorption features typical of previous ULSNe (e.g., C II, Si III), and strong and narrow Mg II and Fe II absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [O II]λ3727 emission line at the same redshift. The equivalent widths of the Fe II λ2600 and Mg II λ2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of τ* ≈ 15-45 Myr and a stellar mass of M * ≈ (1.1-2.6) × 109 M ⊙ (for Z = 0.05-1 Z ⊙). The star formation rate inferred from the UV continuum and [O II]λ3727 emission line is ≈10 M ⊙ yr-1, higher than in previous ULSN hosts. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the ISM in distant galaxies. The depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z ~ 1-2); the future combination of LSST and 30 m class telescopes promises to extend this technique to z ~ 4.

An ultraviolet-optical flare from the tidal disruption of a helium-rich stellar core

The flare of radiation from the tidal disruption and accretion of a star can be used as a marker for supermassive black holes that otherwise lie dormant and undetected in the centres of distant galaxies. Previous candidate flares have had declining light curves in good agreement with expectations, but with poor constraints on the time of disruption and the type of star disrupted, because the rising emission was not observed. Recently, two `relativistic' candidate tidal disruption events were discovered, each of whose extreme X-ray luminosity and synchrotron radio emission were interpreted as the onset of emission from a relativistic jet. Here we report a luminous ultraviolet-optical flare from the nuclear region of an inactive galaxy at a redshift of 0.1696. The observed continuum is cooler than expected for a simple accreting debris disk, but the well-sampled rise and decay of the light curve follow the predicted mass accretion rate and can be modelled to determine the time of disruption to an accuracy of two days. The black hole has a mass of about two million solar masses, modulo a factor dependent on the mass and radius of the star disrupted. On the basis of the spectroscopic signature of ionized helium from the unbound debris, we determine that the disrupted star was a helium-rich stellar core.

Type Ia supernovae from exploding oxygen-neon white dwarfs

Context. The progenitor problem of Type Ia supernovae (SNe Ia) is still unsolved. Most of these events are thought to be explosions of carbon-oxygen (CO) white dwarfs (WDs), but for many of the explosion scenarios, particularly those involving the externally triggered detonation of a sub-Chandrasekhar mass WD (sub-M-Ch, WD), there is also a possibility of having an oxygen-neon (ONe) WD as progenitor.

Aims. We simulate detonations of ONe WDs and calculate synthetic observables from these models. The results are compared with detonations in CO WDs of similar mass and observational data of SNe Ia.

Methods. We perform hydrodynamic explosion simulations of detonations in initially hydrostatic ONe WDs for a range of masses below the Chandrasekhar mass (M-Ch), followed by detailed nucleosynthetic postprocessing with a 384-isotope nuclear reaction network. The results are used to calculate synthetic spectra and light curves, which are then compared with observations of SNe Ia. We also perform binary evolution calculations to determine the number of SNe Ia involving ONe WDs relative to the number of other promising progenitor channels.

Results. The ejecta structures of our simulated detonations in sub-M-Ch, ONe WDs are similar to those from CO WDs. There are, however, small systematic deviations in the mass fractions and the ejecta velocities. These lead to spectral features that are systematically less blueshifted. Nevertheless, the synthetic observables of our ONe WD explosions are similar to those obtained from CO models.

Conclusions. Our binary evolution calculations show that a significant fraction (3-10%) of potential progenitor systems should contain an ONe WD. The comparison of our ONe models with our CO models of comparable mass (similar to 1.2 M-circle dot) shows that the less blueshifted spectral features fit the observations better, although they are too bright for normal SNe Ia.

5.9-keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kα line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two three-dimensional explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two (1.1 and 0.9 M⊙) white dwarfs. Both models are based on solar metallicity zero-age main-sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesizes ˜3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kα line flux of the delayed-detonation model exceeds that of the merger model by a factor of ˜4.5. Since in both models the 5.9-keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. Of the currently existing telescopes, XMM-Newton/pn is the best instrument for close (≲1-2 Mpc), non-background limited SNe Ia because of its large effective area. Due to its low instrumental background, Chandra/ACIS is currently the best choice for SNe Ia at distances above ˜2 Mpc. For the delayed-detonation scenario, a line detection is feasible with Chandra up to ˜3 Mpc for an exposure time of 106 s. We find that it should be possible with currently existing X-ray instruments (with exposure times ≲5 × 105 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of ˜5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (˜6.4 Mpc) and SN 2014J (˜3.5 Mpc).

SN2014J gamma rays from the 56Ni decay chain

Context. The detection and measurement of gamma-ray lines from the decaychain of 56Ni provides unique information about the explosionin supernovae. SN2014J at 3.3 Mpc is a sufficiently-nearby supernova oftype Ia so that such measurements have been feasible with the gamma-rayspectrometer SPI on ESA's INTEGRAL gamma-ray observatory.

Aims:The 56Ni freshly produced in the supernova is understood topower the optical light curve, because it emits gamma rays upon itsradioactive decay first to 56Co and then to 56Fe.Gamma-ray lines from 56Co decay are expected to becomedirectly visible through the white dwarf material several weeks afterthe explosion, as they progressively penetrate the overlying material ofthe supernova envelope, which is diluted as it expands. The lines areexpected to be Doppler-shifted or broadened from the kinematics of the56Ni ejecta. We aim to exploit high-resolution gamma-rayspectroscopy with the SPI spectrometer on INTEGRAL toward constrainingthe 56Ni distribution and kinematics in this supernova.

Methods: We use the observations with the SPI spectrometer onINTEGRAL, together with an improved instrumental background method.

Results: We detect the two main lines from 56Co decay at847 and 1238 keV, which are significantly Doppler-broadened, and atintensities (3.65 ± 1.21) × 10-4 and (2.27± 0.69) × 10-4 ph cm-2s-1, respectively, at their brightness maximum. We measuretheir rise toward a maximum after about 60-100 days and a declinethereafter. The intensity ratio of the two lines is found to beconsistent with expectations from 56Co decay (0.62 ±0.28 at brightness maximum, the expected ratio is 0.68). We find thatthe broad lines seen in the late, gamma-ray transparent phase are notrepresentative of the early gamma-ray emission, and notice instead thatthe emission spectrum is complex and irregular until the supernova isfully transparent to gamma rays, with progressive uncovering of the bulkof 56Ni. We infer that the explosion morphology is notspherically symmetric, both in the distribution of 56Ni andin the unburnt material which occults the 56Co emission.After we compare light curves from different plausible models, theresulting 56Ni mass is determined to be 0.49 ± 0.09M⊙.

The white dwarf's carbon fraction as a secondary parameter of Type Ia supernovae

Context. Binary stellar evolution calculations predict thatChandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radiallyvarying profile for the composition with a carbon depleted core. Manyrecent multi-dimensional simulations of Type Ia supernovae (SNe Ia),however, assume the progenitor WD has a homogeneous chemicalcomposition.
Aims: In this work, we explore the impact ofdifferent initial carbon profiles of the progenitor WD on the explosionphase and on synthetic observables in the Chandrasekhar-mass delayeddetonation model. Spectra and light curves are compared to observationsto judge the validity of the model.
Methods: The explosion phaseis simulated using the finite volume supernova code Leafs, which isextended to treat different compositions of the progenitor WD. Thesynthetic observables are computed with the Monte Carlo radiativetransfer code Artis. Results: Differences in binding energies ofcarbon and oxygen lead to a lower nuclear energy release for carbondepleted material; thus, the burning fronts that develop are weaker andthe total nuclear energy release is smaller. For otherwise identicalconditions, carbon depleted models produce less 56Ni.Comparing different models with similar 56Ni yields showslower kinetic energies in the ejecta for carbon depleted models, butonly small differences in velocity distributions and line velocities inspectra. The light curve width-luminosity relation (WLR) obtained formodels with differing carbon depletion is roughly perpendicular to theobserved WLR, hence the carbon mass fraction is probably only asecondary parameter in the family of SNe Ia.
Tables 3 and 4 are available in electronic form at http://www.aanda.org

Gamma-ray diagnostics of Type Ia supernovae. Predictions of observables from three-dimensional modeling

Context. Although the question of progenitor systems and detailed explosion mechanisms still remains a matter of discussion, it is commonly believed that Type Ia supernovae (SNe Ia) are production sites of large amounts of radioactive nuclei. Even though the gamma-ray emission due to radioactive decays is responsible for powering the light curves of SNe Ia, gamma rays themselves are of particular interest as a diagnostic tool because they directly lead to deeper insight into the nucleosynthesis and the kinematics of these explosion events. Aims: We study the evolution of gamma-ray line and continuum emission of SNe Ia with the objective of analyzing the relevance of observations in this energy range. We seek to investigate the chances for the success of future MeV missions regarding their capabilities for constraining the intrinsic properties and the physical processes of SNe Ia. Methods: Focusing on two of the most broadly discussed SN Ia progenitor scenarios - a delayed detonation in a Chandrasekhar-mass white dwarf (WD) and a violent merger of two WDs - we used three-dimensional explosion models and performed radiative transfer simulations to obtain synthetic gamma-ray spectra. Both chosen models produce the same mass of 56Ni and have similar optical properties that are in reasonable agreement with the recently observed supernova SN 2011fe. We examine the gamma-ray spectra with respect to their distinct features and draw connections to certain characteristics of the explosion models. Applying diagnostics, such as line and hardness ratios, the detection prospects for future gamma-ray missions with higher sensitivities in the MeV energy range are discussed. Results: In contrast to the optical regime, the gamma-ray emission of our two chosen models proves to be quite different. The almost direct connection of the emission of gamma rays to fundamental physical processes occurring in SNe Ia permits additional constraints concerning several explosion model properties that are not easily accessible within other wavelength ranges. Proposed future MeV missions such as GRIPS will resolve all spectral details only for nearby SNe Ia, but hardness ratio and light curve measurements still allow for a distinction of the two different models at 10 Mpc and 16 Mpc for an exposure time of 106 s. The possibility of detecting the strongest line features up to the Virgo distance will offer the opportunity to build up a first sample of SN Ia detections in the gamma-ray energy range and underlines the importance of future space observatories for MeV gamma rays.

Young Stellar Object VARiability (YSOVAR): Long Timescale Variations in the Mid-infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 and 4.5 μm) time series photometry of the Orion Nebula Cluster plus smaller footprints in 11 other star-forming cores (AFGL 490, NGC 1333, Mon R2, GGD 12-15, NGC 2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC 1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs roughly twice per day for ~40 days. We also define a "standard sample of members" consisting of all the IR-selected members and X-ray-selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data—the Stetson index, a χ2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of six to seven years by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data set; out of members and field objects combined, at most 0.02% may have transient IR excesses.

YSOVAR: Mid-infrared Variability in the Star-forming Region Lynds 1688

The emission from young stellar objects (YSOs) in the mid-infrared (mid-IR) is dominated by the inner rim of their circumstellar disks. We present IR data from the Young Stellar Object VARiability (YSOVAR) survey of ~800 objects in the direction of the Lynds 1688 (L1688) star-forming region over four visibility windows spanning 1.6 yr using the Spitzer Space Telescope in its warm mission phase. Among all light curves, 57 sources are cluster members identified based on their spectral energy distribution and X-ray emission. Almost all cluster members show significant variability. The amplitude of the variability is larger in more embedded YSOs. Ten out of 57 cluster members have periodic variations in the light curves with periods typically between three and seven days, but even for those sources, significant variability in addition to the periodic signal can be seen. No period is stable over 1.6 yr. Nonperiodic light curves often still show a preferred timescale of variability that is longer for more embedded sources. About half of all sources exhibit redder colors in a fainter state. This is compatible with time-variable absorption toward the YSO. The other half becomes bluer when fainter. These colors can only be explained with significant changes in the structure of the inner disk. No relation between mid-IR variability and stellar effective temperature or X-ray spectrum is found.

YSOVAR: Mid-infrared Variability of Young Stellar Objects and Their Disks in the Cluster IRAS 20050+2720

We present a time-variability study of young stellar objects (YSOs) in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 μm with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability (YSOVAR) project. We have collected light curves for 181 cluster members over 60 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2–6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color–magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability timescales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer timescales than the X-ray undetected members.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

YSOVAR: Mid-infrared Variability in NGC 1333

As part of the Young Stellar Object VARiability (YSOVAR) program, wemonitored NGC 1333 for ∼35 days at 3.6 and 4.5 μm using theSpitzer Space Telescope. We report here on the mid-infrared variabilityof the point sources in the ∼10‧ × ∼20‧ areacentered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves ineither channel, we find 78 variables over the YSOVAR campaign. Abouthalf of the members are variable. The variable fraction for the mostembedded spectral energy distributions (SEDs) (Class I, flat) is higherthan that for less embedded SEDs (Class II), which is in turn higherthan the star-like SEDs (Class III). A few objects have amplitudes(10–90th percentile brightness) in [3.6] or [4.5] > 0.2 mag; amore typical amplitude is 0.1–0.15 mag. The largest color changeis >0.2 mag. There are 24 periodic objects, with 40% of them beingflat SED class. This may mean that the periodic signal is primarily fromthe disk, not the photosphere, in those cases. We find 9 variableslikely to be “dippers,” where texture in the disk occultsthe central star, and 11 likely to be “bursters,” whereaccretion instabilities create brightness bursts. There are 39 objectsthat have significant trends in [3.6]–[4.5] color over thecampaign, about evenly divided between redder-when-fainter (consistentwith extinction variations) and bluer-when-fainter. About a third of the17 Class 0 and/or jet-driving sources from the literature are variableover the YSOVAR campaign, and a larger fraction (∼half) are variablebetween the YSOVAR campaign and the cryogenic-era Spitzer observations(6–7 years), perhaps because it takes time for the envelope torespond to changes in the central source. The NGC 1333 brown dwarfs donot stand out from the stellar light curves in any way except there is amuch larger fraction of periodic objects (∼60% of variable browndwarfs are periodic, compared to ∼30% of the variables overall).

Where are the mid-sized flares of ultracool M dwarfs?

We propose to observe the M8.5 dwarf SCR J1845-6357 with XMM-Newton EPIC for 60 ks. Very low-mass M dwarfs show a distinct drop in X-ray luminosity compared to slightly more massive M dwarfs. Surprisingly, this does not happen at the mass threshold where M dwarfs become fully convective (M4), but at significantly lower masses (M8). These very low mass stars seem to have a flaring behaviour different from earlier type stars: they display either occasional large flares or a very low-level "flickering" in their X-ray light curves, but not the canonical power-law flare-energy distribution observed for the Sun and other cool stars. Our aim is to collect a long-duration light curve for one of the most nearby ultracool dwarfs to quantify how its flare-energy distribution differs from earlier type stars.

M Dwarf Flare Continuum Variations on One-second Timescales: Calibrating and Modeling of ULTRACAM Flare Color Indices

We present a large data set of high-cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the near-ultraviolet (NUV) and optical continuum spectral evolution on timescales of ≈1 s. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new “color-color” relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T ≈ 104 K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as observed in moderate-sized and large flares alike. We also discuss the high time resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100. Based on observations obtained with the Apache Point Observatory 3.5 m telescope, which is owned and operated by the Astrophysical Research Consortium, based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofsica de Canarias, and observations, and based on observations made with the ESO Telescopes at the La Silla Paranal Observatory under programme ID 085.D-0501(A).