Swope Supernova Survey 2017a (SSS17a), the optical counterpart to a gravitational wave source

Low-latency gravitational wave alert products and their performance at the time of the fourth LIGO-Virgo-KAGRA observing run

Multimessenger searches for binary neutron star (BNS) and neutron star-black hole (NSBH) mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to gravitational wave (GW)s has resumed with ALIGO's, AdVirgo's and KAGRA's fourth observing run (O4). To support this effort, public semiautomated data products are sent in near real-time and include localization and source properties to guide complementary observations. In preparation for O4, we have conducted a study using a simulated population of compact binaries and a mock data challenge (MDC) in the form of a real-time replay to optimize and profile the software infrastructure and scientific deliverables. End-toend performance was tested, including data ingestion, running online search pipelines, performing annotations, and issuing alerts to the astrophysics community. We present an overview of the low-latency infrastructure and the performance of the data products that are now being released during O4 based on the MDC. We report the expected median latency for the preliminary alert of full bandwidth searches (29.5 s) and show consistency and accuracy of released data products using the MDC. We report the expected median latency for triggers from early warning searches (-3.1 s), which are new in O4 and target neutron star mergers during inspiral phase. This paper provides a performance overview for LIGO-Virgo-KAGRA (LVK) low-latency alert infrastructure and data products using theMDCand serves as a useful reference for the interpretation of O4 detections.

Independent evidence in multi-messenger astrophysics

In this paper I discuss the first "multi-messenger" observations of a binary neutron star merger and kilonova. These observations, touted as "revolutionary," included both gravitational-wave and electromagnetic observations of a single source. I draw on analogies between astrophysics and historical sciences (e.g., paleontology) to explain the significance of this for (gravitational-wave) astrophysics. In particular, I argue that having independent lines of evidence about a target system enables the use of argumentative strategies-the "Sherlock Holmes" method and consilience-that help overcome the key challenges astrophysics faces as an observational and historical science.

Heavy-element production in a compact object merger observed by JWST

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)1, sources of high-frequency gravitational waves (GWs)2 and likely production sites for heavy-element nucleosynthesis by means of rapid neutron capture (the r-process)3. Here we present observations of the exceptionally bright GRB 230307A. We show that GRB 230307A belongs to the class of long-duration GRBs associated with compact object mergers4-6 and contains a kilonova similar to AT2017gfo, associated with the GW merger GW170817 (refs. 7-12). We obtained James Webb Space Telescope (JWST) mid-infrared 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-infrared owing 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.

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M⊙ neutron star to be [Formula: see text] km.

Jets from Neutron-Star Merger Remnants and Massive Blue Kilonovae

We perform three-dimensional general-relativistic magnetohydrodynamic simulations with weak interactions of binary neutron-star (BNS) mergers resulting in a long-lived remnant neutron star, with properties typical of galactic BNS and consistent with those inferred for the first observed BNS merger GW170817. We demonstrate self-consistently that within ≲30  ms postmerger magnetized (σ∼5-10) incipient jets emerge with asymptotic Lorentz factor Γ∼5-10, which successfully break out from the merger debris within ≲20  ms. A fast (v≲0.6c), magnetized (σ∼0.1) wind surrounds the jet core and generates a UV/blue kilonova precursor on timescales of hours, similar to the precursor signal due to free neutron decay in fast dynamical ejecta. Postmerger ejecta are quickly dominated by magnetohydrodynamically driven outflows from an accretion disk. We demonstrate that, within only 50 ms postmerger, ≳2×10^{-2}M_{⊙} of lanthanide-free, quasispherical ejecta with velocities ∼0.1-0.2c is launched, yielding a kilonova signal consistent with GW170817 on timescales of ≲5  d.

Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85 million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017  s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Spherical symmetry in the kilonova AT2017gfo/GW170817

The mergers of neutron stars expel a heavy-element enriched fireball that can be observed as a kilonova^1-4. The kilonova's geometry is a key diagnostic of the merger and is dictated by the properties of ultra-dense matter and the energetics of the collapse to a black hole. Current hydrodynamical merger models typically show aspherical ejecta^5-7. Previously, Sr⁺ was identified in the spectrum⁸ of the only well-studied kilonova^9-11 AT2017gfo¹², associated with the gravitational wave event GW170817. Here we combine the strong Sr⁺ P Cygni absorption-emission spectral feature and the blackbody nature of kilonova spectrum to determine that the kilonova is highly spherical at early epochs. Line shape analysis combined with the known inclination angle of the source¹³ also show the same sphericity independently. We conclude that energy injection by radioactive decay is insufficient to make the ejecta spherical. A magnetar wind or jet from the black-hole disk could inject enough energy to induce a more spherical distribution in the overall ejecta; however, an additional process seems necessary to make the element distribution uniform.

A very luminous jet from the disruption of a star by a massive black hole

Tidal disruption events (TDEs) are bursts of electromagnetic energy that are released when supermassive black holes at the centres of galaxies violently disrupt a star that passes too close¹. TDEs provide a window through which to study accretion onto supermassive black holes; in some rare cases, this accretion leads to launching of a relativistic jet^2-9, but the necessary conditions are not fully understood. The best-studied jetted TDE so far is Swift J1644+57, which was discovered in γ-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical detection of AT2022cmc, a rapidly fading source at cosmological distance (redshift z = 1.19325) the unique light curve of which transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-ray, submillimetre and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron 'afterglow', probably launched by a supermassive black hole with spin greater than approximately 0.3. Using four years of Zwicky Transient Facility¹⁰ survey data, we calculate a rate of [Formula: see text] per gigapascals cubed per year for on-axis jetted TDEs on the basis of the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations¹¹. Correcting for the beaming angle effects, this rate confirms that approximately 1 per cent of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs.

A kilonova following a long-duration gamma-ray burst at 350 Mpc

Gamma-ray bursts (GRBs) are divided into two populations^1,2; long GRBs that derive from the core collapse of massive stars (for example, ref. ³) and short GRBs that form in the merger of two compact objects^4,5. Although it is common to divide the two populations at a gamma-ray duration of 2 s, classification based on duration does not always map to the progenitor. Notably, GRBs with short (≲2 s) spikes of prompt gamma-ray emission followed by prolonged, spectrally softer extended emission (EE-SGRBs) have been suggested to arise from compact object mergers^6-8. Compact object mergers are of great astrophysical importance as the only confirmed site of rapid neutron capture (r-process) nucleosynthesis, observed in the form of so-called kilonovae^9-14. Here we report the discovery of a possible kilonova associated with the nearby (350 Mpc), minute-duration GRB 211211A. The kilonova implies that the progenitor is a compact object merger, suggesting that GRBs with long, complex light curves can be spawned from merger events. The kilonova of GRB 211211A has a similar luminosity, duration and colour to that which accompanied the gravitational wave (GW)-detected binary neutron star (BNS) merger GW170817 (ref. ⁴). Further searches for GW signals coincident with long GRBs are a promising route for future multi-messenger astronomy.

Eccentricity of Long Inspiraling Compact Binaries Sheds Light on Dark Sirens

The localization and distance inference of gravitational waves are two crucial factors for dark sirens as precise probes of cosmology, astrophysics, and fundamental physics. In this Letter, for the first time we investigate the parameter estimation of gravitational waves emitted by the eccentric compact binaries in the midfrequency (0.1-10 Hz) band. Based on the configuration of one cluster of DECIGO (B-DECIGO), we simulate five types of typical compact binaries in GWTC-3 with component mass ranging from O(1∼100)  M_{⊙}. For each type of binaries, we assign discrete eccentricities from 0 to 0.4 at 0.1 Hz in 10^{3} random orientations. The multiple harmonics induced by eccentricity can break the degeneracy between parameters. We find that with eccentricity e_{0}=0.4, these typical binaries can achieve O(10^{2}-10^{4}) improvement for the distance inference in the near face-on orientations, compared to the circular case. More importantly, a nonvanishing eccentricity (0.01-0.4) can significantly improve the source localization of the typical binary black holes, most by 1.5-3.5 orders of magnitude. Our result shows the remarkable significance of eccentricity for dark sirens in the midband as precise probes of the Universe.

Incorporating a Radiative Hydrodynamics Scheme in the Numerical-Relativity Code BAM

To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHBΛϕ. Overall, our test simulations have good agreement with those reported in the literature.

Testing Quantum Gravity in the Multi-Messenger Astronomy Era

Quantum gravity (QG) remains elusive despite almost century-long efforts to combine general relativity and quantum mechanics. All the approaches triggered and powered by purely theoretical considerations eventually failed with a prevailing feeling of a complete lack of guidance from the experimental side. Currently, however, this circumstance is beginning to change considerably. We have entered the era of multi-messenger astronomy. The electromagnetic window to the universe—so far the only one—has been tremendously enlarged in the energy range beyond gamma rays up to ultra-high-energy photons and has been complemented by other messengers: high-energy cosmic rays, cosmic neutrinos, and gravitational waves (GWs). This has created a unique environment in which to observationally constrain various phenomenological QG effects. In this paper, we focus on the LIV phenomenology manifested as energy-dependent time-of-flight delays and strong lensing time delays. We review results regarding time-of-flight delays obtained with GRBs. We also recall the idea of energy-dependent lensing time delays, which allow one to constrain LIV models independently of the intrinsic time delay. Lastly, we show how strongly a gravitationally lensed GW signal would place interesting constraints on the LIV.

Multi-Messenger Constraints on the Hubble Constant Through Combination of Gravitational Waves, Gamma-Ray Bursts and Kilonovae from Neutron Star Mergers

The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant H0 that does not rely on a cosmic distance ladder, nor assumes a specific cosmological model. By using gravitational waves as “standard sirens”, this approach holds promise to arbitrate the existing tension between the H0 value inferred from the cosmic microwave background and those obtained from local measurements. However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source led to a H0 value from GW170817 that was not precise enough to resolve the existing tension. In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on H0. We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.

Electromagnetic Precursors of Short Gamma-Ray Bursts as Counterparts of Gravitational Waves

Precursor emissions are found in some short gamma-ray bursts (SGRBs). In this paper, we review the theories and observations of the SGRB precursor and discuss its prospect as an electromagnetic counterpart of the gravitational wave event produced by neutron star (NS) mergers. The observed luminosity, spectrum, and duration of precursors are explained by the magnetospheric interaction model during the inspiral or the cocoon/jet shock breakout model during the jet propagation. In general, these two models predict that the precursor will be weaker than the main GRB, but will be of a larger opening angle, which makes it an advantageous gamma-ray counterpart for NS mergers in the local Universe, especially for NS - black hole mergers with very low mass ratios, in which the main GRBs are not expected. The joint observation of the precursor, SGRB, and gravitational wave will help to reveal the jet launch mechanism and post-merger remnant.

The first 5 years of gravitational-wave astrophysics

Gravitational waves are ripples in spacetime generated by the acceleration of astrophysical objects; a direct consequence of general relativity, they were first directly observed in 2015. Here, I review the first 5 years of gravitational-wave detections. More than 50 gravitational-wave events have been found, emitted by pairs of merging compact objects such as neutron stars and black holes. These signals yield insights into the formation of compact objects and their progenitor stars, enable stringent tests of general relativity, and constrain the behavior of matter at densities higher than that of an atomic nucleus. Mergers that emit both gravitational and electromagnetic waves probe the formation of short gamma-ray bursts and the nucleosynthesis of heavy elements, and they measure the local expansion rate of the Universe.

A Possible Kilonova Powered by Magnetic Wind from a Newborn Black Hole

The merger of binary neutron stars (NS–NS) as the progenitor of short gamma-ray bursts (GRBs) has been confirmed by the discovery of the association of the gravitational-wave (GW) event GW170817 with GRB 170817A. However, the merger product of binary NS remains an open question. An X-ray plateau followed by a steep decay (“internal plateau”) has been found in some short GRBs, implying that a supramassive magnetar operates as the merger remnant and then collapses into a newborn black hole (BH) at the end of the plateau. X-ray bump or second plateau following the “internal plateau” are considered as the expected signature from the fallback accretion onto this newborn BH through the Blandford–Znajek mechanism (BZ). At the same time, a nearly isotropic wind driven by the Blandford–Payne mechanism (BP) from the newborn BH’s disk can produce a bright kilonova. Therefore, the bright kilonova observation for a short GRB with “internal plateau” (and followed by X-ray bump or second plateau) provides further evidence for this scenario. In this paper, we find that GRB 160821B is a candidate of such a case, and the kilonova emission of GRB 160821B is possibly powered by the BP wind from a newborn BH. Future GW detection of GRB 160821B–like events may provide further support to this scenario, enable us to investigate the properties of the magnetar and the newborn BH, and constrain the equation of state of neutron stars.

Search for Gravitational Waves from High-Mass-Ratio Compact-Binary Mergers of Stellar Mass and Subsolar Mass Black Holes

We present the first search for gravitational waves from the coalescence of stellar mass and subsolar mass black holes with masses between 20-100  M_{⊙} and 0.01-1  M_{⊙}(10-10^{3}  M_{J}), respectively. The observation of a single subsolar mass black hole would establish the existence of primordial black holes and a possible component of dark matter. We search the ∼164 day of public LIGO data from 2015-2017 when LIGO-Hanford and LIGO-Livingston were simultaneously observing. We find no significant candidate gravitational-wave signals. Using this nondetection, we place a 90% upper limit on the rate of 30-0.01  M_{⊙} and 30-0.1  M_{⊙} mergers at <1.2×10^{6} and <1.6×10^{4}  Gpc^{-3}  yr^{-1}, respectively. If we consider binary formation through direct gravitational-wave braking, this kind of merger would be exceedingly rare if only the lighter black hole were primordial in origin (<10^{-4}  Gpc^{-3}  yr^{-1}). If both black holes are primordial in origin, we constrain the contribution of 1(0.1)M_{⊙} black holes to dark matter to <0.3(3)%.

Multimessenger constraints on the neutron-star equation of state and the Hubble constant

Observations of neutron-star mergers with distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus and to measure the expansion rate of the Universe as quantified by the Hubble constant. We performed a joint analysis of the gravitational-wave event GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave event GW190425, both originating from neutron-star mergers. We combined these with previous measurements of pulsars using x-ray and radio observations, and nuclear-theory computations using chiral effective field theory, to constrain the neutron-star equation of state. We found that the radius of a 1.4-solar mass neutron star is [Formula: see text] km at 90% confidence and the Hubble constant is [Formula: see text] at 1σ uncertainty.

Measuring the Hubble constant with a sample of kilonovae

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant (H₀). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for H₀ contraints. Here, we show measurement of H₀ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=73.{8}_{-5.8}^{+6.3}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=73.8−5.8+6.3kms−1Mpc−1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${H}_{0}=71.{2}_{-3.1}^{+3.2}\ {\rm{km}}\ {{\rm{s}}}^{-1}\ {{\rm{Mpc}}}^{-1}$$\end{document}H0=71.2−3.1+3.2kms−1Mpc−1 for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs.

Kilonovae observations can be used to out constraints on the Hubble constant (H0). Here, the authors show H0 measurements by combining light curves of four short gamma-ray burts with GW170817 are about a factor of 2-3 more precise than the standard-siren measurements using only gravitational-waves.

An evolutionary system of mineralogy. Part I: Stellar mineralogy (>13 to 4.6 Ga)

Minerals preserve records of the physical, chemical, and biological histories of their origins and subsequent alteration, and thus provide a vivid narrative of the evolution of Earth and other worlds through billions of years of cosmic history. Mineral properties, including trace and minor elements, ratios of isotopes, solid and fluid inclusions, external morphologies, and other idiosyncratic attributes, represent information that points to specific modes of formation and subsequent environmental histories-information essential to understanding the co-evolving geosphere and biosphere. This perspective suggests an opportunity to amplify the existing system of mineral classification, by which minerals are defined solely on idealized end-member chemical compositions and crystal structures. Here we present the first in a series of contributions to explore a complementary evolutionary system of mineralogy-a classification scheme that links mineral species to their paragenetic modes. The earliest stage of mineral evolution commenced with the appearance of the first crystals in the universe at >13 Ga and continues today in the expanding, cooling atmospheres of countless evolved stars, which host the high-temperature (T > 1000 K), low-pressure (P < 10-2 atm) condensation of refractory minerals and amorphous phases. Most stardust is thought to originate in three distinct processes in carbon- and/or oxygen-rich mineral-forming stars: (1) condensation in the cooling, expanding atmospheres of asymptotic giant branch stars; (2) during the catastrophic explosions of supernovae, most commonly core collapse (Type II) supernovae; and (3) classical novae explosions, the consequence of runaway fusion reactions at the surface of a binary white dwarf star. Each stellar environment imparts distinctive isotopic and trace element signatures to the micro- and nanoscale stardust grains that are recovered from meteorites and micrometeorites collected on Earth's surface, by atmospheric sampling, and from asteroids and comets. Although our understanding of the diverse mineral-forming environments of stars is as yet incomplete, we present a preliminary catalog of 41 distinct natural kinds of stellar minerals, representing 22 official International Mineralogical Association (IMA) mineral species, as well as 2 as yet unapproved crystalline phases and 3 kinds of non-crystalline condensed phases not codified by the IMA.

Prospects for observing and localizing gravitational-wave transients with Advanced LIGO, Advanced Virgo and KAGRA

We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is binary neutron star, neutron star-black hole, and binary black hole systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. During O3, the median localization volume (90% credible region) is expected to be on the order of 10 5 , 10 6 , 10 7 Mpc 3 for binary neutron star, neutron star-black hole, and binary black hole systems, respectively. The localization volume in O4 is expected to be about a factor two smaller than in O3. We predict a detection count of 1 - 1 + 12 ( 10 - 10 + 52 ) for binary neutron star mergers, of 0 - 0 + 19 ( 1 - 1 + 91 ) for neutron star-black hole mergers, and 17 - 11 + 22 ( 79 - 44 + 89 ) for binary black hole mergers in a one-calendar-year observing run of the HLV network during O3 (HLVK network during O4). We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

Constraining the Neutron Star Equation of State Using Multiband Independent Measurements of Radii and Tidal Deformabilities

Using a Bayesian approach, we combine measurements of neutron star macroscopic observables obtained by astrophysical and gravitational observations to derive joint constraints on the equation of state (EOS) of matter at supranuclear density. In our analysis, we use two sets of data: (i) the masses and tidal deformabilities measured in the binary neutron star event GW170817, detected by LIGO and Virgo, and (ii) the masses and stellar radii measured from observations of nuclear bursts in accreting low-mass x-ray binaries. Using a phenomenological parametrization of the equation of state, we compute the posterior probability distributions of the EOS parameters, using which we infer the posterior distribution for the radius and the mass of the two neutron stars of GW170817. The constraints we set on the radii are tighter than previous bounds.

Multimessenger Bayesian parameter inference of a binary neutron star merger

The combined detection of a binary neutron star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum – GW170817/AT2017gfo/GRB170817A – marks a breakthrough in the field of multimessenger astronomy. Between the plethora of modelling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, and GRB170817 to obtain truly multimessenger constraints on the tidal deformability $\tilde{\Lambda } \in [302,860]$, total binary mass M ∈ [2.722, 2.751] M⊙, the radius of a 1.4 solar mass neutron star $R \in [11.3,13.5] \,\,\rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of q ≤ 1.27, all at 90 per cent confidence. Our joint novel analysis uses new phenomenological descriptions of the dynamical ejecta, debris disc mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.

Binary Neutron Star (BNS) Merger: What We Learned from Relativistic Ejecta of GW/GRB 170817A

Gravitational Waves (GW) from coalescence of a Binary Neutron Star (BNS) and its accompanying short Gamma-Ray Burst (GRB) GW/GRB 170817A confirmed the presumed origin of these puzzling transients and opened up the way for relating properties of short GRBs to those of their progenitor stars and their surroundings. Here we review an extensive analysis of the prompt gamma-ray and late afterglows of this event. We show that a fraction of polar ejecta from the merger had been accelerated to ultra-relativistic speeds. This structured jet had an initial Lorentz factor of about 260 in our direction, which was O ( 10 ∘ ) from the jet’s axis, and was a few orders of magnitude less dense than in typical short GRBs. At the time of arrival to circum-burst material the ultra-relativistic jet had a close to Gaussian profile and a Lorentz factor ≳ 130 in its core. It had retained in some extent its internal collimation and coherence, but had extended laterally to create mildly relativistic lobes—a cocoon. Its external shocks on the far from center inhomogeneous circum-burst material and low density of colliding shells generated slowly rising afterglows, which peaked more than 100 days after the prompt gamma-ray. The circum-burst material was somehow correlated with the merger. As non-relativistic outflows or tidally ejected material during BNS merger could not have been arrived to the location of the external shocks before the relativistic jet, circum-burst material might have contained recently ejected materials from resumption of internal activities, faulting and mass loss due to deformation and breaking of stars crusts by tidal forces during latest stages of their inspiral but well before their merger. By comparing these findings with the results of relativistic Magneto-Hydro-Dynamics (MHD) simulations and observed gravitational waves we conclude that progenitor neutron stars were most probably old, had close masses and highly reduced magnetic fields.

Observing the Dark Sector

Despite the observational success of the standard model of cosmology, present-day observations do not tightly constrain the nature of dark matter and dark energy and modifications to the theory of general relativity. Here, we will discuss some of the ongoing and upcoming surveys that will revolutionize our understanding of the dark sector.

Collapsars as a major source of r-process elements

