BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Bernd Bruegmann

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.SRastro-ph.COphysics.comp-phastro-ph.GAmath-phmath.MP

Catalog footprint

What is connected

35works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

35 published item(s)

preprint2021arXiv

High-accuracy simulations of highly spinning binary neutron star systems

With an increasing number of expected gravitational-wave detections of binary neutron star mergers, it is essential that gravitational-wave models employed for the analysis of observational data are able to describe generic compact binary systems. This includes systems in which the individual neutron stars are millisecond pulsars for which spin effects become essential. In this work, we perform numerical-relativity simulations of binary neutron stars with aligned and anti-aligned spins within a range of dimensionless spins of $χ\sim [-0.28,0.58]$. The simulations are performed with multiple resolutions, show a clear convergence order and, consequently, can be used to test existing waveform approximants. We find that for very high spins gravitational-wave models that have been employed for the interpretation of GW170817 and GW190425 are not capable of describing our numerical-relativity dataset. We verify through a full parameter estimation study in which clear biases in the estimate of the tidal deformability and effective spin are present. We hope that in preparation of the next gravitational-wave observing run of the Advanced LIGO and Advanced Virgo detectors our new set of numerical-relativity data can be used to support future developments of new gravitational-wave models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Investigating the mass-ratio dependence of the prompt-collapse threshold with numerical-relativity simulations

The next observing runs of advanced gravitational-wave detectors will lead to a variety of binary neutron star detections and numerous possibilities for multi-messenger observations of binary neutron star systems. In this context a clear understanding of the merger process and the possibility of prompt black hole formation after merger is important, as the amount of ejected material strongly depends on the merger dynamics. These dynamics are primarily affected by the total mass of the binary, however, the mass ratio also influences the postmerger evolution. To determine the effect of the mass ratio, we investigate the parameter space around the prompt-collapse threshold with a new set of fully relativistic simulations. The simulations cover three equations of state and seven mass ratios in the range of $1.0 \leq q \leq 1.75$, with five to seven simulations of binary systems of different total mass in each case. The threshold mass is determined through an empirical relation based on the collapse-time, which allows us to investigate effects of the mass-ratio on the threshold mass and also on the properties of the remnant system. Furthermore, we model effects of mass ratio and equation of state on tidal parameters of threshold configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analytical and numerical treatment of perturbed black holes in horizon-penetrating coordinates

The deviations of non-linear perturbations of black holes from the linear case are important in the context of ringdown signals with large signal-to-noise ratio. To facilitate a comparison between the two we derive several results of linear perturbation theory in coordinates which may be adopted in numerical work. Specifically, our results are derived in Kerr-Schild coordinates adjusted by a general height function. In the first part of the paper we address the questions: for an initial configuration of a massless scalar field, what is the amplitude of the excited quasinormal mode (QNM) for any observer outside outside the event horizon, and furthermore what is the resulting tail contribution? This is done by constructing the full Green's function for the problem with exact solutions of the confluent Heun equation satisfying appropriate boundary conditions. In the second part of the paper, we detail new developments to our pseudospectral numerical relativity code bamps to handle scalar fields. In the linear regime we employ precisely the Kerr-Schild coordinates treated by our previous analysis. In particular, we evolve pure QNM type initial data along with several other types of initial data and report on the presence of overtone modes in the signal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The evolution of hyperboloidal data with the dual foliation formalism: Mathematical analysis and wave equation tests

A long-standing problem in numerical relativity is the satisfactory treatment of future null-infinity. We propose an approach for the evolution of hyperboloidal initial data in which the outer boundary of the computational domain is placed at infinity. The main idea is to apply the `dual foliation' formalism in combination with hyperboloidal coordinates and the generalized harmonic gauge formulation. The strength of the present approach is that, following the ideas of Zenginoglu, a hyperboloidal layer can be naturally attached to a central region using standard coordinates of numerical relativity applications. Employing a generalization of the standard hyperboloidal slices, developed by Calabrese et. al., we find that all formally singular terms take a trivial limit as we head to null-infinity. A byproduct is a numerical approach for hyperboloidal evolution of nonlinear wave equations violating the null-condition. The height-function method, used often for fixed background spacetimes, is generalized in such a way that the slices can be dynamically `waggled' to maintain the desired outgoing coordinate lightspeed precisely. This is achieved by dynamically solving the eikonal equation. As a first numerical test of the new approach we solve the 3D flat space scalar wave equation. The simulations, performed with the pseudospectral bamps code, show that outgoing waves are cleanly absorbed at null-infinity and that errors converge away rapidly as resolution is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Pseudospectral Method for Gravitational Wave Collapse

