Search

The recent detection of gravitational waves and electromagnetic counterparts emitted during and after the collision of two neutron stars marks a breakthrough in the field of multi-messenger astronomy. Numerical relativity simulations are the only tool to describe the binary's merger dynamics in the regime when speeds are largest and gravity is strongest. In this work we report state-of-the-art binary neutron star simulations for irrotational (non-spinning) and spinning configurations. The main use of these simulations is to model the gravitational-wave signal. Key numerical requirements are the understanding of the convergence properties of the numerical data and a detailed error budget. The simulations have been performed on different HPC clusters, they use multiple grid resolutions, and are based on eccentricity reduced quasi-circular initial data. We obtain convergent waveforms with phase errors of 0.5-1.5 rad accumulated over approximately 12 orbits to merger. The waveforms have been used for the construction of a phenomenological waveform model which has been applied for the analysis of the recent binary neutron star detection. Additionally, we show that the data can also be used to test other state-of-the-art semi-analytical waveform models., Comment: published as proceeding for the 26th Euromicro International Conference on Parallel, Distributed, and Network-Based Processing. PDP 2018 IEEE Catalog Number: CFP18169

It is well-known that Bowen-York initial data contain spurious radiation. Although this ``junk'' radiation has been seen to be small for non-spinning black-hole binaries in circular orbit, its magnitude increases when the black holes are given spin. It is possible to reduce the spurious radiation by applying the puncture approach to multiple Kerr black holes, as we demonstrate for examples of head-on collisions of equal-mass black-hole binaries., Comment: 10 pages, 2 figures, submitted to special "New Frontiers in Numerical Relativity" issue of Classical and Quantum Gravity

We present a series of test beds for numerical codes designed to find apparent horizons. We consider three apparent horizon finders that use different numerical methods: one of them in axisymmetry, and two fully three-dimensional. We concentrate first on a toy model that has a simple horizon structure, and then go on to study single and multiple black hole data sets. We use our finders to look for apparent horizons in Brill wave initial data where we discover that some results published previously are not correct. For pure wave and multiple black hole spacetimes, we apply our finders to survey parameter space, mapping out properties of interesting data sets for future evolutions., Comment: 19 pages Revtex, 22 figures included

Numerical relativity is finally approaching a state where the evolution of rather general (3+1)-dimensional data sets can be computed in order to solve the Einstein equations. After a general introduction, three topics of current interest are briefly reviewed: binary black hole mergers, the evolution of strong gravitational waves, and shift conditions for neutron star binaries., Comment: 23 pages, 5 figures, review talk for Journees Relativistes 1999

We investigate the stationary end-state obtained by evolving a collapsing spherical star with the gauges routinely adopted to study puncture black holes. We compare the end-state of the collapse with the trumpet solution found in the evolution of a single wormhole slice and show that the two solutions closely agree. We demonstrate that the agreement is caused by the use of the Gamma-driver shift condition, which allows the matter to fall inwards into a region of spacetime that is not resolved by the numerical grid, and which simultaneously finds the stationary coordinates of the trumpet outside the matter.

Solutions to both the diffeomorphism and the hamiltonian constraint of quantum gravity have been found in the loop representation, which is based on Ashtekar's new variables. While the diffeomorphism constraint is easily solved by considering loop functionals which are knot invariants, there remains the puzzle why several of the known knot invariants are also solutions to the hamiltonian constraint. We show how the Jones polynomial gives rise to an infinite set of solutions to all the constraints of quantum gravity thereby illuminating the structure of the space of solutions and suggesting the existance of a deep connection between quantum gravity and knot theory at a dynamical level., 7pp

We quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases., Comment: 17 pages, 9 figures. Version accepted by PRD

We study zoom-whirl behaviour of equal mass, non-spinning black hole binaries in full general relativity. The magnitude of the linear momentum of the initial data is fixed to that of a quasi-circular orbit, and its direction is varied. We find a global maximum in radiated energy for a configuration which completes roughly one orbit. The radiated energy in this case exceeds the value of a quasi-circular binary with the same momentum by 15%. The direction parameter only requires minor tuning for the localization of the maximum. There is non-trivial dependence of the energy radiated on eccentricity (several local maxima and minima). Correlations with orbital dynamics shortly before merger are discussed. While being strongly gauge dependent, these findings are intuitive from a physical point of view and support basic ideas about the efficiency of gravitational radiation from a binary system., Comment: 9 pages, 6 figures, Amaldi8 conference proceedings as published

The success of the moving puncture method for the numerical simulation of black hole systems can be partially explained by the properties of stationary solutions of the 1+log coordinate condition. We compute stationary 1+log slices of the Schwarzschild spacetime in isotropic coordinates in order to investigate the coordinate singularity that the numerical methods have to handle at the puncture. We present an alternative integration method to obtain isotropic coordinates that simplifies numerical integration and that gives direct access to a local expansion in the isotropic radius near the puncture. Numerical results have shown that certain quantities are well approximated by a function linear in the isotropic radius near the puncture, while here we show that in some cases the isotropic radius appears with an exponent that is close to but unequal to one., Comment: 11 pages, 5 figures, pdflatex

Moving puncture simulations of black hole binaries rely on a specific gauge choice that leads to approximately stationary coordinates near each black hole. Part of the shift condition is a damping parameter, which has to be properly chosen for stable evolutions. However, a constant damping parameter does not account for the difference in mass in unequal mass binaries. We introduce a position dependent shift damping that addresses this problem. Although the coordinates change, the changes in the extracted gravitational waves are small., Comment: 15 pages, submitted to CQG for NRDA 2009 conference proceedings

We present a preliminary study of the multipolar structure of gravitational radiation from spinning black hole binary mergers. We consider three different spinning binary configurations: (1) one "hang-up" run, where the black holes have equal masses and large spins initially aligned with the orbital angular momentum; (2) seven "spin-flip" runs, where the holes have a mass ratio q=4, the spins are anti-aligned with the orbital angular momentum, and the initial Kerr parameters of the holes j_1=j_2=j_i are fine-tuned to produce a Schwarzschild remnant after merger; (3) three "super-kick" runs where the mass ratio q=M_1/M_2=1, 2, 4 and the spins of the two holes are initially located on the orbital plane, pointing in opposite directions. For all of these simulations we compute the multipolar energy distribution and the Kerr parameter of the final hole. For the hang-up run, we show that including leading-order spin-orbit and spin-spin terms in a multipolar decomposition of the post-Newtonian waveforms improves the agreement with the numerical simulation., corrected minor typos in Eqs.(2),(3); final version accepted by CQG

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., 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 2006

We describe a modification of a fourth-order accurate ``moving puncture'' evolution code, where by replacing spatial fourth-order accurate differencing operators in the bulk of the grid by a specific choice of sixth-order accurate stencils we gain significant improvements in accuracy. We illustrate the performance of the modified algorithm with an equal-mass simulation covering nine orbits., Comment: 13 pages, 6 figures

Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to higher-mass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors., Minor changes, Submitted to Class. Quantum Grav. (Proc. GWDAW11)

We perform an analytic late time analysis for maximal slicing of the Reissner-Nordstr\"om black hole spacetime. In particular, we discuss the collapse of the lapse in terms of its late time behavior at the throat and at the event horizon for the even and the puncture lapse. In the latter case we also determine the value of the lapse at the puncture. Furthermore, in the limit of late times slice stretching effects are studied as they arise for maximal slicing of puncture evolutions. We perform numerical experiments for a Schwarzschild black hole with puncture lapse and find agreement with the analytical results., Comment: 22 pages, 4 figures, changes in Sec. IV B, as published in Phys. Rev. D

We prove by explicit construction that there exists a maximal slicing of the Schwarzschild spacetime such that the lapse has zero gradient at the puncture. This boundary condition has been observed to hold in numerical evolutions, but in the past it was not clear whether the numerically obtained maximal slices exist analytically. We show that our analytical result agrees with numerical simulation. Given the analytical form for the lapse, we can derive that at late times the value of the lapse at the event horizon approaches the value ${3/16}\sqrt{3} \approx 0.3248$, justifying the numerical estimate of 0.3 that has been used for black hole excision in numerical simulations. We present our results for the non-extremal Reissner-Nordstr\"om metric, generalizing previous constructions of maximal slices., Comment: 21 pages, 9 figures, published version with changes to Sec. VII

The strong-field region inside a black hole needs special attention during numerical simulation. One approach for handling the problem is the moving puncture method, which has become an important tool in numerical relativity since it allows long term simulations of binary black holes. An essential component of this method is the choice of the '1+log'-slicing condition. We present an investigation of this slicing condition in rotating black hole spacetimes. We discuss how the results of the stationary Schwarzschild '1+log'-trumpet change when spin is added. This modification enables a simple and cheap algorithm for determining the spin of a non-moving black hole for this particular slicing condition. Applicability of the algorithm is verified in simulations of single black hole, binary neutron star and mixed binary simulations., conference

