Search

Gravitational-wave observations of binary neutron star mergers and their electromagnetic counterparts provide an independent measurement of the Hubble constant, $H_0$, through the standard-sirens approach. Current methods of determining $H_0$, such as measurements from the early universe and the local distance ladder, are in tension with one another. If gravitational waves are to break this tension a thorough understanding of systematic uncertainty for gravitational-wave observations is required. To accurately estimate the properties of gravitational-wave signals measured by LIGO and Virgo, we need to understand the characteristics of the detectors noise. Non-gaussian transients in the detector data and rapid changes in the instrument, known as non-stationary noise, can both add a systematic uncertainty to inferred results. We investigate how non-stationary noise affects the estimation of the luminosity distance of the source, and therefore of $H_0$. Using a population of 200 simulated binary neutron star signals, we show that non-stationary noise can bias the estimation of the luminosity distance by up to 2.4\%. However, only $\sim$15\% of binary neutron star signals would be affected around their merger time with non-stationary noise at a similar level to that seen in the first half of LIGO-Virgo's third observing run. Comparing the expected bias to other systematic uncertainties, we argue that non-stationary noise in the current generation of detectors will not be a limiting factor in resolving the tension on $H_0$ using standard sirens. Although, evaluating non-stationarity in gravitational-wave data will be crucial to obtain accurate estimates of $H_0$., Comment: 10 pages, 5 figures, 1 table

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate., Comment: Published by JOSS. The package is available at https://github.com/skypyproject/skypy. Comments, issues and pull requests are welcome

Gravitational waves from compact binaries measured by the LIGO and Virgo detectors are routinely analyzed using Markov Chain Monte Carlo sampling algorithms. Because the evaluation of the likelihood function requires evaluating millions of waveform models that link between signal shapes and the source parameters, running Markov chains until convergence is typically expensive and requires days of computation. In this extended abstract, we provide a proof of concept that demonstrates how the latest advances in neural simulation-based inference can speed up the inference time by up to three orders of magnitude -- from days to minutes -- without impairing the performance. Our approach is based on a convolutional neural network modeling the likelihood-to-evidence ratio and entirely amortizes the computation of the posterior. We find that our model correctly estimates credible intervals for the parameters of simulated gravitational waves., Comment: V1: First version; V2: Updated references; V3: Update references and camera-ready version; V4: Correct figure labels; V5: Updated references

The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system. In particular, constraining the inclination of the binary can reduce the distance uncertainty of the source. Here we present the first realistic set of localizations for binary black hole mergers, including different prior constraints on the binaries' inclinations. We find that prior information on the inclination can reduce the localization volume by a factor of 3. We discuss two astrophysical scenarios of interest: (i) follow-up searches for beamed electromagnetic/neutrino counterparts and (ii) mergers in the accretion disks of active galactic nuclei., Comment: Accepted for publication in MNRAS; 5 pages, 2 figures

Gravitational waveforms from numerical simulations are a critical tool to test and analytically calibrate the waveform models used to study the properties of merging compact objects. In this paper, we present a series of high-accuracy waveforms produced with the SpEC code for systems involving at least one neutron star. We provide for the first time waveforms with sub-radian accuracy over more than twenty cycles for low-mass black hole-neutron star binaries, including binaries with non-spinning objects, and binaries with rapidly spinning neutron stars that maximize the impact on the gravitational wave signal of the near-resonant growth of the fundamental excitation mode of the neutron star (f-mode). We also provide for the first time with SpEC a high-accuracy neutron star-neutron star waveform. These waveforms are made publicly available as part of the SxS catalogue. We compare our results to analytical waveform models currently implemented in data analysis pipelines. For most simulations, the models lie outside of the predicted numerical errors in the last few orbits before merger, but do not show systematic deviations from the numerical results: comparing different models appears to provide reasonable estimates of the modeling errors. The sole exception is the equal-mass simulation using a rapidly counter-rotating neutron star to maximize the impact of the excitation of the f-mode, for which all models perform poorly. This is however expected, as even the single model that takes f-mode excitation into account ignores the significant impact of the neutron star spin on the f-mode excitation frequency., Comment: 17p, 10 figures

Gravitational wave (GW) observations of binary neutron star (BNS) mergers can be used to measure luminosity distances and hence, when coupled with estimates for the mergers' host redshifts, infer the Hubble constant, $H_0$. These observations are, however, affected by GW measurement noise, uncertainties in host redshifts and peculiar velocities, and are potentially biased by selection effects and the mis-specification of the cosmological model or the BNS population. The estimation of $H_0$ from samples of BNS mergers with optical counterparts is tested here by using a phenomenological model for the GW strains that captures both the data-driven event selection and the distance-inclination degeneracy, while being simple enough to facilitate large numbers of simulations. A rigorous Bayesian approach to analyzing the data from such simulated BNS merger samples is shown to yield results that are unbiased, have the appropriate uncertainties, and are robust to model mis-specification. Applying such methods to a sample of $N \simeq 50$ BNS merger events, as LIGO+Virgo could produce in the next $\sim 5$ years, should yield robust and accurate Hubble constant estimates that are precise to a level of $\sim 2$ km s$^{-1}$ Mpc$^{-1}$, sufficient to reliably resolve the current tension between local and cosmological measurements of $H_0$., Comment: 21 pages (including eight-page appendices); 9 figures. Matches version accepted by PRD

Synchrotron-radiation angle-resolved and core-level photoemission spectroscopy are used together to investigate the origin of the two-dimensional electron gas on the surface of single crystal CdO (100) films. A reduction in the two-dimensional electron density of the surface state is observed under the synchrotron beam during ARPES, which is shown to be accompanied by a concomitant reduction in the surface-adsorbed species (monitored through the O 1s core level signal). This shows that surface adsorbates donate electrons into the surface accumulation layer. When the surface is cleaned, the surface conduction band state empties. A surface doped with atomic H is also studied. Here, interstitial H increases the two-dimensional electron density at the surface. This demonstrates that reversible donor doping is possible. The surface band bending profiles, 2D electron densities, and effective masses are calculated from subband dispersion simulations., Comment: 8 pages, 6 figures

The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations., Comment: 8 pages

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

The Hubble constant ($H_0$) estimated from the local Cepheid-supernova (SN) distance ladder is in 3-$\sigma$ tension with the value extrapolated from cosmic microwave background (CMB) data assuming the standard cosmological model. Whether this tension represents new physics or systematic effects is the subject of intense debate. Here, we investigate how new, independent $H_0$ estimates can arbitrate this tension, assessing whether the measurements are consistent with being derived from the same model using the posterior predictive distribution (PPD). We show that, with existing data, the inverse distance ladder formed from BOSS baryon acoustic oscillation measurements and the Pantheon SN sample yields an $H_0$ posterior near-identical to the Planck CMB measurement. The observed local distance ladder value is a very unlikely draw from the resulting PPD. Turning to the future, we find that a sample of $\sim50$ binary neutron star "standard sirens" (detectable within the next decade) will be able to adjudicate between the local and CMB estimates., Comment: Eight pages, four figures. v3: matches version accepted by Physical Review Letters. Code available at https://github.com/sfeeney/hh0

The convex feasibility problem consists in finding a point in the intersection of a finite family of closed convex sets. When the intersection is empty, a best compromise is to search for a point that minimizes the sum of the squared distances to the sets. In 2001, de Pierro conjectured that the limit cycles generated by the $\varepsilon$-under-relaxed cyclic projection method converge when $\varepsilon\downarrow 0$ towards a least squares solution. While the conjecture has been confirmed under fairly general conditions, we show that it is false in general by constructing a system of three compact convex sets in $\mathbb{R}^3$ for which the $\varepsilon$-under-relaxed cycles do not converge.

In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at multiple resolutions and with two independent techniques to solve Einstein's equations. For the nonprecessing and precessing simulations, we demonstrate the two techniques agree mode by mode, at a precision substantially in excess of statistical uncertainties in current LIGO's observations. Conversely, we demonstrate our full numerical solutions contain information which is not accurately captured with the approximate phenomenological models commonly used to infer compact binary parameters. To quantify the impact of these differences on parameter inference for GW170104 specifically, we compare the predictions of our simulations and these approximate models to LIGO's observations of GW170104., Comment: 11 figures, 20 pages

The first detection of a binary neutron star merger, GW170817, and an associated short gamma-ray burst confirmed that neutron star mergers are responsible for at least some of these bursts. The prompt gamma ray emission from these events is thought to be highly relativistically beamed. We present a method for inferring limits on the extent of this beaming by comparing the number of short gamma-ray bursts observed electromagnetically to the number of neutron star binary mergers detected in gravitational waves. We demonstrate that an observing run comparable to the expected Advanced LIGO 2016--2017 run would be capable of placing limits on the beaming angle of approximately $\theta \in (2.88^\circ,14.15^\circ)$, given one binary neutron star detection. We anticipate that after a year of observations with Advanced LIGO at design sensitivity in 2020 these constraints would improve to $\theta \in (8.10^\circ,14.95^\circ)$., Comment: 11 pages with 10 figures

Calculations within density functional theory on the basis of the local density approximation are carried out to study the properties of intrinsic point defects in aluminum antimonide. Special care is taken to address finite-size effects, band gap error, and symmetry reduction in the defect structures. The correction of the band gap is based on a set of GW calculations. The most important defects are identified to be the aluminum interstitial $Al_{i,Al}^{1+}$, the antimony antisites $Sb_{Al}^0$ and $Sb_{Al}^{1+}$, and the aluminum vacancy $V_{Al}^{3-}$. The intrinsic defect and charge carrier concentrations in the impurity-free material are calculated by self-consistently solving the charge neutrality equation. The impurity-free material is found to be n-type conducting at finite temperatures., Comment: 11 pages, 5 figures

We investigate the impact of water, a polar solvent, on the optical absorption of prototypical silicon clusters with oxygen passivation. We approach this complex problem by assessing the contributions of three factors: chemical reactivity; thermal equilibration and dielectric screening. We find that the silanone (Si=O) functional group is not chemically stable in the presence of water and exclude this as a source of significant red shift in absorption in aqueous environments. We perform first principles molecular dynamics simulations of the solvation of an oxygenated silicon cluster with explicit water molecules at 300 K. We find a systematic 0.7 eV red shift in the absorption gap of this cluster, which we attribute to thermal strain of the molecular structure. Surprisingly, we find no observable screening impact of the solvent, in contrast with consistent blue shifts observed for similarly sized organic molecules in polar solvents. The predicted red shift is expected to be significantly smaller for larger Si quantum dots produced experimentally, guaranteeing that their vacuum optical properties are preserved even in aqueous environments., Comment: 11 pages, 7 figures, 3 tables. submitted to JACS

We evaluate the optical gap and Stokes shift of several candidate 1 nm silicon nanocrystal structures using density functional and quantum Monte Carlo (QMC) methods. We find that the combination of absorption gap calculations and Stokes shift calculations may be used to determine structures. We find that although absorption gaps calculated within B3LYP and QMC agree for spherical, completely hydrogenated silicon nanocrystals, they disagree in clusters with different surface bonding networks. The nature of the Stokes shift of the ultrabright luminescence is examined by comparing possible relaxation mechanisms. We find that the exciton which reproduces the experimental value of the Stokes shift is most likely a state formed by a collective structural relaxation distributed over the entire cluster.