22 research outputs found
Binary black holes on a budget: Simulations using workstations
Binary black hole simulations have traditionally been computationally very
expensive: current simulations are performed in supercomputers involving dozens
if not hundreds of processors, thus systematic studies of the parameter space
of binary black hole encounters still seem prohibitive with current technology.
Here we show how the multi-layered refinement level code BAM can be used on
dual processor workstations to simulate certain binary black hole systems. BAM,
based on the moving punctures method, provides grid structures composed of
boxes of increasing resolution near the center of the grid. In the case of
binaries, the highest resolution boxes are placed around each black hole and
they track them in their orbits until the final merger when a single set of
levels surrounds the black hole remnant. This is particularly useful when
simulating spinning black holes since the gravitational fields gradients are
larger. We present simulations of binaries with equal mass black holes with
spins parallel to the binary axis and intrinsic magnitude of S/m^2= 0.75. Our
results compare favorably to those of previous simulations of this particular
system. We show that the moving punctures method produces stable simulations at
maximum spatial resolutions up to M/160 and for durations of up to the
equivalent of 20 orbital periods.Comment: 20 pages, 8 figures. Final version, to appear in a special issue of
Class. Quantum Grav. based on the New Frontiers in Numerical Relativity
Conference, Golm, July 200
BSSN in Spherical Symmetry
The BSSN (Baumgarte-Shapiro-Shibata-Nakamura) formulation of the Einstein
evolution equations is written in spherical symmetry. These equations can be
used to address a number of technical and conceptual issues in numerical
relativity in the context of a single Schwarzschild black hole. One of the
benefits of spherical symmetry is that the numerical grid points can be tracked
on a Kruskal--Szekeres diagram. Boundary conditions suitable for puncture
evolution of a Schwarzschild black hole are presented. Several results are
shown for puncture evolution using a fourth--order finite difference
implementation of the equations.Comment: This is the final version to be published in CQG. It contains much
more information and detail than the original versio
Radiation content of Conformally flat initial data
We study the radiation of energy and linear momentum emitted to infinity by
the headon collision of binary black holes, starting from rest at a finite
initial separation, in the extreme mass ratio limit. For these configurations
we identify the radiation produced by the initially conformally flat choice of
the three geometry. This identification suggests that the radiated energy and
momentum of headon collisions will not be dominated by the details of the
initial data for evolution of holes from initial proper separations
. For non-headon orbits, where the amount of radiation is orders of
magnitude larger, the conformally flat initial data may provide a relative even
better approximation.Comment: 4 pages, 4 figure
Initial data for Einstein's equations with superposed gravitational waves
A method is presented to construct initial data for Einstein's equations as a
superposition of a gravitational wave perturbation on an arbitrary stationary
background spacetime. The method combines the conformal thin sandwich formalism
with linear gravitational waves, and allows detailed control over
characteristics of the superposed gravitational wave like shape, location and
propagation direction. It is furthermore fully covariant with respect to
spatial coordinate changes and allows for very large amplitude of the
gravitational wave.Comment: Version accepted by PRD; added convergence plots, expanded
discussion. 9 pages, 9 figure
Status of NINJA: the Numerical INJection Analysis project
The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise
Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration
The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a
joint effort between members of the numerical relativity, analytical relativity
and gravitational-wave data analysis communities. The goal of the NRAR
collaboration is to produce numerical-relativity simulations of compact
binaries and use them to develop accurate analytical templates for the
LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and
extracting astrophysical information from them. We describe the results of the
first stage of the NRAR project, which focused on producing an initial set of
numerical waveforms from binary black holes with moderate mass ratios and
spins, as well as one non-spinning binary configuration which has a mass ratio
of 10. All of the numerical waveforms are analysed in a uniform and consistent
manner, with numerical errors evaluated using an analysis code created by
members of the NRAR collaboration. We compare previously-calibrated,
non-precessing analytical waveforms, notably the effective-one-body (EOB) and
phenomenological template families, to the newly-produced numerical waveforms.
We find that when the binary's total mass is ~100-200 solar masses, current EOB
and phenomenological models of spinning, non-precessing binary waveforms have
overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary
numerical waveforms with mass ratios <= 4, when maximizing over binary
parameters. This implies that the loss of event rate due to modelling error is
below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to
five non-spinning waveforms with mass ratio smaller than 6 have overlaps above
99.7% with the numerical waveform with a mass ratio of 10, without even
maximizing on the binary parameters.Comment: 51 pages, 10 figures; published versio
The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries
The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations
The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries (Addendum to vol 29, 124001, 2012)
This is an Addendum for the article 2012 Class. Quantum Grav. 29 124001 A recent paper (Ajith et al 2012 Class. Quantum Grav. 29 124001) described a catalog of 56 hybrid post-Newtonian/numerical-relativity waveforms modeling the inspiral, merger and ringdown of binary black hole systems spanning a range of mass ratios and spins. This catalog has been created and validated for use in the NINJA-2 project to study the sensitivity of gravitational-wave search and parameter-estimation algorithms. The contents of this catalog are being made available for public use. This addendum describes the public release
Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project
The Numerical INJection Analysis (NINJA) project is a collaborative effort
between members of the numerical relativity and gravitational-wave data
analysis communities. The purpose of NINJA is to study the sensitivity of
existing gravitational-wave search algorithms using numerically generated
waveforms and to foster closer collaboration between the numerical relativity
and data analysis communities. We describe the results of the first NINJA
analysis which focused on gravitational waveforms from binary black hole
coalescence. Ten numerical relativity groups contributed numerical data which
were used to generate a set of gravitational-wave signals. These signals were
injected into a simulated data set, designed to mimic the response of the
Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this
data using search and parameter-estimation pipelines. Matched filter
algorithms, un-modelled-burst searches and Bayesian parameter-estimation and
model-selection algorithms were applied to the data. We report the efficiency
of these search methods in detecting the numerical waveforms and measuring
their parameters. We describe preliminary comparisons between the different
search methods and suggest improvements for future NINJA analyses.Comment: 56 pages, 25 figures; various clarifications; accepted to CQ
Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration
The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters