Search

We present EXplosive TRAnsient Spectral Simulator (extrass), a newly developed code aimed at generating 3D spectra for supernovae in the nebular phase by using modern multidimensional explosion models as input. It is well established that supernovae are asymmetric by nature, and that the morphology is encoded in the line profiles during the nebular phase, months after the explosion. In this work, we use extrass to study one such simulation of a $3.3\, \mathrm{ M}_\odot$ He-core explosion ($M_\text{ejecta}=1.3\, M_\odot$, $E_\text{kin}=1.05\times 10^{51}\,$erg) modelled with the Prometheus-HotB code and evolved to the homologous phase. Our code calculates the energy deposition from the radioactive decay of 56Ni → 56Co → 56Fe and uses this to determine the Non-Local-Thermodynamic-Equilibrium temperature, excitation, and ionization structure across the nebula. From the physical condition solutions we generate the emissivities to construct spectra depending on viewing angles. Our results show large variations in the line profiles with viewing angles, as diagnosed by the first three moments of the line profiles; shifts, widths, and skewness. We compare line profiles from different elements, and study the morphology of line-of-sight slices that determine the flux at each part of a line profile. We find that excitation conditions can sometimes make the momentum vector of the ejecta emitting in the excited states significantly different from that of the bulk of the ejecta of the respective element, thus giving blueshifted lines for bulk receding material, and vice versa. We compare the 3.3 M⊙ He-core model to observations of the Type Ib supernova SN 2007Y.

A binary neutron star merger produces a rapidly evolving transient known as a kilonova (KN), which peaks a few days after merger. Modelling of KNe has often been approached assuming local thermodynamic equilibrium (LTE) conditions in the ejecta. We present the first analysis of non-local thermodynamic equilibrium (NLTE) level populations, using the spectral synthesis code SUMO, and compare these to LTE values. We investigate the importance of the radiation field by conducting NLTE excitation calculations with and without radiative transfer. Level populations, in particular higher lying ones, start to show deviations from LTE a few days after merger. Excitation is lower in NLTE for the majority of ions and states, and this tends to give lower expansion opacities. While the difference is small for the first few days, it grows to factors 2-10 after this. Our results are important both for demonstrating validity of LTE expansion opacities for an initial phase (few days), while highlighting the need for NLTE modelling during later phases. Considering also NLTE ionisation, our results indicate that NLTE can give both higher or lower opacities, depending on composition and wavelength, sometimes by orders of magnitudes., 15 pages, 7 figures in main body, 31 pages, 24 figures with appendices. Submitted to MNRAS

We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features., Accepted for publication in A&A. 27 pages, 12 figures (v4: updated to match published version). The ejecta models and output files from the simulations are available at https://github.com/sn-rad-trans/data1

Superluminous supernovae radiate up to 100 times more energy than normal supernovae. The origin of this energy and the nature of their stellar progenitors are poorly understood. We identify neutral iron lines in the spectrum of one such transient, SN 2006gy, and show that they require a large mass of iron (>~0.3 Msun) expanding at 1500 km/s. We demonstrate that a model of a standard Type Ia supernova hitting a shell of circumstellar material produces a light curve and late-time iron-dominated spectrum that match SN 2006gy. In such a scenario, common envelope evolution of the progenitor system can synchronize envelope ejection and supernova explosion and may explain these bright transients., Author version of paper published in Science on Jan 24 2020

Supernovae provide environments with strong links to laboratory astrophysics. Diverse physical processes spanning from hot gas and semi-relativistic particles down to cold dusty clumps require extensive atomic data and understanding of processes across different physical regimes. The current status of modelling and analyzing supernova spectra is reviewed, with focus on recent results for diagnosing the production of oxygen and nickel.

We thank the anonymous referee for the very useful comments, which contributed to improve the manuscript. AF is partially supported by the PRIN-INAF 2017 with the project Towards the SKA and CTA era: discovery, localisation, and physics of transients sources (P.I. M. Giroletti). These observations made use of the LCO network. DAH, CP, DH, and JB are supported by NSF Grant AST1911225 and NASA Grant 80NSSC19k1639. TMB was funded by the CONICYT PFCHA/DOCTORADOBECAS CHILE/201772180113. MG is supported by the Polish NCN MAESTRO grant 2014/14/A/ST9/00121. TWC acknowledges the funding provided by the Alexander von Humboldt Foundation and the EU Funding under Marie Sklodowska-Curie grant agreement No 842471, and Thomas Kruhler for reducing X-Shooter spectrum. LG was funded by the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 839090. This work has been partially supported by the Spanish grant PGC2018-095317-B-C21 within the European Funds for Regional Development (FEDER). CPG acknowledges support from EU/FP7ERC grant no. [615929]. GL was supported by a research grant (19054) from VILLUM FONDEN. MN is supported by a Royal Astronomical Society Research Fellowship. RL is supported by a Marie Sklodowska-Curie Individual Fellowship within the Horizon 2020 European Union (EU) Framework Programme for Research and Innovation (H2020-MSCA-IF-2017-794467). GT acknowledges partial support by the National Science Foundation under Award No. AST-1909796. Research by SV is supported by NSF grants AST-1813176 and AST-2008108. Some of the observations reported here were obtained at the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution under program 2018A-UAO-G16 (PI Terreran). Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration under program NW440 (PI Fong). The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. W. M. Keck Observatory and MMT Observatory accesswas supported by Northwestern University and the Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA). Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programmes 199.D-0143, 0100.D-0751(B), 0101.D-0199(B), 099.A-9025(A), 0100.A-9099(A)099.A-9099 and 0100.A-9099. This work makes use of observations from the LCO network. Part of the funding for GROND (both hardware as well as personnel) was generously granted from the Leibniz-Prize to Prof. G. Hasinger (DFG grant HA 1850/28-1). The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg, and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen's University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration Grants No.s NNX08AR22G, NNX12AR65G, and NNX14AM74G, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory and the Gordon and Betty Moore foundation. TheATLAS surveys are funded through NASA grants NNX12AR55G. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This research made use of TARDIS, a community-developed software package for spectral synthesis in supernovae (Kerzendorf & Sim 2014). The development of TARDIS received support from the Google Summer of Code initiative and from ESA's Summer of Code in Space program. TARDIS makes extensive use of Astropy and PyNE., This article has been accepted for publication in MNRAS published by Oxford University Press on behalf of the Royal Astronomical Society, We present and discuss the optical spectrophotometric observations of the nearby (z = 0.087) Type I superluminous supernova (SLSN I) SN 2017gci, whose peak K-corrected absolute magnitude reaches M-g = -21.5 mag. Its photometric and spectroscopic evolution includes features of both slow- and of fast-evolving SLSN I, thus favoring a continuum distribution between the two SLSN-I subclasses. In particular, similarly to other SLSNe I, the multiband light curves (LCs) of SN 2017gci show two re-brightenings at about 103 and 142 d after the maximum light. Interestingly, this broadly agrees with a broad emission feature emerging around 6520 angstrom after similar to 51 d from the maximum light, which is followed by a sharp knee in the LC. If we interpret this feature as H alpha, this could support the fact that the bumps are the signature of late interactions of the ejecta with a (hydrogen-rich) circumstellar material. Then we fitted magnetar- and CSM-interaction-powered synthetic LCs on to the bolometric one of SN 2017gci. In the magnetar case, the fit suggests a polar magnetic field B-p similar or equal to 6 x 10(14) G, an initial period of the magnetar P-initial similar or equal to 2.8 ms, an ejecta mass M-ejecta similar or equal to 9M(circle dot) and an ejecta opacity kappa similar or equal to 0.08 cm(2) g(-1). A CSM-interaction scenario would imply a CSM mass similar or equal to 5 M-circle dot and an ejecta mass similar or equal to 12M(circle dot). Finally, the nebular spectrum of phase + 187 d was modeled, deriving a mass of similar or equal to 10 M-circle dot for the ejecta. Our models suggest that either a magnetar or CSM interaction might be the power sources for SN 2017gci and that its progenitor was a massive (40 M-circle dot) star., PRIN-INAF 2017, National Science Foundation (NSF) AST1911225 AST-1813176 AST-2008108, CONICYT PFCHA/DOCTORADOBECAS CHILE/2017 72180113, Polish NCN MAESTRO 2014/14/A/ST9/00121, Alexander von Humboldt Foundation, European Commission 842471, European Commission 839090, European Commission PGC2018-095317-B-C21, VILLUM FONDEN 19054, Royal Astronomical Society Research Fellowship, Marie Sklodowska-Curie Individual Fellowship within the Horizon 2020 European Union (EU) Framework Programme for Research and Innovation H2020-MSCA-IF-2017-794467, National Science Foundation (NSF) AST-1909796 AST-1238877, University of Arizona 2018A-UAO-G16, Smithsonian Institution 2018A-UAO-G16, National Aeronautics & Space Administration (NASA) NW440 NNX08AR22G NNX12AR65G NNX14AM74G, W.M. Keck Foundation, Northwestern University, Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA), European Organisation for Astronomical Research in the Southern Hemisphere 199.D-0143 0100.D-0751(B) 0101.D-0199(B) 099.A-9025(A) 0100.A-9099(A)099.A-9099 0100.A-9099, German Research Foundation (DFG), European Commission HA 1850/28-1, National Aeronautics & Space Administration (NASA) NNX12AR55G 80NSSC19k1639, Google Summer of Code initiative, ESA's Summer of Code in Space program, European Commission 615929

Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.--Full list of authors: Ackley, K.; Amati, L.; Barbieri, C.; Bauer, F. E.; Benetti, S.; Bernardini, M. G.; Bhirombhakdi, K.; Botticella, M. T.; Branchesi, M.; Brocato, E.; Bruun, S. H.; Bulla, M.; Campana, S.; Cappellaro, E.; Castro-Tirado, A. J.; Chambers, K. C.; Chaty, S.; Chen, T. -W.; Ciolfi, R.; Coleiro, A.; Copperwheat, C. M.; Covino, S.; Cutter, R.; D'Ammando, F.; D'Avanzo, P.; De Cesare, G.; D'Elia, V.; Della Valle, M.; Denneau, L.; De Pasquale, M.; Dhillon, V. S.; Dyer, M. J.; Elias-Rosa, N.; Evans, P. A.; Eyles-Ferris, R. A. J.; Fiore, A.; Fraser, M.; Fruchter, A. S.; Fynbo, J. P. U.; Galbany, L.; Gall, C.; Galloway, D. K.; Getman, F. I.; Ghirlanda, G.; Gillanders, J. H.; Gomboc, A.; Gompertz, B. P.; González-Fernández, C.; González-Gaitán, S.; Grado, A.; Greco, G.; Gromadzki, M.; Groot, P. J.; Gutiérrez, C. P.; Heikkilä, T.; Heintz, K. E.; Hjorth, J.; Hu, Y. -D.; Huber, M. E.; Inserra, C.; Izzo, L.; Japelj, J.; Jerkstrand, A.; Jin, Z. P.; Jonker, P. G.; Kankare, E.; Kann, D. A.; Kennedy, M.; Kim, S.; Klose, S.; Kool, E. C.; Kotak, R.; Kuncarayakti, H.; Lamb, G. P.; Leloudas, G.; Levan, A. J.; Longo, F.; Lowe, T. B.; Lyman, J. D.; Magnier, E.; Maguire, K.; Maiorano, E.; Mandel, I.; Mapelli, M.; Mattila, S.; McBrien, O. R.; Melandri, A.; Michałowski, M. J.; Milvang-Jensen, B.; Moran, S.; Nicastro, L.; Nicholl, M.; Nicuesa Guelbenzu, A.; Nuttal, L.; Oates, S. R.; O'Brien, P. T.; Onori, F.; Palazzi, E.; Patricelli, B.; Perego, A.; Torres, M. A. P.; Perley, D. A.; Pian, E.; Pignata, G.; Piranomonte, S.; Poshyachinda, S.; Possenti, A.; Pumo, M. L.; Quirola-Vásquez, J.; Ragosta, F.; Ramsay, G.; Rau, A.; Rest, A.; Reynolds, T. M.; Rosetti, S. S.; Rossi, A.; Rosswog, S.; Sabha, N. B.; Sagués Carracedo, A.; Salafia, O. S.; Salmon, L.; Salvaterra, R.; Savaglio, S.; Sbordone, L.; Schady, P.; Schipani, P.; Schultz, A. S. B.; Schweyer, T.; Smartt, S. J.; Smith, K. W.; Smith, M.; Sollerman, J.; Srivastav, S.; Stanway, E. R.; Starling, R. L. C.; Steeghs, D.; Stratta, G.; Stubbs, C. W.; Tanvir, N. R.; Testa, V.; Thrane, E.; Tonry, J. L.; Turatto, M.; Ulaczyk, K.; van der Horst, A. J.; Vergani, S. D.; Walton, N. A.; Watson, D.; Wiersema, K.; Wiik, K.; Wyrzykowski, Ł.; Yang, S.; Yi, S. -X.; Young, D. R., Context. Gravitational wave (GW) astronomy has rapidly reached maturity, becoming a fundamental observing window for modern astrophysics. The coalescences of a few tens of black hole (BH) binaries have been detected, while the number of events possibly including a neutron star (NS) is still limited to a few. On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. A preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. Aims. In this paper, we present our extensive search campaign aimed at uncovering the potential optical and near infrared electromagnetic counterpart of S190814bv. We found no convincing electromagnetic counterpart in our data. We therefore use our non-detection to place limits on the properties of the putative outflows that could have been produced by the binary during and after the merger. Methods. Thanks to the three-detector observation of S190814bv, and given the characteristics of the signal, the LIGO and Virgo Collaborations delivered a relatively narrow localisation in low latency - a 50% (90%) credible area of 5 deg2 (23 deg2) - despite the relatively large distance of 267 ± 52 Mpc. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical and near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS, and VINROUGE projects also carried out a search on this event. In this paper, we describe the combined observational campaign of these groups. Results. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN), which was possibly generated by this NS-BH merger, and for the strategy of future searches. The typical depth of our wide-field observations, which cover most of the projected sky localisation probability (up to 99.8%, depending on the night and filter considered), is r ∼ 22 (resp. K ∼ 21) in the optical (resp. near infrared). We reach deeper limits in a subset of our galaxy-targeted observations, which cover a total ∼50% of the galaxy-mass-weighted localisation probability. Altogether, our observations allow us to exclude a KN with large ejecta mass M 0.1 M- to a high (> 90%) confidence, and we can exclude much smaller masses in a sub-sample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Conclusions. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv, we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundred megaparsecs will be detected only by large facilities with both a high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event. © K. Ackley et al. 2020., Based on observations collected at the European Southern Observatory under ESO programmes 1102.D-0353(E), 1102.D0353(F), 1102.D-0353(Q), 1102.D-0353(G), 0103.D-0070(A), 0103.D-0070(B), 0103.D-0703(A), 0103.D-0722(A), 0103.A-9099(A), 198.D-2010(D) and 60.A9285(A). ATLAS is primarily funded through NEO NASA grants NN12AR55G, 80NSSC18K0284, and 80NSSC18K1575. The ATLAS science products have been made possible through the contributions of the University of Hawaii IfA, the Queen's University Belfast, the Space Telescope Science Institute, and the South African Astronomical Observatory. PanSTARRS is primarily funded through NEO NASA grants NASA Grants NNX08AR22G, NNX14AM74G. The PanSTARRS science products for LIGO-Virgo follow-up are made possible through the contributions of the University of Hawaii IfA and the Queen's University Belfast. The Gravitational-wave Optical Transient Observer (GOTO) project acknowledges the support of the Monash-Warwick Alliance; Warwick University; Monash University; She ffield University; Leicester University; Armagh Observatory & Planetarium; the National Astronomical Research Institute of Thailand (NARIT); University of Portsmouth; Turku University and the Instituto de Astrofisica de Canarias (IAC). Part of the funding for GROND was generously granted from the Leibniz-Prize to Prof. G. Hasinger (DFG grant HA 1850/28-1). The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council. The WHT and its override programme are operated on the island of La Palma by the Isaac Newton Group of Telescopes in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias; part of these data were taken under program (19A)N3. FEB thanks CONICYT Basal AFB-170002 and Chile's Ministry of Economy fund IC120009. MGB, PDA and AM acknowledge support from ASI grant I/004/11/3. MBr, EC, AP and SPi acknowledge support from MIUR (PRIN 2017 grant 20179ZF5KS). EB, EM and MT acknowledge funding from GRAWITA. SHB is indebted to the Danish National Research Foundation (DNRF132) for support. SCa acknowledges support from grant MAE0065741. EC acknowledges the support of the H2020 OPTICON programme 730890. TWC acknowledges the Humboldt Foundation and Marie Sklodowska-Curie grant 842471. MDP thanks Istanbul University for support. PAE acknowledges UKSA support. RAJEF is supported by an STFC studentship. MF is supported by a Royal Society -SFI University Research Fellowship. LG was funded by the EU H2020 programme under MSCA grant no. 839090. CG, JH and LI were supported by a research grant from VILLUM FONDEN (project 16599). CG and LI were supported by a research grant from VILLUM FONDEN (25501). GGh acknowledges the PRIN MIUR "Figaro" for financial support. AGo acknowledges financial support from the Slovenian Research Agency (grants P1-0031, I0-0033, and J1-8136). BPG, AJL and JDL acknowledge support from ERC grant 725246 (TEDE, PI Levan). SGG acknowledges support by FCT Fundacao para a Ciencia e Tecnologia and by Project PTDC/FIS-AST-31546. GGr acknowledges the ESCAPE H2020 project no. 824064. MG is supported by the Polish NCN MAESTRO grant 2014/14/A/ST9/00121. PJG acknowledges support from NOVA and from the South African NRF SARChI grant 111692. CPG and MS acknowledge support from EU/FP7-ERC grant no. 615929. KEH acknowledges support by a Project Grant from The Icelandic Research Fund. YDH acknowledges support from the China Scholarships Council. JJ acknowledges support from NOVA and NWO-FAPESP grant for instrumentation. AJ acknowledges funding from the European Research Council (ERC). ZPJ was supported by the Foundation for Distinguished Young Scholars of Jiangsu Province (no. BK20180050). PGJ acknowledges funding from the ERC under Consolidator Grant agreement no. 647208. DAK acknowledges Spanish research project RTI2018-098104-J-I00 (GRBPhot). SKl acknowledges support by DFG grant Kl 766/16-3. ECK acknowledges support from the GREAT research environment. GPL acknowledges support from STFC via grant ST/N000757/1. GL was supported by a research grant (19054) from VILLUM FONDEN. KM acknowledges support from the ERC (grant no. 758638). IM is partially supported by OzGrav (ARC project CE17010000). MMacknowledges support from ERC through ERC-2017-CoG no. 770017. MJM acknowledges the National Science Centre, Poland, grant 2018/30/E/ST9/00208. BMJ and DW are supported by Independent Research Fund Denmark grant DFF-7014-00017. MN is supported by a Royal Astronomical Society Research Fellowship. ANG acknowledges support by grant DFG Kl 766/16-3. PTOB acknowledges funding from STFC. SRO gratefully acknowledges the support of the Leverhulme Trust. FO acknowledges the support of the H2020 Hemera program, grant no. 730970. MAPT was supported by grants RYC-2015-17854 and AYA201783216-P. EP aknowledges financial support from INAF. GP is supported by the Millennium Science Initiative through grant IC120009. MLP is partially supported by a "Linea 2" project of the Catania University. JQV acknowledges support from CONICYT folio 21180886. TMR acknowledges the support of the Vilho, Yrjo and Kalle Vaisala Foundation. ARo acknowledges support from Premiale LBT 2013. SR is supported by VR grants 2016-03657_3 and the research environment grant GREAT, Dnr. 2016-06012, and the Swedish National Space board, Dnr. 107/16. OSS acknowledges the Italian Ministry of Research (MIUR) grant 1.05.06.13. LSa acknowledges the Irish Research Council Scholarship no. GOIPG/2017/1525. SJS acknowledges support from STFC Grant ST/P000312/1. ERS and DS acknowledge funding from UK STFC CG ST/P000495/1. RLCS acknowledges funding from STFC. DS acknowledges support from STFC via grant ST/T007184/1. SDV acknowledges the support of the CNES. LWsupported by Polish NCN DAINA 2017/27/L/ST9/03221. The Cosmic DAWN center is funded by the Danish National Research Foundation.

