The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

Ten binary black-hole mergers have already been detected during the first two observing runs of advanced LIGO and Virgo, and many more are expected to be observed in the near future. This opens the possibility for gravitational-wave astronomy to better constrain the properties of black hole binaries, not only as single sources, but as a whole astrophysical population. In this paper, we address the problem of using gravitational-wave measurements to estimate the proportion of merging black holes produced either via isolated binaries or binaries evolving in young star clusters. To this end, we use a Bayesian hierarchical modeling approach applied to catalogs of merging binary black holes generated using state-of-the-art population synthesis and N-body codes. In particular, we show that, although current advanced LIGO/Virgo observations only mildly constrain the mixing fraction $f \in [0,1]$ between the two formation channels, we expect to narrow down the fractional errors on $f$ to $10-20\%$ after a few hundreds of detections.

Light bosons, proposed as a possible solution to various problems in fundamental physics and cosmology, include a broad class of candidates for beyond the Standard Model physics, such as dilatons and moduli, wave dark matter and axion-like particles. If light bosons exist in nature, they will spontaneously form "clouds" by extracting rotational energy from rotating massive black holes through superradiance, a classical wave amplification process that has been studied for decades. The superradiant growth of the cloud sets the geometry of the final black hole, and the black hole geometry determines the shape of the cloud. Hence, both the black hole geometry and the cloud encode information about the light boson. For this reason, measurements of the gravitational field of the black hole/cloud system (as encoded in gravitational waves) are over-determined. We show that a single gravitational wave measurement can be used to verify the existence of light bosons by model selection, rule out alternative explanations for the signal, and measure the boson mass. Such measurements can be done generically for bosons in the mass range $[10^{-16.5},10^{-14}]$ eV using LISA observations of extreme mass-ratio inspirals.

Two of the dominant channels to produce merging stellar-mass black-hole binaries are believed to be the isolated evolution of binary stars in the field and dynamical formation in star clusters. The first black-hole binary event from the third LIGO/Virgo observing run (GW190412) is unusual in that it has unequal masses, nonzero effective spin, and nonzero primary spin at 90% confidence interval. We show that this event should be exceedingly rare in the context of both the field and cluster formation scenarios. Interpreting GW190412 as a remnant of a previous black-hole merger provides a promising route to explain its features. If GW190412 indeed formed hierarchically, we show that the region of the parameter space which is best motivated from an astrophysical standpoint (low natal spins and light clusters) cannot accommodate the observation. We analyze public GW190412 LIGO/Virgo data with a Bayesian prior where the more massive black hole resulted from a previous merger, and find that this scenario is equally supported by the data. If the heavier component of GW190412 is indeed a merger remnant, then its spin magnitude is $\chi_1=0.56_{-0.21}^{+0.19}$, higher than the value previously reported by the LIGO/Virgo collaboration.

Binary neutron-star mergers will predominantly produce black-hole remnants of mass $\sim 3-4\,M_{\odot}$, thus populating the putative \emph{low mass gap} between neutron stars and stellar-mass black holes. If these low-mass black holes are in dense astrophysical environments, mass segregation could lead to "second-generation" compact binaries merging within a Hubble time. In this paper, we investigate possible signatures of such low-mass compact binary mergers in gravitational-wave observations. We show that this unique population of objects, if present, will be uncovered by the third-generation gravitational-wave detectors, such as Cosmic Explorer and Einstein Telescope. Future joint measurements of chirp mass ${\cal M}$ and effective spin $\chi_{\rm eff}$ could clarify the formation scenario of compact objects in the low mass gap. As a case study, we show that the recent detection of GW190425 (along with GW170817) favors a double Gaussian mass model for neutron stars, under the assumption that the primary in GW190425 is a black hole formed from a previous binary neutron star merger.

Gravitational waves from the merger of two neutron stars cannot be easily distinguished from those produced by a comparable-mass mixed binary in which one of the companions is a black hole. Low-mass black holes are interesting because they could form in the aftermath of the coalescence of two neutron stars, from the collapse of massive stars, from matter overdensities in the primordial Universe, or as the outcome of the interaction between neutron stars and dark matter. Gravitational waves carry the imprint of the internal composition of neutron stars via the so-called tidal deformability parameter, which depends on the stellar equation of state and is equal to zero for black holes. We present a new data analysis strategy powered by Bayesian inference and machine learning to identify mixed binaries, hence low-mass black holes, using the distribution of the tidal deformability parameter inferred from gravitational-wave observations.

Two of the dominant channels to produce the black-hole binary mergers observed by LIGO and Virgo are believed to be the isolated evolution of stellar binaries in the field and dynamical formation in star clusters. Their relative efficiency can be characterized by a "mixing fraction." Pair instabilities prevent stellar collapse from generating black holes more massive than about $45 M_\odot$. This "mass gap" only applies to the field formation scenario, and it can be filled by repeated mergers in clusters. A similar reasoning applies to the binary's effective spin. If black holes are born slowly rotating, the high-spin portion of the parameter space (the "spin gap") can only be populated by black hole binaries that were assembled dynamically. Using a semianalytical cluster model, we show that future gravitational-wave events in either the mass gap, the spin gap, or both can be leveraged to infer the mixing fraction between the field and cluster formation channels.

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

LISA can detect higher harmonics of the ringdown gravitational-wave signal from massive black-hole binary mergers with large signal-to-noise ratio. The most massive black-hole binaries are more likely to have electromagnetic counterparts, and the inspiral will contribute little to their signal-to-noise ratio. Here we address the following question: can we extract the binary parameters and localize the source using LISA observations of the ringdown only? Modulations of the amplitude and phase due to LISA's motion around the Sun can be used to disentangle the source location and orientation when we detect the long-lived inspiral signal, but they can not be used for ringdown-dominated signals, which are very short-lived. We show that (i) we can still measure the mass ratio and inclination of high-mass binaries by carefully combining multiple ringdown harmonics, and (ii) we can constrain the sky location and luminosity distance by relying on the relative amplitudes and phases of various harmonics, as measured in different LISA channels.

Black-hole spectroscopy is arguably the most promising tool to test gravity in extreme regimes and to probe the ultimate nature of black holes with unparalleled precision. These tests are currently limited by the lack of a ringdown parametrization that is both robust and accurate. We develop an observable-based parametrization of the ringdown of spinning black holes beyond general relativity, which we dub ParSpec (Parametrized Ringdown Spin Expansion Coefficients). This approach is perturbative in the spin, but it can be made arbitrarily precise (at least in principle) through a high-order expansion. It requires O(10) ringdown detections, which should be routinely available with the planned space mission LISA and with third-generation ground-based detectors. We provide a preliminary analysis of the projected bounds on parametrized ringdown parameters with LISA and with the Einstein Telescope, and discuss extensions of our model that can be straightforwardly included in the future.

We present new numerical-relativity simulations of eccentric merging black holes with initially antiparallel spins lying in the orbital plane (the so-called \emph{superkick} configuration). Binary eccentricity boosts the recoil of the merger remnant by up to $25\%$. The increase in the energy flux is much more modest, and therefore this kick enhancement is mainly due to asymmetry in the binary dynamics. Our findings might have important consequences for the retention of stellar-mass black holes in star clusters and supermassive black holes in galactic hosts.

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

One of the Holy Grails of observational astronomy is to confirm the prediction that black holes in the Universe are described by the Kerr solution of Einstein's field equations of general relativity. This Topical Collection provides a status report of theoretical and experimental progress towards confirming the "Kerr paradigm" through X-ray astronomy, gravitational lensing, stellar tidal disruption events, superradiance, and gravitational-wave observations of black hole binary mergers.

Using simulations performed with the population synthesis code MOBSE, we compute the merger rate densities and detection rates of compact binary mergers formed in isolation for second- and third-generation gravitational-wave detectors. We estimate how rates are affected by uncertainties on key stellar-physics parameters, namely common envelope evolution and natal kicks. We estimate how future upgrades will increase the size of the available catalog of merger events, and we discuss features of the merger rate density that will become accessible with third-generation detectors.

The space-based detector LISA may observe gravitational waves from the early inspiral of stellar-mass black hole binaries, some of which could have significant eccentricity. Current gravitational waveform templates are only valid for small orbital velocities (i.e., in a post-Newtonian expansion) and small initial eccentricity $e_0$ ("post-circular" expansion). We conventionally define $e_0$ as the eccentricity corresponding to an orbital frequency of $5 \text{ mHz}$, and we study the convergence properties of frequency-domain inspiral templates that are accurate up to 2PN and order $e_0^6$ in eccentricity. We compute the so-called "unfaithfulness" between the full template and "reduced" templates obtained by dropping some terms in the phasing series; we investigate the conditions under which systematic errors are negligible with respect to statistical errors, and we study the convergence properties of statistical errors. In general, eccentric waveforms lead to larger statistical errors than circular waveforms due to correlations between the parameters, but the error estimates do not change significantly as long as we include terms of order $e_0^2$ or higher in the phasing.

We study the post-Newtonian dynamics of black hole binaries in Einstein-scalar-Gauss-Bonnet gravity theories. To this aim we build static, spherically symmetric black hole solutions at fourth order in the Gauss-Bonnet coupling $\alpha$. We then "skeletonize" these solutions by reducing them to point particles with scalar field-dependent masses, showing that this procedure amounts to fixing the Wald entropy of the black holes during their slow inspiral. The cosmological value of the scalar field plays a crucial role in the dynamics of the binary. We compute the two-body Lagrangian at first post-Newtonian order and show that no regularization procedure is needed to obtain the Gauss-Bonnet contributions to the fields, which are finite. We illustrate the power of our approach by Pad\'e-resumming the so-called "sensitivities," which measure the coupling of the skeletonized body to the scalar field, for some specific theories of interest.

Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

Unlike traditional electromagnetic measurements, gravitational-wave observations are not affected by crowding and extinction. For this reason, compact object binaries orbiting around a massive black hole can be used as probes of the inner environment of the black hole in regions inaccessible to traditional astronomical measurements. The orbit of the binary's barycenter around the massive black hole will cause a Doppler shift in the gravitational waveform which is in principle measurable by future space-based gravitational-wave interferometers, such as the Laser Interferometer Space Antenna (LISA). We investigate the conditions under which these Doppler shifts are observable by LISA. Our results imply that Doppler shift observations can be used to study the central region of globular clusters in the Milky Way, as well the central environment of extragalactic massive black holes.

Pair instabilities in supernovae might prevent the formation of black holes with masses between $\sim 50 M_\odot$ and $\sim 130 M_\odot$. Multiple generations of black-hole mergers provide a possible way to populate this "mass gap" from below. However this requires an astrophysical environment with a sufficiently large escape speed to retain merger remnants, and prevent them from being ejected by gravitational-wave recoils. We show that, if the mass gap is indeed populated by multiple mergers, the observation of a single black-hole binary component in the mass gap implies that its progenitors grew in an environment with escape speed $v_{\rm esc} \gtrsim 50$ km/s. This is larger than the escape speeds of most globular clusters, requiring denser and heavier environments such as nuclear star clusters or disks-assisted migration in galactic nuclei. A single detection in the upper mass gap would hint at the existence of a much larger population of first-generation events from the same environment, thus providing a tool to disentangle the contribution of different formation channels to the observed merger rate.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA.

In this paper, we present the analysis and results of a direct measurement of the cosmic-ray proton spectrum with the CALET instrument onboard the International Space Station, including the detailed assessment of systematic uncertainties. The observation period used in this analysis is from October 13, 2015 to August 31, 2018 (1054 days). We have achieved the very wide energy range necessary to carry out measurements of the spectrum from 50 GeV to 10 TeV covering, for the first time in space, with a single instrument the whole energy interval previously investigated in most cases in separate subranges by magnetic spectrometers (BESS-TeV, PAMELA, and AMS-02) and calorimetric instruments (ATIC, CREAM, and NUCLEON). The observed spectrum is consistent with AMS-02 but extends to nearly an order of magnitude higher energy, showing a very smooth transition of the power-law spectral index from -2.81 +- 0.03 (50--500 GeV) neglecting solar modulation effects (or -2.87 +- 0.06 including solar modulation effects in the lower energy region) to -2.56 +- 0.04 (1--10 TeV), thereby confirming the existence of spectral hardening and providing evidence of a deviation from a single power law by more than 3 sigma.

Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected number of multiband detections first by using information from current LIGO/Virgo data, and then through population synthesis simulations of binary stars. We estimate that few to tens of LISA detections can be used to predict mergers detectable on the ground. Conversely, hundreds of events could potentially be extracted from the LISA data stream using prior information from ground detections. In general, the merger signal of binaries observable by LISA is strong enough to be unambiguously identified by both current and future ground-based detectors. Therefore third-generation detectors will not increase the number of multiband detections compared to LIGO/Virgo. We use population synthesis simulations of isolated binary stars to explore some of the stellar physics that could be constrained with multiband events, and we show that specific formation pathways might be overrepresented in multiband events compared to ground-only detections.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

Within the framework of 2PN black-hole binary spin precession, we explore configurations where one of the two spins oscillates from being completely aligned with the orbital angular momentum to being completely anti-aligned with it during a single precession cycle. This "wide nutation" is the extreme limit of the generic phenomenon of spin nutation in black-hole binaries. Crucially, wide nutation happens on the short precession time scale and it is not a secular effect due to gravitational-wave radiation reaction. The spins of these binaries, therefore, flip repeatedly as one of these special configurations is entered. Binaries with total mass $M$, mass ratio $q$, and dimensionless spin $\chi_1$ ($\chi_2$) of the more (less) massive black hole are allowed to undergo wide nutation at binary separations $r \leq r_{\rm wide} \equiv [(q \chi_2 - \chi_1)/(1-q)]^2 M$. Sources that are more likely to nutate widely have similar masses and effective spins close to zero.

A space-based interferometer operating in the previously unexplored mHz gravitational band has tremendous discovery potential. If history is any guide, the opening of a new spectral band will lead to the discovery of entirely new sources and phenomena. The mHz band is ideally suited to exploring beyond standard model processes in the early universe, and with the sensitivities that can be reached with current technologies, the discovery space for exotic astrophysical systems is vast.

Scalar-tensor theories of gravity where a new scalar degree of freedom couples to the Gauss-Bonnet invariant can exhibit the phenomenon of spontaneous black hole scalarization. These theories admit both the classic black holes predicted by general relativity as well as novel hairy black hole solutions. The stability of the hairy black holes is strongly dependent on the precise form of the scalar-gravity coupling. A radial stability investigation revealed that all scalarized black hole solutions are unstable when the coupling between the scalar field and the Gauss-Bonnet invariant is quadratic in the scalar, whereas stable solutions exist for exponential couplings. Here we elucidate this behavior. We demonstrate that, while the quadratic term controls the onset of the tachyonic instability that gives rise to the black hole hair, the higher-order coupling terms control the nonlinearities that quench that instability, and hence also control the stability of the hairy black hole solutions.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

The CALorimetric Electron Telescope (CALET) is a high-energy astroparticle physics space experiment installed on the International Space Station (ISS), developed and operated by Japan in collaboration with Italy and the United States. The CALET mission goals include the investigation of possible nearby sources of high-energy electrons, of the details of galactic particle acceleration and propagation, and of potential signatures of dark matter. CALET measures the cosmic-ray electron + positron flux up to 20 TeV, gamma-rays up to 10 TeV, and nuclei with Z=1 to 40 up to 1,000 TeV for the more abundant elements during a long-term observation aboard the ISS. Starting science operation in mid-October 2015, CALET performed continuous observation without major interruption with close to 20 million triggered events over 10 GeV per month. Based on the data taken during the first two-years, we present an overview of CALET observations: uses w/o major interruption 1) Electron + positron energy spectrum, 2) Nuclei analysis, 3) Gamma-ray observation including a characterization of on-orbit performance. Results of the electromagnetic counterpart search for LIGO/Virgo gravitational wave events are discussed as well.

It has recently been shown that nontrivial couplings between a scalar and the Gauss-Bonnet invariant can give rise to black hole spontaneous scalarization. Theories that exhibit this phenomenon are among the leading candidates for testing gravity with upcoming black hole observations. All models considered so far have focused on specific forms for the coupling, neglecting scalar self-interactions. In this work, we take the first steps towards placing this phenomenon on a more robust theoretical footing by considering the leading-order scalar self-interactions as well as the scalar-Gauss-Bonnet coupling. Our approach is consistent with the principles of effective field theory and yields the simplest and most natural model. We find that a mass term for the scalar alters the threshold for the onset of scalarization, and we study the mass range over which scalarized black hole solutions exist. We also demonstrate that the quartic self-coupling is sufficient to produce scalarized solutions that are stable against radial perturbations, without the need to resort to higher-order terms in the Gauss-Bonnet coupling function. Our model therefore represents a canonical model that can be studied further, with the ultimate aim of developing falsifiable tests of black hole scalarization.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

The LIGO/Virgo gravitational-wave (GW) interferometers have to-date detected ten merging black hole (BH) binaries, some with masses considerably larger than had been anticipated. Stellar-mass BH binaries at the high end of the observed mass range (with "chirp mass" ${\cal M} \gtrsim 25 M_{\odot}$) should be detectable by a space-based GW observatory years before those binaries become visible to ground-based GW detectors. This white paper discusses some of the synergies that result when the same binaries are observed by instruments in space and on the ground. We consider intermediate-mass black hole binaries (with total mass $M \sim 10^2 -10^4 M_{\odot}$) as well as stellar-mass black hole binaries. We illustrate how combining space-based and ground-based data sets can break degeneracies and thereby improve our understanding of the binary's physical parameters. While early work focused on how space-based observatories can forecast precisely when some mergers will be observed on the ground, the reverse is also important: ground-based detections will allow us to "dig deeper" into archived, space-based data to confidently identify black hole inspirals whose signal-to-noise ratios were originally sub-threshold, increasing the number of binaries observed in both bands by a factor of $\sim 4 - 7$.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

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.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

Cunha et al. (2018) recently reexamined the possibility of detecting gravitational waves from exoplanets, claiming that three ultra-short period systems would be observable by LISA. We revisit their analysis and conclude that the currently known exoplanetary systems are unlikely to be detectable, even assuming a LISA observation time $T_{\rm obs}=4$ yrs. Conclusive statements on the detectability of one of these systems, GP Com b, will require better knowledge of the system's properties, as well as more careful modeling of both LISA's response and the galactic confusion noise. Still, the possibility of exoplanet detection with LISA is interesting enough to warrant further study, as gravitational waves could yield dynamical properties that are difficult to constrain with electromagnetic observations.

We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $\chi_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $\chi_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars.

The first two LIGO/Virgo observing runs have detected several black hole binary mergers. One of the most exciting prospects of future observing runs is the possibility to identify the remnants of these mergers as Kerr black holes by measuring their (complex) quasinormal mode frequencies. This idea - similar to the identification of atomic elements through their spectral lines - is sometimes called "black hole spectroscopy". Third-generation Earth-based detectors and the space-based interferometer LISA could measure multiple spectral lines from different multipolar components of the radiation, and therefore provide qualitatively better tests of the Kerr hypothesis. In this paper we quantify the redshift out to which the various modes would be detectable (or, conversely, the number of detectable modes at any given redshift) as a function of the intrinsic parameters of the merging binary. LISA could detect so many modes that current numerical relativity simulations would not be sufficient to extract all available science from the data.

The Laser Interferometer Space Antenna (LISA) gravitational-wave (GW) observatory will be limited in its ability to detect mergers of binary black holes (BBHs) in the stellar-mass range. A future ground-based detector network, meanwhile, will achieve by the LISA launch date a sensitivity that ensures complete detection of all mergers within a volume $>\!\mathcal{O}(10)\,{\rm Gpc}^{3}$. We propose a method to use the information from the ground to revisit the LISA data in search for sub-threshold events. By discarding spurious triggers that do not overlap with the ground-based catalogue, we show that the signal-to-noise threshold $\rho_{\rm LISA}$ employed in LISA can be significantly lowered, greatly boosting the detection rate. The efficiency of this method depends predominantly on the rate of false-alarm increase when the threshold is lowered and on the uncertainty in the parameter estimation for the LISA events. As an example, we demonstrate that while all current LIGO BBH-merger detections would have evaded detection by LISA when employing a standard $\rho_{\rm LISA}=8$ threshold, this method will allow us to easily (possibly) detect an event similar to GW150914 (GW170814) in LISA. Overall, we estimate that the total rate of stellar-mass BBH mergers detected by LISA can be boosted by a factor $\sim\!4$ ($\gtrsim\!8$) under conservative (optimistic) assumptions. This will enable new tests using multi-band GW observations, significantly aided by the greatly increased lever arm in frequency.

We present results on searches for gamma-ray counterparts of the LIGO/Virgo gravitational-wave events using CALorimetric Electron Telescope ({\sl CALET}) observations. The main instrument of {\sl CALET}, CALorimeter (CAL), observes gamma-rays from $\sim1$ GeV up to 10 TeV with a field of view of nearly 2 sr. In addition, the {\sl CALET} gamma-ray burst monitor (CGBM) views $\sim$3 sr and $\sim2\pi$ sr of the sky in the 7 keV -- 1 MeV and the 40 keV -- 20 MeV bands, respectively, by using two different crystal scintillators. The {\sl CALET} observations on the International Space Station started in October 2015, and here we report analyses of events associated with the following gravitational wave events: GW151226, GW170104, GW170608, GW170814 and GW170817. Although only upper limits on gamma-ray emission are obtained, they correspond to a luminosity of $10^{49}\sim10^{53}$ erg s$^{-1}$ in the GeV energy band depending on the distance and the assumed time duration of each event, which is approximately the order of luminosity of typical short gamma-ray bursts. This implies there will be a favorable opportunity to detect high-energy gamma-ray emission in further observations if additional gravitational wave events with favorable geometry will occur within our field-of-view. We also show the sensitivity of {\sl CALET} for gamma-ray transient events which is the order of $10^{-7}$~erg\,cm$^{-2}$\,s$^{-1}$ for an observation of 100~s duration.

Extended results on the cosmic-ray electron + positron spectrum from 11 GeV to 4.8 TeV are presented based on observations with the Calorimetric Electron Telescope (CALET) on the International Space Station utilizing the data up to November 2017. The analysis uses the full detector acceptance at high energies, approximately doubling the statistics compared to the previous result. CALET is an all-calorimetric instrument with a total thickness of 30 $X_0$ at normal incidence and fine imaging capability, designed to achieve large proton rejection and excellent energy resolution well into the TeV energy region. The observed energy spectrum in the region below 1 TeV shows good agreement with Alpha Magnetic Spectrometer (AMS-02) data. In the energy region below $\sim$300 GeV, CALET's spectral index is found to be consistent with the AMS-02, Fermi Large Area Telescope (Fermi-LAT) and Dark Matter Particle Explorer (DAMPE), while from 300 to 600 GeV the spectrum is significantly softer than the spectra from the latter two experiments. The absolute flux of CALET is consistent with other experiments at around a few tens of GeV. However, it is lower than those of DAMPE and Fermi-LAT with the difference increasing up to several hundred GeV. The observed energy spectrum above $\sim$1 TeV suggests a flux suppression consistent within the errors with the results of DAMPE, while CALET does not observe any significant evidence for a narrow spectral feature in the energy region around 1.4 TeV. Our measured all-electron flux, including statistical errors and a detailed breakdown of the systematic errors, is tabulated in the Supplemental Material in order to allow more refined spectral analyses based on our data.

We compare binary evolution models with different assumptions about black-hole natal kicks to the first gravitational-wave observations performed by the LIGO detectors. Our comparisons attempt to reconcile merger rate, masses, spins, and spin-orbit misalignments of all current observations with state-of-the-art formation scenarios of binary black holes formed in isolation. We estimate that black holes (BHs) should receive natal kicks at birth of the order of $\sigma\simeq 200$ (50) km/s if tidal processes do (not) realign stellar spins. Our estimate is driven by two simple factors. The natal kick dispersion $\sigma$ is bounded from above because large kicks disrupt too many binaries (reducing the merger rate below the observed value). Conversely, the natal kick distribution is bounded from below because modest kicks are needed to produce a range of spin-orbit misalignments. A distribution of misalignments increases our models' compatibility with LIGO's observations, if all BHs are likely to have natal spins. Unlike related work which adopts a concrete BH natal spin prescription, we explore a range of possible BH natal spin distributions. Within the context of our models, for all of the choices of $\sigma$ used here and within the context of one simple fiducial parameterized spin distribution, observations favor low BH natal spin.

We infer the mass distribution of neutron stars in binary systems using a flexible Gaussian mixture model and use Bayesian model selection to explore evidence for multi-modality and a sharp cut-off in the mass distribution. We find overwhelming evidence for a bimodal distribution, in agreement with previous literature, and report for the first time positive evidence for a sharp cut-off at a maximum neutron star mass. We measure the maximum mass to be $2.0M_\odot < m_\mathrm{max} < 2.2M_\odot$ (68\%), $2.0M_\odot < m_\mathrm{max}< 2.6M_\odot$ (90\%), and evidence for a cut-off is robust against the choice of model for the mass distribution and to removing the most extreme (highest mass) neutron stars from the dataset. If this sharp cut-off is interpreted as the maximum stable neutron star mass allowed by the equation of state of dense matter, our measurement puts constraints on the equation of state. For a set of realistic equations of state that support $>2M_\odot$ neutron stars, our inference of $m_\mathrm{max}$ is able to distinguish between models at odds ratios of up to $12:1$, whilst under a flexible piecewise polytropic equation of state model our maximum mass measurement improves constraints on the pressure at $3-7\times$ the nuclear saturation density by $\sim 30-50\%$ compared to simply requiring $m_\mathrm{max}> 2M_\odot$. We obtain a lower bound on the maximum sound speed attained inside the neutron star of $c_s^\mathrm{max} > 0.63c$ (99.8\%), ruling out $c_s^\mathrm{max} < c/\sqrt{3}$ at high significance. Our constraints on the maximum neutron star mass strengthen the case for neutron star-neutron star mergers as the primary source of short gamma-ray bursts.

The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better.

The LIGO/Virgo detections of binary black hole mergers marked a watershed moment in astronomy, ushering in the era of precision tests of Kerr dynamics. We review theoretical and experimental challenges that must be overcome to carry out black hole spectroscopy with present and future gravitational wave detectors. Among other topics, we discuss quasinormal mode excitation in binary mergers, astrophysical event rates, tests of black hole dynamics in modified theories of gravity, parameterized "post-Kerr" ringdown tests, exotic compact objects, and proposed data analysis methods to improve spectroscopic tests of Kerr dynamics by stacking multiple events.

The observation of the inspiral and merger of compact binaries by the LIGO/Virgo collaboration ushered in a new era in the study of strong-field gravity. We review current and future tests of strong gravity and of the Kerr paradigm with gravitational-wave interferometers, both within a theory-agnostic framework (the parametrized post-Einsteinian formalism) and in the context of specific modified theories of gravity (scalar-tensor, Einstein-dilaton-Gauss-Bonnet, dynamical Chern-Simons, Lorentz-violating, and extra dimensional theories). In this contribution we focus on (i) the information carried by the inspiral radiation, and (ii) recent progress in numerical simulations of compact binary mergers in modified gravity.

We identify a class of scalar-tensor theories with coupling between the scalar and the Gauss-Bonnet invariant that exhibit spontaneous scalarization for both black holes and compact stars. In particular, these theories formally admit all of the stationary solutions of general relativity, but these are not dynamically preferred if certain conditions are satisfied. Remarkably, black holes exhibit scalarization if their mass lies within one of many narrow bands. We find evidence that scalarization can occur in neutron stars as well.

We prove the following theorem: axisymmetric, stationary solutions of the Einstein field equations formed from classical gravitational collapse of matter obeying the null energy condition, that are everywhere smooth and ultracompact (i.e., they have a light ring) must have at least two light rings, and one of them is stable. It has been argued that stable light rings generally lead to nonlinear spacetime instabilities. Our result implies that smooth, physically and dynamically reasonable ultracompact objects are not viable as observational alternatives to black holes whenever these instabilities occur on astrophysically short time scales. The proof of the theorem has two parts: (i) We show that light rings always come in pairs, one being a saddle point and the other a local extremum of an effective potential. This result follows from a topological argument based on the Brouwer degree of a continuous map, with no assumptions on the spacetime dynamics, and hence it is applicable to any metric gravity theory where photons follow null geodesics. (ii) Assuming Einstein's equations, we show that the extremum is a local minimum of the potential (i.e., a stable light ring) if the energy-momentum tensor satisfies the null energy condition.

The relaxation of a distorted black hole to its final state provides important tests of general relativity within the reach of current and upcoming gravitational wave facilities. In black hole perturbation theory, this phase consists of a simple linear superposition of exponentially damped sinusoids (the quasinormal modes) and of a power-law tail. How many quasinormal modes are necessary to describe waveforms with a prescribed precision? What error do we incur by only including quasinormal modes, and not tails? What other systematic effects are present in current state-of-the-art numerical waveforms? These issues, which are basic to testing fundamental physics with distorted black holes, have hardly been addressed in the literature. We use numerical relativity waveforms and accurate evolutions within black hole perturbation theory to provide some answers. We show that (i) a determination of the fundamental $l=m=2$ quasinormal mode to within $1\%$ or better requires the inclusion of at least the first overtone, and preferably of the first two or three overtones; (ii) a determination of the black hole mass and spin with precision better than $1\%$ requires the inclusion of at least two quasinormal modes for any given angular harmonic mode $(\ell,\,m)$. We also improve on previous estimates and fits for the ringdown energy radiated in the various multipoles. These results are important to quantify theoretical (as opposed to instrumental) limits in parameter estimation accuracy and tests of general relativity allowed by ringdown measurements with high signal-to-noise ratio gravitational wave detectors.

Ultralight scalar fields around spinning black holes can trigger superradiant instabilities, forming a long-lived bosonic condensate outside the horizon. We use numerical solutions of the perturbed field equations and astrophysical models of massive and stellar-mass black hole populations to compute, for the first time, the stochastic gravitational-wave background from these sources. In optimistic scenarios the background is observable by Advanced LIGO and LISA for field masses $m_s$ in the range $\sim [2\times 10^{-13}, 10^{-12}]\,{\rm eV}$ and $\sim 5\times[ 10^{-19}, 10^{-16}]\,{\rm eV}$, respectively, and it can affect the detectability of resolvable sources. Our estimates suggest that an analysis of the stochastic background limits from LIGO O1 might already be used to marginally exclude axions with mass $\sim 10^{-12.5}{\rm eV}$. Semicoherent searches with Advanced LIGO (LISA) should detect $\sim 15~(5)$ to $200~(40)$ resolvable sources for scalar field masses $3\times 10^{-13}$ ($10^{-17}$) eV. LISA measurements of massive BH spins could either rule out bosons in the range $\sim [10^{-18}, 2\times 10^{-13}]$ eV, or measure $m_s$ with ten percent accuracy in the range $\sim[10^{-17}, 10^{-13}]$ eV.