The production of elements by rapid neutron capture (r-process) in neutron-star mergers is expected theoretically and is supported by multimessenger observations^1-3 of gravitational-wave event GW170817: this production route is in principle sufficient to account for most of the r-process elements in the Universe⁴. Analysis of the kilonova that accompanied GW170817 identified^5,6 delayed outflows from a remnant accretion disk formed around the newly born black hole^7-10 as the dominant source of heavy r-process material from that event^9,11. Similar accretion disks are expected to form in collapsars (the supernova-triggering collapse of rapidly rotating massive stars), which have previously been speculated to produce r-process elements^12,13. Recent observations of stars rich in such elements in the dwarf galaxy Reticulum II¹⁴, as well as the Galactic chemical enrichment of europium relative to iron over longer timescales^15,16, are more consistent with rare supernovae acting at low stellar metallicities than with neutron-star mergers. Here we report simulations that show that collapsar accretion disks yield sufficient r-process elements to explain observed abundances in the Universe. Although these supernovae are rarer than neutron-star mergers, the larger amount of material ejected per event compensates for the lower rate of occurrence. We calculate that collapsars may supply more than 80 per cent of the r-process content of the Universe.

Compact radio emission indicates a structured jet was produced by a binary neutron star merger

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High-spatial-resolution measurements of the source size and displacement can discriminate between these scenarios. We present very-long-baseline interferometry observations, performed 207.4 days after the merger by using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milli-arc seconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

Fingerprints of Heavy-Element Nucleosynthesis in the Late-Time Lightcurves of Kilonovae

The kilonova emission observed following the binary neutron star merger event GW170817 provided the first direct evidence for the synthesis of heavy nuclei through the rapid neutron capture process (r process). The late-time transition in the spectral energy distribution to near-infrared wavelengths was interpreted as indicating the production of lanthanide nuclei, with atomic mass number A≳140. However, compelling evidence for the presence of even heavier third-peak (A≈195) r-process elements (e.g., gold, platinum) or translead nuclei remains elusive. At early times (∼days) most of the r-process heating arises from a large statistical ensemble of β decays, which thermalize efficiently while the ejecta is still dense, generating a heating rate that is reasonably approximated by a single power law. However, at later times of weeks to months, the decay energy input can also possibly be dominated by a discrete number of α decays, ^{223}Ra (half-life t_{1/2}=11.43  d), ^{225}Ac (t_{1/2}=10.0  d, following the β decay of ^{225}Ra with t_{1/2}=14.9  d), and the fissioning isotope ^{254}Cf (t_{1/2}=60.5  d), which liberate more energy per decay and thermalize with greater efficiency than β-decay products. Late-time nebular observations of kilonovae which constrain the radioactive power provide the potential to identify signatures of these individual isotopes, thus confirming the production of heavy nuclei. In order to constrain the bolometric light to the required accuracy, multiepoch and wideband observations are required with sensitive instruments like the James Webb Space Telescope. In addition, by comparing the nuclear heating rate obtained with an abundance distribution that follows the solar r abundance pattern, to the bolometric lightcurve of AT2017gfo, we find that the yet-uncertain r abundance of ^{72}Ge plays a decisive role in powering the lightcurve, if one assumes that GW170817 has produced a full range of the solar r abundances down to mass number A∼70.

High-energy emissions from neutron star mergers

In 2017, LIGO-Virgo collaborations reported detection of the first neutron star merger event, GW170817, which is accompanied by electromagnetic counterparts from radio to gamma rays. Although high-energy neutrinos were not detected from this event, mergers of neutron stars are expected to produce such high-energy particles. Relativistic jets are launched when neutron stars merge. If the jets contain protons, they can emit high-energy neutrinos through photomeson production. In addition, neutron star mergers produce massive and fast ejecta, which can be a source of Galactic high-energy cosmic rays above the knee. We briefly review what we learned from the multi-messenger event, GW170817, and discuss prospects for multi-messenger detections and hadronic cosmic-ray production related to the neutron star mergers.

Multi-messenger Astrophysics at Ultra-High Energy with the Pierre Auger Observatory

The study of correlations between observations of fundamentally different nature from extreme cosmic sources promises extraordinary physical insights into the Universe. With the Pierre Auger Observatory, we can significantly contribute to multi-messenger astrophysics by searching for ultra-high energy particles, particularly neutrinos and photons which, being electrically neutral, point back to their origin. Using Pierre Auger Observatory data, stringent limits at EeV energies have been established on the photon and neutrino fluxes from a large fraction of the sky, probing the production mechanisms of ultra-high energy cosmic rays. The good angular resolution and the neutrino identification capabilities of the Observatory at EeV energies allow the follow-up of events detected in gravitational waves, such as the binary mergers observed with the Advanced LIGO/Virgo detectors, or from other energetic sources of particles.