Search

The ongoing discrepancy in the Hubble constant ($H_0$) estimates obtained through local distance ladder methods and early universe observations poses a significant challenge to the $\Lambda$CDM model, suggesting potential new physics. Type II supernovae (SNe II) offer a promising technique for determining $H_0$ in the local universe independently of the traditional distance ladder approach, opening up a complimentary path for testing this discrepancy. We aim to provide the first $H_0$ estimate using the tailored expanding photosphere method (EPM) applied to SNe II, made possible by recent advancements in spectral modelling that enhance its precision and efficiency. Our tailored EPM measurement utilizes a spectral emulator to interpolate between radiative transfer models calculated with TARDIS, allowing us to fit supernova spectra efficiently and derive self-consistent values for luminosity-related parameters. We apply the method on public data for ten SNe II at redshifts between 0.01 and 0.04. Our analysis demonstrates that the tailored EPM allows for $H_0$ measurements with precision comparable to the most competitive established techniques, even when applied to literature data not designed for cosmological applications. We find an independent $H_0$ value of $74.9\pm1.9$ (stat) km/s/Mpc, which is consistent with most current local measurements. Considering dominant sources of systematic effects, we conclude that our systematic uncertainty is comparable to or less than the current statistical uncertainty. This proof-of-principle study highlights the potential of the tailored EPM as a robust and precise tool for investigating the Hubble tension independently of the local distance ladder. Observations of SNe II tailored to $H_0$ estimation can make this an even more powerful tool by improving the precision and by allowing us to better understand and control systematic uncertainties., Comment: 40 pages, 57 figures, 4 tables; submitted to A&A

Neutron star (NS) mergers are known to produce heavy elements through rapid neutron capture (r-process) nucleosynthesis. Actinides are expected to be created solely by the r-process in the most neutron rich environments. Confirming if NS mergers provide the requisite conditions for actinide creation is therefore central to determining their origin in the Universe. Actinide signatures in kilonova (KN) spectra may yield an answer, provided adequate models are available in order to interpret observational data. In this study, we investigate actinide signatures in neutron rich merger ejecta. We use three ejecta models with different compositions and radioactive power, generated by nucleosynthesis calculations using the same initial electron fraction ($Y_e = 0.15$) but with different nuclear physics inputs and thermodynamic expansion history. These are evolved from 10 - 100 days after merger using the SUMO non-local thermodynamic equilibrium (NLTE) radiative transfer code. We highlight how uncertainties in nuclear properties, as well as choices in thermodynamic trajectory, may yield entirely different outputs for equal values of $Y_e$. We consider an actinide-free model and two actinide-rich models, and find that the emergent spectra and lightcurve evolution are significantly different depending on the amount of actinides present, and the overall decay properties of the models. We also present potential key actinide spectral signatures, of which doubly ionized $_{89}$Ac and $_{90}$Th may be particularly interesting as spectral indicators of actinide presence in KN ejecta., Comment: 21 pages with appendices, 11 figures, published in MNRAS

Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent)., Comment: Accepted by MNRAS

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $\alpha$ and $\beta^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned treatments, a line-by-line treatment enables us include fluorescence effects and associate spectral features with the emitting and absorbing lines of individual elements. We find variations in the synthetic observables with both the polar and azimuthal viewing angles. The spectra exhibit blended features with strong interactions by Ce III, Sr II, Y II, and Zr II that vary with time and viewing direction. We demonstrate the importance of wavelength-calibration of atomic data using a model with calibrated Sr, Y, and Zr data, and find major differences in the resulting spectra, including a better agreement with AT2017gfo. The synthetic spectra for near-polar inclination show a feature at around 8000 A, similar to AT2017gfo. However, they evolve on a more rapid timescale, likely due to the low ejecta mass (0.005 M$_\odot$) as we take into account only the early ejecta. The comparatively featureless spectra for equatorial observers gives a tentative prediction that future observations of edge-on kilonovae will appear substantially different from AT2017gfo. We also show that 1D models obtained by spherically averaging the 3D ejecta lead to dramatically different direction-integrated luminosities and spectra compared to full 3D calculations., Comment: 12 pages, 5 figures. Accepted by ApJL

