GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral

Gravitational wave astronomy with TianQin

The opening of the gravitational wave window has significantly enhanced our capacity to explore the Universe's most extreme and dynamic sector. In the mHz frequency range, a diverse range of compact objects, from the most massive black holes at the farthest reaches of the Universe to the lightest white dwarfs in our cosmic backyard, generate a complex and dynamic symphony of gravitational wave signals. Once recorded by gravitational wave detectors, these unique fingerprints have the potential to decipher the birth and growth of cosmic structures over a wide range of scales, from stellar binaries and stellar clusters to galaxies and large-scale structures. The TianQin space-borne gravitational wave mission is scheduled for launch in the 2030s, with an operational lifespan of five years. It will facilitate pivotal insights into the history of our Universe. This document presents a concise overview of the detectable sources of TianQin, outlining their characteristics, the challenges they present, and the expected impact of the TianQin observatory on our understanding of them.

Impact of Newly Measured β-Delayed Neutron Emitters around ^{78}Ni on Light Element Nucleosynthesis in the Neutrino Wind Following a Neutron Star Merger

Neutron emission probabilities and half-lives of 37 β-delayed neutron emitters from ^{75}Ni to ^{92}Br were measured at the RIKEN Nishina Center in Japan, including 11 one-neutron and 13 two-neutron emission probabilities and six half-lives for the first time that supersede theoretical estimates. These nuclei lie in the path of the weak r process occurring in neutrino-driven winds from the accretion disk formed after the merger of two neutron stars synthesizing elements in the A∼80 abundance peak. The presence of such elements dominates the accompanying kilonova emission over the first few days and have been identified in the AT2017gfo event, associated to the gravitational wave detection GW170817. Abundance calculations based on over 17 000 simulated trajectories describing the evolution of matter properties in the merger outflows show that the new data lead to an increase of 50%-70% in the abundance of Y, Zr, Nb, and Mo. This enhancement is large compared to the scatter of relative abundances observed in old very metal poor stars and thus is significant in the comparison with other possible astrophysical processes contributing to the light-element production. These results underline the importance of including experimental decay data for very neutron-rich β-delayed neutron emitters into r-process models.

Multi-messenger gravitational lensing

We introduce the rapidly emerging field of multi-messenger gravitational lensing-the discovery and science of gravitationally lensed phenomena in the distant universe through the combination of multiple messengers. This is framed by gravitational lensing phenomenology that has grown since the first discoveries in the twentieth century, messengers that span 30 orders of magnitude in energy from high-energy neutrinos to gravitational waves, and powerful 'survey facilities' that are capable of continually scanning the sky for transient and variable sources. Within this context, the main focus is on discoveries and science that are feasible in the next 5-10 years with current and imminent technology including the LIGO-Virgo-KAGRA network of gravitational wave detectors, the Vera C. Rubin Observatory and contemporaneous gamma/X-ray satellites and radio surveys. The scientific impact of even one multi-messenger gravitational lensing discovery will be transformational and reach across fundamental physics, cosmology and astrophysics. We describe these scientific opportunities and the key challenges along the path to achieving them. This article therefore describes the consensus that emerged at the eponymous Theo Murphy meeting in March 2024, and also serves as an introduction to this Theo Murphy meeting issue.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 2)'.

Kilonova simulations: connecting observations with the underlying physics

Kilonova observations contain information about heavy-element r-process nucleosynthesis and the behaviour of high-density matter. However, interpreting what these observations tell us about the underlying physics requires detailed modelling. We outline recent kilonova radiative transfer simulations that are based on hydrodynamical models of neutron star merger ejecta. The simulated spectra in the polar directions show a remarkably similar evolution to the observations of AT2017gfo. Using these simulations, we show the importance of accurate atomic data for kilonova modelling, as well as the importance of three-dimensional simulations. By improving radiative transfer simulations and by extending this study to consider a range of theoretical equations of state, simulations will be able to connect observations to the underlying merger physics and place constraints on the high-density equation of state and r-process nucleosynthesis.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 2)'.

False positives for gravitational lensing: the gravitational-wave perspective

For the first detection of a novel astrophysical phenomenon, scientific standards are particularly high. Especially in a multi-messenger context, there are also opportunity costs to follow-up observations on any detection claims. So in searching for the still-elusive lensed gravitational waves (GWs), care needs to be taken in controlling false positives. In particular, many methods for identifying strong lensing rely on some form of parameter similarity or waveform consistency, which under rapidly growing catalogue sizes can expose them to false positives from coincident but unlensed events if proper care is not taken. Searches for waveform deformations in all lensing regimes are subject to degeneracies; we need to mitigate between lensing, intrinsic parameters, insufficiently modelled effects such as orbital eccentricity, or even deviations from general relativity. Robust lensing studies also require understanding and mitigation of glitches and non-stationarities in the detector data. This article reviews sources of possible false positives (and their flip side: false negatives) in GW lensing searches and the main approaches the community is pursuing to mitigate them.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 2)'.

Gravitational lensing: towards combining the multi-messengers

The next generation of gravitational wave (GW) detectors and electromagnetic telescopes are beckoning the onset of the multi-messenger era and the exciting science that lies ahead. Multi-messenger strong gravitational lensing will help probe some of the most important questions of the Universe in an unprecedented manner. In particular, understanding the nature of GW sources, the underlying physical processes and mechanisms that produce emissions well before or right until the time of the merger, their associations to the seemingly distinct populations of gamma-ray bursts, fast radio bursts and kilonovae. Not to mention the fact that multi-messenger lensing will offer unique probes of test of gravity models and constraints on cosmological parameters complementary to other probes. Enabling multi-messenger science calls for concerted follow-up efforts and the development of new and shared resources in the community.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 2)'.

Accessing Universal Relations of Binary Neutron Star Waveforms in Massive Scalar-Tensor Theory

We investigate how the quasiuniversal relations connecting tidal deformability with gravitational waveform characteristics and/or properties of individual neutron stars that were proposed in the literature within general relativity would be influenced in the massive Damour-Esposito-Farese-type scalar-tensor gravity. For this purpose, we systematically perform numerical relativity simulations of ∼120 binary neutron-star mergers with varying scalar coupling constants. Although only three neutron-star equations of state are adopted, a clear breach of universality can be observed in the datasets. In addition to presenting difficulties in constructing quasiuniversal relations in alternative gravity theories, we also briefly compare the impacts of non-general-relativity physics on the waveform features and those due to the first order or cross-over quantum chromodynamical phase transition.

Electromagnetic follow-up of gravitational waves: review and lessons learned

The detection of gravitational waves (GWs) has provided a new tool to study the Universe, with the scientific return enriched when combined with established probes: electromagnetic (EM) radiation and energetic particles. Since the groundbreaking detection in 2017 of merging neutron stars producing GW emission, a gamma-ray burst and an optical 'kilonova', the field has grown rapidly. At present, no additional neutron star mergers have been jointly detected in GW and EM radiation, but with upgrades in EM and GW facilities now is a chance to take stock of almost a decade of observations. We discuss the motivations for following up GW sources and the basic challenges of searching large areas for a rapidly evolving EM signal. We examine how the kilonova counterpart to GW170817 was discovered and the association confirmed, and outline some of the key physics enabled by this discovery. We then review the status of EM searches since 2017, highlighting areas where more information (in GW alerts or catalogs) can improve efficiency, and discuss what we have learned about kilonovae despite the lack of further multi-messenger detections. We discuss upcoming facilities and the many lessons learned, considering also how these could inform searches for lensed mergers.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 1)'.

Finding black holes: an unconventional multi-messenger

A rather clear problem has remained for astrophysicists studying black holes: localizing black holes. One of the recent theoretical ways proposed to identify black hole mergers' hosts is through multi-messenger gravitational lensing: matching the properties of a lensed host galaxy with those of a lensed gravitational wave (GW). This paper reviews the most recent literature and introduces some of the ongoing work on the localization of binary black holes (BBHs) and their host galaxies through lensing of GWs and their electromagnetically bright hosts.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 1)'.

Gravitational lensing in gamma-ray bursts

Gravitationally lensed gamma-ray bursts (GRBs) offer critical advantages over other lensed sources. They can be detected via continuously operating detectors covering most of the sky. They offer extremely high time resolution to determine lensing delays and find short time delays accurately. They are detectable across most of the visible Universe. However, they are also rare and frequently poorly localized. In this article, we review searches for gravitational lensing in GRBs and comment on promising avenues for the future. We note that the highly structured jets in GRBs can show variations on sufficiently small scales that, unlike lensing of most transient sources, gravitational lensing in GRBs may not be achromatic. Such behaviour would weaken the stringent requirements for identifying lensed bursts but would also make robust identification of lensing more challenging. A continuously running search that could identify candidate lensed events in near real-time would enable afterglow searches with current and near-future wide-field optical/infrared surveys that could yield the first unambiguous detection of a lensed GRB. The new generation of sensitive X-ray and γ-ray detectors, such as the Einstein Probe and SVOM, will complement Swift and significantly enhance the number of well-localized γ-ray and X-ray transients. Tuned strategies could dramatically improve the probability of observing a lensed GRB.This article is part of the Theo Murphy meeting issue 'Multi-messenger gravitational lensing (Part 1)'.

Conservative Spin-Magnitude Change in Orbital Evolution in General Relativity

We show that physical scattering observables for compact spinning objects in general relativity can depend on additional degrees of freedom in the spin tensor beyond those described by the spin vector alone. The impulse, spin kick, and leading-order waveforms exhibit such a nontrivial dependence. A signal of this additional structure is the change in the magnitude of the spin vector under conservative Hamiltonian evolution, similar to our previous studies in electrodynamics. These additional degrees of freedom describe dynamical mass multipoles of compact objects and decouple for black holes. We also show that the conservative impulse, spin kick, and change of the additional degrees of freedom are encoded in the eikonal phase.

Real-time inference for binary neutron star mergers using machine learning

Mergers of binary neutron stars emit signals in both the gravitational-wave (GW) and electromagnetic spectra. Famously, the 2017 multi-messenger observation of GW170817 (refs. 1,2) led to scientific discoveries across cosmology3, nuclear physics4-6 and gravity7. Central to these results were the sky localization and distance obtained from the GW data, which, in the case of GW170817, helped to identify the associated electromagnetic transient, AT 2017gfo (ref. 8), 11 h after the GW signal. Fast analysis of GW data is critical for directing time-sensitive electromagnetic observations. However, owing to challenges arising from the length and complexity of signals, it is often necessary to make approximations that sacrifice accuracy. Here we present a machine-learning framework that performs complete binary neutron star inference in just 1 s without making any such approximations. Our approach enhances multi-messenger observations by providing: (1) accurate localization even before the merger; (2) improved localization precision by around 30% compared to approximate low-latency methods; and (3) detailed information on luminosity distance, inclination and masses, which can be used to prioritize expensive telescope time. Additionally, the flexibility and reduced cost of our method open new opportunities for equation-of-state studies. Finally, we demonstrate that our method scales to long signals, up to an hour in length, thus serving as a blueprint for data analysis for next-generation ground- and space-based detectors.

Impact of Bulk Viscosity on the Postmerger Gravitational-Wave Signal from Merging Neutron Stars

In the violent postmerger of binary neutron-star mergers strong oscillations are present that impact the emitted gravitational-wave (GW) signal. The frequencies, temperatures, and densities involved in these oscillations allow for violations of the chemical equilibrium promoted by weak interactions, thus leading to a nonzero bulk viscosity that can impact dynamics and GW signals. We present the first simulations of binary neutron-star mergers employing the self-consistent and second-order formulation of the equations of relativistic hydrodynamics for dissipative fluids proposed by Müller, Israel, and Stewart. With the spirit of obtaining a first assessment of the impact of bulk viscosity on the structure and radiative efficiency of the merger remnant we adopt a simplified but realistic approach for the viscosity, which we assume to be determined by direct and modified Urca reactions and hence to vary within the stars. At the same time, to compensate for the lack of a precise knowledge about the strength of bulk viscosity, we explore the possible behaviors by considering three different scenarios of low, medium, and high bulk viscosity. In this way, we find that large values of the bulk viscosities damp the collision-and-bounce oscillations that characterize the dynamics of the stellar cores right after the merger. At the same time, large viscosities tend to preserve the m=2 deformations in the remnant, thus leading to a comparatively more efficient GW emission and to changes in the postmerger spectrum that can be up to 100 Hz in the case of the most extreme configurations. Overall, our self-consistent results indicate that bulk viscosity increases the energy radiated in GWs soon after the merger by ≲2% in the (realistic) scenario of small viscosity, and by ≲30% in the (unrealistic) scenario of large viscosity.

Direct Mass Measurements of Neutron-Rich Zinc and Gallium Isotopes: An Investigation of the Formation of the First r-Process Peak

The prediction of isotopic abundances resulting from the rapid neutron capture process (r process) requires high-precision mass measurements. Using TITAN's on-line time-of-flight spectrometer, first time mass measurements are performed for ^{83}Zn and ^{86}Ga. These measurements reduced uncertainties, and are used to calculate isotopic abundances near the first r-process abundance peak using astrophysical conditions present during a binary neutron star (BNS) merger. Good agreement in abundance across a range of trajectories is found when comparing to several metal-poor stars while also strongly deviating from the solar r-process pattern. These findings point to a high degree of sensitivity to the electron fraction of a BNS merger on the final elemental abundance pattern for certain elements near the first r-process peak while others display universality. We find that small changes in electron fraction can produce distinct abundance patterns that match those of metal-poor stars with different classifications.

Constraining the equation of state in neutron-star cores via the long-ringdown signal

Multimessenger signals from binary neutron star (BNS) mergers are promising tools to infer the properties of nuclear matter at densities inaccessible to laboratory experiments. Gravitational waves (GWs) from BNS merger remnants can constrain the neutron-star equation of state (EOS) complementing constraints from late inspiral, direct mass-radius measurements, and ab-initio calculations. We perform a series of general-relativistic simulations of BNS systems with EOSs constructed to comprehensively cover the high-density regime. We identify a tight correlation between the ratio of the energy and angular-momentum losses in the late-time portion of the post-merger signal, called the long ringdown, and the EOS at the highest pressures and densities in neutron-star cores. Applying this correlation to post-merger GW signals significantly reduces EOS uncertainty at densities several times the nuclear saturation density, where no direct constraints are currently available. Hence, the long ringdown can provide stringent constraints on material properties of neutron stars cores.

Testing Mirror Symmetry in the Universe with LIGO-Virgo Black-Hole Mergers

Certain precessing black-hole mergers produce gravitational waves with net circular polarization, understood as an imbalance between right- and left-handed amplitudes. According to the cosmological principle, such emission must average to zero across all binary mergers in our Universe to preserve mirror-reflection symmetry at very large scales. We present a new independent gravitational-wave test of this hypothesis. Using a novel observable based on the Chern-Pontryagin pseudoscalar, we measure the emission of net circular polarization across 47 black-hole mergers recently analyzed by [T. Islam et al., arXiv:2309.14473.] with a state-of-the art model for precessing black-hole mergers in general relativity. The average value obtained is consistent with zero. Remarkably, however, we find that at least 82% of the analyzed sources must have produced net circular polarization. Of these, GW200129 shows strong evidence for mirror asymmetry, with a Bayes factor of 12.6 or, equivalently, 93.1% probability. We obtain consistent (although stronger) results of 97.5% and 94.3%, respectively, using public results on this event from [M. Hannam et al., Nature (London) 610, 652 (2022).NATUAS0028-083610.1038/s41586-022-05212-z] and performing our own parameter inference. This finding further implies evidence of astrophysical sources that can spontaneously emit circularly polarized photons by quantum effects. Forthcoming black-hole merger detections will enable stronger constraints on large-scale mirror asymmetry and the cosmological principle.

Criteria for identifying and evaluating locations that could potentially host the Cosmic Explorer observatories

Cosmic Explorer is a next-generation ground-based gravitational-wave observatory that is being designed in the 2020s and is envisioned to begin operations in the 2030s together with the Einstein Telescope in Europe. The Cosmic Explorer concept currently consists of two widely separated L-shaped observatories in the United States, one with 40 km-long arms and the other with 20 km-long arms. This order of magnitude increase in scale with respect to the LIGO-Virgo-KAGRA observatories will, together with technological improvements, deliver an order of magnitude greater astronomical reach, allowing access to gravitational waves from remnants of the first stars and opening a wide discovery aperture to the novel and unknown. In addition to pushing the reach of gravitational-wave astronomy, Cosmic Explorer endeavors to approach the lifecycle of large scientific facilities in a way that prioritizes mutually beneficial relationships with local and Indigenous communities. This article describes the (scientific, cost and access, and social) criteria that will be used to identify and evaluate locations that could potentially host the Cosmic Explorer observatories.

Low-vibration cryogenic test facility for next generation of ground-based gravitational-wave observatories

We present the design and commissioning of a cryogenic low-vibration test facility that measures displacement noise from a gram-scale silicon cantilever at the level of 10-16m/Hz at 1 kHz. This sensitivity is necessary for future tests of thermal noise models on cross sections of silicon suspension samples proposed for future gravitational-wave detectors. A volume of ∼36 l is enclosed by radiation shields cooling an optical test cavity that is suspended from a multi-stage pendulum chain providing isolation from acoustic and environmental noise. This 3 kg test cavity housing a crystalline silicon cantilever is radiatively cooled to 123 K in 41 h and held at that temperature over many months with a relative temperature stability of ±1 mK. The facility housing the test cavity is sensitive to cavity length changes, which can resolve thermal fluctuations at the desired sensitivity. It is capable of interferometrically measuring temperature-dependent broadband displacement noise directly between 50 Hz and 10 kHz, where current and future ground-based gravitational wave observatories are the most sensitive. With a suitable cantilever design, the cryogenic facility we describe here will allow for the measurement of broadband thermal noise in crystalline silicon at 123 K. This will guide the design of suspensions in planned future cryogenic ground-based gravitational-wave detectors such as LIGO Voyager and may have implications for suspensions in the Einstein Telescope. This facility is also suitable for the testing of new mirror coatings at cryogenic temperatures.

Accurate Standard Siren Cosmology with Joint Gravitational-Wave and γ-Ray Burst Observations

Joint gravitational-wave and γ-ray burst (GRB) observations are among the best prospects for standard siren cosmology. However, the strong selection effect for the coincident GRB detection, which is possible only for sources with small inclination angles, induces a systematic uncertainty that is currently not accounted for. We show that this severe source of bias can be removed by inferring the a priori unknown electromagnetic detection probability directly from multimessenger data. This leads at the same time to an unbiased measurement of the Hubble constant, to constrain the properties of GRB emission, and to accurately measure the viewing angle of each source. Our inference scheme is applicable to real data already in the small-statistics regime, a scenario that might become reality in the near future. Additionally, we introduce a novel likelihood approximant for gravitational-wave events which treats the dependence on distance and inclination as exact.

Core Payload of the Space Gravitational Wave Observatory: Inertial Sensor and Its Critical Technologies

Since Einstein's prediction regarding the existence of gravitational waves was directly verified by the ground-based detector Advanced LIGO, research on gravitational wave detection has garnered increasing attention. To overcome limitations imposed by ground vibrations and interference at arm's length, a space-based gravitational wave detection initiative was proposed, which focuses on analyzing a large number of waves within the frequency range below 1 Hz. Due to the weak signal intensity, the TMs must move along their geodesic orbit with a residual acceleration less than 10-15 m/s2/Hz1/2. Consequently, the core payload-inertial sensor was designed to shield against stray force noise while maintaining the high-precision motion of the test mass through a drag-free control system, providing an ultra-stable inertial reference for laser interferometry. To meet these requirements, the inertial sensor integrates a series of unit settings and innovative designs, involving numerous subsystems and technologies. This article provides a comprehensive overview of these critical technologies used in the development of inertial sensors for space gravitational wave detection and discusses future trends and potential applications for these sensors.

Post-Minkowskian Theory Meets the Spinning Effective-One-Body Approach for Bound-Orbit Waveforms

Driven by advances in scattering amplitudes and worldline-based methods, recent years have seen significant progress in our ability to calculate gravitational two-body scattering observables. These observables effectively encapsulate the gravitational two-body problem in the weak-field and high-velocity regime [post-Minkowskian (PM)], with applications to the bound two-body problem and gravitational-wave modeling. We leverage PM data to construct a complete inspiral-merger-ringdown waveform model for nonprecessing spinning black holes within the effective-one-body (EOB) formalism SEOBNR-PM. This model is closely based on the highly successful SEOBNRv5 model, used by the LIGO-Virgo-KAGRA Collaboration, with its key new feature being an EOB Hamiltonian derived by matching the two-body scattering angle in a perturbative PM expansion. The model performs remarkably well, showing a median mismatch against 441 numerical-relativity (NR) simulations that is somewhat lower than a similarly calibrated version of SEOBNRv5. Comparisons of the binding energy with NR also demonstrate better agreement than SEOBNRv5, despite the latter containing additional calibration to NR simulations.

Enhanced Extreme Mass Ratio Inspiral Rates and Intermediate Mass Black Holes

Extreme mass ratio inspirals (EMRIs) occur when stellar-mass compact objects begin a gravitational wave (GW) driven inspiral into massive black holes. EMRI waveforms can precisely map the surrounding spacetime, making them a key target for future space-based GW interferometers such as LISA, but their event rates and parameters are massively uncertain. One of the largest uncertainties is the ratio of true EMRIs (which spend at least thousands of orbits in the LISA band) and direct plunges, which are in-band for at most a handful of orbits and are not detectable in practice. In this Letter, we show that the traditional dichotomy between EMRIs and plunges-EMRIs originate from small semimajor axes, plunges from large-does not hold for intermediate-mass black holes with masses M_{•}≲10^{5}M_{⊙}. In this low-mass regime, a plunge always has an O(1) probability of failing and transitioning into a novel "cliffhanger" EMRI. Cliffhanger EMRIs are more easily produced for larger stellar-mass compact objects, and are less likely for smaller ones. This new EMRI production channel can dominate volumetric EMRI rates n[over ˙]_{EMRI} if intermediate-mass black holes are common in dwarf galactic nuclei, potentially increasing n[over ˙]_{EMRI} by an order of magnitude.

On networks of space-based gravitational-wave detectors

The space-based laser interferometers, LISA, Taiji and TianQin, are targeting to observe milliHz gravitational waves (GWs) in the 2030s. The joint observations from multiple space-based detectors yield significant advantages. In this work, we recap the studies and investigations for the joint space-based GW detector networks to highlight: 1) the high precision of sky localization for the massive binary black hole (BBH) coalescences and the GW sirens in the cosmological implication, 2) the effectiveness to test the parity violation in the stochastic GW background observations, 3) the efficiency of subtracting galactic foreground, 4) the improvement in stellar-mass BBH observations. We inspect alternative networks by trading off massive BBH observations and stochastic GW background observation.

Detecting single gravitons with quantum sensing

The quantization of gravity is widely believed to result in gravitons - particles of discrete energy that form gravitational waves. But their detection has so far been considered impossible. Here we show that signatures of single graviton exchange can be observed in laboratory experiments. We show that stimulated and spontaneous single-graviton processes can become relevant for massive quantum acoustic resonators and that stimulated absorption can be resolved through continuous sensing of quantum jumps. We analyze the feasibility of observing the exchange of single energy quanta between matter and gravitational waves. Our results show that single graviton signatures are within reach of experiments. In analogy to the discovery of the photo-electric effect for photons, such signatures can provide the first experimental clue of the quantization of gravity.

Estimate for the Bulk Viscosity of Strongly Coupled Quark Matter Using Perturbative QCD and Holography

Modern hydrodynamic simulations of core-collapse supernovae and neutron-star mergers require knowledge not only of the equilibrium properties of strongly interacting matter, but also of the system's response to perturbations, encoded in various transport coefficients. Using perturbative and holographic tools, we derive here an improved weak-coupling and a new strong-coupling result for the most important transport coefficient of unpaired quark matter, its bulk viscosity. These results are combined in a simple analytic pocket formula for the quantity that is rooted in perturbative quantum chromodynamics at high densities but takes into account nonperturbative holographic input at neutron-star densities, where the system is strongly coupled. This expression can be used in the modeling of unpaired quark matter at astrophysically relevant temperatures and densities.

Imprints of massive black-hole binaries on neighbouring decihertz gravitational-wave sources

The most massive black holes in our Universe form binaries at the centre of merging galaxies. The recent evidence for a gravitational-wave (GW) background from pulsar timing may constitute the first observation that these supermassive black-hole binaries (SMBHBs) merge. Yet, the most massive SMBHBs are out of reach of interferometric GW detectors and are exceedingly difficult to resolve individually with pulsar timing. These limitations call for unexplored strategies to detect individual SMBHBs in the uncharted frequency band ≲10-5 Hz to establish their abundance and decipher the coevolution with their host galaxies. Here we show that SMBHBs imprint detectable long-term modulations on GWs from stellar-mass binaries residing in the same galaxy at a distance d ≲ 1 kpc. We determine that proposed decihertz GW interferometers sensitive to numerous stellar-mass binaries could uncover modulations from ~O(10-1-104) SMBHBs with masses ~O(107-108) M⊙out to redshift z ≈ 3.5. This offers a unique opportunity to map the population of SMBHBs through cosmic time, which might remain inaccessible otherwise.

Methodological notes on gauge invariance in the treatment of waves and oscillations in plasmas via the Einstein-Vlasov-Maxwell system: Fundamental equations

The theory of gauge transformations in linearized gravitation is investigated. After a brief discussion of the fundamentals of the kinetic theory in curved spacetime, the Einstein-Vlasov-Maxwell (EVM) system of equations in terms of gauge-invariant quantities is established without neglecting the equations of motion associated with the dynamics of the nonradiative components of the metric tensor. The established theory is applied to a noncollisional electron-positron plasma, leading to a dispersion relation for gravitational waves in this model system. The problem of Landau damping is addressed and some attention is given to the issue of the energy exchanges between the plasma and the gravitational wave. In a future paper, a more complete set of approximate dispersion relations for waves and oscillations in plasmas will be presented, including the dynamics of nonradiative components of the metric tensor, with special attention to the problems of the Landau damping and of the energy exchanges between matter, the electromagnetic field and the gravitational field.

A Field-Theory Approach for Modeling Dissipative Relativistic Fluids

We develop an action principle for producing a single-fluid two-constituent system with dissipation in general relativity. The two constituents in the model are particles and entropy. The particle flux creation rate is taken to be zero, while the entropy creation rate is non-zero. Building on previous work, it is demonstrated that a new term (the proper time derivative of the matter space "metric") is required in the Lagrangian in order to produce terms typically associated with bulk and shear viscosity. Equations of motion, entropy creation rate, and energy-momentum-stress tensor are derived. Using an Onsager approach of identifying thermodynamic "forces" and "fluxes", a model is produced which delivers the same entropy creation rate as the standard, relativistic Navier-Stokes equations. This result is then contrasted with a model generated in the spirit of the action principle, which takes as its starting point a specific Lagrangian and then produces the equations of motion, entropy creation rate, and energy-momentum-stress tensor. Unlike the equations derived from Onsager reasoning, where the analogs of the bulk and shear viscosity coefficients are prescribed "externally", we find that the forms of the coefficients in the second example are a direct result of the specified Lagrangian. Furthermore, the coefficients are shown to satisfy evolution equations along the fluid worldline, also a product of the specific Lagrangian.

Photon strength functions and nuclear level densities: invaluable input for nucleosynthesis

The pivotal role of nuclear physics in nucleosynthesis processes is being investigated, in particular the intricate influence of photon strength functions (PSFs) and nuclear level densities (NLDs) on shaping the outcomes of the i-, r- and p-processes. Exploring diverse NLD and PSF model combinations uncovers large uncertainties for (p,[Formula: see text]), (n,[Formula: see text]) and ([Formula: see text],[Formula: see text]) rates across many regions of the nuclear chart. These lead to potentially significant abundance variations of the nucleosynthesis processes and highlight the importance of accurate experimental nuclear data. Theoretical insights and advanced experimental techniques lay the ground work for profound understanding that can be gained of nucleosynthesis mechanisms and the origin of the elements. Recent results further underscore the effect of PSF and NLD data and its contribution to understanding abundance distributions and refining knowledge of the intricate nucleosynthesis processes. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

Nuclear structure opportunities with GeV radioactive beams at FAIR

The Facility for Antiproton and Ion Research (FAIR) is in its final construction stage next to the campus of the Gesellschaft für Schwerionenforschung Helmholtzzentrum for heavy-ion research in Darmstadt, Germany. Once it starts its operation, it will be the main nuclear physics research facility in many basic sciences and their applications in Europe for the coming decades. Owing to the ability of the new fragment separator, Super-FRagment Separator, to produce high-intensity radioactive ion beams in the energy range up to about 2 GeV/nucleon, these can be used in various nuclear reactions. This opens a unique opportunity for various nuclear structure studies across a range of fields and scales: from low-energy physics via the investigation of multi-neutron systems and halos to high-density nuclear matter and the equation of state, following heavy-ion collisions, fission and study of short-range correlations in nuclei and hypernuclei. The newly developed reactions with relativistic radioactive beams (R3B) set up at FAIR would be the most suitable and versatile for such studies. An overview of highlighted physics cases foreseen at R3B is given, along with possible future opportunities, at FAIR. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

A decomposition of light's spin angular momentum density

Light carries intrinsic spin angular momentum (SAM) when the electric or magnetic field vector rotates over time. A familiar vector equation calculates the direction of light's SAM density using the right-hand rule with reference to the electric and magnetic polarisation ellipses. Using Maxwell's equations, this vector equation can be decomposed into a sum of two distinct terms, akin to the well-known Poynting vector decomposition into orbital and spin currents. We present the first general study of this spin decomposition, showing that the two terms, which we call canonical and Poynting spin, are chiral analogies to the canonical and spin momenta of light in its interaction with matter. Like canonical momentum, canonical spin is directly measurable. Both canonical and Poynting spin incorporate spatial variation of the electric and magnetic fields and are influenced by optical vortices. The decomposition allows us to show that a linearly polarised vortex beam, which has no total SAM, can nevertheless exert longitudinal chiral pressure due to equal and opposite canonical and Poynting spins.

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J0348+0432, PSR J1624-2230, PSR J0740+6620, GW170817). At baryon number chemical potential μ=2.6  GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. This constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

Conservative Binary Dynamics at Order α^{5} in Electrodynamics

We compute the potential-photon contributions to the classical relativistic scattering angle of two charged nonspinning bodies in electrodynamics through fifth order in the coupling. We use the scattering amplitudes framework, effective field theory, and multiloop integration techniques based on integration by parts and differential equations. At fifth order, the result is expressed in terms of cyclotomic polylogarithms. Our calculation demonstrates the feasibility of the corresponding calculations in general relativity, including the evaluation of the encountered four-loop integrals.

Conservative Black Hole Scattering at Fifth Post-Minkowskian and First Self-Force Order

We compute the fifth post-Minkowskian (5PM) order contributions to the scattering angle and impulse of classical black hole scattering in the conservative sector at first self-force order using the worldline quantum field theory formalism. This challenging four-loop computation required the use of advanced integration-by-parts and differential equation technology implemented on high-performance computing systems. Use of partial fraction identities allowed us to render the complete integrand in a fully planar form. The resulting function space is simpler than expected: In the scattering angle, we see only multiple polylogarithms up to weight three and a total absence of the elliptic integrals that appeared at 4PM order. All checks on our result, both internal-cancellation of dimensional regularization poles and preservation of the on-shell condition-and external-matching the slow-velocity limit with the post-Newtonian (PN) literature up to 5PN order and matching the tail terms to the 4PM loss of energy-are passed.

Abundances of Neutron-capture Elements in 62 Stars in the Globular Cluster Messier 15

M15 is a globular cluster with a known spread in neutron-capture elements. This paper presents abundances of neutron-capture elements for 62 stars in M15. Spectra were obtained with the Michigan/Magellan Fiber System spectrograph, covering a wavelength range from ∼4430 to 4630 Å. Spectral lines from Fe i, Fe ii, Sr i, Zr ii, Ba ii, La ii, Ce ii, Nd ii, Sm ii, Eu ii, and Dy ii were measured, enabling classifications and neutron-capture abundance patterns for the stars. Of the 62 targets, 44 are found to be highly Eu-enhanced r-II stars, another 17 are moderately Eu-enhanced r-I stars, and one star is found to have an s-process signature. The neutron-capture patterns indicate that the majority of the stars are consistent with enrichment by the r-process. The 62 target stars are found to show significant star-to-star spreads in Sr, Zr, Ba, La, Ce, Nd, Sm, Eu, and Dy, but no significant spread in Fe. The neutron-capture abundances are further found to have slight correlations with sodium abundances from the literature, unlike what has been previously found; follow-up studies are needed to verify this result. The findings in this paper suggest that the Eu-enhanced stars in M15 were enhanced by the same process, that the nucleosynthetic source of this Eu pollution was the r-process, and that the r-process source occurred as the first generation of cluster stars was forming.

Calabi-Yau Meets Gravity: A Calabi-Yau Threefold at Fifth Post-Minkowskian Order

We study geometries occurring in Feynman integrals that contribute to the scattering of black holes in the post-Minkowskian (PM) expansion. These geometries become relevant to gravitational-wave production from binary mergers through the classical conservative potential. At 4PM, a K3 surface is known to occur in a three-loop integral, leading to elliptic integrals in the result. In this Letter, we identify a Calabi-Yau threefold in a four-loop integral, contributing at 5PM. The presence of this Calabi-Yau geometry indicates that completely new functions occur in the full analytical results.

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.

QCD challenges from pp to AA collisions: 4th edition

This paper is a write-up of the ideas that were presented, developed and discussed at the fourth International Workshop on QCD Challenges from pp to AA, which took place in February 2023 in Padua, Italy. The goal of the workshop was to focus on some of the open questions in the field of high-energy heavy-ion physics and to stimulate the formulation of concrete suggestions for making progresses on both the experimental and theoretical sides. The paper gives a brief introduction to each topic and then summarizes the primary results.

Nuclear Charge Radii of Silicon Isotopes

The nuclear charge radius of ^{32}Si was determined using collinear laser spectroscopy. The experimental result was confronted with ab initio nuclear lattice effective field theory, valence-space in-medium similarity renormalization group, and mean field calculations, highlighting important achievements and challenges of modern many-body methods. The charge radius of ^{32}Si completes the radii of the mirror pair ^{32}Ar-^{32}Si, whose difference was correlated to the slope L of the symmetry energy in the nuclear equation of state. Our result suggests L≤60  MeV, which agrees with complementary observables.

Kerr-Enhanced Optical Spring

We propose and experimentally demonstrate the generation of enhanced optical springs using the optical Kerr effect. A nonlinear optical crystal is inserted into a Fabry-Perot cavity with a movable mirror, and a chain of second-order nonlinear optical effects in the phase-mismatched condition induces the Kerr effect. The optical spring constant is enhanced by a factor of 1.6±0.1 over linear theory. To our knowledge, this is the first realization of optomechanical coupling enhancement using a nonlinear optical effect, which has been theoretically investigated to overcome the performance limitations of linear optomechanical systems. The tunable nonlinearity of demonstrated system has a wide range of potential applications, from observing gravitational waves emitted by binary neutron star postmerger remnants to cooling macroscopic oscillators to their quantum ground state.

Pulsar glitches from quantum vortex networks

Neutron stars or pulsars are very rapidly rotating compact stars with extremely high density. One of the unsolved long-standing problems of these enigmatic celestial bodies is the origin of pulsars' glitches, i.e., the sudden rapid deceleration in the rotation speed of neutron stars. Although many glitch events have been reported, there is no consensus on the microscopic mechanism responsible for them. One of the important characterizations of the glitches is the scaling law P ( E ) ∼ E - α of the probability distribution for a glitch with energy E. Here, we reanalyse the accumulated up-to-date observation data to obtain the exponent α ≈ 0.88 for the scaling law, and propose a simple microscopic model that naturally deduces this scaling law without any free parameters. Our model explains the appearance of these glitches in terms of the presence of quantum vortex networks arising at the interface of two different kinds of superfluids in the core of neutron stars; a p-wave neutron superfluid in the inner core which interfaces with the s-wave neutron superfluid in the outer core, where each integer vortex in the s-wave superfluid connects to two half-quantized vortices in the p-wave superfluid through structures called "boojums".

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.

Absolute Ranging with Time Delay Interferometry for Space-Borne Gravitational Wave Detection

In future space-borne gravitational wave (GW) detectors, time delay interferometry (TDI) will be utilized to reduce the overwhelming noise, including the laser frequency noise and the clock noise etc., by time shifting and recombining the data streams in post-processing. The successful operation of TDI relies on absolute inter-satellite ranging with meter-level precision. In this work, we numerically and experimentally demonstrate a strategy for inter-satellite distance measurement. The distances can be coarsely determined using the technique of arm-locking ranging with a large non-ambiguity range, and subsequently TDI can be used for precise distance measurement (TDI ranging) by finding the minimum value of the power of the residual noises. The measurement principle is introduced. We carry out the numerical simulations, and the results show millimeter-level precision. Further, we perform the experimental verifications based on the fiber link, and the distances can be measured with better than 0.05 m uncertainty, which can well satisfy the requirement of time delay interferometry.

Momentum Shell in Quarkyonic Matter from Explicit Duality: A Dual Model for Cold, Dense QCD

We present a model of cold QCD matter that bridges nuclear and quark matter through the duality relation between quarks and baryons. The baryon number and energy densities are expressed as functionals of either the baryon momentum distribution, f_{B}, or the quark distribution, f_{Q}, which are subject to the constraints on fermions, 0≤f_{B,Q}≤1. The theory is ideal in the sense that the confinement of quarks into baryons is reflected in the duality relation between f_{Q} and f_{B}, while other possible interactions among quarks and baryons are all neglected. The variational problem with the duality constraints is formulated and we explicitly construct analytic solutions, finding two distinct regimes: a nuclear matter regime at low density and a quarkyonic regime at high density. In the quarkyonic regime, baryons underoccupy states at low momenta but form a momentum shell with f_{B}=1 on top of a quark Fermi sea. Such a theory describes a rapid transition from a soft nuclear equation of state to a stiff quarkyonic equation of state. At this transition, there is a rapid increase in the pressure.

Multimessenger Constraints on Radiatively Decaying Axions from GW170817

The metastable hypermassive neutron star produced in the coalescence of two neutron stars can copiously produce axions that radiatively decay into O(100)  MeV photons. These photons can form a fireball with characteristic temperature smaller than 1 MeV. By relying on x-ray observations of GW170817/GRB 170817A with CALET CGBM, Konus-Wind, and Insight-HXMT/HE, we present new bounds on the axion-photon coupling for axion masses in the range 1-400 MeV. We exclude couplings down to 5×10^{-11}  GeV^{-1}, complementing and surpassing existing constraints. Our approach can be extended to any feebly interacting particle decaying into photons.

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.

Binary Coalescences as Sources of Ultrahigh-Energy Cosmic Rays

Binary coalescences are known sources of gravitational waves (GWs) and they encompass combinations of black holes (BHs) and neutron stars (NSs). Here we show that when BHs are embedded in magnetic fields (B's) larger than approximately 10^{10}  G, charged particles colliding around their event horizons can easily have center-of-mass energies in the range of ultrahigh energies (≳10^{18}  eV) and become more likely to escape. Such B-embedding and high-energy particles can take place in BH-NS binaries, or even in BH-BH binaries with one of the BHs being charged (with charge-to-mass ratios as small as 10^{-5}, which do not change GW waveforms) and having a residual accretion disk. Ultrahigh center-of-mass energies for particle collisions arise for basically any rotation parameter of the BH when B≳10^{10}  G, meaning that it should be a common aspect in binaries, especially in BH-NS ones given the natural presence of a B onto the BH and charged particles due to the magnetosphere of the NS. We estimate that the number of ultrahigh center-of-mass collisions ranges from a few up to millions before the merger of binary compact systems. Thus, binary coalescences may also be efficient sources of ultrahigh energy cosmic rays (UHECRs) and constraints to NS/BH parameters would be possible if UHECRs are detected along with GWs.

Numerical investigation of effective nonlinear coefficient model for coupled third harmonic generation

In this paper, the optimal solution of effective nonlinear coefficient of quasi-phase-matching (QPM) crystals for coupled third harmonic generation (CTHG) was numerically investigated. The effective nonlinear coefficient of CTHG was converted to an Ising model for optimizing domain distributions of aperiodically poled lithium niobate (APPLN) crystals with lengths as 0.5 mm and 1 mm, and fundamental wavelengths ranging from 1000 nm to 6000 nm. A method for reconstructing crystal domain poling weight curve of coupled nonlinear processes was also proposed, which demonstrated the optimal conversion ratio between two coupled nonlinear processes at each place along the crystal. In addition, by applying the semidefinite programming, the upper bound on the effective nonlinear coefficients deff for different fundamental wavelengths were calculated. The research can be extended to any coupled dual χ(2) process and will help us to understand better the dynamics of coupled nonlinear interactions based on QPM crystals.