Search

Natal supernova kicks, the linear momentum compact remnants receive during their formation, are an essential part of binary population synthesis (BPS) models. Although these kicks are well-supported by evidence, their underlying distributions and incorporation into BPS models is uncertain. In this work, we investigate the nature of natal kicks using a previously proposed analytical prescription where the strength of the kick is given by $v_\text{k}=\alpha\frac{m_\text{ejecta}}{m_\text{remnant}}+\beta~\text{km s}^{-1}$, for free parameters $\alpha$ and $\beta$. We vary the free parameters over large ranges of possible values, comparing these synthetic populations simultaneously against four constraints: the merger rate of compact binary neutron star (BNS) systems, the period-eccentricity distribution of galactic BNSs, the velocity distribution of single-star pulsars, and the likelihood for low-ejecta mass supernovae to produce low-velocity kicks. We find that different samples of the parameter space satisfy each tests, and only 1 per cent of the models satisfy all four constraints simultaneously. Although we cannot identify a single best kick model, we report $\alpha=115^{+40}_{-55}~\text{km s}^{-1}, \beta=15^{+10}_{-15}~\text{km s}^{-1}$ as the center of the region of the parameter space that fulfils all of our constraints, and expect $\beta\geq0~\text{km s}^{-1}$ as a further constraint. We also suggest further observations that will enable future refinement of the kick model. A sensitive test for the kick model will be the redshift evolution of the BNS merger rate since this is effectively a direct measure of the delay-time distribution for mergers. For our best fitting values, we find that the peak of the BNS merger rate is the present-day., Comment: Accepted by MNRAS

S. J. Brennan et al., We present the bolometric light curve, identification and analysis of the progenitor candidate, and preliminary modelling of AT 2016jbu (Gaia16cfr). We find a progenitor consistent with a ∼ 22–25 M⊙ yellow hypergiant surrounded by a dusty circumstellar shell, in agreement with what has been previously reported. We see evidence for significant photometric variability in the progenitor, as well as strong Hα emission consistent with pre-existing circumstellar material. The age of the environment, as well as the resolved stellar population surrounding AT 2016jbu, supports a progenitor age of >10 Myr, consistent with a progenitor mass of ∼22 M⊙. A joint analysis of the velocity evolution of AT 2016jbu and the photospheric radius inferred from the bolometric light curve shows the transient is consistent with two successive outbursts/explosions. The first outburst ejected material with velocity ∼650 km s−1, while the second, more energetic event ejected material at ∼4500 km s−1. Whether the latter is the core collapse of the progenitor remains uncertain. We place a limit on the ejected 56Ni mass of, SJB acknowledges support from Science Foundation Ireland and the Royal Society (RS-EA/3471). MF is supported by a Royal Society–Science Foundation Ireland University Research Fellowship. TMB was funded by CONICYT PFCHA/DOCTORADOBECAS CHILE/2017-72180113 and acknowledges financial support from the Spanish Ministerio de Ciencia e Innovación (MCIN), the Agencia Estatal de Investigación (AEI) 10.13039/501100011033 under the PID2020-115253GA-I00 HOSTFLOWS project, and Centro Superior de Investigaciones Científicas (CSIC) under the PIE project 20215AT016. KM is funded by the EU H2020 ERC grant no. 758638. TWC acknowledges EU Funding under Marie Skłodowska–Curie grant H2020-MSCA-IF-2018-842471, and thanks Thomas Krühler for GROND data reduction. MN is supported by a Royal Astronomical Society Research Fellowship. BJS is supported by NSF grants AST-1908952, AST-1920392, AST-1911074, and NASA award 80NSSC19K1717. MS is supported by generous grants from Villum FONDEN (13261,28021) and by a project grant (8021-00170B) from the Independent Research Fund Denmark. LH acknowledges support for Watcher from Science Foundation Ireland grant 07/RFP/PHYF295. Time-domain research by DJS is supported by NSF grants AST-1821987, 1813466, and 1908972, and by the Heising–Simons Foundation under grant #2020-1864. NER acknowledges support from MIUR, PRIN 2017 (grant 20179ZF5KS). Support for JLP is provided in part by ANID through the Fondecyt regular grant 1191038 and through the Millennium Science Initiative grant ICN12_009, awarded to The Millennium Institute of Astrophysics, MAS. LG acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MICIU) under the 2019 Ramón y Cajal program RYC2019-027683 and from the Spanish MICIU project PID2020-115253GA-I00. DAH and DH are supported by AST-1911151, AST19-11225, and NASA Swift grant 80NSSC19K1639. GP acknowledges support by the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. LT acknowledges support from MIUR (PRIN 2017 grant 20179ZF5KS). Support for TW-SH was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51458.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. HK was funded by the Academy of Finland projects 324504 and 328898. This research made use of ASTROPY,8 a community-developed core PYTHON package for Astronomy (Astropy Collaboration 2013; Price-Whelan et al. 2018). This research made use of data provided by Astrometry.net.9 Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in Three Dimensions (ASTRO 3D), through project number CE170100013. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We acknowledge the Telescope Access Program (TAP) funded by the NAOC, CAS, and the Special Fund for Astronomy from the Ministry of Finance. This work was partially supported by Polish NCN grants: Harmonia No. 2018/30/M/ST9/00311 and Daina No. 2017/27/L/ST9/03221. This work made use of v2.2.1 of the Binary Population and Spectral Synthesis (BPASS) models as described in Eldridge et al. (2017) and Stanway & Eldridge (2018). This research is based on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program 15645. Observations were also obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESAC/ESA), and the Canadian Astronomy Data Centre (CADC/NRC/CSA). This research has made use of the SVO Filter Profile Service10 supported from the Spanish MINECO through grant AYA2017-84089.

Machine learning has become widely used in astronomy. Gaussian Process (GP) regression in particular has been employed a number of times to fit or re-sample supernova (SN) light-curves, however by their nature typical GP models are not suited to fit SN photometric data and they will be prone to over-fitting. Recently GP re-sampling was used in the context of studying the morphologies of type II and IIb SNe and they were found to be clearly distinct with respect to four parameters: the rise time (t$_{\rm rise}$), the magnitude difference between 40 and 30 days post explosion ($\Delta m_{\rm 40-30}$), the earliest maximum (post-peak) of the first derivative (dm1) and minimum of the second derivative (dm2). Here we take a close look at GP regression and its limitations in the context of SN light-curves in general, and we also discuss the uncertainties on these specific parameters, finding that dm1 and dm2 cannot give reliable astrophysical information. We do reproduce the clustering in t$_{\rm rise}$--$\Delta m_{\rm 40-30}$ space although it is not as clear cut as previously presented. The best strategy to accurately populate the t$_{\rm rise}$-- $\Delta m_{\rm 40-30}$ space will be to use an expanded sample of high quality light-curves (such as those in the ATLAS transient survey) and analytical fitting methods. Finally, using the BPASS fiducial models, we predict that future photometric studies will reveal clear clustering of the type IIb and II light curve morphologies with a distinct continuum of transitional events., Comment: 13 pages, 11 figures, 2 tables, Accepted for publication in MNRAS

Brennan, S. J. et al (ASTROPY Collaboration), We present the results from a high-cadence, multiwavelength observation campaign of AT 2016jbu (aka Gaia16cfr), an interacting transient. This data set complements the current literature by adding higher cadence as well as extended coverage of the light-curve evolution and late-time spectroscopic evolution. Photometric coverage reveals that AT 2016jbu underwent significant photometric variability followed by two luminous events, the latter of which reached an absolute magnitude of MV ∼-18.5 mag. This is similar to the transient SN 2009ip whose nature is still debated. Spectra are dominated by narrow emission lines and show a blue continuum during the peak of the second event. AT 2016jbu shows signatures of a complex, non-homogeneous circumstellar material (CSM). We see slowly evolving asymmetric hydrogen line profiles, with velocities of 500 km s-1 seen in narrow emission features from a slow-moving CSM, and up to 10 000 km s-1 seen in broad absorption from some high-velocity material. Late-time spectra (∼+1 yr) show a lack of forbidden emission lines expected from a core-collapse supernova and are dominated by strong emission from H, He i, and Ca ii. Strong asymmetric emission features, a bumpy light curve, and continually evolving spectra suggest an inhibit nebular phase. We compare the evolution of H α among SN 2009ip-like transients and find possible evidence for orientation angle effects. The light-curve evolution of AT 2016jbu suggests similar, but not identical, circumstellar environments to other SN 2009ip-like transients., SJB acknowledges support from Science Foundation Ireland and the Royal Society (RS-EA/3471). MF is supported by a Royal Society - Science Foundation Ireland University Research Fellowship. TMB was funded by the CONICYT PFCHA / DOCTORADOBECAS CHILE/2017-72180113 and acknowledges their financial support from the Spanish Ministerio de Ciencia e Innovación (MCIN), the Agencia Estatal de Investigación (AEI) 10.13039/501100011033 under the PID2020-115253GA-I00 HOSTFLOWS project, and from Centro Superior de Investigaciones Científicas (CSIC) under the PIE project 20215AT016. KM is funded by the EU H2020 ERC grant no. 758638. TWC acknowledges the EU Funding under Marie Skłodowska-Curie grant H2020-MSCA-IF-2018-842471, and thanks to Thomas Krühler for GROND data reduction. MN is supported by a Royal Astronomical Society Research Fellowship. BJS is supported by NSF grants AST-1908952, AST-1920392, AST-1911074, and NASA award 80NSSC19K1717. MS is supported by generous grants from Villum FONDEN (13261, 28021) and by a project grant (8021-00170B) from the Independent Research Fund Denmark. LH acknowledges support for Watcher from Science Foundation Ireland grant 07/RFP/PHYF295. Time domain research by DJS is supported by NSF grants AST-1821987, 1813466, and 1908972, and by the Heising-Simons Foundation under grant #2020-1864. NER acknowledges support from MIUR, PRIN 2017 (grant 20179ZF5KS). LG acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MICIU) under the 2019 Ramón y Cajal program RYC2019-027683 and from the Spanish MICIU project PID2020-115253GA-I00. Support for TW-SH was provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51458.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. Support for JLP is provided in part by ANID through the Fondecyt regular grant 1191038 and through the Millennium Science Initiative grant ICN12_009, awarded to The Millennium Institute of Astrophysics, MAS. DAH and DH are supported by AST-1911151, AST19-11225, and NASA Swift grant 80NSSC19K1639. GP acknowledge support by the Ministry of Economy, Development, and Tourism’s Millennium Science Initiative through grant IC120009, awarded to The Millennium Institute of Astrophysics, MAS. LT acknowledges support from MIUR (PRIN 2017 grant 20179ZF5KS). HK was funded by the Academy of Finland projects 324504 and 328898.

