Dynamical ejecta of neutron star mergers with nucleonic weak processes II: Kilonova emission

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo.Comment: 23 pages, 15 figures, 1 table, accepted to MNRAS; updated for correct version of Fig. A1, left panel, and corrected Eqs. 7 and

Self-consistent 3D radiative transfer for kilonovae: directional spectra from merger simulations

We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $\alpha$ and $\beta^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned treatments, a line-by-line treatment enables us include fluorescence effects and associate spectral features with the emitting and absorbing lines of individual elements. We find variations in the synthetic observables with both the polar and azimuthal viewing angles. The spectra exhibit blended features with strong interactions by Ce III, Sr II, Y II, and Zr II that vary with time and viewing direction. We demonstrate the importance of wavelength-calibration of atomic data using a model with calibrated Sr, Y, and Zr data, and find major differences in the resulting spectra, including a better agreement with AT2017gfo. The synthetic spectra for near-polar inclination show a feature at around 8000 A, similar to AT2017gfo. However, they evolve on a more rapid timescale, likely due to the low ejecta mass (0.005 M$_\odot$) as we take into account only the early ejecta. The comparatively featureless spectra for equatorial observers gives a tentative prediction that future observations of edge-on kilonovae will appear substantially different from AT2017gfo. We also show that 1D models obtained by spherically averaging the 3D ejecta lead to dramatically different direction-integrated luminosities and spectra compared to full 3D calculations.Comment: 12 pages, 5 figures. Accepted by ApJ

Quark-Novae, cosmic reionization, and early r-process element production

We examine the case for Quark-Novae (QNe) as possible sources for the reionization and early metal enrichment of the universe. Quark-Novae are predicted to arise from the explosive collapse (and conversion) of sufficiently massive neutron stars into quark stars. A Quark-Nova (QN) can occur over a range of time scales following the supernova event. For QNe that arise days to weeks after the supernovae, we show that dual-shock that arises as the QN ejecta encounter the supernova ejecta can produce enough photons to reionize hydrogen in most of the Inter-Galactic medium (IGM) by z ~ 6. Such events can explain the large optical depth tau_e ~ 0.1 as measured by WMAP, if the clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a normal initial mass function (IMF) for the oldest stars can be reconciled with a large optical depth as well as the mean metallicity of the early IGM post reionization. We find that QN also make a contribution to r-process element abundances for atomic numbers A > 130. We predict that the main cosmological signatures of Quark-Novae are the gamma-ray bursts that announce their birth. These will be clustered at redshifts in the range z ~ 7-8 in our model.Comment: Accepted for publication in the Astrophysical Journal. 9 journal pages, 2 figures. This version includes more discussion, 3 new appendices and extended literatur

Towards inferring the geometry of kilonovae

Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent).Comment: 11 pages, submitted to MNRA

Strangeness in Astrophysics and Cosmology

Some recent developments concerning the role of strange quark matter for astrophysical systems and the QCD phase transition in the early universe are addressed. Causality constraints of the soft nuclear equation of state as extracted from subthreshold kaon production in heavy-ion collisions are used to derive an upper mass limit for compact stars. The interplay between the viscosity of strange quark matter and the gravitational wave emission from rotation-powered pulsars are outlined. The flux of strange quark matter nuggets in cosmic rays is put in perspective with a detailed numerical investigation of the merger of two strange stars. Finally, we discuss a novel scenario for the QCD phase transition in the early universe, which allows for a small inflationary period due to a pronounced first order phase transition at large baryochemical potential.Comment: 8 pages, invited talk given at the International Conference on Strangeness in Quark Matter (SQM2009), Buzios, Brasil, September 28 - October 2, 200

Identification of strontium in the merger of two neutron stars.

Half of all of the elements in the Universe that are heavier than iron were created by rapid neutron capture. The theory underlying this astrophysical r-process was worked out six decades ago, and requires an enormous neutron flux to make the bulk of the elements1. Where this happens is still debated2. A key piece of evidence would be the discovery of freshly synthesized r-process elements in an astrophysical site. Existing models3-5 and circumstantial evidence6 point to neutron-star mergers as a probable r-process site; the optical/infrared transient known as a 'kilonova' that emerges in the days after a merger is a likely place to detect the spectral signatures of newly created neutron-capture elements7-9. The kilonova AT2017gfo-which was found following the discovery of the neutron-star merger GW170817 by gravitational-wave detectors10-was the first kilonova for which detailed spectra were recorded. When these spectra were first reported11,12, it was argued that they were broadly consistent with an outflow of radioactive heavy elements; however, there was no robust identification of any one element. Here we report the identification of the neutron-capture element strontium in a reanalysis of these spectra. The detection of a neutron-capture element associated with the collision of two extreme-density stars establishes the origin of r-process elements in neutron-star mergers, and shows that neutron stars are made of neutron-rich matter13

Challenges and opportunities of gravitational-wave searches at MHz to GHz frequencies

The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational-wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts that have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop “Challenges and opportunities of high-frequency gravitational wave detection” held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative

Evidence for quark-matter cores in massive neutron stars

The theory governing the strong nuclear force-quantum chromodynamics-predicts that at sufficiently high energy densities, hadronic nuclear matter undergoes a deconfinement transition to a new phase of quarks and gluons(1). Although this has been observed in ultrarelativistic heavy-ion collisions(2,3), it is currently an open question whether quark matter exists inside neutron stars(4). By combining astrophysical observations and theoretical ab initio calculations in a model-independent way, we find that the inferred properties of matter in the cores of neutron stars with mass corresponding to 1.4 solar masses (M-circle dot) are compatible with nuclear model calculations. However, the matter in the interior of maximally massive stable neutron stars exhibits characteristics of the deconfined phase, which we interpret as evidence for the presence of quark-matter cores. For the heaviest reliably observed neutron stars(5,6) with mass M approximate to 2M(circle dot), the presence of quark matter is found to be linked to the behaviour of the speed of sound c(s) in strongly interacting matter. If the conformal bound cs2Peer reviewe

Exploring the sensitivity of gravitational wave detectors to neutron star physics

The physics of neutron stars can be studied with gravitational waves emitted from coalescing binary systems. Tidal effects become significant during the last few orbits and can be visible in the gravitational-wave spectrum above 500 Hz. After the merger, the neutron star remnant oscillates at frequencies above 1 kHz and can collapse into a black hole. Gravitational-wave detectors with a sensitivity of ~10^{-24} strain/sqHz at 2-4 kHz can observe these oscillations from a source which is ~100 Mpc away. The current observatories, such as LIGO and Virgo, are limited by shot noise at high frequencies and have a sensitivity of > 2 * 10^{-23} strain/sqHz at 3 kHz. In this paper, we propose an optical configuration of gravitational-wave detectors which can be set up in present facilities using the current interferometer topology. This scheme has a potential to reach 7 * 10^{-25} strain/sqHz at 2.5 kHz without compromising the detector sensitivity to black hole binaries. We argue that the proposed instruments have a potential to detect similar amount of post-merger neutron star oscillations as the next generation detectors, such as Cosmic Explorer and Einstein Telescope. We also optimise the arm length of the future detectors for neutron star physics and find that the optimal arm length is ~20 km. These instruments have the potential to observe neutron star post-merger oscillations at a rate of ~30 events per year with a signal-to-noise ratio of 5 or more

Horizons: nuclear astrophysics in the 2020s and beyond

Nuclear astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities