Source author record

Emanuele Berti

Emanuele Berti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

108works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

108 published item(s)

preprint2026arXiv

Extreme mass ratio inspirals in dark matter halos: dynamics and distinguishability of halo models

The gravitational wave (GW) signals from extreme mass-ratio inspirals (EMRIs), a key target for the Laser Interferometer Space Antenna (LISA), will be affected in the presence of dark matter (DM) halos. In this paper we explore whether the effects of DM are detectable by LISA within a fully relativistic framework. We model the massive EMRI component as a nonrotating black hole (BH) surrounded by a DM halo. We compute axial and polar GW fluxes for circular orbits at linear order in the mass ratio for DM density profiles with varying mass and compactness. By comparing the phase evolution with vacuum systems, we find that DM halos can induce dephasings of tens to hundreds of radians over a one-year observation period. We demonstrate that even highly diluted DM distributions can significantly affect the emitted waveforms, and that the resulting GW signals can usually be distinguished from each other. While it is important to generalize these findings to more generic orbits and to spinning BHs, our results suggest that LISA could not only reveal the presence of DM halos, but also discriminate between different halo models.

preprint2022arXiv

Analysis of Ringdown Overtones in GW150914

We analyze GW150914 post-merger data to understand if ringdown overtone detection claims are robust. We find no evidence in favor of an overtone in the data after the waveform peak. Around the peak, the Bayes factor does not indicate the presence of an overtone, while the support for a nonzero amplitude is sensitive to changes in the starting time much smaller than the overtone damping time. This suggests that claims of an overtone detection are noise-dominated. We perform GW150914-like injections in neighboring segments of the real detector noise, and we show that noise can indeed induce artificial evidence for an overtone.

preprint2022arXiv

Black hole sensitivities in Einstein-scalar-Gauss-Bonnet gravity

The post-Newtonian dynamics of black hole binaries in Einstein-scalar-Gauss-Bonnet theories of gravity depends on the so-called "sensitivities", quantities which characterize a black hole's adiabatic response to the time-dependent scalar field environment sourced by its companion. In this work, we calculate numerically the sensitivities of nonrotating black holes, including spontaneously scalarized ones, in three classes of Einstein-scalar-Gauss-Bonnet gravity: the shift-symmetric, dilatonic and Gaussian theories. When possible, we compare our results against perturbative analytical results, finding excellent agreement. Unlike their general relativistic counterparts, black holes in Einstein-scalar-Gauss-Bonnet gravity only exist in a restricted parameter space controlled by the theory's coupling constant. A preliminary study of the role played by the sensitivities in black hole binaries suggests that, in principle, black holes can be driven outside of their domain of existence during the inspiral, for binary parameters which we determine.

preprint2022arXiv

Compact Binary Foreground Subtraction in Next-Generation Ground-Based Observatories

The stochastic gravitational-wave backgrounds (SGWBs) for current detectors are dominated by binary black-hole (BBH) and binary neutron-star (BNS) coalescences. The sensitivity of current networks of gravitational-wave (GW) detectors allows only a small fraction of BBHs and BNSs to be resolved and subtracted, but previous work indicated that the situation should significantly improve with next-generation (XG) observatories. We revisit these conclusions by taking into account waveform-modeling uncertainties, updated astrophysical models, and (crucially) the full set of parameters that must be estimated to remove the resolved sources. Compared to previous studies, we find that the residual background from BBHs and BNSs is large even with XG detector networks. New data analysis methods will thus be required to observe the SGWB from cosmic supernovae or contributions from early-Universe phenomena like cosmic strings, stiff post-inflation fluids, or axion inflation.

preprint2022arXiv

Constraining high-redshift stellar-mass primordial black holes with next-generation ground-based gravitational-wave detectors

The possible existence of primordial black holes in the stellar mass window has received considerable attention because their mergers may contribute to current and future gravitational-wave detections. Primordial black hole mergers, together with mergers of black holes originating from Population~III stars, are expected to dominate at high redshifts ($z\gtrsim 10$). However the primordial black hole merger rate density is expected to rise monotonically with redshift, while Population~III mergers can only occur after the birth of the first stars. Next-generation gravitational-wave detectors such as Cosmic Explorer~(CE) and Einstein Telescope~(ET) can access this distinctive feature in the merger rates as functions of redshift, allowing for a direct measurement of the abundance of the two populations, and hence for robust constraints on the abundance of primordial black holes. We simulate four-months worth of data observed by a CE-ET detector network and perform hierarchical Bayesian analysis to recover the merger rate densities. We find that if the Universe has no primordial black holes with masses of $\mathcal{O}(10M_{\odot})$, the projected upper limit on their abundance $f_{\rm PBH}$ as a fraction of dark matter energy density may be as low as $f_{\rm PBH}\sim \mathcal{O}({10^{-5}})$, about two orders of magnitude lower than current upper limits in this mass range. If instead $f_{\rm PBH}\gtrsim 10^{-4}$, future gravitational wave observations would exclude $f_{\rm PBH}=0$ at the 95\% credible interval.

preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

preprint2022arXiv

Destabilizing the Fundamental Mode of Black Holes: The Elephant and the Flea

Recent work applying the notion of pseudospectrum to gravitational physics showed that the quasinormal mode spectrum of black holes is unstable, with the possible exception of the longest-lived (fundamental) mode. The fundamental mode dominates the expected signal in gravitational wave astronomy, and there is no reason why it should have privileged status. We compute the quasinormal mode spectrum of two model problems where the Schwarzschild potential is perturbed by a small "bump" consisting of either a Pöschl-Teller potential or a Gaussian, and we show that the fundamental mode is destabilized under generic perturbations. We present phase diagrams and study a simple double-barrier toy problem to clarify the conditions under which the spectral instability occurs.

preprint2022arXiv

How to assess the primordial origin of single gravitational-wave events with mass, spin, eccentricity, and deformability measurements

A population of primordial black holes formed in the early Universe could contribute to at least a fraction of the black-hole merger events detectable by current and future gravitational-wave interferometers. With the ever-increasing number of detections, an important open problem is how to discriminate whether a given event is of primordial or astrophysical origin. We systematically present a comprehensive and interconnected list of discriminators that would allow us to rule out, or potentially claim, the primordial origin of a binary by measuring different parameters, including redshift, masses, spins, eccentricity, and tidal deformability. We estimate how accurately future detectors (such as the Einstein Telescope and LISA) could measure these quantities, and we quantify the constraining power of each discriminator for current interferometers. We apply this strategy to the GWTC-3 catalog of compact binary mergers. We show that current measurement uncertainties do not allow us to draw solid conclusions on the primordial origin of individual events, but this may become possible with next-generation ground-based detectors.

preprint2022arXiv

Hunting intermediate-mass black holes with LISA binary radial velocity measurements

Despite their potential role as massive seeds for quasars, in dwarf galaxy feedback, and in tidal disruption events, the observational evidence for intermediate-mass black holes (IMBHs) is scarce. LISA may observe stellar-mass black hole binaries orbiting Galactic IMBHs, and reveal the presence of the IMBH by measuring the Doppler shift in the gravitational waveform induced by the binary's radial velocity. We estimate the number of detectable Doppler shift events from the Milky Way globular clusters (assuming they host IMBHs) and we find that it decreases with the IMBH mass. A few Galactic globular clusters (including M22 and $ω$ Centauri) may produce at least one event detectable by LISA. Even in more pessimistic scenarios, one could still expect $\sim$ 1 event overall in the Milky Way. We also estimate the number of Doppler shift events for IMBHs wandering in the Milky Way as a result of the disruption of their parent clusters. If there is at least one binary black hole orbiting around each wandering IMBH, LISA may detect up to a few tens of Doppler shift events from this elusive IMBH population. Under more pessimistic assumptions, LISA may still detect $\sim 1$ wandering IMBH that would hardly be observable otherwise.

preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

preprint2022arXiv

Searching for a subpopulation of primordial black holes in LIGO/Virgo gravitational-wave data

With several dozen binary black hole events detected by LIGO/Virgo to date and many more expected in the next few years, gravitational-wave astronomy is shifting from individual-event analyses to population studies. Using the GWTC-2 catalog, we perform a hierarchical Bayesian analysis that for the first time combines several state-of-the-art astrophysical formation models with a population of primordial black holes (PBHs) and constrains the fraction of a putative subpopulation of PBHs in the data. We find that this fraction depends significantly on the set of assumed astrophysical models. While a primordial population is statistically favored against certain competitive astrophysical channels, such as globular clusters and nuclear stellar clusters, a dominant contribution from the stable-mass-transfer isolated formation channel drastically reduces the need for PBHs, except for explaining the rate of mass-gap events like GW190521. The tantalizing possibility that black holes formed after inflation are contributing to LIGO/Virgo observations could only be verified by further reducing uncertainties in astrophysical and primordial formation models, and it may ultimately be confirmed by third-generation interferometers.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Fundamental Physics and Beyond the Standard Model

