Observations Of SN 1987A Research Articles

Insights into the Production of 44Ti and Nickel Isotopes in Core-collapse Supernovae

We report nucleosynthetic results for both 44Ti and nickel isotopes for 18 three-dimensional (3D) core-collapse supernova (CCSN) simulations extended to ∼20 s after bounce. We find that many of our long-term models are able to achieve 44Ti/56Ni ratios similar to that observed in Cassiopeia A, and modern supernova models can synthesize up to 2 × 10−4 M ⊙ of 44Ti. Neutrino-driven winds and the fact that there can be simultaneous accretion and explosion in 3D models of CCSNe play central roles in its production. We conclude that the 44Ti underproduction problem in previous CCSN models is no longer an issue. In addition, we discuss the production of both 57Ni and stable nickel/iron ratios and compare our results to observations of SN 1987A and the Crab.

Investigating the Time Evolution of the Thermal Emission from the Putative Neutron Star in SN 1987A for 50+ Years

Observations collected with the Atacama Large Millimeter/submillimeter Array (ALMA) and analysis of broadband X-ray spectra have recently suggested the presence of a central compact object (CCO) in SN 1987A. However, no direct evidence of the CCO has been found yet. Here we analyze Chandra X-ray observations of SN 1987A collected in 2007 and 2018, and synthesize 2027 Chandra and 2037 Lynx spectra of the faint inner region of SN 1987A. We estimate the temporal evolution of the upper limits of the intrinsic luminosity of the putative CCO in three epochs (2018, 2027, and 2037). We find that these upper limits are higher for higher neutron star (NS) kick velocities due to increased absorption from the surrounding cold ejecta. We compare NS cooling models with both the intrinsic luminosity limits obtained from the X-ray spectra and the ALMA constraints with the assumption that the observed blob of SN 1987A is primarily heated by thermal emission. We find that the synthetic Lynx spectra are crucial to constrain the physical properties of the CCO, which will be confirmed by future observations in the 2040s. We draw our conclusions based on two scenarios, namely the nondetection and detection of the NS by Lynx. If the NS is not detected, its kick velocity should be ≃700 km s−1. Furthermore, nondetection of the NS would suggest rapid cooling processes at the age of 40 yr, implying strong crust superfluidity. Conversely, in the case of NS detection, the mass of the NS envelope must be high.

Supernova 1987A: 3D Mixing and Light Curves for Explosion Models Based on Binary-merger Progenitors

Abstract Six binary-merger progenitors of supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak −69°202 are exploded by neutrino heating and evolved until long after shock breakout in 3D and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize 56Ni masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km s−1 compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive 56Ni mixing depends mainly on high growth factors of Rayleigh–Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak −69°202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first ∼7 days. Models that mix about 0.5 M ⊙ of hydrogen into the He-shell and exhibit strong outward mixing of 56Ni with maximum velocities exceeding the 3000 km s−1 observed for the bulk of ejected 56Ni have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one.

Testing the ALP-photon coupling with polarization measurements of Sagittarius A⋆

Ultra-light bosons such as axions or axion-like particles (ALPs), are promising candidates to solve the dark matter problem. A unique way to detect such ALPs is to search for the periodic oscillation feature of the position angles of linearly polarized photons emitted from the galaxy center. In this work, we use the high-resolution polarimetric measurements of the radiation near the super-massive black hole (SMBH) in the center of the Milky Way, i.e., Sagittarius A⋆ (Sgr A⋆), by a sub-array of the Event Horizon Telescope to search for the ultra-light ALPs. We derive upper limits on the ALP-photon coupling of ∼ 10-12 GeV-1 for ALP masses of ma∼ (10-19-10-18) eV, with a solitonic core + NFW dark matter density profile. Our results are stronger than that derived from the observations of SN1987A and a population of supernovae in the mass window of (10-19-10-17) eV. Improved polarimetric measurements with the full Event Horizon Telescope can further strengthen the constraints.

An Alternative to Dark Matter? Part 1: The Early Universe (<i>t<sub>p</sub></i> to 10<sup>-9</sup> s), Energy Creation the Alphaton, Baryogenesis

A cosmological model was developed using the equation of state of photon gas, as well as cosmic time. The primary objective of this model is to see if determining the observed rotation speed of galactic matter is possible, without using dark matter (halo) as a parameter. To do so, a numerical application of the evolution of variables in accordance with cosmic time and a new state equation was developed to determine precise, realistic values for a number of cosmological parameters, such as the energy of the universe U, cosmological constant Λ, the curvature of space k, energy density ρΛe, age of the universe tΩ etc. The development of the state equation highlights the importance of not neglecting any of the differential terms given the very large amounts in play that can counterbalance the infinitesimals. Some assumptions were put forth in order to solve these equations. The current version of the model partially explains several of the observed phenomena that raise questions. Numerical application of the model has yielded the following results, among others: Initially, during the Planck era, at the very beginning of Planck time, tp, the universe contained a single photon at Planck temperature TP, almost Planck energy EP in the Planck volume. During the photon inflation phase (before characteristic time ~10-9 [s]), the number of original photons (alphatons) increased at each unit of Planck time tp and geometrical progression~n3, where n is the quotient of cosmic time over Planck time t/tp. Then, the primordial number of photons reached a maximum of N~1089, where it remained constant. These primordial photons (alphatons) are still present today and represent the essential of the energy contained in the universe via the cosmological constant expressed in the form of energy EΛ. Such geometric growth in the number of photons can bring a solution to the horizon problem through γγ exchange and a photon energy volume that is in phase with that of the volume energy of the universe. The predicted total mass (p, n, e, and ν), based on the Maxwell-Juttner relativistic statistical distribution, is ~7 × 1050 [kg]. The predicted cosmic neutrino mass is ≤8.69 × 10-32 [kg] (≤48.7 [keV·c-2]) if based on observations of SN1987A. The temperature variation of the cosmic microwave background (CMB), as measured by Planck, can be said to be partially due to energy variations in the universe (ΔU/U) during the primordial baryon synthesis (energy jump from the creation of protons and neutrons).

Properties of gamma-ray decay lines in 3D core-collapse supernova models, with application to SN 1987A and Cas A

ABSTRACT Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometres per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For MZAMS = 15−20 M⊙ progenitors exploding with ∼1.5 × 1051 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ∼400 km s−1 and a bulk velocity of at least 1500 km s−1. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600 d, an ejecta mass around 14 M⊙ is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint VNS > Vredshift, and applying this to SN 1987A gives a NS kick of at least 500 km s−1. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

Gamma-ray spectral modulations induced by photon-ALP-dark photon oscillations

Recently it has been noticed that the Fermi-LAT data of gamma-rays from some galactic pulsars and supernova remnants reveal spectral modulations that might be explained by the conversion of photons to ALPs (axion-like particles) induced by the conventional ALP coupling to photon in the presence of galactic magnetic fields. However the corresponding ALP mass and coupling are in a severe tension with the observational constraints from CAST, SN1987A, and other gamma-ray observations. Motivated by this, we examine an alternative possibility that those spectral modulations are explained by other type of ALP coupling involving both the ordinary photon and a massless dark photon, when nonzero background dark photon gauge fields are assumed. We find that our scheme results in oscillations among the photon, ALP, and dark photon, which can explain the gamma-ray spectral modulations of galactic pulsars or supernova remnants, while satisfying the known observational constraints.

Final Spitzer IRAC Observations of the Rise and Fall of SN 1987A

Abstract Spitzer's final Infrared Array Camera observations of SN 1987A show the 3.6 and 4.5 μm emission from the equatorial ring (ER) continues a period of steady decline. Deconvolution of the images reveals that the emission is dominated by the ring, not the ejecta, and is brightest on the west side. Decomposition of the marginally resolved emission also confirms this, and shows that the west side of the ER has been brightening relative to the other portions of the ER. The infrared morphological changes resemble those seen in both the soft X-ray emission and the optical emission. The integrated ER light curves at 3.6 and 4.5 μm are more similar to the optical light curves than the soft X-ray light curve, though differences would be expected if dust is responsible for this emission and its destruction is rapid. Future observations with the James Webb Space Telescope will continue to monitor the ER evolution, and will reveal the true spectrum and nature of the material responsible for the broadband emission at 3.6 and 4.5 μm. The present observations also serendipitously reveal a nearby variable source, subsequently identified as a Be star, that has gone through a multiyear outburst during the course of these observations.

Observing Supernova Neutrino Light Curves with Super-Kamiokande: Expected Event Number over 10 s

Abstract Supernova neutrinos are crucially important to probe the final phases of massive star evolution. As is well known from observations of SN 1987A, neutrinos provide information on the physical conditions responsible for neutron star formation and on the supernova explosion mechanism. However, there is still no complete understanding of the long-term evolution of neutrino emission in supernova explosions, although there are a number of modern simulations of neutrino radiation hydrodynamics, which study neutrino emission at times less than one second after the bounce. In the present work we systematically calculate the number of neutrinos that can be observed in Super-Kamiokande over periods longer than 10 seconds using the database of Nakazato et al. anticipating that neutrinos from a Galactic supernova can be detected for several tens of seconds. We find that for a supernova at a distance of 10 kpc, neutrinos remain observable for longer than 30 s for a low-mass neutron star (1.20 M ⊙ gravitational mass) and even longer than 100 s for a high-mass neutron star (2.05 M ⊙). These scenarios are much longer than the observations of SN 1987A and longer than the duration of existing numerical simulations. We propose a new analysis method based on the cumulative neutrino event distribution as a function of reverse time from the last observed event, as a useful probe of the neutron star mass. Our result demonstrates the importance of complete modeling of neutrino light curves in order to extract physical quantities essential for understanding supernova explosion mechanisms, such as the mass and radius of the resulting neutron star.

Observable signatures of dark photons from supernovae

A dark photon is a well-motivated new particle which, as a component of an associated dark sector, could explain dark matter. One strong limit on dark photons arises from excessive cooling of supernovae. We point out that even at couplings where too few dark photons are produced in supernovae to violate the cooling bound, they can be observed directly through their decays. Supernovae produce dark photons which decay to positrons, giving a signal in the 511 keV annihilation line observed by SPI/INTEGRAL. Further, prompt gamma-ray emission by these decaying dark photons gives a signal for gamma-ray telescopes. Existing GRS observations of SN1987a already constrain this, and a future nearby SN could provide a detection. Finally, dark photon decays from extragalactic SN would produce a diffuse flux of gamma rays observable by detectors such as SMM and HEAO-1. Together these observations can probe dark photon couplings several orders of magnitude beyond current constraints for masses of roughly 1-100 MeV.

Very Deep inside the SN 1987A Core Ejecta: Molecular Structures Seen in 3D

Abstract Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outward through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks (“nickel heating”). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained.

High-Resolution Observations of Dust in SN 1987A

AbstractThe dust produced by supernovae is an important topic for understanding supernova physics and the chemical evolution of galaxies. Recent ALMA observations of SN 1987A have allowed us to peer into the inner ejecta to the cool dust, with spatial resolution from 0.″3 at ~300 GHz down to 0.″09 at ~680 GHz – an improvement over the previous 300 GHz Cycle 0 observations at 0.″69. Comparison of the dust location and morphology with other multiwavelength emission presents an interesting picture of the role dust plays in the ejecta. The mm-FIR SED is compared to radiative models to study the dust composition 30 years after the initial explosion. Fits to the ring emission also probe the drift of the center of the system over time.

ALMA observations of Molecules in Supernova 1987A

AbstractSupernova (SN) 1987A has provided a unique opportunity to study how SN ejecta evolve in 30 years time scale. We report our ALMA spectral observations of SN 1987A, taken in 2014, 2015 and 2016, with detections of CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO.We find a dip in the SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities causes mixing of gas, with heavier elements much more disturbed, making more elongated structure.Using 28SiO and its isotopologues, Si isotope ratios were estimated for the first time in SN 1987A. The estimated ratios appear to be consistent with theoretical predictions of inefficient formation of neutron rich atoms at lower metallicity, such as observed in the Large Magellanic Cloud (about half a solar metallicity).The deduced large HCO+ mass and small SiS mass, which are inconsistent to the predictions of chemical model, might be explained by some mixing of elements immediately after the explosion. The mixing might have made some hydrogen from the envelope to sink into carbon and oxygen-rich zone during early days after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may penetrate into silicon and sulphur zone, suppressing formation of SiS.Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae.

THREE-DIMENSIONAL DISTRIBUTION OF EJECTA IN SUPERNOVA 1987A AT 10,000 DAYS

ABSTRACT Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ∼10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of Hα to date, the first 3D maps for [Ca ii] , [O i] , and Mg ii , as well as new maps for [Si i]+[Fe ii] and He i 2.058 μm. A comparison with previous observations shows that the [Si i]+[Fe ii] flux and morphology have not changed significantly during the past ten years, providing evidence that this line is powered by 44Ti. The time evolution of Hα shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north–south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, Hα and [Si i]+[Fe ii] 1.644 μm, show substructures at the level of ∼200–1000 km s–1and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.

CHANDRA OBSERVES THE END OF AN ERA IN SN 1987A