Low-luminosity type II supernovae (LL SNe~II) make up the low explosion energy end of core-collapse SNe, but their study and physical understanding remain limited. We present SN\,2016aqf, a LL SN~II with extensive spectral and photometric coverage. We measure a $V$-band peak magnitude of $-14.58$\,mag, a plateau duration of $\sim$100\,days, and an inferred $^{56}$Ni mass of $0.008 \pm 0.002$\,\msun. The peak bolometric luminosity, L$_{\rm bol} \approx 10^{41.4}$\,erg\,s$^{-1}$, and its spectral evolution is typical of other SNe in the class. Using our late-time spectra, we measure the [\ion{O}{i}] $\lambda\lambda6300, 6364$ lines, which we compare against SN II spectral synthesis models to constrain the progenitor zero-age main-sequence mass. We find this to be 12 $\pm$ 3\,\msun. Our extensive late-time spectral coverage of the [\ion{Fe}{ii}] $\lambda7155$ and [\ion{Ni}{ii}] $\lambda7378$ lines permits a measurement of the Ni/Fe abundance ratio, a parameter sensitive to the inner progenitor structure and explosion mechanism dynamics. We measure a constant abundance ratio evolution of $0.081^{+0.009}_{-0.010}$, and argue that the best epochs to measure the ratio are at $\sim$200 -- 300\,days after explosion. We place this measurement in the context of a large sample of SNe II and compare against various physical, light-curve and spectral parameters, in search of trends which might allow indirect ways of constraining this ratio. We do not find correlations predicted by theoretical models; however, this may be the result of the exact choice of parameters and explosion mechanism in the models, the simplicity of them and/or primordial contamination in the measured abundance ratio., Comment: Accepted for publication in MNRAS

Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mecha ...