Gravitational wave detectors are formidable tools to explore strong-field gravity, especially black holes and neutron stars. These compact objects are extraordinarily efficient at producing electromagnetic and gravitational radiation. As such, they are ideal laboratories for fundamental physics and have an immense discovery potential. The detection of black hole binaries by third-generation Earth-based detectors, space-based detectors and pulsar timing arrays will provide exquisite tests of general relativity. Loud "golden" events and extreme mass-ratio inspirals can strengthen the observational evidence for horizons by mapping the exterior spacetime geometry, inform us on possible near-horizon modifications, and perhaps reveal a breakdown of Einstein's gravity. Measurements of the black-hole spin distribution and continuous gravitational-wave searches can turn black holes into efficient detectors of ultralight bosons across ten or more orders of magnitude in mass. A precise monitoring of the phase of inspiralling binaries can constrain the existence of additional propagating fields and characterize the environment in which the binaries live, bounding the local dark matter density and properties. Gravitational waves from compact binaries will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address fundamental issues in our current understanding of the cosmos.

preprint2022arXiv

Superradiance in massive vector fields with spatially varying mass

Superradiance is a process by which massive bosonic particles can extract energy from spinning black holes, leading to the build up of a "cloud" if the particle has a Compton wavelength comparable to the black hole's Schwarzschild radius. One interesting possibility is that superradiance may occur for photons in a diffuse plasma, where they gain a small effective mass. Studies of the spin-0 case have indicated that such a build up is suppressed by a spatially varying effective mass, supposed to mimic the photons' interaction with a physically realistic plasma density profile. We carry out relativistic simulations of a massive Proca field evolving on a Kerr background, with modifications to account for the spatially varying effective mass. This allows us to treat the spin-1 case directly relevant to photons, and to study the effect of thinner disk profiles in the plasma. We find similar qualitative results to the scalar case, and so support the conclusions of that work: either a constant asymptotic mass or a shell-like plasma structure is required for superradiant growth to occur. We study thin disks and find a leakage of the superradiant cloud that suppresses its growth, concluding that thick disks are more likely to support the instability.

preprint2022arXiv

The Effect of Mission Duration on LISA Science Objectives

The science objectives of the LISA mission have been defined under the implicit assumption of a 4 yr continuous data stream. Based on the performance of LISA Pathfinder, it is now expected that LISA will have a duty cycle of $\approx 0.75$, which would reduce the effective span of usable data to 3 yr. This paper reports the results of a study by the LISA Science Group, which was charged with assessing the additional science return of increasing the mission lifetime. We explore various observational scenarios to assess the impact of mission duration on the main science objectives of the mission. We find that the science investigations most affected by mission duration concern the search for seed black holes at cosmic dawn, as well as the study of stellar-origin black holes and of their formation channels via multi-band and multi-messenger observations. We conclude that an extension to 6 yr of mission operations is recommended.

preprint2021arXiv

Constraining the primordial black hole scenario with Bayesian inference and machine learning: the GWTC-2 gravitational wave catalog

Primordial black holes (PBHs) might be formed in the early Universe and could comprise at least a fraction of the dark matter. Using the recently released GWTC-2 dataset from the third observing run of the LIGO-Virgo Collaboration, we investigate whether current observations are compatible with the hypothesis that all black hole mergers detected so far are of primordial origin. We constrain PBH formation models within a hierarchical Bayesian inference framework based on deep learning techniques, finding best-fit values for distinctive features of these models, including the PBH initial mass function, the fraction of PBHs in dark matter, and the accretion efficiency. The presence of several spinning binaries in the GWTC-2 dataset favors a scenario in which PBHs accrete and spin up. Our results indicate that PBHs may comprise only a fraction smaller than $0.3 \%$ of the total dark matter, and that the predicted PBH abundance is still compatible with other constraints.

preprint2021arXiv

Spin-induced black hole scalarization in Einstein-scalar-Gauss-Bonnet theory

We construct black hole solutions with spin-induced scalarization in a class of models where a scalar field is quadratically coupled to the topological Gauss-Bonnet term. Starting from the tachyonically unstable Kerr solutions, we obtain families of scalarized black holes such that the scalar field has either even or odd parity, and we investigate their domain of existence. The scalarized black holes can violate the Kerr rotation bound. We identify "critical" families of scalarized black hole solutions such that the expansion of the metric functions and of the scalar field at the horizon no longer allows for real coefficients. For the quadratic coupling considered here, solutions with spin-induced scalarization are entropically favored over Kerr solutions with the same mass and angular momentum.

preprint2020arXiv

$d+1$ formalism in Einstein-scalar-Gauss-Bonnet gravity

We present the $d+1$ formulation of Einstein-scalar-Gauss-Bonnet (ESGB) theories in dimension $D=d+1$ and for arbitrary (spacelike or timelike) slicings. We first build an action which generalizes those of Gibbons-Hawking-York and Myers to ESGB theories, showing that they can be described by a Dirichlet variational principle. We then generalize the Arnowitt-Deser-Misner (ADM) Lagrangian and Hamiltonian to ESGB theories, as well as the resulting $d+1$ decomposition of the equations of motion. Unlike general relativity, the canonical momenta of ESGB theories are nonlinear in the extrinsic curvature. This has two main implications: (i) the ADM Hamiltonian is generically multivalued, and the associated Hamiltonian evolution is not predictable; (ii) the "$d+1$" equations of motion are quasilinear, and they may break down in strongly curved, highly dynamical regimes. Our results should be useful to guide future developments of numerical relativity for ESGB gravity in the nonperturbative regime.

preprint2020arXiv

Astrophysical implications of GW190412 as a remnant of a previous black-hole merger

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 reported 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 that 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 interpretation is equally supported by the data. If the heavier component of GW190412 is indeed a merger remnant, then its spin magnitude is $χ_1=0.56_{-0.21}^{+0.19}$, which is higher than the value previously reported by the LIGO/Virgo collaboration.

preprint2020arXiv

Black holes in the low mass gap: Implications for gravitational wave observations

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 $χ_{\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.

preprint2020arXiv

Constraining the fraction of binary black holes formed in isolation and young star clusters with gravitational-wave data

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.

preprint2020arXiv

Distinguishing double neutron star from neutron star-black hole binary populations with gravitational wave observations

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.

preprint2020arXiv

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA.

preprint2020arXiv

Gravitational-wave signal-to-noise interpolation via neural networks

Computing signal-to-noise ratios (SNRs) is one of the most common tasks in gravitational-wave data analysis. While a single SNR evaluation is generally fast, computing SNRs for an entire population of merger events could be time consuming. We compute SNRs for aligned-spin binary black-hole mergers as a function of the (detector-frame) total mass, mass ratio and spin magnitudes using selected waveform models and detector noise curves, then we interpolate the SNRs in this four-dimensional parameter space with a simple neural network (a multilayer perceptron). The trained network can evaluate $10^6$ SNRs on a 4-core CPU within a minute with a median fractional error below $10^{-3}$. This corresponds to average speed-ups by factors in the range $[120,\,7.5\times10^4]$, depending on the underlying waveform model. Our trained network (and source code) is publicly available at https://github.com/kazewong/NeuralSNR, and it can be easily adapted to similar multidimensional interpolation problems.

preprint2020arXiv

LISA parameter estimation and source localization with higher harmonics of the ringdown

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.

preprint2020arXiv

Parametrized ringdown spin expansion coefficients: a data-analysis framework for black-hole spectroscopy with multiple events

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.

preprint2020arXiv

Probing the existence of ultralight bosons with a single gravitational-wave measurement

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.

preprint2020arXiv

Prospects for Fundamental Physics with LISA

In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.

preprint2020arXiv

The mass gap, the spin gap, and the origin of merging binary black holes

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.

preprint2019arXiv

Spinning and excited black holes in Einstein-scalar-Gauss-Bonnet theory

We construct rotating black holes in Einstein-scalar-Gauss-Bonnet theory with a quadratic coupling function. We map the domain of existence of the rotating fundamental solutions, we construct radially excited rotating black holes (including their existence lines), and we show that there are angularly excited rotating black holes. The bifurcation points of the radially and angularly excited solutions branching out of the Schwarzschild solution follow a regular pattern.

preprint2016arXiv

Astrophysical applications of the post-Tolman-Oppenheimer-Volkoff formalism

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 $μ_1$ and $χ$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

preprint2016arXiv

Compact Binary Merger Rates: Comparison with LIGO/Virgo Upper Limits

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.

preprint2016arXiv

Compact objects in Horndeski gravity

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.

preprint2016arXiv

Distinguishing black-hole spin-orbit resonances by their gravitational wave signatures. II: Full parameter estimation

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" $ξ$ and resonant family $ΔΦ=0,π$ (which uniquely label the source), the inclination $θ_{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.

preprint2016arXiv

eLISA eccentricity measurements as tracers of binary black hole formation

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.

preprint2016arXiv

Low-mass neutron stars: universal relations, the nuclear symmetry energy and gravitational radiation

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 $η\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 $η$ and of the central density. These fits may be used to constrain $η$, either via observations of binary pulsars in the electromagnetic spectrum, or via near-future observations of inspiralling compact binaries in the gravitational-wave spectrum.

preprint2016arXiv

Neutron stars in Horndeski gravity

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 $Λ$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.

preprint2016arXiv

Numerical Relativity and High Energy Physics: Recent Developments

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.

preprint2016arXiv

Science with the space-based interferometer eLISA. I: Supermassive black hole binaries

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.

preprint2016arXiv

Spectroscopy of Kerr black holes with Earth- and space-based interferometers

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.

preprint2015arXiv

Aligned spin neutron star-black hole mergers: a gravitational waveform amplitude model

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.

preprint2015arXiv

Effective potentials and morphological transitions for binary black-hole spin precession

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.

preprint2015arXiv

Gravitational-wave cutoff frequencies of tidally disruptive neutron star-black hole binary mergers

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 $χ$. 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 $χ$. 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.

preprint2015arXiv

Gravity-dominated unequal-mass black hole collisions

We continue our series of studies of high-energy collisions of black holes investigating unequal-mass, boosted head-on collisions in four dimensions. We show that the fraction of the center-of-mass energy radiated as gravitational waves becomes independent of mass ratio and approximately equal to $13\%$ at large energies. We support this conclusion with calculations using black hole perturbation theory and Smarr's zero-frequency limit approximation. These results lend strong support to the conjecture that the detailed structure of the colliding objects is irrelevant at high energies.

preprint2015arXiv

Massive Black Hole Science with eLISA

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.

preprint2015arXiv

Multi-timescale analysis of phase transitions in precessing black-hole binaries

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.

preprint2015arXiv

Post-Tolman-Oppenheimer-Volkoff formalism for relativistic stars

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.

preprint2015arXiv

Precessional instability in binary black holes with aligned spins

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{χ_1} \pm \sqrt{q χ_2})^4 (1-q)^{-2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is the mass ratio, and $χ_1$ ($χ_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.

preprint2015arXiv

Reply to "On the branching of quasinormal resonances of near-extremal Kerr black holes" by Shahar Hod

In a study of the quasinormal mode frequencies of nearly extremal black holes, we pointed out a bifurcation of the mode spectrum into modes with finite decay and modes with vanishing decay in the extremal limit. We provided analytic and semi-analytic results identifying which families of modes bifurcated, and when modes with finite decay rates exist when approaching the extremal limit. In a recent note (arXiv:1510.05604), Hod suggests that additional modes asymptote to finite decay at extremely high spin parameter, based on past work by Detweiler. We search for these suggested modes and find no evidence of their existence. In addition, we point out an inconsistency in the derivation of the proposed modes, which further indicates that these damped modes do not exist.

preprint2015arXiv

Slowly Rotating Anisotropic Neutron Stars in General Relativity and Scalar-Tensor Theory

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.

preprint2015arXiv

Slowly rotating black hole solutions in Horndeski gravity

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.

preprint2015arXiv

Tensor-multi-scalar theories: relativistic stars and 3+1 decomposition

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.

preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

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.

preprint2015arXiv

Ultra-high-energy debris from the collisional Penrose process

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 $η$ between the energy of post-collision particles detected at infinity and the energy of the colliding particles should be modest ($η\lesssim 1.5$). Schnittman has shown that these studies underestimated the maximum efficiency by about one order of magnitude (i.e., $η\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.

preprint2014arXiv

A Black-Hole Primer: Particles, Waves, Critical Phenomena and Superradiant Instabilities

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.

preprint2014arXiv

Distinguishing black-hole spin-orbit resonances by their gravitational-wave signatures

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 $ρ\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.

preprint2014arXiv

I-Love-Q, Spontaneously: Slowly Rotating Neutron Stars in Scalar-Tensor Theories

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.

preprint2014arXiv

Mixing of spherical and spheroidal modes in perturbed Kerr black holes

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 $ω_{\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.

preprint2014arXiv

Post-Circular Expansion of Eccentric Binary Inspirals: Fourier-Domain Waveforms in the Stationary Phase Approximation

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.

preprint2014arXiv

Torsional oscillations of neutron stars in scalar-tensor theory of gravity

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.

preprint2013arXiv

Astrophysical black holes as natural laboratories for fundamental physics and strong-field gravity

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.

preprint2013arXiv

Double Compact Objects II: Cosmological Merger Rates

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.

preprint2013arXiv

Gravito-Electromagnetic Perturbations of Kerr-Newman Black Holes: Stability and Isospectrality in the Slow-Rotation Limit

The most general stationary black-hole solution of Einstein-Maxwell theory in vacuum is the Kerr-Newman metric, specified by three parameters: mass M, spin J and charge Q. Within classical general relativity, the most important and challenging open problem in black-hole perturbation theory is the study of gravitational and electromagnetic fields in the Kerr-Newman geometry, because of the indissoluble coupling of the perturbation functions. Here we circumvent this long-standing problem by working in the slow-rotation limit. We compute the quasinormal modes up to linear order in J for any value of Q and provide the first, fully-consistent stability analysis of the Kerr-Newman metric. For scalar perturbations the quasinormal modes can be computed exactly, and we demonstrate that the method is accurate within 3% for spins J/Jmax<~0.5, where Jmax is the maximum allowed spin for any value of Q. Quite remarkably, we find numerical evidence that the axial and polar sectors of the gravito-electromagnetic perturbations are isospectral to linear order in the spin. The extension of our results to nonasymptotically flat space-times could be useful in the context of gauge/gravity dualities and string theory.

preprint2013arXiv

Nonspinning black hole-neutron star mergers: a model for the amplitude of gravitational waveforms

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.

preprint2013arXiv

Numerical simulations of single and binary black holes in scalar-tensor theories: circumventing the no-hair theorem

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.

preprint2013arXiv

Quasinormal modes of nearly extremal Kerr spacetimes: spectrum bifurcation and power-law ringdown

We provide an in-depth investigation of quasinormal-mode oscillations of Kerr black holes with nearly extremal angular momenta. We first discuss in greater detail the two distinct types of quasinormal mode frequencies presented in a recent paper (arXiv:1212.3271). One set of modes, that we call "zero-damping modes", has vanishing imaginary part in the extremal limit, and exists for all corotating perturbations (i.e., modes with azimuthal index m being nonnegative). The other set (the "damped modes") retains a finite decay rate even for extremal Kerr black holes, and exists only for a subset of corotating modes. As the angular momentum approaches its extremal value, the frequency spectrum bifurcates into these two distinct branches when both types of modes are present. We discuss the physical reason for the mode branching by developing and using a bound-state formulation for the perturbations of generic Kerr black holes. We also numerically explore the specific case of the fundamental l=2 modes, which have the greatest astrophysical interest. Using the results of these investigations, we compute the quasinormal mode response of a nearly extremal Kerr black hole to perturbations. We show that many superimposed overtones result in a slow power-law decay of the quasinormal ringing at early times, which later gives way to exponential decay. This exceptional early-time power-law decay implies that the ringdown phase is long-lived for black holes with large angular momentum, which could provide a promising strong source for gravitational-wave detectors.

preprint2013arXiv

Quasinormal ringing of Kerr black holes. II. Excitation by particles falling radially with arbitrary energy

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.

preprint2013arXiv

Relations Between Neutron-Star Parameters in the Hartle-Thorne Approximation

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.

preprint2013arXiv

Resonant-plane locking and spin alignment in stellar-mass black-hole binaries: a diagnostic of compact-binary formation

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 ΔΦ=0° or ΔΦ=+/-180°, where ΔΦis the angle between the components of the black-hole spins in the plane orthogonal to the orbital angular momentum. If the angles ΔΦ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.

preprint2013arXiv

Scalar, Electromagnetic and Gravitational Perturbations of Kerr-Newman Black Holes in the Slow-Rotation Limit

In Einstein-Maxwell theory, according to classic uniqueness theorems, the most general stationary black-hole solution is the axisymmetric Kerr-Newman metric, which is defined by three parameters: mass, spin and electric charge. The radial and angular dependence of gravitational and electromagnetic perturbations in the Kerr-Newman geometry do not seem to be separable. In this paper we circumvent this problem by studying scalar, electromagnetic and gravitational perturbations of Kerr-Newman black holes in the slow-rotation limit. We extend (and provide details of) the analysis presented in a recent Letter [arXiv:1304.1160]. Working at linear order in the spin, we present the first detailed derivation of the axial and polar perturbation equations in the gravito-electromagnetic case, and we compute the corresponding quasinormal modes for any value of the electric charge. Our study is the first self-consistent stability analysis of the Kerr-Newman metric, and in principle it can be extended to any order in the small rotation parameter. We find numerical evidence that the axial and polar sectors are isospectral at first order in the spin, and speculate on the possible implications of this result.

preprint2013arXiv

The Transient Gravitational-Wave Sky

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.

preprint2013arXiv

Universality, maximum radiation and absorption in high-energy collisions of black holes with spin

We explore the impact of black hole spins on the dynamics of high-energy black hole collisions. We report results from numerical simulations with gamma-factors up to 2.49 and dimensionless spin parameter +0.85, +0.6, 0, -0.6, -0.85. We find that the scattering threshold becomes independent of spin at large center-of-mass energies, confirming previous conjectures that structure does not matter in ultrarelativistic collisions. It has further been argued that in this limit all of the kinetic energy of the system may be radiated by fine tuning the impact parameter to threshold. On the contrary, we find that only about 60% of the kinetic energy is radiated for gamma=2.49. By monitoring apparent horizons before and after scattering events we show that the "missing energy" is absorbed by the individual black holes in the encounter, and moreover the individual black-hole spins change significantly. We support this conclusion with perturbative calculations. An extrapolation of our results to the limit gamma->infinity suggests that about half of the center-of-mass energy of the system can be emitted in gravitational radiation, while the rest must be converted into rest-mass and spin energy.

preprint2012arXiv

Black-Hole Bombs and Photon-Mass Bounds

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.

preprint2012arXiv

Branching of quasinormal modes for nearly extremal Kerr black holes

We show that nearly extremal Kerr black holes have two distinct sets of quasinormal modes, which we call zero-damping modes (ZDMs) and damped modes (DMs). The ZDMs exist for all harmonic indices $l$ and $m \ge 0$, and their frequencies cluster onto the real axis in the extremal limit. The DMs have nonzero damping for all black hole spins; they exist for all counterrotating modes ($m<0$) and for corotating modes with $0\leq μ\lesssim μ_c=0.74$ (in the eikonal limit), where $μ\equiv m/(l+1/2)$. When the two families coexist, ZDMs and DMs merge to form a single set of quasinormal modes as the black hole spin decreases. Using the effective potential for perturbations of the Kerr spacetime, we give intuitive explanations for the absence of DMs in certain areas of the spectrum and for the branching of the spectrum into ZDMs and DMs at large spins.

preprint2012arXiv

Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates

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.

preprint2012arXiv

Effects of post-Newtonian Spin Alignment on the Distribution of Black-Hole Recoils

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.

preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

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.

preprint2012arXiv

Gravitational radiation from compact binary systems in the massive Brans-Dicke theory of gravity

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 ω_{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 ω_{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 ω_{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 ω_{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.

preprint2012arXiv

Light scalar field constraints from gravitational-wave observations of compact binaries

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 ω_{BD}. We show that gravitational-wave observations of nonspinning neutron star-black hole binary inspirals can be used to set lower bounds on ω_{BD} and upper bounds on the combination m_s/\sqrt{ω_{BD}}$. We estimate via a Fisher matrix analysis that individual observations with signal-to-noise ratio ρwould yield (m_s/\sqrt{ω_{BD}})(ρ/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.

preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

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.

preprint2012arXiv

Perturbations of slowly rotating black holes: massive vector fields in the Kerr metric

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.

preprint2011arXiv

Accuracy of the Post-Newtonian Approximation: Optimal Asymptotic Expansion for Quasi-Circular, Extreme-Mass Ratio Inspirals

We study the accuracy of the post-Newtonian (PN) approximation and its formal region of validity, by investigating its optimal asymptotic expansion for the quasi-circular, adiabatic inspiral of a point particle into a Schwarzschild black hole. By comparing the PN expansion of the energy flux to numerical calculations in the perturbative Teukolsky formalism, we show that (i) the inclusion of higher multipoles is necessary to establish the accuracy of high-order PN terms, and (ii) the region of validity of PN theory is largest at relative O(1/c^6) (3PN order). The latter result suggests that the series diverges beyond 3PN order, at least in the extreme mass-ratio limit, probably due to the appearance of logarithmic terms in the energy flux. The study presented here is a first formal attempt to determine the region of validity of the PN approximation using asymptotic analysis. Therefore, it should serve as a template to perform similar studies on other systems, such as comparable-mass quasi-circular inspirals computed by high-accuracy numerical relativistic simulations.

preprint2011arXiv

Accuracy of the post-Newtonian approximation. II. Optimal asymptotic expansion of the energy flux for quasicircular, extreme mass-ratio inspirals into a Kerr black hole

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.

preprint2011arXiv

Compact stars in alternative theories of gravity. Einstein-Dilaton-Gauss-Bonnet gravity

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.

preprint2011arXiv

Floating and sinking: the imprint of massive scalars around rotating black holes

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.

preprint2011arXiv

Graviton mass bounds from space-based gravitational-wave observations of massive black hole populations

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 λ_g=h/(m_g c) by observing multiple inspiralling black hole binaries. We show that while observations of individual inspirals will yield mean bounds λ_g~3x10^15 km, the combined bound from observing ~50 events in a two-year mission is about ten times better: λ_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 λ_g applies to any theory whose deviations from general relativity are parametrized by a set of global parameters.

preprint2011arXiv

The vacuum revealed: the final state of vacuum instabilities in compact stars

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.

preprint2010arXiv

Constraining properties of the black hole population using LISA

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.

preprint2010arXiv

Final spins from the merger of precessing binary 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.

preprint2010arXiv

Gravitational instabilities of superspinars

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.

preprint2010arXiv

Gravitational-wave signatures of the absence of an event horizon. II. Extreme mass ratio inspirals in the spacetime of a thin-shell gravastar

We study gravitational-wave emission from the quasi-circular, extreme mass ratio inspiral of compact objects of mass m0 into massive objects of mass M>>m0 whose external metric is identical to the Schwarzschild metric, except for the absence of an event horizon. To be specific we consider one of the simplest realizations of such an object: a nonrotating thin-shell gravastar. The power radiated in gravitational waves during the inspiral shows distinctive peaks corresponding to the excitation of the polar oscillation modes of the gravastar. For ultra-compact gravastars the frequency of these peaks depends mildly on the gravastar compactness. For masses M~10^6Msun the peaks typically lie within the optimal sensitivity bandwidth of LISA, potentially providing a unique signature of the horizonless nature of the central object. For relatively modest values of the gravastar compactness the radiated power has even more peculiar features, carrying the signature of the microscopic properties of the physical surface replacing the event horizon.

preprint2010arXiv

Reconstructing the massive black hole cosmic history through gravitational waves

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.

preprint2010arXiv

Relativistic Suppression of Black Hole Recoils

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.

preprint2010arXiv

Semianalytical estimates of scattering thresholds and gravitational radiation in ultrarelativistic black hole encounters

Ultrarelativistic collisions of black holes are ideal gedanken experiments to study the nonlinearities of general relativity. In this paper we use semianalytical tools to better understand the nature of these collisions and the emitted gravitational radiation. We explain many features of the energy spectra extracted from numerical relativity simulations using two complementary semianalytical calculations. In the first calculation we estimate the radiation by a "zero-frequency limit" analysis of the collision of two point particles with finite impact parameter. In the second calculation we replace one of the black holes by a point particle plunging with arbitrary energy and impact parameter into a Schwarzschild black hole, and we explore the multipolar structure of the radiation paying particular attention to the near-critical regime. We also use a geodesic analogy to provide qualitative estimates of the dependence of the scattering threshold on the black hole spin and on the dimensionality of the spacetime.

preprint2010arXiv

Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called "superkicks". In astrophysical, quasicircular inspirals, such kicks can be as large as ~3,000 km/s; here, we find configurations that exceed ~15,000 km/s. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the "bobbing" motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.

preprint2010arXiv

Up to eleven: radiation from particles with arbitrary energy falling into higher-dimensional black holes

We consider point particles with arbitrary energy per unit mass E that fall radially into a higher-dimensional, nonrotating, asymptotically flat black hole. We compute the energy and linear momentum radiated in this process as functions of E and of the spacetime dimensionality D=n+2 for n=2,...,9 (in some cases we go up to 11). We find that the total energy radiated increases with n for particles falling from rest (E=1). For fixed particle energies 1<E<=2 we show explicitly that the radiation has a local minimum at some critical value of n, and then it increases with n. We conjecture that such a minimum exists also for higher particle energies. The present point-particle calculation breaks down when n=11, because then the radiated energy becomes larger than the particle mass. Quite interestingly, for n=11 the radiated energy predicted by our calculation would also violate Hawking's area bound. This hints at a qualitative change in gravitational radiation emission for n>11. Our results are in very good agreement with numerical simulations of low-energy, unequal-mass black hole collisions in D=5 (that will be reported elsewhere) and they are a useful benchmark for future nonlinear evolutions of the higher-dimensional Einstein equations.

preprint2009arXiv

Geodesic stability, Lyapunov exponents and quasinormal modes

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.

preprint2009arXiv

Gravitational wave signatures of the absence of an event horizon. I. Nonradial oscillations of a thin-shell gravastar

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.

preprint2009arXiv

Internal dynamics and dynamical friction effects in the dwarf spheroidal galaxy in Fornax

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.

preprint2009arXiv

Massive Black Hole Binary Inspirals: Results from the LISA Parameter Estimation Taskforce

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.

preprint2009arXiv

Quasinormal modes of black holes and black branes

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.

preprint2008arXiv

Cosmological black hole spin evolution by mergers and accretion

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.

preprint2006arXiv

Eigenvalues and eigenfunctions of spin-weighted spheroidal harmonics in four and higher dimensions

Spin-weighted spheroidal harmonics are useful in a variety of physical situations, including light scattering, nuclear modeling, signal processing, electromagnetic wave propagation, black hole perturbation theory in four and higher dimensions, quantum field theory in curved space-time and studies of D-branes. We first review analytic and numerical calculations of their eigenvalues and eigenfunctions in four dimensions, filling gaps in the existing literature when necessary. Then we compute the angular dependence of the spin-weighted spheroidal harmonics corresponding to slowly-damped quasinormal mode frequencies of the Kerr black hole, providing numerical tables and approximate formulas for their scalar products. Finally we present an exhaustive analytic and numerical study of scalar spheroidal harmonics in (n+4) dimensions.

preprint2005arXiv

Rotating neutron stars: an invariant comparison of approximate and numerical spacetime models

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.

preprint2005arXiv

What we (don't) know about black hole formation in high-energy collisions

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.

preprint2004arXiv

Quasinormal modes and classical wave propagation in analogue black holes

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.

preprint2003arXiv

Stability of five-dimensional rotating black holes projected on the brane

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.