We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational-wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between ...1.5 and 9.5 days post-merger, using the Southern Astrophysical Research and Magellan telescopes; the UV spectrum was obtained with the Hubble Space Telescope at 5.5 days. Our data reveal a rapidly fading blue component ( T 5500 K at 1.5 days) that quickly reddens; spectra later than 4.5 days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at ∼7900 at t 4.5 days. The colors, rapid evolution, and featureless spectrum are consistent with a "blue" kilonova from polar ejecta comprised mainly of light r-process nuclei with atomic mass number A 140 . This indicates a sightline within θ obs 45 ° of the orbital axis. Comparison to models suggests ∼0.03 M of blue ejecta, with a velocity of ∼ 0.3 c . The required lanthanide fraction is ∼ 10 − 4 , but this drops to < 10 − 5 in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of 12 km. This mass also supports the idea that neutron star mergers are a major contributor to r-process nucleosynthesis.
We present Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA) radio observations of GW170817, the first Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo ...gravitational wave (GW) event from a binary neutron star merger and the first GW event with an electromagnetic (EM) counterpart. Our data include the first observations following the discovery of the optical transient at both the centimeter (13.7 hr post-merger) and millimeter (2.41 days post-merger) bands. We detect faint emission at 6 GHz at 19.47 and 39.23 days after the merger, but not in an earlier observation at 2.46 days. We do not detect cm/mm emission at the position of the optical counterpart at frequencies of 10-97.5 GHz at times ranging from 0.6 to 30 days post-merger, ruling out an on-axis short gamma-ray burst (SGRB) for energies 10 48 erg. For fiducial SGRB parameters, our limits require an observer viewer angle of 20°. The radio and X-ray data can be jointly explained as the afterglow emission from an SGRB with a jet energy of ∼ 10 49 - 10 50 erg that exploded in a uniform density environment with n ∼ 10 − 4 - 10 − 2 cm−3, viewed at an angle of ∼20°-40° from the jet axis. Using the results of our light curve and spectral modeling, in conjunction with the inference of the circumbinary density, we predict the emergence of late-time radio emission from the deceleration of the kilonova (KN) ejecta on a timescale of ∼5-10 years that will remain detectable for decades with next-generation radio facilities, making GW170817 a compelling target for long-term radio monitoring.
Abstract
A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining ...cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models—including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is 2.5× larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range 2000–11,000 Å (1800 Å redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-
z I
and
iz
photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-
z
Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties, we compare distances of low (0.01 <
z
< 0.2) and high (0.4 <
z
< 0.6) redshift SNe in the training compilation, finding an insignificant 3 ± 14 mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at
https://saltshaker.readthedocs.io/en/latest/
, and the model is integrated into the
sncosmo
and
SNANA
software packages.
We present the properties of NGC 4993, the host galaxy of GW170817, the first gravitational-wave (GW) event from the merger of a binary neutron star (BNS) system and the first with an electromagnetic ...(EM) counterpart. We use both archival photometry and new optical/near-IR imaging and spectroscopy, together with stellar population synthesis models to infer the global properties of the host galaxy. We infer a star formation history peaked at ago, with subsequent exponential decline leading to a low current star formation rate of 0.01 yr−1, which we convert into a binary merger timescale probability distribution. We find a median merger timescale of Gyr, with a 90% confidence range of . This in turn indicates an initial binary separation of , comparable to the inferred values for Galactic BNS systems. We also use new and archival Hubble Space Telescope images to measure a projected offset of the optical counterpart of 2.1 kpc (0.64re) from the center of NGC 4993 and to place a limit of mag on any pre-existing emission, which rules out the brighter half of the globular cluster luminosity function. Finally, the age and offset of the system indicates it experienced a modest natal kick with an upper limit of ∼200 km s−1. Future GW-EM observations of BNS mergers will enable measurement of their population delay time distribution, which will directly inform their viability as the dominant source of r-process enrichment in the universe.
We present a comprehensive comparison of the properties of the radio through X-ray counterpart of GW170817 and the properties of short-duration gamma-ray bursts (GRBs). For this effort, we utilize a ...sample of 36 short GRBs spanning a redshift range of z 0.12 - 2.6 discovered over 2004-2017. We find that the counterpart to GW170817 has an isotropic-equivalent luminosity that is 3000 times less than the median value of on-axis short GRB X-ray afterglows, and 104 times less than that for detected short GRB radio afterglows. Moreover, the allowed jet energies and particle densities inferred from the radio and X-ray counterparts to GW170817 and on-axis short GRB afterglows are remarkably similar, suggesting that viewing angle effects are the dominant, and perhaps only, difference in their observed radio and X-ray behavior. From comparison to previous claimed kilonovae following short GRBs, we find that the optical and near-infrared (NIR) counterpart to GW170817 is comparatively under-luminous by a factor of 3 - 5 , indicating a range of kilonova luminosities and timescales. A comparison of the optical limits following short GRBs on 1 day timescales also rules out a "blue" kilonova of comparable optical isotropic-equivalent luminosity in one previous short GRB. Finally, we investigate the host galaxy of GW170817, NGC 4993, in the context of short GRB host galaxy stellar population properties. We find that NGC 4993 is superlative in terms of its large luminosity, old stellar population age, and low star formation rate compared to previous short GRB hosts. Additional events within the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo volume will be crucial in delineating the properties of the host galaxies of neutron star-neutron star (NS-NS) mergers, and connecting them to their cosmological counterparts.
The observed tension (∼9% difference) between the local distance ladder measurement of the Hubble constant, H0, and its value inferred from the cosmic microwave background could hint at new, exotic, ...cosmological physics. We test the impact of the assumption about the expansion history of the universe ( ) on the local distance ladder estimate of H0. In the fiducial analysis, the Hubble flow Type Ia supernova (SN Ia) sample is truncated to z < 0.15, and the deceleration parameter (q0) is fixed to −0.55. We create realistic simulations of the calibrator and Pantheon samples, and account for a full systematics covariance between these two sets. We fit several physically motivated dark-energy models, and derive combined constraints from calibrator and Pantheon SNe Ia and simultaneously infer H0 and dark-energy properties. We find that the assumption on the dark-energy model does not significantly change the local distance ladder value of H0, with a maximum difference (ΔH0) between the inferred value for different models of 0.47 km , i.e., a 0.6% shift in H0, significantly smaller than the observed tension. Additional freedom in the dark-energy models does not increase the error in the inferred value of H0. Including systematics covariance between the calibrators, low-redshift SNe, and high-redshift SNe can induce small shifts in the inferred value for H0. The SN Ia systematics in this study contribute 0.8% to the total uncertainty of H0.
ABSTRACT
We present a revised Spectral Adaptive Light Curve Template (SALT2) surface (SALT2-2021) for fitting the light curves of Type Ia supernovae (SNe Ia), which incorporates new measurements of ...zero-point calibration offsets and Milky Way extinction. The most notable change in the new surface occurs in the UV region. This new surface alters the distance measurements of SNe Ia, which can be used to investigate the nature of dark energy by probing the expansion history of the Universe. Using the revised SALT2 surface on public data from the first 3 yr of the Dark Energy Survey Supernova Program (combined with an external low-z SNe Ia sample) and combining with cosmic microwave background constraints, we find a change in the dark energy equation-of-state parameter, Δw = 0.015 ± 0.004. This result highlights the continued importance of controlling and reducing systematic uncertainties, particularly with the next generation of supernova analyses aiming to improve constraints on dark energy properties.
We present a near-infrared spectral sequence of the electromagnetic counterpart to the binary neutron star merger GW170817 detected by Advanced Laser Interferometer Gravitational-wave Observatory ...(LIGO)/Virgo. Our data set comprises seven epochs of J+H spectra taken with FLAMINGOS-2 on Gemini-South between 1.5 and 10.5 days after the merger. In the initial epoch, the spectrum is dominated by a smooth blue continuum due to a high-velocity, lanthanide-poor blue kilonova component. Starting the following night, all of the subsequent spectra instead show features that are similar to those predicted in model spectra of material with a high concentration of lanthanides, including spectral peaks near 1.07 and 1.55 m. Our fiducial model with 0.04 M of ejecta, an ejection velocity of v = 0.1c, and a lanthanide concentration of Xlan = 10−2 provides a good match to the spectra taken in the first five days, although it over-predicts the late-time fluxes. We also explore models with multiple fitting components, in each case finding that a significant abundance of lanthanide elements is necessary to match the broad spectral peaks that we observe starting at 2.5 days after the merger. These data provide direct evidence that binary neutron star mergers are significant production sites of even the heaviest r-process elements.
ABSTRACT
Type Ia supernovae (SNe Ia) are powerful tools for measuring the expansion history of the Universe, but the impact of dust around SNe Ia remains unknown and is a critical systematic ...uncertainty. One way to improve our empirical description of dust is to analyse highly reddened SNe Ia E(B − V) > 0.4, roughly equivalent to the fitted SALT2 light-curve parameter c > 0.3. With the recently released Pantheon+ sample, there are 57 SNe Ia that were removed because of their high colour alone (with colours up to c = 1.61), which can provide enormous leverage on understanding line-of-sight RV. Previous studies have claimed that RV decreases with redder colour, though it is unclear if this is due to limited statistics, selection effects, or an alternative explanation. To test this claim, we fit two separate colour–luminosity relationships, one for the main cosmological sample (c < 0.3) and one for highly reddened (c > 0.3) SNe Ia. We find the change in the colour–luminosity coefficient to be consistent with zero. Additionally, we compare the data to simulations with different colour models, and find that the data prefers a model with a flat dependence of RV on colour over a declining dependence. Finally, our results strongly support that line-of-sight RV to SNe Ia is not a single value, but forms a distribution.
We present griz sub(P1) light curves of 146 spectroscopically confirmed Type Ia supernovae (SNe Ia; 0.03 < z < 0.65) discovered during the first 1.5 yr of the Pan-STARRS1 Medium Deep Survey. The ...Pan-STARRS1 natural photometric system is determined by a combination of on-site measurements of the instrument response function and observations of spectrophotometric standard stars. We find that the systematic uncertainties in the photometric system are currently 1.2% without accounting for the uncertainty in the Hubble Space Telescope Calspec definition of the AB system. A Hubble diagram is constructed with a subset of 113 out of 146 SNe Ia that pass our light curve quality cuts. The cosmological fit to 310 SNe Ia (113 PSI SNe Ia + 222 light curves from 197 low-z SNe Ia), using only supernovae (SNe) and assuming a constant dark energy equation of state and flatness, yields w = -1.120 super(+0.360) sub(-0.206)(Stat) super(+0.269) sub(-0.291)(Sys). When combined with BAO+CMB(Planck)+H sub(0), the analysis yields Omega sub(M) = 0.280 super(0.013) sub(-0.012) and w = 1.166 super(+0.072) sub(-0.069) including all identified systematics. The value of w is inconsistent with the cosmological constant value of -1 at the 2.3sigma level. Tension endures after removing either the baryon acoustic oscillation (BAO) or the H sub(0) constraint, though it is strongest when including the H sub(0) constraint. If we include WMAP9 cosmic microwave background (CMB) constraints instead of those from Planck, we find w = -1.124 super(+0.083) sub(-0.065) which diminishes the discord to <2sigma. We cannot conclude whether the tension with flat ACDM is a feature of dark energy, new physics, or a combination of chance and systematic errors. The full Pan-STARRS1 SN sample with ~three times as many SNe should provide more conclusive results.