ABSTRACT
The radioactively powered transient following a binary neutron star merger, known as a kilonova (KN), is expected to enter the steady-state nebular phase a few days after merger. ...Steady-state holds until thermal reprocessing time-scales become long, at which point the temperature and ionization states need to be evolved time-dependently. We study the onset and significance of time-dependent effects using the non-local thermodynamic equilibrium spectral synthesis code sumo. We employ a simple single-zone model with an elemental composition of Te, Ce, Pt, and Th, scaled to their respective solar abundances. The atomic data are generated using the Flexible Atomic Code (fac), and consist of energy levels and radiative transitions, including highly forbidden lines. We explore the KN evolution from 5 to 100 d after merger, varying ejecta mass and velocity. We also consider variations in the degree of electron magnetic field trapping, as well as radioactive power generation for alpha and beta decay (but omitting fission products). We find that the transition time, and magnitude of steady-state deviations are highly sensitive to these parameters. For typical KN ejecta, the deviations are minor within the time-frame studied. However, low density ejecta with low energy deposition show significant differences from ∼10 d. Important deviation of the ionization structure solution impacts the temperature by altering the overall line cooling. Adiabatic cooling becomes important at t ≥ 60 d which, in addition to the temperature and ionization effects, lead to the bolometric light-curve deviating from the instantaneous radioactive power deposited.
NLTE effects on kilonova expansion opacities Pognan, Quentin; Jerkstrand, Anders; Grumer, Jon
Monthly notices of the Royal Astronomical Society,
05/2022, Letnik:
513, Številka:
4
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
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 several 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 (less than a week), while highlighting the need for NLTE modelling during later phases. Considering also NLTE ionization, our results indicate that NLTE can give both higher or lower opacities, depending on composition and wavelength, sometimes by orders of magnitudes.
ABSTRACT The observational effects of the "Infrared Catastrophe" are discussed in view of the very late observations of the Type Ia SN 2011fe. Our model spectra at 1000 days take non-local radiative ...transfer into account and find that this has a crucial impact on the spectral formation. Although rapid cooling of the ejecta to a few 100 K occurs also in these models, the late-time optical/NIR flux is brighter by 1-2 mag due to redistribution of UV emissivity, resulting from non-thermal excitation and ionization. This effect brings models into better agreement with late-time observations of SN 2011fe, and other SNe Ia, and offers a solution to the long-standing discrepancy between models and observations. The models show that spectral formation shifts from Fe ii and Fe iii at 300 days to Fe i at 1000 days, which explains the apparent wavelength shifts seen in SN 2011fe. We discuss the effects of time dependence and energy input from 57Co, finding both to be important at 1000 days.
Superluminous supernovae radiate up to 100 times more energy than normal supernovae. The origin of this energy and the nature of the stellar progenitors of these transients are poorly understood. We ...identify neutral iron lines in the spectrum of one such supernova, SN 2006gy, and show that they require a large mass of iron (≳0.3 solar masses) expanding at 1500 kilometers per second. By modeling a standard type Ia supernova hitting a shell of circumstellar material, we produce a light curve and late-time iron-dominated spectrum that match the observations of SN 2006gy. In such a scenario, common envelope evolution of a progenitor binary system can synchronize envelope ejection and supernova explosion and may explain these bright transients.
ABSTRACT
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.
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 MV = −16 mag and has H expansion velocities under 1350 km s−1 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.
Despite more than 30 years of searching, 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 ( erg s−1 cm−2 Hz−1) at 213 GHz, 1 L ( erg s−1 cm−2 Hz−1) in the optical if our line of sight is free of ejecta dust, and 1036 erg s−1 ( 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. The allowed bolometric luminosity of the compact object is 22 L if our line of sight is free of ejecta dust, or 138 L if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star (NS) to -8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency is limited to M yr−1, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength (B) for a given spin period (P) is G s−2, which firmly excludes pulsars comparable to the Crab. By combining information about radiation reprocessing and geometry, we infer that the compact object is a dust-obscured thermally emitting NS, which may appear as a region of higher-temperature ejecta dust emission.
We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of Mr = −17.98, a peak bolometric ...luminosity of ∼1.4 × 1043 erg s−1, and a rapid decline rate of 0.48 0.03 mag day−1 in the r band. The bolometric luminosity evolves too quickly to be explained by pure 56Ni 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 Mej = 0.10 0.05 M . Our spectral modeling suggests the ejecta is composed of intermediate mass elements including O, Si, and Mg and trace amounts of Fe-peak elements, which disfavors a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition, and high likelihood of additional powering-the core collapse of an ultra-stripped progenitor, the accretion induced collapse (AIC) of a white dwarf, and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favors either the AIC with magnetar powering or a white dwarf-neutron star merger with energy from disk wind shocks.
Abstract
Hydrogen-rich, core-collapse supernovae are typically divided into four classes: IIP, IIL, IIn, and IIb. Recent hydrodynamic modelling shows that circumstellar material is required to ...produce the early light curves of most IIP/IIL supernovae. In this scenario, IIL supernovae experience large amounts of mass-loss before exploding. We test this hypothesis on ASASSN-15oz, a Type IIL supernova. With extensive follow-up in the X-ray, UV, optical, IR, and radio, we present our search for signs of interaction and the mass-loss history indicated by their detection. We find evidence of short-lived intense mass-loss just prior to explosion from light-curve modelling, amounting in 1.5 M⊙ of material within 1800 R⊙ of the progenitor. We also detect the supernova in the radio, indicating mass-loss rates of 10−6 to 10−7 M⊙ yr−1 prior to the extreme mass-loss period. Our failure to detect the supernova in the X-ray and the lack of narrow emission lines in the UV, optical, and NIR do not contradict this picture and place an upper limit on the mass-loss rate outside the extreme period of <10−4 M⊙ yr−1. This paper highlights the importance gathering comprehensive data on more Type II supernovae to enable detailed modelling of the progenitor and supernova which can elucidate their mass-loss histories and envelope structures and thus inform stellar evolution models.
ABSTRACT
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, an ...LL SN II with extensive spectral and photometric coverage. We measure a V-band peak magnitude of −14.58 mag, a plateau duration of ∼100 d, and an inferred 56Ni mass of 0.008 ± 0.002 M⊙. The peak bolometric luminosity, Lbol ≈ 1041.4 erg s−1, and its spectral evolution are typical of other SNe in the class. Using our late-time spectra, we measure the O i λλ6300, 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 ± 3 M⊙. Our extensive late-time spectral coverage of the Fe ii λ7155 and Ni ii λ7378 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 ∼200–300 d 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 that 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.