We present a new pseudospectral code, bamps, for numerical relativity written with the evolution of collapsing gravitational waves in mind. We employ the first order generalized harmonic gauge formulation. The relevant theory is reviewed and the numerical method is critically examined and specialized for the task at hand. In particular we investigate formulation parameters, gauge and constraint preserving boundary conditions well-suited to non-vanishing gauge source functions. Different types of axisymmetric twist-free moment of time symmetry gravitational wave initial data are discussed. A treatment of the axisymmetric apparent horizon condition is presented with careful attention to regularity on axis. Our apparent horizon finder is then evaluated in a number of test cases. Moving on to evolutions, we investigate modifications to the generalized harmonic gauge constraint damping scheme to improve conservation in the strong field regime. We demonstrate strong-scaling of our pseudospectral penalty code. We employ the Cartoon method to efficiently evolve axisymmetric data in our 3+1 dimensional code. We perform test evolutions of Schwarzschild perturbed by gravitational waves and by gauge pulses, both to demonstrate the use of our blackhole excision scheme and for comparison with earlier results. Finally numerical evolutions of supercritical Brill waves are presented to demonstrate durability of the excision scheme for the dynamical formation of a blackhole.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Binary Neutron Stars with Generic Spin, Eccentricity, Mass ratio, and Compactness - Quasi-equilibrium Sequences and First Evolutions

Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $\sim10^{-2}M_\odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

We study equal and unequal-mass neutron star mergers by means of new numerical relativity simulations in which the general relativistic hydrodynamics solver employs an algorithm that guarantees mass conservation across the refinement levels of the computational mesh. We consider eight binary configurations with total mass $M=2.7\,M_\odot$, mass-ratios $q=1$ and $q=1.16$, and four different equation of states (EOSs), and one configuration with a stiff EOS, $M=2.5M_\odot$ and $q=1.5$. We focus on the post-merger dynamics and study the merger remnant, dynamical ejecta and the postmerger gravitational wave spectrum. Although most of the merger remnants form a hypermassive neutron star collapsing to a black hole+disk system on dynamical timescales, stiff EOSs can eventually produce a stable massive neutron star. Ejecta are mostly emitted around the orbital plane; favored by large mass ratios and softer EOS. The postmerger wave spectrum is mainly characterized by non-axisymmetric oscillations of the remnant. The stiff EOS configuration consisting of a $1.5M_\odot$ and a $1.0M_\odot$ neutron star shows a rather peculiar dynamics. During merger the companion star is very deformed; about~$\sim0.03M_\odot$ of rest-mass becomes unbound from the tidal tail due torque; and the merger remnant forms stable neutron star surrounded by a massive accretion disk $\sim0.3M_\odot$. Similar configurations might be particularly interesting for electromagnetic counterparts. Comparing results obtained with and without the conservative mesh refinement algorithm, we find that post-merger simulations can be affected by systematic errors if mass conservation is not enforced in the mesh refinement strategy. However, mass conservation also depends on grid details and on the artificial atmosphere setup. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solving 3D relativistic hydrodynamical problems with WENO discontinuous Galerkin methods

Discontinuous Galerkin (DG) methods coupled to WENO algorithms allow high order convergence for smooth problems and for the simulation of discontinuities and shocks. In this work, we investigate WENO-DG algorithms in the context of numerical general relativity, in particular for general relativistic hydrodynamics. We implement the standard WENO method at different orders, a compact (simple) WENO scheme, as well as an alternative subcell evolution algorithm. To evaluate the performance of the different numerical schemes, we study non-relativistic, special relativistic, and general relativistic testbeds. We present the first three-dimensional simulations of general relativistic hydrodynamics, albeit for a fixed spacetime background, within the framework of WENO-DG methods. The most important testbed is a single TOV-star in three dimensions, showing that long term stable simulations of single isolated neutron stars can be obtained with WENO-DG methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Initial data for binary neutron stars with adjustable eccentricity

Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipse's inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mergers of binary neutron stars with realistic spin