We present a new technique for the numerical simulation of axisymmetric systems. This technique avoids the coordinate singularities which often arise when cylindrical or polar-spherical coordinate finite difference grids are used, particularly in simulating tensor partial differential equations like those of 3+1 numerical relativity. For a system axisymmetric about the z axis, the basic idea is to use a 3-dimensional Cartesian (x,y,z) coordinate grid which covers (say) the y=0 plane, but is only one finite-difference-molecule--width thick in the y direction. The field variables in the central y=0 grid plane can be updated using normal (x,y,z)--coordinate finite differencing, while those in the y \neq 0 grid planes can be computed from those in the central plane by using the axisymmetry assumption and interpolation. We demonstrate the effectiveness of the approach on a set of fully nonlinear test computations in 3+1 numerical general relativity, involving both black holes and collapsing gravitational waves., Comment: 17 pages, 4 figures

We present the first results in a new program intended to make the best use of all available technologies to provide an effective understanding of waves from inspiralling black hole binaries in time for imminent observations. In particular, we address the problem of combining the close-limit approximation describing ringing black holes and full numerical relativity, required for essentially nonlinear interactions. We demonstrate the effectiveness of our approach using general methods for a model problem, the head-on collision of black holes. Our method allows a more direct physical understanding of these collisions indicating clearly when non-linear methods are important. The success of this method supports our expectation that this unified approach will be able to provide astrophysically relevant results for black hole binaries in time to assist gravitational wave observations., Comment: 4 pages, 3 eps figures, Revtex

The standard approach to the numerical evolution of black hole data using the ADM formulation with maximal slicing and vanishing shift is extended to non-symmetric black hole data containing black holes with linear momentum and spin by using a time-independent conformal rescaling based on the puncture representation of the black holes. We give an example for a concrete three dimensional numerical implementation. The main result of the simulations is that this approach allows for the first time to evolve through a brief period of the merger phase of the black hole inspiral., 8 pages, 9 figures, REVTeX; expanded discussion, results unchanged

In this work we present a formulation of the Teukolsky equation for generic spin perturbations on the hyperboloidal and horizon penetrating foliation of Kerr recently proposed by Racz and Toth. An additional, spin-dependent rescaling of the field variable can be used to achieve stable, long-term, and accurate time-domain evolutions of generic spin perturbations. As an application (and a severe numerical test), we investigate the late-time decays of electromagnetic and gravitational perturbations at the horizon and future null infinity by means of 2+1 evolutions. As initial data we consider four combinations of (non-)stationary and (non-)compact-support initial data with a pure spin-weighted spherical harmonic profile. We present an extensive study of late time decays of axisymmetric perturbations. We verify the power-law decay rates predicted analytically, together with a certain "splitting" behaviour of the power-law exponent. We also present results for non-axisymmetric perturbations. In particular, our approach allows to study the behaviour of the late time decays of gravitational fields for nearly extremal and extremal black holes. For rapid rotation we observe a very prolonged, weakly damped, quasi-normal-mode phase. For extremal rotation the field at future null infinity shows an oscillatory behaviour decaying as the inverse power of time, while at the horizon it is amplified by several orders of magnitude over long time scales. This behaviour can be understood in terms of the superradiance cavity argument.

Although an important issue in canonical quantization, the problem of representing the constraint algebra in the loop representation of quantum gravity has received little attention. The only explicit computation was performed by Gambini, Garat and Pullin for a formal point-splitting regularization of the diffeomorphism and Hamiltonian constraints. It is shown that the calculation of the algebra simplifies considerably when the constraints are expressed not in terms of generic area derivatives but rather as the specific shift operators that reflect the geometric meaning of the constraints., Comment: 17 pages, LaTeX, 2 figures included with epsf

We find a consistent formulation of the constraints of Quantum Gravity with a cosmological constant in terms of the Ashtekar new variables in the connection representation, including the existence of a state that is a solution to all the constraints. This state is related to the Chern-Simons form constructed from the Ashtekar connection and has an associated metric in spacetime that is everywhere nondegenerate. We then transform this state to the loop representation and find solutions to all the constraint equations for intersecting loops. These states are given by suitable generalizations of the Jones knot polynomial for the case of intersecting knots. These are the first physical states of Quantum Gravity for which an explicit form is known both in the connection and loop representations. Implications of this result are also discussed., Comment: 23pp