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'.

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.

JWST detection of a supernova associated with GRB 221009A without an r-process signature

Identifying the sites of r-process nucleosynthesis, a primary mechanism of heavy element production, is a key goal of astrophysics. The discovery of the brightest gamma-ray burst (GRB) to date, GRB 221009A, presented an opportunity to spectroscopically test the idea that r-process elements are produced following the collapse of rapidly rotating massive stars. Here we present James Webb Space Telescope observations of GRB 221009A obtained +168 and +170 rest-frame days after the gamma-ray trigger, and demonstrate that they are well described by a SN 1998bw-like supernova (SN) and power-law afterglow, with no evidence for a component from r-process emission. The SN, with a nickel mass of approximately 0.09 M ⊙, is only slightly fainter than the brightness of SN 1998bw at this phase, which indicates that the SN is not an unusual GRB-SN. This demonstrates that the GRB and SN mechanisms are decoupled and that highly energetic GRBs are not likely to produce significant quantities of r-process material, which leaves open the question of whether explosions of massive stars are key sources of r-process elements. Moreover, the host galaxy of GRB 221009A has a very low metallicity of approximately 0.12 Z ⊙ and strong H2 emission at the explosion site, which is consistent with recent star formation, hinting that environmental factors are responsible for its extreme energetics.

Thallium-208: A Beacon of In Situ Neutron Capture Nucleosynthesis

We demonstrate that the well-known 2.6 MeV gamma-ray emission line from thallium-208 could serve as a real-time indicator of astrophysical heavy element production, with both rapid (r) and intermediate (i) neutron capture processes capable of its synthesis. We consider the r process in a Galactic neutron star merger and show Tl-208 to be detectable from ∼12 hours to ∼ten days, and again ∼1-20  years postevent. Detection of Tl-208 represents the only identified prospect for a direct signal of lead production (implying gold synthesis), arguing for the importance of future MeV telescope missions which aim to detect Galactic events but may also be able to reach some nearby galaxies in the Local Group.

Minutes-duration optical flares with supernova luminosities

In recent years, certain luminous extragalactic optical transients have been observed to last only a few days1. Their short observed duration implies a different powering mechanism from the most common luminous extragalactic transients (supernovae), whose timescale is weeks2. Some short-duration transients, most notably AT2018cow (ref. 3), show blue optical colours and bright radio and X-ray emission4. Several AT2018cow-like transients have shown hints of a long-lived embedded energy source5, such as X-ray variability6,7, prolonged ultraviolet emission8, a tentative X-ray quasiperiodic oscillation9,10 and large energies coupled to fast (but subrelativistic) radio-emitting ejecta11,12. Here we report observations of minutes-duration optical flares in the aftermath of an AT2018cow-like transient, AT2022tsd (the 'Tasmanian Devil'). The flares occur over a period of months, are highly energetic and are probably nonthermal, implying that they arise from a near-relativistic outflow or jet. Our observations confirm that, in some AT2018cow-like transients, the embedded energy source is a compact object, either a magnetar or an accreting black hole.

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.

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.

Thinking Outside the Box: Numerical Relativity with Particles

The observation of gravitational waves from compact objects has now become an active part of observational astronomy. For a sound interpretation, one needs to compare such observations against detailed Numerical Relativity simulations, which are essential tools to explore the dynamics and physics of compact binary mergers. To date, essentially all simulation codes that solve the full set of Einstein’s equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations, we focus on perfectly symmetric binary systems that are irrotational and have 1.3 M⊙ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle–mesh mapping), and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the Ψ4 formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times 10−3M⊙ as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity (∼0.5…0.7c) ejecta component of ∼10−4M⊙ that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission, and it could also potentially cause a kilonova afterglow.

Shell Model Applications in Nuclear Astrophysics

In recent years, shell model studies have significantly contributed in improving the nuclear input, required in simulations of the dynamics of astrophysical objects and their associated nucleosynthesis. This review highlights a few examples such as electron capture rates and neutrino-nucleus cross sections, important for the evolution and nucleosynthesis of supernovae. For simulations of rapid neutron-capture (r-process) nucleosynthesis, shell model studies have contributed to an improved understanding of half lives of neutron-rich nuclei with magic neutron numbers and of the nuclear level densities and γ-strength functions that are both relevant for neutron capture rates.

Dark Matter Effects on the Compact Star Properties

The neutron star properties are generally determined by the equation of state of β-equilibrated dense matter. In this work, we consider the interaction of fermionic dark matter (DM) particles with nucleons via Higgs exchange and investigate the effect on the neutron star properties with the relativistic mean-field model equation of state coupled with DM. We deduce that DM significantly affects the neutron star properties, such as considerably reducing the maximum mass of the star, which depends on the percentage of the DM considered inside the neutron star. The tidal Love numbers both for electric and magnetic cases and surficial Love numbers are also studied for DM admixed NS. We observed that the magnitude of tidal and surficial Love numbers increases with a greater DM percentage. Further, we present post-Newtonian tidal corrections to gravitational waves decreased by increasing the DM percentage. The DM effect on the GW signal is significant during the late inspiral and merger stages of binary evolution for GW frequencies >500 Hz.

Probing Gamma-Ray Burst VHE Emission with the Southern Wide-Field-of-View Gamma-Ray Observatory

Recent observations have confirmed that Gamma-Ray Burst (GRB) afterglows produce Very High-Energy radiation (VHE, E>100GeV). This highly anticipated discovery opens new scenarios in the interpretation of GRBs and in their role as probes of Extragalactic Background Light (EBL) and Lorentz Invariance Violation (LIV). However, some fundamental questions about the actual nature of VHE emission in GRBs and its evolution during the burst are still unsolved. These questions will be difficult to address, even with future imaging Cherenkov telescopes, such as the Cherenkov Telescope Array (CTA). Here we investigate the prospects of gamma-ray sky monitoring with Extensive Air Showers arrays (EAS) to address these problems. We discuss the theoretical aspects connected with VHE radiation emission and the implications that its temporal evolution properties have on the interpretation of GRBs. By revisiting the high-energy properties of some Fermi-LAT detected GRBs, we estimate the typical fluxes expected in the VHE band and compare them with a range of foreseeable instrument performances, based on the Southern Wide Field-of-view Gamma-ray Observatory concept (SWGO). We focus our analysis on how different instrument capabilities affect the chances to explore the burst onset and early evolution in VHE, providing invaluable complementary information with respect to Cherenkov telescope observations. We show that under the assumption of conditions already observed in historical events, the next-generation ground monitoring detectors can actually contribute to answer several key questions.

A Neutron Star Is Born

A neutron star was first detected as a pulsar in 1967. It is one of the most mysterious compact objects in the universe, with a radius of the order of 10 km and masses that can reach two solar masses. In fact, neutron stars are star remnants, a kind of stellar zombie (they die, but do not disappear). In the last decades, astronomical observations yielded various contraints for neutron star masses, and finally, in 2017, a gravitational wave was detected (GW170817). Its source was identified as the merger of two neutron stars coming from NGC 4993, a galaxy 140 million light years away from us. The very same event was detected in γ-ray, X-ray, UV, IR, radio frequency and even in the optical region of the electromagnetic spectrum, starting the new era of multi-messenger astronomy. To understand and describe neutron stars, an appropriate equation of state that satisfies bulk nuclear matter properties is necessary. GW170817 detection contributed with extra constraints to determine it. On the other hand, magnetars are the same sort of compact object, but bearing much stronger magnetic fields that can reach up to 1015 G on the surface as compared with the usual 1012 G present in ordinary pulsars. While the description of ordinary pulsars is not completely established, describing magnetars poses extra challenges. In this paper, I give an overview on the history of neutron stars and on the development of nuclear models and show how the description of the tiny world of the nuclear physics can help the understanding of the cosmos, especially of the neutron stars.

Nuclear Physics and Astrophysics Constraints on the High Density Matter Equation of State

(1) This review has been written in memory of Steven Moszkowski who unexpectedly passed away in December 2020. It has been inspired by our many years of discussions. Steven’s enthusiasm, drive and determination to understand atomic nuclei in simple terms of basic laws of physics was infectious. He sought the fundamental origin of nuclear forces in free space, and their saturation and modification in nuclear medium. His untimely departure left our job unfinished but his legacy lives on. (2) Focusing on the nuclear force acting in nuclear matter of astrophysical interest and its equation of state (EoS), we take several typical snapshots of evolution of the theory of nuclear forces. We start from original ideas in the 1930s moving through to its overwhelming diversity today. The development is supported by modern observational and terrestrial data and their inference in the multimessenger era, as well as by novel mathematical techniques and computer power. (3) We find that, despite the admirable effort both in theory and measurement, we are facing multiple models dependent on a large number of variable correlated parameters which cannot be constrained by data, which are not yet accurate, nor sensitive enough, to identify the theory closest to reality. The role of microphysics in the theories is severely limited or neglected, mostly deemed to be too difficult to tackle. (4) Taking the EoS of high-density matter as an example, we propose to develop models, based, as much as currently possible, on the microphysics of the nuclear force, with a minimal set of parameters, chosen under clear physical guidance. Still somewhat phenomenological, such models could pave the way to realistic predictions, not tracing the measurement, but leading it.

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.

129I and 247Cm in meteorites constrain the last astrophysical source of solar r-process elements

The composition of the early Solar System can be inferred from meteorites. Many elements heavier than iron were formed by the rapid neutron capture process (r-process), but the astrophysical sources where this occurred remain poorly understood. We demonstrate that the near-identical half-lives [Formula: see text] of the radioactive r-process nuclei iodine-129 and curium-247 preserve their ratio, irrespective of the time between production and incorporation into the Solar System. We constrain the last r-process source by comparing the measured meteoritic ratio ¹²⁹I/²⁴⁷Cm = 438 ± 184 with nucleosynthesis calculations based on neutron star merger and magneto-rotational supernova simulations. Moderately neutron-rich conditions, often found in merger disk ejecta simulations, are most consistent with the meteoritic value. Uncertain nuclear physics data limit our confidence in this conclusion.

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.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Neutron star merger remnants

Binary neutron star mergers observations are a unique way to constrain fundamental physics and astrophysics at the extreme. The interpretation of gravitational-wave events and their electromagnetic counterparts crucially relies on general-relativistic models of the merger remnants. Quantitative models can be obtained only by means of numerical relativity simulations in 3 + 1 dimensions including detailed input physics for the nuclear matter, electromagnetic and weak interactions. This review summarizes the current understanding of merger remnants focusing on some of the aspects that are relevant for multimessenger observations.

The origin of the elements: a century of progress

This review assesses the current state of knowledge of how the elements were produced in the Big Bang, in stellar lives and deaths, and by interactions in interstellar gas. We begin with statements of fact and discuss the evidence that convinced astronomers that the Sun is fusing hydrogen, that low-mass stars produce heavy elements through neutron capture, that massive stars can explode as supernovae and that supernovae of all types produce new elements. Nucleosynthesis in the Big Bang, through cosmic ray spallation, and in exploding white dwarfs is only ranked below the above facts in certainty because the evidence, while overwhelming, is so far circumstantial. Next, we highlight the flaws in our current understanding of the predictions for lithium production in the Big Bang and/or its destruction in stars and for the production of the elements with atomic number [Formula: see text]. While the theory that neutron star mergers produce elements through neutron-capture has powerful circumstantial evidence, we are unconvinced that they produce all of the elements past nickel. Also in dispute is the exact mechanism or mechanisms that cause the white dwarfs to explode. It is difficult to determine the origin of rare isotopes because signatures of their production are weak. We are uncertain about the production sites of some lithium and nitrogen isotopes and proton-rich heavy nuclei. Finally, Betelgeuse is probably not the next star to become a supernovae in the Milky Way, in part because Betelgeuse may collapse directly to a black hole instead. The accumulated evidence in this review shows that we understand the major production sites for the elements, but islands of uncertainty in the periodic table exist. Resolving these uncertainties requires in particular understanding explosive events with compact objects and understanding the nature of the first stars and is therefore primed for new discoveries in the next decades. This article is part of the theme issue 'Mendeleev and the periodic table'.

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron-star merger GW170817 remains the only confirmed electromagnetic counterpart to a gravitational-wave event^1,2. Observations of the late-time electromagnetic emission, however, do not agree with the expectations from standard neutron-star merger models. Although the large measured ejecta mass^3,4 could be explained by a progenitor system that is asymmetric in terms of the stellar component masses (that is, with a mass ratio q of 0.7 to 0.8)⁵, the known Galactic population of merging double neutron-star systems (that is, those that will coalesce within billions of years or less) has until now consisted only of nearly equal-mass (q > 0.9) binaries⁶. The pulsar PSR J1913+1102 is a double system in a five-hour, low-eccentricity (0.09) orbit, with an orbital separation of 1.8 solar radii⁷, and the two neutron stars are predicted to coalesce in [Formula: see text] million years owing to gravitational-wave emission. Here we report that the masses of the pulsar and the companion neutron star, as measured by a dedicated pulsar timing campaign, are 1.62 ± 0.03 and 1.27 ± 0.03 solar masses, respectively. With a measured mass ratio of q = 0.78 ± 0.03, this is the most asymmetric merging system reported so far. On the basis of this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30 per cent (90 per cent confidence) of the total population of merging binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

The periodic table and the physics that drives it

Mendeleev's introduction of the periodic table of elements is one of the most important milestones in the history of chemistry, as it brought order into the known chemical and physical behaviour of the elements. The periodic table can be seen as parallel to the Standard Model in particle physics, in which the elementary particles known today can be ordered according to their intrinsic properties. The underlying fundamental theory to describe the interactions between particles comes from quantum theory or, more specifically, from quantum field theory and its inherent symmetries. In the periodic table, the elements are placed into a certain period and group based on electronic configurations that originate from the Pauli and Aufbau principles for the electrons surrounding a positively charged nucleus. This order enables us to approximately predict the chemical and physical properties of elements. Apparent anomalies can arise from relativistic effects, partial-screening phenomena (of type lanthanide contraction) and the compact size of the first shell of every l-value. Further, ambiguities in electron configurations and the breakdown of assigning a dominant configuration, owing to configuration mixing and dense spectra for the heaviest elements in the periodic table. For the short-lived transactinides, the nuclear stability becomes an important factor in chemical studies. Nuclear stability, decay rates, spectra and reaction cross sections are also important for predicting the astrophysical origin of the elements, including the production of the heavy elements beyond iron in supernova explosions or neutron-star mergers. In this Perspective, we critically analyse the periodic table of elements and the current status of theoretical predictions and origins for the heaviest elements, which combine both quantum chemistry and physics.

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.

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.

A nearby neutron-star merger explains the actinide abundances in the early Solar System

A growing body of evidence indicates that binary neutron-star mergers are the primary origin of heavy elements produced exclusively through rapid neutron capture^1-4 (the 'r-process'). As neutron-star mergers occur infrequently, their deposition of radioactive isotopes into the pre-solar nebula could have been dominated by a few nearby events. Although short-lived r-process isotopes-with half-lives shorter than 100 million years-are no longer present in the Solar System, their abundances in the early Solar System are known because their daughter products were preserved in high-temperature condensates found in meteorites⁵. Here we report that abundances of short-lived r-process isotopes in the early Solar System point to their origin in neutron-star mergers, and indicate substantial deposition by a single nearby merger event. By comparing numerical simulations with the early Solar System abundance ratios of actinides produced exclusively through the r-process, we constrain the rate of occurrence of their Galactic production sites to within about 1-100 per million years. This is consistent with observational estimates of neutron-star merger rates^6-8, but rules out supernovae and stellar sources. We further find that there was probably a single nearby merger that produced much of the curium and a substantial fraction of the plutonium present in the early Solar System. Such an event may have occurred about 300 parsecs away from the pre-solar nebula, approximately 80 million years before the formation of the Solar System.

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.

Numerical relativity of compact binaries in the 21st century

We review the dramatic progress in the simulations of compact objects and compact-object binaries that has taken place in the first two decades of the twenty-first century. This includes simulations of the inspirals and violent mergers of binaries containing black holes and neutron stars, as well as simulations of black-hole formation through failed supernovae and high-mass neutron star-neutron star mergers. Modeling such events requires numerical integration of the field equations of general relativity in three spatial dimensions, coupled, in the case of neutron-star containing binaries, with increasingly sophisticated treatment of fluids, electromagnetic fields, and neutrino radiation. However, it was not until 2005 that accurate long-term evolutions of binaries containing black holes were even possible (Pretorius 2005 Phys. Rev. Lett. 95 121101, Campanelli et al 2006 Phys. Rev. Lett. 96 111101, Baker et al 2006 Phys. Rev. Lett. 96 111102). Since then, there has been an explosion of new results and insights into the physics of strongly-gravitating system. Particular emphasis has been placed on understanding the gravitational wave and electromagnetic signatures from these extreme events. And with the recent dramatic discoveries of gravitational waves from merging black holes by the Laser Interferometric Gravitational Wave Observatory and Virgo, and the subsequent discovery of both electromagnetic and gravitational wave signals from a merging neutron star-neutron star binary, numerical relativity became an indispensable tool for the new field of multimessenger astronomy.