Search

The LIGO-Virgo-KAGRA (LVK) Collaboration has made breakthrough discoveries in gravitational-wave astronomy, a new field that provides a different means of observing our Universe. Gravitational-wave discoveries are possible thanks to the work of thousands of people from across the globe working together. In this article, we discuss the range of engagement activities used to communicate LVK gravitational-wave discoveries and the stories of the people behind the science, using the activities surrounding the release of the third Gravitational-Wave Transient Catalog as a case study., Comment: 15 pages, 7 figures, published in JCOM: https://doi.org/10.22323/2.23070803

A key component of automated algorithm selection and configuration, which in most cases are performed using supervised machine learning (ML) methods is a good-performing predictive model. The predictive model uses the feature representation of a set of problem instances as input data and predicts the algorithm performance achieved on them. Common machine learning models struggle to make predictions for instances with feature representations not covered by the training data, resulting in poor generalization to unseen problems. In this study, we propose a workflow to estimate the generalizability of a predictive model for algorithm performance, trained on one benchmark suite to another. The workflow has been tested by training predictive models across benchmark suites and the results show that generalizability patterns in the landscape feature space are reflected in the performance space., Comment: To appear at GECCO 2023

Fair algorithm evaluation is conditioned on the existence of high-quality benchmark datasets that are non-redundant and are representative of typical optimization scenarios. In this paper, we evaluate three heuristics for selecting diverse problem instances which should be involved in the comparison of optimization algorithms in order to ensure robust statistical algorithm performance analysis. The first approach employs clustering to identify similar groups of problem instances and subsequent sampling from each cluster to construct new benchmarks, while the other two approaches use graph algorithms for identifying dominating and maximal independent sets of nodes. We demonstrate the applicability of the proposed heuristics by performing a statistical performance analysis of five portfolios consisting of three optimization algorithms on five of the most commonly used optimization benchmarks. The results indicate that the statistical analyses of the algorithms' performance, conducted on each benchmark separately, produce conflicting outcomes, which can be used to give a false indication of the superiority of one algorithm over another. On the other hand, when the analysis is conducted on the problem instances selected with the proposed heuristics, which uniformly cover the problem landscape, the statistical outcomes are robust and consistent., Comment: 10 pages, 6 figures

After their successful first observing run (September 12, 2015 - January 12, 2016), the Advanced LIGO detectors were upgraded to increase their sensitivity for the second observing run (November 30, 2016 - August 26, 2017). The Advanced Virgo detector joined the second observing run on August 1, 2017. We discuss the updates that happened during this period in the GstLAL-based inspiral pipeline, which is used to detect gravitational waves from the coalescence of compact binaries both in low latency and an offline configuration. These updates include deployment of a zero-latency whitening filter to reduce the over-all latency of the pipeline by up to 32 seconds, incorporation of the Virgo data stream in the analysis, introduction of a single-detector search to analyze data from the periods when only one of the detectors is running, addition of new parameters to the likelihood ratio ranking statistic, increase in the parameter space of the search, and introduction of a template mass-dependent glitch-excision thresholding method., Comment: 12 pages, 7 figures, to be submitted to Phys. Rev. D, comments welcome

We describe the methods used to construct the aligned-spin template bank of gravitational waveforms used by the GstLAL-based inspiral pipeline to analyze data from the second observing run of Advanced LIGO and the first observing run of advanced Virgo. The bank expands upon the parameter space covered during the first observing run, including coverage for merging compact binary systems with total mass between 2 $\mathrm{M}_{\odot}$ and 400 $\mathrm{M}_{\odot}$ and mass ratios between 1 and 97.988. Thus the systems targeted include merging neutron star-neutron star systems, neutron star-black hole binaries, and black hole-black hole binaries expanding into the intermediate-mass range. Component masses less than 2 $\mathrm{M}_{\odot}$ have allowed (anti-)aligned spins between $\pm0.05$ while component masses greater than 2 $\mathrm{M}_{\odot}$ have allowed (anti-)aligned between $\pm0.999$. The bank placement technique combines a stochastic method with a new grid-bank method to better isolate noisy templates, resulting in a total of 677,000 templates., Comment: 9 pages, 13 figures

We describe a stream-based analysis pipeline to detect gravitational waves from the merger of binary neutron stars, binary black holes, and neutron-star-black-hole binaries within ~ 1 minute of the arrival of the merger signal at Earth. Such low-latency detection is crucial for the prompt response by electromagnetic facilities in order to observe any fading electromagnetic counterparts that might be produced by mergers involving at least one neutron star. Even for systems expected not to produce counterparts, low-latency analysis of the data is useful for deciding when not to point telescopes, and as feedback to observatory operations. Analysts using this pipeline were the first to identify GW151226, the second gravitational-wave event ever detected. The pipeline also operates in an offline mode, in which it incorporates more refined information about data quality and employs acausal methods that are inapplicable to the online mode. The pipeline's offline mode was used in the detection of the first two gravitational-wave events, GW150914 and GW151226, as well as the identification of a third candidate, LVT151012.

We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, $n=1$ polytropes and are magnetized. We explore two types of magnetic-field geometries: one where each star is endowed with a dipole magnetic field extending from the interior into the exterior, as in a pulsar, and the other where the dipole field is initially confined to the interior. In both cases the adopted magnetic fields are initially dynamically unimportant. The merger outcome is a hypermassive neutron star that undergoes delayed collapse to a black hole (spin parameter $a/M_{\rm BH} \sim 0.74$) immersed in a magnetized accretion disk. About $4000M \sim 60(M_{\rm NS}/1.625M_\odot)$ ms following merger, the region above the black hole poles becomes strongly magnetized, and a collimated, mildly relativistic outflow --- an incipient jet --- is launched. The lifetime of the accretion disk, which likely equals the lifetime of the jet, is $\Delta t \sim 0.1 (M_{\rm NS}/1.625M_\odot)$ s. In contrast to black hole--neutron star mergers, we find that incipient jets are launched even when the initial magnetic field is confined to the interior of the stars., Comment: 6 pages, 3 figures, 1 table, matches published version

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime., Comment: 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

We derive the scalar waveform generated by a binary of nonspinning compact objects (black holes or neutron stars) in a general class of scalar-tensor theories of gravity. The waveform is accurate to 1.5 post-Newtonian order [$O((v/c)^3)$] beyond the leading-order tensor gravitational waves (the "Newtonian quadrupole"). To solve the scalar-tensor field equations, we adapt the direct integration of the relaxed Einstein equations formalism developed by Will, Wiseman, and Pati. The internal gravity of the compact objects is treated with an approach developed by Eardley. We find that the scalar waves are described by the same small set of parameters which describes the equations of motion and tensor waves. For black hole--black hole binaries, the scalar waveform vanishes, as expected from previous results which show that these systems in scalar-tensor theory are indistinguishable from their general relativistic counterparts. For black hole--neutron star binaries, the scalar waveform simplifies considerably from the generic case, essentially depending on only a single parameter up to first post-Newtonian order. With both the tensor and scalar waveforms in hand, we calculate the total energy flux carried by the outgoing waves. This quantity is computed to first post-Newtonian order relative to the "quadrupole formula" and agrees with previous, lower order calculations., Comment: 35 pages

General relativity (GR) has been extensively tested in the solar system and in binary pulsars, but never in the strong-field, dynamical regime. Soon, gravitational-wave (GW) detectors like Advanced LIGO and eLISA will be able to probe this regime by measuring GWs from inspiraling and merging compact binaries. One particularly interesting alternative to GR is scalar-tensor gravity. We present progress in the calculation of second post-Newtonian (2PN) gravitational waveforms for inspiraling compact binaries in a general class of scalar-tensor theories. The waveforms are constructed using a standard GR method known as "direct integration of the relaxed Einstein equations," appropriately adapted to the scalar-tensor case. We find that differences from general relativity can be characterized by a reasonably small number of parameters. Among the differences are new hereditary terms which depend on the past history of the source. In one special case, binary black hole systems, we find that the waveform is indistinguishable from that of general relativity. In another, mixed black hole-neutron star systems, all differences from GR can be characterized by only a single parameter., Comment: 7 pages, submitted to the proceedings of the 10th International LISA Symposium

We derive the tensor gravitational waveform generated by a binary of nonspinning compact objects (black holes or neutron stars) in a general class of scalar-tensor theories of gravity. The waveform is accurate to second post-Newtonian order beyond the leading order quadrupole approximation. We use the direct integration of the relaxed Einstein equations formalism, appropriately adapted to scalar-tensor theories, along with previous results for the equations of motion in these theories. The self-gravity of the compact objects is treated with an approach developed by Eardley. The scalar field causes deviations from the general relativistic waveform that depend only on a small number of parameters. Among the effects of the scalar field are new hereditary terms which depend on the past history of the source. One of these, a dipole-dipole coupling, produces a zero-frequency "gravitational-wave memory" equivalent to the Christodoulou memory of general relativity. In the special case of two black holes, the waveform reduces to the general relativistic waveform. For a mixed (black hole-neutron star) system, the waveform is identical to that of Einstein's theory to first post-Newtonian order, with deviations at higher order depending only on a single parameter. The behavior in these cases matches that found for the equations of motion., Comment: 44 pages

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime., Comment: 20 pages, 8 figures, proceedings of the 9th Amaldi Conference on Gravitational Waves. Final journal version. For a longer exposition of the eLISA science case, see http://arxiv.org/abs/1201.3621

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits., Comment: Submitted to GW Notes; abstract abridged due to arXiv requirements