Context: Type II supernovae provide a direct way to estimate distances through the expanding photosphere method, which is independent of the cosmic distance ladder. A recently introduced Gaussian process-based method allows for a fast and precise modelling of spectral time series, which puts accurate and computationally cheap Type II-based absolute distance determinations within reach. Aims: The goal of the paper is to assess the internal consistency of this new modelling technique coupled with the distance estimation empirically, using the spectral time series of supernova siblings, i.e. supernovae that exploded in the same host galaxy. Methods: We use a recently developed spectral emulator code, which is trained on \textsc{Tardis} radiative transfer models and is capable of a fast maximum likelihood parameter estimation and spectral fitting. After calculating the relevant physical parameters of supernovae we apply the expanding photosphere method to estimate their distances. Finally, we test the consistency of the obtained values by applying the formalism of Bayes factors. Results: The distances to four different host galaxies were estimated based on two supernovae in each. The distance estimates are not only consistent within the errors for each of the supernova sibling pairs, but in the case of two hosts they are precise to better than 5\%. Conclusions: Even though the literature data we used was not tailored for the requirements of our analysis, the agreement of the final estimates shows that the method is robust and is capable of inferring both precise and consistent distances. By using high-quality spectral time series, this method can provide precise distance estimates independent of the distance ladder, which are of high value for cosmology., Comment: 20 pages, 20 figures, 6 tables, Accepted in A&A

Even though the electromagnetic counterpart AT2017gfo to the binary neutron star merger GW170817 is powered by the radioactive decay of r-process nuclei, only few tentative identifications of light r-process elements have been made so far. One of the major limitations for the identification of heavy nuclei is incomplete or missing atomic data. While substantial progress has been made on lanthanide atomic data over the last few years, for actinides there has been less emphasis, with the first complete set of opacity data only recently published. We perform atomic structure calculations of neodymium $(Z=60)$ as well as the corresponding actinide uranium $(Z=92)$. Using two different codes (FAC and HFR) for the calculation of the atomic data, we investigate the accuracy of the calculated data (energy levels and electric dipole transitions) and their effect on kilonova opacities. For the FAC calculations, we optimise the local central potential and the number of included configurations and use a dedicated calibration technique to improve the agreement between theoretical and available experimental atomic energy levels (AELs). For ions with vast amounts of experimental data available, the presented opacities agree quite well with previous estimations. On the other hand, the optimisation and calibration method cannot be used for ions with only few available AELs. For these cases, where no experimental nor benchmarked calculations are available, a large spread in the opacities estimated from the atomic data obtained with the various atomic structure codes is observed.We find that the opacity of uranium is almost double the neodymium opacity., Comment: 20 pages, 13 figures. Accepted by 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., Comment: 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