Ultralight bosons can induce superradiant instabilities in spinning black holes, tapping their rotational energy to trigger the growth of a bosonic condensate. Possible observational imprints of these boson clouds include (i) direct detection of the nearly monochromatic (resolvable or stochastic) gravitational waves emitted by the condensate, and (ii) statistically significant evidence for the formation of "holes" at large spins in the spin versus mass plane (sometimes also referred to as "Regge plane") of astrophysical black holes. In this work, we focus on the prospects of LISA and LIGO detecting or constraining scalars with mass in the range $m_s\in [10^{-19},\,10^{-15}]$ eV and $m_s\in [10^{-14},\,10^{-11}]$ eV, respectively. Using astrophysical models of black-hole populations calibrated to observations and black-hole perturbation theory calculations of the gravitational emission, we find that, in optimistic scenarios, LIGO could observe a stochastic background of gravitational radiation in the range $m_s\in [2\times 10^{-13}, 10^{-12}]$ eV, and up to $10^4$ resolvable events in a $4$-year search if $m_s\sim 3\times 10^{-13}\,{\rm eV}$. LISA could observe a stochastic background for boson masses in the range $m_s\in [5\times 10^{-19}, 5\times 10^{-16}]$, and up to $\sim 10^3$ resolvable events in a $4$-year search if $m_s\sim 10^{-17}\,{\rm eV}$. LISA could further measure spins for black-hole binaries with component masses in the range $[10^3, 10^7]~M_\odot$, which is not probed by traditional spin-measurement techniques. A statistical analysis of the spin distribution of these binaries could either rule out scalar fields in the mass range $\sim [4 \times 10^{-18}, 10^{-14}]$ eV, or measure $m_s$ with ten percent accuracy if light scalars in the mass range $\sim [10^{-17}, 10^{-13}]$ eV exist.

Advanced LIGO detectors at Hanford and Livingston made two confirmed and one marginal detection of binary black holes during their first observing run. The first event, GW150914, was from the merger of two black holes much heavier that those whose masses have been estimated so far, indicating a formation scenario that might differ from "ordinary" stellar evolution. One possibility is that these heavy black holes resulted from a previous merger. When the progenitors of a black hole binary merger result from previous mergers, they should (on average) merge later, be more massive, and have spin magnitudes clustered around a dimensionless spin ~0.7. Here we ask the following question: can gravitational-wave observations determine whether merging black holes were born from the collapse of massive stars ("first generation"), rather than being the end product of earlier mergers ("second generation")? We construct simple, observationally motivated populations of black hole binaries, and we use Bayesian model selection to show that measurements of the masses, luminosity distance (or redshift), and "effective spin" of black hole binaries can indeed distinguish between these different formation scenarios.

One of the central goals of the newborn field of gravitational wave astronomy is to test gravity in the highly nonlinear, strong field regime characterizing the spacetime of black holes. In particular, "black hole spectroscopy" (the observation and identification of black hole quasinormal mode frequencies in the gravitational wave signal) is expected to become one of the main tools for probing the structure and dynamics of Kerr black holes. In this paper we take a significant step towards that goal by constructing a "post-Kerr" quasinormal mode formalism. The formalism incorporates a parametrized but general perturbative deviation from the Kerr metric and exploits the well-established connection between the properties of the spacetime's circular null geodesics and the fundamental quasinormal mode to provide approximate, eikonal limit formulae for the modes' complex frequencies. The resulting algebraic toolkit can be used in waveform templates for ringing black holes with the purpose of measuring deviations from the Kerr metric. As a first illustrative application of our framework, we consider the Johannsen-Psaltis deformed Kerr metric and compute the resulting deviation in the quasinormal mode frequency relative to the known Kerr result.

One of the key astrophysical sources for the Laser Interferometer Space Antenna (LISA) are the inspirals of stellar-origin compact objects into massive black holes in the centres of galaxies. These extreme-mass-ratio inspirals (EMRIs) have great potential for astrophysics, cosmology and fundamental physics. In this paper we describe the likely numbers and properties of EMRI events that LISA will observe. We present the first results computed for the 2.5 Gm interferometer that was the new baseline mission submitted in January 2017 in response to the ESA L3 mission call. In addition, we attempt to quantify the astrophysical uncertainties in EMRI event rate estimates by considering a range of different models for the astrophysical population. We present both likely event rates and estimates for the precision with which the parameters of the observed sources could be measured. We finish by discussing the implications of these results for science using EMRIs.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

We study the contribution of binary black hole (BH-BH) mergers from the first, metal-free stars in the Universe (Pop III) to gravitational wave detection rates. Our study combines initial conditions for the formation of Pop III stars based on N-body simulations of binary formation (including rates, binary fraction, initial mass function, orbital separation and eccentricity distributions) with an updated model of stellar evolution specific for Pop III stars. We find that the merger rate of these Pop III BH-BH systems is relatively small (< 0.1 Gpc^-3 yr^-1) at low redshifts (z<2), where it can be compared with the LIGO empirical estimate of 9-240 Gpc^-3 yr^-1 (Abbott et al. 2016). The predicted rates are even smaller for Pop III double neutron star and black hole neutron star mergers. Our rates are compatible with those of Hartwig et al. (2016), but significantly smaller than those found in previous work (Bond & Carr 1984; Belczynski et al. 2004; Kinugawa et al. 2014, 2016). We explain the reasons for this discrepancy by means of detailed model comparisons and point out that (i) identification of Pop III BH-BH mergers may not be possible by advanced LIGO, and (ii) the level of stochastic gravitational wave background from Pop III mergers may be lower than recently estimated (Kowalska et al. 2012; Inayoshi et al. 2016; Dvorkin et al. 2016). We further estimate gravitational wave detection rates for third-generation interferometric detectors. Our calculations are relevant for low to moderately rotating Pop III stars. We can now exclude significant (> 1 per cent) contribution of these stars to low-redshift BH-BH mergers. However, it remains to be tested whether (and at what level) rapidly spinning Pop III stars (homogeneous evolution) can contribute to BH-BH mergers in the local Universe.

We estimate the potential of present and future interferometric gravitational-wave detectors to test the Kerr nature of black holes through "gravitational spectroscopy," i.e. the measurement of multiple quasinormal mode frequencies from the remnant of a black hole merger. Using population synthesis models of the formation and evolution of stellar-mass black hole binaries, we find that Voyager-class interferometers will be necessary to perform these tests. Gravitational spectroscopy in the local Universe may become routine with the Einstein Telescope, but a 40-km facility like Cosmic Explorer is necessary to go beyond $z\sim 3$. In contrast, eLISA-like detectors should carry out a few - or even hundreds - of these tests every year, depending on uncertainties in massive black hole formation models. Many space-based spectroscopical measurements will occur at high redshift, testing the strong gravity dynamics of Kerr black holes in domains where cosmological corrections to general relativity (if they occur in nature) must be significant.

Up to hundreds of black hole binaries individually resolvable by eLISA will coalesce in the Advanced LIGO/Virgo band within ten years, allowing for multi-band gravitational wave observations. Binaries formed via dynamical interactions in dense star clusters are expected to have eccentricities $e_0\sim 10^{-3}$-$10^{-1}$ at the frequencies $f_0=10^{-2}$ Hz where eLISA is most sensitive, while binaries formed in the field should have negligible eccentricity in both frequency bands. We estimate that eLISA should always be able to detect a nonzero $e_0$ whenever $e_0\gtrsim 10^{-2}$; if $e_0\sim 10^{-3}$, eLISA should detect nonzero eccentricity for a fraction $\sim 90\%$ ($\sim 25\%$) of binaries when the observation time is $T_{\rm obs}=5$ ($2$) years, respectively. Therefore eLISA observations of BH binaries have the potential to distinguish between field and cluster formation scenarios.

Mergers of two stellar origin black holes are a prime source of gravitational waves and are under intensive investigations. One crucial ingredient in their modeling has so far been neglected. Pair-instability pulsation supernovae with associated severe mass loss may suppress formation of massive black holes, decreasing black hole merger rates for the highest black hole masses. The mass loss associated with pair-instability pulsation supernovae limits the Population I/II stellar-origin black hole mass to 50 Msun, in tension with earlier predictions that the maximum black hole mass could be as high as 100 Msun. Suppression of double black hole merger rates by pair-instability pulsation supernovae is negligible for our evolutionary channel. Our standard evolutionary model with inclusion of pair-instability pulsation supernovae and with pair-instability supernovae is fully consistent with the LIGO detections of black hole mergers: GW150914, GW151226 and LVT151012. The LIGO observations seem to exclude high (>400 km/s) BH natal kicks. We predict the detection of several up to 60 BH-BH mergers with total mass 10--150 Msun (most likely range: 20--80 Msun) in the forthcoming 60 effective days of LIGO O2 observations.

A space-based interferometer such as eLISA could observe few to few thousands progenitors of black hole binaries (BHBs) similar to those recently detected by Advanced LIGO. Gravitational radiation circularizes the orbit during inspiral, but some BHBs retain a measurable eccentricity at the low frequencies where eLISA is most sensitive. The eccentricity of a BHB carries precious information about its formation channel: BHBs formed in the field, in globular clusters, or close to a massive black hole (MBH) have distinct eccentricity distributions in the eLISA band. We generate mock eLISA observations, folding in measurement errors, and using Bayesian model selection we study whether eLISA measurements can identify the BHB formation channel. We find that a handful of observations would suffice to tell whether BHBs were formed in the gravitational field of a MBH. Conversely, several tens of observations are needed to tell apart field formation from globular cluster formation. A five-year eLISA mission with the longest possible armlength is desirable to shed light on BHB formation scenarios.

Horndeski gravity holds a special position as the most general extension of Einstein's theory of general relativity with a single scalar degree of freedom and second-order field equations. Because of these features, Horndeski gravity is an attractive phenomenological playground to investigate the consequences of modifications of general relativity in cosmology and astrophysics. We present a review of the progress made so far in the study of compact objects (black holes and neutron stars) within Horndeski gravity. In particular, we review our recent work on slowly rotating black holes and present some new results on slowly rotating neutron stars.

Horndeski's theory of gravity is the most general scalar-tensor theory with a single scalar whose equations of motion contain at most second-order derivatives. A subsector of Horndeski's theory known as "Fab Four" gravity allows for dynamical self-tuning of the quantum vacuum energy, and therefore it has received particular attention in cosmology as a possible alternative to the $\Lambda$CDM model. Here we study compact stars in Fab Four gravity, which includes as special cases general relativity ("George"), Einstein-dilaton-Gauss-Bonnet gravity ("Ringo"), theories with a nonminimal coupling with the Einstein tensor ("John") and theories involving the double-dual of the Riemann tensor ("Paul"). We generalize and extend previous results in theories of the John class and were not able to find realistic compact stars in theories involving the Paul class.

The bulk properties of spherically symmetric stars in general relativity can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations. In previous work we developed a "post-TOV" formalism - inspired by parametrized post-Newtonian theory - which allows us to classify in a parametrized, phenomenological form all possible perturbative deviations from the structure of compact stars in general relativity that may be induced by modified gravity at second post-Newtonian order. In this paper we extend the formalism to deal with the stellar exterior, and we compute several potential astrophysical observables within the post-TOV formalism: the surface redshift $z_s$, the apparent radius $R_{\rm app}$, the Eddington luminosity at infinity $L_{\rm E}^\infty$ and the orbital frequencies. We show that, at leading order, all of these quantities depend on just two post-TOV parameters $\mu_1$ and $\chi$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

The lowest neutron star masses currently measured are in the range $1.0-1.1~M_\odot$, but these measurement have either large uncertainties or refer to isolated neutron stars. The recent claim of a precisely measured mass $M/M_{\odot} = 1.174 \pm 0.004$ by Martinez et al [Astrophys. J. 812, 143 (2015)] in a double neutron star system suggests that low-mass neutron stars may be an interesting target for gravitational-wave detectors. Furthermore, Sotani et al [PTEP 2014, 051E01 (2014)] recently found empirical formulas relating the mass and surface redshift of nonrotating neutron stars to the star's central density and to the parameter $\eta\equiv (K_0 L^2)^{1/3}$, where $K_0$ is the incompressibility of symmetric nuclear matter and $L$ is the slope of the symmetry energy at saturation density. Motivated by these considerations, we extend the work by Sotani et al to slowly rotating and tidally deformed neutron stars. We compute the moment of inertia, quadrupole moment, quadrupole ellipticity, tidal and rotational Love number and apsidal constant of slowly rotating neutron stars by integrating the Hartle-Thorne equations at second order in rotation, and we fit all of these quantities as functions of $\eta$ and of the central density. These fits may be used to constrain $\eta$, either via observations of binary pulsars in the electromagnetic spectrum, or via near-future observations of inspiralling compact binaries in the gravitational-wave spectrum.

We compare evolutionary predictions of double compact object merger rate densities with initial and forthcoming LIGO/Virgo upper limits. We find that: (i) Due to the cosmological reach of advanced detectors, current conversion methods of population synthesis predictions into merger rate densities are insufficient. (ii) Our optimistic models are a factor of 18 below the initial LIGO/Virgo upper limits for BH-BH systems, indicating that a modest increase in observational sensitivity (by a factor of 2.5) may bring the first detections or first gravitational wave constraints on binary evolution. (iii) Stellar-origin massive BH-BH mergers should dominate event rates in advanced LIGO/Virgo and can be detected out to redshift z=2 with templates including inspiral, merger, and ringdown. Normal stars (<150 Msun) can produce such mergers with total redshifted mass up to 400 Msun. (iv) High black hole natal kicks can severely limit the formation of massive BH-BH systems (both in isolated binary and in dynamical dense cluster evolution), and thus would eliminate detection of these systems even at full advanced LIGO/Virgo sensitivity. We find that low and high black hole natal kicks are allowed by current observational electromagnetic constraints. (v) The majority of our models yield detections of all types of mergers (NS-NS, BH-NS, BH-BH) with advanced detectors. Numerous massive BH-BH merger detections will indicate small (if any) natal kicks for massive BHs.

We review recent progress in the application of numerical relativity techniques to astrophysics and high-energy physics. We focus on some developments that took place within the "Numerical Relativity and High Energy Physics" network, a Marie Curie IRSES action that we coordinated, namely: spin evolution in black hole binaries, high-energy black hole collisions, compact object solutions in scalar-tensor gravity, superradiant instabilities and hairy black hole solutions in Einstein's gravity coupled to fundamental fields, and the possibility to gain insight into these phenomena using analog gravity models.

Gravitational waves emitted by the merger of two black holes have been detected, setting the course for a new era of observational astrophysics.

Gravitational waves from coalescing binary black holes encode the evolution of their spins prior to merger. In the post-Newtonian regime and on the precession timescale, this evolution has one of three morphologies, with the spins either librating around one of two fixed points ("resonances") or circulating freely. In this work we perform full parameter estimation on resonant binaries with fixed masses and spin magnitudes, changing three parameters: a conserved "projected effective spin" $\xi$ and resonant family $\Delta\Phi=0,\pi$ (which uniquely label the source), the inclination $\theta_{JN}$ of the binary's total angular momentum with respect to the line of sight (which determines the strength of precessional effects in the waveform), and the signal amplitude. We demonstrate that resonances can be distinguished for a wide range of binaries, except for highly symmetric configurations where precessional effects are suppressed. Motivated by new insight into double-spin evolution, we introduce new variables to characterize precessing black hole binaries which naturally reflects the timescale separation of the system and therefore better encode the dynamical information carried by gravitational waves.

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers.

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.

We study black hole solutions at first order in the Hartle-Thorne slow-rotation approximation in Horndeski gravity theories. We derive the equations of motion including also cases where the scalar depends linearly on time. In the Hartle-Thorne formalism, all first-order rotational corrections are described by a single frame-dragging function. We show that the frame-dragging function is exactly the same as in general relativity for all known black hole solutions in shift symmetric Horndeski theories, with the exception of theories with a linear coupling to the Gauss-Bonnet invariant. Our results extend previous no-hair theorems for a broad class of Horndeski gravity theories.

The GAMMA-400 gamma-ray telescope with excellent angular and energy resolutions is designed to search for signatures of dark matter in the fluxes of gamma-ray emission and electrons + positrons. Precision investigations of gamma-ray emission from Galactic Center, Crab, Vela, Cygnus, Geminga, and other regions will be performed, as well as diffuse gamma-ray emission, along with measurements of high-energy electron + positron and nuclei fluxes. Furthermore, it will study gamma-ray bursts and gamma-ray emission from the Sun during periods of solar activity. The energy range of GAMMA-400 is expected to be from ~20 MeV up to TeV energies for gamma rays, up to 20 TeV for electrons + positrons, and up to 10E15 eV for cosmic-ray nuclei. For high-energy gamma rays with energy from 10 to 100 GeV, the GAMMA-400 angular resolution improves from 0.1{\deg} to ~0.01{\deg} and energy resolution from 3% to ~1%; the proton rejection factor is ~5x10E5. GAMMA-400 will be installed onboard the Russian space observatory.

Binary black holes on quasicircular orbits with spins aligned with their orbital angular momentum have been testbeds for analytic and numerical relativity for decades, not least because symmetry ensures that such configurations are equilibrium solutions to the spin-precession equations. In this work, we show that these solutions can be unstable when the spin of the higher-mass black hole is aligned with the orbital angular momentum and the spin of the lower-mass black hole is anti-aligned. Spins in these configurations are unstable to precession to large misalignment when the binary separation $r$ is between the values $r_{\rm ud\pm}= (\sqrt{\chi_1} \pm \sqrt{q \chi_2})^4 (1-q)^{-2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is the mass ratio, and $\chi_1$ ($\chi_2$) is the dimensionless spin of the more (less) massive black hole. This instability exists for a wide range of spin magnitudes and mass ratios and can occur in the strong-field regime near merger. We describe the origin and nature of the instability using recently developed analytical techniques to characterize fully generic spin precession. This instability provides a channel to circumvent astrophysical spin alignment at large binary separations, allowing significant spin precession prior to merger affecting both gravitational-wave and electromagnetic signatures of stellar-mass and supermassive binary black holes.

Tidal disruption has a dramatic impact on the outcome of neutron star-black hole mergers. The phenomenology of these systems can be divided in three classes: nondisruptive, mildly disruptive or disruptive. The cutoff frequency of the gravitational radiation produced during the merger (which is potentially measurable by interferometric detectors) is very different in each regime, and when the merger is disuptive it carries information on the neutron star equation of state. Here we use semianalytical tools to derive a formula for the critical binary mass ratio $Q=M_{\rm BH}/M_{\rm NS}$ below which mergers are disruptive as a function of the stellar compactness $\mathcal{C}=M_{\rm NS}/R_{\rm NS}$ and the dimensionless black hole spin $\chi$. We then employ a new gravitational waveform amplitude model, calibrated to $134$ general relativistic numerical simulations of binaries with black hole spin (anti-)aligned with the orbital angular momentum, to obtain a fit to the gravitational-wave cutoff frequency in the disruptive regime as a function of $\mathcal{C}$, $Q$ and $\chi$. Our findings are important to build gravitational wave template banks, to determine whether neutron star-black hole mergers can emit electromagnetic radiation (thus helping multimessenger searches), and to improve event rate calculations for these systems.

The dynamics of precessing binary black holes (BBHs) in the post-Newtonian regime has a strong timescale hierarchy: the orbital timescale is very short compared to the spin-precession timescale which, in turn, is much shorter than the radiation-reaction timescale on which the orbit is shrinking due to gravitational-wave emission. We exploit this timescale hierarchy to develop a multi-scale analysis of BBH dynamics elaborating on the analysis of Kesden et al. (2015). We solve the spin-precession equations analytically on the precession time and then implement a quasi-adiabatic approach to evolve these solutions on the longer radiation-reaction time. This procedure leads to an innovative "precession-averaged" post-Newtonian approach to studying precessing BBHs. We use our new solutions to classify BBH spin precession into three distinct morphologies, then investigate phase transitions between these morphologies as BBHs inspiral. These precession-averaged post-Newtonian inspirals can be efficiently calculated from arbitrarily large separations, thus making progress towards bridging the gap between astrophysics and numerical relativity.

Gravitational theories with multiple scalar fields coupled to the metric and each other --- a natural extension of the well studied single-scalar-tensor theories --- are interesting phenomenological frameworks to describe deviations from general relativity in the strong-field regime. In these theories, the $N$-tuple of scalar fields takes values in a coordinate patch of an $N$-dimensional Riemannian target-space manifold whose properties are poorly constrained by weak-field observations. Here we introduce for simplicity a non-trivial model with two scalar fields and a maximally symmetric target-space manifold. Within this model we present a preliminary investigation of spontaneous scalarization for relativistic, perfect fluid stellar models in spherical symmetry. We find that the scalarization threshold is determined by the eigenvalues of a symmetric scalar-matter coupling matrix, and that the properties of strongly scalarized stellar configurations additionally depend on the target-space curvature radius. In preparation for numerical relativity simulations, we also write down the $3+1$ decomposition of the field equations for generic tensor-multi-scalar theories.

The gravitational radiation emitted during the merger of a black hole with a neutron star is rather similar to the radiation from the merger of two black holes when the neutron star is not tidally disrupted. When tidal disruption occurs, gravitational waveforms can be broadly classified in two groups, depending on the spatial extent of the disrupted material. Extending previous work by some of us, here we present a phenomenological model for the gravitational waveform amplitude in the frequency domain encompassing the three possible outcomes of the merger: no tidal disruption, "mild" and "strong" tidal disruption. The model is calibrated to 134 general-relativistic numerical simulations of binaries where the black hole spin is either aligned or antialigned with the orbital angular momentum. All simulations were produced using the SACRA code and piecewise polytropic neutron star equations of state. The present model can be used to determine when black-hole binary waveforms are sufficient for gravitational-wave detection, to extract information on the equation of state from future gravitational-wave observations, to obtain more accurate estimates of black hole-neutron star merger event rates, and to determine the conditions under which these systems are plausible candidates as central engines of gamma-ray bursts, macronovae and kilonovae.

Besides their astrophysical interest, compact stars also provide an arena for understanding the properties of theories of gravity that differ from Einstein's general relativity. Numerous studies have shown that different modified theories of gravity can modify the bulk properties (such as mass and radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. Motivated by this observation, in this paper we take an alternative route and construct a stellar structure formalism which, without adhering to any particular theory of gravity, describes in a simple parametrized form the departure from compact stars in general relativity. This "post-Tolman-Oppenheimer-Volkoff (TOV)" formalism for spherical static stars is inspired by the well-known parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable correction terms to the fully relativistic TOV equations. We show how neutron star properties are modified within our formalism, paying special attention to the effect of each correction term. We also show that the formalism is equivalent to general relativity with an "effective" (gravity-modified) equation of state.

The unprecedented range of second-generation gravitational-wave (GW) observatories calls for refining the predictions of potential sources and detection rates. The coalescence of double compact objects (DCOs)---i.e., neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS), and black hole-black hole (BH-BH) binary systems---is the most promising source of GWs for these detectors. We compute detection rates of coalescing DCOs in second-generation GW detectors using the latest models for their cosmological evolution, and implementing inspiral-merger-ringdown (IMR) gravitational waveform models in our signal-to-noise ratio calculations. We find that: (1) the inclusion of the merger/ringdown portion of the signal does not significantly affect rates for NS-NS and BH-NS systems, but it boosts rates by a factor $\sim 1.5$ for BH-BH systems; (2) in almost all of our models BH-BH systems yield by far the largest rates, followed by NS-NS and BH-NS systems, respectively, and (3) a majority of the detectable BH-BH systems were formed in the early Universe in low-metallicity environments. We make predictions for the distributions of detected binaries and discuss what the first GW detections will teach us about the astrophysics underlying binary formation and evolution.

Some models (such as the Skyrme model, a low-energy effective field theory for QCD) suggest that the high-density matter prevailing in neutron star interiors may be significantly anisotropic. Anisotropy is known to affect the bulk properties of nonrotating neutron stars in General Relativity. In this paper we study the effects of anisotropy on slowly rotating stars in General Relativity. We also consider one of the most popular extensions of Einstein's theory, namely scalar-tensor theories allowing for spontaneous scalarization (a phase transition similar to spontaneous magnetization in ferromagnetic materials). Anisotropy affects the moment of inertia of neutron stars (a quantity that could potentially be measured in binary pulsar systems) in both theories. We find that the effects of scalarization increase (decrease) when the tangential pressure is bigger (smaller) than the radial pressure, and we present a simple criterion to determine the onset of scalarization by linearizing the scalar-field equation. Our calculations suggest that binary pulsar observations may constrain the degree of anisotropy or even, more optimistically, provide evidence for anisotropy in neutron star cores.

Soon after the discovery of the Kerr metric, Penrose realized that superradiance can be exploited to extract energy from black holes. The original idea (involving the breakup of a single particle) yields only modest energy gains. A variant of the Penrose process consists of particle collisions in the ergoregion. The collisional Penrose process has been explored recently in the context of dark matter searches, with the conclusion that the ratio $\eta$ between the energy of post-collision particles detected at infinity and the energy of the colliding particles should be modest ($\eta \lesssim 1.5$). Schnittman has shown that these studies underestimated the maximum efficiency by about one order of magnitude (i.e., $\eta \lesssim 15$). In this work we show that particle collisions in the vicinity of rapidly rotating black holes can produce high-energy ejecta and result in high efficiencies under much more generic conditions. The astrophysical likelihood of these events deserves further scrutiny, but our study hints at the tantalizing possibility that the collisional Penrose process may power gamma rays and ultra-high-energy cosmic rays.

The GAMMA-400 gamma-ray telescope is intended to measure the fluxes of gamma rays and cosmic-ray electrons and positrons in the energy range from 100 MeV to several TeV. Such measurements concern with the following scientific goals: search for signatures of dark matter, investigation of gamma-ray point and extended sources, studies of the energy spectra of Galactic and extragalactic diffuse emission, studies of gamma-ray bursts and gamma-ray emission from the active Sun, as well as high-precision measurements of spectra of high-energy electrons and positrons, protons, and nuclei up to the knee. The main components of cosmic rays are protons and helium nuclei, whereas the part of lepton component in the total flux is ~10E-3 for high energies. In present paper, the capability of the GAMMA-400 gamma-ray telescope to distinguish electrons and positrons from protons in cosmic rays is investigated. The individual contribution to the proton rejection is studied for each detector system of the GAMMA-400 gamma-ray telescope. Using combined information from all detector systems allow us to provide the proton rejection from electrons with a factor of ~4x10E5 for vertical incident particles and ~3x10E5 for particles with initial inclination of 30 degrees. The calculations were performed for the electron energy range from 50 GeV to 1 TeV.

GAMMA-400 is a new space mission, designed as a dual experiment, capable to study both high energy gamma rays (from $\sim$100 MeV to few TeV) and cosmic rays (electrons up to 20 TeV and nuclei up to $\sim$10$^{15}$ eV). The full simulation framework of GAMMA-400 is based on the Geant4 toolkit. The details of the gamma-ray reconstruction pipeline in the pre-shower and calorimeter will be outlined. The performance of GAMMA-400 (PSF, effective area) have been obtained using this framework. The most updated results on them will be shown.

We derive an effective potential for binary black-hole (BBH) spin precession at second post-Newtonian order. This effective potential allows us to solve the orbit-averaged spin-precession equations analytically for arbitrary mass ratios and spins. These solutions are quasiperiodic functions of time: after a fixed period the BBH spins return to their initial relative orientations and jointly precess about the total angular momentum by a fixed angle. Using these solutions, we classify BBH spin precession into three distinct morphologies between which BBHs can transition during their inspiral. We also derive a precession-averaged evolution equation for the total angular momentum that can be integrated on the radiation-reaction time and identify a new class of spin-orbit resonances that can tilt the direction of the total angular momentum during the inspiral. Our new results will help efforts to model and interpret gravitational waves from generic BBH mergers and predict the distributions of final spins and gravitational recoils.

GAMMA-400 is a new space mission which will be installed on board the Russian space platform Navigator. It is scheduled to be launched at the beginning of the next decade. GAMMA-400 is designed to study simultaneously gamma rays (up to 3 TeV) and cosmic rays (electrons and positrons from 1 GeV to 20 TeV, nuclei up to 10$^{15}$-10$^{16}$ eV). Being a dual-purpose mission, GAMMA-400 will be able to address some of the most impelling science topics, such as search for signatures of dark matter, cosmic-rays origin and propagation, and the nature of transients. GAMMA-400 will try to solve the unanswered questions on these topics by high-precision measurements of the Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission and the spectra of cosmic-ray electrons + positrons and nuclei, thanks to excellent energy and angular resolutions.

The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the parameters of individual binaries (such as their masses, spins and luminosity distance), and a deep understanding of the underlying cosmic massive black hole parent population. This wealth of unprecedented information can lead to breakthroughs in many areas of physics, including astrophysics, cosmology and fundamental physics. We review the current status of the field, recent progress and future challenges.

The present design of the new space observatory GAMMA-400 is presented in this paper. The instrument has been designed for the optimal detection of gamma rays in a broad energy range (from ~100 MeV up to 3 TeV), with excellent angular and energy resolution. The observatory will also allow precise and high statistic studies of the electron component in the cosmic rays up to the multi TeV region, as well as protons and nuclei spectra up to the knee region. The GAMMA-400 observatory will allow to address a broad range of science topics, like search for signatures of dark matter, studies of Galactic and extragalactic gamma-ray sources, Galactic and extragalactic diffuse emission, gamma-ray bursts and charged cosmic rays acceleration and diffusion mechanism up to the knee.

The measurements of gamma-ray fluxes and cosmic-ray electrons and positrons in the energy range from 100 MeV to several TeV, which will be implemented by the specially designed GAMMA-400 gamma-ray telescope, concern with the following broad range of science topics. Searching for signatures of dark matter, surveying the celestial sphere in order to study gamma-ray point and extended sources, measuring the energy spectra of Galactic and extragalactic diffuse gamma-ray emission, studying gamma-ray bursts and gamma-ray emission from the Sun, as well as high precision measuring spectra of high-energy electrons and positrons, protons and nuclei up to the knee. To clarify these scientific problems with the new experimental data the GAMMA-400 gamma-ray telescope possesses unique physical characteristics comparing with previous and present experiments. For gamma-ray energies more than 100 GeV GAMMA-400 provides the energy resolution of ~1% and angular resolution better than 0.02 deg. The methods developed to reconstruct the direction of incident gamma photon are presented in this paper, as well as, the capability of the GAMMA-400 gamma-ray telescope to distinguish electrons and positrons from protons in cosmic rays is investigated.

These notes were prepared for a lecture on black holes delivered at the DPG Physics School "General Relativity @ 99" (Physikzentrum Bad Honnef, Germany, September 2014). The common thread of the lecture is the relation between geodesic stability and black-hole perturbations in the geometric optics limit. Chapter 1 establishes notation and discusses a common misconception on Michell's "Newtonian black holes". Chapters 2 and 3 deal with particle dynamics and wave dynamics in black-hole spacetimes, respectively. All calculations should be simple enough that they can be done with pen and paper. Chapter 4 builds on this introduction to discuss two exciting topics in current research: critical phenomena in black-hole mergers and the black-hole bomb instability.

We study torsional oscillations of neutron stars in the scalar-tensor theory of gravity using the relativistic Cowling approximation. We compute unperturbed neutron star models adopting realistic equations of state for the neutron star's core and crust. For scalar-tensor theories that allow for spontaneous scalarization, the crust thickness can be significantly smaller than in general relativity. We derive the perturbation equation describing torsional oscillations in scalar-tensor theory, and we solve the corresponding eigenvalue problem to find the oscillation frequencies. The fundamental mode (overtone) frequencies become smaller (larger) than in general relativity for scalarized stellar models. Torsional oscillation frequencies may yield information on the crust microphysics if microphysics effects are not degenerate with strong-gravity effects, such as those due to scalarization. To address this issue, we consider two different models for the equation of state of the crust and we look at the effects of electron screening. The effect of scalarization on torsional oscillation frequencies turns out to be smaller than uncertainties in the microphysics for all spontaneous scalarization models allowed by binary pulsar observations. Our study shows that the observation of quasi-periodic oscillations (QPOs) following giant flares can be used to constrain the microphysics of neutron star crusts, whether spontaneous scalarization occurs or not.

The angular dependence of the gravitational radiation emitted in compact binary mergers and gravitational collapse is usually separated using spin-weighted spherical harmonics ${}_sY_{\ell m}$ of spin weight $s$, that reduce to the ordinary spherical harmonics $Y_{\ell m}$ when $s=0$. Teukolsky first showed that the perturbations of the Kerr black hole that may be produced as a result of these events are separable in terms of a different set of angular functions: the spin-weighted spheroidal harmonics ${}_sS_{\ell m n}$, where $n$ denotes the "overtone index" of the corresponding Kerr quasinormal mode frequency $\omega_{\ell m n}$. In this paper we compute the complex-valued scalar products of the ${}_sS_{\ell m n}$'s with the ${}_sY_{\ell m}$'s ("spherical-spheroidal mixing coefficients") and with themselves ("spheroidal-spheroidal mixing coefficients") as functions of the dimensionless Kerr parameter $j$. Tables of these coefficients and analytical fits of their dependence on $j$ are available online for use in gravitational-wave source modeling and in other applications of black-hole perturbation theory.

We construct models of slowly rotating, perfect-fluid neutron stars by extending the classical Hartle-Thorne formalism to generic scalar-tensor theories of gravity. Working at second order in the dimensionless angular momentum, we compute the mass M, radius R, scalar charge q, moment of inertia I and spin-induced quadrupole moment Q, as well as the tidal and rotational Love numbers. Our formalism applies to generic scalar-tensor theories, but we focus in particular on theories that allow for spontaneous scalarization. It was recently discovered that the moment of inertia, quadrupole moment and Love numbers are connected by approximately universal (i.e., equation-of-state independent) "I-Love-Q" relations. We find that similar relations hold also for spontaneously scalarized stars. More interestingly, the I-Love-Q relations in scalar-tensor theories coincide with the general relativistic ones within less than a few percent, even for spontaneously scalarized stars with the largest couplings allowed by current binary-pulsar constraints. This implies that astrophysical measurements of these parameters cannot be used to discriminate between general relativity and scalar-tensor theories, even if spontaneous scalarization occurs in nature. Because of the well known equivalence between f(R) theories and scalar-tensor theories, the theoretical framework developed in this paper can be used to construct rotating compact stellar models in f(R) gravity. Our slow-rotation expansion can also be used as a benchmark for numerical calculations of rapidly spinning neutron stars in generic scalar-tensor theories.

If binary black holes form following the successive core collapses of sufficiently massive binary stars, precessional dynamics may align their spins $\mathbf S_1$ and $\mathbf S_2$ and the orbital angular momentum $\mathbf L$ into a plane in which they jointly precess about the total angular momentum $\mathbf J$. These spin orientations are known as spin-orbit resonances since $\mathbf S_1$, $\mathbf S_2$, and $\mathbf L$ all precess at the same frequency to maintain their planar configuration. Two families of such spin-orbit resonances exist, alike in dignity but differentiated by whether the components of the two spins in the orbital plane are either aligned or antialigned. The fraction of binary black holes in each family is determined by the stellar evolution of their progenitors, so if gravitational-wave detectors could measure this fraction they could provide important insights into astrophysical formation scenarios for binary black holes. In this paper, we show that even under the conservative assumption that binary black holes are observed along the direction of $\mathbf J$ (where precession-induced modulations to the gravitational waveforms are minimized), the waveforms of many members of each resonant family can be distinguished from all members of the opposing family in events with signal-to-noise ratios $\rho \simeq 10$, typical of those expected for the first detections with Advanced LIGO/Virgo. We hope that our preliminary findings inspire a greater appreciation of the capability of gravitational-wave detectors to constrain stellar astrophysics and stimulate further studies of the distinguishability of spin-orbit resonant families in more expanded regions of binary black-hole parameter space.

We lay the foundations for the construction of analytic expressions for Fourier-domain gravitational waveforms produced by eccentric, inspiraling compact binaries in a post-circular or small-eccentricity approximation. The time-dependent, "plus" and "cross" polarizations are expanded in Bessel functions, which are then self-consistently re-expanded in a power series about zero initial eccentricity to eighth order. The stationary phase approximation is then employed to obtain explicit analytic expressions for the Fourier transform of the post-circular expanded, time-domain signal. We exemplify this framework by considering Newtonian-accurate waveforms, which in the post-circular scheme give rise to higher harmonics of the orbital phase and amplitude corrections both to the amplitude and the phase of the Fourier domain waveform. Such higher harmonics lead to an effective increase in the inspiral mass reach of a detector as a function of the binary's eccentricity e_0 at the time when the binary enters the detector sensitivity band. Using the largest initial eccentricity allowed by our approximations (e_0 < 0.4), the mass reach is found to be enhanced up to factors of approximately 5 relative to that of circular binaries for Advanced LIGO, LISA, and the proposed Einstein Telescope at a signal-to-noise ratio of ten. A post-Newtonian generalization of the post circular scheme is also discussed, which holds the promise to provide "ready-to-use" Fourier-domain waveforms for data analysis of eccentric inspirals.

The development of advanced gravitational wave (GW) observatories, such as Advanced LIGO and Advanced Virgo, provides impetus to refine theoretical predictions for what these instruments might detect. In particular, with the range increasing by an order of magnitude, the search for GW sources is extending beyond the "local" Universe and out to cosmological distances. Double compact objects (neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS) and black hole-black hole (BH-BH) systems) are considered to be the most promising gravitational wave sources. In addition, NS-NS and/or BH-NS systems are thought to be the progenitors of gamma ray bursts (GRBs), and may also be associated with kilonovae. In this paper we present the merger event rates of these objects as a function of cosmological redshift. We provide the results for four cases, each one investigating a different important evolution parameter of binary stars. Each case is also presented for two metallicity evolution scenarios. We find that (i) in most cases NS-NS systems dominate the merger rates in the local Universe, while BH-BH mergers dominate at high redshift; (ii) BH-NS mergers are less frequent than other sources per unit volume, for all time; and (iii) natal kicks may alter the observable properties of populations in a significant way, allowing the underlying models of binary evolution and compact object formation to be easily distinguished. This is the second paper in a series of three. The third paper will focus on calculating the detection rates of mergers by gravitational wave telescopes.

Black hole-neutron star binary mergers display a much richer phenomenology than black hole-black hole mergers, even in the relatively simple case - considered in this paper - in which both the black hole and the neutron star are nonspinning. When the neutron star is tidally disrupted, the gravitational wave emission is radically different from the black hole-black hole case and it can be broadly classified in two groups, depending on the spatial extent of the disrupted material. We present a phenomenological model for the gravitational waveform amplitude in the frequency domain that encompasses the three possible outcomes of the merger: no tidal disruption, "mild", and "strong" tidal disruption. The model is calibrated to general relativistic numerical simulations using piecewise polytropic neutron star equations of state. It should prove useful to extract information on the nuclear equation of state from future gravitational-wave observations, and also to obtain more accurate estimates of black hole-neutron star merger event rates in second- and third-generation interferometric gravitational-wave detectors. We plan to extend and improve the model as longer and more accurate gravitational waveforms become available, and we will make it publicly available online as a Mathematica package. We also present in appendix analytical fits of the projected KAGRA noise spectral density, that should be useful in data analysis applications.

Using stellar structure calculations in the Hartle-Thorne approximation, we derive analytic expressions connecting the ellipticity of the stellar surface to the compactness, the spin angular momentum, and the quadrupole moment of the spacetime. We also obtain empirical relations between the compactness, the spin angular momentum, and the spacetime quadrupole. Our formulae reproduce the results of numerical calculations to within a few percent and help reduce the number of parameters necessary to model the observational appearance of moderately spinning neutron stars. This is sufficient for comparing theoretical spectroscopic and timing models to observations that aim to measure the masses and radii of neutron stars and to determine the equation of state prevailing in their interiors.

Scalar-tensor theories are a compelling alternative to general relativity and one of the most accepted extensions of Einstein's theory. Black holes in these theories have no hair, but could grow "wigs" supported by time-dependent boundary conditions or spatial gradients. Time-dependent or spatially varying fields lead in general to nontrivial black hole dynamics, with potentially interesting experimental consequences. We carry out a numerical investigation of the dynamics of single and binary black holes in the presence of scalar fields. In particular we study gravitational and scalar radiation from black-hole binaries in a constant scalar-field gradient, and we compare our numerical findings to analytical models. In the single black hole case we find that, after a short transient, the scalar field relaxes to static configurations, in agreement with perturbative calculations. Furthermore we predict analytically (and verify numerically) that accelerated black holes in a scalar-field gradient emit scalar radiation. For a quasicircular black-hole binary, our analytical and numerical calculations show that the dominant component of the scalar radiation is emitted at twice the binary's orbital frequency.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

The analytical understanding of quasinormal mode ringing requires an accurate knowledge of the Green's function describing the response of the black hole to external perturbations. We carry out a comprehensive study of quasinormal mode excitation for Kerr black holes. Relying on the formalism developed by Mano, Suzuki and Takasugi, we improve and extend previous calculations of the quasinormal mode residues in the complex frequency plane ("excitation factors" B_q). Using these results we compute the "excitation coefficients" C_q (essentially the mode amplitudes) in the special case where the source of the perturbations is a particle falling into the black hole along the symmetry axis. We compare this calculation with numerical integrations of the perturbation equations, and we show quantitatively how the addition of higher overtones improves the agreement with the numerical waveforms. Our results should find applications in models of the ringdown stage and in the construction of semianalytical template banks for gravitational-wave detectors, especially for binaries with large mass ratios and/or fast-spinning black holes.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

We study the influence of astrophysical formation scenarios on the precessional dynamics of spinning black-hole binaries by the time they enter the observational window of second- and third-generation gravitational-wave detectors, such as Advanced LIGO/Virgo, LIGO-India, KAGRA and the Einstein Telescope. Under the plausible assumption that tidal interactions are efficient at aligning the spins of few-solar mass black-hole progenitors with the orbital angular momentum, we find that black-hole spins should be expected to preferentially lie in a plane when they become detectable by gravitational-wave interferometers. This "resonant plane" is identified by the conditions \Delta\Phi=0{\deg} or \Delta\Phi=+/-180{\deg}, where \Delta\Phi is the angle between the components of the black-hole spins in the plane orthogonal to the orbital angular momentum. If the angles \Delta \Phi can be accurately measured for a large sample of gravitational-wave detections, their distribution will constrain models of compact binary formation. In particular, it will tell us whether tidal interactions are efficient and whether a mechanism such as mass transfer, stellar winds, or supernovae can induce a mass-ratio reversal (so that the heavier black hole is produced by the initially lighter stellar progenitor). Therefore our model offers a concrete observational link between gravitational-wave measurements and astrophysics. We also hope that it will stimulate further studies of precessional dynamics, gravitational-wave template placement and parameter estimation for binaries locked in the resonant plane.

Astrophysical tests of general relativity belong to two categories: 1) "internal", i.e. consistency tests within the theory (for example, tests that astrophysical black holes are indeed described by the Kerr solution and its perturbations), or 2) "external", i.e. tests of the many proposed extensions of the theory. I review some ways in which astrophysical black holes can be used as natural laboratories for both "internal" and "external" tests of general relativity. The examples provided here (ringdown tests of the black hole "no-hair" theorem, bosonic superradiant instabilities in rotating black holes and gravitational-wave tests of massive scalar-tensor theories) are shamelessly biased towards recent research by myself and my collaborators. Hopefully this colloquial introduction aimed mainly at astrophysicists will convince skeptics (if there are any) that space-based detectors will be crucial to study fundamental physics through gravitational-wave observations.

We discuss a general method to study linear perturbations of slowly rotating black holes which is valid for any perturbation field, and particularly advantageous when the field equations are not separable. As an illustration of the method we investigate massive vector (Proca) perturbations in the Kerr metric, which do not appear to be separable in the standard Teukolsky formalism. Working in a perturbative scheme, we discuss two important effects induced by rotation: a Zeeman-like shift of nonaxisymmetric quasinormal modes and bound states with different azimuthal number m, and the coupling between axial and polar modes with different multipolar index l. We explicitly compute the perturbation equations up to second order in rotation, but in principle the method can be extended to any order. Working at first order in rotation we show that polar and axial Proca modes can be computed by solving two decoupled sets of equations, and we derive a single master equation describing axial perturbations of spin s=0 and s=+-1. By extending the calculation to second order we can study the superradiant regime of Proca perturbations in a self-consistent way. For the first time we show that Proca fields around Kerr black holes exhibit a superradiant instability, which is significantly stronger than for massive scalar fields. Because of this instability, astrophysical observations of spinning black holes provide the tightest upper limit on the mass of the photon: mv<4x10^-20 eV under our most conservative assumptions. Spin measurements for the largest black holes could reduce this bound to mv<10^-22 eV or lower.

Generic extensions of the standard model predict the existence of ultralight bosonic degrees of freedom. Several ongoing experiments are aimed at detecting these particles or constraining their mass range. Here we show that massive vector fields around rotating black holes can give rise to a strong superradiant instability which extracts angular momentum from the hole. The observation of supermassive spinning black holes imposes limits on this mechanism. We show that current supermassive black hole spin estimates provide the tightest upper limits on the mass of the photon (mv<4x10^{-20} eV according to our most conservative estimate), and that spin measurements for the largest known supermassive black holes could further lower this bound to mv<10^{-22} eV. Our analysis relies on a novel framework to study perturbations of rotating Kerr black holes in the slow-rotation regime, that we developed up to second order in rotation, and that can be extended to other spacetime metrics and other theories.

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.

The last decade of observational and theoretical developments in stellar and binary evolution provides an opportunity to incorporate major improvements to the predictions from populations synthesis models. We compute the Galactic merger rates for NS-NS, BH-NS, and BH-BH mergers with the StarTrack code. The most important revisions include: updated wind mass loss rates (allowing for stellar mass black holes up to $80 \msun$), a realistic treatment of the common envelope phase (a process that can affect merger rates by 2--3 orders of magnitude), and a qualitatively new neutron star/black hole mass distribution (consistent with the observed "mass gap"). Our findings include: (i) The binding energy of the envelope plays a pivotal role in determining whether a binary merges within a Hubble time. (ii) Our description of natal kicks from supernovae plays an important role, especially for the formation of BH-BH systems. (iii) The masses of BH-BH systems can be substantially increased in the case of low metallicities or weak winds. (iv) Certain combinations of parameters underpredict the Galactic NS-NS merger rate, and can be ruled out. {\em (v)} Models incorporating delayed supernovae do not agree with the observed NS/BH "mass gap", in accordance with our previous work. This is the first in a series of three papers. The second paper will study the merger rates of double compact objects as a function of redshift, star formation rate, and metallicity. In the third paper we will present the detection rates for gravitational wave observatories, using up-to-date signal waveforms and sensitivity curves.

We review recent progress in numerical relativity simulations of black-hole (BH) spacetimes. Following a brief summary of the methods employed in the modeling, we summarize the key results in three major areas of BH physics: (i) BHs as sources of gravitational waves (GWs), (ii) astrophysical systems involving BHs, and (iii) BHs in high-energy physics. We conclude with a list of the most urgent tasks for numerical relativity in these three areas.

Scalar-tensor theories are among the simplest extensions of general relativity. In theories with light scalars, deviations from Einstein's theory of gravity are determined by the scalar mass m_s and by a Brans-Dicke-like coupling parameter \omega_{BD}. We show that gravitational-wave observations of nonspinning neutron star-black hole binary inspirals can be used to set lower bounds on \omega_{BD} and upper bounds on the combination m_s/\sqrt{\omega_{BD}}$. We estimate via a Fisher matrix analysis that individual observations with signal-to-noise ratio \rho would yield (m_s/\sqrt{\omega_{BD}})(\rho/10)<10^{-15}, 10^{-16} and 10^{-19} eV for Advanced LIGO, ET and eLISA, respectively. A statistical combination of multiple observations may further improve these bounds.

Recent numerical relativity simulations have shown that the final black hole produced in a binary merger can recoil with a velocity as large as 5,000 km/s. Because of enhanced gravitational-wave emission in the so-called "hang-up" configurations, this maximum recoil occurs when the black-hole spins are partially aligned with the orbital angular momentum. We revisit our previous statistical analysis of post-Newtonian evolutions of black-hole binaries in the light of these new findings. We demonstrate that despite these new configurations with enhanced recoil velocities, spin alignment during the post-Newtonian stage of the inspiral will still significantly suppress (or enhance) kick magnitudes when the initial spin of the more massive black hole is more (or less) closely aligned with the orbital angular momentum than that of the smaller hole. We present a preliminary study of how this post-Newtonian spin alignment affects the ejection probabilities of supermassive black holes from their host galaxies with astrophysically motivated mass ratio and initial spin distributions. We find that spin alignment suppresses (enhances) ejection probabilities by ~ 40% (20%) for an observationally motivated mass-dependent galactic escape velocity, and by an even greater amount for a constant escape velocity of 1,000 km/s. Kick suppression is thus at least a factor two more efficient than enhancement.

We derive the equations of motion, the periastron shift, and the gravitational radiation damping for quasicircular compact binaries in a massive variant of the Brans-Dicke theory of gravity. We also study the Shapiro time delay and the Nordtvedt effect in this theory. By comparing with recent observational data, we put bounds on the two parameters of the theory: the Brans-Dicke coupling parameter \omega_{BD} and the scalar mass m_s. We find that the most stringent bounds come from Cassini measurements of the Shapiro time delay in the Solar System, that yield a lower bound \omega_{BD}>40000 for scalar masses m_s<2.5x10^{-20} eV, to 95% confidence. In comparison, observations of the Nordtvedt effect using Lunar Laser Ranging (LLR) experiments yield \omega_{BD}>1000 for m_s<2.5x10^{-20} eV. Observations of the orbital period derivative of the quasicircular white dwarf-neutron star binary PSR J1012+5307 yield \omega_{BD}>1250 for m_s<10^{-20} eV. A first estimate suggests that bounds comparable to the Shapiro time delay may come from observations of radiation damping in the eccentric white dwarf-neutron star binary PSR J1141-6545, but a quantitative prediction requires the extension of our work to eccentric orbits.

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.

We study the coupling of massive scalar fields to matter in orbit around rotating black holes. It is generally expected that orbiting bodies will lose energy in gravitational waves, slowly inspiralling into the black hole. Instead, we show that the coupling of the field to matter leads to a surprising effect: because of superradiance, matter can hover into "floating orbits" for which the net gravitational energy loss at infinity is entirely provided by the black hole's rotational energy. Orbiting bodies remain floating until they extract sufficient angular momentum from the black hole, or until perturbations or nonlinear effects disrupt the orbit. For slowly rotating and nonrotating black holes floating orbits are unlikely to exist, but resonances at orbital frequencies corresponding to quasibound states of the scalar field can speed up the inspiral, so that the orbiting body "sinks". These effects could be a smoking gun of deviations from general relativity.

Space-based gravitational-wave detectors, such as LISA or a similar ESA-led mission, will offer unique opportunities to test general relativity. We study the bounds that space-based detectors could place on the graviton Compton wavelength \lambda_g=h/(m_g c) by observing multiple inspiralling black hole binaries. We show that while observations of individual inspirals will yield mean bounds \lambda_g~3x10^15 km, the combined bound from observing ~50 events in a two-year mission is about ten times better: \lambda_g~3x10^16 km (m_g~4x10^-26 eV). The bound improves faster than the square root of the number of observed events, because typically a few sources provide constraints as much as three times better than the mean. This result is only mildly dependent on details of black hole formation and detector characteristics. The bound achievable in practice should be one order of magnitude better than this figure (and hence almost competitive with the static, model-dependent bounds from gravitational effects on cosmological scales), because our calculations ignore the merger/ringdown portion of the waveform. The observation that an ensemble of events can sensibly improve the bounds that individual binaries set on \lambda_g applies to any theory whose deviations from general relativity are parametrized by a set of global parameters.

We develop a theoretical framework to study slowly rotating compact stars in a rather general class of alternative theories of gravity, with the ultimate goal of investigating constraints on alternative theories from electromagnetic and gravitational-wave observations of compact stars. Our Lagrangian includes as special cases scalar-tensor theories (and indirectly f(R) theories) as well as models with a scalar field coupled to quadratic curvature invariants. As a first application of the formalism, we discuss (for the first time in the literature) compact stars in Einstein-Dilaton-Gauss-Bonnet gravity. We show that compact objects with central densities typical of neutron stars cannot exist for certain values of the coupling constants of the theory. In fact, the existence and stability of compact stars sets more stringent constraints on the theory than the existence of black hole solutions. This work is a first step in a program to systematically rule out (possibly using Bayesian model selection) theories that are incompatible with astrophysical observations of compact stars.

We study the effect of black hole spin on the accuracy of the post-Newtonian approximation. We focus on the gravitational energy flux for the quasicircular, equatorial, extreme mass-ratio inspiral of a compact object into a Kerr black hole of mass M and spin J. For a given dimensionless spin a=J/M^2 (in geometrical units), the energy flux depends only on the orbital velocity v or (equivalently) on the Boyer-Lindquist orbital radius r. We investigate the formal region of validity of the Taylor post-Newtonian expansion of the energy flux (which is known up to order v^8 beyond the quadrupole formula), generalizing previous work by two of us. The "error function" used to determine the region of validity of the post-Newtonian expansion can have two qualitatively different kinds of behavior, and we deal with these two cases separately. We find that, at any fixed post-Newtonian order, the edge of the region of validity (as measured by v/v_{ISCO}, where v_{ISCO} is the orbital velocity at the innermost stable circular orbit) is only weakly dependent on a. Unlike in the nonspinning case, the lack of sufficiently high order terms does not allow us to determine if there is a convergent to divergent transition at order v^6. Independently of a, the inclusion of angular multipoles up to and including l=5 in the numerical flux is necessary to achieve the level of accuracy of the best-known (N=8) PN expansion of the energy flux.

Quantum fields in compact stars can be amplified due to a semiclassical instability. This generic feature of scalar fields coupled to curvature may affect the birth and the equilibrium structure of relativistic stars. We point out that the semiclassical instability has a classical counterpart, which occurs exactly in the same region of the parameter space. For negative values of the coupling parameter the instability is equivalent to the well-known "spontaneous scalarization" effect: the plausible end-state of the instability is a static, asymptotically flat equilibrium configuration with nonzero expectation value for the quantum fields, which is compatible with experiments in the weak-field regime and energetically favored over stellar solutions in general relativity. For positive values of the coupling parameter the new configurations are energetically disfavored, and the end-point of the instability remains an open and interesting issue. The vacuum instability may provide a natural mechanism to produce spontaneous scalarization, leading to new experimental opportunities to probe the nature of vacuum energy via astrophysical observations of compact stars.

LISA should detect gravitational waves from tens to hundreds of systems containing black holes with mass in the range from 10 thousand to 10 million solar masses. Black holes in this mass range are not well constrained by current electromagnetic observations, so LISA could significantly enhance our understanding of the astrophysics of such systems. In this paper, we describe a framework for combining LISA observations to make statements about massive black hole populations. We summarise the constraints that LISA observations of extreme-mass-ratio inspirals might be able to place on the mass function of black holes in the LISA range. We also describe how LISA observations can be used to choose between different models for the hierarchical growth of structure in the early Universe. We consider four models that differ in their prescription for the initial mass distribution of black hole seeds, and in the efficiency of accretion onto the black holes. We show that with as little as 3 months of LISA data we can clearly distinguish between these models, even under relatively pessimistic assumptions about the performance of the detector and our knowledge of the gravitational waveforms.

The massive black holes we observe in galaxies today are the natural end-product of a complex evolutionary path, in which black holes seeded in proto-galaxies at high redshift grow through cosmic history via a sequence of mergers and accretion episodes. Electromagnetic observations probe a small subset of the population of massive black holes (namely, those that are active or those that are very close to us), but planned space-based gravitational-wave observatories such as the Laser Interferometer Space Antenna (LISA) can measure the parameters of ``electromagnetically invisible'' massive black holes out to high redshift. In this paper we introduce a Bayesian framework to analyze the information that can be gathered from a set of such measurements. Our goal is to connect a set of massive black hole binary merger observations to the underlying model of massive black hole formation. In other words, given a set of observed massive black hole coalescences, we assess what information can be extracted about the underlying massive black hole population model. For concreteness we consider ten specific models of massive black hole formation, chosen to probe four important (and largely unconstrained) aspects of the input physics used in structure formation simulations: seed formation, metallicity ``feedback'', accretion efficiency and accretion geometry. For the first time we allow for the possibility of ``model mixing'', by drawing the observed population from some combination of the ``pure'' models that have been simulated. A Bayesian analysis allows us to recover a posterior probability distribution for the ``mixing parameters'' that characterize the fractions of each model represented in the observed distribution. Our work shows that LISA has enormous potential to probe the underlying physics of structure formation.

Superspinars are ultracompact objects whose mass M and angular momentum J violate the Kerr bound (cJ/GM^2>1). Recent studies analyzed the observable consequences of gravitational lensing and accretion around superspinars in astrophysical scenarios. In this paper we investigate the dynamical stability of superspinars to gravitational perturbations, considering either purely reflecting or perfectly absorbing boundary conditions at the "surface" of the superspinar. We find that these objects are unstable independently of the boundary conditions, and that the instability is strongest for relatively small values of the spin. Also, we give a physical interpretation of the various instabilities that we find. Our results (together with the well-known fact that accretion tends to spin superspinars down) imply that superspinars are very unlikely astrophysical alternatives to black holes.

Numerical-relativity simulations indicate that the black hole produced in a binary merger can recoil with a velocity up to v_max ~ 4,000 km/s with respect to the center of mass of the initial binary. This challenges the paradigm that most galaxies form through hierarchical mergers, yet retain supermassive black holes at their centers despite having escape velocities much less than v_max. Interaction with a circumbinary disk can align the binary black hole spins with their orbital angular momentum, reducing the recoil velocity of the final black hole produced in the subsequent merger. However, the effectiveness of this alignment depends on highly uncertain accretion flows near the binary black holes. In this Letter, we show that if the spin S_1 of the more massive binary black hole is even partially aligned with the orbital angular momentum L, relativistic spin precession on sub-parsec scales can align the binary black hole spins with each other. This alignment significantly reduces the recoil velocity even in the absence of gas. For example, if the angle between S_1 and L at large separations is 10 degrees while the second spin S_2 is isotropically distributed, the spin alignment discussed in this paper reduces the median recoil from 864 km/s to 273 km/s for maximally spinning black holes with a mass ratio of 9/11. This reduction will greatly increase the fraction of galaxies retaining their supermassive black holes.

The inspiral of binary black holes is governed by gravitational radiation reaction at binary separations r < 1000 M, yet it is too computationally expensive to begin numerical-relativity simulations with initial separations r > 10 M. Fortunately, binary evolution between these separations is well described by post-Newtonian equations of motion. We examine how this post-Newtonian evolution affects the distribution of spin orientations at separations r ~ 10 M where numerical-relativity simulations typically begin. Although isotropic spin distributions at r ~ 1000 M remain isotropic at r ~ 10 M, distributions that are initially partially aligned with the orbital angular momentum can be significantly distorted during the post-Newtonian inspiral. Spin precession tends to align (anti-align) the binary black hole spins with each other if the spin of the more massive black hole is initially partially aligned (anti-aligned) with the orbital angular momentum, thus increasing (decreasing) the average final spin. Spin precession is stronger for comparable-mass binaries, and could produce significant spin alignment before merger for both supermassive and stellar-mass black hole binaries. We also point out that precession induces an intrinsic accuracy limitation (< 0.03 in the dimensionless spin magnitude, < 20 degrees in the direction) in predicting the final spin resulting from the merger of widely separated binaries.

Gravitational waves from compact objects provide information about their structure, probing deep into strong-gravity regions. Here we illustrate how the presence or absence of an event horizon can produce qualitative differences in the gravitational waves emitted by ultra-compact objects. In order to set up a straw-man ultra-compact object with no event horizon, but which is otherwise almost identical to a black hole, we consider a nonrotating thin-shell model inspired by Mazur and Mottola's gravastar, which has a Schwarzschild exterior, a de Sitter interior and an infinitely thin shell with finite tension separating the two regions. As viewed from the external space-time, the shell can be located arbitrarily close to the Schwarzschild radius, so a gravastar might seem indistinguishable from a black hole when tests are only performed on its external metric. We study the linearized dynamics of the system, and in particular the junction conditions connecting internal and external gravitational perturbations. As a first application of the formalism we compute polar and axial oscillation modes of a thin-shell gravastar. We show that the quasinormal mode spectrum is completely different from that of a black hole, even in the limit when the surface redshift becomes infinite. Polar QNMs depend on the equation of state of matter on the shell and can be used to distinguish between different gravastar models. Our calculations suggest that low-compactness gravastars could be unstable when the sound speed on the shell vs/c>0.92.

It has been suggested that rotating black holes could serve as particle colliders with arbitrarily high center-of-mass energy. Astrophysical limitations on the maximal spin, back-reaction effects and sensitivity to the initial conditions impose severe limits on the likelihood of such collisions.

Quasinormal modes are eigenmodes of dissipative systems. Perturbations of classical gravitational backgrounds involving black holes or branes naturally lead to quasinormal modes. The analysis and classification of the quasinormal spectra requires solving non-Hermitian eigenvalue problems for the associated linear differential equations. Within the recently developed gauge-gravity duality, these modes serve as an important tool for determining the near-equilibrium properties of strongly coupled quantum field theories, in particular their transport coefficients, such as viscosity, conductivity and diffusion constants. In astrophysics, the detection of quasinormal modes in gravitational wave experiments would allow precise measurements of the mass and spin of black holes as well as new tests of general relativity. This review is meant as an introduction to the subject, with a focus on the recent developments in the field.

In the Fornax dwarf spheroidal galaxy the globular clusters are distributed widely, without any significant central concentration. Oh et al. pointed out that such a distribution is paradoxical: dynamical friction effects estimated using single-component King models would have forced the globular clusters to spiral down to the center of the galaxy well within a Hubble time. This paper is devoted to a discussion of this paradox. We describe a model in which the stars of the dwarf spheroidal galaxy are embedded in a cloud of dark matter, and each of these components is specified by its own phase space distribution function. This model allows us to fit self-consistently the observed luminosity profile and the spatial variation of the velocity dispersion of the stars. This fitting yields two basic parameters, related to the central density and velocity dispersion, that characterize the phase space distribution of dark matter. The dynamical friction effects calculated on the basis of this self-consistent model are small enough that the observed spatial distribution of the globular clusters poses no difficulty, and the apparent paradox is resolved. Thus we have at hand a model for Fornax that reproduces the main observed features of this dwarf spheroidal galaxy.

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.

Geodesic motion determines important features of spacetimes. Null unstable geodesics are closely related to the appearance of compact objects to external observers and have been associated with the characteristic modes of black holes. By computing the Lyapunov exponent, which is the inverse of the instability timescale associated with this geodesic motion, we show that, in the eikonal limit, quasinormal modes of black holes in any dimensions are determined by the parameters of the circular null geodesics. This result is independent of the field equations and only assumes a stationary, spherically symmetric and asymptotically flat line element, but it does not seem to be easily extendable to anti-de Sitter spacetimes. We further show that (i) in spacetime dimensions greater than four, equatorial circular timelike geodesics in a Myers-Perry black hole background are unstable, and (ii) the instability timescale of equatorial null geodesics in Myers-Perry spacetimes has a local minimum for spacetimes of dimension d > 5.

We study the transition from inspiral to plunge in general relativity by computing gravitational waveforms of non-spinning, equal-mass black-hole binaries. We consider three sequences of simulations, starting with a quasi-circular inspiral completing 1.5, 2.3 and 9.6 orbits, respectively, prior to coalescence of the holes. For each sequence, the binding energy of the system is kept constant and the orbital angular momentum is progressively reduced, producing orbits of increasing eccentricity and eventually a head-on collision. We analyze in detail the radiation of energy and angular momentum in gravitational waves, the contribution of different multipolar components and the final spin of the remnant. We find that the motion transitions from inspiral to plunge when the orbital angular momentum L=L_crit is about 0.8M^2. For L<L_crit the radiated energy drops very rapidly. Orbits with L of about L_crit produce our largest dimensionless Kerr parameter for the remnant, j=J/M^2=0.724. Generalizing a model recently proposed by Buonanno, Kidder and Lehner to eccentric binaries, we conjecture that (1) j=0.724 is the maximal Kerr parameter that can be obtained by any merger of non-spinning holes, and (2) no binary merger (even if the binary members are extremal Kerr black holes with spins aligned to the orbital angular momentum, and the inspiral is highly eccentric) can violate the cosmic censorship conjecture.

We study the transformation of multipolar decompositions of gravitational radiation under rotations and boosts. Rotations to the remnant black hole's frame simplify the waveforms from the merger of generic spinning black hole binaries. Boosts may be important to get an accurate gravitational-wave phasing, especially for configurations leading to large recoil velocities of the remnant. As a test of our formalism we revisit the classic problem of point particles falling into a Schwarzschild black hole. Then we highlight by specific examples the importance of choosing the right frame in numerical simulations of unequal-mass, spinning binary black-hole mergers.

Using recent results from numerical relativity simulations of black hole mergers, we revisit previous studies of cosmological black hole spin evolution. We show that mergers are very unlikely to yield large spins, unless alignment of the spins of the merging holes with the orbital angular momentum is very efficient. We analyze the spin evolution in three specific scenarios: (1) spin evolves only through mergers, (2) spin evolves through mergers and prolonged accretion episodes, (3) spin evolves through mergers and short-lived (chaotic) accretion episodes. We study how different diagnostics can distinguish between these evolutionary scenarios, assessing the discriminating power of gravitational-wave measurements and X-ray spectroscopy. Gravitational radiation can produce three different types of spin measurements, yielding respectively the spins of the two black holes in a binary inspiral prior to merger, the spin of the merger remnant (as encoded in the ringdown waves), and the spin of ``single'' black holes during the extreme mass-ratio inspiral (EMRI) of compact objects. The latter spin population is also accessible to iron-line measurements. We compute and compare the spin distributions relevant for these different observations. If iron-line measurements and gravitational-wave observations of EMRIs only yield dimensionless spins j=J/M^2>0.9, then prolonged accretion should be responsible for spin-up, and chaotic accretion scenarios would be very unlikely. If only a fraction of the whole population of low-redshift black holes spins rapidly, spin-alignment during binary mergers (rather than prolonged accretion) could be responsible for spin-ups.

Using recent results from numerical relativity simulations of non-spinning binary black hole mergers we revisit the problem of detecting ringdown waveforms and of estimating the source parameters, considering both LISA and Earth-based interferometers. We find that Advanced LIGO and EGO could detect intermediate-mass black holes of mass up to about 1000 solar masses out to a luminosity distance of a few Gpc. For typical multipolar energy distributions, we show that the single-mode ringdown templates presently used for ringdown searches in the LIGO data stream can produce a significant event loss (> 10% for all detectors in a large interval of black hole masses) and very large parameter estimation errors on the black hole's mass and spin. We estimate that more than 10^6 templates would be needed for a single-stage multi-mode search. Therefore, we recommend a "two stage" search to save on computational costs: single-mode templates can be used for detection, but multi-mode templates or Prony methods should be used to estimate parameters once a detection has been made. We update estimates of the critical signal-to-noise ratio required to test the hypothesis that two or more modes are present in the signal and to resolve their frequencies, showing that second-generation Earth-based detectors and LISA have the potential to perform no-hair tests.

We study the inspiral, merger and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q=M2/M1=1 to q=4). We compare numerical and Post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l,m). We find that the Post-Newtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l,m), and that the convergence of the Post-Newtonian series to the numerical results is non-monotonic. To leading order the total energy emitted in the merger phase scales like eta^2 and the spin of the final black hole scales like eta, where eta=q/(1+q)^2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By cross-correlating information from different multipoles we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations.

We use a post-Newtonian diagnostic tool to examine numerically generated quasiequilibrium initial data sets for non-spinning double neutron star and neutron star-black hole binary systems. The PN equations include the effects of tidal interactions, parametrized by the compactness of the neutron stars and by suitable values of ``apsidal'' constants, which measure the degree of distortion of stars subjected to tidal forces. We find that the post-Newtonian diagnostic agrees well with the double neutron star initial data, typically to better than half a percent except where tidal distortions are becoming extreme. We show that the differences could be interpreted as representing small residual eccentricity in the initial orbits. In comparing the diagnostic with preliminary numerical data on neutron star-black hole binaries, we find less agreement.

The ringdown phase following a binary black hole merger is usually assumed to be well described by a linear superposition of complex exponentials (quasinormal modes). In the strong-field conditions typical of a binary black hole merger, non-linear effects may produce mode coupling. Artificial mode coupling can also be induced by the black hole's rotation, if the radiation field is expanded in terms of spin-weighted spherical (rather than spheroidal) harmonics. Observing deviations from linear black hole perturbation theory requires optimal fitting techniques to extract ringdown parameters from numerical waveforms, which are inevitably affected by errors. So far, non-linear least-squares fitting methods have been used as the standard workhorse to extract frequencies from ringdown waveforms. These methods are known not to be optimal for estimating parameters of complex exponentials. Furthermore, different fitting methods have different performance in the presence of noise. The main purpose of this paper is to introduce the gravitational wave community to modern variations of a linear parameter estimation technique first devised in 1795 by Prony: the Kumaresan-Tufts and matrix pencil methods. Using "test" damped sinusoidal signals in Gaussian white noise we illustrate the advantages of these methods, showing that they have variance and bias at least comparable to standard non-linear least-squares techniques. Then we compare the performance of different methods on unequal-mass binary black hole merger waveforms. The methods we discuss should be useful both theoretically (to monitor errors and search for non-linearities in numerical relativity simulations) and experimentally (for parameter estimation from ringdown signals after a gravitational wave detection).

Distorted black holes radiate gravitational waves. In the so-called ringdown phase radiation is emitted in a discrete set of complex quasinormal frequencies, whose values depend only on the black hole's mass and angular momentum. Ringdown radiation could be detectable with large signal-to-noise ratio by the Laser Interferometer Space Antenna LISA. If more than one mode is detected, tests of the black hole nature of the source become possible. The detectability of different modes depends on their relative excitation, which in turn depends on the cause of the perturbation (i.e. on the initial data). A ``universal'', initial data-independent measure of the relative mode excitation is encoded in the poles of the Green's function that propagates small perturbations of the geometry (``excitation factors''). We compute for the first time the excitation factors for general-spin perturbations of Kerr black holes. We find that for corotating modes with $l=m$ the excitation factors tend to zero in the extremal limit, and that the contribution of the overtones should be more significant when the black hole is fast rotating. We also present the first analytical calculation of the large-damping asymptotics of the excitation factors for static black holes, including the Schwarzschild and Reissner-Nordstrom metrics. This is an important step to determine the convergence properties of the quasinormal mode expansion.

We present numerical results from three-dimensional evolutions of scalar perturbations of Kerr black holes. Our simulations make use of a high-order accurate multi-block code which naturally allows for fixed adaptivity and smooth inner (excision) and outer boundaries. We focus on the quasinormal ringing phase, presenting a systematic method for extraction of the quasinormal mode frequencies and amplitudes and comparing our results against perturbation theory. The amplitude of each mode depends exponentially on the starting time of the quasinormal regime, which is not defined unambiguously. We show that this time-shift problem can be circumvented by looking at appropriately chosen relative mode amplitudes. From our simulations we extract the quasinormal frequencies and the relative and absolute amplitudes of corotating and counterrotating modes (including overtones in the corotating case). We study the dependence of these amplitudes on the shape of the initial perturbation, the angular dependence of the mode and the black hole spin, comparing against results from perturbation theory in the so-called asymptotic approximation. We also compare the quasinormal frequencies from our numerical simulations with predictions from perturbation theory, finding excellent agreement. Finally we study under what conditions the relative amplitude between given pairs of modes gets maximally excited and present a quantitative analysis of rotational mode-mode coupling. The main conclusions and techniques of our analysis are quite general and, as such, should be of interest in the study of ringdown gravitational waves produced by astrophysical gravitational wave sources.

The evidence for supermassive Kerr black holes in galactic centers is strong and growing, but only the detection of gravitational waves will convincingly rule out other possibilities to explain the observations. The Kerr spacetime is completely specified by the first two multipole moments: mass and angular momentum. This is usually referred to as the ``no-hair theorem'', but it is really a ``two-hair'' theorem. If general relativity is the correct theory of gravity, the most plausible alternative to a supermassive Kerr black hole is a rotating boson star. Numerical calculations indicate that the spacetime of rotating boson stars is determined by the first three multipole moments (``three-hair theorem''). LISA could accurately measure the oscillation frequencies of these supermassive objects. We propose to use these measurements to ``count their hair'', unambiguously determining their nature and properties.

The observability of gravitational waves from supermassive and intermediate-mass black holes by the forecoming Laser Interferometer Space Antenna (LISA), and the physics we can learn from the observations, will depend on two basic factors: the event rates for massive black hole mergers occurring in the LISA best sensitivity window, and our theoretical knowledge of the gravitational waveforms. We first provide a concise review of the literature on LISA event rates for massive black hole mergers, as predicted by different formation scenarios. Then we discuss what (in our view) are the most urgent issues to address in terms of waveform modelling. For massive black hole binary inspiral these include spin precession, eccentricity, the effect of high-order Post-Newtonian terms in the amplitude and phase, and an accurate prediction of the transition from inspiral to plunge. For black hole ringdown, numerical relativity will ultimately be required to determine the relative quasinormal mode excitation, and to reduce the dimensionality of the template space in matched filtering.

Newly formed black holes are expected to emit characteristic radiation in the form of quasi-normal modes, called ringdown waves, with discrete frequencies. LISA should be able to detect the ringdown waves emitted by oscillating supermassive black holes throughout the observable Universe. We develop a multi-mode formalism, applicable to any interferometric detectors, for detecting ringdown signals, for estimating black hole parameters from those signals, and for testing the no-hair theorem of general relativity. Focusing on LISA, we use current models of its sensitivity to compute the expected signal-to-noise ratio for ringdown events, the relative parameter estimation accuracy, and the resolvability of different modes. We also discuss the extent to which uncertainties on physical parameters, such as the black hole spin and the energy emitted in each mode, will affect our ability to do black hole spectroscopy.

We provide some considerations on the excitation of black hole quasinormal modes (QNMs) in different physical scenarios. Considering a simple model in which a stream of particles accretes onto a black hole, we show that resonant QNM excitation by hyperaccretion requires a significant amount of fine-tuning, and is quite unlikely to occur in nature. Then we summarize and discuss present estimates of black hole QNM excitation from gravitational collapse, distorted black holes and head-on black hole collisions. We emphasize the areas that, in our opinion, are in urgent need of further investigation from the point of view of gravitational wave source modeling.

Observations of binary inspirals with LISA will allow us to place bounds on alternative theories of gravity and to study the merger history of massive black holes (MBH). These possibilities rely on LISA's parameter estimation accuracy. We update previous studies of parameter estimation including non-precessional spin effects. We work both in Einstein's theory and in alternative theories of gravity of the scalar-tensor and massive-graviton types. Inclusion of non-precessional spin terms in MBH binaries has little effect on the angular resolution or on distance determination accuracy, but it degrades the estimation of the chirp mass and reduced mass by between one and two orders of magnitude. The bound on the coupling parameter of scalar-tensor gravity is significantly reduced by the presence of spin couplings, while the reduction in the graviton-mass bound is milder. LISA will measure the luminosity distance of MBHs to better than ~10% out to z~4 for a (10^6+10^6) Msun binary, and out to z~2 for a (10^7+10^7) Msun binary. The chirp mass of a MBH binary can always be determined with excellent accuracy. Ignoring spin effects, the reduced mass can be measured within ~1% out to z=10 and beyond for a (10^6+10^6) Msun binary, but only out to z~2 for a (10^7+10^7) Msun binary. Present-day MBH coalescence rate calculations indicate that most detectable events should originate at z~2-6: at these redshifts LISA can be used to measure the two black hole masses and their luminosity distance with sufficient accuracy to probe the merger history of MBHs. If the low-frequency LISA noise can only be trusted down to 10^-4 Hz, parameter estimation for MBHs (and LISA's ability to perform reliable cosmological observations) will be significantly degraded.

We compute numerically the quasinormal modes of Kerr-Newman black holes in the scalar case, for which the perturbation equations are separable. Then we study different approximations to decouple electromagnetic and gravitational perturbations of the Kerr-Newman metric, computing the corresponding quasinormal modes. Our results suggest that the Teukolsky-like equation derived by Dudley and Finley gives a good approximation to the dynamics of a rotating charged black hole for Q<M/2. Though insufficient to deal with Kerr-Newman based models of elementary particles, the Dudley-Finley equation should be adequate for astrophysical applications.

Higher-dimensional scenarios allow for the formation of mini-black holes from TeV-scale particle collisions. The purpose of this paper is to review and compare different methods for the estimate of the total gravitational energy emitted in this process. To date, black hole formation has mainly been studied using an apparent horizon search technique. This approach yields only an upper bound on the gravitational energy emitted during black hole formation. Alternative calculations based on instantaneous collisions of point particles and black hole perturbation theory suggest that the emitted gravitational energy may be smaller. New and more refined methods may be necessary to accurately describe black hole formation in high-energy particle collisions.

We investigate the effect of spin-orbit and spin-spin couplings on the estimation of parameters for inspiralling compact binaries of massive black holes, and for neutron stars inspiralling into intermediate-mass black holes, using hypothetical data from the proposed Laser Interferometer Space Antenna (LISA). We work both in Einstein's theory and in alternative theories of gravity of the scalar-tensor and massive-graviton types. We restrict the analysis to non-precessing spinning binaries, i.e. to cases where the spins are aligned normal to the orbital plane. We find that the accuracy with which intrinsic binary parameters such as chirp mass and reduced mass can be estimated within general relativity is degraded by between one and two orders of magnitude. We find that the bound on the coupling parameter omega_BD of scalar-tensor gravity is significantly reduced by the presence of spin couplings, while the reduction in the graviton-mass bound is milder. Using fast Monte-Carlo simulations of 10^4 binaries, we show that inclusion of spin terms in massive black-hole binaries has little effect on the angular resolution or on distance determination accuracy. For stellar mass inspirals into intermediate-mass black holes, the angular resolution and the distance are determined only poorly, in all cases considered. We also show that, if LISA's low-frequency noise sensitivity can be extrapolated from 10^-4 Hz to as low as 10^-5 Hz, the accuracy of determining both extrinsic parameters (distance, sky location) and intrinsic parameters (chirp mass, reduced mass) of massive binaries may be greatly improved.

We compare three different models of rotating neutron star spacetimes: the Hartle-Thorne (HT) slow-rotation approximation at second order in rotation, the exact analytic vacuum solution of Manko et al. and a numerical solution of the full Einstein equations. We integrate the HT structure equations for five representative equations of state. Then we match the HT models to numerical solutions of the Einstein equations, imposing that the mass and angular momentum of the models be the same. We estimate the limits of validity of the HT expansion computing relative errors in the spacetime's quadrupole moment Q and in the ISCO radii. We find that ISCO radii computed in the HT approximation are accurate to better than 1%, even for the fastest observed ms pulsar. At the same rotational rates the accuracy on Q is of order 20%. In the second part of the paper we focus on the exterior vacuum spacetimes. We introduce a physically motivated `quasi-Kinnersley' Newman-Penrose frame. In this frame we evaluate the speciality index S, a coordinate-independent quantity measuring the deviation of each model from Petrov Type D. On the equatorial plane this deviation is smaller than 5%, even for the fastest rotating models. Our main conclusion is that the HT approximation is very reliable for most astrophysical applications.

Many properties of black holes can be studied using acoustic analogues in the laboratory through the propagation of sound waves. We investigate in detail sound wave propagation in a rotating acoustic (2+1)-dimensional black hole, which corresponds to the ``draining bathtub'' fluid flow. We compute the quasinormal mode frequencies of this system and discuss late-time power-law tails. Due to the presence of an ergoregion, waves in a rotating acoustic black hole can be superradiantly amplified. We compute superradiant reflection coefficients and instability timescales for the acoustic black hole bomb, the equivalent of the Press-Teukolsky black hole bomb. Finally we discuss quasinormal modes and late-time tails in a non-rotating canonical acoustic black hole, corresponding to an incompressible, spherically symmetric (3+1)-dimensional fluid flow.

This is a short review of the quasinormal mode spectrum of Schwarzschild, Reissner-Nordstrom and Kerr black holes. The summary includes previously unpublished calculations of i) the eigenvalues of spin-weighted spheroidal harmonics, and ii) quasinormal frequencies of extremal Reissner-Nordstrom black holes.

We investigate the properties of a closed-form analytic solution recently found by Manko et al. (2000) for the exterior spacetime of rapidly rotating neutron stars. For selected equations of state we numerically solve the full Einstein equations to determine the neutron star spacetime along constant rest mass sequences. The analytic solution is then matched to the numerical solutions by imposing the condition that the quadrupole moment of the numerical and analytic spacetimes be the same. For the analytic solution we consider, such a matching condition can be satisfied only for very rapidly rotating stars. When solutions to the matching condition exist, they belong to one of two branches. For one branch the current octupole moment of the analytic solution is very close to the current octupole moment of the numerical spacetime; the other branch is more similar to the Kerr solution. We present an extensive comparison of the radii of innermost stable circular orbits (ISCOs) obtained with a) the analytic solution, b) the Kerr metric, c) an analytic series expansion derived by Shibata and Sasaki (1998) and d) a highly accurate numerical code. In most cases where a corotating ISCO exists, the analytic solution has an accuracy consistently better than the Shibata-Sasaki expansion. The numerical code is used for tabulating the mass-quadrupole and current-octupole moments for several sequences of constant rest mass.

We compute for the first time very highly damped quasinormal modes of the (rotating) Kerr black hole. Our numerical technique is based on a decoupling of the radial and angular equations, performed using a large-frequency expansion for the angular separation constant_{s}A_{l m}. This allows us to go much further in overtone number than ever before. We find that the real part of the quasinormal frequencies approaches a non-zero constant value which does not depend on the spin s of the perturbing field and on the angular index l: \omega_R=m\varpi(a). We numerically compute \varpi(a). Leading-order corrections to the asymptotic frequency are likely to be of order 1/\omega_I. The imaginary part grows without bound, the spacing between consecutive modes being a monotonic function of a.

Motivated by recent suggestions that highly damped black hole quasinormal modes (QNM's) may provide a link between classical general relativity and quantum gravity, we present an extensive computation of highly damped QNM's of Kerr black holes. We do not limit our attention to gravitational modes, thus filling some gaps in the existing literature. The frequency of gravitational modes with l=m=2 tends to \omega_R=2 \Omega, \Omega being the angular velocity of the black hole horizon. If Hod's conjecture is valid, this asymptotic behaviour is related to reversible black hole transformations. Other highly damped modes with m>0 that we computed do not show a similar behaviour. The real part of modes with l=2 and m<0 seems to asymptotically approach a constant value \omega_R\simeq -m\varpi, \varpi\simeq 0.12 being (almost) independent of a. For any perturbing field, trajectories in the complex plane of QNM's with m=0 show a spiralling behaviour, similar to the one observed for Reissner-Nordstrom (RN) black holes. Finally, for any perturbing field, the asymptotic separation in the imaginary part of consecutive modes with m>0 is given by 2\pi T_H (T_H being the black hole temperature). We conjecture that for all values of l and m>0 there is an infinity of modes tending to the critical frequency for superradiance (\omega_R=m) in the extremal limit. Finally, we study in some detail modes branching off the so--called ``algebraically special frequency'' of Schwarzschild black holes. For the first time we find numerically that QNM multiplets emerge from the algebraically special Schwarzschild modes, confirming a recent speculation.

We study the stability of five-dimensional Myers-Perry black holes with a single angular momentum under linear perturbations, and we compute the quasinormal modes (QNM's) of the black hole metric projected on the brane, using Leaver's continued fraction method. In our numerical search we do not find unstable modes. The damping time of modes having l=m=2 and l=m=1 tends to infinity as the black hole spin tends to the extremal value, showing a behaviour reminiscent of the one observed for ordinary 4-dimensional Kerr black holes.

According to a recent proposal, the so-called Barbero-Immirzi parameter of Loop Quantum Gravity can be fixed, using Bohr's correspondence principle, from a knowledge of highly-damped black hole oscillation frequencies. Such frequencies are rather difficult to compute, even for Schwarzschild black holes. However, it is now quite likely that they may provide a fundamental link between classical general relativity and quantum theories of gravity. Here we carry out the first numerical computation of very highly damped quasinormal modes (QNM's) for charged and rotating black holes. In the Reissner-Nordstr\"om case QNM frequencies and damping times show an oscillatory behaviour as a function of charge. The oscillations become faster as the mode order increases. At fixed mode order, QNM's describe spirals in the complex plane as the charge is increased, tending towards a well defined limit as the hole becomes extremal. Kerr QNM's have a similar oscillatory behaviour when the angular index $m=0$. For $l=m=2$ the real part of Kerr QNM frequencies tends to $2\Omega$, $\Omega$ being the angular velocity of the black hole horizon, while the asymptotic spacing of the imaginary parts is given by $2\pi T_H$.

We study scalar, electromagnetic and gravitational perturbations of a Reissner-Nordstr\"om-anti-de Sitter (RN-AdS) spacetime, and compute its quasinormal modes (QNM's). We confirm and extend results previously found for Schwarzschild-anti-de Sitter (S-AdS) black holes. For ``large'' black holes, whose horizon is much larger than the AdS radius, different classes of perturbations are almost exactly {\it isospectral}; this isospectrality is broken when the black hole's horizon radius is comparable to the AdS radius. We provide very accurate fitting formulas for the QNM's, which are valid for black holes of any size and charge $Q<Q_{ext}/3$. Electromagnetic and axial perturbations of large black holes are characterized by the existence of pure-imaginary (purely damped) modes. The damping of these modes tends to infinity as the black hole charge approaches the extremal value; if the corresponding mode amplitude does not tend to zero in the same limit, this implies that {\it extremally charged RN-AdS black holes are marginally unstable}. This result is relevant in view of the AdS/CFT conjecture, since, according to it, the AdS QNM's give the timescales for approach to equilibrium in the corresponding conformal field theory.

We compute the overlap function between Post-Newtonian (PN) templates and gravitational signals emitted by binary systems composed of one neutron star and one point mass, obtained by a perturbative approach. The calculations are performed for different stellar models and for different detectors, to estimate how effectual and faithful the PN templates are, and to establish whether effects related to the internal structure of neutron stars may possibly be extracted by the matched filtering technique.

A phase transition occurring in the inner core of a neutron star could be associated to a density discontinuity that would affect the frequency spectrum of the non-radial oscillation modes in two ways. Firstly, it would produce a softening of the equation of state, leading to more compact equilibrium configurations and changing the frequency of the fundamental and pressure modes of the neutron star. Secondly, a new non-zero frequency g-- mode would appear, associated to each discontinuity. These discontinuity g--modes have typical frequencies larger than those of g--modes previously studied in the literature (thermal, core g-- modes, or g--modes due to chemical inhomogeneities in the outer layers), and smaller than that of the fundamental mode; therefore they should be distinguishable from the other modes of non radial oscillation. In this paper we investigate how high density discontinuities change the frequency spectrum of the non-radial oscillations, in the framework of the general relativistic theory of stellar perturbations. Our purpose is to understand whether a gravitational signal, emitted at the frequencies of the quasi normal modes, may give some clear information on the equation of state of the neutron star and, in particular, on the parameters that characterize the density discontinuity. We discuss some astrophysical processes that may be associated to the excitation of these modes, and estimate how much gravitational energy should the modes convey to produce a signal detectable by high frequency gravitational detectors.

The effects that the structure of a neutron star would have on the gravitational emission of a binary system are studied in a perturbative regime, and in the frequency domain. Assuming that a neutron star is perturbed by a point mass moving on a close, circular orbit, we solve the equations of stellar perturbations in general relativity to evaluate the energy lost by the system in gravitational waves. We compare the energy output obtained for different stellar models with that found by assuming that the perturbed object is a black hole with the same mass, and we discuss the role played by the excitation of the stellar modes. Ouresults indicate that the stellar structure begins to affect the emitted power when the orbital velocity is v >0.2c (about 185 Hz for a binary system composed of two canonical neutron stars). We show that the differences between different stellar models and a black hole are due mainly to the excitation of the quasinormal modes of the star. Finally, we discuss to what extent and up to which distance the perturbative approach can be used to describe the interaction of a star and a pointlike massive body.

We compute the energy spectra of the gravitational signals emitted when a pointlike mass moves on a closed orbit around a non rotating neutron star, inducing a perturbation of its gravitational field and its internal structure. The Einstein equations and the hydrodynamical equations are perturbed and numerically integrated in the frequency domain. The results are compared with the energy spectra computed by the quadrupole formalism which assumes that both masses are pointlike, and accounts only for the radiation emitted because the orbital motion produces a time dependent quadrupole moment. The results of our perturbative approach show that, in general, the quadrupole formalism overestimates the amount of emitted radiation, especially when the two masses are close. However, if the pointlike mass is allowed to move on an orbit so tight that the keplerian orbital frequency resonates with the frequency of the fundamental quasi-normal mode of the star (2w_K=w_f), this mode can be excited and the emitted radiation can be considerably larger than that computed by the quadrupole approach.

In this paper we consider the Extra-solar Planetary Systems recently discovered in our Galaxy as potential sources of gravitational waves. We estimate the frequency and characteristic amplitude of the radiation they emit due to the orbital motion, using the quadrupole formalism. In addition, we check whether the conditions needed for the resonant excitation of the f- and g-modes of the central star can be fulfilled. By a Roche-lobe analysis, we show that there could exist systems in which the low-order g-modes could be excited, although this does not happen in the systems discovered up to now.

The possibility of exciting the g-modes of a solar-type star as a consequence of the gravitational interaction with a close companion (a planet or a brown dwarf) is studied by a perturbative approach. The amplitude of the emitted gravitational wave is computed and compared with the quadrupole emission of the system, showing that in some cases it can be considerably larger. The effects of radiation reaction are considered to evaluate the timescale of the emission process, and a Roche lobe analysis is used to establish the region where the companion can orbit without being disrupted by tidal interactions with the star.