The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics including as the progenitors of gammaray bursts GRBsSUP1SUP sources of highfrequency gravitational waves GWsSUP2SUP and likely production sites for heavyelement nucleosynthesis by means of rapid neutron capture the rprocessSUP3SUP. Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of longduration GRBs associated with compact object mergersSUP46SUP and contains a kilonova similar to AT2017gfo associated with the GW merger GW170817 refs. SUP712SUP. We obtained James Webb Space Telescope JWST midinfrared imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium atomic mass A  130 and a very red source emitting most of its light in the midinfrared owing to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create rprocess elements across a broad atomic mass range and play a central role in heavyelement nucleosynthesis across the Universe.
['Astrophysics - High Energy Astrophysical Phenomena', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Heavyelement production in a compact object merger observed by JWST
{'JWST detection of heavy neutron capture elements in a compact object merger': "Andrew Levan 1,2* , Benjamin P. Gompertz 3 , Om Sharan Salafia 4,5 , Mattia Bulla 6,7,8 , Eric Burns 9 , Kenta Hotokezaka 10,11 , Luca Izzo 12,13 , Gavin P. Lamb 14,15 , Daniele B. Malesani 1,16,17 , Samantha R. Oates 3 , Maria Edvige Ravasio 1,4 , Alicia Rouco Escorial 18 , Benjamin Schneider 19 , Nikhil Sarin 20,21 , Steve Schulze 21 , Nial R. Tanvir 15 , Kendall Ackley 2 , Gemma Anderson 22 , Gabriel B. Brammer 16,17 , Lise Christensen 16,17 , Vikram S. Dhillon 23,24 , Phil A. Evans 15 , Michael Fausnaugh 19,25 , Wen-fai Fong 26,27 , Andrew S. Fruchter 28 , Chris Fryer 29,30,31,32 , Johan P. U. Fynbo 16,17 , Nicola Gaspari 1 , Kasper E. Heintz 16,17 , Jens Hjorth 12 , Jamie A. Kennea 33 , Mark R. Kennedy 34,35 , Tanmoy Laskar 1,36 , Giorgos Leloudas 37 , Ilya Mandel 38,39 , Antonio Martin-Carrillo 40 , Brian D. Metzger 41,42 , Matt Nicholl 43 , Anya Nugent 26,27 , Jesse T. Palmerio 44 , Giovanna Pugliese 45 , Jillian Rastinejad 26,27 , Lauren Rhodes 46 , Andrea Rossi 47 , Stephen J. Smartt 43,46 , Heloise F. Stevance 46,48 , Aaron Tohuvavohu 49 , Alexander van der Horst 32 , Susanna D. Vergani 44 , Darach Watson 16,17 , Thomas Barclay 50 , Kornpob Bhirombhakdi 28 , Elm'e Breedt 51 , Alice A. Breeveld 52 , Alexander J. Brown 23 , Sergio Campana 4 , Ashley A. Chrimes 1 , Paolo D'Avanzo 4 , Valerio D'Elia 53,54 , Massimiliano De Pasquale 55 , Martin J. Dyer 23 , Duncan K. Galloway 38,39 , James A. Garbutt 23 , Matthew J. Green 56 , Dieter H. Hartmann 57 , P'all Jakobsson 58 , Paul Kerry 23 , Danial Langeroodi 12 , James K. Leung 59,60,39 , Stuart P. Littlefair 23 , James Munday 2,61 , Paul O'Brien 15 , Steven G. Parsons 23 , Ingrid Pelisoli 2 , Andrea Saccardi 44 , David I. Sahman 23 , Ruben Salvaterra 62 , Boris Sbarufatti 4 , Danny Steeghs 2,39 , Gianpiero Tagliaferri 4 , Christina C. Thone 63 , Antonio de Ugarte Postigo 64 , David Alexander Kann 65  \n- 1 Department of Astrophysics/IMAPP, Radboud University, 6525 AJ Nijmegen, The Netherlands.\n- 2 Department of Physics, University of Warwick, Coventry, CV4 7AL, UK. 3 Institute for Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Birmingham, B15 2TT, UK.\n- INAF - Osservatorio Astronomico di Brera, Via E. Bianchi 46, I-23807,\n- 4 Merate (LC), Italy.\n- 5 INFN - Sezione di Milano-Bicocca, Piazza della Scienza 2, I-20146, Milano (MI), Italy.\n- 6 Department of Physics and Earth Science, University of Ferrara, via Saragat 1, I-44122 Ferrara, Italy.  \n7  \n- INFN, Sezione di Ferrara, via Saragat 1, I-44122 Ferrara, Italy. 8 INAF, Osservatorio Astronomico d'Abruzzo, via Mentore Maggini snc,\n- 64100 Teramo, Italy.\n- 9 Department of Physics & Astronomy, Louisiana State University, Baton Rouge, LA 70803, USA.\n- 10 Research Center for the Early Universe, Graduate School of Science, The University of Tokyo, Bunkyo, Tokyo 113-0033, Japan.\n- 11 Kavli IPMU (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan.\n- 12 DARK, Niels Bohr Institute, University of Copenhagen, Jagtvej 128, 2200 Copenhagen N, Denmark.\n- 13 INAF-Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, 80131, Napoli, Italy.\n- 14\n- Astrophysics Research Institute, Liverpool John Moores University, IC2 Liverpool Science Park, Liverpool, L3 5RF, Liverpool, UK. 15 School of Physics & Astronomy, University of Leicester, University Road, Leicester, LE1 7RH, UK.  \n16  \n- Cosmic Dawn Center (DAWN), Denmark.\n- 17 Niels Bohr Institute, University of Copenhagen, Jagtvej 128, 2200 Copenhagen N, Denmark.\n- 18 European Space Agency (ESA), European Space Astronomy Centre (ESAC), Camino Bajo del Castillo s/n, 28692 Villanueva de la Ca˜nada, Madrid, Spain.  \n19 Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA 02139, USA.  \n20 Nordita, Stockholm University and KTH Royal Institute of Technology Hannes Alfv'ens vag 12, SE-106 91 Stockholm, Sweden.  \n- 21 The Oskar Klein Centre, Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden.\n- 22 International Centre for Radio Astronomy Research, Curtin University, GPO Box U1987, Perth, WA 6845, Australia.\n- 23 Department of Physics and Astronomy, University of Sheffield, Sheffield, S3 7RH, United Kingdom.\n- 24 Instituto de Astrof'ısica de Canarias, E-38205 La Laguna, Tenerife, Spain.\n- 25 Department of Physics & Astronomy, Texas Tech University, Lubbock TX, 79410-1051, USA.\n- 26 Center for Interdisciplinary Exploration and Research in Astrophysics, Northwestern University, 1800 Sherman Ave., Evanston, 60208, IL, USA. 27 Department of Physics and Astronomy, Northwestern University, 2145\n- Sheridan Road, Evanston, 60208-3112, IL, USA.\n- 28 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218.\n- 29 Center for Theoretical Astrophysics, Los Alamos National Laboratory, Los Alamos, NM 87545.\n- 30 Department of Astronomy, The University of Arizona, Tucson, AZ 85721.\n- 31 Department of Physics and Astronomy, The University of New Mexico, Albuquerque, NM 87131.\n- 32 Department of Physics, The George Washington University, Washington, DC 20052.\n- 33 Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Lab, University Park, PA 16802, USA. 34 School of Physics, Kane Building, University College Cork, Cork, Ireland.\n- 35 Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, The University of Manchester, M13 9PL, UK.\n- 36\n- Department of Physics & Astronomy, University of Utah, Salt Lake City, UT 84112, USA.\n- 37 DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark. School of Physics and Astronomy, Monash University, Clayton,\n- 38 Victoria 3800, Australia.  \n- 61 Isaac Newton Group of Telescopes, Apartado de Correos 368, E-38700 Santa Cruz de La Palma, Spain. 62 INAF/IASF-MI, Via Alfonso Corti 12, I-20133, Milano, Italy. 63 Astronomical Institute of the Czech Academy of Sciences, Friˇcova 298, 251 65 Ondˇrejov, Czech Republic. 64 Artemis, Observatoire de la Cˆote d'Azur, Universit'e Cˆote d'Azur, Boulevard de l'Observatoire, F-06304 Nice, France. 65 Hessian Research Cluster ELEMENTS, Giersch Science Center,\n- Max-von-Laue-Straße 12, Goethe University Frankfurt, Campus Riedberg, D-60438 Frankfurt am Main, Germany.  \n*Corresponding author(s). E-mail(s): a.levan@astro.ru.nl;", 'Abstract': "The mergers of binary compact objects such as neutron stars and black holes are of central interest to several areas of astrophysics, including as the progenitors of gamma-ray bursts (GRBs), sources of high-frequency gravitational waves and likely production sites for heavy element nucleosynthesis via rapid neutron capture (the r -process). These heavy elements include some of great geophysical, biological and cultural importance, such as thorium, iodine and gold. Here we present observations of the exceptionally bright gamma-ray burst GRB 230307A. We show that GRB 230307A belongs to the class of long-duration gamma-ray bursts associated with compact object mergers, and contains a kilonova similar to AT2017gfo, associated with the gravitational-wave merger GW170817. We obtained James Webb Space Telescope mid-infrared (mid-IR) imaging and spectroscopy 29 and 61 days after the burst. The spectroscopy shows an emission line at 2.15 microns which we interpret as tellurium (atomic mass A=130), and a very red source, emitting most of its light in the mid-IR due to the production of lanthanides. These observations demonstrate that nucleosynthesis in GRBs can create r -process elements across a broad atomic mass range and play a central role in heavy element nucleosynthesis across the Universe.  \nKeywords: Gamma-ray burst, Nucleosynthesis, Neutron star merger  \nGRB 230307A was first detected by the Fermi Gamma-ray Burst Monitor (GBM) at 15:44:06 UT on 7 Mar 2023 [1]. It was an exceptionally bright gamma-ray burst with a duration of T 90 ∼ 35 s and a prompt fluence of (2.951 ± 0.004) × 10 -3 erg cm -2 in the 10-1000 keV band [2] (Figure 1). These properties place this event at the peak of the distribution of the class of 'long-soft' GRBs. The measured fluence makes it the second-brightest GRB ever detected [3].  \nIn addition to Fermi , the burst was also detected by several other high-energy instruments (see Methods). The multiple detections enabled source triangulation by  \nthe InterPlanetary Network (IPN). The Neil Gehrels Swift Observatory (hereafter Swift ) tiled the initial ∼ 2 sq. deg IPN region [4] which revealed one candidate X-ray afterglow [5] consistent with the final 8 sq. arcmin IPN localization [6].  \nWe obtained optical observations of the field of GRB 230307A with the ULTRACAM instrument, mounted on the 3.5m New Technology Telescope (NTT) at La Silla, Chile. Visual inspection of the images compared to those obtained with the Legacy Survey [7] revealed a new source coincident with the Swift X-ray source, and we identified it as the optical afterglow of GRB 230307A [8]. This identification was subsequently confirmed via imaging by several additional observatories [9-11]. The location was also serendipitously imaged by the Transiting Exoplanet Survey Satellite (TESS) from 3 days before to 3 days after the GRB [12].  \nFollowing the precise localisation, we obtained extensive follow-up observations, in the optical and near infrared with the Gemini South Telescope and the Very Large Telescope (VLT); in the X-ray with the Swift /XRT and the Chandra X-ray observatory; and in the radio with the Australia Telescope Compact Array (ATCA) and MeerKAT in South Africa. This campaign included spectroscopy with the VLT Xshooter instrument, as well as the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph. The latter provides redshift information for many galaxies in the field, including, in particular, a bright spiral galaxy at z = 0 . 0646 ± 0 . 0001 offset 30.2 arcseconds (38.9 kiloparsec in projection) from the burst position (Figure 2). Of the optically detected galaxies in the field, this is the one with the lowest probability of being an unrelated chance alignment, and hence is most likely to be the host (see also [13]).  \nOur ground-based campaign included imaging from 1.4 to 41 days after the burst (see Supplementary Information Table 1). At 11 days, infrared observations demonstrated a transition from an early blue spectral slope, to a much redder one with i -K > 2 . 9(AB). This extremely red colour appeared similar to the expectations for a kilonova, powered by the decay of unstable isotopes of heavy elements synthesised by rapid neutron capture within the ejecta produced during the merger of a neutron star and another compact object [14-16]. Based on this detection, we requested James Webb Space Telescope (JWST) observations (GO 4434, 4445, PI Levan), which were initiated on 5 April 2023. At the first epoch (+28.9 days after GRB), we took 6colour observations with the Near Infrared Camera (NIRCam) (Figure 2), as well as a spectrum with the Near Infrared Spectrograph (NIRSpec) covering 0 . 5 -5 . 5 microns (Figure 3).  \nThe NIRCam observations reveal an extremely red source that is only weakly detected in the bluer bands, where F150W(AB) = 28 . 11 ± 0 . 12 mag, but rises sharply through the mid-IR to F444W(AB) = 24 . 4 ± 0 . 01 mag. The NIRSpec observations also exhibit this steep rise. A faint galaxy is detected in these data at z = 3 . 87, offset approximately 0.3 arcseconds from the burst position. However, the burst's properties are inconsistent with an origin at this redshift, in particular because the implied isotropic equivalent energy release would exceed all known GRBs by an order of magnitude or more (see Supplementary Information). A second epoch of JWST observations was obtained approximately 61 days after the burst. These observations showed that the source had faded by 2.4 magnitudes in F444W, demonstrating a rapid  \ndecay expected in the kilonova scenario and effectively ruling out alternatives (see Supplementary Information).  \nSome of the burst properties are remarkably similar to those of the bright GRB 211211A, which was also shown to be accompanied by a kilonova [17-19]. In particular, the prompt emission consists of a hard pulse lasting for ∼ 19 seconds, followed by much softer emission (Figure 1). The prompt emission spectrum is well modelled by a double broken power-law with two spectral breaks moving rapidly through the gammaray band (see Methods), suggesting a synchrotron origin of the emission [20]. The X-ray afterglow is exceptionally faint, much fainter than most bursts when scaled by the prompt GRB fluence (see Figure 1 and Supplementary Information). The development of the optical and IR counterpart is also similar to GRB 211211A, with an early blue colour and a subsequent transition to red on a timescale of a few days. In Figure 4, we plot the evolution of the counterpart compared with the kilonova AT2017gfo [2128], identified in association with the gravitational-wave detected binary neutron star merger, GW170817 [29]. AT2017gfo is the most rapidly evolving thermal transient ever observed; much more rapid than supernovae or even, for example, fast blue optical transients [30]. The counterpart of GRB 230307A appears to show near identical decline rates to AT2017gfo both at early times in the optical and IR, and later in the mid-IR [31]. These similarities are confirmed by a joint fit of afterglow and kilonova models to our multi-wavelength data (see Supplementary Information). Our JWST observations rule out supernovae: for any redshift z < 1, a supernova would need to be more than 100 times fainter than the canonical GRB-supernova, SN 1998bw, to be compatible with our observations. Therefore, we conclude that GRB 230307A is a long-duration GRB formed from a compact object merger. This falls into a class that includes GRB 211211A [17-19, 32], GRB 060614 [33], GRB 111005A [34] and GRB 191019A [35], among others.  \nThe JWST observations provide a view of a kilonova in the mid-IR with high spatial resolution and sensitivity. On timescales of ∼ 30 days, it is apparent that the kilonova emits almost all of its light in the mid-infrared, beyond the limits of sensitive groundbased observations (effectively limited to below 2.5 microns). This is consistent with previous model predictions [36], but has not previously been observationally confirmed. Late-time studies of such emission in the nebular phase must therefore be conducted in the mid-IR. Strikingly, despite its powerful and long-lived prompt emission that stands in stark contrast to GRB 170817A, the GRB 230307A kilonova is remarkably similar to AT 2017gfo. This was also the case for the long-lived merger GRB 211211A [1719, 32], and suggests the kilonova signal, particularly the red component, is relatively insensitive to the GRB.  \nOur NIRSpec spectrum shows a broad emission feature with a central wavelength of 2.15 microns, visible in both epochs of JWST spectroscopy. At longer wavelengths, the spectrum displays a slowly rising continuum up to 4.5 microns followed by either an additional feature or change of spectral slope. The colours of the counterpart at this time can be readily explained by kilonova models (see Supplementary Information Section 6).  \nA similar emission-like feature is also visible in the later epochs of X-shooter observations of AT2017gfo [23], measured at 2.1 microns by [37]. This further strengthens both the kilonova interpretation and the redshift measurement of GRB 230307A (Figure 3). We interpret this feature as arising from the forbidden [Te III] transition between the ground level and the first fine structure level of tellurium, with an experimentally-determined wavelength of 2.1044 microns [38]. The presence of tellurium is plausible, as it lies at the second peak in the r -process abundance pattern, which occurs at atomic masses around A ≈ 130 [39]. It should therefore be abundantly produced in kilonovae, as is seen in hydrodynamical simulations of binary neutron star mergers with nucleosynthetic compositions similar to those favoured for AT2017gfo [40]. Furthermore, the typical ionization state of Te in kilonova ejecta is expected to be Te III at this epoch because of the efficient radioactive ionization [41]. Tellurium has recently been suggested as the origin of the same feature in the spectrum of AT2017gfo [42]. A previous study [37] also identified this tellurium transition and noted that the observed feature is most likely two blended emission lines. However, alternative transitions from heavy r -process elements have been considered for this feature [e.g. Ce III; 37]. Tellurium can also be produced via the slower capture of neutrons in the s -process. Indeed, this line is also seen in planetary nebulae [43]. The detection of [Te III] 2 . 1 µ m extends on the earlier detection of strontium, a first r -process peak element, in the early time photospheric emission of AT2017gfo [44]. The mass of Te III estimated from the observed line flux is ∼ 10 -3 M ⊙ (see Supplementary Information 6.2). Although weaker, we also note that the spectral feature visible at 4.5 microns is approximately consistent with the expected location of the first peak element selenium and the near-third peak element tungsten [45].  \nDetailed spectral fitting at late epochs is challenging because of the breakdown of the assumptions regarding local thermodynamic equilibrium (LTE) which are used to predict kilonova spectra at earlier ages, as well as fundamental uncertainties in the atomic physics and electron transitions in the highly complex electron shells of r -process elements. However, these observations provide a calibration sample for informing future models. The red continuum emission indicates a large opacity in the mid-IR at low temperatures, e.g., ∼ 10 cm 2 g -1 at ∼ 700 K. Since the Planck mean opacity under LTE is expected to be less than 1 cm 2 g -1 at ≲ 1000 K [46], this large opacity may suggest that lanthanides are abundant in the ejecta. Indeed, it has been shown that non-LTE effects can increase the lanthanide opacity in mid-IR at late times [47]. Therefore, systematic studies of non-LTE opacity are necessary to answer the question whether lanthanides are the origin of the red emission at this epoch.  \nA fit to our combined MUSE + JWST data for the host galaxy (Supplementary Information) suggests a relatively low mass ( ∼ 2 . 5 × 10 9 M ⊙ ) dominated by an older ∼ 10 Gyr population, but with a second more recent burst of star formation. These properties are entirely consistent with the population of short GRB host galaxies [48, 49]. The host normalised offset of the burst from the host galaxy places it in the top 10% of those seen for short GRBs [48, 50]. The offset could readily be achieved by a binary with a velocity of a few hundred km s -1 and a merger time of hundreds of millions of years. Alternatively, in the second epoch of JWST observations, a faint source is detected in the F150W images at the location of the transient. This may  \nrepresent continued emission from the transient. However, its absolute magnitude of M F 150 W ∼ -8 . 5 is comparable to the absolute magnitude of globular clusters in which dynamical interactions could be at play to create merging systems at enhanced rates [51]. Future observations should readily distinguish between a fading afterglow or underlying cluster.  \nIt is striking that GRB 230307A is an extremely bright GRB, with only the exceptional GRB 221009A being brighter [52]. Of the ten most fluent Fermi /GBM GRBs, two are now associated with kilonovae (230307A and 211211A), three are associated with supernovae, and the nature of the remaining five appears unclear (see Table 4). For bright GRBs, there may be a significant contribution from mergers. Indeed, such a conclusion can also be reached by considering the energetics. Both GRB 230307A and GRB 211211A have isotropic equivalent energies of E iso > 10 51 erg. The majority of local GRBs for which the connection between GRBs and core-collapse supernovae has been established are much less energetic (typically E iso > 10 49 -50 erg) and it has been suggested that they represent a separate population powered via different processes [53]. For more energetic bursts in the local Universe (where supernovae can still readily be discovered) the fraction of long GRBs with and without supernovae appear similar (see Supplementary Information). If a substantial number of long GRBs are associated with compact object mergers, they provide an essential complement to gravitationalwave (GW) detections. Firstly, joint GW-GRB detections, including long GRBs, can push the effective horizons of GW detectors to greater distances and provide much smaller localisations. In the case of GRB 230307A, the distance of 300 Mpc could have provided a robust detection in gravitational waves for the relevant O4 sensitivity [17, 54]. Secondly, long GRBs can be detected without GW detectors, providing a valuable route to enhancing the number of kilonova detections. Thirdly, JWST can detect kilonova emission at redshifts substantially beyond the horizons of the current generation of GW detectors, enabling the study of kilonovae across a greater volume of the Universe.  \nThe duration of the prompt γ -ray emission in these mergers remains a challenge to explain. In particular, the natural timescales for emission in compact object mergers are much shorter than the measured duration of GRB 230307A. Previously suggested models that may also explain GRB 230307A include magnetars [55, 56], black hole - neutron star mergers [57, 58], or even neutron star - white dwarf systems [19, 59]. Recent results have also shown that the jet timescale does not directly track the accretion timescale in compact object mergers, and that long GRBs may be created from very short lived engines [60], and hence from binary neutron star mergers without magnetars.  \nThere is evidence that the kilonova in GRB 230307A produced elements across a wide range of atomic mass. The detection of second peak elements in the spectrum of a kilonova demonstrates that nuclei with atomic masses around A ∼ 130 are being created in the mergers of compact objects. Many second peak elements have important biological roles. For example iodine is essential for mammals and may have been used by the single cell Last Universal Common Ancestor [61]. The creation of these elements in compact object mergers, which can have long delay times, may have important consequences for the time at which certain evolutionary channels become plausible.  \nFig. 1 The high energy properties of GRB 230307A. The left panel shows the light curve of the GRB at 64 ms time resolution with the Fermi /GBM. The shaded region indicates the region where saturation may be an issue. The burst begins very hard, with the count rate dominated by photons in the hardest (100-900 keV) band, but rapidly softens, with the count rate in the hard band being progressively overtaken by softer bands (e.g. 8-25 keV and 25-100 keV) beyond ∼ 20s. This strong hard-to-soft evolution is reminiscent of GRB 211211A [20] and is caused by the motion of two spectral breaks through the γ -ray regime (see Methods). The right panel shows the X-ray light curves of GRBs from the Swift X-ray telescope. These have been divided by the prompt fluence of the burst, which broadly scales with the X-ray light curve luminosity, resulting in a modest spread of afterglows. The greyscale background represents the ensemble of long GRBs. GRB 230307A is an extreme outlier of the > 1000 Swift -GRBs, with an extremely faint afterglow for the brightness of its prompt emission. Other merger GRBs from long bursts occupy a similar region of parameter space. This suggests the prompt to afterglow fluence could be a valuable tool in distinguishing long GRBs from mergers and those from supernovae.  \n<!-- image -->  \nFig. 2 JWST images of GRB 230307A at 28.5 days post burst. The upper panel shows the wide field image combining the F115W, F150W and F444W images. The putative host is the bright face-on spiral galaxy, while the afterglow appears at a 30-arcsecond offset, within the white box. The lower panel shows cut-outs of the NIRCAM data around the GRB afterglow location. The source is faint and barely detected in the bluer bands but very bright and well detected in the red. In the red bands, a faint galaxy is present northeast of the transient position. This galaxy has a redshift of z = 3 . 87, but we consider it to be a background object unrelated to the GRB (see Supplementary Information).  \n<!-- image -->  \nFig. 3 JWST/NIRSpec spectroscopy of the counterpart of GRB 230307A taken on 5 April 2023. The top panel shows the 2-D spectrum rectified to a common wavelength scale. The transient is well detected beyond 2 microns but not short ward, indicative of an extremely red source. Emission lines from the nearby galaxy at z = 3 . 87 can also be seen offset from the afterglow trace. The lower panel shows the 1D extraction of the spectrum in comparison with the latest (10-day) AT2017gfo epoch and different kilonova models. A clear emission feature can be seen at ∼ 2.15 microns at both 29 and 61 days. This feature is consistent with the expected location of [Te III], while redder features are compatible with lines from [Se III] and [W III]. This line is also clearly visible in the late time spectrum of AT2017gfo [37, 42]  \n<!-- image -->  \nFig. 4 A comparison of the counterpart of GRB 230307A with AT2017gfo (the kilonova associated with GW170817) at z = 0 . 065. Beyond ∼ 2 days, the kilonova dominates the counterpart (see Supplementary Information). The decay rates in both the optical and IR are very similar to those in AT2017gfo. These are too rapid for any plausible afterglow model (e.g. as a power-law, they decay faster than t -3 . 5 over a prolonged period). There is also good agreement in the late time slope between the measurements made at 4.4 microns with JWST and 4.5 microns for AT2017gfo with Spitzer [31].  \n<!-- image -->", '2 Observations': 'Below we outline the observational data that were used in this paper. Magnitudes are given in the AB system unless stated otherwise. We utilize cosmology resulting from the Planck observations [62]. All uncertainties are given at the 1 σ level unless explicitly stated.', '2.1 γ -ray observations': 'GRB 230307A was first detected by Fermi /GBM and GECAM at 15:44:06 UT on 7 Mar 2023 [1]. It had a duration of T 90 ∼ 35s and an exceptionally bright prompt fluence of (2.951 ± 0.004) × 10 -3 erg cm -2 [2]. The burst fell outside of the coded field-of-view of the Swift Burst Alert Telescope (BAT), and so did not receive a subdegree localisation despite a strong detection. However, detections by Swift , GECAM [63], STIX on the Solar Orbiter [64], AGILE [65], ASTROSAT [66], GRBalpha [67], VZLUSAT [68], Konus-WIND [69] and ASO-HXI [70] enabled an enhanced position via the InterPlanetary Network to increasingly precise localisations of 1.948 deg 2 [71], 30 arcmin 2 [72], and ultimately to 8 arcmin 2 [6]. This was sufficiently small to enable tiling with Swift and ground-based telescopes.', '2.1.1 Fermi /GBM data analysis': 'In Figure 1, we plot the light curve of GRB 230307A as seen by the Fermi /GBM in several bands, built by selecting Time Tagged Event (TTE) data, binned with a time resolution of 64 ms. The highlighted time interval of 3-7 s after trigger are affected by data loss due to the bandwidth limit for TTE data [73].  \nFor the spectral analysis, we made use of the CSPEC data, which have 1024 ms time resolution. Data files were obtained from the online archive 1 . Following the suggestion reported by the Fermi Collaboration [73], we analysed the data detected by NaI 10 and BGO 1, which had a source viewing angle less than 60 · , and excluded the time intervals affected by pulse pile-up issues (from 2.5 s to 7.5 s). The data extraction was performed with the public software gtburst , while data were analysed with Xspec . The background, whose time intervals have been selected before and after the source, was modelled with a polynomial function whose order is automatically found by gtburst and manually checked. In the fitting procedure, we used intercalibration factors among the detectors, scaled to the only NaI analysed and free to vary within 30%. We used the PG-Statistic, valid for Poisson data with a Gaussian background. The best-fit parameters and their uncertainties were estimated through a Markov Chain Monte Carlo (MCMC) approach. We selected the time intervals before and after the excluded period of 2.5-7.5 s due to instrumental effects. In particular, we extracted 2 time intervals from 0 to 2.5 s (1.25 s each) and 14 time intervals from 7.5 s to 40.5 (bin width of 2 s, except the last two with integration of 5 s to increase the signal-to-noise ratio), for a total of 16 time intervals. We fitted the corresponding  \nspectra with the two smoothly broken power-law (2SBPL) function [74, 75], which has been shown to successfully model the synchrotron-like spectral shape of bright long GRBs, including the merger-driven GRB 211211A [20].  \nFrom our spectral analysis we found that all spectra up to ∼ 20 s are well modelled by the 2SBPL function, namely they are described by the presence of two spectral breaks inside the GBM band (8 keV-40 MeV). In particular, in the time intervals between 7.5 and 19.5 s, the low-energy break E break is coherently decreasing from 304 . 3 +5 . 2 -2 . 6 keV down to 52 . 1 +4 . 3 -5 . 1 keV, and the typical νF ν peak energy E peak is also becoming softer, moving from ∼ 1 MeV to 450 keV. The spectral indices of the two power-laws below and above the low-energy break are distributed around the values of -0.82 and -1.72, which are similar to the predictions for synchrotron emission in marginally fast-cooling regime (i.e. -2 / 3 and -3 / 2). This is consistent with what has been found in GRB 211211A [20]. We notice, however, that in all spectra the highenergy power-law above E peak is characterised by a much softer index (with a mean value of -4 . 10 ± 0 . 24) with respect to the value of ∼ -2.5 typically found in Fermi GRBs. This suggests that the spectral data might require a cut-off at high energy, although further investigations are needed to support this. From 19.5 s until 40.5 s (the last time interval analysed), all the break energies are found to be below 20 keV, close to the GBM low energy threshold. In the same time intervals, the peak energy E peak decreases from 682 . 4 +3 . 2 -6 . 1 to 123 . 1 +5 . 4 -4 . 9 keV, and the index of the power-law below the peak energy is fully consistent (mean value of -1 . 45 ± 0 . 06) with the synchrotron predicted value of -1.5.', '2.2.1 NTT - Afterglow discovery': 'Following the refinement of the IPN error box to an area of 30 arcmin 2 [72], we obtained observations of the field of GRB 230307A with the ULTRACAM instrument [76], mounted on the 3.5m New Technology Telescope (NTT) at La Silla, Chile. The instrument obtains images in 3 simultaneous bands, and is optimised for short exposure, low dead-time observations [76]. We obtained 10 × 20 s exposures in two pointings in each of the Super SDSS u , g and r , bands (where the Super SDSS bands match the wavelength range of the traditional SDSS filters but with a higher throughput; Dhillon et al. 77). The observations began at 01:53:21 UT on 2023-03-09, approximately 34 hr after the GRB. The images were reduced via the HIPERCAM pipeline [77] using bias and flat frames taken on the same night. Visual inspection of the images compared to those obtained with the Legacy Survey [7] revealed a new source coincident with an X-ray source identified via Swift /XRT observations [5], and we identified it as the likely optical afterglow of GRB 230307A [8]. The best available optical position of this source (ultimately measured from our JWST observations, see below) is RA(J2000) = 04:03:26.02, Dec(J2000) = -75:22:42.76, with an uncertainty of 0.05 arcseconds in each axis. The IPN error box and the footprint of the ULTRACAM observations are shown in Figure 5.  \nThis identification was subsequently confirmed via observations from a number of additional observatories, including [9-11, 78-80]. We acquired two further epochs of  \nobservations with ULTRACAM on the following nights with 10 × 20s exposures in the Super SDSS u , g and i , bands. Aperture photometry of the source is reported in Table 1, and is reported relative to the Legacy survey for the gri bands, and to SkyMapper for the u -band.', '2.2.2 TESS': 'The prompt and afterglow emission of GRB 230307A was detected by TESS, which observed the field continuously from 3 days before the Fermi trigger to 3 days after at a cadence of 200 s [12]. A reference image was subtracted from the observations to obtain GRB-only flux over this period. The measured flux in the broad TESS filter (600nm - 1000nm) is corrected for Galactic extinction and converted to the I c band assuming a power-law spectrum with F ∝ ν -0 . 8 . We then bin the light curve logarithmically, taking the mean flux of the observations in each bin and converting to AB magnitudes. A systematic error of 0.1 magnitudes was added in quadrature to the measured statistical errors to account for the uncertainties in the data processing. These data are presented in Table 1.', '2.2.3 Swift /UVOT': 'The Swift Ultra-violet/Optical Telescope [UVOT; 81] began observing the field of GRB 230307A ∼ 84 . 6 ks after the Fermi /GBM trigger [1]. The source counts were extracted using a source region of 5 arcsec radius. Background counts were extracted using a circular region of 20 arcsec radius located in a source-free part of the sky. The count rates were obtained from the image lists using the Swift tool uvotsource . A faint catalogued unrelated source also falls within the 5 arcsec radius, this will affect the photometry, particularly at late times. We, therefore, requested a deep template image in white in order to estimate the level of contamination. We extracted the count rate in the template image using the same 5 arcsec radii aperture. This was subtracted from the source count rates to obtain the afterglow count rates. The afterglow count rates were converted to magnitudes using the UVOT photometric zero points [82, 83].', '2.2.4 Gemini': 'We obtained three epochs of K-band observations using the FLAMINGOS-2 instrument on the Gemini-South telescope. These observations were reduced through the DRAGONS pipeline to produce dark and sky-subtracted and flat-fielded images [84]. At the location of the optical counterpart to GRB 230307A, we identify a relatively bright K-band source in the first and second epochs, with only a upper limit in epoch 3. We report our photometry, performed relative to secondary standards in the VISTA hemisphere survey [85], in Table 1.', '2.2.5 VLT imaging': "We carried out observations of the GRB 230307A field with the 8.2-m VLT telescopes located in Cerro Paranal, Chile. The observations were obtained with the FORS2 camera (mounted on the Unit Telescope 1, UT1, ANTU) in B , R , I and z bands at  \nmultiple epochs, and with the HAWK-I instrument (mounted on the Unit Telescope 4, UT4, Yepun) in the K band at one epoch. All images were reduced using the standard ESO (European Southern Observatory) Reflex pipeline [86]. The source was detected in the FORS2 z -band image at ∼ 6.4 days after the Fermi /GBM detection. A single r ' -band observation of the GRB 230307A was also executed with the 2.6m VLT Survey Telescope (VST) after 2.37 days from the GRB discovery. In later observations the source was not detected (see bottom right panels of Fig. 5) and the upper limit values at 3 σ level are reported in Table 1.", '2.2.6 VLT spectroscopy': 'To attempt to measure the redshift of GRB 230307A and the nearby candidate host galaxies, we obtained spectroscopy with the VLT utilising both the X-shooter and MUSE instruments, mounted respectively on the Unit Telescope 2 (UT2, Kueyen) and on UT4 (Yepun).  \nX-shooter spectroscopy, covering the wavelength range 3000-22000 ˚ A was undertaken on 2023-03-15. Observations were taken at a fixed position angle, with the slit centred on a nearby bright star. X-shooter data have been reduced with standard esorex recipes. Given that only two of the four nod exposures were covering the GRB position, resulting in a total exposure time of 2400s on-source, we have reduced each single exposure using the stare mode data reduction. Then, we have stacked the two 2D frames covering the GRB position using dedicated post-processing tools developed in a python framework [87]. We further obtained observations with the MUSE integral field unit on 2023-03-23. The MUSE observations cover multiple galaxies within the field, as well as the GRB position, and cover the wavelength range 4750-9350 ˚ A. MUSE data have been reduced using standard esorex recipes embedded within a single python script that performs the entire data reduction procedure. Later, the resulting datacube has been corrected for sky emission residuals using ZAP [88]. The MUSE observations reveal the redshifts for a large number of galaxies in the field, including a prominent spiral G1 at z = 0 . 0646 [see also 80] and a group of galaxies, G2, G3 and G4 at z = 0 . 263.', '2.3 X-ray afterglow': "Swift began tiled observations of the IPN localisation region with its X-ray Telescope [XRT; 89] at 12:56:42 on 8 Mar 2023 2 [90]. XRT made the first reported detection of the afterglow (initially identified as 'Source 2') with a count rate of 0 . 019 ± 0 . 004 cts -1 [91] and later confirmed it to be fading with a temporal power-law index of 1 . 1 +0 . 6 -0 . 5 [92]. XRT data were downloaded from the UK Swift Science Data Centre [UKSSDC; 93, 94].  \nWe further obtained observations with the Chandra X-ray observatory (programme ID 402458: PI Fong/Gompertz). A total of 50.26 ks (49.67 ks of effective exposure) of data were obtained in three visits between 31 March 2023 and 2 April 2023. The source was placed at the default aim point on the S3 chip of the ACIS detector. At the location of the optical and X-ray afterglow of GRB 230307A, we detect a total of  \n12 counts, with an expected background of ∼ 1, corresponding to a detection of the afterglow at > 5 σ based on the photon statistics of [95]. To obtain fluxes, we performed a joint spectral fit of the Chandra and Swift /XRT data. The best fitting spectrum, adopting uniform priors on all parameters, is a power law with a photon index of Γ = 2 . 50 +0 . 30 -0 . 29 when fitting with a Galactic N H = 1 . 26 × 10 21 cm -2 [96] and zero intrinsic absorption (neither XRT nor Chandra spectra have sufficient signal to noise to constrain any intrinsic absorption component). The resultant flux in the 0.3 - 10 keV band is F X (1 . 7 d) = 4 . 91 +0 . 89 -0 . 79 × 10 -13 erg cm -2 s -1 during the XRT observation and F X (24 . 8 d) = 1 . 19 +0 . 87 -0 . 62 × 10 -14 erg cm -2 s -1 during the Chandra observation. Due to the low count number, the Chandra flux posterior support extends to considerably below the reported median, with the 5th percentile being as low as F X , 5th = 3 × 10 -15 erg cm -2 s -1 . If a uniform-in-the-logarithm prior on the flux were adopted, this would extend to even lower values. Chandra and XRT fluxes are converted to 1 keV flux densities using the best fit spectrum (Table 2).", '2.4 ATCA': 'Following the identification of the optical afterglow [97], we requested target of opportunity observations of GRB 230307A (proposal identification CX529) with the Australia Telescope Compact Array (ATCA) to search for a radio counterpart. These data were processed using Miriad [98], which is the native reduction software package for ATCA data using standard techniques. Flux and bandpass calibration were performed using PKS 1934-638, with phase calibration using interleaved observations of 0454-810.  \nThe first observation took place on 2023-03-12 at 4.46 d post-burst, which was conducted using the 4 cm dual receiver with frequencies centered at 5.5 and 9 GHz, each with a 2 GHz bandwidth. The array was in the 750C configuration 3 with a maximum baseline of 6 km. A radio source was detected at the position of the optical afterglow at 9 GHz with a flux density of 92 ± 22 µ Jy but went undetected at 5.5 GHz (3 σ upper limit of 84 µ Jy). A further two follow-up observations were also obtained swapping between the 4 cm and 15 mm dual receivers (the latter with central frequencies of 16.7 and 21.2 GHz, each with a 2 GHz bandwidth). During our second epoch at 10.66 d we detected the radio counterpart again, having become detectable at 5.5 GHz with marginal fading at 9 GHz. By the third epoch, the radio afterglow had faded below detectability. We did not detect the radio transient at 16.7 or 21.2 GHz in either epoch. All ATCA flux densities are listed in Table 2.', '2.5 MeerKAT': "We were awarded time to observe the position of GRB 230307A with the MeerKAT radio telescope via a successful Director's Discretionary Time proposal (PI: Rhodes, DDT-20230313-LR-01). The MeerKAT radio telescope is a 64-dish interferometer based in the Karoo Desert, Northern Cape, South Africa [99]. Each dish is 12 m in diameter and the longest baseline is ∼ 8 km allowing for an angular resolution of  \n∼ 7 arccsec and a field of view of 1 deg 2 . The observations we were awarded were made at both L and S-band.  \nGRB 230307A was observed over three separate epochs between seven and 41 days post-burst. The first two observations were made at both L and S4-band (the highest frequency of the five S-band sub-bands), centred at 1.28 and 3.06 GHz with a bandwidth of 0.856 and 0.875 GHz, respectively. Each observation spent two hours at L-band and 20 minutes at S4-band. The final observation was made only at S4-band with one hour on target. Please see the paper by MPIfR for further details on the new MeerKAT S-band receiver.  \nEach observation was processed using oxkat , a series of semi-automated Python scripts designed specifically to deal with MeerKAT imaging data [100]. The scripts average the data and perform flagging on the calibrators from which delay, bandpass and gain corrections are calculated and then applied to the target. The sources J04086545 and J0252-7104 were used at the flux and complex gain calibrator, respectively. Flagging and imaging of the target field are performed. We also perform a single round of phase-only self-calibration. We do not detect a radio counterpart in any epoch in either band. The rms noise in the field was measured using an empty region of the sky and used to calculate 3 σ upper limits which are given in Table 2.", '2.6 JWST observations': 'We obtained two epochs of observations of the location of GRB 230307A with JWST. The first on 5 April 2023, with observations beginning at 00:16 UT (MJD=60039.01), 28.4 days after the burst (under programme GO 4434, PI Levan), and the second on 8 May 2023, 61.5 days after the burst (programme 4445, PI Levan). The observations were at a post-peak epoch because the source was not in the JWST field of regard at the time of the burst and only entered it on 2 April 2023.  \nAt the first epoch, we obtained observations in the F070W, F115W, F150W, F277W, F356W and F444W filters of NIRCam [101], as well as a prism spectrum with NIRSpec [102]. In the second epoch we obtained NIRCam observations in F115W, F150W, F277W and F444W and a further NIRSpec prism observation. However, in the second epoch the prism observation is contaminated by light from the diffraction spike of a nearby star and is of limited use, in particular at the blue end of the spectrum. We therefore use only light redwad of 1.8 microns. However, even here we should be cautious in interpreting the overall spectral shape. However, the feature at 2.15 microns is visible in both the 29 and 61 day spectra.  \nWe reprocessed and re-drizzled the NIRCam data products to remove 1 /f striping and aid point spread function recovery, with the final images having plate scales of 0.02 arcsec/pixel (blue channel) and 0.04 arcsec/pixel (red channel).  \nIn the NIRCam imaging we detect a source at the location of the optical counterpart of GRB 230307A. This source is weakly detected in all three bluer filters (F070W, F115W and F150W), but is at high signal-to-noise ratio in the redder channels (see Figure 2). The source is compact. We also identify a second source offset (H1) approximately 0.3 arcseconds from the burst location. This source is also weakly, or non-detected in the bluer bands, and is brightest in the F277W filter.  \nBecause of the proximity of the nearby star and a contribution from diffraction spikes close to the afterglow position we model point spread functions for the appropriate bands using WebbPSF [103], and then scale and subtract these from star position. Photometry is measured in small (0.05 arcsec (blue) and 0.1 arcsec (red)) apertures and then corrected using tabulated encircled energy corrections. In addition to the direct photometry of the NIRCam images we also report a a K-band point based on folding the NIRSpec spectrum (see below), through a 2MASS, Ks filter. This both provides a better broadband SED and a direct comparison with ground based K-band observations. Details of photometric measurements are shown in Table 1.  \nFor NIRSpec, we utilise the available archive-processed level 3 two-dimensional spectrum (Figure 3). In this spectrum we clearly identify the trace of the optical counterpart, which appears effectively undetected until 2 microns and then rises rapidly. We also identify two likely emission lines which are offset from the burst position. These are consistent with the identification with H α and [O iii ] (4959/5007) at a redshift of z = 3 . 87. Both of these lines lie within the F277W filter in NIRCam and support the identification of the nearby source as the origin of these lines.  \nWe extract the spectrum in two small (2 pixel) apertures. One of these is centred on the transient position, while the second is centred on the location of the emission lines. Since the offset between these two locations is only ∼ 0 . 2 arcseconds there is naturally some contamination of each spectrum with light from both sources, but this is minimised by the use of small extraction apertures. The resulting 1D spectra are shown in Figure 3. The counterpart is very red, with a sharp break at 2 microns and an apparent emission feature at 2.15 microns. The spectrum then continues to rise to a possible second feature (or a change in the associated spectral slope) at around 4.5 microns.', '4.1 Prompt emission': "GRB 230307A is an exceptionally bright GRB. It has the second highest fluence of any GRB observed in more than 50 years of GRB observations [52]. While it remains a factor of 50 less fluent than GRB 221009A, it is still a factor ∼ 2 brighter than GRB 130427A, the third brightest burst. Bursts with these extreme fluences are rare. In Figure 6, we plot the distribution of observed fluence for Fermi /GBM detected bursts. At the brighter end, the slope of the distribution is consistent with the expected -3 / 2 slope for a uniform distribution of sources. The extrapolation of this relation suggests that bursts like GRB 230307A should occur once every several decades. Notably, three bursts well above the extrapolation (GRB 130427A, GRB 230307A, GRB 221009A) may indicate that bright bursts arise more frequently than expected. However, observationally it is clear that GRB 230307A is, at least, a once-per-decade event.  \nThe prompt light curve of GRB 230307A (Figure 1) shows two distinct emission features: an initial episode of hard emission from the trigger until ≈ 18 s, then a softer episode from ≈ 19 s onwards. These distinct episodes of hard and soft emission are strongly reminiscent of the long-duration merger GRB 211211A, but the initial pulse complex is ∼ 50 per cent longer in GRB 230307A when compared to the ∼ 12 s duration seen in GRB 211211A [17]. The relative durations of the initial pulse complex in the two GRBs bear a striking resemblance to their relative time-averaged peak energies [936 ± 3 keV vs 647 ± 8 keV; 2, 104]. In GRB 211211A, substantial spectral evolution was seen to drive the light curve, and the underlying radiation mechanism was identified as fast-cooling synchrotron emission [20]. The coherent development of the hardness ratio (Figure 1, lower) indicated similar spectral evolution in GRB 230307A, which the spectral analysis confirmed. Indeed, as described in Section 2.1.1, the time-resolved spectral analysis of the prompt emission revealed the presence of two spectral breaks in the GBM band, E break and E peak , coherently becoming softer from 7.5 s up to 19.5 s. Also, in this case, the spectral indices indicate synchrotron emission in the marginally fast-cooling regime. From 19.5 s onwards (approximately when the softer and dimmer emission episode starts), the low-energy break E break is continuously approaching the lower limit of the GBM band (8 keV), presumably crossing it to enter the X-ray regime. Unfortunately, the lack of simultaneous observations in X-rays with another telescope, e.g. Swift /XRT, prevents us from fully tracing the evolution of the spectral break down to X-rays at later times, as was done for GRB 211211A.  \nThe time-averaged Fermi /GBM spectrum of GRB 230307A across the T 90 interval is best fit with a cutoff power-law with α = 1 . 07 ± 0 . 01 and cutoff energy 936 ± 3 keV [2]. From this, we calculate a hardness ratio (the ratio of the 50 - 300 keV photon flux to the 10 - 50 keV photon flux) of 0 . 88 +0 . 01 -0 . 02 . This is higher than the value for 211211A (0 . 57) but comfortably within the 1 σ distribution of hardness ratios for canonical long GRBs (i.e. with T 90 > 2 s) in the Fermi catalogue, which we calculate to be 0 . 66 +0 . 51 -0 . 29 from the data in von Kienlin et al. [105]. Like GRB 211211A before it, GRB  \n230307A appears to have 'typical' long GRB properties in terms of its time-averaged hardness ratio and its T 90 . This strengthens the case for a significant number [17, 18] of long-duration GRBs having been mistakenly identified as stellar collapse events.  \nHowever, in some ways, GRB 230307A differs significantly from several of the other brightest GRBs. For example, the afterglow was relatively faint, while the burst was very bright. In Figure 7, we plot the prompt fluence in the 15-150 keV band against the X-ray afterglow brightness at 11 hours [updated from 106, 107]. The general trend between the afterglow brightness and fluence is seen; the best-fit slope to this relation is approximately one. So, while there is substantial scatter, there is a direct proportionality between the fluence and the afterglow brightness. Notably, while the afterglow and prompt emission of GRB 221009A were exceptionally bright (after correcting for the heavy foreground extinction), they were in keeping with this relatively broad relationship. GRB 230307A is different. Here we extrapolate the X-ray flux to 11 hours based on the measured X-ray flux at ∼ 1 day and the decay slope. We also re-calculate the GRB 230307A fluence in the relevant 15-150 keV energy band for comparison to Swift /BAT. This burst is a notable outlier in the relation, with a faint X-ray flux for its extraordinary prompt brightness. The afterglow brightness depends both on the energy of the burst and the density of the interstellar medium; it is, therefore, possible that the location in this fluence - afterglow brightness plane is indicative of a lowdensity medium, which would be consistent with expectations for such a large GRB host offset.  \nIt is also of interest that another burst in a similar location is GRB 211211A. This long burst has a clear signature of kilonova emission within its light curve. If GRB 230307A is a similar event, faint afterglows (relative to the prompt emission) may be an effective route for disentangling mergers from collapsars.  \nTo further compare the ratio between the X-ray brightness and the γ -ray fluence, we retrieve the X-ray light curve of all Swift -detected GRBs from the Swift Burst Analyser [108] and limit the sample to 985 long GRBs and 55 short GRBs with at least two XRT detections and measured BAT fluence. The fluences are taken from the Swift /BAT Gamma-Ray Burst Catalog 4 [109] and represent the measurements from 15 to 150 keV integrated over the total burst duration. We add to this sample the GRBs 170817A [off-axis GRB; 110] and 221009A [brightest GRB detected to date; 52]. For the former, we retrieve the X-ray light curve from Hajela et al. [111] and use a γ -ray fluence of 2 . 4 × 10 -7 erg cm -2 [110]. For the latter, we take the X-ray light curve from the Swift Burst Analyser and assume a fluence of 0 . 007 erg cm -2 (corrected from the 1-10000 keV fluence in [52] to the 15-150 keV band). Following [112], we resample the X-ray light curves and normalise them by the γ -ray fluence on a grid defined by the observed F X / Fluence ratios and the time-span probed by the data. If no data are available at a specific time of the grid, we linearly interpolate between adjacent observations but do not extrapolate any data. Hence the paucity of observations at later times reflects the last time at which sources were detected by the Swift /XRT.  \nShort and long GRBs occupy the same part of the F X / Fluence vs time parameter space (Figure 1). In contrast, GRB 230307A has an unprecedentedly low F X / Fluence ratio that is almost 10-fold lower than the faintest GRBs at the same time. To  \nemphasise the uniqueness of GRB 230307A, we also show in the same figure the Swift /BAT-detected GRBs 050925, 051105A, 070209, 070810B, 100628A, 130313A, 170112A that evaded detection with Swift /XRT. The limits on their F X / Fluence ratio (shown by downward pointing triangles in that figure) are consistent with the observed range of F X / Fluence ratios, ruling out a selection bias against GRBs with lower than usual F X / Fluence. Intriguingly, GRBs 080503, 191019A and 211211A had markedly low F X / Fluence ratios during the shallow decline phase of their X-ray light curves. Furthermore, GRB 211211A reached a value of 1 . 2 × 10 -9 s -1 at 120 ks, comparable to GRB 230307A.", '4.2 Counterpart Evolution': 'Although the afterglow of GRB 230307A was promptly detected thanks to TESS, this data was not available to the community for several days. Further follow-up was, therefore, much slower, and the counterpart was not discovered until the localisation was narrowed down to several sq. arcminutes, approximately 24 hours after the burst. The result is that the counterpart is poorly sampled (particularly in colour) during the early phases, while later observations suffer from typically modest signal-to-noise.  \nThe TESS observations detected a relatively bright (though not exceptional given the fluence of the burst) outburst, coincident with the prompt emission, likely peaking at I < 15 [12]. The afterglow was much fainter, apparently no brighter than I = 18 in the minutes to hours after the burst was detected. It was relatively flat during this period, with a power-law through the first to last TESS observations decaying as F ( t ) ∝ t -0 . 2 . The TESS and ground-based observations can be consistently modelled with a forward shock afterglow + kilonova (see Section 6.1).  \nThere are no simultaneous colours at the time of the first ground-based afterglow detections (1.4 days), although extrapolation of the r -band detection with ULTRACAM to the WHITE detection with the Swift /UVOT suggests a relatively red colour (WHITE-r = 1.6 ± 0.4) . However, such an interpretation is difficult due both to the large photometric errors and the width of the WHITE filter on the Swift /UVOT.  \nOptical observations obtained multiple colours at an epoch ∼ 2 . 4 days post-burst. These show the afterglow to have a blue colour with g = 22 . 35 ± 0 . 26, i = 21 . 68 ± 0 . 09 and z = 21 . 8 [80]. This is consistent with GRB afterglows in general(i.e F ν ∝ ν -β gives β ≈ 1). Observations in the near infrared (NIR) were not undertaken until ∼ 10.4 days post-burst. However, these reveal a relatively bright K-band source. The inferred i-K(AB) > 2.9 at this epoch is very red. Interpreted as a change in the spectral slope, it is β ≈ 2 . 5. The K-band light hence appears to be in significant excess with respect to the afterglow expectations based both on optical data and on the X-ray light curve.  \nIt is relevant to consider if such an excess could arise via extinction. However, this is not straightforward to explain. For a generic β = 1 slope we expect i -K ( AB ) ≈ 1 . 1. At z = 0, to obtain i -K ( AB ) = 2 . 9 would correspond to a foreground extinction of A V ≈ 4. However, this would also predict g -i ≈ 3, which is entirely inconsistent with the earlier observations. This problem becomes more acute for higher redshifts, where the bluer bands probe increasingly into the UV.  \nThe IR excess becomes extremely prominent by the time of the JWST observations. At 28.5 days, the source is detected in all bands but is very faint in the NIRCam blue  \nchannel (F070W, F115W, F150W) and rises rapidly (in F ν ) through the redder bands (F277W, F356W, F444W). Expressed as a power-law, this is β ≈ 3 . 1 in the 2-5 micron region, and β ≈ 1 between 0.7-1.5 microns. This does not match the expectations for any plausible spectral break in a GRB afterglow or any plausible extinction (where one would expect the slope to steepen towards the blue). This strongly implies that the red excess seen in the K-band at ten days and with JWST at 28.5 days is some additional component. Indeed, in the JWST observations, the other component, beginning at around 2 microns, is very clearly visible in both photometry and spectroscopy.  \nThis component evolves exceptionally rapidly. In the K-band, the inferred decay rate from 11.5 to 28.5 days is ∼ t -3 . 5 expressed as a power-law or ∼ 0 . 25 mag per day, if exponential. This is much faster than observed in GRB afterglows or supernovae. It is, however, consistent with the expectations for kilonovae. As shown in Figure 4, the overall evolution shows substantial similarity with AT2017gfo. To constrain the temporal and spectral evolution within a plausible physical model more accurately, we fit the multi-band photometry with afterglow and kilonova models. The outputs of these models are described in detail in section 6.', '4.3 Identification of the host galaxy': "Deep optical imaging of the field identifies several relatively bright galaxies in the vicinity of the sky position of GRB 230307A. Our preferred host galaxy is the brightest of these, which we denote as G1. It lies at z = 0 . 065 and is offset 30 arcseconds (40 kpc in projection) from the location of the afterglow. Following the method of [113] this galaxy has a probability of chance alignment of P chance ∼ 0 . 09 (see also [80]). Although this is not extremely low, and so is only suggestive of a connection to the transient, we note that i) the luminosity of the late time counterpart at this redshift is very similar to AT2017gfo and ii) the spectral feature seen at 2.1 microns in AT2017gfo matches with the emission feature seen in the JWST spectroscopy of GRB 230307A. This is a broad line, but assuming they have the same physical origin, they fix the redshift to the range 0 . 04 < z < 0 . 08. G1 is the only galaxy within this range in the field. The physical properties of this galaxy are outlined in section 4.4.  \nOur MUSE observations provide redshifts for this galaxy and several others, also identifying a small group of galaxies (G2, G3, G4) at a common redshift of z = 0 . 263. All of these galaxies have P chance values substantially greater than our preferred host. Furthermore, because of the larger redshifts, the implied offsets from GRB 230307A are ≫ 100kpc. This is larger than seen for any short GRB with a firmly identified host. We, therefore, disfavour these as plausible host galaxies for GRB 230307A.  \nDeep JWST observations reveal no evidence of a directly underlying host galaxy for GRB 230307A, as would be expected if it had a collapsar origin. In particular, at late times, the faint source at the counterpart's location is consistent with a pointsource (i.e. a subtraction of the PSF constructed by WebbPSF yields no significant residuals). However, we identify a faint galaxy, undetected in the blue and with F277W =27 . 9 ± 0 . 1, offset only 0.3 arcseconds from the burst position. We designate this galaxy H1.  \nOur NIRSpec observations provide a redshift of z = 3 . 87 for H1 based on the detection of [O III] (5007) and H α . At this redshift, the offset is only ∼ 1 . 3 kpc.  \nAlthough many z ∼ 4 galaxies are extremely compact [114], it seems likely that some stellar population from this galaxy does extend under the burst position, and there may be marginal evidence for extension in this direction in the F444W image. However, this region is neither UV-bright nor an emission line region where one may expect to observe massive stars.  \nThe galaxy photometry, performed in 0.1-arcsecond apertures and subsequently corrected for encircled energy assuming point-source curves is F070W > 29.0, F115W=28.4 ± 0.3, F150W=28.6 ± 0.4, F277W=27.9 ± 0.1, F444W=28.3 ± 0.1, and the galaxy is only robustly detected in the redder bands (see Figure 2). We note that because of the proximity of the afterglow, we use a smaller aperture than may be optimal, although the galaxy is also compact.  \nWe can estimate the probability of chance alignment of this source with the GRB position via various routes. In principle, one can use number counts of galaxies on the sky in the multiple bands. These have recently been updated based on the first observations with JWST to provide number counts in appropriate bands [115]. We find that P chance , following the approach of [113] to be in the range ∼ 3-6 % for F277W and F444W (with no bound in the filters where the galaxy is undetected). Alternatively, we also estimate the probability directly from the data. We extract sources within the field via Source Extractor to create a mask of objects within the field. In the brightest detection (F277W) approximately 5% of the image is covered with objects of equal or brighter magnitude to H1, and we note that the burst position is not contained within this mask. This suggests that in this particular field, P chance > 5%.  \nThe absolute magnitude of H1 is M i ∼ -17 . 7, and the H α star formation rate is approximately 1 M ⊙ yr -1 . The half-light radius of the galaxy is approximately 0.1 arcseconds (700 pc). Although limited information is available, these values are generally consistent with those of the long GRB population. The burst offset from its host galaxy is ∼ 2 . 5 half-light radii. This is large but within the range seen for long-GRBs [116].  \nIn our X-shooter and MUSE observations there is no trace visible in 1D or 2D extractions at the source position, although a weak continuum is seen in the X-shooter spectrum when heavily binned. This is consistent with its faint magnitude at the time of the observations. At the location of Ly α at z = 3 . 87 we place limits of F < 2 . 5 × 10 -17 erg s -1 cm -2 assuming an unresolved line.  \nWe also examined both spectra for any emission lines at other redshifts. This is worthwhile given the strong emission lines often seen in long GRB hosts [117], which may make emission line redshifts possible, even if the host itself is undetected. However, despite deep observations, there are no visible emission lines consistent with no directly underlying host galaxy, consistent with a compact object merger, but not a collapsar.  \nUnsurprisingly, there are also numerous faint galaxies in the JWST images. However, all of these have large P chance values, and we do not consider them plausible host galaxies.  \nTaken a face value, the probability of chance alignment for G1 (our preferred host) and H1 ( z = 3 . 87) is similar. However, the luminosity, lightcurve evolution and spectroscopic feature at the redshift of G1 offer strong support for it as the host galaxy of GRB 230307A. Furthermore, there is no straightforward, reasonably viable  \nphysical model that could explain the burst's extreme properties at z = 3 . 87. This scenario would require extreme energetics, exceptionally rapid evolution and yields unphysical outcomes in standard GRB or supernovae scenarios. We outline this in detail in section 7.1.", '4.4 Host galaxy properties': 'To better understand the properties of G1, the likely GRB host galaxy, we performed a fit to both the MUSE spectrum and photometric measurements from the far-UV to the mid-IR. For the photometric measurements, we retrieved science-ready coadded images from the Galaxy Evolution Explorer (GALEX) general release 6/7 [118], DESI Legacy Imaging Surveys (LS; [119]) data release 9, and re-processed WISE images [120] from the unWISE archive [121] 5 . The unWISE images are based on the public WISE data and include images from the ongoing NEOWISE-Reactivation mission R7 [122, 123]. We measured the brightness of the galaxy G1 using the Lambda Adaptive Multi-Band Deblending Algorithm in R ( LAMBDAR [124]) and the methods described in [125]. We augment the SED with Swift /UVOT photometry in the u band and our 6-band JWST/NIRCAM photometry. The photometry on the UVOT images was done with uvotsource in HEASoft and an aperture encircling the entire galaxy. For JWST photometry, we used a 6-arcsec circular aperture, which allows us to gather all the observed light observed in JWST filters from the host galaxy. All measurements are summarised in Table 7.4.  \nTo derive the main physical properties of the host galaxy, such as its stellar mass, we employ two separate methodologies based on the photometric and spectroscopic data available for the host, and finally compare the results to assess the robustness of our conclusions. We first fit the multi-wavelength (0.1-4.4 µ m) dataset using the prospector python package [126], which allows us to model the host galaxy spectrum starting from its main constituents, namely a set of stellar population base spectra, built from the Flexible Stellar Population Synthesis (FSPS) package [127], and combined with a specific star-formation history (SFH) model. Moreover, we have also considered a fixed attenuation model based on the Calzetti [128] attenuation curve, and an additional nebular model originating from the gas component, which is built using the Cloudy photo-ionization code [129], and considering the FSPS stellar population as ionising sources. We have adopted a parametric SFH model, which is described by a delayed-exponential model where the star-formation rate varies as a function of time t = t age -t lt , with t lt being the lookback time [126], as SFR ∝ ( t/τ ) exp( -t/τ ), with τ being the e -folding time. We finally have used the dynesty Speagle [130] ensemble sampler to reconstruct the posterior distribution.  \nThe results of the prospector analysis are shown in Fig. 8. We obtain a mass value of the living stars of M ∗ = 2 . 37(+0 . 24 , -0 . 35) × 10 9 M ⊙ yr -1 . The mass of all stars ever formed is 0 . 20(+0 . 02 , -0 . 04) dex larger. The light-weighted stellar age resulting from the fit is 1 . 13(+1 . 49 , -0 . 36) Gyr.  \nAn alternative to parametric SED fits is to use synthetic stellar population SEDs as templates and combine them to fit the galactic spectra (the underlying assumption being instantaneous star formations rather than continuous functions of time). We  \ncan use the spectral synthesis from the BPASSv2.2.2[131, 132] binary populations and create templates with hoki [133] that are compatible with the ppxf fitting package[134], as described in [135]. Because SED fitting has a high level of degeneracies (see [134]), at first we do not fit all 13 BPASS metallicities at once with ppxf , as this can result in unphysical results (see discussion in Stevance et al. 135); instead we fit the metallicities individually to find which ones result in the best fits on their own. We find that a low Z (0.001) population and solar metallicity population (Z=0.014) result in decent fits, but the low metallicity population fails to predict a young stellar component that is seen in the images, whilst the solar metallicity fit fails to accurately match the H β and neighbouring absorption features in the blue part of the spectrum. So we then fit the galaxy simultaneously with Z=0.001 and Z=0.014 templates, and retrieve a good fit shown in Figure 9 alongside the recovered SFH.  \nWe find evidence of three main stellar populations: > 95% of the mass is found in lower metallicity (Z=0.001) stars with ages ranging from a few Gyr to 10 Gyr, with a peak of star formation around 5 Gyr; > 4.7% of the mass originates from a solar metallicity population (Z=0.014) that formed around 400 Myr ago; finally a small fraction ( < 0.05%) of the stellar mass in the host originates from the star-forming regions with ages a few Myrs.  \nThe details of the age distributions and exact metallicity values can be model dependent so we also fit the integrated galaxy spectrum with the single stellar population synthesis code STARLIGHT [136], which uses stellar populations based on 25 different ages and six metallicity values [137], and a Chabrier IMF [138]. The SFH retrieved by this method is more complex and would require odd configurations (including some high metallicties at old ages and low metallicities around 100 Myr, which is counter-intuitive, unless inflow from pristine gas will trigger a burst of SFR), but it also finds that overall the galaxy is dominated by an old population with lower metallicity and has a younger component at higher metallicity. In Figure 9 we show a comparison of the STARLIGHT and BPASS fits in the bottom left panel and see that they are very similar, despite STARLIGHT containing 6 different metallicities and assuming solely single star populations. This highlights the level of degeneracy we face when performing galaxy SED fits. We leave further comparisons to a follow-up study dedicated to the host and the progenitor populations of GRB 230307A, where we will also present detailed, specially resolved, fits to the datacube including its kinematics.  \nFor now we use the BPASS integrated fits to infer the stellar mass and the star formation rate of the host of GRB 230307A, as the fit and SFH is more convincing that the one obtained with STARLIGHT. We find that there are currently M ∗ = 1 . 65 × 10 9 M ⊙ of living stars (corresponding to 3.1 × 10 9 M ⊙ at ZAMS) in G1. Using the nebular component retrieved from subtracting the fit of the stellar component to the observed data, we can also estimate the star formation rate and metallicity. From the H α feature we estimate that the SFR is 5.47 ± 0 . 30 × 10 -1 M ⊙ yr -1 using the Kennicutt formulation [139], and using the N2 index, in the CALIFA formulation [140], we infer an oxygen abundance of 8.20 ± 0.16 (12 + log(O/H)).  \nThere are qualitative similarities between the host of GRB 230307A and NGC 4993[135], the host of the first confirmed kilonova (they are both dominated by an older stellar populations and include a younger more metal rich component), but there  \nare some key differences: NGC 4993 was a lenticular galaxy without a clear young component, whereas the host of GRB 230307A shows clear spiral arms and star forming regions. Another major difference is that the metallicity of the old population in this galaxy is 10 times lower than that of NGC 4993 (Z=0.001 compared to Z=0.010), which will influence the stellar evolution of potential progenitors. Finally, NGC 4993 had a large stellar (and presumably dark halo) mass M ∗ ≈ 10 11 M ⊙ [135], a factor of > 50 larger than the host of GRB 230307A.  \nThe location of GRB 230307A relative to its host galaxy is consistent with these properties. In particular, the low mass of the galaxy suggests a modest gravitational potential such that binaries with velocities of a few hundred km s -1 can readily escape. The large offset also suggests that the binary is formed from the older stellar population.', '4.5 Properties of the brightest GRBs': 'GRB 230307A is the second brightest 6 burst observed in over 50 years of observations [52]. If it arises from a compact object merger, this implies that such bright bursts can be created in mergers. Indeed, such a picture appears likely based on GRB 211211A [17-19], the sixth brightest burst. Of the ten brightest bursts observed by the Fermi /GBM, and subsequently localised at the arcsecond level, three have apparently secure associations with supernovae (GRB 130427A, GRB 171010A, GRB 190114C), and two (GRB 211211A, GRB 230307A) are associated with kilonovae, and hence mergers. Of the remaining five, one lies at z = 1 . 4 and has energetics which suggest a collapsar; three have no redshift information, although one of these (GRB 160821A) lies in proximity to several galaxies at z = 0 . 19; and one is GRB 221009A whose associated with a supernova remains unclear [97, 141, 142], although recent observations suggest a collapsar with an associated supernova is most likely [143]. Within this very bright population, collapsars are likely as common as mergers.', '5 Event rates': "One key question of interest is the likely event rate for such merger GRBs. A simple estimate of the event rate associated to a single event is given by  \nR = 1 Ω tV max . (1)  \nHere Ω reflects the fraction of the sky covered by the detection mission, t the effective mission duration (accounting for the duty cycle) and V max the maximum co-moving volume within which a burst with the same properties could be identified.  \nFor GRB 230307A, Ω = 0 . 65 (average for the Fermi /GBM) and t ≈ 15 years. V max is more complicated: as shown in Figure 6, the fluence distribution for GBM bursts extends to ∼ 10 -8 erg cm -2 and is likely complete to around 10 -6 erg cm -2 . Given the extreme brightness of GRB 230307A, it would likely have been recovered to a distance ∼ 50 times greater than its observed distance. If at z = 0 . 065 the inferred z max = 2 . 03  \nor V max = 630 Gpc -3 . In this case, the inferred rate of such bursts becomes extremely small, R ≈ 1 . 6 × 10 -4 Gpc -3 yr -1 . However, in practice, such bursts would not readily be identified at such redshifts since neither supernova nor kilonova signatures could be observed. A more realistic estimate would correspond to the distance at which associated supernovae can be either identified or ruled out with moderate confidence. In this case z max = 0 . 5 (also adopted by [144]), V max = 29 Gpc 3 , and R ≈ 3 . 5 × 10 -3 Gpc -3 yr -1 .  \nThese rate estimates also assume that GRB 230307A is the only merger-GRB to have occurred within the 15-year lifetime of the Fermi /GBM. This is almost certainly not the case. Indeed, GRB 211211A was also identified by Fermi /GBM and has rather similar estimates of the intrinisc rate [5 . 7 × 10 -3 Gpc -3 yr -1 , 144].  \nHowever, even the interpretation of ∼ 2 events is problematic. In particular, the V/V max for GRB 230307A is 0 . 004, and for GRB 211211A = 0 . 005 (again assuming z max = 0 . 5). For a sample average of uniformly distributed sources of comparable energy or luminosity, we expect V/V max ∼ 0 . 5). That the initial identification of such a population should arise from bursts with such extreme V/V max values is surprising, but may reflect that these bursts are the brightest, which likely encouraged a detailed follow-up. However, it is improbable that they represent the only such bursts observed, and we should expect a much larger population.  \nTo better quantify this, we extend our analysis to the Swift bursts and utilise the fluence of GRB 230307A converted to a 15-150 keV equivalent fluence using the observed spectral parameters. At z = 0 . 065, E iso (15-150) keV ∼ 7 × 10 51 erg, and for GRB 211211A E iso (15-150) keV = 2 × 10 51 erg. As expected, low energy events dominate the low redshift GRB population. However, at z < 0 . 5, there are 12 (out of 42) bursts with E iso ≳ 10 51 erg. This includes some further supernova-less GRBs, in particular GRB 060614 ( E iso = 9 × 10 50 erg), GRB 191019A ( E iso = 2 . 0 × 10 51 erg), and some bursts for which supernova searches have not been reported (e.g. GRB 150727A, GRB 061021, and the 'ultra-long' GRB 130925A). This sets an upper limit on the number of bursts at low redshift, which may be associated with mergers. In practice, selection effects would support a scenario where mergers generate a larger fraction of these bursts. In particular, the afterglows of GRB 230307A and GRB 211211A appear to be faint, despite the bright prompt emission. Such afterglows are difficult to find and may evade detection. In these cases redshifts may only be obtained from host galaxies. The associations may not be obvious if the bursts are offset from host galaxies at moderate redshifts. Such follow-up may occur late after the burst, or optical afterglow non-detections may lead to a lack of optical/IR follow-up because of uncertainty regarding the optical brightness of the event or suggestions it may be optically dark because of host galaxy extinction. Finally, given the afterglow brightness issues, it is possible that the small fraction of bright GRBs without redshift measurements may arise from a similar channel. These observations would imply that between 30-70% bursts at z < 0 . 5 and E iso ≳ 10 51 erg could arise from mergers, although it is likely less. A modest number of events at higher redshift is consistent with the observations, and would alleviate concerns regarding V/V max for GRB 211211A and GRB 230307A.  \nThis fraction is surprisingly high given the strong evidence that long GRBs arise from broad-lined type Ic supernovae and short GRBs from compact object mergers.  \nHowever, the dominant contributors to the long-GRB supernova connection occur at low energy, and belong to a population of low luminosity GRBs (LLGRBs) [145]. In a significant number of these, we may observe a energy source in the prompt emission separate from the highly relativistic jet seen in on-axis, energetic bursts. For example, the long-lived, soft nature of some bursts suggests a contribution from shock breakout or cocoon emission. If, for this reason, the luminosity function of collapsar GRBs is steeper at low luminosity than that of merger-GRBs, it is possible that at low luminosity the long GRB population is dominated by collapsars, while at high luminosity the contribution of mergers is significant. Such an interpretation is not without problems, given the star-forming nature of long-GRB hosts and their typically small offsets from their host galaxies. However, it is a logical investigation for future work.", '6.1 Light curve modelling': "In order to shed light onto the properties of the jet and, even more importantly, to separate the contribution of the kilonova from that of the jet afterglow in the UVOIR bands, we modelled the multi-wavelength light curves from radio to X-rays as a superposition of synchrotron emission from the forward shock driven by the jet into the interstellar medium (ISM), following [32, 146], and blackbody emission from the photophere of a kilonova, using the simple single-component model of [147].  \nThe forward shock synchrotron emission model has eight parameters, namely the isotropic-equivalent kinetic energy in the jet E K , its initial bulk Lorentz factor Γ 0 , its half-opening angle θ j , the ISM number density n , the fraction ξ N of ISM electrons that undergo diffusive shock acceleration in the forward shock, the fraction ϵ e of the shock downstream internal energy that is shared by such electrons, the slope p of the power law d N e / d γ ∝ γ -p that describes the Lorentz factor (as measured in the shock downstream comoving frame) distribution of the accelerated electrons as they leave the acceleration region, and the fraction ϵ B of the shock downstream internal energy that is shared by a small-scale, turbulence-driven, random magnetic field. The shock hydrodynamics is computed from energy conservation and accounts for the lateral expansion of the shock [146]. The effective electron energy distribution is computed accounting for the cooling induced by synchrotron and synchrotron-self-Compton emission, including an approximate treatment of the Klein-Nishina suppression of the Thomson cross section [146]. In computing flux densities, the synchrotron surface brightness of the shock is integrated over equal-arrival-time surfaces to account for the effects of relativistic aberration and latitude-dependent retarded times on the spectral shape [148].  \nThe kilonova model [147] assumes spherical ejecta expanding homologously, v = r/t , and featuring a power law density profile ρ ( r, t ) ∝ t -3 v -δ between a minimum and a maximum velocity, v ej ≤ v ≤ v ej , max . The density normalization is set by the total ejecta mass M ej . In general, the model allows for the ejecta opacity (assumed grey) κ to be piecewise-constant within the profile, but here we assume a uniform opacity across the ejecta for simplicity. The model divides the ejecta into 100 small  \nshells and computes the heating rate and thermalization efficiency within each. This allows for the derivation of the internal energy evolution in each shell and eventually the computation of the photospheric luminosity L KN in the diffusion approximation. The fixed ejecta opacity also allows for the computation of the optical depth and hence for the identification of a photospheric radius, which then sets the effective temperature T KN by the Stefan-Boltzmann law. In our modelling of GRB 230307A, we computed the flux density by simply assuming pure blackbody emission with the given luminosity and effective temperature at each given time. We fixed v max = 0 . 6 c and left M ej , v ej , κ and δ as free parameters.  \nTo carry out the model fitting, we defined an asymmetric Gaussian log-likelihood term for the i -th datapoint, which corresponds to an observation at time t i and in a band whose central frequency is ν i , as  \nln L i = -1 2 ( F ν, m ( ν i , t i ) -F ν, obs , i ) 2 σ 2 i + f 2 sys F 2 ν, m -ln [ √ 2 π ( σ 2 l ,i + f 2 sys F 2 ν, m ) + √ 2 π ( σ 2 h ,i + f 2 sys F 2 ν, m ) ] , (2)  \nwhere F ν, m ( ν, t ) is the flux density predicted by the model, F ν, obs ,i is the measured flux density, the one-sigma error reflects the potentially asymmetric error bars  \nσ i = { σ l ,i if F ν, m ( ν i , t i ) ≤ F ν, obs ,i σ h ,i if F ν, m ( ν i , t i ) > F ν, obs ,i , (3)  \nand we introduced a fractional systematic error contribution f sys , which we take as an additional nuisance parameter, to account for potential inter-calibration uncertainties between different instruments and for the fact that error bars typically only account for statistical uncertainties. For X-ray detections, we fit the integrated flux and the spectral index independently, with an analogous term for each (but with no systematic error contribution for the spectral index). Upper limits were treated simply by setting F ν, obs ,i equal to the reported upper limit, σ h ,i = F ν, obs ,i / 10 and σ l ,i = 10 F ν, obs ,i . The final log-likelihood was taken as the sum of these terms.  \nIn order to derive a posterior probability density on our 13-dimensional parameter space, we assumed the priors reported in Table 6 and we sampled the posterior with a Markov Chain Monte Carlo approach using the emcee python module [149], which implements the Goodman and Weare [150] affine-invariant ensemple sampler. The medians and 90% credible intervals of the marginalised posteriors on each parameter obtained in this way are reported in Table 6. The posterior is visualised by means of corner plots in Figures 10 (jet afterglow parameters), 11 (kilonova parameters) and 12 (all parameters).  \nThe left-hand panel in Figure 13 shows the observed light curve data (markers) along with the best-fitting model (solid lines). Dashed lines single out the contribution of the kilonova. The right-hand panel in the same figure shows some selected spectra, showing in particular the good agreement of the first JWST epoch with the blackbody plus power law spectrum implied by our model at those times.  \nWhile the best-fit model demonstrates a relatively good agreement with most of the measurements, some discrepancies stand out, most prominently with the 61.5 d  \nJWST data and with the 28.5 d Chandra detection. The former is not too surprising, as the assumptions in the kilonova model (in particular that of blackbody photopsheric emission, which is particularly rough in such a nebular phase, and that of constant and uniform grey opacity, due to recombination of at least some species) are expected to break down at such late epochs. The latter is linked to the steepening ('jet break') apparent at around 2 days in the model X-ray light curve, which in turn is mainly driven by the need to not exceed the optical and near-infrared fluxes implied by observations at around one week and beyond. In absence of these constraints, the fit would have accommodated a larger jet half-opening angle, postponing the jet break and hence allowing for a better match with the best-fit Chandra flux. On the other hand, as noted in Methods, this flux is rather uncertain, with the low-end uncertainty possibly extending to fluxes lower by one order of magnitude or more, depending on the adopted prior in the spectral analysis (see Methods). Still, such a discrepancy might indicate the presence of additional X-ray emission that is not accounted for by the model, as has been seen previously in e.g. [151, 152].", '6.2 Spectral analysis modeling': "The JWST/NIRSpec spectrum taken on 5 April 2023 exhibits a red continuum component with emission line features. The most distinctive feature is a broad emission line at 2 . 15 microns (in the rest frame, assuming z = 0 . 065). This may be a blend (visibly split in Figure 3) and a simultaneous fit of two Gaussians provides measured centroids of 20285 ± 10 ˚ A and 22062 ± 10 ˚ A. The line widths are both consistent at v FWHM = 19100kms -1 (0.064c). This 2 . 1 micron feature is quite similar in strength and width to the 2 . 07 micron feature in AT2017gfo at 10.5 days after merger [discussed in 37]. The AT2017gfo line also appears to be better fit as a blend of two features rather than a single transition, with line velocities of v FWHM = 38900kms -1 . While the average line centre is reasonably consistent between the two, the components inferred for AT2017gfo and the kilonova of GRB 230307A are each quite different. Reference [37] finds them at 20590 ˚ A and 21350 ˚ A and there is no consistent velocity shift that could be applied to match AT2017gfo with our JWST spectrum. Nevertheless, the similarity in their average line centroids, velocities and equivalent widths is striking, as demonstrated in Figure 3.  \nWith a Doppler broadening parameter of ≲ 0 . 1 c , it is unlikely that the continuum component is formed as a result of the superposition of emission lines. Because kilonova radiation transfer at such late times is not yet fully understood, here we attempt to model the spectrum with the assumption that the emission consists of blackbody radiation from the photosphere and forbidden emission lines of heavy elements formed outside the photosphere.  \nIf the continuum is described with blackbody radiation, the temperature and photospheric velocity are ≈ 670 K and ≈ 0 . 08 c , respectively. The continuum luminosity is estimated as ∼ 2 × 10 39 erg / s in the NIRSpec band and ∼ 5 × 10 39 erg / s if the blackbody emission extends to much longer wavelengths. Assuming this emission is entirely powered by radioactivity of r -process nuclei, these correspond to an ejecta mass of ∼ 0 . 03-0 . 07 M ⊙ [147]. With the ejecta mass and velocity, the opacity is required to be ≳ 5 cm 2 / g in order to keep the ejecta optically thick at 30 day. It is worth noting  \nthat such a high opacity in the mid-IR indicates that the inner part of the ejecta is lanthanide rich [153-155].  \nForbidden emission lines in the infrared are expected to arise from fine structure transitions of low-lying energy levels of heavy elements. Most abundant ions are expected to produce the strongest lines. We attribute the strongest observed line at 2.15 microns to tellurium (Te) III from an M1 line list of heavy elements presented in [45], where the line wavelengths are experimentally calibrated according to the NIST database [156]. Te belongs to the second r -process peak. With the M1 line list, we model kilonova emission line spectra under the assumption that photons from forbidden lines produced outside the photosphere freely escape from the ejecta. The collision strengths of Te III are taken from an R-matrix calculation [43] and those of other ions are obtained by using an atomic structure code HULLAC[157]. The abundance pattern is chosen to be the solar r -process but we separate 'light' and 'heavy' elements at an atomic mass of 85 and introduce a parameter, the abundance ratio of the two (see figure 14). The ionization fractions are fixed to be ( Y +1 , Y +2 , Y +3 ) = (0 . 2 , 0 . 5 , 0 . 3) motivated by the Te ionization evolution in kilonova ejecta [41]. The line shape is approximated by a Gaussian with a line broadening velocity of 0 . 08 c , which is the same as the photospheric velocity. The mass in the line forming region is estimated by assuming that the observed line luminosity, 5 × 10 38 erg s -1 , is locally generated by radioactivity of r -process nuclei, corresponding to ∼ 0 . 02 M ⊙ . Given the abundance pattern and ionization state, the mass of Te III in the line forming region is ≈ 8 · 10 -4 M ⊙ . The electron temperature of the line forming region is then determined such that the total line luminosity agrees with the observed one. The estimated electron temperature is ∼ 3000 K, which is slightly higher than that derived from the pure neodymium nebular modeling [158]. This is because the cooling by tellurium ions is more efficient than neodymium.  \nWe find that [Te III] 2 . 10 µ m line is indeed the most outstanding emission line around 2 microns. Several weaker lines also contribute to the flux around 3-4 microns. There is another potential line feature around 4 . 5 microns in the NIRSpec spectrum. The location of this feature is consistent with [Se III] 4 . 55 µ m and [W III] 4 . 43 µ m as pointed out by [45] for the kilonova AT2017gfo. From the spectral modeling, we obtain the total ejecta mass of ∼ 0 . 05-0 . 1 M ⊙ , which agrees with the one obtained from the light curve modeling ∼ 0 . 1 M ⊙ .  \nHere we show a brief estimate of the Te III mass from the observed line at 2 . 15 microns ( 3 P 0 -3 P 1 ). The collisional excitation rate per Te III ion from the ground level ( 3 P 0 ) to the first exited level ( 3 P 1 ) is given by  \nk 01 = 8 . 63 · 10 -6 n e √ T e Ω 01 g 0 e -E 01 /kT e s -1 , (4)  \nwhere n e and T e are the thermal electron density and temperature, Ω 01 ≈ 5 . 8 is the collision strength [43], E 01 ≈ 0 . 6 eV is the excitation energy, g 0 is the statistical weight of the ground level. Assuming that the ejecta mass in the line forming region is 0 . 02 M ⊙ expanding with 0 . 08 c and the ions are typically doubly ionised, we estimate  \nn e ∼ 3 · 10 5 cm -3 , and thus, the line emissivity per Te III ion is  \nϵ 10 ≈ 2 . 5 · 10 -14 ( n e 3 · 10 5 cm -3 ) erg / s , (5)  \nwhere T e = 3000K is used. Combining the line emissivity with the observed line luminosity in 2 . 25 ± 0 . 23 µ m, L line ≈ 3 · 10 38 erg / s, we obtain  \nM (Te III) ≈ 10 -3 M ⊙ ( n e 3 · 10 5 cm -3 ) -1 ( L line 3 · 10 38 erg / s ) . (6)  \nThe mass estimated from the line is somewhat dependent on T e and n e . However, we emphasise that, with T e ≈ 3000 K and n e ≈ 3 · 10 5 cm -3 , the line luminosity is consistent with the radioactive power in the line forming region. It is also interesting to note that the Te III mass of 10 -3 M ⊙ is in good agreement with the one obtained based on the same line seen in AT2017gfo at 10.5 day [42].  \nWhile we conclude that the observed line feature at 2.1 microns is most likely attributed to Te III, it is important to note that there are caveats associated with our modeling. One obvious caveat is that the model does not include E1 lines. Lanthanides and actinides have E1 transitions between low-lying levels in the mid-IR [discussed in 37]. Due to their lower abundances, these lines are expected to be weaker compared to the Te line if collisional excitation dominates the excitation processes. However, as we make an implicit assumption that their E1 lines contribute to the opacity in the mid-IR, they may produce P-Cygni like features, see, e.g., [37, 159]. For example, Ce III has a strong line at 2 . 07 µ m with log gf = -1 . 67 [159]. We estimate that its line optical depth is ≲ 0 . 1 at 700 K with ∼ 0 . 05 M ⊙ and ∼ 0 . 1 c even if Ce is purely in Ce III. However, more careful analyses including non-LTE effects are needed to quantify it. Another caveat is that the opacity of lanthanides is expected to have some wavelength dependence. Including this effect may also affect the spectral modelings.", '7 Alternative progenitor possibilities': 'Our interpretation of GRB 230307A provides a self-consistent model for the source in which the temporal and spectral evolution, as well as the source location, can be readily explained. The kilonova has marked similarities with AT2017gfo providing a robust indication of its origin, and we do not need to postulate new and unseen phenomena to explain it. However, it is also relevant to consider alternative possibilities. In particular, given the location of the galaxy at z = 3 . 87, it is important to consider if the burst could originate at that redshift.', '7.1 GRB 230307A as a high redshift, highly energetic GRB': "The nearby galaxy H1 (F277W(AB)=27.5 ± 0.1 , r proj = 0 . 3 arcsec) with a spectroscopic redshift of z=3.87 has a relatively low probability of occurring by chance ( ∼ 5 -10%, see section 4.3). This galaxy has a comparable P chance to G1 (the z = 0 . 065 galaxy). The host-normalised offset for H1 is ∼ 2 . 5, which is large but not unprecedented for long GRBs [116]. However, assuming the late time light at the GRB position  \nis all from the transient, it does not lie on the stellar field of this galaxy, which is unusual, for example, in the samples of [116, 160, 161], there is only one (of > 100) sub-arcsecond localised GRB not on the stellar field of its host.  \nAt z = 3 . 87, the inferred isotropic energy release and luminosity of GRB 230307A would be E iso = 1 . 2 × 10 56 erg and L iso = 1 . 7 × 10 56 erg s -1 (using a 64 ms peak flux). This is approximately an order of magnitude more energetic and two orders of magnitude more luminous than any other previously identified GRB [52].  \nIf at z = 3 . 87, we can have some confidence that GRB 230307A would be the most energetic burst ever detected by Swift or Fermi , including those without redshift or even afterglow identifications. In Table 4, we tabulate the most fluent GRBs observed by Fermi . Most of these have either redshifts or optical detections, which constrain z < 6 via the detection of the source in the optical band. This leads to a set of measured or maximum E iso values. For events without any redshift information, we can place a conservative upper redshift limit of z = 16. No GRBs detected by Fermi without a redshift can have energy over 10 56 erg unless they lie beyond z ∼ 20. Hence, GRB 230307A is sufficiently rare, if at z = 3 . 87, that events like it occur less frequently than once per decade across the Universe (i.e. no more than one in the combined lifetimes of Swift and Fermi ).  \nThe energetics of the burst at this redshift lie would lie far beyond those of the general GRB population and beyond those suggested as the upper limit for GRB energetics [162]. The only population of core-collapse GRBs whose energetics have been suggested to approach this value are those from first-generation population III stars [163-165]. It is not expected that such stars should exist at z ∼ 4. However, while a pop-III origin may alleviate energy concerns, the properties of the GRB and its optical/IR counterpart do not resemble the predictions for pop-III stars. In particular, pop III GRBs are suggested to have particularly long durations given the mass and radii of their progenitors [165], and so require extremely long durations of engine activity to enable jet breakout. However, the ∼ 35 s duration of GRB 230307A and its rapid variability do not readily fit this expectation.  \nIf one ascribes the GRB to a stellar collapse event, considering the afterglow's properties and associated supernovae is also relevant. Firstly, at 28.5 days, the JWST spectral observations are inconsistent synchrotron emission, suggesting that the counterpart must be dominated by another source in the mid-IR and the earlier K-band points. This excess, which in our preferred model is explained by a kilonova, would have to be due to the supernova or shock breakout if at z = 3 . 87. The K-band (rest frame B-band) light would reach a peak of M B ( AB ) < -23 . 5 on a timeframe of < 2-days (rest-frame) before decaying at a rate of > 1 mag day -1 for the next four days, or as a power-law decay, a rate of approximately t -4 . This appears too rapid for radioactively powered transients, at least based on the sample to date (note that at z = 3 . 87, the timescales are a factor of ∼ 5 faster than in the z = 0 . 065 scenario due to cosmological time dilation). The most likely option for such emission would be shock-breakout, which may begin blue but rapidly cool. There are simulations for the shock breakout associated with pop-III supernovae, which show an early peak [166]. However, this emission peaks in the UV to soft X-ray regime and at luminosities below that seen in GRB 230307A. Indeed, taking 28.5 days as a baseline; for a plausible  \nmaximum Pop III radius (e.g. 2000 R ⊙ , [167]) and a luminosity of L bol ∼ 10 43 erg s -1 the inferred temperature is T ∼ 300 , 000 K. This is incompatible with the spectral shape seen in the counterpart to GRB 230307A which peaks at > 4 . 5 microns ( T < 3000K at z = 3 . 87). Dust or metal line blanketing could alleviate this discrepancy to some degree, but it would be extreme to explain the observations. It would also come at the cost of an even higher intrinsic luminosity. Conversely, the radius at which the luminosity and temperature would be consistent is extremely large ( ∼ 0 . 3 pc) and indeed would require super-luminal expansion to reach from a single explosion within the time since burst. These constraints become even more extreme for the Kband observations at 11.5 days, where the luminosity is > 50 times higher. However, we lack detailed information regarding the spectral shape at this time. We can conclude that a thermal transient launched at the time of GRB 230307A cannot explain the observed source at z = 3 . 87.  \nWe should consider if GRB 230307A could be related to an explosion which bears little to no similarity to long-GRB progenitors. Given the inferred energetics, this is not an unreasonable proposition. However, the emission is too bright and too fast for, for example, the fast blue optical transients (e.g. [168]), the fastest of which have halftimes of ∼ 4 days [169, c.f. < 1 day for GRB 230307A at z = 3 . 87]. A further alternative may be a relativistic tidal disruption event. This would face significant challenges with the rapid variability timescales seen in the prompt emission and the non-nuclearity of the source within the galaxy at z = 3 . 87. Putting aside these concerns, the peak optical/IR luminosity is comparable to AT2022cmc [170, 171], but the evolution is too rapid and the dynamic range too large.  \nFinally, it is possible that the red excess seen at later times is not directly related to the progenitor or the transient but is a result of the re-processing of the GRB radiation by material within the host galaxy. In particular, for GRB 211211A [172] suggest that an alternative explanation for the emission could be the heating of dust. However, this model also encounters significant issues at z = 3 . 87. In particular, the observed K -band excess is a rest-frame B-band excess, much bluer than expected for dust heating. If this represented the peak of the thermal spectrum, it would be above the sublimation temperature of the dust. Alternatively, if it were the blue tail of a much cooler black body, the luminosity would be extremely high.  \nA final challenge to the high -z scenario is that the afterglow is detected in the UVOT-white and ground-based g -bands. These observations all have substantial sensitivity blue-ward of Ly α at z=3.87. A typical column from the intergalactic medium should attenuate ∼ 50% of the light in these bands, inconsistent with observations. Indeed, for a typical β = 1 spectrum, we expect to observe white -i = 2.8 and g -i = 1 . 9, approximately 3 and 4 σ away from the observed colours. There is significant variation in the absorption strength as a function of the line of sight, so a low (or near zero) absorption column would alleviate this tension. The sample of [173] implies that at z ∼ 4, perhaps 10% of galaxies have such low absorption sight lines.  \nHence, while the proximate galaxy could indicate a high redshift, there are few other indications in the transient properties that would support this interpretation. In particular, neither standard thermal nor non-thermal emission can explain the observed counterpart properties. If the burst does arise from z = 3 . 87, it requires a  \nnew kind of explosion, unlike any seen until now. In practice, such explosions could be extremely rare: the volumetric rate of GRBs with E iso > 10 56 erg is minimal, but postulating them is unnecessary when a robust, physically motivated explanation can be obtained for a lower redshift solution.", '7.2 Other cosmological scenarios': "We should also consider a further option which is that GRB 230307A does not reside at either z = 0 . 065 or z = 3 . 87 but is a chance super-position with both galaxies. In this scenario, the actual host is undetected or is one of the other galaxies within the field. The absence of direct redshift measurements makes placing constraints on this scenario challenging. However, we can use the non-detection of the late-time JWST magnitudes to limit the brightness of a supernova component at any redshift.  \nTo quantify the exclusion of 'normal' long GRBs at intermediate redshift we utilize model light curves for SN 1998bw from MOSFiT [174] calculated at a range of redshift from 0 . 05 < z < 4 . 0 (we take z = 4 as an upper limit for the redshift the GRB based on the observed g-band detections together with the detection of continuum emission to 5300 ˚ A [13] and similar faint trace seen to ∼ 5100 ˚ Ain our X-shooter spectrum). At each redshift, we compare our observed JWST photometry in each band to the model predictions at that time and report the most constraining limit (e.g. the lowest ratio of F obs /F 98 bw ). This is shown in Figure 15. At all plausible redshifts, any supernova must be at least a factor ∼ 3 fainter than SN 1998bw at similar times. For any redshift where the burst energetics fall within the range seen in the bulk GRB population ( z < 1 . 2), any supernova must be a factor > 100 fainter than SN 1998bw. Hence, there is no route in which GRB 230307A can arise from a classical long GRB ( E iso < 10 55 erg with an associated broad-lined Ic supernova. The strength of this rejection is predominantly based on the faintness of the source in the bluer bands, whereas at the first epoch, the redder bands are substantially brighter than SN 1998bw at z = 3 . 87.  \nIf the burst lies at an intermediate redshift, dust models may become more appealing. In particular, for a low to moderate redshift, the luminosty and timescales may be a suitable match (e.g. for GRB 211211A [172] find plausible explanations at z ∼ 0 . 5). However, in this case, the lack of a supernova to extremely deep limits would be surprising, as would the non-detection of the host galaxy.", '7.3 Galactic objects': 'In the absence of a robust absorption redshift, it is also necessary to consider if GRB 230307A could arise from a Galactic system.  \nThe very faint magnitudes and extreme red colours observed at late times can effectively rule out X-ray binary outbursts. For example, with an M-dwarf companion (absolute magnitude 9), the distance to the source would be ∼ 100 kpc, and larger for any more massive star, while the late-time colours of the source are not stellar. In practice, given that the source is transient, and we may also expect some contribution from an accretion disc, the overall properties cannot be remedied within the accreting binary framework.  \nThe inferred energetics for Galactic systems are E iso ∼ 5 × 10 43 ( d/ 10 kpc) 2 . This energetic output is within the bounds of giant outbursts from magnetars. For example, the giant flare of SGR 1806-20 had an inferred isotropic energy of E iso ∼ 2 × 10 46 erg (see [175], although subsequent downward revisions in its distance lower this somewhat [176]). However, GRB 230307A does not appear to meet the requirements of a magnetar. In particular, it is at high Galactic latitude (-36 degrees) and far from any plausible star formation that could give rise to a young neutron star. Furthermore, the emission in all magnetar outbursts is dominated by a very short pulse followed by decaying emission in which the pulse period of the neutron star is visible. This is not the case for GRB 230307A.  \nGRBs have previously been suggested to arise from the tidal disruption and accretion of rocky material onto a neutron star [177], and such events are seen in the case of white dwarfs [178]. However, for accretion onto a neutron star, we would usually expect to observe a relatively soft outburst (e.g. in the model of [177] the temperature is ∼ 10 keV). The spectrum we observe for GRB 230307A consists of evolving synchrotron emission which does not contain a thermal component and is inconsistent with directly observing heating accreting material. Indeed, even for near direct accretion, it is not clear how such a spectrum would be formed in an accreting neutron star scenario. Indeed, in this case, the evolution to very low temperatures on a timescale of ∼ 60 days would also not be natural.  \nHence, we conclude that no known Galactic systems could explain the observed properties of GRB 230307A.', '7.4 GRB 230307A as a white-dwarf - neutron star merger': "A final alternative is that GRB 230307A is related to the merger of a white dwarf and a neutron star. Although this is still a 'compact object merger', such mergers are very different from those of neutron stars with another neutron star or a black hole. In particular, simulations show that no r -process material is produced [179], and so we should not expect the very red emission. Although there are suggestions that GRB 211211A could have been produced by a WD-NS merger [19, 59], it is unclear if these could readily explain the detailed spectrophotometric evolution of GRB 230307A.  \nWhite dwarf neutron star mergers are appealing, because the long-duration of the gamma-ray emission could suggest a less compact remnant merger event. The wider separation of the binary at the disruption of the white dwarf produces disk accretions times from 100 to 1000s, matching the long duration of these bursts [180]. However, we expect the accretion rate in the mergers to be low, producing less-powerful, and hence, less luminous GRBs. Current WD/NS merger simulations predict a range of light curves that span the emission from GRB 230307A. A Ca feature does exist that is close to the observed line feature, but the models are too blue to explain the shape of the spectra [179]. This is because WD/NS mergers do not produce elements much heavier than the iron peak elements. As we have noticed in matching kilonova models, these elements do not have strong lines beyond ∼ 20 , 000 ˚ A. The subsequent emission above these wavelengths is very weak in these models. In addition, WD/NS mergers are expected to have fairly week kicks to ensure that the binary remains bound and these mergers are expected to have much lower offsets than neutron star mergers. However,  \nthe mass of the best-candidate host galaxy is sufficiently low that the observed offset for this burst can be attained [181].  \nTable 1 Optical and IR observations of the optical counterpart of GRB 230307A. Limits are given at the 3 σ level.  \nTable 2 X-ray and radio observations of the afterglow of GRB 230307A.  \nTable 3 Properties of possible host galaxies for GRB 230307A. *Formally, because the galaxy is undetected in the r -band, P chance is unbounded. This probability is based on the magnitudes measured at other wavelengths.  \nTable 4 The properties of the brightest 10 bursts observed by Fermi /GBM. Eight of these have afterglow identifications which place limits on the redshift, while it is likely that all originate from z < 16. As expected for the (observationally) brightest bursts there is a preference for low redshift, with the bursts all arising at much less than the mean for Swift GRBs [182]. In 3-4 cases a supernova has been identified, meaning the bursts arise from core collapse (with GRB 221009A still ambiguous at the time of writing). In 2 cases a kilonova origin appears more likely, while for the remaining 4 the absence of afterglows or the redshifts render such diagnostics impossible. However, it would appear that samples of bright, long-GRBs may contain significant number of compact object mergers.  \nTable 5 Photometry (in the AB photometric system) of the host galaxy of GRB 230307A. No reddening correction was applied.  \nTable 6 Light curve model parameters, priors and posterior medians and 90% credible intervals.  \nFig. 5 Ground based optical and IR imaging of the afterglow of GRB 230307A. Panel (a) shows a legacy survey image of the field, the green trapezoid shows the IPN error box at the time of afterlgow discovery, with the blue dotted lines showing its ultimate refinement. The shaded regions represent observations taken with ULTRACAM 1.4 days post burst. Panel (b) shows the same image, zoomed in to a region around the afterglow. The red source in the centre is an unrelated foreground star (confirmed both by its measured proper motion and spectrum), and several other galaxies can also be seen. Panel (c) shows an ULTRACAM image, where the afterglow can be seen to the north-east of the star. The remaining panels (d,e,f) show the afterglow imaged on ∼ 10 -20 day timescales. At 10 days the source is undetected in deep FORS2 i-band imaging, but well detected in the K-band with Gemini-South. This very red colour was suggestive of an additional component over and above any afterglow, and motivated further follow-up.  \n<!-- image -->", 'Author contribution statements': 'AJL led the project, including the location of the afterglow and kilonova and the JWST observations. BPG first identified the source as a likely compact object merger, was co-PI of the Chandra observations, and contributed to analysis and writing. OS contributed to afterglow and kilonova modelling and led the writing of these sections. MB was involved in kilonova modelling, EB contributed to interpretation, placing the burst in context and high energy properties. KH was involved in kilonova spectral modelling and identified the 2.15 micron feature. LI reduced VLT/MUSE and X-shooter observations and led the host analysis. GPL contributed to afterglow and kilonova modelling. ARE analysed the Chandra observations, BS reduced and analysed VLT observations. NS contributed to afterglow and kilonova modelling. SS was responsible for placing the burst afterglow in context and demonstrating its faintness. NRT  \ncontributed to analysis, interpretation and writing. KA was involved in the ULTRACAM observations and interpretation. GA led the ATCA observations. GB reduced the JWST NIRCAM data. LC processed and analysed the MUSE observations. VSD is the ULTRACAM PI. JPUF contributed to the interpretation. WF was the PI on the Chandra observations and contributed to discussion. CF contributed to the theoretical interpretation. NG was involved in host analysis. KEH, GP, AR, SDV, SC, PDA, DH, MDP, CCT, AdUP and DA contributed to ESO observations and discussion. DW contributed to spectral and progenitor modelling and discussion. MJD, PK, SP, JM, SGP, IP, DIS contributed to the ULTRACAM observations. GL investigated potential similarities with other transients. AT, PAE, BS and JAK contributed to the Swift observations. MF extracted and flux-calibrated the TESS light curve. SJS analysed the JWST spectral lines, and contributed to interpretation and writing. HFS performed the BPASS-hoki-ppxf fits to the integrated MUSE flux and contributed the associate figure and text. All authors contributed to manuscript preparation through contributions to concept development, discussion and text.', 'Acknowledgements': "We dedicate this paper to David Alexander Kann, who passed on March 10. The final messages he sent were regarding follow-up of GRB 230307A, and we hope it would satisfy his curiosity to know the final conclusions.  \nThis work is based on observations made with the NASA/ESA/CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #4434 and 4445. Support for Program numbers 4434 and 4445 was provided through grants from the STScI under NASA contract NAS5- 03127.  \nThis paper is partly based on observations collected at the European Southern Observatory under ESO programme 110.24CF (PI Tanvir), and on observations obtained at the international Gemini Observatory (program IDs GS-2023A-DD-106), a program of NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini Observatory partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigaci'on y Desarrollo (Chile), Ministerio de Ciencia, Tecnolog'ıa e Innovaci'on (Argentina), Minist'erio da Ciˆencia, Tecnologia, Inova¸c˜oes e Comunica¸c˜oes (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). Processed using the Gemini IRAF package and DRAGONS (Data Reduction for Astronomy from Gemini Observatory North and South).  \nAJL, DBM and NRT were supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 725246).  \nN. Sarin is supported by a Nordita fellowship. Nordita is supported in part by NordForsk. B. Metzger is supported in part by the NSF (grant AST-2002577).  \nJ.H. and D.L. were supported by a VILLUM FONDEN Investigator grant (project number 16599). G.P.L. is supported by a Royal Society Dorothy Hodgkin Fellowship (grant Nos. DHF-R1-221175 and DHF-ERE-221005). G.L. was supported by a research grant (19054) from VILLUM FONDEN. K. H. is supported by JST FOREST Program (JPMJFR2136) and the JSPS Grant-in-Aid for Scientific Research (20H05639, 20H00158, 23H01169, 20K14513). KEH acknowledges support from the Carlsberg Foundation Reintegration Fellowship Grant CF21-0103. S. Schulze acknowledges support from the G.R.E.A.T. research environment, funded by Vetenskapsr˚adet , the Swedish Research Council, project number 2016-06012. The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation under grant No. 140. JPUF is supported by the Independent Research Fund Denmark (DFF-409000079) and thanks the Carlsberg Foundation for support. SJS acknowledges funding from STFC Grant ST/X006506/1 and ST/T000198/1. VSD and ULTRACAM are funded by STFC grant ST/V000853/1 AAB acknowledges funding from the UK Space Agency. MN is supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 948381) and by UK Space Agency Grant No. ST/Y000692/1. H.F.S is supported by the Eric and Wendy Schmidt AI in Science Postdoctoral Fellowship, a Schmidt Futures program. DS acknowledges funding from STFC grants ST/T000406/1, ST/T003103/1, ST/X001121/1. MER acknowledges support from the research programme Athena with project number 184.034.002, which is financed by the Dutch Research Council (NWO). POB acknowledges funding from STFC grant ST/W000857/1. D.K.G acknowledges support from the Australian Research Council Centre of Excellence for Gravitational Wave Discovery (OzGrav), through project number CE170100004.", 'Data availability': 'JWST data are directly available from the MAST archive. Chandra and Swift data are also in the public domain. ESO and Gemini data are stored in their respective archives and will be available to all once the proprietary period expires. Data can be obtained from the corresponding author between the date of publication and the end of the proprietary period. This research has made use of Fermi data which are publicly available and can be obtained through the High Energy Astrophysics Science Archive Research Center (HEASARC) website at https://heasarc.gsfc.nasa.gov/W3Browse/ fermi/fermigbrst.html', 'Code availability': "Much analysis for this paper has been undertaken with publically available codes and the details required to reproduce the analysis are contained within the manuscript.  \nFig. 6 The distribution of measure fluence from the Fermi /GBM catalogue [105]. The solid line shows an expected slope of -3 / 2 for a uniform distribution. The faint end deviates from this line because of incompleteness. At the brighter end, there are three bursts which appear to be extremely rare, GRB 130427A, GRB 221009A [52] and GRB 230307A. To indicate the apparent rarity we also plot lines representing the expected frequency of events under the assumption of a -3 / 2 slope. We would expect to observe bursts akin to GRB 230307A only once per several decades.  \n<!-- image -->  \nFig. 7 A comparison between the prompt fluence (in the 15-150 keV band) and the X-ray flux at 11 hours for Swift GRBs with GRB 221009A and GRB 230307A added, updated from [106, 107]. The general 1:1 trend between the prompt fluence and X-ray brightness can clearly be seen, although it has a significant scatter, although a very rare event, GRB 221009A apparently lies on the same relation. However, GRB 211211A and GRB 230307A are clearly outliers to this relation with GRB 230307A occupying a region devoid of other GRBs. This very faint afterglow compared to the prompt emission may be related to a location at large projected offset from its host galaxy in which the density of the ambient medium is very low.  \n<!-- image -->  \nFig. 8 Spectral energy distribution (SED) of the galaxy G1 from 1000 to 60,000 ˚ A (black data points) and its best fit with the prospector SED fitting code (grey shaded curve). In the top right, we also report the values of the model parameters and their 1 σ uncertainties. The red squares represent the model-predicted magnitudes. The fitting parameters are shown in the upper-left corner. The abbreviation 'n.o.f.' stands for the number of filters.  \n<!-- image -->  \nFig. 9 SED fit (top left panel) and Star Formation History (right hand panel of the host galaxy obtained with BPASSv.2.2.2 -hoki templates fit with ppxf . We also include a comparison of the final fits obtained with BPASS (2 metallicities) and STALRIGHT (6 metallicites) in the bottom left panel to highlight how similar both fits are.  \n<!-- image -->  \nFig. 10 Corner plot of posterior probability density from multi-wavelength light curve fitting, limited to the relativistic jet afterglow parameter space. Histograms on the diagonal show the marginalized posterior probability densities on the parameters constructed from our MCMC posterior samples. Dashed black lines show the 90% credible interval, while red lines show the medians. The remaining plots show the one, two and three sigma equivalent contours of the joint posterior probability densities of parameter pairs. Red lines and squares mark the medians.  \n<!-- image -->  \nFig. 11 Corner plot of posterior probability density from multi-wavelength light curve fitting, limited to the kilonova parameter space. Similar to figure 10, but for the kilonova model parameters.  \n<!-- image -->  \nFig. 12 Corner plot of posterior probability density from multi-wavelength light curve fitting. Similar to figures 10 and 11, but showing all model parameters, including the nuisance parameter f sys .  \n<!-- image -->  \n1 2  \n3  \n4  \n5  \nFig. 13 Multi-wavelength light curves and model predictions. Markers in the figure show the observed flux density at the position of GRB 230307A in various bands (see legend in the left-hand panel) and at various times. Downward-facing triangles represent upper limits. The optical and near infrared flux densities are multiplied by the numbers reported in the legend for presentation purposes. The butterfly-shaped filled regions in the right-hand panel encompass flux densities consistent at one, two and three sigma (progressively lighter shades) with the Swift /XRT and Chandra detections in the 0.3-10 keV band, according to our analysis and adopting a uniform prior on the flux. Solid lines of the corresponding colours show the predicted light curves (left-hand panel) and spectra (right-hand panel) of our afterglow (forward shock only) plus kilonova model at the central frequencies of the bands. Dashed lines single out the contribution of the kilonova.  \n<!-- image -->  \nFig. 14 Abundance used in the spectral modeling. The abundance is chosen based on the solar rprocess residuals. The abundance of the 'light' elements ( A < 85) is reduced relative to the solar pattern. The locations of selenium (Se), tellurium (Te) and tungsten (W) are marked.  \n<!-- image -->  \nFig. 15 Limits on supernova similar to SN 1998bw as a function of redshift, based on the most constraining detection with JWST . At any redshift for which GRB 230307A would not be the most energetic GRB ever observed, any supernova is at least a factor of ∼ 100 fainter than SN 1998bw.  \n<!-- image -->", 'References': "- [1] Fermi GBM Team: GRB 230307A: Fermi GBM Final Real-time Localization. GRB Coordinates Network 33405 , 1 (2023)\n- [2] Dalessi, S., Roberts, O.J., Meegan, C., Fermi GBM Team: GRB 230307A: Fermi GBM Observation of a very bright burst. GRB Coordinates Network 33411 , 1 (2023)\n- [3] Burns, E., Goldstein, A., Lesage, S., Dalessi, S.: Grb 230307a: possibly the second highest grb energy fluence ever identified. GRB Coordinates Network 33414 , 1 (2023)\n- [4] Svinkin, D., Frederiks, D., Ridnaia, A., Lysenko, A., Cline, T., Konus-Wind Team, von Kienlin, A., Zhang, X., Rau, A., Savchenko, V., Bozzo, E., Ferrigno, C., INTEGRAL SPI-ACS Team, Goldstein, A., Briggs, M.S., Wilson-Hodge, C., Fermi Gbm Team, Barthelmy, S., Cummings, J., Krimm, H., Palmer, D., Tohuvavohu, A., Swift-Bat Team: IPN triangulation of extremely bright GRB 221009A. GRB Coordinates Network 32641 (2022)\n- [5] Burrows, D.N., Gropp, J.D., Osborne, J.P., Page, K.L., D'Elia, V., Sbarufatti, B., D'Ai, A., Dichiara, S., Evans, P.A., Swift-XRT Team: GRB 230307A: SwiftXRT observations. GRB Coordinates Network 33429 , 1 (2023)\n- [6] Kozyrev, A.S., Golovin, D.V., Litvak, M.L., Mitrofanov, I.G., Sanin, A.B., Mgns/Bepicolombo, Hend/Mars Odyssey Teams, Benkhoff, J., Bepicolombo Team, Svinkin, D., Fredericks, D., Ridnaia, A., Lysenko, A., Cline, T., Konus-Wind Team, Goldstein, A., Briggs, M.S., Wilson-Hodge, C., Burns, E., Fermi Gbm Team, Bozzo, E., Ferrigno, C., INTEGRAL SPI-ACS Grb Team, Barthelmy, S., Cummings, J., Krimm, H., Palmer, D., Tohuvavohu, A., Swift-Bat Team, Boynton, W., Fellows, C., Harshman, K., Enos, H., Starr, R., Gardner, A.S., Grs-Odyssey Grb Team: Further improved IPN localization for GRB 230307A. GRB Coordinates Network 33461 , 1 (2023)\n- [7] Dey, A., Schlegel, D.J., Lang, D., Blum, R., Burleigh, K., Fan, X., Findlay, J.R., Finkbeiner, D., Herrera, D., Juneau, S., Landriau, M., Levi, M., McGreer, I., Meisner, A., Myers, A.D., Moustakas, J., Nugent, P., Patej, A., Schlafly, E.F., Walker, A.R., Valdes, F., Weaver, B.A., Y'eche, C., Zou, H., Zhou, X., Abareshi, B., Abbott, T.M.C., Abolfathi, B., Aguilera, C., Alam, S., Allen, L., Alvarez, A., Annis, J., Ansarinejad, B., Aubert, M., Beechert, J., Bell, E.F., BenZvi, S.Y., Beutler, F., Bielby, R.M., Bolton, A.S., Brice˜no, C., BuckleyGeer, E.J., Butler, K., Calamida, A., Carlberg, R.G., Carter, P., Casas, R., Castander, F.J., Choi, Y., Comparat, J., Cukanovaite, E., Delubac, T., DeVries, K., Dey, S., Dhungana, G., Dickinson, M., Ding, Z., Donaldson, J.B., Duan, Y., Duckworth, C.J., Eftekharzadeh, S., Eisenstein, D.J., Etourneau, T., Fagrelius, P.A., Farihi, J., Fitzpatrick, M., Font-Ribera, A., Fulmer, L., Gansicke, B.T., Gaztanaga, E., George, K., Gerdes, D.W., Gontcho, S.G.A., Gorgoni, C., Green,  \n- G., Guy, J., Harmer, D., Hernandez, M., Honscheid, K., Huang, L.W., James, D.J., Jannuzi, B.T., Jiang, L., Joyce, R., Karcher, A., Karkar, S., Kehoe, R., Kneib, J.-P., Kueter-Young, A., Lan, T.-W., Lauer, T.R., Le Guillou, L., Le Van Suu, A., Lee, J.H., Lesser, M., Perreault Levasseur, L., Li, T.S., Mann, J.L., Marshall, R., Mart'ınez-V'azquez, C.E., Martini, P., du Mas des Bourboux, H., McManus, S., Meier, T.G., M'enard, B., Metcalfe, N., Mu˜noz-Guti'errez, A., Najita, J., Napier, K., Narayan, G., Newman, J.A., Nie, J., Nord, B., Norman, D.J., Olsen, K.A.G., Paat, A., Palanque-Delabrouille, N., Peng, X., Poppett, C.L., Poremba, M.R., Prakash, A., Rabinowitz, D., Raichoor, A., Rezaie, M., Robertson, A.N., Roe, N.A., Ross, A.J., Ross, N.P., Rudnick, G., Safonova, S., Saha, A., S'anchez, F.J., Savary, E., Schweiker, H., Scott, A., Seo, H.-J., Shan, H., Silva, D.R., Slepian, Z., Soto, C., Sprayberry, D., Staten, R., Stillman, C.M., Stupak, R.J., Summers, D.L., Sien Tie, S., Tirado, H., Vargas-Maga˜na, M., Vivas, A.K., Wechsler, R.H., Williams, D., Yang, J., Yang, Q., Yapici, T., Zaritsky, D., Zenteno, A., Zhang, K., Zhang, T., Zhou, R., Zhou, Z.: Overview of the DESI Legacy Imaging Surveys. Astron. J. 157 (5), 168 (2019) https://doi. org/10.3847/1538-3881/ab089d arXiv:1804.08657 [astro-ph.IM]\n- [8] Levan, A.J., Gompertz, B., Ackley, K., Kennedy, M., Dhillon, V.: GRB 230307A: Possible optical afterglow. GRB Coordinates Network 33439 , 1 (2023)\n- [9] Lipunov, V., Buckley, D.A.H., Chasovnikov, A., Kechin, Y., Kuznetsov, A., Tiurina, N., Gress, O., Gorbovskoy, E., Antipov, G., Balanutsa, P., Zhirkov, K., Vlasenko, D., Senik, V., Kuvshinov, D., Topolev, V., Tselik, Y., Cheryasov, D., Gorbunov, I., Kechin, Y., Francile, C., Podesta, R., Podesta, F., Rebolo, R., Serra, M., Carrasco, L., Valdes, J.R., Chavushyan, V., Patino Alvarez, V.M., Martinez, J., Corella, A., Rodriguez, L.H., Budnev, N.M., Tlatov, A., Dormidontov, D., Gabovich, A., Yurkov, V.: GRB 230204B: MASTER detection of an object near NTT Optical Afterglow Candidate. GRB Coordinates Network 33441 , 1 (2023)\n- [10] O'Connor, B., Dichiara, S., Troja, E., Gillanders, J., Cenko, S.B., Kouveliotou, C.: GRB 230307A: Gemini-South Confirmation of the Optical Afterglow. GRB Coordinates Network 33447 , 1 (2023)\n- [11] Im, M., Paek, G.S.-H., Jeong, M., Chang, S.-W., Choi, H., Lee, C.-U.: GRB 230307A: GECKO optical observations and detections of the afterglow. GRB Coordinates Network 33449 , 1 (2023)\n- [12] Fausnaugh, M.M., Jayaraman, R., Vanderspek, R., Ricker, G.R., Burke, C.J., Colon, K.D., Fleming, S.W., Lewis, H.M., Mullally, S., Youngblood, A., Barclay, T., Burns, E., Latham, D.W., Seager, S., Winn, J.N., Jenkins, J.M.: Observations of GRB 230307A by TESS. arXiv e-prints, 2303-07319 (2023) https://doi.org/ 10.48550/arXiv.2303.07319 arXiv:2303.07319 [astro-ph.HE]  \n- [13] Gillanders, J., O'Connor, B., Dichiara, S., Troja, E.: GRB 230307A: Continued Gemini-South observations confirm rapid optical fading. GRB Coordinates Network 33485 , 1 (2023)\n- [14] Lattimer, J.M., Mackie, F., Ravenhall, D.G., Schramm, D.N.: The decompression of cold neutron star matter. Astrophys. J. 213 , 225-233 (1977) https: //doi.org/10.1086/155148\n- [15] Li, L.-X., Paczy'nski, B.: Transient Events from Neutron Star Mergers. Astrophys. J. Lett. 507 (1), 59-62 (1998) https://doi.org/10.1086/311680\n- [16] Metzger, B.D.: Kilonovae. Living Reviews in Relativity 23 (1), 1 (2019) https: //doi.org/10.1007/s41114-019-0024-0\n- [17] Rastinejad, J.C., Gompertz, B.P., Levan, A.J., Fong, W.-f., Nicholl, M., Lamb, G.P., Malesani, D.B., Nugent, A.E., Oates, S.R., Tanvir, N.R., de Ugarte Postigo, A., Kilpatrick, C.D., Moore, C.J., Metzger, B.D., Ravasio, M.E., Rossi, A., Schroeder, G., Jencson, J., Sand, D.J., Smith, N., Agu'ı Fern'andez, J.F., Berger, E., Blanchard, P.K., Chornock, R., Cobb, B.E., De Pasquale, M., Fynbo, J.P.U., Izzo, L., Kann, D.A., Laskar, T., Marini, E., Paterson, K., Escorial, A.R., Sears, H.M., Thone, C.C.: A kilonova following a long-duration gammaray burst at 350 Mpc. Nature 612 (7939), 223-227 (2022) https://doi.org/10. 1038/s41586-022-05390-w arXiv:2204.10864 [astro-ph.HE]\n- [18] Troja, E., Fryer, C.L., O'Connor, B., Ryan, G., Dichiara, S., Kumar, A., Ito, N., Gupta, R., Wollaeger, R.T., Norris, J.P., Kawai, N., Butler, N.R., Aryan, A., Misra, K., Hosokawa, R., Murata, K.L., Niwano, M., Pandey, S.B., Kutyrev, A., van Eerten, H.J., Chase, E.A., Hu, Y.-D., Caballero-Garcia, M.D., CastroTirado, A.J.: A nearby long gamma-ray burst from a merger of compact objects. Nature 612 (7939), 228-231 (2022) https://doi.org/10.1038/s41586-022-05327-3 arXiv:2209.03363 [astro-ph.HE]\n- [19] Yang, J., Ai, S., Zhang, B.-B., Zhang, B., Liu, Z.-K., Wang, X.I., Yang, Y.H., Yin, Y.-H., Li, Y., Lu, H.-J.: A long-duration gamma-ray burst with a peculiar origin. Nature 612 (7939), 232-235 (2022) https://doi.org/10.1038/ s41586-022-05403-8 arXiv:2204.12771 [astro-ph.HE]\n- [20] Gompertz, B.P., Ravasio, M.E., Nicholl, M., Levan, A.J., Metzger, B.D., Oates, S.R., Lamb, G.P., Fong, W.-f., Malesani, D.B., Rastinejad, J.C., Tanvir, N.R., Evans, P.A., Jonker, P.G., Page, K.L., Pe'er, A.: The case for a minutelong merger-driven gamma-ray burst from fast-cooling synchrotron emission. Nature Astronomy 7 , 67-79 (2023) https://doi.org/10.1038/s41550-022-01819-4 arXiv:2205.05008 [astro-ph.HE]\n- [21] Arcavi, I., Hosseinzadeh, G., Howell, D.A., McCully, C., Poznanski, D., Kasen, D., Barnes, J., Zaltzman, M., Vasylyev, S., Maoz, D., Valenti, S.: Optical emission from a kilonova following a gravitational-wave-detected neutron-star  \n[22] Coulter, D.A., Foley, R.J., Kilpatrick, C.D., Drout, M.R., Piro, A.L., Shappee, B.J., Siebert, M.R., Simon, J.D., Ulloa, N., Kasen, D., Madore, B.F., MurguiaBerthier, A., Pan, Y.-C., Prochaska, J.X., Ramirez-Ruiz, E., Rest, A., RojasBravo, C.: Swope Supernova Survey 2017a (SSS17a), the optical counterpart to a gravitational wave source. Science 358 (6370), 1556-1558 (2017) https://doi. org/10.1126/science.aap9811  \n- [23] Pian, E., D'Avanzo, P., Benetti, S., Branchesi, M., Brocato, E., Campana, S., Cappellaro, E., Covino, S., D'Elia, V., Fynbo, J.P.U., Getman, F., Ghirlanda, G., Ghisellini, G., Grado, A., Greco, G., Hjorth, J., Kouveliotou, C., Levan, A., Limatola, L., Malesani, D., Mazzali, P.A., Melandri, A., Møller, P., Nicastro, L., Palazzi, E., Piranomonte, S., Rossi, A., Salafia, O.S., Selsing, J., Stratta, G., Tanaka, M., Tanvir, N.R., Tomasella, L., Watson, D., Yang, S., Amati, L., Antonelli, L.A., Ascenzi, S., Bernardini, M.G., Boer, M., Bufano, F., Bulgarelli, A., Capaccioli, M., Casella, P., Castro-Tirado, A.J., Chassande-Mottin, E., Ciolfi, R., Copperwheat, C.M., Dadina, M., De Cesare, G., di Paola, A., Fan, Y.Z., Gendre, B., Giuffrida, G., Giunta, A., Hunt, L.K., Israel, G.L., Jin, Z.-P., Kasliwal, M.M., Klose, S., Lisi, M., Longo, F., Maiorano, E., Mapelli, M., Masetti, N., Nava, L., Patricelli, B., Perley, D., Pescalli, A., Piran, T., Possenti, A., Pulone, L., Razzano, M., Salvaterra, R., Schipani, P., Spera, M., Stamerra, A., Stella, L., Tagliaferri, G., Testa, V., Troja, E., Turatto, M., Vergani, S.D., Vergani, D.: Spectroscopic identification of r-process nucleosynthesis in a double neutron-star merger. Nature 551 (7678), 67-70 (2017) https://doi.org/10.1038/nature24298  \n[24] Smartt, S.J., Chen, T.-W., Jerkstrand, A., Coughlin, M., Kankare, E., Sim, S.A., Fraser, M., Inserra, C., Maguire, K., Chambers, K.C., Huber, M.E., Kruhler, T., Leloudas, G., Magee, M., Shingles, L.J., Smith, K.W., Young, D.R., Tonry, J., Kotak, R., Gal-Yam, A., Lyman, J.D., Homan, D.S., Agliozzo, C., Anderson, J.P., Angus, C.R., Ashall, C., Barbarino, C., Bauer, F.E., Berton, M., Botticella, M.T., Bulla, M., Bulger, J., Cannizzaro, G., Cano, Z., Cartier, R., Cikota, A., Clark, P., De Cia, A., Della Valle, M., Denneau, L., Dennefeld, M., Dessart, L., Dimitriadis, G., Elias-Rosa, N., Firth, R.E., Flewelling, H., Flors, A., Franckowiak, A., Frohmaier, C., Galbany, L., Gonz'alez-Gait'an, S., Greiner, J., Gromadzki, M., Guelbenzu, A.N., Guti'errez, C.P., Hamanowicz, A., Hanlon, L., Harmanen, J., Heintz, K.E., Heinze, A., Hernandez, M.-S., Hodgkin, S.T., Hook, I.M., Izzo, L., James, P.A., Jonker, P.G., Kerzendorf, W.E., Klose, S., Kostrzewa-Rutkowska, Z., Kowalski, M., Kromer, M., Kuncarayakti, H., Lawrence, A., Lowe, T.B., Magnier, E.A., Manulis, I., Martin-Carrillo, A., Mattila, S., McBrien, O., Muller, A., Nordin, J., O'Neill, D., Onori, F., Palmerio, J.T., Pastorello, A., Patat, F., Pignata, G., Podsiadlowski, P., Pumo, M.L., Prentice, S.J., Rau, A., Razza, A., Rest, A., Reynolds, T., Roy, R., Ruiter, A.J., Rybicki, K.A., Salmon, L., Schady, P., Schultz, A.S.B., Schweyer, T., Seitenzahl, I.R., Smith, M., Sollerman, J., Stalder, B., Stubbs, C.W., Sullivan, M., Szegedi,  \nH., Taddia, F., Taubenberger, S., Terreran, G., van Soelen, B., Vos, J., Wainscoat, R.J., Walton, N.A., Waters, C., Weiland, H., Willman, M., Wiseman, P., Wright, D.E., Wyrzykowski, glyph[suppress]L., Yaron, O.: A kilonova as the electromagnetic counterpart to a gravitational-wave source. Nature 551 (7678), 75-79 (2017) https://doi.org/10.1038/nature24303  \n[25] Tanvir, N.R., Levan, A.J., Gonz'alez-Fern'andez, C., Korobkin, O., Mandel, I., Rosswog, S., Hjorth, J., D'Avanzo, P., Fruchter, A.S., Fryer, C.L., Kangas, T., Milvang-Jensen, B., Rosetti, S., Steeghs, D., Wollaeger, R.T., Cano, Z., Copperwheat, C.M., Covino, S., D'Elia, V., de Ugarte Postigo, A., Evans, P.A., Even, W.P., Fairhurst, S., Figuera Jaimes, R., Fontes, C.J., Fujii, Y.I., Fynbo, J.P.U., Gompertz, B.P., Greiner, J., Hodosan, G., Irwin, M.J., Jakobsson, P., Jørgensen, U.G., Kann, D.A., Lyman, J.D., Malesani, D., McMahon, R.G., Melandri, A., O'Brien, P.T., Osborne, J.P., Palazzi, E., Perley, D.A., Pian, E., Piranomonte, S., Rabus, M., Rol, E., Rowlinson, A., Schulze, S., Sutton, P., Thone, C.C., Ulaczyk, K., Watson, D., Wiersema, K., Wijers, R.A.M.J.: The Emergence of a Lanthanide-rich Kilonova Following the Merger of Two Neutron Stars. Astrophys. J. Lett. 848 (2), 27 (2017) https://doi.org/10.3847/2041-8213/ aa90b6  \n[26] Soares-Santos, M., Holz, D.E., Annis, J., Chornock, R., Herner, K., Berger, E., Brout, D., Chen, H.-Y., Kessler, R., Sako, M., Allam, S., Tucker, D.L., Butler, R.E., Palmese, A., Doctor, Z., Diehl, H.T., Frieman, J., Yanny, B., Lin, H., Scolnic, D., Cowperthwaite, P., Neilsen, E., Marriner, J., Kuropatkin, N., Hartley, W.G., Paz-Chinch'on, F., Alexander, K.D., Balbinot, E., Blanchard, P., Brown, D.A., Carlin, J.L., Conselice, C., Cook, E.R., Drlica-Wagner, A., Drout, M.R., Durret, F., Eftekhari, T., Farr, B., Finley, D.A., Foley, R.J., Fong, W., Fryer, C.L., Garc'ıa-Bellido, J., Gill, M.S.S., Gruendl, R.A., Hanna, C., Kasen, D., Li, T.S., Lopes, P.A.A., Louren¸co, A.C.C., Margutti, R., Marshall, J.L., Matheson, T., Medina, G.E., Metzger, B.D., Mu˜noz, R.R., Muir, J., Nicholl, M., Quataert, E., Rest, A., Sauseda, M., Schlegel, D.J., Secco, L.F., Sobreira, F., Stebbins, A., Villar, V.A., Vivas, K., Walker, A.R., Wester, W., Williams, P.K.G., Zenteno, A., Zhang, Y., Abbott, T.M.C., Abdalla, F.B., Banerji, M., Bechtol, K., Benoit-L'evy, A., Bertin, E., Brooks, D., Buckley-Geer, E., Burke, D.L., Carnero Rosell, A., Carrasco Kind, M., Carretero, J., Castander, F.J., Crocce, M., Cunha, C.E., D'Andrea, C.B., da Costa, L.N., Davis, C., Desai, S., Dietrich, J.P., Doel, P., Eifler, T.F., Fernandez, E., Flaugher, B., Fosalba, P., Gaztanaga, E., Gerdes, D.W., Giannantonio, T., Goldstein, D.A., Gruen, D., Gschwend, J., Gutierrez, G., Honscheid, K., Jain, B., James, D.J., Jeltema, T., Johnson, M.W.G., Johnson, M.D., Kent, S., Krause, E., Kron, R., Kuehn, K., Kuhlmann, S., Lahav, O., Lima, M., Maia, M.A.G., March, M., McMahon, R.G., Menanteau, F., Miquel, R., Mohr, J.J., Nichol, R.C., Nord, B., Ogando, R.L.C., Petravick, D., Plazas, A.A., Romer, A.K., Roodman, A., Rykoff, E.S., Sanchez, E., Scarpine, V., Schubnell, M., Sevilla-Noarbe, I., Smith, M., Smith, R.C., Suchyta, E., Swanson, M.E.C., Tarle, G., Thomas, D., Thomas, R.C., Troxel, M.A., Vikram, V., Wechsler, R.H., Weller, J., Dark Energy Survey, Dark  \nEnergy Camera GW-EM Collaboration: The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. I. Discovery of the Optical Counterpart Using the Dark Energy Camera. Astrophys. J. Lett. 848 (2), 16 (2017) https://doi.org/10.3847/2041-8213/aa9059  \n- [27] Valenti, S., Sand, D.J., Yang, S., Cappellaro, E., Tartaglia, L., Corsi, A., Jha, S.W., Reichart, D.E., Haislip, J., Kouprianov, V.: The Discovery of the Electromagnetic Counterpart of GW170817: Kilonova AT 2017gfo/DLT17ck. Astrophys. J. Lett. 848 (2), 24 (2017) https://doi.org/10.3847/2041-8213/aa8edf\n- [28] Villar, V.A., Guillochon, J., Berger, E., Metzger, B.D., Cowperthwaite, P.S., Nicholl, M., Alexander, K.D., Blanchard, P.K., Chornock, R., Eftekhari, T., Fong, W., Margutti, R., Williams, P.K.G.: The Combined Ultraviolet, Optical, and Near-infrared Light Curves of the Kilonova Associated with the Binary Neutron Star Merger GW170817: Unified Data Set, Analytic Models, and Physical Implications. Astrophys. J. Lett. 851 (1), 21 (2017) https://doi.org/10.3847/ 2041-8213/aa9c84\n- [29] Abbott, B.P., Abbott, R., Abbott, T.D., et al. : GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 119 (16), 161101 (2017) https://doi.org/10.1103/PhysRevLett.119.161101 arXiv:1710.05832 [gr-qc]\n- [30] Drout, M.R., Chornock, R., Soderberg, A.M., Sanders, N.E., McKinnon, R., Rest, A., Foley, R.J., Milisavljevic, D., Margutti, R., Berger, E., Calkins, M., Fong, W., Gezari, S., Huber, M.E., Kankare, E., Kirshner, R.P., Leibler, C., Lunnan, R., Mattila, S., Marion, G.H., Narayan, G., Riess, A.G., Roth, K.C., Scolnic, D., Smartt, S.J., Tonry, J.L., Burgett, W.S., Chambers, K.C., Hodapp, K.W., Jedicke, R., Kaiser, N., Magnier, E.A., Metcalfe, N., Morgan, J.S., Price, P.A., Waters, C.: Rapidly Evolving and Luminous Transients from Pan-STARRS1. Astrophys. J. 794 (1), 23 (2014) https://doi.org/10.1088/0004-637X/794/1/23\n- [31] Kasliwal, M.M., Kasen, D., Lau, R.M., Perley, D.A., Rosswog, S., Ofek, E.O., Hotokezaka, K., Chary, R.-R., Sollerman, J., Goobar, A., Kaplan, D.L.: Spitzer mid-infrared detections of neutron star merger GW170817 suggests synthesis of the heaviest elements. Mon. Not. R. Astron. Soc. 510 (1), 7-12 (2022) https: //doi.org/10.1093/mnrasl/slz007 arXiv:1812.08708 [astro-ph.HE]\n- [32] Mei, A., Banerjee, B., Oganesyan, G., Salafia, O.S., Giarratana, S., Branchesi, M., D'Avanzo, P., Campana, S., Ghirlanda, G., Ronchini, S., Shukla, A., Tiwari, P.: Gigaelectronvolt emission from a compact binary merger. Nature 612 (7939), 236-239 (2022) https://doi.org/10.1038/s41586-022-05404-7 arXiv:2205.08566 [astro-ph.HE]\n- [33] Gehrels, N., Norris, J.P., Barthelmy, S.D., Granot, J., Kaneko, Y., Kouveliotou, C., Markwardt, C.B., M'esz'aros, P., Nakar, E., Nousek, J.A., O'Brien, P.T., Page, M., Palmer, D.M., Parsons, A.M., Roming, P.W.A., Sakamoto, T., Sarazin, C.L.,  \nSchady, P., Stamatikos, M., Woosley, S.E.: A new γ -ray burst classification scheme from GRB060614. Nature 444 (7122), 1044-1046 (2006) https://doi.org/ 10.1038/nature05376 arXiv:astro-ph/0610635 [astro-ph]  \n- [34] Michaglyph[suppress]lowskI, M.J., Xu, D., Stevens, J., Levan, A., Yang, J., Paragi, Z., Kamble, A., Tsai, A.-L., Dannerbauer, H., van der Horst, A.J., Shao, L., Crosby, D., Gentile, G., Stanway, E., Wiersema, K., Fynbo, J.P.U., Tanvir, N.R., Kamphuis, P., Garrett, M., Bartczak, P.: The second-closest gamma-ray burst: sub-luminous GRB 111005A with no supernova in a super-solar metallicity environment. Astron. Astrophys. 616 , 169 (2018) https://doi.org/10.1051/ 0004-6361/201629942 arXiv:1610.06928 [astro-ph.HE]\n- [35] Levan, A.J., Malesani, D.B., Gompertz, B.P., Nugent, A.E., Nicholl, M., Oates, S., Perley, D.A., Rastinejad, J., Metzger, B.D., Schulze, S., Stanway, E.R., Inkenhaag, A., Zafar, T., Agui Fernandez, J.F., Chrimes, A., Bhirombhakdi, K., de Ugarte Postigo, A., Fong, W.-f., Fruchter, A.S., Fragione, G., Fynbo, J.P.U., Gaspari, N., Heintz, K.E., Hjorth, J., Jakobsson, P., Jonker, P.G., Lamb, G.P., Mandel, I., Mandhai, S., Ravasio, M.E., Sollerman, J., Tanvir, N.R.: A longduration gamma-ray burst of dynamical origin from the nucleus of an ancient galaxy. arXiv e-prints, 2303-12912 (2023) https://doi.org/10.48550/arXiv.2303. 12912 arXiv:2303.12912 [astro-ph.HE]\n- [36] Wollaeger, R.T., Korobkin, O., Fontes, C.J., Rosswog, S.K., Even, W.P., Fryer, C.L., Sollerman, J., Hungerford, A.L., van Rossum, D.R., Wollaber, A.B.: Impact of ejecta morphology and composition on the electromagnetic signatures of neutron star mergers. Mon. Not. R. Astron. Soc. 478 (3), 3298-3334 (2018) https://doi.org/10.1093/mnras/sty1018 arXiv:1705.07084 [astro-ph.HE]\n- [37] Gillanders, J.H., Sim, S.A., Smartt, S.J., Goriely, S., Bauswein, A.: Modelling the spectra of the kilonova AT2017gfo - II: Beyond the photospheric epochs. arXiv e-prints, 2306-15055 (2023) https://doi.org/10.48550/arXiv.2306.15055 arXiv:2306.15055 [astro-ph.HE]\n- [38] Joshi, Y.N., Tauheed, A., Davison, I.G.: The analysis of the 5s 2 5p 2 , 5s5p 3 5s 2 5p5d, and 5s 2 5p6s configurations of Te III. Canadian Journal of Physics 70 , 740 (1992) https://doi.org/10.1139/p92-117\n- [39] Hotokezaka, K., Beniamini, P., Piran, T.: Neutron star mergers as sites of r-process nucleosynthesis and short gamma-ray bursts. International Journal of Modern Physics D 27 (13), 1842005 (2018) https://doi.org/10.1142/ S0218271818420051 arXiv:1801.01141 [astro-ph.HE]\n- [40] Gillanders, J.H., Smartt, S.J., Sim, S.A., Bauswein, A., Goriely, S.: Modelling the spectra of the kilonova AT2017gfo - I. The photospheric epochs. Mon. Not. R. Astron. Soc. 515 (1), 631-651 (2022) https://doi.org/10.1093/mnras/stac1258 arXiv:2202.01786 [astro-ph.HE]  \n- [50] Fong, W.-f., Nugent, A.E., Dong, Y., Berger, E., Paterson, K., Chornock, R.,  \nLevan, A., Blanchard, P., Alexander, K.D., Andrews, J., Cobb, B.E., Cucchiara, A., Fox, D., Fryer, C.L., Gordon, A.C., Kilpatrick, C.D., Lunnan, R., Margutti, R., Miller, A., Milne, P., Nicholl, M., Perley, D., Rastinejad, J., Escorial, A.R., Schroeder, G., Smith, N., Tanvir, N., Terreran, G.: Short GRB Host Galaxies. I. Photometric and Spectroscopic Catalogs, Host Associations, and Galactocentric Offsets. Astrophys. J. 940 (1), 56 (2022) https://doi.org/10.3847/1538-4357/ ac91d0 arXiv:2206.01763 [astro-ph.GA]  \n- [51] Grindlay, J., Portegies Zwart, S., McMillan, S.: Short gamma-ray bursts from binary neutron star mergers in globular clusters. Nature Physics 2 (2), 116-119 (2006) https://doi.org/10.1038/nphys214 arXiv:astro-ph/0512654 [astro-ph]\n- [52] Burns, E., Svinkin, D., Fenimore, E., Kann, D.A., Feliciano Agu'ı Fern'andez, J., Frederiks, D., Hamburg, R., Lesage, S., Temiraev, Y., Tsvetkova, A., Bissaldi, E., Briggs, M.S., Fletcher, C., Goldstein, A., Hui, C.M., Hristov, B.A., Kocevski, D., Lysenko, A.L., Mailyan, B., Racusin, J., Ridnaia, A., Roberts, O.J., Ulanov, M., Veres, P., Wilson-Hodge, C.A., Wood, J.: GRB 221009A, The BOAT. arXiv e-prints, 2302-14037 (2023) https://doi.org/10.48550/arXiv.2302.14037 arXiv:2302.14037 [astro-ph.HE]\n- [53] Virgili, F.J., Liang, E.-W., Zhang, B.: Low-luminosity gamma-ray bursts as a distinct GRB population: a firmer case from multiple criteria constraints. Mon. Not. R. Astron. Soc. 392 (1), 91-103 (2009) https://doi.org/10.1111/j. 1365-2966.2008.14063.x arXiv:0801.4751 [astro-ph]\n- [54] Sarin, N., Lasky, P.D., Nathan, R.S.: Missed opportunities: GRB 211211A and the case for continual gravitational-wave coverage with a single observatory. Mon. Not. R. Astron. Soc. 518 (4), 5483-5489 (2023) https://doi.org/10.1093/ mnras/stac3441 arXiv:2210.14938 [astro-ph.HE]\n- [55] Metzger, B.D., Quataert, E., Thompson, T.A.: Short-duration gamma-ray bursts with extended emission from protomagnetar spin-down. Mon. Not. R. Astron. Soc. 385 (3), 1455-1460 (2008) https://doi.org/10.1111/j.1365-2966.2008.12923. x arXiv:0712.1233 [astro-ph]\n- [56] Bucciantini, N., Metzger, B.D., Thompson, T.A., Quataert, E.: Short gamma-ray bursts with extended emission from magnetar birth: jet formation and collimation. Mon. Not. R. Astron. Soc. 419 (2), 1537-1545 (2012) https://doi.org/10. 1111/j.1365-2966.2011.19810.x arXiv:1106.4668 [astro-ph.HE]\n- [57] Rosswog, S.: Fallback accretion in the aftermath of a compact binary merger. Mon. Not. R. Astron. Soc. 376 (1), 48-51 (2007) https://doi.org/10.1111/j. 1745-3933.2007.00284.x arXiv:astro-ph/0611440 [astro-ph]\n- [58] Desai, D., Metzger, B.D., Foucart, F.: Imprints of r-process heating on fallback accretion: distinguishing black hole-neutron star from double neutron star mergers. Mon. Not. R. Astron. Soc. 485 (3), 4404-4412 (2019) https://doi.org/", '10.1093/mnras/stz644 arXiv:1812.04641 [astro-ph.HE]': "- [59] Zhong, S.-Q., Li, L., Dai, Z.-G.: GRB 211211A: A Neutron Star-White Dwarf Merger? Astrophys. J. Lett. 947 (2), 21 (2023) https://doi.org/10.3847/ 2041-8213/acca83 arXiv:2304.04009 [astro-ph.HE]\n- [60] Gottlieb, O., Issa, D., Jacquemin-Ide, J., Liska, M., Foucart, F., Tchekhovskoy, A., Metzger, B.D., Quataert, E., Perna, R., Kasen, D., Duez, M.D., Kidder, L.E., Pfeiffer, H.P., Scheeli, M.A.: Large-scale Evolution of Seconds-long Relativistic Jets from Black Hole-Neutron Star Mergers. arXiv e-prints, 2306-14947 (2023) https://doi.org/10.48550/arXiv.2306.14947 arXiv:2306.14947 [astro-ph.HE]\n- [61] Crockford, S.J.: Evolutionary Roots of Iodine and Thyroid Hormones in Cell-Cell Signaling. Integrative and Comparative Biology 49 (7785), 155-166 (2009) https://doi.org/https://www.jstor.org/stable/40306080\n- [62] Planck Collaboration, Aghanim, N., Akrami, Y., Ashdown, M., Aumont, J., Baccigalupi, C., Ballardini, M., Banday, A.J., Barreiro, R.B., Bartolo, N., Basak, S., Battye, R., Benabed, K., Bernard, J.-P., Bersanelli, M., Bielewicz, P., Bock, J.J., Bond, J.R., Borrill, J., Bouchet, F.R., Boulanger, F., Bucher, M., Burigana, C., Butler, R.C., Calabrese, E., Cardoso, J.-F., Carron, J., Challinor, A., Chiang, H.C., Chluba, J., Colombo, L.P.L., Combet, C., Contreras, D., Crill, B.P., Cuttaia, F., de Bernardis, P., de Zotti, G., Delabrouille, J., Delouis, J.-M., Di Valentino, E., Diego, J.M., Dor'e, O., Douspis, M., Ducout, A., Dupac, X., Dusini, S., Efstathiou, G., Elsner, F., Enßlin, T.A., Eriksen, H.K., Fantaye, Y., Farhang, M., Fergusson, J., Fernandez-Cobos, R., Finelli, F., Forastieri, F., Frailis, M., Fraisse, A.A., Franceschi, E., Frolov, A., Galeotta, S., Galli, S., Ganga, K., G'enova-Santos, R.T., Gerbino, M., Ghosh, T., Gonz'alez-Nuevo, J., G'orski, K.M., Gratton, S., Gruppuso, A., Gudmundsson, J.E., Hamann, J., Handley, W., Hansen, F.K., Herranz, D., Hildebrandt, S.R., Hivon, E., Huang, Z., Jaffe, A.H., Jones, W.C., Karakci, A., Keihanen, E., Keskitalo, R., Kiiveri, K., Kim, J., Kisner, T.S., Knox, L., Krachmalnicoff, N., Kunz, M., Kurki-Suonio, H., Lagache, G., Lamarre, J.-M., Lasenby, A., Lattanzi, M., Lawrence, C.R., Le Jeune, M., Lemos, P., Lesgourgues, J., Levrier, F., Lewis, A., Liguori, M., Lilje, P.B., Lilley, M., Lindholm, V., L'opez-Caniego, M., Lubin, P.M., Ma, Y.-Z., Mac'ıas-P'erez, J.F., Maggio, G., Maino, D., Mandolesi, N., Mangilli, A., Marcos-Caballero, A., Maris, M., Martin, P.G., Martinelli, M., Mart'ınez-Gonz'alez, E., Matarrese, S., Mauri, N., McEwen, J.D., Meinhold, P.R., Melchiorri, A., Mennella, A., Migliaccio, M., Millea, M., Mitra, S., MivilleDeschˆenes, M.-A., Molinari, D., Montier, L., Morgante, G., Moss, A., Natoli, P., Nørgaard-Nielsen, H.U., Pagano, L., Paoletti, D., Partridge, B., Patanchon, G., Peiris, H.V., Perrotta, F., Pettorino, V., Piacentini, F., Polastri, L., Polenta, G., Puget, J.-L., Rachen, J.P., Reinecke, M., Remazeilles, M., Renzi, A., Rocha, G., Rosset, C., Roudier, G., Rubi˜no-Mart'ın, J.A., Ruiz-Granados, B., Salvati, L., Sandri, M., Savelainen, M., Scott, D., Shellard, E.P.S., Sirignano, C., Sirri, G., Spencer, L.D., Sunyaev, R., Suur-Uski, A.-S., Tauber, J.A., Tavagnacco,  \n- D., Tenti, M., Toffolatti, L., Tomasi, M., Trombetti, T., Valenziano, L., Valiviita, J., Van Tent, B., Vibert, L., Vielva, P., Villa, F., Vittorio, N., Wandelt, B.D., Wehus, I.K., White, M., White, S.D.M., Zacchei, A., Zonca, A.: Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys. 641 , 6 (2020) https://doi.org/10.1051/0004-6361/201833910 arXiv:1807.06209 [astro-ph.CO]\n- [63] Xiong, S., Wang, C., Huang, Y., Gecam Team: GRB 230307A: GECAM detection of an extremely bright burst. GRB Coordinates Network 33406 , 1 (2023)\n- [64] Xiao, H., Krucker, S.: Solar Orbiter STIX observation of GRB 230307A. GRB Coordinates Network 33410 , 1 (2023)\n- [65] Casentini, C., Tavani, M., Pittori, C., Lucarelli, F., Verrecchia, F., Argan, A., Cardillo, M., Evangelista, Y., Foffano, L., Piano, G., Addis, A., Baroncelli, L., Bulgarelli, A., di Piano, A., Fioretti, V., Panebianco, G., Parmiggiani, N., Longo, F., Romani, M., Marisaldi, M., Pilia, M., Trois, A., Donnarumma, I., Menegoni, E., Ursi, A., Giuliani, A., Tempesta, P., Agile Team: GRB 230307A: AGILE/MCAL detection. GRB Coordinates Network 33412 , 1 (2023)\n- [66] Navaneeth, P.K., Waratkar, G., Vibhute, A., Bhalerao, V., Bhattacharya, D., Rao, A.R., Vadawale, S., AstroSat CZTI Collaboration: GRB 230307A: AstroSat CZTI detection. GRB Coordinates Network 33415 , 1 (2023)\n- [67] Dafcikova, M., Ripa, J., Pal, A., Werner, N., Ohno, M., Takahashi, H., Meszaros, L., Csak, B., Husarikova, N., Munz, F., Topinka, M., Kolar, M., Breuer, J.-P., Hroch, F., Urbanec, T., Kasal, M., Povalac, A., Hudec, J., Kapus, J., Frajt, M., Laszlo, R., Koleda, M., Smelko, M., Hanak, P., Lipovsky, P., Galgoczi, G., Uchida, Y., Poon, H., Matake, H., Uchida, N., Bozoki, T., Dalya, G., Enoto, T., Frei, Z., Friss, G., Fukazawa, Y., Hirose, K., Hisadomi, S., Ichinohe, Y., Kapas, K., Kiss, L.L., Mizuno, T., Nakazawa, K., Odaka, H., Takatsy, J., Torigoe, K., Kogiso, N., Yoneyama, M., Moritaki, M., Kano, T., GRBAlpha Collaboration.: GRB 230307A: Detection by GRBAlpha. GRB Coordinates Network 33418 , 1 (2023)\n- [68] Ripa, J., Dafcikova, M., Pal, A., Werner, N., Ohno, M., Meszaros, L., Csak, B., Takahashi, H., Munz, F., Topinka, M., Hroch, F., Husarikova, N., Breuer, J.-P., Hudec, J., Kapus, J., Frajt, M., Rezenov, M., Laszlo, R., Galgoczi, G., Uchida, N., Enoto, T., Frei, Z., Fukazawa, Y., Hirose, K., Matake, H., Hisadomi, S., Ichinohe, Y., Kiss, L.L., Mizuno, T., Nakazawa, K., Odaka, H., Torigoe, K., Svoboda, P., Daniel, V., Dudas, J., Junas, M., Gromes, J., Vertat, I., Vzlusat-2/Grb Payload Collaboration.: GRB 230307A: VZLUSAT-2 detection. GRB Coordinates Network 33424 , 1 (2023)\n- [69] Svinkin, D., Frederiks, D., Ulanov, M., Tsvetkova, A., Lysenko, A., Ridnaia, A., Cline, T., Konus-Wind Team: Konus-Wind detection of GRB 230307A. GRB Coordinates Network 33427 , 1 (2023)  \n- [70] Li, Z., Su, Y., HXI Team: GRB 230307A: ASO-S/HXI detection. GRB Coordinates Network 33438 , 1 (2023)\n- [71] Kozyrev, A.S., Golovin, D.V., Litvak, M.L., Mitrofanov, I.G., Sanin, A.B., Hend/Mars Odyssey Team, Svinkin, D., Lysenko, A., Ridnaia, A., Ipn, Goldstein, A., Briggs, M.S., Wilson-Hodge, C., Burns, E., Fermi Gbm Team, Bozzo, E., Ferrigno, C., INTEGRAL SPI-ACS Grb Team, Barthelmy, S., Cummings, J., Krimm, H., Palmer, D., Tohuvavohu, A., Swift-Bat Team, Boynton, W., Fellows, C., Harshman, K., Enos, H., Starr, R., Gardner, A.S., Grs-Odyssey Grb Team: IPN triangulation of GRB 230307A (long/very bright). GRB Coordinates Network 33413 , 1 (2023)\n- [72] Kozyrev, A.S., Golovin, D.V., Litvak, M.L., Mitrofanov, I.G., Sanin, A.B., Hend/Mars Odyssey Team, Svinkin, D., Frederiks, D., Ridnaia, A., Lysenko, A., Cline, T., Konus-Wind Team, Goldstein, A., Briggs, M.S., Wilson-Hodge, C., E. Burns, E., Fermi Gbm Team, Bozzo, E., Ferrigno, C., The INTEGRAL SPI-ACS Grb Team, Barthelmy, S., Cummings, J., H. Krimm, H., Palmer, D., Tohuvavohu, A., Swift-Bat Team, Boynton, W., Fellows, C., Harshman, K., Enos, H., Starr, R., Gardner, A.S., Grs-Odyssey Grb Team: Improved IPN localization for GRB 230307A. GRB Coordinates Network 33425 , 1 (2023)\n- [73] Dalessi, S., Fermi GBM Team: GRB 230307A: Bad Time Intervals for Fermi GBM data. GRB Coordinates Network 33551 , 1 (2023)\n- [74] Ravasio, M.E., Oganesyan, G., Ghirlanda, G., Nava, L., Ghisellini, G., Pescalli, A., Celotti, A.: Consistency with synchrotron emission in the bright GRB 160625B observed by Fermi. Astron. Astrophys. 613 , 16 (2018) https://doi.org/ 10.1051/0004-6361/201732245 arXiv:1711.03106 [astro-ph.HE]\n- [75] Ravasio, M.E., Ghirlanda, G., Nava, L., Ghisellini, G.: Evidence of two spectral breaks in the prompt emission of gamma-ray bursts. Astron. Astrophys. 625 , 60 (2019) https://doi.org/10.1051/0004-6361/201834987 arXiv:1903.02555 [astroph.HE]\n- [76] Dhillon, V.S., Marsh, T.R., Stevenson, M.J., Atkinson, D.C., Kerry, P., Peacocke, P.T., Vick, A.J.A., Beard, S.M., Ives, D.J., Lunney, D.W., McLay, S.A., Tierney, C.J., Kelly, J., Littlefair, S.P., Nicholson, R., Pashley, R., Harlaftis, E.T., O'Brien, K.: ULTRACAM: an ultrafast, triple-beam CCD camera for high-speed astrophysics. Mon. Not. R. Astron. Soc. 378 (3), 825-840 (2007) https://doi.org/10.1111/j.1365-2966.2007.11881.x arXiv:0704.2557 [astro-ph]\n- [77] Dhillon, V.S., Bezawada, N., Black, M., Dixon, S.D., Gamble, T., Gao, X., Henry, D.M., Kerry, P., Littlefair, S.P., Lunney, D.W., Marsh, T.R., Miller, C., Parsons, S.G., Ashley, R.P., Breedt, E., Brown, A., Dyer, M.J., Green, M.J., Pelisoli, I., Sahman, D.I., Wild, J., Ives, D.J., Mehrgan, L., Stegmeier, J., Dubbeldam, C.M., Morris, T.J., Osborn, J., Wilson, R.W., Casares, J., Mu˜nozDarias, T., Pall'e, E., Rodr'ıguez-Gil, P., Shahbaz, T., Torres, M.A.P., de Ugarte  \n- [85] McMahon, R.G., Banerji, M., Gonzalez, E., Koposov, S.E., Bejar, V.J., Lodieu,  \n- M., Perri, M., Gehrels, N., Romano, P.: An online repository of Swift/XRT light curves of γ -ray bursts. Astron. Astrophys. 469 (1), 379-385 (2007) https: //doi.org/10.1051/0004-6361:20077530 arXiv:0704.0128 [astro-ph]\n- [94] Evans, P.A., Beardmore, A.P., Page, K.L., Osborne, J.P., O'Brien, P.T., Willingale, R., Starling, R.L.C., Burrows, D.N., Godet, O., Vetere, L., Racusin, J., Goad, M.R., Wiersema, K., Angelini, L., Capalbi, M., Chincarini, G., Gehrels, N., Kennea, J.A., Margutti, R., Morris, D.C., Mountford, C.J., Pagani, C., Perri, M., Romano, P., Tanvir, N.: Methods and results of an automatic analysis of a complete sample of Swift-XRT observations of GRBs. Mon. Not. R. Astron. Soc. 397 (3), 1177-1201 (2009) https://doi.org/10.1111/j.1365-2966.2009.14913. x arXiv:0812.3662 [astro-ph]\n- [95] Kraft, R.P., Burrows, D.N., Nousek, J.A.: Determination of Confidence Limits for Experiments with Low Numbers of Counts. Astrophys. J. 374 , 344 (1991) https://doi.org/10.1086/170124\n- [96] Willingale, R., Starling, R.L.C., Beardmore, A.P., Tanvir, N.R., O'Brien, P.T.: Calibration of X-ray absorption in our Galaxy. Mon. Not. R. Astron. Soc. 431 (1), 394-404 (2013) https://doi.org/10.1093/mnras/stt175 arXiv:1303.0843 [astro-ph.HE]\n- [97] Levan, A.J., Lamb, G.P., Schneider, B., Hjorth, J., Zafar, T., de Ugarte Postigo, A., Sargent, B., Mullally, S.E., Izzo, L., D'Avanzo, P., Burns, E., Agu'ı Fern'andez, J.F., Barclay, T., Bernardini, M.G., Bhirombhakdi, K., Bremer, M., Brivio, R., Campana, S., Chrimes, A.A., D'Elia, V., De Pasquale, M., Ferro, M., Fong, W., Fruchter, A.S., Fynbo, J.P.U., Gaspari, N., Gompertz, B.P., Hartmann, D.H., Hedges, C.L., Heintz, K.E., Hotokezaka, K., Jakobsson, P., Kann, D.A., Kennea, J.A., Laskar, T., Le Floc'h, E., Malesani, D.B., Melandri, A., Metzger, B.D., Oates, S.R., Pian, E., Piranomonte, S., Pugliese, G., Racusin, J.L., Rastinejad, J.C., Ravasio, M.E., Rossi, A., Saccardi, A., Salvaterra, R., Sbarufatti, B., Starling, R.L.C., Tanvir, N.R., Thone, C.C., Vergani, S.D., Watson, D., Wiersema, K., Xu, D.: The first JWST spectrum of a GRB afterglow: No bright supernova in observations of the brightest GRB of all time, GRB 221009A. arXiv e-prints, 2302-07761 (2023) arXiv:2302.07761 [astro-ph.HE]\n- [98] Sault, R.J., Teuben, P.J., Wright, M.C.H.: A Retrospective View of MIRIAD. In: Shaw, R.A., Payne, H.E., Hayes, J.J.E. (eds.) Astronomical Data Analysis Software and Systems IV. Astronomical Society of the Pacific Conference Series, vol. 77, p. 433 (1995)\n- [99] Jonas, J., MeerKAT Team: The MeerKAT Radio Telescope. In: MeerKAT Science: On the Pathway to the SKA, p. 1 (2016). https://doi.org/10.22323/1.277. 0001\n- [100] Heywood, I.: oxkat: Semi-automated imaging of MeerKAT observations. Astrophysics Source Code Library, record ascl:2009.003 (2020)  \n- [101] Beichman, C.A., Rieke, M., Eisenstein, D., Greene, T.P., Krist, J., McCarthy, D., Meyer, M., Stansberry, J.: Science opportunities with the near-IR camera (NIRCam) on the James Webb Space Telescope (JWST). In: Clampin, M.C., Fazio, G.G., MacEwen, H.A., Oschmann, J. Jacobus M. (eds.) Space Telescopes and Instrumentation 2012: Optical, Infrared, and Millimeter Wave. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 8442, p. 84422 (2012). https://doi.org/10.1117/12.925447\n- [102] Jakobsen, P., Ferruit, P., Alves de Oliveira, C., Arribas, S., Bagnasco, G., Barho, R., Beck, T.L., Birkmann, S., Boker, T., Bunker, A.J., Charlot, S., de Jong, P., de Marchi, G., Ehrenwinkler, R., Falcolini, M., Fels, R., Franx, M., Franz, D., Funke, M., Giardino, G., Gnata, X., Holota, W., Honnen, K., Jensen, P.L., Jentsch, M., Johnson, T., Jollet, D., Karl, H., Kling, G., Kohler, J., Kolm, M.G., Kumari, N., Lander, M.E., Lemke, R., L'opez-Caniego, M., Lutzgendorf, N., Maiolino, R., Manjavacas, E., Marston, A., Maschmann, M., Maurer, R., Messerschmidt, B., Moseley, S.H., Mosner, P., Mott, D.B., Muzerolle, J., Pirzkal, N., Pittet, J.-F., Plitzke, A., Posselt, W., Rapp, B., Rauscher, B.J., Rawle, T., Rix, H.-W., Rodel, A., Rumler, P., Sabbi, E., Salvignol, J.-C., Schmid, T., Sirianni, M., Smith, C., Strada, P., te Plate, M., Valenti, J., Wettemann, T., Wiehe, T., Wiesmayer, M., Willott, C.J., Wright, R., Zeidler, P., Zincke, C.: The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope. I. Overview of the instrument and its capabilities. Astron. Astrophys. 661 , 80 (2022) https://doi.org/10.1051/0004-6361/202142663 arXiv:2202.03305 [astroph.IM]\n- [103] Perrin, M.D., Sivaramakrishnan, A., Lajoie, C.-P., Elliott, E., Pueyo, L., Ravindranath, S., Albert, L.: Updated point spread function simulations for JWST with WebbPSF. In: Oschmann, J. Jacobus M., Clampin, M., Fazio, G.G., MacEwen, H.A. (eds.) Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave. Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, vol. 9143, p. 91433 (2014). https://doi.org/10. 1117/12.2056689\n- [104] Mangan, J., Dunwoody, R., Meegan, C., Fermi GBM Team: GRB 211211A: Fermi GBM observation. GRB Coordinates Network 31210 , 1 (2021)  \n[105] von Kienlin, A., Meegan, C.A., Paciesas, W.S., Bhat, P.N., Bissaldi, E., Briggs, M.S., Burns, E., Cleveland, W.H., Gibby, M.H., Giles, M.M., Goldstein, A., Hamburg, R., Hui, C.M., Kocevski, D., Mailyan, B., Malacaria, C., Poolakkil, S., Preece, R.D., Roberts, O.J., Veres, P., Wilson-Hodge, C.A.: The Fourth FermiGBM Gamma-Ray Burst Catalog: A Decade of Data. Astrophys. J. 893 (1), 46 (2020) https://doi.org/10.3847/1538-4357/ab7a18 arXiv:2002.11460 [astroph.HE]  \n[106] Gehrels, N., Ramirez-Ruiz, E., Fox, D.B.: Gamma-Ray Bursts in the Swift Era. Annu. Rev. Astron. Astrophys. 47 (1), 567-617 (2009) https://doi.org/10.1146/", 'annurev.astro.46.060407.145147 arXiv:0909.1531 [astro-ph.HE]': "- [107] Nysewander, M., Fruchter, A.S., Pe'er, A.: A Comparison of the Afterglows of Short- and Long-duration Gamma-ray Bursts. Astrophys. J. 701 (1), 824-836 (2009) https://doi.org/10.1088/0004-637X/701/1/824 arXiv:0806.3607 [astroph]\n- [108] Evans, P.A., Willingale, R., Osborne, J.P., O'Brien, P.T., Page, K.L., Markwardt, C.B., Barthelmy, S.D., Beardmore, A.P., Burrows, D.N., Pagani, C., Starling, R.L.C., Gehrels, N., Romano, P.: The Swift Burst Analyser. I. BAT and XRT spectral and flux evolution of gamma ray bursts. Astron. Astrophys. 519 , 102 (2010) https://doi.org/10.1051/0004-6361/201014819 arXiv:1004.3208 [astro-ph.IM]\n- [109] Lien, A., Sakamoto, T., Barthelmy, S.D., Baumgartner, W.H., Cannizzo, J.K., Chen, K., Collins, N.R., Cummings, J.R., Gehrels, N., Krimm, H.A., Markwardt, C.B., Palmer, D.M., Stamatikos, M., Troja, E., Ukwatta, T.N.: The Third Swift Burst Alert Telescope Gamma-Ray Burst Catalog. Astrophys. J. 829 (1), 7 (2016) https://doi.org/10.3847/0004-637X/829/1/7 arXiv:1606.01956 [astro-ph.HE]\n- [110] Goldstein, A., Veres, P., Burns, E., Briggs, M.S., Hamburg, R., Kocevski, D., Wilson-Hodge, C.A., Preece, R.D., Poolakkil, S., Roberts, O.J., Hui, C.M., Connaughton, V., Racusin, J., von Kienlin, A., Dal Canton, T., Christensen, N., Littenberg, T., Siellez, K., Blackburn, L., Broida, J., Bissaldi, E., Cleveland, W.H., Gibby, M.H., Giles, M.M., Kippen, R.M., McBreen, S., McEnery, J., Meegan, C.A., Paciesas, W.S., Stanbro, M.: An Ordinary Short Gamma-Ray Burst with Extraordinary Implications: Fermi-GBM Detection of GRB 170817A. Astrophys. J. Lett. 848 (2), 14 (2017) https://doi.org/10.3847/2041-8213/aa8f41 arXiv:1710.05446 [astro-ph.HE]  \n[111] Hajela, A., Margutti, R., Bright, J.S., Alexander, K.D., Metzger, B.D., Nedora, V., Kathirgamaraju, A., Margalit, B., Radice, D., Guidorzi, C., Berger, E., MacFadyen, A., Giannios, D., Chornock, R., Heywood, I., Sironi, L., Gottlieb, O., Coppejans, D., Laskar, T., Cendes, Y., Duran, R.B., Eftekhari, T., Fong, W., McDowell, A., Nicholl, M., Xie, X., Zrake, J., Bernuzzi, S., Broekgaarden, F.S., Kilpatrick, C.D., Terreran, G., Villar, V.A., Blanchard, P.K., Gomez, S., Hosseinzadeh, G., Matthews, D.J., Rastinejad, J.C.: Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW170817: A New Emission Component. Astrophys. J. Lett. 927 (1), 17 (2022) https://doi.org/10.3847/2041-8213/ac504a arXiv:2104.02070 [astro-ph.HE]  \n[112] Schulze, S., Malesani, D., Cucchiara, A., Tanvir, N.R., Kruhler, T., de Ugarte Postigo, A., Leloudas, G., Lyman, J., Bersier, D., Wiersema, K., Perley, D.A., Schady, P., Gorosabel, J., Anderson, J.P., Castro-Tirado, A.J., Cenko, S.B., De Cia, A., Ellerbroek, L.E., Fynbo, J.P.U., Greiner, J., Hjorth, J., Kann, D.A.,  \nKaper, L., Klose, S., Levan, A.J., Mart'ın, S., O'Brien, P.T., Page, K.L., Pignata, G., Rapaport, S., S'anchez-Ram'ırez, R., Sollerman, J., Smith, I.A., Sparre, M., Thone, C.C., Watson, D.J., Xu, D., Bauer, F.E., Bayliss, M., Bjornsson, G., Bremer, M., Cano, Z., Covino, S., D'Elia, V., Frail, D.A., Geier, S., Goldoni, P., Hartoog, O.E., Jakobsson, P., Korhonen, H., Lee, K.Y., Milvang-Jensen, B., Nardini, M., Nicuesa Guelbenzu, A., Oguri, M., Pandey, S.B., Petitpas, G., Rossi, A., Sandberg, A., Schmidl, S., Tagliaferri, G., Tilanus, R.P.J., Winters, J.M., Wright, D., Wuyts, E.: GRB 120422A/SN 2012bz: Bridging the gap between low- and high-luminosity gamma-ray bursts. Astron. Astrophys. 566 , 102 (2014) https://doi.org/10.1051/0004-6361/201423387 arXiv:1401.3774 [astro-ph.HE]  \n[113] Bloom, J.S., Kulkarni, S.R., Djorgovski, S.G.: The Observed Offset Distribution of Gamma-Ray Bursts from Their Host Galaxies: A Robust Clue to the Nature of the Progenitors. Astron. J. 123 (3), 1111-1148 (2002) https://doi.org/10.1086/ 338893 arXiv:astro-ph/0010176 [astro-ph]  \n[114] Straatman, C.M.S., Labb'e, I., Spitler, L.R., Glazebrook, K., Tomczak, A., Allen, R., Brammer, G.B., Cowley, M., van Dokkum, P., Kacprzak, G.G., Kawinwanichakij, L., Mehrtens, N., Nanayakkara, T., Papovich, C., Persson, S.E., Quadri, R.F., Rees, G., Tilvi, V., Tran, K.-V.H., Whitaker, K.E.: The Sizes of Massive Quiescent and Star-forming Galaxies at z ∼ 4 with ZFOURGE and CANDELS. Astrophys. J. Lett. 808 (1), 29 (2015) https://doi.org/10.1088/ 2041-8205/808/1/L29 arXiv:1506.01380 [astro-ph.GA]  \n[115] Windhorst, R.A., Cohen, S.H., Jansen, R.A., Summers, J., Tompkins, S., Conselice, C.J., Driver, S.P., Yan, H., Coe, D., Frye, B., Grogin, N., Koekemoer, A., Marshall, M.A., O'Brien, R., Pirzkal, N., Robotham, A., Ryan, R.E., Willmer, C.N.A., Carleton, T., Diego, J.M., Keel, W.C., Porto, P., Redshaw, C., Scheller, S., Wilkins, S.M., Willner, S.P., Zitrin, A., Adams, N.J., Austin, D., Arendt, R.G., Beacom, J.F., Bhatawdekar, R.A., Bradley, L.D., Broadhurst, T., Cheng, C., Civano, F., Dai, L., Dole, H., D'Silva, J.C.J., Duncan, K.J., Fazio, G.G., Ferrami, G., Ferreira, L., Finkelstein, S.L., Furtak, L.J., Gim, H.B., Griffiths, A., Hammel, H.B., Harrington, K.C., Hathi, N.P., Holwerda, B.W., Honor, R., Huang, J.-S., Hyun, M., Im, M., Joshi, B.A., Kamieneski, P.S., Kelly, P., Larson, R.L., Li, J., Lim, J., Ma, Z., Maksym, P., Manzoni, G., Meena, A.K., Milam, S.N., Nonino, M., Pascale, M., Petric, A., Pierel, J.D.R., del Carmen Polletta, M., Rottgering, H.J.A., Rutkowski, M.J., Smail, I., Straughn, A.N., Strolger, L.-G., Swirbul, A., Trussler, J.A.A., Wang, L., Welch, B., B. Wyithe, J.S., Yun, M., Zackrisson, E., Zhang, J., Zhao, X.: JWST PEARLS. Prime Extragalactic Areas for Reionization and Lensing Science: Project Overview and First Results. Astron. J. 165 (1), 13 (2023) https://doi.org/10.3847/1538-3881/aca163 arXiv:2209.04119 [astro-ph.CO]  \n[116] Blanchard, P.K., Berger, E., Fong, W.-f.: The Offset and Host Light Distributions of Long Gamma-Ray Bursts: A New View From HST Observations of Swift Bursts. Astrophys. J. 817 (2), 144 (2016) https://doi.org/10.3847/0004-637X/", '817/2/144 arXiv:1509.07866 [astro-ph.HE]': "[117] Kruhler, T., Malesani, D., Fynbo, J.P.U., Hartoog, O.E., Hjorth, J., Jakobsson, P., Perley, D.A., Rossi, A., Schady, P., Schulze, S., Tanvir, N.R., Vergani, S.D., Wiersema, K., Afonso, P.M.J., Bolmer, J., Cano, Z., Covino, S., D'Elia, V., de Ugarte Postigo, A., Filgas, R., Friis, M., Graham, J.F., Greiner, J., Goldoni, P., Gomboc, A., Hammer, F., Japelj, J., Kann, D.A., Kaper, L., Klose, S., Levan, A.J., Leloudas, G., Milvang-Jensen, B., Nicuesa Guelbenzu, A., Palazzi, E., Pian, E., Piranomonte, S., S'anchez-Ram'ırez, R., Savaglio, S., Selsing, J., Tagliaferri, G., Vreeswijk, P.M., Watson, D.J., Xu, D.: GRB hosts through cosmic time. VLT/X-Shooter emission-line spectroscopy of 96 γ -ray-burst-selected galaxies at 0.1 ¡z ¡ 3.6. Astron. Astrophys. 581 , 125 (2015) https://doi.org/10.1051/ 0004-6361/201425561 arXiv:1505.06743 [astro-ph.GA]  \n[118] Martin, D.C., Fanson, J., Schiminovich, D., Morrissey, P., Friedman, P.G., Barlow, T.A., Conrow, T., Grange, R., Jelinsky, P.N., Milliard, B., Siegmund, O.H.W., Bianchi, L., Byun, Y.-I., Donas, J., Forster, K., Heckman, T.M., Lee, Y.-W., Madore, B.F., Malina, R.F., Neff, S.G., Rich, R.M., Small, T., Surber, F., Szalay, A.S., Welsh, B., Wyder, T.K.: The Galaxy Evolution Explorer: A Space Ultraviolet Survey Mission. Astrophys. J. 619 , 1 (2005) https://doi.org/ 10.1086/426387 astro-ph/0411302  \n[119] Dey, A., Schlegel, D.J., Lang, D., Blum, R., Burleigh, K., Fan, X., Findlay, J.R., Finkbeiner, D., Herrera, D., Juneau, S., Landriau, M., Levi, M., McGreer, I., Meisner, A., Myers, A.D., Moustakas, J., Nugent, P., Patej, A., Schlafly, E.F., Walker, A.R., Valdes, F., Weaver, B.A., Y'eche, C., Zou, H., Zhou, X., Abareshi, B., Abbott, T.M.C., Abolfathi, B., Aguilera, C., Alam, S., Allen, L., Alvarez, A., Annis, J., Ansarinejad, B., Aubert, M., Beechert, J., Bell, E.F., BenZvi, S.Y., Beutler, F., Bielby, R.M., Bolton, A.S., Brice˜no, C., BuckleyGeer, E.J., Butler, K., Calamida, A., Carlberg, R.G., Carter, P., Casas, R., Castander, F.J., Choi, Y., Comparat, J., Cukanovaite, E., Delubac, T., DeVries, K., Dey, S., Dhungana, G., Dickinson, M., Ding, Z., Donaldson, J.B., Duan, Y., Duckworth, C.J., Eftekharzadeh, S., Eisenstein, D.J., Etourneau, T., Fagrelius, P.A., Farihi, J., Fitzpatrick, M., Font-Ribera, A., Fulmer, L., Gansicke, B.T., Gaztanaga, E., George, K., Gerdes, D.W., Gontcho, S.G.A., Gorgoni, C., Green, G., Guy, J., Harmer, D., Hernand ez, M., Honscheid, K., Huang, L.W., James, D.J., Jannuzi, B.T., Jiang, L., Joyce, R., Karcher, A., Karkar, S., Kehoe, R., Kneib, J.-P., Kueter-Young, A., Lan, T.-W., Lauer, T.R., Le Guillou, L., Le Van Suu, A., Lee, J.H., Lesser, M., Perreault Levasseur, L., Li, T.S., Mann, J.L., Marshall, R., Mart'ınez-V'azquez, C.E., Martini, P., du Mas des Bourboux, H., McManus, S., Meier, T.G., M'enard, B., Metcalfe, N., Mu˜noz-Guti'errez, A., Najita, J., Napier, K., Narayan, G., Newman, J.A., Nie, J., Nord, B., Norman, D.J., Olsen, K.A.G., Paat, A., Palanque-Delabrouille, N., Peng, X., Poppett, C.L., Poremba, M.R., Prakash, A., Rabinowitz, D., Raichoor, A., Rezaie, M., Robertson, A.N., Roe, N.A., Ross, A.J., Ross, N.P., Rudnick, G., Safonova, S., Saha, A., S'anchez, F.J., Savary, E., Schweiker, H., Scott, A., Seo, H.-J.,  \nShan, H., Silva, D.R., Slepian, Z., Soto, C., Sprayberry, D., Staten, R., Stillman, C.M., Stupak, R.J., Summers, D.L., Sien Tie, S., Tirado, H., Vargas-Maga˜na, M., Vivas, A.K., Wechsler, R.H., Williams, D., Yang, J., Yang, Q., Yapici, T., Zaritsky, D., Zenteno, A., Zhang, K., Zhang, T., Zhou, R., Zhou, Z.: Overview of the DESI Legacy Imaging Surveys. Astron. J. 157 (5), 168 (2019) https://doi. org/10.3847/1538-3881/ab089d arXiv:1804.08657 [astro-ph.IM]  \n- [120] Wright, E.L., Eisenhardt, P.R.M., Mainzer, A.K., Ressler, M.E., Cutri, R.M., Jarrett, T., Kirkpatrick, J.D., Padgett, D., McMillan, R.S., Skrutskie, M., Stanford, S.A., Cohen, M., Walker, R.G., Mather, J.C., Leisawitz, D., Gautier, T.N. III, McLean, I., Benford, D., Lonsdale, C.J., Blain, A., Mendez, B., Irace, W.R., Duval, V., Liu, F., Royer, D., Heinrichsen, I., Howard, J., Shannon, M., Kendall, M., Walsh, A.L., Larsen, M., Cardon, J.G., Schick, S., Schwalm, M., Abid, M., Fabinsky, B., Naes, L., Tsai, C.-W.: The Wide-field Infrared Survey Explorer (WISE): Mission Description and Initial On-orbit Performance. Astron. J. 140 , 1868-1881 (2010) https://doi.org/10.1088/0004-6256/140/6/ 1868 arXiv:1008.0031 [astro-ph.IM]\n- [121] Lang, D.: unWISE: Unblurred Coadds of the WISE Imaging. Astron. J. 147 , 108 (2014) https://doi.org/10.1088/0004-6256/147/5/108 arXiv:1405.0308 [astro-ph.IM]  \n[122] Mainzer, A., Bauer, J., Cutri, R.M., Grav, T., Masiero, J., Beck, R., Clarkson, P., Conrow, T., Dailey, J., Eisenhardt, P., Fabinsky, B., Fajardo-Acosta, S., Fowler, J., Gelino, C., Grillmair, C., Heinrichsen, I., Kendall, M., Kirkpatrick, J.D., Liu, F., Masci, F., McCallon, H., Nugent, C.R., Papin, M., Rice, E., Royer, D., Ryan, T., Sevilla, P., Sonnett, S., Stevenson, R., Thompson, D.B., Wheelock, S., Wiemer, D., Wittman, M., Wright, E., Yan, L.: Initial Performance of the NEOWISE Reactivation Mission. Astrophys. J. 792 (1), 30 (2014) https://doi. org/10.1088/0004-637X/792/1/30 arXiv:1406.6025 [astro-ph.EP]  \n[123] Meisner, A.M., Lang, D., Schlegel, D.J.: Full-depth Coadds of the WISE and First-year NEOWISE-reactivation Images. Astron. J. 153 , 38 (2017) https:// doi.org/10.3847/1538-3881/153/1/38 arXiv:1603.05664 [astro-ph.IM]  \n[124] Wright, A.H., Robotham, A.S.G., Bourne, N., Driver, S.P., Dunne, L., Maddox, S.J., Alpaslan, M., Andrews, S.K., Bauer, A.E., Bland-Hawthorn, J., Brough, S., Brown, M.J.I., Clarke, C., Cluver, M., Davies, L.J.M., Grootes, M.W., Holwerda, B.W., Hopkins, A.M., Jarrett, T.H., Kafle, P.R., Lange, R., Liske, J., Loveday, J., Moffett, A.J., Norberg, P., Popescu, C.C., Smith, M., Taylor, E.N., Tuffs, R.J., Wang, L., Wilkins, S.M.: Galaxy And Mass Assembly: accurate panchromatic photometry from optical priors using LAMBDAR. Mon. Not. R. Astron. Soc. 460 , 765 (2016) https://doi.org/10.1093/mnras/stw832 arXiv:1604.01923  \n[125] Schulze, S., Yaron, O., Sollerman, J., Leloudas, G., Gal, A., Wright, A.H., Lunnan, R., Gal-Yam, A., Ofek, E.O., Perley, D.A., Filippenko, A.V., Kasliwal,  \n- M.M., Kulkarni, S.R., Neill, J.D., Nugent, P.E., Quimby, R.M., Sullivan, M., Strotjohann, N.L., Arcavi, I., Ben-Ami, S., Bianco, F., Bloom, J.S., De, K., Fraser, M., Fremling, C.U., Horesh, A., Johansson, J., Kelly, P.L., Kneˇzevi'c, N., Kneˇzevi'c, S., Maguire, K., Nyholm, A., Papadogiannakis, S., Petrushevska, T., Rubin, A., Yan, L., Yang, Y., Adams, S.M., Bufano, F., Clubb, K.I., Foley, R.J., Green, Y., Harmanen, J., Ho, A.Y.Q., Hook, I.M., Hosseinzadeh, G., Howell, D.A., Kong, A.K.H., Kotak, R., Matheson, T., McCully, C., Milisavljevic, D., Pan, Y.-C., Poznanski, D., Shivvers, I., van Velzen, S., Verbeek, K.K.: The Palomar Transient Factory Core-collapse Supernova Host-galaxy Sample. I. Host-galaxy Distribution Functions and Environment Dependence of Core-collapse Supernovae. Astrophys. J. Suppl. Ser. 255 (2), 29 (2021) https: //doi.org/10.3847/1538-4365/abff5e arXiv:2008.05988 [astro-ph.GA]\n- [126] Johnson, B.D., Leja, J., Conroy, C., Speagle, J.S.: Stellar Population Inference with Prospector. Astrophys. J. Suppl. Ser. 254 (2), 22 (2021) https://doi.org/ 10.3847/1538-4365/abef67 arXiv:2012.01426 [astro-ph.GA]\n- [127] Conroy, C., Gunn, J.E.: FSPS: Flexible Stellar Population Synthesis. Astrophysics Source Code Library, record ascl:1010.043 (2010)\n- [128] Calzetti, D., Armus, L., Bohlin, R.C., Kinney, A.L., Koornneef, J., StorchiBergmann, T.: The Dust Content and Opacity of Actively Star-forming Galaxies. Astrophys. J. 533 (2), 682-695 (2000) https://doi.org/10.1086/308692 arXiv:astro-ph/9911459 [astro-ph]\n- [129] Ferland, G.J., Chatzikos, M., Guzm'an, F., Lykins, M.L., van Hoof, P.A.M., Williams, R.J.R., Abel, N.P., Badnell, N.R., Keenan, F.P., Porter, R.L., Stancil, P.C.: The 2017 Release Cloudy. Rev. Mexicana Astron. Astrofis. 53 , 385-438 (2017) https://doi.org/10.48550/arXiv.1705.10877 arXiv:1705.10877 [astro-ph.GA]\n- [130] Speagle, J.S.: DYNESTY: a dynamic nested sampling package for estimating Bayesian posteriors and evidences. Mon. Not. R. Astron. Soc. 493 (3), 3132-3158 (2020) https://doi.org/10.1093/mnras/staa278 arXiv:1904.02180 [astro-ph.IM]\n- [131] Eldridge, J.J., Stanway, E.R., Xiao, L., McClelland, L.A.S., Taylor, G., Ng, M., Greis, S.M.L., Bray, J.C.: Binary Population and Spectral Synthesis Version 2.1: Construction, Observational Verification, and New Results. Publ. Astron. Soc. Aust. 34 , 058 (2017) https://doi.org/10.1017/pasa.2017.51 arXiv:1710.02154 [astro-ph.SR]\n- [132] Stanway, E.R., Eldridge, J.J.: Re-evaluating old stellar populations. Mon. Not. R. Astron. Soc. 479 (1), 75-93 (2018) https://doi.org/10.1093/mnras/sty1353 arXiv:1805.08784 [astro-ph.GA]\n- [133] Stevance, H., Eldridge, J., Stanway, E.: Hoki: Making BPASS accessible through Python. The Journal of Open Source Software 5 (45), 1987 (2020) https://doi.", 'org/10.21105/joss.01987 arXiv:2001.11069 [astro-ph.SR]': "- [134] Cappellari, M.: Improving the full spectrum fitting method: accurate convolution with Gauss-Hermite functions. Mon. Not. R. Astron. Soc. 466 (1), 798-811 (2017) https://doi.org/10.1093/mnras/stw3020 arXiv:1607.08538 [astro-ph.GA]\n- [135] Stevance, H.F., Eldridge, J.J., Stanway, E.R., Lyman, J., McLeod, A.F., Levan, A.J.: End-to-end study of the host galaxy and genealogy of the first binary neutron star merger. Nature Astronomy 7 , 444-450 (2023) https://doi.org/10. 1038/s41550-022-01873-y arXiv:2301.05236 [astro-ph.HE]\n- [136] Cid Fernandes, R., Mateus, A., Sodr'e, L., Stasi'nska, G., Gomes, J.M.: Semiempirical analysis of Sloan Digital Sky Survey galaxies - I. Spectral synthesis method. Mon. Not. R. Astron. Soc. 358 (2), 363-378 (2005) https://doi.org/10. 1111/j.1365-2966.2005.08752.x arXiv:astro-ph/0412481 [astro-ph]\n- [137] Bruzual, G., Charlot, S.: Stellar population synthesis at the resolution of 2003. Mon. Not. R. Astron. Soc. 344 (4), 1000-1028 (2003) https://doi.org/10.1046/j. 1365-8711.2003.06897.x arXiv:astro-ph/0309134 [astro-ph]\n- [138] Chabrier, G.: Galactic Stellar and Substellar Initial Mass Function. Publ. Astron. Soc. Pac. 115 (809), 763-795 (2003) https://doi.org/10.1086/376392 arXiv:astro-ph/0304382 [astro-ph]\n- [139] Kennicutt, J. Robert C.: The Global Schmidt Law in Star-forming Galaxies. Astrophys. J. 498 (2), 541-552 (1998) https://doi.org/10.1086/305588 arXiv:astro-ph/9712213 [astro-ph]\n- [140] Marino, R.A., Rosales-Ortega, F.F., S'anchez, S.F., Gil de Paz, A., V'ılchez, J., Miralles-Caballero, D., Kehrig, C., P'erez-Montero, E., Stanishev, V., IglesiasP'aramo, J., D'ıaz, A.I., Castillo-Morales, A., Kennicutt, R., L'opez-S'anchez, A.R., Galbany, L., Garc'ıa-Benito, R., Mast, D., Mendez-Abreu, J., MonrealIbero, A., Husemann, B., Walcher, C.J., Garc'ıa-Lorenzo, B., Masegosa, J., Del Olmo Orozco, A., Mour˜ao, A.M., Ziegler, B., Moll'a, M., Papaderos, P., S'anchez-Bl'azquez, P., Gonz'alez Delgado, R.M., Falc'on-Barroso, J., Roth, M.M., van de Ven, G., CALIFA Team: The O3N2 and N2 abundance indicators revisited: improved calibrations based on CALIFA and T e -based literature data. Astron. Astrophys. 559 , 114 (2013) https://doi.org/10.1051/0004-6361/ 201321956 arXiv:1307.5316 [astro-ph.CO]  \n[141] Fulton, M.D., Smartt, S.J., Rhodes, L., Huber, M.E., Villar, A.V., Moore, T., Srivastav, S., Schultz, A.S.B., Chambers, K.C., Izzo, L., Hjorth, J., Chen, T.-W., Nicholl, M., Foley, R.J., Rest, A., Smith, K.W., Young, D.R., Sim, S.A., Bright, J., Zenati, Y., de Boer, T., Bulger, J., Fairlamb, J., Gao, H., Lin, C.-C., Lowe, T., Magnier, E.A., Smith, I.A., Wainscoat, R., Coulter, D.A., Jones, D.O., Kilpatrick, C.D., McGill, P., Ramirez-Ruiz, E., Lee, K.S., Narayan, G., Ramakrishnan, V., Ridden-Harper, R., Singh, A., Wang, Q.,  \nKong, A.K.H., Ngeow, C.-C., Pan, Y.-C., Yang, S., Davis, K.W., Piro, A.L., Rojas-Bravo, C., Sommer, J., Yadavalli, S.K.: The optical light curve of GRB 221009A: the afterglow and detection of the emerging supernova SN 2022xiw. arXiv e-prints, 2301-11170 (2023) https://doi.org/10.48550/arXiv.2301.11170 arXiv:2301.11170 [astro-ph.HE]  \n- [142] Shrestha, M., Sand, D.J., Alexander, K.D., Bostroem, K.A., Hosseinzadeh, G., Pearson, J., Aghakhanloo, M., Vink'o, J., Andrews, J.E., Jencson, J.E., Lundquist, M.J., Wyatt, S., Howell, D.A., McCully, C., Padilla Gonzalez, E., Pellegrino, C., Terreran, G., Hiramatsu, D., Newsome, M., Farah, J., Jha, S.W., Smith, N., Wheeler, J.C., Mart'ınez-V'azquez, C., Carballo-Bello, J.A., DrlicaWagner, A., James, D.J., Mutlu-Pakdil, B., Stringfellow, G.S., Sakowska, J.D., Noel, N.E.D., Bom, C.R., Kuehn, K.: Lack of Bright Supernova Emission in the Brightest Gamma-ray Burst, GRB ˜ 221009A. arXiv e-prints, 2302-03829 (2023) https://doi.org/10.48550/arXiv.2302.03829 arXiv:2302.03829 [astro-ph.HE]\n- [143] Blanchard, P.K., Villar, V.A., Chornock, R., Sears, H., LeBaron, N., Yadavalli, S.K., Laskar, T., Alexander, K.D., Margutti, R., Berger, E., Barnes, J., Siegel, D., Metzger, B., Kasen, D., Cendes, Y., Eftekhari, T., Leja, J.: GRB 221009A: Late-time JWST Observations and Detection of Supernova Emission. GRB Coordinates Network 33676 , 1 (2023)\n- [144] Yin, Y.-H.I., Zhang, B.-B., Sun, H., Yang, J., Kang, Y., Shao, L., Yang, Y.H., Zhang, B.: GRB 211211A-like Events and How Gravitational Waves May Tell Their Origin. arXiv e-prints, 2304-06581 (2023) https://doi.org/10.48550/ arXiv.2304.06581 arXiv:2304.06581 [astro-ph.HE]\n- [145] Liang, E., Zhang, B., Virgili, F., Dai, Z.G.: Low-Luminosity Gamma-Ray Bursts as a Unique Population: Luminosity Function, Local Rate, and Beaming Factor. Astrophys. J. 662 (2), 1111-1118 (2007) https://doi.org/10.1086/517959 arXiv:astro-ph/0605200 [astro-ph]\n- [146] Salafia, O.S., Ravasio, M.E., Yang, J., An, T., Orienti, M., Ghirlanda, G., Nava, L., Giroletti, M., Mohan, P., Spinelli, R., Zhang, Y., Marcote, B., Cim'o, G., Wu, X., Li, Z.: Multiwavelength View of the Close-by GRB 190829A Sheds Light on Gamma-Ray Burst Physics. Astrophys. J. Lett. 931 (2), 19 (2022) https: //doi.org/10.3847/2041-8213/ac6c28 arXiv:2106.07169 [astro-ph.HE]\n- [147] Hotokezaka, K., Nakar, E.: Radioactive Heating Rate of r-process Elements and Macronova Light Curve. Astrophys. J. 891 (2), 152 (2020) https://doi.org/10. 3847/1538-4357/ab6a98\n- [148] Panaitescu, A., Kumar, P.: Analytic Light Curves of Gamma-Ray Burst Afterglows: Homogeneous versus Wind External Media. Astrophys. J. 543 (1), 66-76 (2000) https://doi.org/10.1086/317090\n- [149] Foreman-Mackey, D., Hogg, D.W., Lang, D., Goodman, J.: emcee: The mcmc  \n[167] Ohkubo, T., Nomoto, K., Umeda, H., Yoshida, N., Tsuruta, S.: Evolution of Very Massive Population III Stars with Mass Accretion from Pre-main Sequence to Collapse. Astrophys. J. 706 (2), 1184-1193 (2009) https://doi.org/10.1088/ 0004-637X/706/2/1184 arXiv:0902.4573 [astro-ph.SR]  \n[168] Perley, D.A., Mazzali, P.A., Yan, L., Cenko, S.B., Gezari, S., Taggart, K., Blagorodnova, N., Fremling, C., Mockler, B., Singh, A., Tominaga, N., Tanaka, M., Watson, A.M., Ahumada, T., Anupama, G.C., Ashall, C., Becerra, R.L., Bersier, D., Bhalerao, V., Bloom, J.S., Butler, N.R., Copperwheat, C., Coughlin, M.W., De, K., Drake, A.J., Duev, D.A., Frederick, S., Gonz'alez, J.J., Goobar, A., Heida, M., Ho, A.Y.Q., Horst, J., Hung, T., Itoh, R., Jencson, J.E., Kasliwal, M.M., Kawai, N., Khanam, T., Kulkarni, S.R., Kumar, B., Kumar, H., Kutyrev, A.S., Lee, W.H., Maeda, K., Mahabal, A., Murata, K.L., Neill, J.D., Ngeow, C.-C., Penprase, B., Pian, E., Quimby, R., Ramirez-Ruiz, E., Richer, M.G., Rom'an-Z'u˜niga, C.G., Sahu, D.K., Srivastav, S., Socia, Q., Sollerman, J., Tachibana, Y., Taddia, F., Tinyanont, S., Troja, E., Ward, C., Wee, J., Yu, P.-C.: The fast, luminous ultraviolet transient AT2018cow: extreme supernova, or disruption of a star by an intermediate-mass black hole? Mon. Not. R. Astron. Soc. 484 (1), 1031-1049 (2019) https://doi.org/10.1093/mnras/sty3420 arXiv:1808.00969 [astro-ph.HE]  \n[169] Ho, A.Y.Q., Perley, D.A., Gal-Yam, A., Lunnan, R., Sollerman, J., Schulze, S., Das, K.K., Dobie, D., Yao, Y., Fremling, C., Adams, S., Anand, S., Andreoni, I., Bellm, E.C., Bruch, R.J., Burdge, K.B., Castro-Tirado, A.J., Dahiwale, A., De, K., Dekany, R., Drake, A.J., Duev, D.A., Graham, M.J., Helou, G., Kaplan, D.L., Karambelkar, V., Kasliwal, M.M., Kool, E.C., Kulkarni, S.R., Mahabal, A.A., Medford, M.S., Miller, A.A., Nordin, J., Ofek, E., Petitpas, G., Riddle, R., Sharma, Y., Smith, R., Stewart, A.J., Taggart, K., Tartaglia, L., Tzanidakis, A., Winters, J.M.: A Search for Extragalactic Fast Blue Optical Transients in ZTF and the Rate of AT2018cow-like Transients. Astrophys. J. 949 (2), 120 (2023) https://doi.org/10.3847/1538-4357/acc533  \n[170] Andreoni, I., Coughlin, M.W., Perley, D.A., Yao, Y., Lu, W., Cenko, S.B., Kumar, H., Anand, S., Ho, A.Y.Q., Kasliwal, M.M., de Ugarte Postigo, A., Sagu'es-Carracedo, A., Schulze, S., Kann, D.A., Kulkarni, S.R., Sollerman, J., Tanvir, N., Rest, A., Izzo, L., Somalwar, J.J., Kaplan, D.L., Ahumada, T., Anupama, G.C., Auchettl, K., Barway, S., Bellm, E.C., Bhalerao, V., Bloom, J.S., Bremer, M., Bulla, M., Burns, E., Campana, S., Chandra, P., Charalampopoulos, P., Cooke, J., D'Elia, V., Das, K.K., Dobie, D., Agu'ı Fern'andez, J.F., Freeburn, J., Fremling, C., Gezari, S., Goode, S., Graham, M.J., Hammerstein, E., Karambelkar, V.R., Kilpatrick, C.D., Kool, E.C., Krips, M., Laher, R.R., Leloudas, G., Levan, A., Lundquist, M.J., Mahabal, A.A., Medford, M.S., Miller, M.C., Moller, A., Mooley, K.P., Nayana, A.J., Nir, G., Pang, P.T.H., Paraskeva, E., Perley, R.A., Petitpas, G., Pursiainen, M., Ravi, V., RiddenHarper, R., Riddle, R., Rigault, M., Rodriguez, A.C., Rusholme, B., Sharma, Y., Smith, I.A., Stein, R.D., Thone, C., Tohuvavohu, A., Valdes, F., van Roestel,  \n- J., Vergani, S.D., Wang, Q., Zhang, J.: A very luminous jet from the disruption of a star by a massive black hole. Nature 612 (7940), 430-434 (2022) https://doi.org/10.1038/s41586-022-05465-8 arXiv:2211.16530 [astro-ph.HE]\n- [171] Pasham, D.R., Lucchini, M., Laskar, T., Gompertz, B.P., Srivastav, S., Nicholl, M., Smartt, S.J., Miller-Jones, J.C.A., Alexander, K.D., Fender, R., Smith, G.P., Fulton, M., Dewangan, G., Gendreau, K., Coughlin, E.R., Rhodes, L., Horesh, A., van Velzen, S., Sfaradi, I., Guolo, M., Castro Segura, N., Aamer, A., Anderson, J.P., Arcavi, I., Brennan, S.J., Chambers, K., Charalampopoulos, P., Chen, T.-W., Clocchiatti, A., de Boer, T., Dennefeld, M., Ferrara, E., Galbany, L., Gao, H., Gillanders, J.H., Goodwin, A., Gromadzki, M., Huber, M., Jonker, P.G., Joshi, M., Kara, E., Killestein, T.L., Kosec, P., Kocevski, D., Leloudas, G., Lin, C.-C., Margutti, R., Mattila, S., Moore, T., Muller-Bravo, T., Ngeow, C.-C., Oates, S., Onori, F., Pan, Y.-C., Perez-Torres, M., Rani, P., Remillard, R., Ridley, E.J., Schulze, S., Sheng, X., Shingles, L., Smith, K.W., Steiner, J.F., Wainscoat, R., Wevers, T., Yang, S.: The Birth of a Relativistic Jet Following the Disruption of a Star by a Cosmological Black Hole. Nature Astronomy 7 , 88-104 (2023) https://doi.org/10.1038/s41550-022-01820-x arXiv:2211.16537 [astro-ph.HE]\n- [172] Waxman, E., Ofek, E.O., Kushnir, D.: Strong NIR emission following the long duration GRB 211211A: Dust heating as an alternative to a kilonova. arXiv e-prints, 2206-10710 (2022) https://doi.org/10.48550/arXiv.2206.10710 arXiv:2206.10710 [astro-ph.HE]\n- [173] Thomas, R., Pentericci, L., Le Fevre, O., Zamorani, G., Schaerer, D., Amorin, R., Castellano, M., Carnall, A.C., Cristiani, S., Cullen, F., Finkelstein, S.L., Fontanot, F., Guaita, L., Hibon, P., Hathi, N., Fynbo, J.P.U., Khusanova, Y., Koekemoer, A.M., McLeod, D., McLure, R.J., Marchi, F., Pozzetti, L., Saxena, A., Talia, M., Bolzonella, M.: The intergalactic medium transmission towards z ≳ 4 galaxies with VANDELS and the impact of dust attenuation. Astron. Astrophys. 634 , 110 (2020) https://doi.org/10.1051/0004-6361/ 201935925 arXiv:1911.12532 [astro-ph.CO]\n- [174] Nicholl, M., Guillochon, J., Berger, E.: The Magnetar Model for Type I Superluminous Supernovae. I. Bayesian Analysis of the Full Multicolor Light-curve Sample with MOSFiT. Astrophys. J. 850 (1), 55 (2017) https://doi.org/10.3847/ 1538-4357/aa9334 arXiv:1706.00825 [astro-ph.HE]\n- [175] Palmer, D.M., Barthelmy, S., Gehrels, N., Kippen, R.M., Cayton, T., Kouveliotou, C., Eichler, D., Wijers, R.A.M.J., Woods, P.M., Granot, J., Lyubarsky, Y.E., Ramirez-Ruiz, E., Barbier, L., Chester, M., Cummings, J., Fenimore, E.E., Finger, M.H., Gaensler, B.M., Hullinger, D., Krimm, H., Markwardt, C.B., Nousek, J.A., Parsons, A., Patel, S., Sakamoto, T., Sato, G., Suzuki, M., Tueller, J.: A giant γ -ray flare from the magnetar SGR 1806 - 20. Nature 434 (7037), 1107-1109 (2005) https://doi.org/10.1038/nature03525 arXiv:astro-ph/0503030"}