Simulations of binary neutron stars have seen great advances in terms of physical detail and numerical quality. However, the spin of the neutron stars, one of the simplest global parameters of binaries, remains mostly unstudied. We present the first, fully nonlinear general relativistic dynamical evolutions of the last three orbits for constraint satisfying initial data of spinning neutron star binaries, with astrophysically realistic spins aligned and anti-aligned to the orbital angular momentum. The initial data is computed with the constant rotational velocity approach. The dynamics of the systems is analyzed in terms of gauge-invariant binding energy vs. orbital angular momentum curves. By comparing to a binary black hole configuration we can estimate the different tidal and spin contributions to the binding energy for the first time. First results on the gravitational wave forms are presented. The phase evolution during the orbital motion is significantly affected by spin-orbit interactions, leading to delayed or early mergers. Furthermore, a frequency shift in the main emission mode of the hyper massive neutron star is observed. Our results suggest that a detailed modeling of merger waveforms requires the inclusion of spin, even for the moderate magnitudes observed in binary neutron star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spinning black hole in the puncture method: Numerical experiments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Collapse of Nonlinear Gravitational Waves in Moving-Puncture Coordinates

We study numerical evolutions of nonlinear gravitational waves in moving-puncture coordinates. We adopt two different types of initial data -- Brill and Teukolsky waves -- and evolve them with two independent codes producing consistent results. We find that Brill data fail to produce long-term evolutions for common choices of coordinates and parameters, unless the initial amplitude is small, while Teukolsky wave initial data lead to stable evolutions, at least for amplitudes sufficiently far from criticality. The critical amplitude separates initial data whose evolutions leave behind flat space from those that lead to a black hole. For the latter we follow the interaction of the wave, the formation of a horizon, and the settling down into a time-independent trumpet geometry. We explore the differences between Brill and Teukolsky data and show that for less common choices of the parameters -- in particular negative amplitudes -- Brill data can be evolved with moving-puncture coordinates, and behave similarly to Teukolsky waves.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Eccentric binary neutron star mergers

Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circular. In addition to gravitational waves generated by the orbital motion, we find that the signal also contains a strong component due to stellar oscillations (f-modes) induced by tidal forces, extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on the order of 10% of the total initial mass.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical solution of the 2+1 Teukolsky equation on a hyperboloidal and horizon penetrating foliation of Kerr and application to late-time decays

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solving the Hamiltonian constraint for 1+log trumpets

The puncture method specifies black hole data on a hypersurface with the aid of a conformal rescaling of the metric that exhibits a coordinate singularity at the puncture point. When constructing puncture initial data by solving the Hamiltonian constraint for the conformal factor, the coordinate singularity requires special attention. The standard way to treat the pole singularity occurring in wormhole puncture data is not generally applicable to trumpet puncture data. We investigate a new approach based on inverse powers of the conformal factor and present numerical examples for single punctures of the wormhole and 1+log-trumpet type. Additionally, we describe a method to solve the Hamiltonian constraint for two 1+log trumpets for a given extrinsic curvature with non-vanishing trace. We investigate properties of this constructed initial data during binary black hole evolutions and find that the initial gauge dynamics is reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accuracy of numerical relativity waveforms from binary neutron star mergers and their comparison with post-Newtonian waveforms