ABSTRACT Updated imaging and photometric results from Chandra observations of SN 1987A, covering the last 16 years, are presented. We find that the 0.5–2 keV light curve has remained constant at ∼8 × 10−12 erg s−1 cm−2 since 9500 days, with the 3–8 keV light curve continuing to increase until at least 10,000 days. The expansion rate of the ring is found to be energy dependent, such that after day 6000 the ring expands faster in the 2–10 keV band than it does at energies <2 keV. Images show a reversal of the east–west asymmetry between 7000 and 8000 days after the explosion. The latest images suggest the southeastern side of the equatorial ring (ER) is beginning to fade. Consistent with the latest optical and infrared results, our Chandra analysis indicates the blast wave is now leaving the dense ER, which marks the beginning of a major change in the evolutionary phase of the supernova remnant 1987A.

INFRARED CONTINUUM AND LINE EVOLUTION OF THE EQUATORIAL RING AROUND SN 1987A

ABSTRACT Spitzer observations of SN 1987A have now spanned more than a decade. Since day ∼4000, mid-infrared (mid-IR) emission has been dominated by that from shock-heated dust in the equatorial ring (ER). From 6000 to 8000 days after the explosion, Spitzer observations included broadband photometry at 3.6–24 μm, and low and moderate resolution spectroscopy at 5–35 μm. Here we present later Spitzer observations, through day 10,377, which include only the broadband measurements at 3.6 and 4.5 μm. These data show that the 3.6 and 4.5 μm brightness has clearly begun to fade after day ∼8500, and no longer tracks the X-ray emission as well as it did at earlier epochs. This can be explained by the destruction of the dust in the ER on timescales shorter than the cooling time for the shocked gas. We find that the evolution of the late time IR emission is also similar to the now fading optical emission. We provide the complete record of the IR emission lines, as seen by Spitzer prior to day 8000. The past evolution of the gas as seen by the IR emission lines seems largely consistent with the optical emission, although the IR [Fe ii] and [Si ii] lines show different, peculiar velocity structures.

THE INTERPLAY BETWEEN CHEMISTRY AND NUCLEATION IN THE FORMATION OF CARBONACEOUS DUST IN SUPERNOVA EJECTA

ABSTRACT Core-collapse supernovae (CCSNe) are considered to be important contributors to the primitive dust enrichment of the interstellar medium in the high-redshift universe. Theoretical models of dust formation in stellar explosions have so far provided controversial results and a generally poor fit to the observations of dust formation in local supernovae. We present a new methodology for the calculation of carbonaceous dust formation in young supernova remnants. Our new technique uses both the nucleation theory and a chemical reaction network to allow us to compute the dust growth beyond the molecular level as well as consider the chemical erosion of the forming grains. We find that carbonaceous dust forms efficiently in the core of the ejecta, but takes several years to condensate, longer than previously estimated. It forms unevenly and remains concentrated in the inner part of the remnant. These results support the role of CCSNe as dust factories and provide new insight into the observations of SN 1987A, in which large amounts of dust have been detected to form on a timescale of years after core-collapse.

Supernova 1987A: neutrino-driven explosions in three dimensions and light curves

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

The extinction law inside the 30 Doradus nebula★

We have studied the interstellar extinction in a field of ~3' x 3' at the core of the 30 Doradus nebula, including the central R136 cluster, in the Large Magellanic Cloud. Observations at optical and near-infrared wavelengths, obtained with the WFC3 camera on board the Hubble Space Telescope, show that the stars belonging to the red giant clump are spread across the colour-magnitude diagrams because of the considerable and uneven levels of extinction in this region. Since these stars share very similar physical properties and are all at the same distance, they allow us to derive the absolute extinction in a straightforward and reliable way. Thus we have measured the extinction towards about 180 objects and the extinction law in the range 0.3 - 1.6 micron. At optical wavelengths, the extinction curve is almost parallel to that of the diffuse Galactic interstellar medium. Taking the latter as a template, the value of Rv = 4.5 +/- 0.2 that we measure indicates that in the optical there is an extra grey component due to a larger fraction of large grains. At wavelengths longer than ~1 micron, the contribution of this additional component tapers off as lambda^-1.5, like in the Milky Way, suggesting that the nature of the grains is otherwise similar to those in our Galaxy, but with a ~2.2 times higher fraction of large grains. These results are consistent with the addition of "fresh" large grains by supernova explosions, as recently revealed by Herschel and ALMA observations of SN 1987A.

MATTER MIXING IN ASPHERICAL CORE-COLLAPSE SUPERNOVAE: A SEARCH FOR POSSIBLE CONDITIONS FOR CONVEYING56Ni INTO HIGH VELOCITY REGIONS

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.