We present the results from multi-wavelength observations of a transient discovered during the follow-up of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN2019wxt, a young transient in a galaxy whose sky position (in the 80\% GW contour) and distance ($\sim$150\,Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient's tightly constrained age, its relatively faint peak magnitude ($M_i \sim -16.7$\,mag) and the $r-$band decline rate of $\sim 1$\,mag per 5\,days appeared suggestive of a compact binary merger. However, SN2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of $\sim 0.1\,M_\odot$, with $^{56}$Ni comprising $\sim 20\%$ of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitors that could give rise to the observed properties of SN2019wxt, and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling electromagnetic counterparts to GW events from transients such as SN2019wxt is challenging: in a bid to characterise the level of contamination, we estimated the rate of events with properties comparable to those of SN2019wxt and found that $\sim 1$ such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500\,Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns., Comment: By the ENGRAVE collaboration (engrave-eso.org). 35 pages, 20 figures, final version accepted by A&A

A thermonuclear explosion triggered by a helium-shell detonation on a carbon-oxygen white dwarf core has been predicted to have strong UV line blanketing at early times due to the iron-group elements produced during helium-shell burning. We present the photometric and spectroscopic observations of SN 2016dsg, a sub-luminous peculiar Type I SN consistent with a thermonuclear explosion involving a thick He shell. With a redshift of 0.04, the $i$-band peak absolute magnitude is derived to be around -17.5. The object is located far away from its host, an early-type galaxy, suggesting it originated from an old stellar population. The spectra collected after the peak are unusually red, show strong UV line blanketing and weak O I $\lambda$7773 absorption lines, and do not evolve significantly over 30 days. An absorption line around 9700-10500 \AA is detected in the near-infrared spectrum and is likely from the unburnt helium in the ejecta. The spectroscopic evolution is consistent with the thermonuclear explosion models for a sub-Chandrasekhar mass white dwarf with a thick helium shell, while the photometric evolution is not well described by existing models., Comment: accepted for publication in ApJ

The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping., Comment: 13 pages. Accepted by MNRAS

We test the adequacy of ultraviolet (UV) spectra for characterizing the outer structure of Type Ia supernova (SN) ejecta. For this purpose, we perform spectroscopic analysis for ASASSN-14lp, a normal SN Ia showing low continuum in the mid-UV regime. To explain the strong UV suppression, two possible origins have been investigated by mapping the chemical profiles over a significant part of their ejecta. We fit the spectral time series with mid-UV coverage obtained before and around maximum light by HST, supplemented with ground-based optical observations for the earliest epochs. The synthetic spectra are calculated with the one dimensional MC radiative-transfer code TARDIS from self-consistent ejecta models. Among several physical parameters, we constrain the abundance profiles of nine chemical elements. We find that a distribution of $^{56}$Ni (and other iron-group elements) that extends toward the highest velocities reproduces the observed UV flux well. The presence of radioactive material in the outer layers of the ejecta, if confirmed, implies strong constraints on the possible explosion scenarios. We investigate the impact of the inferred $^{56}$Ni distribution on the early light curves with the radiative transfer code TURTLS, and confront the results with the observed light curves of ASASSN-14lp. The inferred abundances are not in conflict with the observed photometry. We also test whether the UV suppression can be reproduced if the radiation at the photosphere is significantly lower in the UV regime than the pure Planck function. In this case, solar metallicity might be sufficient enough at the highest velocities to reproduce the UV suppression., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

Manual fits to spectral times series of Type Ia supernovae have provided a method of reconstructing the explosion from a parametric model but due to lack of information about model uncertainties or parameter degeneracies direct comparison between theory and observation is difficult. In order to mitigate this important problem we present a new way to probabilistically reconstruct the outer ejecta of the normal Type Ia supernova SN 2002bo. A single epoch spectrum, taken 10 days before maximum light, is fit by a 13-parameter model describing the elemental composition of the ejecta and the explosion physics (density, temperature, velocity, and explosion epoch). Model evaluation is performed through the application of a novel rapid spectral synthesis technique in which the radiative transfer code, TARDIS, is accelerated by a machine-learning framework. Analysis of the posterior distribution reveals a complex and degenerate parameter space and allows direct comparison to various hydrodynamic models. Our analysis favors detonation over deflagration scenarios and we find that our technique offers a novel way to compare simulation to observation., Comment: 12 Pages, 4 Figures, submitted to AAS Journals, comments and constructive criticism welcome

A non-local-thermodynamic-equilibrium (NLTE) level population model of the first and second ionisation stages of iron, nickel and cobalt is used to fit a sample of XShooter optical + near-infrared (NIR) spectra of Type Ia supernovae (SNe Ia). From the ratio of the NIR lines to the optical lines limits can be placed on the temperature and density of the emission region. We find a similar evolution of these parameters across our sample. Using the evolution of the Fe II 12$\,$570$\,\mathring{A}\,$to 7$\,$155$\,\mathring{A}\,$line as a prior in fits of spectra covering only the optical wavelengths we show that the 7200$\,\mathring{A}\,$feature is fully explained by [Fe II] and [Ni II] alone. This approach allows us to determine the abundance of Ni II$\,$/$\,$Fe II for a large sample of 130 optical spectra of 58 SNe Ia with uncertainties small enough to distinguish between Chandrasekhar mass (M$_{\text{Ch}}$) and sub-Chandrasekhar mass (sub-M$_{\text{Ch}}$) explosion models. We conclude that the majority (85$\%$) of normal SNe Ia have a Ni/Fe abundance that is in agreement with predictions of sub-M$_{\text{Ch}}$ explosion simulations of $\sim Z_\odot$ progenitors. Only a small fraction (11$\%$) of objects in the sample have a Ni/Fe abundance in agreement with M$_{\text{Ch}}$ explosion models., Comment: 19 pages, 11 figures, Accepted for publication in MNRAS

This work presents the observations and analysis of ATLAS19dqr/SN 2019bkc, an extraordinary rapidly evolving transient event located in an isolated environment, tens of kiloparsecs from any likely host. Its light curves rise to maximum light in $5-6$ d and then display a decline of $\Delta m_{15} \sim5$ mag. With such a pronounced decay, it has one of the most rapidly evolving light curves known for a stellar explosion. The early spectra show similarities to normal and `ultra-stripped' type Ic SNe, but the early nebular phase spectra, which were reached just over two weeks after explosion, display prominent calcium lines, marking SN 2019bkc as a Ca-rich transient. The Ca emission lines at this phase show an unprecedented and unexplained blueshift of 10,000 -- 12,000 km/s. Modelling of the light curve and the early spectra suggests that the transient had a low ejecta mass of $0.2 - 0.4$ M$_\odot$ and a low kinetic energy of $ (2-4)\times 10^{50}$ erg, giving a specific kinetic energy $\sim1$ [$10^{51}$ erg]/M$_\odot$. The origin of this event cannot be unambiguously defined. While the abundance distribution used to model the spectra marginally favours a progenitor of white dwarf origin through the tentative identification of \ArII, the specific kinetic energy, which is defined by the explosion mechanism, is found to be more similar to an ultra-stripped core-collapse events. SN 2019bkc adds to the diverse range of physical properties shown by Ca-rich events., Comment: Accepted for publication in Astronomy & Astrophysics. Minor changes to section 4.3, some minor discussion added regarding opacities and line identification

We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks., Comment: 11 pages, 5 figures, 3 tables Accepted by ApJL on 2019 October 15

As a candidate 'super-Chandrasekhar' or 09dc-like Type Ia supernova (SN Ia), SN 2012dn shares many characteristics with other members of this remarkable class of objects but lacks their extraordinary luminosity. Here, we present and discuss the most comprehensive optical data set of this SN to date, comprised of a densely sampled series of early-time spectra obtained within the Nearby Supernova Factory project, plus photometry and spectroscopy obtained at the VLT about 1 yr after the explosion. The light curves, colour curves, spectral time series and ejecta velocities of SN 2012dn are compared with those of other 09dc-like and normal SNe Ia, the overall variety within the class of 09dc-like SNe Ia is discussed, and new criteria for 09dc-likeness are proposed. Particular attention is directed to additional insight that the late-phase data provide. The nebular spectra show forbidden lines of oxygen and calcium, elements that are usually not seen in late-time spectra of SNe Ia, while the ionisation state of the emitting iron plasma is low, pointing to low ejecta temperatures and high densities. The optical light curves are characterised by an enhanced fading starting ~60 d after maximum and very low luminosities in the nebular phase, which is most readily explained by unusually early formation of clumpy dust in the ejecta. Taken together, these effects suggest a strongly perturbed ejecta density profile, which might lend support to the idea that 09dc-like characteristics arise from a brief episode of interaction with a hydrogen-deficient envelope during the first hours or days after the explosion., Comment: 16 pages, 11 figures, MNRAS in press, missing line added in Table 1

We present the results of our photometric and spectroscopic follow-up of the intermediate-luminosity optical transient AT 2017jfs. At peak, the object reaches an absolute magnitude of Mg=-15.46+-0.15 mag and a bolometric luminosity of 5.5x10^41 erg/s. Its light curve has the double-peak shape typical of Luminous Red Novae (LRNe), with a narrow first peak bright in the blue bands, while the second peak is longer lasting and more luminous in the red and near-infrared (NIR) bands. During the first peak, the spectrum shows a blue continuum with narrow emission lines of H and Fe II. During the second peak, the spectrum becomes cooler, resembling that of a K-type star, and the emission lines are replaced by a forest of narrow lines in absorption. About 5 months later, while the optical light curves are characterized by a fast linear decline, the NIR ones show a moderate rebrightening, observed until the transient disappeared in solar conjunction. At these late epochs, the spectrum becomes reminiscent of that of M-type stars, with prominent molecular absorption bands. The late-time properties suggest the formation of some dust in the expanding common envelope or an IR echo from foreground pre-existing dust. We propose that the object is a common-envelope transient, possibly the outcome of a merging event in a massive binary, similar to NGC4490-2011OT1., Comment: 7 pages, 4 figures, 2 Tables (incl. one on-line). Published in A&A Letters

We obtained optical and near-infrared spectra of Type$\,$Ia supernovae (SNe$\,$Ia) at epochs ranging from 224 to 496 days after the explosion. The spectra show emission lines from forbidden transitions of singly ionised iron and cobalt atoms. We used non-local thermodynamic equilibrium (NLTE) modelling of the first and second ionisation stages of iron, nickel, and cobalt to fit the spectra using a sampling algorithm allowing us to probe a broad parameter space. We derive velocity shifts, line widths, and abundance ratios for iron and cobalt. The measured line widths and velocity shifts of the singly ionised ions suggest a shared emitting region. Our data are fully compatible with radioactive $^{56}$Ni decay as the origin for cobalt and iron. We compare the measured abundance ratios of iron and cobalt to theoretical predictions of various SN$\,$Ia explosion models. These models include, in addition to $^{56}$Ni, different amounts of $^{57}$Ni and stable $^{54,56}$Fe. We can exclude models that produced only $^{54,56}$Fe or only $^{57}$Ni in addition to $^{56}$Ni. If we consider a model that has $^{56}$Ni, $^{57}$Ni, and $^{54,56}$Fe then our data imply that these ratios are $^{54,56}$Fe / $^{56}$Ni $=0.272\pm0.086$ and $^{57}$Ni / $^{56}$Ni $=0.032\pm0.011$., Comment: 10 pages, 7 figures, Accepted for publication in A&A

We present near infrared (NIR) spectroscopy of the nearby supernova 2014J obtained $\sim$450 d after explosion. We detect the [Ni II] 1.939 $\mu$m line in the spectra indicating the presence of stable $^{58}$Ni in the ejecta. The stable nickel is not centrally concentrated but rather distributed as the iron. The spectra are dominated by forbidden [Fe II] and [Co II] lines. We use lines, in the NIR spectra, arising from the same upper energy levels to place constraints on the extinction from host galaxy dust. We find that that our data are in agreement with the high $A_V$ and low $R_V$ found in earlier studies from data near maximum light. Using a $^{56}$Ni mass prior from near maximum light $\gamma$-ray observations, we find $\sim$0.05 M$_\odot$ of stable nickel to be present in the ejecta. We find that the iron group features are redshifted from the host galaxy rest frame by $\sim$600 km s$^{-1}$., Comment: 6 pages, 4 figures, submitted to A&A

We investigate the light-curve properties of a sample of 26 spectroscopically confirmed hydrogen-poor superluminous supernovae (SLSNe-I) in the Palomar Transient Factory (PTF) survey. These events are brighter than SNe Ib/c and SNe Ic-BL, on average, by about 4 and 2~mag, respectively. The peak absolute magnitudes of SLSNe-I in rest-frame $g$ band span $-22\lesssim M_g \lesssim-20$~mag, and these peaks are not powered by radioactive $^{56}$Ni, unless strong asymmetries are at play. The rise timescales are longer for SLSNe than for normal SNe Ib/c, by roughly 10 days, for events with similar decay times. Thus, SLSNe-I can be considered as a separate population based on photometric properties. After peak, SLSNe-I decay with a wide range of slopes, with no obvious gap between rapidly declining and slowly declining events. The latter events show more irregularities (bumps) in the light curves at all times. At late times, the SLSN-I light curves slow down and cluster around the $^{56}$Co radioactive decay rate. Powering the late-time light curves with radioactive decay would require between 1 and 10${\rm M}_\odot$ of Ni masses. Alternatively, a simple magnetar model can reasonably fit the majority of SLSNe-I light curves, with four exceptions, and can mimic the radioactive decay of $^{56}$Co, up to $\sim400$ days from explosion. The resulting spin values do not correlate with the host-galaxy metallicities. Finally, the analysis of our sample cannot strengthen the case for using SLSNe-I for cosmology., Comment: 120 pages, 48 figures, 78 tables. ApJ in press