14332015
Binary black hole merger gravitational waves and recoil in the large mass ratio limit
Spectacular breakthroughs in numerical relativity now make it possible to compute spacetime dynamics in almost complete generality, allowing us to model the coalescence and merger of binary black holes with essentially no approximations. The primary limitation of these calculations is now computational. In particular, it is difficult to model systems with large mass ratio and large spins, since one must accurately resolve the multiple length scales that play a role in such systems. Perturbation theory can play an important role in extending the reach of computational modeling for binary systems. In this paper, we present first results of a code that allows us to model the gravitational waves generated by the inspiral, merger, and ringdown of a binary system in which one member of the binary is much more massive than the other. This allows us to accurately calibrate binary dynamics in the large mass ratio regime. We focus in this analysis on the recoil imparted to the merged remnant by these waves. We closely examine the ''antikick,'' an antiphase cancellation of the recoil arising from the plunge and ringdown waves, described in detail by Schnittman et al. We find that, for orbits aligned with the black hole spin, themore » antikick grows as a function of spin. The total recoil is smallest for prograde coalescence into a rapidly rotating black hole, and largest for retrograde coalescence. Amusingly, this completely reverses the predicted trend for kick versus spin from analyses that only include inspiral information.« less
119020610
Towards adiabatic waveforms for inspiral into Kerr black holes : A new model of the source for the time domain perturbation equation
We revisit the problem of the emission of gravitational waves from a test mass orbiting and thus perturbing a Kerr black hole. The source term of the Teukolsky perturbation equation contains a Dirac delta function which represents a point particle. We present a technique to effectively model the delta function and its derivatives using as few as four points on a numerical grid. The source term is then incorporated into a code that evolves the Teukolsky equation in the time domain as a (2+1) dimensional partial differential equation. The waveforms and energy fluxes are extracted far from the black hole. Our comparisons with earlier work show an order of magnitude gain in performance (speed) and numerical errors less than 1% for a large fraction of parameter space. As a first application of this code, we analyze the effect of finite extraction radius on the energy fluxes. This paper is the first in a series whose goal is to develop adiabatic waveforms describing the inspiral of a small compact body into a massive Kerr black hole.
8108905
Untangling the merger history of massive black holes with LISA
Binary black hole coalescences emit gravitational waves that will be measurable by the space-based detector LISA to large redshifts. This suggests that LISA may be able to observe black holes grow and evolve as the Universe evolves, mapping the distribution of black hole masses as a function of redshift. An immediate difficulty with this idea is that LISA measures certain redshifted combinations of masses with good accuracy: if a system has some mass parameter m, then LISA measures (1 + z)m. This mass-redshift degeneracy makes it difficult to follow the mass evolution. In many cases, LISA will also measure the luminosity distance D of a coalescence accurately. Since cosmological parameters (particularly the mean density, the cosmological constant and the Hubble constant) are now known with moderate precision, we can obtain z from D and break the degeneracy. This makes it possible to untangle the mass and redshift and to study the mass and merger history of black holes. Mapping the black hole mass distribution could open a window on to an early epoch of structure formation.
36773341
Circular and non-circular nearly horizon-skimming orbits in Kerr spacetimes
We have performed a detailed analysis of orbital motion in the vicinity of a nearly extremal Kerr black hole. For very rapidly rotating black holes—spin parameter aJ=M > 0:9524M —we have found a class of very strong-field eccentric orbits whose orbital angular momentum Lz increases with the orbit's inclination with respect to the equatorial plane, while keeping latus rectum and eccentricity fixed. This behavior is in contrast with Newtonian intuition, and is in fact opposite to the normal behavior of black hole orbits. Such behavior was noted previously for circular orbits; since it only applies to orbits very close to the black hole, they were named ''nearly horizon-skimming orbits.'' Our current analysis generalizes this result, mapping out the full generic (inclined and eccentric) family of nearly horizon- skimming orbits. The earlier work on circular orbits reported that, under gravitational radiation emission, nearly horizon-skimming orbits exhibit unusual inspiral, tending to evolve to smaller orbit inclination, toward prograde equatorial configuration. Normal orbits, by contrast, always demonstrate slowly growing orbit inclination —orbits evolve toward the retrograde equatorial configuration. Using up-to-date Teukolsky-based fluxes, we have concluded that the earlier result was incorrect —all circular orbits, including nearly horizon-skimming ones, exhibit growing orbit inclination under radiative backreaction. Using kludge fluxes based on a Post-Newtonian expansion corrected with fits to circular and to equatorial Teukolsky-based fluxes, we argue that the inclination grows also for eccentric nearly horizon-skimming orbits. We also find that the inclination change is, in any case, very small. As such, we conclude that these orbits are not likely to have a clear and peculiar imprint on the gravitational waveforms expected to be measured by the space-based detector LISA.
16126239
Seismic gravity-gradient noise in interferometric gravitational-wave detectors
When ambient seismic waves pass near and under an interferometric gravitational-wave detector, they induce density perturbations in the earth, which in turn produce fluctuating gravitational forces on the interferometer’s test masses. These forces mimic a stochastic background of gravitational waves and thus constitute a noise source. This seismic gravitygradient noise has been estimated and discussed previously by Saulson using a simple model of the earth’s ambient seismic motions. In this paper, we develop a more sophisticated model of these motions, based on the theory of multimode Rayleigh and Love waves propagating in a multilayer medium that approximates the geological strata at the LIGO sites (Tables II–IV), and we use this model to revisit seismic gravity gradients. We characterize the seismic gravity-gradient noise by a transfer function, T (f) ≡ x̃(f)/W̃ (f), from the spectrum of rms seismic displacements averaged over vertical and horizontal directions, W̃ (f), to the spectrum of interferometric test-mass motions, x̃(f) ≡ Lh̃(f); here L is the interferometer arm length, h̃(f) is the gravitational-wave noise spectrum, and f is frequency. Our model predicts a transfer function with essentially the same functional form as that derived by Saulson, T ' 4πGρ(2πf)β(f), where ρ is the density of the earth near the test masses, G is Newton’s constant, and β(f) ≡ γ(f)Γ(f)β(f) is a dimensionless reduced transfer function whose components γ ' 1 and Γ ' 1 account for a weak correlation between the interferometer’s two test masses (Fig. 1) and a slight reduction of the noise due to the height of the test masses above the earth’s surface. This paper’s primary foci are (i) a study of how β(f) ' β(f) depends on the various Rayleigh and Love modes that are present in the seismic spectrum (Figs. 4–11 and Table I), (ii) an attempt to estimate which modes are actually present at the two LIGO sites at quiet times and at noisy times, and (iii) a corresponding estimate of the magnitude of β(f) at quiet and noisy times. We conclude that at quiet times β ' 0.35–0.6 at the LIGO sites, and at noisy times β ' 0.15–1.4. (For comparison, Saulson’s simple model gave β = β = 1/ √ 3 = 0.58.) By folding our resulting transfer function into the “standard LIGO seismic spectrum” [Eq. (1.29)], which approximates W̃ (f) at typical times, we obtain the gravity-gradient noise spectra shown in Fig. 2. At quiet times this noise is below the benchmark noise level of “advanced LIGO interferometers” at all frequencies (though not by much at ∼ 10 Hz); at noisy times it may significantly exceed the advanced noise level near 10 Hz. The lower edge of our quiet-time noise constitutes a limit, beyond which there would be little gain from further improvements in vibration isolation and thermal noise, unless one can also reduce the seismic gravity gradient noise. Two methods of such reduction are briefly discussed: monitoring the earth’s density perturbations near each test mass, computing the gravitational forces they produce, and correcting the data for those forces; and constructing narrow moats around the interferometers’ corner and end stations to shield out the fundamental-mode Rayleigh waves, which we suspect dominate at quiet times.
1891241
How black holes get their kicks: Radiation recoil in binary black hole mergers
Gravitational waves from the coalescence of binary black holes carry linear momentum, causing center of mass recoil. This ``radiation rocket'' has important implications for systems with escape speeds of order the recoil velocity. We describe new recoil calculations using high precision black hole perturbation theory to estimate the magnitude of the recoil for the slow ``inspiral'' coalescence phase; coupled with a cruder calculation for the final ``plunge'', we estimate the total recoil imparted to a merged black hole. We find that velocities of many tens to a few hundred km/sec can be achieved fairly easily. The recoil probably never exceeds about 500 km/sec.
119858324
Anyons from bosons
Abstract We consider systems containing two or more distinct species of particles in two spatial dimensions. In quantizations of these systems, the statistics of composites containing more than one type of particle are not completely determined by the statistics of the constituents. In particular there exist quantum theories in which two bosons can combine to form an anyon with any desired statistical angle. We demonstrate these results using the topological approach to quantum kinematics, which in this case leads to a generalization of ordinary braid theory in two dimensions. Comparisons are made to three-dimensional systems.
121457285
INITIAL STEREO ANALYSIS OF MRK 421 FROM THE VERITAS TELESCOPES
VERITAS (Very Energetic Radiation Imaging Telescope Array System), an array of ground-based gamma-ray telescopes in southern Arizona, USA, has been taking data in hardware stereo mode since March, 2006. The April–May 2006 dark run provided a large set of data from two telescopes on the known blazar Markarian (Mrk) 421. An initial analysis produced a light curve and preliminary cuts showing the two telescope array's angular resolution to be 0.19°. The remaining two VERITAS telescopes will be brought online by January, 2007.
119377823
The Hubble Space Telescope Extragalactic Distance Scale Key Project. XXIII. The Discovery of Cepheids in NGC 3319
The distance to NGC 3319 has been determined from Cepheid variable stars as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen and four epochs of observations, using filters F555W (V) and F814W (I), respectively, were made with the Wide Field Planetary Camera 2. Thirty-three Cepheid variables between periods of 8 and 47 days were discovered. Adopting a Large Magellanic Cloud distance modulus of 18.50 ± 0.10 mag and extinction of E(V-I) = 0.13 mag, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.78 ± 0.14 (random) ± 0.10 (systematic) mag and the extinction of E(V-I) = 0.06 mag were determined for NGC 3319. This galaxy is the last galaxy observed for the HST H0 Key Project.
197531313
The HST Key Project on the Extragalactic Distance Scale. XXVIII. Combining the Constraints on the Hubble Constant
Since the launch of the Hubble Space Telescope nine years ago, Cepheid distances to 25 galaxies have been determined for the purpose of calibrating secondary distance indicators. A variety of these can now be calibrated, and the accompanying papers by Sakai, Kelson, Ferrarese, and Gibson employ the full set of 25 galaxies to consider the Tully-Fisher relation, the fundamental plane of elliptical galaxies, Type Ia supernovae, and surface brightness fluctuations. 
When calibrated with Cepheid distances, each of these methods yields a measurement of the Hubble constant and a corresponding measurement uncertainty. We combine these measurements in this paper, together with a model of the velocity field, to yield the best available estimate of the value of H_0 within the range of these secondary distance indicators and its uncertainty. 
The result is H_0 = 71 +/- 6 km/sec/Mpc. The largest contributor to the uncertainty of this 67% confidence level result is the distance of the Large Magellanic Cloud, which has been assumed to be 50 +/- 3 kpc.
6418107
A new search for primordial black hole evaporations using the Whipple gamma-ray telescope
Stephen Hawking's prediction that black holes should radiate like black bodies has several important consequences, including the possibility of the detection of small (~1015 g) black holes created in the very early universe. The detection of such primordial black holes (PBHs) would be an important discovery, not only confirming Hawking's theory, but also providing valuable insights into the history of the early universe. A search through 5.5 years of archival data from the Whipple Atmospheric Cerenkov Telescope is made for TeV gamma-ray bursts on 1, 3, and 5 s timescales. On the basis of a null result from this direct search for PBH evaporations, an upper limit of 1.08 × 106 pc−3 yr−1 (99% CL) is set on the PBH evaporation rate in the local region of the galaxy, assuming the Standard Model of particle physics. This is more than a factor of two better than the previous limit at this energy range and includes longer timescales than have previously been explored. Comparison of this result with previous limits on the fraction of the critical density comprised by PBHs, Ωpbh, depends strongly on assumptions made about PBH clustering; in models predicting strong PBH clustering, the limit in this work could be as many as ten orders of magnitude more stringently than those set by diffuse MeV gamma-ray observations.