Researcher profile

Bernd Bruegmann

Bernd Bruegmann contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 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.

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.

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.

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.

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.

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.

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%.

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.

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.

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.

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.

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.