Quasinormal modes of Kerr-Newman black holes: coupling of electromagnetic and gravitational perturbations

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

Explaining LIGO's observations via isolated binary evolution with natal kicks

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

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.Comment: 8 pages, 2 figures. Updated to match the version published on PRD. Corrected a small typo (which does not affect the results) in equation (6) of the published pape

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.Comment: 20 pages, 16 figures, new PN terms, submitted to PR

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.Comment: 15 pages, 7 figures, 1 tabl

Atmospheric circulation patterns, cloud-to-ground lightning, and locally intense convective rainfall associated with debris flow initiation in the Dolomite Alps of northeastern Italy

The Dolomite Alps of northeastern Italy experience debris flows with great frequency during the summer months. An ample supply of unconsolidated material on steep slopes and a summer season climate regime characterized by recurrent thunderstorms combine to produce an abundance of these destructive hydro-geologic events. In the past, debris flow events have been studied primarily in the context of their geologic and geomorphic characteristics. The atmospheric contribution to these mass-wasting events has been limited to recording rainfall and developing intensity thresholds for debris mobilization. This study aims to expand the examination of atmospheric processes that preceded both locally intense convective rainfall (LICR) and debris flows in the Dolomite region. 500 hPa pressure level plots of geopotential heights were constructed for a period of 3 days prior to debris flow events to gain insight into the synoptic-scale processes which provide an environment conducive to LICR in the Dolomites. Cloud-to-ground (CG) lightning flash data recorded at the meso-scale were incorporated to assess the convective environment proximal to debris flow source regions. Twelve events were analyzed and from this analysis three common synoptic-scale circulation patterns were identified. Evaluation of CG flashes at smaller spatial and temporal scales illustrated that convective processes vary in their production of CF flashes (total number) and the spatial distribution of flashes can also be quite different between events over longer periods. During the 60 min interval immediately preceding debris flow a majority of cases exhibited spatial and temporal colocation of LICR and CG flashes. Also a number of CG flash parameters were found to be significantly correlated to rainfall intensity prior to debris flow initiation

On gravitational-wave spectroscopy of massive black holes with the space interferometer LISA

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

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.Comment: 12 pages, 3 figures, submitted to Class. Quantum Grav. for proceedings of 8th LISA Symposium; v2 minor changes for consistency with accepted versio

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.Comment: 24 pages, 16 figures, submitted to Phys. Rev.

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

This is the author accepted manuscript. The final version is available from APS via http://dx.doi.org/10.1103/PhysRevD.92.064016The 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.D.G. is supported by the UK STFC and the Isaac Newton Studentship of the University of Cambridge. M.K. is supported by Alfred P. Sloan Foundation grant FG-2015-65299. R.O'S. is supported by NSF grants PHY-0970074 and PHY-1307429. E.B. is sup- ported by NSF CAREER Grant PHY-1055103 and by FCT contract IF/00797/2014/CP1214/CT0012 under the IF2014 Programme. U.S. is supported by FP7- PEOPLE-2011-CIG Grant No. 293412, FP7-PEOPLE- 2011-IRSES Grant No.295189, H2020 ERC Consolida- tor Grant Agreement No. MaGRaTh-646597, SDSC and TACC through XSEDE Grant No. PHY-090003 by the NSF, Finis Terrae through Grant No. ICTS- CESGA-249, STFC Roller Grant No. ST/L000636/1 and DiRAC's Cosmos Shared Memory system through BIS Grant No. ST/J005673/1 and STFC Grant Nos. ST/H008586/1, ST/K00333X/1