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.

First Constraints on the Photon Coupling of Axionlike Particles from Multimessenger Studies of the Neutron Star Merger GW170817

We use multimessenger observations of the neutron star merger event GW170817 to derive new constraints on axionlike particles (ALPs) coupling to photons. ALPs are produced via Primakoff and photon coalescence processes in the merger, escape the remnant, and decay back into two photons, giving rise to a photon signal approximately along the line of sight to the merger. We analyze the spectral and temporal information of the ALP-induced photon signal and use the Fermi Large Area Telescope (Fermi-LAT) observations of GW170817 to derive our new ALP constraints. We also show the improved prospects with future MeV γ-ray missions, taking the spectral and temporal coverage of Fermi-LAT as an example.

The space cold atom interferometer for testing the equivalence principle in the China Space Station

The precision of the weak equivalence principle (WEP) test using atom interferometers (AIs) is expected to be extremely high in microgravity environment. The microgravity scientific laboratory cabinet (MSLC) in the China Space Station (CSS) can provide a higher-level microgravity than the CSS itself, which provides a good experimental environment for scientific experiments that require high microgravity. We designed and realized a payload of a dual-species cold rubidium atom interferometer. The payload is highly integrated and has a size of [Formula: see text]. It will be installed in the MSLC to carry out high-precision WEP test experiment. In this article, we introduce the constraints and guidelines of the payload design, the compositions and functions of the scientific payload, the expected test precision in space, and some results of the ground test experiments.

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.

Reconnection-Powered Fast Radio Transients from Coalescing Neutron Star Binaries

It is an open question whether and how gravitational wave events involving neutron stars can be preceded by electromagnetic counterparts. This Letter shows that the collision of two neutron stars with magnetic fields well below magnetar-level strengths can produce millisecond fast-radio-burst-like transients. Using global force-free electrodynamics simulations, we identify the coherent emission mechanism that might operate in the common magnetosphere of a binary neutron star system prior to merger. We predict that the emission show have frequencies in the range of 10-20 GHz for magnetic fields of B^{*}=10^{11}  G at the surfaces of the stars.

Two Classes of Gamma-ray Bursts Distinguished within the First Second of Their Prompt Emission

Studies of Gamma-Ray Burst (GRB) properties, such as duration and spectral hardness, have found evidence for additional classes, beyond the short/hard and long/soft prototypes, using model-dependent methods. In this paper, a model-independent approach was used to analyse the gamma-ray light curves of large samples of GRBs detected by BATSE, Swift/BAT and Fermi/GBM. All the features were extracted from the GRB time profiles in four energy bands using the Stationary Wavelet Transform and Principal Component Analysis. t-distributed Stochastic Neighbourhood Embedding (t-SNE) visualisation of the features revealed two distinct groups of Swift/BAT bursts using the T100 interval with 64 ms resolution data. When the same analysis was applied to 4 ms resolution data, two groups were seen to emerge within the first second (T1) post-trigger. These two groups primarily consisted of short/hard (Group 1) and long/soft (Group 2) bursts, and were 95% consistent with the groups identified using the T100 64 ms resolution data. Kilonova candidates, arising from compact object mergers, were found to belong to Group 1, while those events with associated supernovae fell into Group 2. Differences in cumulative counts between the two groups in the first second, and in the minimum variability timescale, identifiable only with the 4 ms resolution data, may account for this result. Short GRBs have particular significance for multi-messenger science as a distinctive EM signature of a binary merger, which may be discovered by its gravitational wave emissions. Incorporating the T1 interval into classification algorithms may support the rapid classification of GRBs, allowing for an improved prioritisation of targets for follow-up observations.

Critical Tests of Leading Gamma Ray Burst Theories

It has been observationally established that supernovae (SNe) of Type Ic produce long duration gamma-ray bursts (GRBs) and that neutron star mergers generate short hard GRBs. SN-Less GRBs presumably originate in a phase transition of a neutron star in a high mass X-ray binary. How these phenomena actually generate GRBs is debated. The fireball and cannonball models of GRBs and their afterglows have been widely confronted with the huge observational data, with their defenders claiming success. The claims, however, may reflect multiple choices and the use of many adjustable parameters, rather than the validity of the models. Only a confrontation of key falsifiable predictions of the models with solid observational data can test their validity. Such critical tests are reviewed in this report.

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.

Breakthrough Multi-Messenger Astrophysics with the THESEUS Space Mission

The mission concept THESEUS (Transient High Energy Sky and Early Universe Surveyor) aims at exploiting Gamma-Ray Bursts (GRB) to explore the early Universe, as well as becoming a cornerstone of multi-messenger and time-domain astrophysics. To achieve these goals, a key feature is the capability to survey the soft X-ray transient sky and to detect the faint and soft GRB population so far poorly explored. Among the expected transients there will be high-redshift GRBs, nearby low-luminosity, X-ray Flashes and short GRBs. Our understanding of the physics governing the GRB prompt emission will benefit from the 0.3 keV–10 MeV simultaneous observations for an unprecedented large number of hundreds of events per year. In particular the mission will provide the identification, accurate sky localisation and characterization of electromagnetic counterparts to sources of gravitational wave and neutrino sources, which will be routinely detected during the 2030s by the upgraded second generation and third generation Gravitational Wave (GW) interferometers and next generation neutrino detectors.

A compact instrument for gamma-ray burst detection on a CubeSat platform II: Detailed design, assembly and validation

The Gamma-ray Module, GMOD, is a miniaturised novel gamma-ray detector which will be the primary scientific payload on the Educational Irish Research Satellite (EIRSAT-1) 2U CubeSat mission. GMOD comprises a compact (25 mm × 25 mm × 40 mm) cerium bromide scintillator coupled to a tiled array of 4 × 4 silicon photomultipliers, with front-end readout provided by the IDE3380 SIPHRA. This paper presents the detailed GMOD design and the accommodation of the instrument within the restrictive CubeSat form factor. The electronic and mechanical interfaces are compatible with many off-the-shelf CubeSat systems and structures. The energy response of the GMOD engineering qualification model has been determined using radioactive sources, and an energy resolution of 5.4% at 662 keV has been measured. EIRSAT-1 will perform on-board processing of GMOD data. Trigger results, including light-curves and spectra, will be incorporated into the spacecraft beacon and transmitted continuously. Inexpensive hardware can be used to decode the beacon signal, making the data accessible to a wide community. GMOD will have scientific capability for the detection of gamma-ray bursts, in addition to the educational and technology demonstration goals of the EIRSAT-1 mission. The detailed design and measurements to date demonstrate the capability of GMOD in low Earth orbit, the scalability of the design for larger CubeSats and as an element of future large gamma-ray missions.

PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems

The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.

New Relativistic Theory for Modified Newtonian Dynamics

We propose a relativistic gravitational theory leading to modified Newtonian dynamics, a paradigm that explains the observed universal galactic acceleration scale and related phenomenology. We discuss phenomenological requirements leading to its construction and demonstrate its agreement with the observed cosmic microwave background and matter power spectra on linear cosmological scales. We show that its action expanded to second order is free of ghost instabilities and discuss its possible embedding in a more fundamental theory.

Advanced Virgo: Status of the Detector, Latest Results and Future Prospects

The Virgo detector, based at the EGO (European Gravitational Observatory) and located in Cascina (Pisa), played a significant role in the development of the gravitational-wave astronomy. From its first scientific run in 2007, the Virgo detector has constantly been upgraded over the years; since 2017, with the Advanced Virgo project, the detector reached a high sensitivity that allowed the detection of several classes of sources and to investigate new physics. This work reports the main hardware upgrades of the detector and the main astrophysical results from the latest five years; future prospects for the Virgo detector are also presented.

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.

Dark Matters on the Scale of Galaxies

The cold dark-matter model successfully explains both the emergence and evolution of cosmic structures on large scales and, when we include a cosmological constant, the properties of the homogeneous and isotropic Universe. However, the cold dark-matter model faces persistent challenges on the scales of galaxies. Indeed, N-body simulations predict some galaxy properties that are at odds with the observations. These discrepancies are primarily related to the dark-matter distribution in the innermost regions of the halos of galaxies and to the dynamical properties of dwarf galaxies. They may have three different origins: (1) the baryonic physics affecting galaxy formation is still poorly understood and it is thus not properly included in the model; (2) the actual properties of dark matter differs from those of the conventional cold dark matter; (3) the theory of gravity departs from General Relativity. Solving these discrepancies is a rapidly evolving research field. We illustrate some of the solutions proposed within the cold dark-matter model, and solutions when including warm dark matter, self-interacting dark matter, axion-like particles, or fuzzy dark matter. We also illustrate some modifications of the theory of gravity: Modified Newtonian Dynamics (MOND), MOdified Gravity (MOG), and f(R) gravity.