We present numerical relativity simulations of nine-orbit equal-mass binary neutron star covering the quasicircular late inspiral and merger. The extracted gravitational waveforms are analyzed for convergence and accuracy. Second order convergence is observed up to contact, i.e. about 3-4 cycles to merger; error estimates can be made up to this point. The uncertainties on the phase and the amplitude are dominated by truncation errors and can be minimized to 0.13 rad and less then 1%, respectively, by using several simulations and extrapolating in resolution. In the latter case finite-radius extraction uncertainties become a source of error of the same order and have to be taken into account. The waveforms are tested against accuracy standards for data analysis. The uncertainties on the waveforms are such that accuracy standards are generically not met for signal-to-noise ratios relevant for detection, except for some best cases using extrapolation from several runs. A detailed analysis of the errors is thus imperative for the use of numerical relativity waveforms from binary neutron stars in quantitative studies. The waveforms are compared with the post-Newtonian Taylor T4 approximants both for point-particle and including the analytically known tidal corrections. The T4 approximants accumulate significant phase differences of 2 rad at contact and 4 rad at merger, underestimating the influence of finite size effects. Tidal signatures in the waveforms are thus important at least during the last six orbits of the merger process.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Compact binary evolutions with the Z4c formulation

Numerical relativity simulations of compact binaries with the Z4c and BSSNOK formulations are compared. The Z4c formulation is advantageous in every case considered. In simulations of non-vacuum spacetimes the constraint violations due to truncation errors are between one and three orders of magnitude lower in the Z4c evolutions. Improvements are also found in the accuracy of the computed gravitational radiation. For equal-mass irrotational binary neutron star evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of four smaller. The quality of the Z4c numerical data is also demonstrated by a remarkably accurate computation of the ADM mass from surface integrals. For equal-mass non-spinning binary puncture black hole evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of two smaller. In the same evolutions we find that away from the punctures the Hamiltonian constraint violation is reduced by between one and two orders of magnitude. Furthermore, the utility of gravitational radiation controlling, constraint preserving boundary conditions for the Z4c formulation is demonstrated. The evolution of spacetimes containing a single compact object confirm earlier results in spherical symmetry. The boundary conditions avoid spurious and non-convergent effects present in high resolution runs with either formulation with a more naive boundary treatment. We conclude that Z4c is preferable to BSSNOK for the numerical solution of the 3+1 Einstein equations with the puncture gauge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Eccentric black hole mergers and zoom-whirl behavior from elliptic inspirals to hyperbolic encounters

We perform a parameter study of non-spinning, equal and unequal mass black hole binaries on generic, eccentric orbits in numerical relativity. The linear momentum considered ranges from that of a circular orbit to ten times that value. We discuss the different manifestations of zoom-whirl behavior in the hyperbolic and the elliptic regime. The hyperbolic data set applies to dynamical capture scenarios (e.g. in globular clusters). Evolutions in the elliptic regime correspond to possible end states of supermassive black hole binaries. We spot zoom-whirl behavior for eccentricities as low as $e\sim0.5$, i.e. within the expected range of eccentricities in massive black hole binaries from galaxy mergers and binaries near galactic centers. The resulting gravitational waveforms reveal a rich structure, which will effectively break degeneracies in parameter space improving parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal effects in binary neutron star coalescence

We compare dynamics and waveforms from binary neutron star coalescence as computed by new long-term ($\sim 10 $ orbits) numerical relativity simulations and by the tidal effective-one-body (EOB) model including analytical tidal corrections up to second post-Newtonian order (2PN). The current analytical knowledge encoded in the tidal EOB model is found to be sufficient to reproduce the numerical data up to contact and within their uncertainties. Remarkably, no calibration of any tidal EOB free parameters is required, beside those already fitted to binary black holes data. The inclusion of 2PN tidal corrections minimizes the differences with the numerical data, but it is not possible to significantly distinguish them from the leading-order tidal contribution. The presence of a relevant amplification of tidal effects is likely to be excluded, although it can appear as a consequence of numerical inaccuracies. We conclude that the tidally-completed effective-one-body model provides nowadays the most advanced and accurate tool for modelling gravitational waveforms from binary neutron star inspiral up to contact. This work also points out the importance of extensive tests to assess the uncertainties of the numerical data, and the potential need of new numerical strategies to perform accurate simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A pseudospectral matrix method for time-dependent tensor fields on a spherical shell

We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test case we consider the evolution of a single black hole in numerical general relativity. A natural strategy would be the expansion in tensor spherical harmonics in spherical coordinates. Instead, we consider the simpler and potentially more efficient possibility of a double Fourier expansion on the sphere for tensors in Cartesian coordinates. As usual for the double Fourier method, we employ a filter to address time-step limitations and certain stability issues. We find that a tensor filter based on spin-weighted spherical harmonics is successful, while two simplified, non-spin-weighted filters do not lead to stable evolutions. The derivatives and the filter are implemented by matrix multiplication for efficiency. A key technical point is the construction of a matrix multiplication method for the spin-weighted spherical harmonic filter. As example for the efficient parallelization of the double Fourier, spin-weighted filter method we discuss an implementation on a GPU, which achieves a speed-up of up to a factor of 20 compared to a single core CPU implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Black-hole binaries go to eleven orbits

We analyse an eleven-orbit inspiral of a non-spinning black-hole binary with mass ratio q=M1/M2=4. The numerically obtained gravitational waveforms are compared with post-Newtonian (PN) predictions including several sub-dominant multipoles up to multipolar indices (l=5,m=5). We find that (i) numerical and post-Newtonian predictions of the phase of the (2,2) mode accumulate a phase difference of about 0.35 rad at the PN cut off frequency 0.1 for the Taylor T1 approximant when numerical and PN waveforms are matched over a window in the early inspiral phase; (ii) in contrast to previous studies of equal-mass and specific spinning binaries, we find the Taylor T4 approximant to agree less well with numerical results, provided the latter are extrapolated to infinite extraction radius; (iii) extrapolation of gravitational waveforms to infinite extraction radius is particularly important for subdominant multipoles with l unequal m; (iv) 3PN terms in post-Newtonian multipole expansions significantly improve the agreement with numerical predictions for sub-dominant multipoles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Characterization of the gravitational wave emission of three black holes

We study the gravitational wave emission of three compact objects using post-Newtonian (PN) equations of motion derived from the Arnowitt-Deser-Misner Hamiltonian formulation, where we include (for the first time in this context) terms up to 2.5 PN order. We perform numerical simulations of a hierarchical configuration of three compact bodies in which a binary system is perturbed by a third, lighter body initially far away from the binary. The relative importance of the different PN orders is examined. We compute the waveform in the linear regime considering mass quadrupole, current quadrupole and mass octupole contributions. Performing a spherical harmonic decomposition of the waveforms we find that from the l = 3 modes it is possible to extract information about the third body, in particular, the period, eccentricity of its orbit, and the inclination angle between the inner and outer binary orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical relativity simulations of binary neutron stars

We present a new numerical relativity code designed for simulations of compact binaries involving matter. The code is an upgrade of the BAM code to include general relativistic hydrodynamics and implements state-of-the-art high-resolution-shock-capturing schemes on a hierarchy of mesh refined Cartesian grids with moving boxes. We test and validate the code in a series of standard experiments involving single neutron star spacetimes. We present test evolutions of quasi-equilibrium equal-mass irrotational binary neutron star configurations in quasi-circular orbits which describe the late inspiral to merger phases. Neutron star matter is modeled as a zero-temperature fluid; thermal effects can be included by means of a simple ideal-gas prescription. We analyze the impact that the use of different values of damping parameter in the Gamma-driver shift condition has on the dynamics of the system. The use of different reconstruction schemes and their impact in the post-merger dynamics is investigated. We compute and characterize the gravitational radiation emitted by the system. Self-convergence of the waves is tested, and we consistently estimate error-bars on the numerically generated waveforms in the inspiral phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The trumpet solution from spherical gravitational collapse with puncture gauges

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical shift condition for unequal mass black hole binaries

Certain numerical frameworks used for the evolution of binary black holes make use of a gamma driver, which includes a damping factor. Such simulations typically use a constant value for damping. However, it has been found that very specific values of the damping factor are needed for the calculation of unequal mass binaries. We examine carefully the role this damping plays, and provide two explicit, non-constant forms for the damping to be used with mass-ratios further from one. Our analysis of the resultant waveforms compares well against the constant damping case.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical black hole initial data with low eccentricity based on post-Newtonian orbital parameters