With the remarkable success of the LVK consortium in detecting binary black hole mergers, it has become possible to use the population properties to constrain our understanding of the progenitor stars' evolution. The most striking features of the observed primary black hole mass distributions are the extended tail up to 100M$_\odot$ and an excess of masses at 35M$_\odot$. Currently, isolated binary population synthesis have difficulty explaining these features. Using the well-tested BPASS detailed stellar binary evolution models to determine mass transfer stability, accretion rates, and remnant masses, we postulate that stable mass transfer with super-Eddington accretion is responsible for the extended tail. Furthermore, that the excess is not due to pulsation-pair instability, as previously thought, but due to stable mass transfer. These systems are able to merge within the Hubble time due to more stable mass transfer with extreme mass ratios that allows the orbits to shrink sufficiently to allow for a merger. These finding are at odds with those from other population synthesis codes but in agreement with other recent studies using detailed binary evolution models., Comment: Submitted to MNRAS, comments welcome. 22 pages, 18 figures, 9 pages supplementary material

The recent discovery of an unambiguous quiescent BH and main sequence O star companion in VFTS 243 opens the door to new constraints on theoretical stellar evolution and population models looking to reproduce the progenitors of black hole - black hole binaries. Here we show that the Binary Population and Spectral Synthesis fiducial models (BPASSv2.2.1) natively predict VFTS243-like systems: We find that VFTS 243 likely originated from a binary system in a about 15 day orbit with primary mass ranging from 40 to 50 M$_{\odot}$ and secondary star with initial mass 24-25M$_{\odot}$. Additionally we find that the death of the primary star must have resulted in a low energy explosion E < 10$^{50}$ ergs.With a uniform prior we find that the kick velocity of the new-born black hole was, 4 figures, 1 table, Published to MNRAS

We present an observational study of the luminous red nova (LRN) AT 2021biy in the nearby galaxy NGC 4631. The field of the object was routinely imaged during the pre-eruptive stage by synoptic surveys, but the transient was detected only at a few epochs from ∼231 days before maximum brightness. The LRN outburst was monitored with unprecedented cadence both photometrically and spectroscopically. AT 2021biy shows a short-duration blue peak, with a bolometric luminosity of ∼1.6 × 1041 erg s−1, followed by the longest plateau among LRNe to date, with a duration of 210 days. A late-time hump in the light curve was also observed, possibly produced by a shell-shell collision. AT 2021biy exhibits the typical spectral evolution of LRNe. Early-time spectra are characterised by a blue continuum and prominent H emission lines. Then, the continuum becomes redder, resembling that of a K-type star with a forest of metal absorption lines during the plateau phase. Finally, late-time spectra show a very red continuum (TBB ≈ 2050 K) with molecular features (e.g., TiO) resembling those of M-type stars. Spectropolarimetric analysis indicates that AT 2021biy has local dust properties similar to those of V838 Mon in the Milky Way Galaxy. Inspection of archival Hubble Space Telescope data taken on 2003 August 3 reveals a ∼20 M⊙ progenitor candidate with log (L/L⊙) = 5.0 dex and Teff = 5900 K at solar metallicity. The above luminosity and colour match those of a luminous yellow supergiant. Most likely, this source is a close binary, with a 17–24 M⊙ primary component.