Investigation of Infrasound Background Noise at Mátra Gravitational and Geophysical Laboratory (MGGL)

Infrasound and seismic waves are supposed to be the main contributors to the gravity-gradient noise (Newtonian noise) of the third-generation subterranean gravitational wave detectors. This noise will limit the sensitivity of the instrument at frequencies below 20 Hz. Investigation of its origin and the possible methods of mitigation have top priority during the designing period of the detectors. Therefore, long-term site characterizing measurements are needed at several subterranean sites. However, at some sites, mining activities can occur. These activities can cause sudden changes (transients) in the measured signal, and increase the continuous background noise, too. We have developed an algorithm based on discrete Haar transform to find these transients in the infrasound signal. We found that eliminating the transients decreases the variation of the noise spectra, and therefore results a more accurate characterization of the continuous background noise. We carried out experiments for controlling the continuous noise. Machines operating at the mine were turned on and off systematically in order to see their effect on the noise spectra. These experiments showed that the main contributor of the continuous noise is the ventilation system of the mine. We also estimated the contribution of infrasound Newtonian noise at MGGL to the strain noise of a subterranean GW detector similar to Einstein Telescope.

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.

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.

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.

Dark Energy after GW170817 Revisited

We revisit the status of scalar-tensor theories with applications to dark energy in the aftermath of the gravitational wave signal GW170817 and its optical counterpart GRB170817A. At the level of the cosmological background, we identify a class of theories, previously declared unviable in this context, whose anomalous gravitational wave speed is proportional to the scalar equation of motion. As long as the scalar field is assumed not to couple directly to matter, this raises the possibility of compatibility with the gravitational wave data, for any cosmological sources, thanks to the scalar dynamics. This newly "rescued" class of theories includes examples of generalized quintic Galileons from Horndeski theories. Despite the promise of this leading order result, we show that the loophole ultimately fails when we include the effect of large scale inhomogeneities.

Search for Neutrinos in Super-Kamiokande Associated with Gravitational Wave Events

We report the results from a search in Super-Kamiokande for neutrino signals coincident with gravitational-wave events using a neutrino energy range from 3.5 MeV–100 PeV. We searched for coincident neutrino events within a time window of ± 500 s around the gravitational-wave detection time. In this work, we report the number of events within the search-window and the 90 % confidence level upper limits on the neutrino fluence for each gravitational-wave event.

Astrophysical Sites that Can Produce Enantiomeric Amino Acids

Recent work has produced theoretical evidence for two sites, colliding neutron stars and neutron-star–Wolf–Rayet binary systems, which might produce amino acids with the left-handed chirality preference found in meteorites. The Supernova Neutrino Amino Acid Processing (SNAAP) model uses electron antineutrinos and the magnetic field from source objects such as neutron stars to preferentially destroy one enantiomer over another. Large enantiomeric excesses are predicted for isovaline and alanine; although based on an earlier study, similar results are expected for the others. Isotopic abundances of 13 C and 15 O in meteorites provide a new test of the SNAAP model. This presents implications for the origins of life.

Testing general relativity in cosmology

We review recent developments and results in testing general relativity (GR) at cosmological scales. The subject has witnessed rapid growth during the last two decades with the aim of addressing the question of cosmic acceleration and the dark energy associated with it. However, with the advent of precision cosmology, it has also become a well-motivated endeavor by itself to test gravitational physics at cosmic scales. We overview cosmological probes of gravity, formalisms and parameterizations for testing deviations from GR at cosmological scales, selected modified gravity (MG) theories, gravitational screening mechanisms, and computer codes developed for these tests. We then provide summaries of recent cosmological constraints on MG parameters and selected MG models. We supplement these cosmological constraints with a summary of implications from the recent binary neutron star merger event. Next, we summarize some results on MG parameter forecasts with and without astrophysical systematics that will dominate the uncertainties. The review aims at providing an overall picture of the subject and an entry point to students and researchers interested in joining the field. It can also serve as a quick reference to recent results and constraints on testing gravity at cosmological scales.