We thank the anonymous referee for their many suggestions that improved this paper. We thank Miguel Perez-Torres for helpful discussions. MN is supported by a Royal Astronomical Society Research Fellowship. TW is funded in part by European Research Council grant 320360 and by European Commission grant 730980. PGJ and GC acknowledge support from European Research Council Consolidator Grant 647208. GL and PC are supported by a research grant (19054) from Villum Fonden. MG is supported by the Polish NCN MAESTRO Grant 2014/14/A/ST9/00121. NI is partially supported by Polish NCN DAINA Grant 2017/27/L/ST9/03221. IA is a CIFAR Azrieli Global Scholar in the Gravity and the Extreme Universe Program and acknowledges support from that program, from the Israel Science Foundation (grant numbers 2108/18 and 2752/19), from the United States-Israel Binational Science Foundation (BSF), and from the Israeli Council forHigher Education Alon Fellowship. JB, DH, and CP were supported by NASA Grant 80NSSC18K0577. TWCacknowledges the EU Funding underMarie Sklodowska-Curie grant agreement No. 842471. LG was funded by the European Union's Horizon 2020 Framework Programme under theMarie Sklodowska-Curie grant agreement No. 839090. Thiswork has been partially supported by the Spanish grant PGC2018-095317-B-C21 within the European Funds for Regional Development (FEDER). SGGacknowledges support by FCT under Projects CRISP PTDC/FIS-AST-31546 and UIDB/00099/2020. IM is a recipient of the Australian Research Council Future Fellowship FT190100574. TMBwas funded by the CONICYT PFCHA/DOCTORADOBECAS CHILE/2017-72180113. KDA acknowledges support provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51403.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. This work is based on data collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile, under ESO programmes 1103.D-0328 and 0104.B-0709 and as part of ePESSTO+(the advanced Public ESO Spectroscopic Survey for Transient Objects Survey), observations from the Las Cumbres Observatory network, the 6.5 meterMagellan Telescopes located at Las Campanas Observatory, Chile, the MMT Observatory, a joint facility of the University of Arizona and the Smithsonian Institution and the WilliamHerschel Telescope (programme W19B/P7). Swift data were supplied by the UK Swift Science Data Centre at the University of Leicester. The Liverpool Telescope andWilliam Herschel Telescope are operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council., At 66 Mpc, AT2019qiz is the closest optical tidal disruption event (TDE) to date, with a luminosity intermediate between the bulk of the population and the faint-and-fast event iPTF16fnl. Its proximity allowed a very early detection and triggering of multiwavelength and spectroscopic follow-up well before maximum light. The velocity dispersion of the host galaxy and fits to the TDE light curve indicate a black hole mass ≈106M , disrupting a star of ≈1M . By analysing our comprehensive UV, optical, and X-ray data, we show that the early optical emission is dominated by an outflow, with a luminosity evolution L ∝ t2, consistent with a photosphere expanding at constant velocity ( 2000 km s−1), and a line-forming region producing initially blueshifted H and He II profiles with v = 3000–10 000 km s−1. The fastest optical ejecta approach the velocity inferred from radio detections (modelled in a forthcoming companion paper from K. D. Alexander et al.), thus the same outflow may be responsible for both the fast optical rise and the radio emission – the first time this connection has been observed in a TDE. The light-curve rise begins 29 ± 2 d before maximum light, peaking when the photosphere reaches the radius where optical photons can escape. The photosphere then undergoes a sudden transition, first cooling at constant radius then contracting at constant temperature. At the same time, the blueshifts disappear from the spectrum and Bowen fluorescence lines (N III) become prominent, implying a source of far-UV photons, while the X-ray light curve peaks at ≈1041 erg s−1. Assuming that these X-rays are from prompt accretion, the size and mass of the outflow are consistent with the reprocessing layer needed to explain the large optical to X-ray ratio in this and other optical TDEs, possibly favouring accretion-powered over collision-powered outflow models., Royal Astronomical Society Research Fellowship, European Research Council (ERC) 320360 647208, European Commission European Commission Joint Research Centre 730980, Villum Fonden 19054, Polish NCN MAESTRO Grant 2014/14/A/ST9/00121, Polish NCN DAINA Grant 2017/27/L/ST9/03221, Israel Science Foundation 2108/18 2752/19, US-Israel Binational Science Foundation, Israeli Council forHigher Education Alon Fellowship, National Aeronautics & Space Administration (NASA) 80NSSC18K0577 NAS5-26555, European Union (EU) 842471 839090, European Union (EU) PGC2018-095317-B-C21, Portuguese Foundation for Science and Technology CRISP PTDC/FIS-AST-31546 UIDB/00099/2020, Australian Research Council FT190100574, CONICYT PFCHA/DOCTORADOBECAS CHILE 2017-72180113, NASA through the NASA Hubble Fellowship grant - Space Telescope Science Institute HST-HF2-51403.001-A, European Organisation for Astronomical Research in the Southern Hemisphere, Chile, under ESO programmes 1103.D-0328 0104.B-0709, WilliamHerschel Telescope W19B/P7, Science & Technology Facilities Council (STFC)

Spectral observations of the type-IIb supernova (SN) 2016gkg at 300-800 days are reported. The spectra show nebular characteristics, revealing emission from the progenitor star's metal-rich core and providing clues to the kinematics and physical conditions of the explosion. The nebular spectra are dominated by emission lines of [O i] λλ6300, 6364 and [Ca ii] λλ7292, 7324. Other notable, albeit weaker, emission lines include Mg I] λ4571, [Fe ii] λ7155, O I λ7774, Ca II triplet, and a broad, boxy feature at the location of Hα. Unlike in other stripped-envelope SNe, the [O i] doublet is clearly resolved due to the presence of strong narrow components. The doublet shows an unprecedented emission line profile consisting of at least three components for each [O i]λ6300, 6364 line: a broad component (width ∼2000 km s⁻¹), and a pair of narrow blue and red components (width ∼300 km s⁻¹) mirrored against the rest velocity. The narrow component appears also in other lines, and is conspicuous in [O i]. This indicates the presence of multiple distinct kinematic components of material at low and high velocities. The low-velocity components are likely to be produced by a dense, slow-moving emitting region near the center, while the broad components are emitted over a larger volume. These observations suggest an asymmetric explosion, supporting the idea of two-component ejecta that influence the resulting late-time spectra and light curves. SN 2016gkg thus presents striking evidence for significant asymmetry in a standard-energy SN explosion. The presence of material at low velocity, which is not predicted in 1D simulations, emphasizes the importance of multidimensional explosion modeling of SNe., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

During the first few hundred days after the explosion, core-collapse supernovae (SNe) emit down-scattered X-rays and gamma-rays originating from radioactive line emissions, primarily from the $^{56}$Ni $\rightarrow$ $^{56}$Co $\rightarrow$ $^{56}$Fe chain. We use SN models based on three-dimensional neutrino-driven explosion simulations of single stars and mergers to compute this emission and compare the predictions with observations of SN 1987A. A number of models are clearly excluded, showing that high-energy emission is a powerful way of discriminating between models. The best models are almost consistent with the observations, but differences that cannot be matched by a suitable choice of viewing angle are evident. Therefore, our self-consistent models suggest that neutrino-driven explosions are able to produce, in principle, sufficient mixing, although remaining discrepancies may require small changes to the progenitor structures. The soft X-ray cutoff is primarily determined by the metallicity of the progenitor envelope. The main effect of asymmetries is to vary the flux level by a factor of ${\sim}$3. For the more asymmetric models, the shapes of the light curves also change. In addition to the models of SN 1987A, we investigate two models of Type II-P SNe and one model of a stripped-envelope Type IIb SN. The Type II-P models have similar observables as the models of SN 1987A, but the stripped-envelope SN model is significantly more luminous and evolves faster. Finally, we make simple predictions for future observations of nearby SNe., 27 pages, 12 figures. Accepted for publication in ApJ. Updated to match the accepted version; added two paragraphs in Section 2 on differences between the models/the level of mixing and rephrased the NuSTAR predictions slightly