Black hole binaries on non-eccentric orbits form an important subclass of gravitational wave sources, but it is a non-trivial issue to construct numerical initial data with minimal initial eccentricity for numerical simulations. We compute post-Newtonian orbital parameters for quasi-spherical orbits using the method of Buonanno, Chen and Damour (2006) and examine the resulting eccentricity in numerical simulations. Four different methods are studied resulting from the choice of Taylor-expanded or effective-one-body Hamiltonians, and from two choices for the energy flux. The eccentricity increases for unequal masses and for spinning black holes, but remains smaller than that obtained from previous post-Newtonian approaches. The effective-one-body Hamiltonian offers advantages for decreasing initial separation as expected, but in the context of this study also performs significantly better than the Taylor-expanded Hamiltonian for binaries with spin.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical evolution of multiple black holes with accurate initial data

We present numerical evolutions of three equal-mass black holes using the moving puncture approach. We calculate puncture initial data for three black holes solving the constraint equations by means of a high-order multigrid elliptic solver. Using these initial data, we show the results for three black hole evolutions with sixth-order waveform convergence. We compare results obtained with the BAM and AMSS-NCKU codes with previous results. The approximate analytic solution to the Hamiltonian constraint used in previous simulations of three black holes leads to different dynamics and waveforms. We present some numerical experiments showing the evolution of four black holes and the resulting gravitational waveform.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiation from low-momentum zoom-whirl orbits

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations of black-hole binaries with unequal masses or non-precessing spins: accuracy, physical properties, and comparison with post-Newtonian results

We present gravitational waveforms for the last orbits and merger of black-hole-binary (BBH) systems along two branches of the BBH parameter space: equal-mass binaries with equal non-precessing spins, and nonspinning unequal-mass binaries. The waveforms are calculated from numerical solutions of Einstein's equations for black-hole binaries that complete between six and ten orbits before merger. Along the equal-mass spinning branch, the spin parameter of each BH is $χ_i = S_i/M_i^2 \in [-0.85,0.85]$, and along the unequal-mass branch the mass ratio is $q =M_2/M_1 \in [1,4]$. We discuss the construction of low-eccentricity puncture initial data for these cases, the properties of the final merged BH, and compare the last 8-10 GW cycles up to $Mω= 0.1$ with the phase and amplitude predicted by standard post-Newtonian (PN) approximants. As in previous studies, we find that the phase from the 3.5PN TaylorT4 approximant is most accurate for nonspinning binaries. For equal-mass spinning binaries the 3.5PN TaylorT1 approximant (including spin terms up to only 2.5PN order) gives the most robust performance, but it is possible to treat TaylorT4 in such a way that it gives the best accuracy for spins $χ_i > -0.75$. When high-order amplitude corrections are included, the PN amplitude of the $(\ell=2,m=\pm2)$ modes is larger than the NR amplitude by between 2-4%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symplectic Integration of Post-Newtonian Equations of Motion with Spin

We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In the context of the splitting, it is possible to integrate the individual terms of the spin-orbit and spin-spin Hamiltonians analytically, exploiting the special structure of the underlying equations of motion. The outcome is a symplectic, time-reversible integrator, which can be raised to arbitrary order by composition. A fourth-order version is shown to give excellent behavior concerning error growth and conservation of energy and angular momentum in long-term simulations. Favorable properties of the integrator are retained in the presence of weak dissipative forces due to radiation damping in the full post-Newtonian equations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Black-hole binary simulations: the mass ratio 10:1

We present the first numerical simulations of an initially non-spinning black-hole binary with mass ratio as large as 10:1 in full general relativity. The binary completes approximately 3 orbits prior to merger and radiates about 0.415% of the total energy and 12.48% of the initial angular momentum in the form of gravitational waves. The single black hole resulting from the merger acquires a kick of about 66.7 km/s relative to the original center of mass frame. The resulting gravitational waveforms are used to validate existing formulas for the recoil, final spin and radiated energy over a wider range of the symmetric mass ratio parameter eta=M1*M2/(M1+M2)^2 than previously possible. The contributions of l > 2 multipoles are found to visibly influence the gravitational wave signal obtained at fixed inclination angles.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Schwarzschild black hole as moving puncture in isotropic coordinates

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Samurai Project: verifying the consistency of black-hole-binary waveforms for gravitational-wave detection

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Toward a dynamical shift condition for unequal mass black hole binary simulations

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.

0 citations0 reviewsNo code linkedOpen access