We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning black hole binaries, combining for the first time the entire inspiral-merger-ringdown signal, the effect of the LISA orbits, and the complete three-channel LISA response. We consider an ensemble of systems near the peak of LISA's sensitivity band, with total rest mass of 2\times10^6 M\odot, a redshift of z = 1, and randomly chosen orientations and sky positions. We find median sky localization errors of approximately \sim3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, such as the James Webb Space Telescope, that could be used to search for electromagnetic signals associated with merging massive black holes. We investigate the way in which parameter errors decrease with measurement time, focusing specifically on the additional information provided during the merger-ringdown segment of the signal. We find that this information improves all parameter estimates directly, rather than through diminishing correlations with any subset of well- determined parameters. Although we have employed the baseline LISA design for this study, many of our conclusions regarding the information provided by mergers will be applicable to alternative mission designs as well., Comment: 9 pages, 5 figures, submitted to Phys. Rev. D

The future space-based gravitational wave detector LISA will be able to measure parameters of coalescing massive black hole binaries, often to extremely high accuracy. Previous work has demonstrated that the black hole spins can have a strong impact on the accuracy of parameter measurement. Relativistic spin-induced precession modulates the waveform in a manner which can break degeneracies between parameters, in principle significantly improving how well they are measured. Recent studies have indicated, however, that spin precession may be weak for an important subset of astrophysical binary black holes: those in which the spins are aligned due to interactions with gas. In this paper, we examine how well a binary's parameters can be measured when its spins are partially aligned and compare results using waveforms that include higher post-Newtonian harmonics to those that are truncated at leading quadrupole order. We find that the weakened precession can substantially degrade parameter estimation. This degradation is particularly devastating for the extrinsic parameters sky position and distance. Absent higher harmonics, LISA typically localizes the sky position of a nearly aligned binary a factor of $\sim 6$ less accurately than for one in which the spin orientations are random. Our knowledge of a source's sky position will thus be worst for the gas-rich systems which are most likely to produce electromagnetic counterparts. Fortunately, higher harmonics of the waveform can make up for this degradation. By including harmonics beyond the quadrupole in our waveform model, we find that the accuracy with which most of the binary's parameters are measured can be substantially improved. In some cases, parameters can be measured as well in partially aligned binaries as they can be when the binary spins are random., Comment: 18 pages, 16 figures, version accepted by PRD (with improved distributions of partially aligned spins)

The early part of the gravitational wave signal of binary neutron star inspirals can potentially yield robust information on the nuclear equation of state. The influence of a star's internal structure on the waveform is characterized by a single parameter: the tidal deformability lambda, which measures the star's quadrupole deformation in response to the companion's perturbing tidal field. We calculate lambda for a wide range of equations of state and find that the value of lambda spans an order of magnitude for the range of equation of state models considered. An analysis of the feasibility of discriminating between neutron star equations of state with gravitational wave observations of the early part of the inspiral reveals that the measurement error in lambda increases steeply with the total mass of the binary. Comparing the errors with the expected range of lambda, we find that Advanced LIGO observations of binaries at a distance of 100 Mpc will probe only unusually stiff equations of state, while the proposed Einstein Telescope is likely to see a clean tidal signature., Comment: 12 pages, submitted to PRD

The LISA Parameter Estimation (LISAPE) Taskforce was formed in September 2007 to provide the LISA Project with vetted codes, source distribution models, and results related to parameter estimation. The Taskforce's goal is to be able to quickly calculate the impact of any mission design changes on LISA's science capabilities, based on reasonable estimates of the distribution of astrophysical sources in the universe. This paper describes our Taskforce's work on massive black-hole binaries (MBHBs). Given present uncertainties in the formation history of MBHBs, we adopt four different population models, based on (i) whether the initial black-hole seeds are small or large, and (ii) whether accretion is efficient or inefficient at spinning up the holes. We compare four largely independent codes for calculating LISA's parameter-estimation capabilities. All codes are based on the Fisher-matrix approximation, but in the past they used somewhat different signal models, source parametrizations and noise curves. We show that once these differences are removed, the four codes give results in extremely close agreement with each other. Using a code that includes both spin precession and higher harmonics in the gravitational-wave signal, we carry out Monte Carlo simulations and determine the number of events that can be detected and accurately localized in our four population models., Comment: 14 pages, 2 figures, 5 tables, minor changes to match version to be published in the proceedings of the 7th LISA Symposium. For more information see the Taskforce's wiki at http://www.tapir.caltech.edu/dokuwiki/lisape:home

The coalescence of massive black holes is one of the primary sources of gravitational waves (GWs) for LISA. Measurements of the GWs can localize the source on the sky to an ellipse with a major axis of a few tens of arcminutes to a few degrees, depending on source redshift, and a minor axis which is 2--4 times smaller. The distance (and thus an approximate redshift) can be determined to better than a per cent for the closest sources we consider, although weak lensing degrades this performance. It will be of great interest to search this three-dimensional `pixel' for an electromagnetic counterpart to the GW event. The presence of a counterpart allows unique studies which combine electromagnetic and GW information, especially if the counterpart is found prior to final merger of the holes. To understand the feasibility of early counterpart detection, we calculate the evolution of the GW pixel with time. We find that the greatest improvement in pixel size occurs in the final day before merger, when spin precession effects are maximal. The source can be localized to within 10 square degrees as early as a month before merger at $z = 1$; for higher redshifts, this accuracy is only possible in the last few days., Comment: 11 pages, 4 figures, version published in Classical and Quantum Gravity (special issue for proceedings of 7th International LISA Symposium)

Massive black hole binary coalescences are prime targets for space-based gravitational wave (GW) observatories such as {\it LISA}. GW measurements can localize the position of a coalescing binary on the sky to an ellipse with a major axis of a few tens of arcminutes to a few degrees, depending on source redshift, and a minor axis which is $2 - 4$ times smaller. Neglecting weak gravitational lensing, the GWs would also determine the source's luminosity distance to better than percent accuracy for close sources, degrading to several percent for more distant sources. Weak lensing cannot, in fact, be neglected and is expected to limit the accuracy with which distances can be fixed to errors no less than a few percent. Assuming a well-measured cosmology, the source's redshift could be inferred with similar accuracy. GWs alone can thus pinpoint a binary to a three-dimensional ``pixel'' which can help guide searches for the hosts of these events. We examine the time evolution of this pixel, studying it at merger and at several intervals before merger. One day before merger, the major axis of the error ellipse is typically larger than its final value by a factor of $\sim 1.5-6$. The minor axis is larger by a factor of $\sim 2-9$, and, neglecting lensing, the error in the luminosity distance is larger by a factor of $\sim 1.5-7$. This large change over a short period of time is due to spin-induced precession, which is strongest in the final days before merger. The evolution is slower as we go back further in time. For $z = 1$, we find that GWs will localize a coalescing binary to within $\sim 10\ \mathrm{deg}^2$ as early as a month prior to merger and determine distance (and hence redshift) to several percent., Comment: 30 pages, 10 figures, 5 tables. Version published in ApJ

The coalescence of massive black holes generates gravitational waves (GWs) that will be measurable by space-based detectors such as LISA to large redshifts. The spins of a binary's black holes have an important impact on its waveform. Specifically, geodetic and gravitomagnetic effects cause the spins to precess; this precession then modulates the waveform, adding periodic structure which encodes useful information about the binary's members. Following pioneering work by Vecchio, we examine the impact upon GW measurements of including these precession-induced modulations in the waveform model. We find that the additional periodicity due to spin precession breaks degeneracies among certain parameters, greatly improving the accuracy with which they may be measured. In particular, mass measurements are improved tremendously, by one to several orders of magnitude. Localization of the source on the sky is also improved, though not as much -- low redshift systems can be localized to an ellipse which is roughly ${a few} \times 10$ arcminutes in the long direction and a factor of 2-4 smaller in the short direction. Though not a drastic improvement relative to analyses which neglect spin precession, even modest gains in source localization will greatly facilitate searches for electromagnetic counterparts to GW events. Determination of distance to the source is likewise improved: We find that relative error in measured luminosity distance is commonly $\sim 0.2%-0.7%$ at $z \sim 1$. Finally, with the inclusion of precession, we find that the magnitude of the spins themselves can typically be determined for low redshift systems with an accuracy of about $0.1%-10 %$, depending on the spin value, allowing accurate surveys of mass and spin evolution over cosmic time., Comment: 22 pages, 10 figures. Corrects small typo in equations; results unaffected