Context. The inclusion of molecular physics is an important piece that tends to be missing from the puzzle when modeling the spectra of supernovae (SNe). Molecules have both a direct impact on the spectra, particularly in the infrared, and an indirect one as a result of their influence on certain physical conditions, such as temperature. Aims. In this paper, we aim to investigate molecular formation and non-local thermodynamic equilibrium (NLTE) cooling, with a particular focus on CO, the most commonly detected molecule in supernovae. We also aim to determine the dependency of supernova chemistry on physical parameters and the relative sensitivity to rate uncertainties. Methods. We implemented a chemical kinetic description of the destruction and formation of molecules into the SN spectral synthesis code SUMO. In addition, selected molecules were coupled into the full NLTE level population framework and, thus, we incorporated molecular NLTE cooling into the temperature equation. We produced a test model of the CO formation in SN 1987A between 150 and 600 days and investigated the sensitivity of the resulting molecular masses to the input parameters. Results. We find that there is a close inter-dependency between the thermal evolution and the amount of CO formed, mainly through an important temperature-sensitive CO destruction process with O+. After a few hundred days, CO completely dominates the cooling of the oxygen-carbon zone of the supernova which, therefore, contributes little optical emission. The uncertainty of the calculated CO mass scales approximately linearly with the typical uncertainty factor for individual rates. We demonstrate how molecular masses can potentially be used to constrain various physical parameters of the supernova.

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z = 0.1086 implies a peak brightness of M = -21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ∼2000 km s wide Hα and Hβ emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ∼50-60 km s based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ∼2000 km s wide Hα emission persisting at high luminosity (∼3 × 10 erg s) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability. © 2018. The American Astronomical Society. All rights reserved., T.W.C. acknowledges Thomas. Kruhler for the X-Shooter data reduction, Lin Yan and Claes Fransson for providing comparison spectra, Jason Spyromilio for useful discussions, Chien-Hsiu. Lee and You-Hua. Chu for coordinating observational resources, and funding from the Alexander von Humboldt Foundation. M.F. acknowledges the support of a Royal Society-Science Foundation Ireland University Research Fellowship. P.S. acknowledges support through the Sofia Kovalevskaja Award (Alexander von Humboldt Foundation). A.V.F. is grateful for the support of the TABASGO Foundation, the Christopher R. Redlich fund, and the Miller Institute for Basic Research in Science (U.C. Berkeley). A.J.R. and I.R.S. are supported by the Australian Research Council through grants FT170100243 and FT160100028, respectively. F.T. and J.S. acknowledge support from the KAW Foundation. S.J.S. acknowledges funding from the European Research Council Grant agreement #291222 and STFC grant ST/P000312/1. M.G. is supported by Polish National Science Centre grant OPUS 2015/17/B/ST9/03167. K.M. acknowledges support from the UK STFC through an Ernest Rutherford Fellowship and from a Horizon 2020 ERC Starting Grant (#758638). L.G. was supported in part by US NSF grant AST-1311862. C.P.G. acknowledges support from EU/FP7-ERC grant #615929. Z.Y.L., C.C.N., and P.C.Y. are grateful for funding from MoST (Taiwan) under grants 105-2112-M-008-002-MY3, 104-2923-M-008-004-MY5, and 106-2112-M-008-007. A.P. and S.B. are partially supported by PRIN-INAF 2017 "Toward the SKA and CTA era: discovery, localization, and physics of transient sources" (P.I.: Giroletti). A.G.-Y. is supported by the EU via ERC grant No. 725161, the Quantum Universe I-Core program, the ISF, the BSF Transformative program, and a Kimmel award. Part of the funding for GROND was generously granted from the Leibniz Prize to Prof. G. Hasinger (DFG grant HA 1850/28-1). Some observations were made with the Nordic Optical Telescope using ALFOSC. This publication has made use of data collected at Lulin Observatory, partly supported by MoST grant 105-2112-M-008-024-MY3. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and NASA; the observatory was made possible by the generous financial support of the W. M. Keck Foundation.

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe., Accepted for publication by Astronomy & Astrophysics

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to, 35 pages, 10 figures. Accepted for publication in ApJ. Updated to match accepted version; slightly shortened and clarified the discussion

The detection of the binary neutron star merger GW170817 together with the observation of electromagnetic counterparts across the entire spectrum inaugurated a new era of multimessenger astronomy. In this study, we incorporate wavelength-dependent opacities and emissivities calculated from atomic-structure data enabling us to model both the measured light curves and spectra of the electromagnetic transient AT2017gfo. Best fits of the observational data are obtained by Gaussian Process Regression, which allows us to present posterior samples for the kilonova and source properties connected to GW170817. Incorporating constraints obtained from the gravitational wave signal measured by the LIGO-Virgo Scientific Collaboration, we present a $90{{\ \rm per\ cent}}$ upper bound on the mass ratio q ≲ 1.38 and a lower bound on the tidal deformability of $\tilde{\Lambda } \gtrsim 197$, which rules out sufficiently soft equations of state. Our analysis is a path-finder for more realistic kilonova models and shows how the combination of gravitational wave and electromagnetic measurements allow for stringent constraints on the source parameters and the supranuclear equation of state.

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

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Parameterizing the light curve shape through rise and decline timescales shows that the two are highly correlated. Magnetar-powered models can reproduce the correlation, with the diversity in rise and decline rates driven by the diffusion timescale. Circumstellar interaction models can exhibit a similar rise-decline relation, but only for a narrow range of densities, which may be problematic for these models. We find that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times broader than SNe Ibc, but that the intrinsic shapes are similar. There are a number of SLSNe with particularly broad light curves, possibly indicating two progenitor channels, but statistical tests do not cleanly separate two populations. The general spectral evolution is also presented. Velocities measured from Fe II are similar for SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by a factor 2-3. Assuming $\kappa=0.1\,$cm$^2\,$g$^{-1}$, we estimate a mean (median) SLSN ejecta mass of 10$\,$M$_\odot$ (6$\,$M$_\odot$), with a range of 3-30$\,$M$_\odot$. Doubling the assumed opacity brings the masses closer to normal SNe Ibc, but with a high-mass tail. The most probable mechanism for generating SLSNe seems to be the core-collapse of a very massive hydrogen-poor star, forming a millisecond magnetar., Comment: 28 pages, 22 figs, 4 tables -- Updated on 2016-01-13 to fix typo in Table 1

Super-luminous supernovae of type Ic have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. PTF12dam is one of the closest and best studied super-luminous explosions that has a broad and slowly fading lightcurve similar to SN 2007bi. Here we present new photometry and spectroscopy for PTF12dam from 200-500 days (rest-frame) after peak and a detailed analysis of the host galaxy (SDSS J142446.21+461348.6 at z = 0.107). Using deep templates and image subtraction we show that the full lightcurve can be fit with a magnetar model if escape of high-energy gamma rays is taken into account. The full bolometric lightcurve from -53 to +399 days (with respect to peak) cannot be fit satisfactorily with the pair-instability models. An alternative model of interaction with a dense CSM produces a good fit to the data although this requires a very large mass (~ 13 M_sun) of hydrogen free CSM. The host galaxy is a compact dwarf (physical size ~ 1.9 kpc) and with M_g = -19.33 +/- 0.10, it is the brightest nearby SLSN Ic host discovered so far. The host is a low mass system (2.8 x 10^8 M_sun) with a star-formation rate (5.0 M_sun/year), which implies a very high specific star-formation rate (17.9 Gyr^-1). The remarkably strong nebular lines provide detections of the [O III] \lambda 4363 and [O II] \lambda\lambda 7320,7330 auroral lines and an accurate oxygen abundance of 12 + log(O/H) = 8.05 +/- 0.09. We show here that they are at the extreme end of the metallicity distribution of dwarf galaxies and propose that low metallicity is a requirement to produce these rare and peculiar supernovae., Comment: 20 pages, 12 figures, 8 tables, accepted for publication to MNRAS

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric-phase observational coverage. Here, we present Very Large Telescope and PESSTO observations of this SN in the nebular phase. We model the nebular [O I] 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 Msun. SN 2012ec is unique among hydrogen-rich SNe in showing a distinct and unblended line of stable nickel [Ni II] 7378. This line is produced by 58Ni, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] 7378 and [Fe II] 7155 (the progenitor of which is 56Ni) to derive a Ni/Fe production ratio of 0.20pm0.07 (by mass), which is a factor 3.4pm1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 micron line. This is the third reported case of a core-collapse supernova producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models., Published version

We present photometric and spectroscopic observations of the interacting transient SN 2009ip taken during the 2013 and 2014 observing seasons. We characterise the photometric evolution as a steady and smooth decline in all bands, with a decline rate that is slower than expected for a solely $^{56}$Co-powered supernova at late phases. No further outbursts or eruptions were seen over a two year period from 2012 December until 2014 December. SN 2009ip remains brighter than its historic minimum from pre-discovery images. Spectroscopically, SN 2009ip continues to be dominated by strong, narrow ($\lesssim$2000 km~s$^{-1}$) emission lines of H, He, Ca, and Fe. While we make tenuous detections of [Fe~{\sc ii}] $\lambda$7155 and [O~{\sc i}] $\lambda\lambda$6300,6364 lines at the end of 2013 June and the start of 2013 October respectively, we see no strong broad nebular emission lines that could point to a core-collapse origin. In general, the lines appear relatively symmetric, with the exception of our final spectrum in 2014 May, when we observe the appearance of a redshifted shoulder of emission at +550 km~s$^{-1}$. The lines are not blue-shifted, and we see no significant near- or mid-infrared excess. From the spectroscopic and photometric evolution of SN 2009ip until 820 days after the start of the 2012a event, we still see no conclusive evidence for core-collapse, although whether any such signs could be masked by ongoing interaction is unclear., Comment: Submitted to MNRAS

While interaction with circumstellar material is known to play an important role in Type IIn supernovae (SNe), analyses of the more common SNe IIP and IIL have not traditionally included interaction as a significant power source. However, recent campaigns to observe SNe within days of explosion have revealed narrow emission lines of high-ionization species in the earliest spectra of luminous SNe II of all subclasses. These "flash ionization" features indicate the presence of a confined shell of material around the progenitor star. Here we present the first low-luminosity (LL) SN to show flash ionization features, SN 2016bkv. This SN peaked at $M_V = -16$ mag and has H{\alpha} expansion velocities under 1350 km/s around maximum light, placing it at the faint/slow end of the distribution of SNe IIP (similar to SN 2005cs). The light-curve shape of SN 2016bkv is also extreme among SNe IIP. A very strong initial peak could indicate additional luminosity from circumstellar interaction. A very small fall from the plateau to the nickel tail indicates unusually large production of radioactive nickel compared to other LL SNe IIP. A comparison between nebular spectra of SN 2016bkv and models raises the possibility that SN 2016bkv is an electron-capture supernova., Comment: Updated to match published version in ApJ. Photometry table and SYN++ input file available at right

We report on the first clear detection of a kilonova; ejected radioactive debris from a neutron star merger. The transient AT2017gfo was detected 12h after GW170817 and associated short GRB170817A. It reached a peak UVOIR luminosity of $10^{42}$ erg/s on a time-scale of $\sim$1d. Modelling of the bolometric light curve suggests ejection of about 0.04 $M_\odot$ of ejecta at 0.2c. The ejecta has low opacity, indicating little lanthanide content. Spectroscopic modelling suggests identification of r-process elements Cs and Te. The combination of 2s delay between the GW and the GRB, a large ejecta mass, and low opacity all suggest dominance by a wind component and survival of a hypermassive neutron star for $t\sim 2$s.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa, Nature, in press, DOI 10.1038/nature24303. Data files will be made available at http://www.pessto.org

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

We present observations of the interacting transient SN 2009ip, from the start of the outburst in October 2012 until the end of the 2012 observing season. The transient reached a peak of $M_V$=-17.7 mag before fading rapidly, with a total integrated luminosity of 1.9$\times10^{49}$ erg over the period of August-December 2012. The optical and near infrared spectra are dominated by narrow emission lines, signaling a dense circumstellar environment, together with multiple components of broad emission and absorption in H and He at velocities between 0.5-1.2$\times10^4$ km s$^{-1}$\. We see no evidence for nucleosynthesized material in SN 2009ip, even in late-time pseudo-nebular spectra. We set a limit of $, Comment: 28 pages, submitted to MNRAS. Abstract abridged for arXiv

Nebular-phase observations and spectral models of Type Ic superluminous supernovae are presented. LSQ14an and SN 2015bn both display late-time spectra similar to galaxy-subtracted spectra of SN 2007bi, and the class shows strong similarity with broad-lined Type Ic SNe such as SN 1998bw. Near-infrared observations of SN 2015bn show a strong Ca II triplet, O I 9263, O I 1.13 um and Mg I 1.50 um, but no distinct He, Si, or S emission. The high Ca II NIR/[Ca II] 7291,7323 ratio of ~2 indicates a high electron density of n_e >~ 10^8 cm^{-3}. Spectral models of oxygen-zone emission are investigated to put constraints on the emitting region. Models require M(O-zone) >~ 10 Msun to produce enough [O I] 6300,6364 luminosity, irrespective of the powering situation and the density. The high oxygen-zone mass, supported by high estimated magnesium masses, points to explosions of massive CO cores, requiring M_ZAMS >~ 40 Msun. Collisions of pair-instability pulsations do not provide enough mass to account for the emission. [O II] and [O III] lines emerge naturally in many models, which strengthens the identification of broad [O II] 7320,7330, [O III] 4363, and [O III] 4959,5007 in some spectra. A small filling factor f, Comment: Published version

When supernovae enter the nebular phase after a few months, they reveal spectral fingerprints of their deep interiors, glowing by radioactivity produced in the explosion. We are given a unique opportunity to see what an exploded star looks like inside. The line profiles and luminosities encode information about physical conditions, explosive and hydrostatic nucleosynthesis, and ejecta morphology, which link to the progenitor properties and the explosion mechanism. Here, the fundamental properties of spectral formation of supernovae in the nebular phase are reviewed. The formalism between ejecta morphology and line profile shapes is derived, including effects of scattering and absorption. Line luminosity expressions are derived in various physical limits, with examples of applications from the literature. The physical processes at work in the supernova ejecta, including gamma-ray deposition, non-thermal electron degradation, ionization and excitation, and radiative transfer are described and linked to the computation and application of advanced spectral models. Some of the results derived so far from nebular-phase supernova analysis are discussed.

Results of calculations for inelastic e+Mg effective collision strengths for the lowest 25 physical states of Mg I (up to 3s6p 1P), and thus 300 transitions, from the convergent close-coupling (CCC) and the B-spline R-matrix (BSR) methods are presented. At temperatures of interest, ~5000 K, the results of the two calculations differ on average by only 4%, with a scatter of 27%. As the methods are independent, this suggests that the calculations provide datasets for e+Mg collisions accurate to this level. Comparison with the commonly used dataset compiled by Mauas et al. (1988), covering 25 transitions among 12 states, suggests the Mauas et al. data are on average ~57% too low, and with a very large scatter of a factor of ~6.5. In particular the collision strength for the transition corresponding to the Mg I intercombination line at 457 nm is significantly underestimated by Mauas et al., which has consequences for models that employ this dataset. In giant stars the new data leads to a stronger line compared to previous non-LTE calculations, and thus a reduction in the non-LTE abundance correction by ~0.1 dex (~25%). A non-LTE calculation in a supernova ejecta model shows this line becomes significantly stronger, by a factor of around two, alleviating the discrepancy where the 457 nm line in typical models with Mg/O ratios close to solar tended to be too weak compared to observations., Comment: Accepted for publication in Astronomy and Astrophysics

A small group of the newly discovered superluminous supernovae show broad and slowly evolving light curves. Here we present extensive observational data for the slow-evolving superluminous supernova LSQ14an, which brings this group of transients to four in total in the low redshift Universe (z$, Comment: Paper accepted by MNRAS on 31/03/2017

We present the one-year long observing campaign of SN 2012A which exploded in the nearby (9.8 Mpc) irregular galaxy NGC 3239. The photometric evolution is that of a normal type IIP supernova. The absolute maximum magnitude, with MB = -16.23 +- 0.16 mag. SN2012A reached a peak luminosity of about 2X10**42 erg/s, which is brighter than those of other SNe with a similar 56Ni mass. The latter was estimated from the luminosity in the exponential tail of the light curve and found to be M(56Ni) = 0.011 +-0.004 Msun. The spectral evolution of SN 2012A is also typical of SN IIP, from the early spectra dominated by a blue continuum and very broad (~10**4 km/s) Balmer lines, to the late-photospheric spectra characterized by prominent P-Cygni features of metal lines (Fe II, Sc II, Ba II, Ti II, Ca II, Na ID). The photospheric velocity is moderately low, ~3X10**3 km/s at 50 days, for the low optical depth metal lines. The nebular spectrum obtained 394 days after the shock breakout shows the typical features of SNe IIP and the strength of the [O I] doublet suggests a progenitor of intermediate mass, similar to SN 2004et (~15 Msun). A candidate progenitor for SN 2012A has been identified in deep, pre-explosion K'-band Gemini North (NIRI) images, and found to be consistent with a star with a bolometric magnitude -7.08+-0.36 (log L/Lsun = 4.73 +- 0.14$ dex). The magnitude of the recovered progenitor in archival images points toward a moderate-mass 10.5 (-2/+4.5) Msun star as the precursor of SN 2012A. The explosion parameters and progenitor mass were also estimated by means of a hydrodynamical model, fitting the bolometric light curve, the velocity and the temperature evolution. We found a best fit for a kinetic energy of 0.48 foe, an initial radius of 1.8X10**13 cm and ejecta mass of 12.5 Msun., 25 pages, 17 figures. Revised version (minor revision) submitted to MNRAS 2013

We present and analyse an extensive dataset of the superluminous supernova (SLSN) LSQ14mo (z = 0.256), consisting of a multi-colour lightcurve from -30 d to +70 d in the rest-frame and a series of 6 spectra from PESSTO covering -7 d to +50 d. This is among the densest spectroscopic coverage, and best-constrained rising lightcurve, for a fast-declining hydrogen-poor SLSN. The bolometric lightcurve can be reproduced with a millisecond magnetar model with ~ 4 M_sol ejecta mass, and the temperature and velocity evolution is also suggestive of a magnetar as the power source. Spectral modelling indicates that the SN ejected ~ 6 M_sol of CO-rich material with a kinetic energy of ~ 7 x 10^51 erg, and suggests a partially thermalised additional source of luminosity between -2 d and +22 d. This may be due to interaction with a shell of material originating from pre-explosion mass loss. We further present a detailed analysis of the host galaxy system of LSQ14mo. PESSTO and GROND imaging show three spatially resolved bright regions, and we used the VLT and FORS2 to obtain a deep (five-hour exposure) spectra of the SN position and the three star-forming regions, which are at a similar redshift. The FORS spectrum at +300 days shows no trace of SN emission lines and we place limits on the strength of [O I] from comparisons with other Ic SNe. The deep spectra provides a unique chance to investigate spatial variations in the host star-formation activity and metallicity. The specific star-formation rate is similar in all three components, as is the presence of a young stellar population. However, the position of LSQ14mo exhibits a lower metallicity, with 12 + log(O/H) = 8.2 in both the R23 and N2 scales (corresponding to ~ 0.3 Z_sol). We propose that the three bright regions in the host system are interacting, which thus triggers star-formation and forms young stellar populations., Accepted for publication in A&A

We investigate two stripped-envelope supernovae (SNe) discovered in the nearby galaxy NGC 5806 by the (i)PTF. These SNe, designated PTF12os/SN 2012P and iPTF13bvn, exploded at a similar distance from the host-galaxy center. We classify PTF12os as a Type IIb SN based on our spectral sequence; iPTF13bvn has previously been classified as Type Ib having a likely progenitor with zero age main sequence (ZAMS) mass below ~17 solar masses. Our main objective is to constrain the explosion parameters of iPTF12os and iPTF13bvn, and to put constraints on the SN progenitors. We present comprehensive datasets on the SNe, and introduce a new reference-subtraction pipeline (FPipe) currently in use by the iPTF. We perform a detailed study of the light curves (LCs) and spectral evolution of the SNe. The bolometric LCs are modeled using the hydrodynamical code HYDE. We use nebular models and late-time spectra to constrain the ZAMS mass of the progenitors. We perform image registration of ground-based images of PTF12os to archival HST images of NGC 5806 to identify a potential progenitor candidate. Our nebular spectra of iPTF13bvn indicate a low ZAMS mass of ~12 solar masses for the progenitor. The late-time spectra of PTF12os are consistent with a ZAMS mass of ~15 solar masses. We successfully identify a progenitor candidate to PTF12os using archival HST images. This source is consistent with being a cluster of massive stars. Our hydrodynamical modeling suggests that the progenitor of PTF12os had a compact He core with a mass of 3.25 solar masses, and that 0.063 solar masses of strongly mixed 56Ni was synthesized. Spectral comparisons to the Type IIb SN 2011dh indicate that the progenitor of PTF12os was surrounded by a hydrogen envelope with a mass lower than 0.02 solar masses. We also find tentative evidence that the progenitor of iPTF13bvn could have been surrounded by a small amount of hydrogen., 34 pages, 18 figures, accepted by A&A

We present nebular-phase imaging and spectroscopy for the hydrogen-poor superluminous supernova SN 2015bn, at redshift z=0.1136, spanning +250-400 d after maximum light. The light curve exhibits a steepening in the decline rate from 1.4 mag/(100 d) to 1.7 mag/(100 d), suggestive of a significant decrease in the opacity. This change is accompanied by a transition from a blue continuum superposed with photospheric absorption lines to a nebular spectrum dominated by emission lines of oxygen, calcium and magnesium. There are no obvious signatures of circumstellar interaction or large nickel mass. We show that the spectrum at +400 d is virtually identical to a number of energetic Type Ic supernovae such as SN 1997dq, SN 2012au, and SN 1998bw, indicating similar core conditions and strengthening the link between `hypernovae'/long gamma-ray bursts and superluminous supernovae. A single explosion mechanism may unify these events that span absolute magnitudes of -22 < M_B < -17. Both the light curve and spectrum of SN 2015bn are consistent with an engine-driven explosion ejecting 7-30 M$_\odot$ of oxygen-dominated ejecta (for reasonable choices in temperature and opacity). A strong and relatively narrow O II $��$7774 line, seen in a number of these energetic events but not in normal supernovae, may point to an inner shell that is the signature of a central engine., Updated to match published version

High-quality collections of Type II supernova (SN) light curves are scarce because they evolve for hundreds of days, making follow-up observations time consuming and often extending over multiple observing seasons. In light of these difficulties, the diversity of SNe II is not fully understood. Here we present ultraviolet and optical photometry of 12 SNe II monitored by the Las Cumbres Observatory Global Telescope Network (LCOGT) during 2013-2014, and compare them with previously studied SNe having well-sampled light curves. We explore SN II diversity by searching for correlations between the slope of the linear light-curve decay after maximum light (historically used to divide SNe II into IIL and IIP) and other measured physical properties. While SNe IIL are found to be on average more luminous than SNe IIP, SNe IIL do not appear to synthesize more 56Ni than SNe IIP. Finally, optical nebular spectra obtained for several SNe in our sample are found to be consistent with models of red supergiant progenitors in the 12-16 Msun range. Consequently, SNe IIL appear not to account for the deficit of massive red supergiants as SN II progenitors., 27 pages, 28 figures, submitted to MNRAS, data in supplementary table, accepted in MNRAS (Accepted 2016 April 12)

We present a study of the nebular phase spectra of a sample of Type II-Plateau supernovae with identified progenitors or restrictive limits. The evolution of line fluxes, shapes, and velocities are compared within the sample, and interpreted by the use of a spectral synthesis code. The small diversity within the dataset can be explained by strong mixing occurring during the explosion, and by recognising that most lines have significant contributions from primordial metals in the H envelope, which dominates the total ejecta mass in these type of objects. In particular, when using the [O I] 6300, 6364 Angstrom doublet for estimating the core mass of the star, care has to be taken to account for emission from primordial O in the envelope. Finally, a correlation between the H-alpha line width and the mass of 56Ni is presented, suggesting that higher energy explosions are associated with higher 56Ni production.

The interpretation of supernova (SN) spectra is essential for deriving SN ejecta properties such as density and composition, which in turn can tell us about their progenitors and the explosion mechanism. A very large number of atomic processes are important for spectrum formation. Several tools for calculating SN spectra exist, but they mainly focus on the very early or late epochs. The intermediate phase, which requires a NLTE treatment of radiation transport has rarely been studied. In this paper we present a new SN radiation transport code, NERO, which can look at those epochs. All the atomic processes are treated in full NLTE, under a steady-state assumption. This is a valid approach between roughly 50 and 500 days after the explosion depending on SN type. This covers the post-maximum photospheric and the early and the intermediate nebular phase. As a test, we compare NERO to the radiation transport code of Jerkstrand et al. (2011) and to the nebular code of Mazzali et al. (2001). All three codes have been developed independently and a comparison provides a valuable opportunity to investigate their reliability. Currently, NERO is one-dimensional and can be used for predicting spectra of synthetic explosion models or for deriving SN properties by spectral modelling. To demonstrate this, we study the spectra of the 'normal' SN Ia 2005cf between 50 and 350 days after the explosion and identify most of the common SN Ia line features at post maximum epochs.

We present optical and near-infrared photometric and spectroscopic observations of SN 2013ej, in galaxy M74, from 1 to 450 days after the explosion. SN 2013ej is a hydrogen-rich supernova, classified as a Type IIL due to its relatively fast decline following the initial peak. It has a relatively high peak luminosity (absolute magnitude M$_\rm{V}$ = -17.6) but a small $^{56}$Ni production of ~0.023 M$_\odot$. Its photospheric evolution is similar to other Type II SNe, with shallow absorption in the H{\alpha} profile typical for a Type IIL. During transition to the radioactive decay tail at ~100 days, we find the SN to grow bluer in B - V colour, in contrast to some other Type II supernovae. At late times, the bolometric light curve declined faster than expected from $^{56}$Co decay and we observed unusually broad and asymmetric nebular emission lines. Based on comparison of nebular emission lines most sensitive to the progenitor core mass, we find our observations are best matched to synthesized spectral models with a M$_\rm{ZAMS}$ = 12 - 15 M$_\odot$ progenitor. The derived mass range is similar to but not higher than the mass estimated for Type IIP progenitors. This is against the idea that Type IIL are from more massive stars. Observations are consistent with the SN having a progenitor with a relatively low-mass envelope., Comment: 23 pages, 18 figures, submitted to MNRAS

We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches., Submitted to ApJL on 6th June 2016. Accepted on 8th August 2016. This is the accepted version after referee report and minor revisions

Both CO and SiO have been observed at early and late phases in SN 1987A. H_2 was predicted to form at roughly the same time as these molecules, but was not detected at early epochs. Here we report the detection of NIR lines from H_2 at 2.12 mu and 2.40 mu in VLT/SINFONI spectra obtained between days 6489 and 10,120. The emission is concentrated to the core of the supernova in contrast to H-alpha and approximately coincides with the [Si I]/[Fe II] emission detected previously in the ejecta. Different excitation mechanisms and power sources of the emission are discussed. From the nearly constant H_2 luminosities we favour excitation resulting from the 44Ti decay., Comment: 5 pages 3 figures, accepted for ApJ Letters

If very massive stars (M >~ 100 Msun) can form and avoid too strong mass loss during their evolution, they are predicted to explode as pair-instability supernovae (PISNe). One critical test for candidate events is whether their nucleosynthesis yields and internal ejecta structure, being revealed through nebular-phase spectra at t >~ 1 yr, match those of model predictions. Here we compute theoretical spectra based on model PISN ejecta at 1-3 years post-explosion to allow quantitative comparison with observations. The high column densities of PISNe lead to complete gamma-ray trapping for t >~ 2 years which, combined with fulfilled conditions of steady state, leads to bolometric supernova luminosities matching the 56Co decay. Most of the gamma-rays are absorbed by the deep-lying iron and silicon/sulphur layers. The ionization balance shows a predominantly neutral gas state, which leads to emission lines of Fe I, Si I, and S I. For low-mass PISNe the metal core expands slowly enough to produce a forest of distinct lines, whereas high-mass PISNe expand faster and produce more featureless spectra. Line blocking is complete below ~5000 A for several years, and the model spectra are red. The strongest line is typically [Ca II] 7291,7323, one of few lines from ionized species. We compare our models with proposed PISN candidates SN 2007bi and PTF12dam, finding discrepancies for several key observables and thus no support for a PISN interpretation. We discuss distinct spectral features predicted by the models, and the possibility of detecting pair-instability explosions among non-superluminous supernovae., 25 pages, accepted for publication in MNRAS

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P., Comment: 19 pages, 10 figures, accepted by A&A

We present our findings on a supernova (SN) impostor, SNHunt248, based on optical and near-IR data spanning $\sim$15 yrs before discovery, to $\sim$1 yr post-discovery. The light curve displays three distinct peaks, the brightest of which is at $M_{R} \sim -15.0$ mag. The post-discovery evolution is consistent with the ejecta from the outburst interacting with two distinct regions of circumstellar material. The 0.5 - 2.2 $\mu$m spectral energy distribution at -740 d is well-matched by a single 6700 K blackbody with $\log(L/L_\odot) \sim 6.1$. This temperature and luminosity support previous suggestions of a yellow hypergiant progenitor; however, we find it to be brighter than the brightest and most massive Galactic late-F to early-G spectral type hypergiants. Overall the historical light curve displays variability of up to $\sim \pm1$ mag. At current epochs ($\sim$1 yr post-outburst), the absolute magnitude ($M_{R} \sim -9$ mag) is just below the faintest observed historical absolute magnitude $\sim$10 yrs before discovery., Comment: 7 pages, 3 figures, accepted for A&A Letters

We present data for LSQ14bdq, a hydrogen-poor super-luminous supernova (SLSN) discovered by the La Silla QUEST survey and classified by the Public ESO Spectroscopic Survey of Transient Objects. The spectrum and light curve are very similar to slow-declining SLSNe such as PTF12dam. However, detections within ~1 day after explosion show a bright and relatively fast initial peak, lasting for ~15 days, prior to the usual slow rise to maximum light. The broader, main peak can be fit with either central engine or circumstellar interaction models. We discuss the implications of the precursor peak in the context of these models. It is too bright and narrow to be explained as a normal 56Ni-powered SN, and we suggest that interaction models may struggle to fit the two peaks simultaneously. We propose that the initial peak may arise from the post-shock cooling of extended stellar material, and reheating by a central engine drives the second peak. In this picture, we show that an explosion energy of ~2x10^{52} erg and a progenitor radius of a few hundred solar radii would be required to power the early emission. The competing engine models involve rapidly spinning magnetars (neutron stars) or fall-back accretion onto a central black hole. The prompt energy required may favour the black hole scenario. The bright initial peak effectively rules out a compact Wolf-Rayet star as a progenitor, since the inferred energies and ejected masses become unphysical., Accepted in ApJ Letters

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies., ApJ accepted, 8 pages