Source author record

Matthew Giesler

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2022arXiv

Multipole moments on the common horizon in a binary-black-hole simulation

We construct the covariantly defined multipole moments on the common horizon of an equal-mass, non-spinning, quasicircular binary-black-hole system. We see a strong correlation between these multipole moments and the gravitational waveform. We find that the multipole moments are well described by the fundamental quasinormal modes at sufficiently late times. For each multipole moment, at least two fundamental modes of different $\ell$ are detectable in the best model. These models provide faithful estimates of the true mass and spin of the remnant black hole. We also show that by including overtones, the $\ell=m=2$ mass multipole moment admits an excellent quasinormal-mode description at all times after the merger. This demonstrates the perhaps surprising power of perturbation theory near the merger.

preprint2022arXiv

Redshift factor and the small mass-ratio limit in binary black hole simulations

We present a calculation of the Detweiler redshift factor in binary black hole simulations based on its relation to the surface gravity. The redshift factor has far-reaching applications in analytic approximations, gravitational self-force calculations, and conservative two-body dynamics. By specializing to non-spinning, quasi-circular binaries with mass ratios ranging from $m_A/m_B = 1$ to $m_A/m_B = 9.5$ we are able to recover the leading small-mass-ratio (SMR) prediction with relative differences of order $10^{-5}$ from simulations alone. The next-to-leading order term that we extract agrees with the SMR prediction arising from self-force calculations, with differences of a few percent. These deviations from the first-order conservative prediction are consistent with non-adiabatic effects that can be accommodated in an SMR expansion. This fact is also supported by a comparison to the conservative post-Newtonian prediction of the redshifts. For the individual redshifts, a re-expansion in terms of the symmetric mass ratio $ν$ does not improve the convergence of the series. However we find that when looking at the sum of the redshift factors of both back holes, $z_A + z_B$, which is symmetric under the exchange of the masses, a re-expansion in $ν$ accelerates its convergence. Our work provides further evidence of the surprising effectiveness of SMR approximations in modeling even comparable mass binary black holes.

preprint2022arXiv

Simulating magnetized neutron stars with discontinuous Galerkin methods

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$ΛΠ^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method.

preprint2022arXiv

Targeted large mass ratio numerical relativity surrogate waveform model for GW190814

Gravitational wave observations of large mass ratio compact binary mergers like GW190814 highlight the need for reliable, high-accuracy waveform templates for such systems. We present NRHybSur2dq15, a new surrogate model trained on hybridized numerical relativity (NR) waveforms with mass ratios $q\leq15$, and aligned spins $|χ_{1z}|\leq0.5$ and $χ_{2z}=0$. We target the parameter space of GW190814-like events as large mass ratio NR simulations are very expensive. The model includes the (2,2), (2,1), (3,3), (4,4), and (5,5) spin-weighted spherical harmonic modes, and spans the entire LIGO bandwidth (with $f_{\mathrm{low}}=20$ Hz) for total masses $M \gtrsim 9.5 \, M_{\odot}$. NRHybSur2dq15 accurately reproduces the hybrid waveforms, with mismatches below $\sim 2 \times 10^{-3}$ for total masses $10 \, M_{\odot} \leq M \leq 300 \, M_{\odot}$. This is at least an order of magnitude improvement over existing semi-analytical models for GW190814-like systems. Finally, we reanalyze GW190814 with the new model and obtain source parameter constraints consistent with previous work.

preprint2020arXiv

Black hole ringdown: the importance of overtones

It is possible to infer the mass and spin of the remnant black hole from binary black hole mergers by comparing the ringdown gravitational wave signal to results from studies of perturbed Kerr spacetimes. Typically these studies are based on the fundamental quasinormal mode of the dominant $\ell=m=2$ harmonic. By modeling the ringdown of accurate numerical relativity simulations, we find that the fundamental mode alone is insufficient to recover the true underlying mass and spin, unless the analysis is started very late in the ringdown. Including higher overtones associated with this $\ell=m=2$ harmonic resolves this issue, and provides an unbiased estimate of the true remnant parameters. Further, including overtones allows for the modeling of the ringdown signal for all times beyond the peak strain amplitude, indicating that the linear quasinormal regime starts much sooner than previously expected. A model for the ringdown beginning at the peak strain amplitude can exploit the higher signal-to-noise ratio in detectors, reducing uncertainties in the extracted remnant quantities. Tests of the no-hair theorem should consider incorporating overtones in the analysis.

preprint2020arXiv

Spin and Eccentricity Evolution in Triple Systems: from the Lidov-Kozai Interaction to the Final Merger of the Inner Binary

We study the spin and eccentricity evolution of black-hole (BH) binaries that are perturbed by tertiary masses and experience the Lidov-Kozai (LK) excitation. We focus on three aspects. Firstly, we study the spin-orbit alignment of the inner binary following the approach outlined by Antonini et al. [MNRAS 480, L58 (2018)] and Liu and Lai [ApJ 863, 68 (2018)], yet allowing the spins to have random initial orientations. We confirm the existence of a dynamical attractor that drives the spin-orbit angle at the end of the LK evolution to a value given by the initial angle between the spin and the outer orbital angular momentum (instead of to a specific value of the effective spin). Secondly, we follow the (inner) binary's evolution further to the merger to study the final spin-spin alignment. We generalize the effective potential theory to include orbital eccentricity, which allows us to efficiently evolve the system in the early inspiral stages. We further find that the spin-spin and spin-orbit alignments are correlated and the correlation is determined by the initial spin-orbit angle. For systems with the spin vectors initially in the orbital plane, the final spins strongly disfavor an aligned configuration and could thus lead to a greater value of the GW recoil than a uniform spin-spin alignment would predict. Lastly, we study the maximum eccentricity excitation that can be achieved during the LK process, including the effects of gravitational-wave radiation. We find that when the tertiary mass is a super-massive BH and the inner binary is massive, then even with the maximum LK excitation, the residual eccentricity is typically less than 0.1 when the binary's orbital frequency reaches 10 Hz, and a decihertz detector would be necessary to follow such a system's orbital evolution.

preprint2015arXiv

Nearly extremal apparent horizons in simulations of merging black holes

The spin angular momentum $S$ of an isolated Kerr black hole is bounded by the surface area $A$ of its apparent horizon: $8πS \le A$, with equality for extremal black holes. In this paper, we explore the extremality of individual and common apparent horizons for merging, rapidly spinning binary black holes. We consider simulations of merging black holes with equal masses $M$ and initial spin angular momenta aligned with the orbital angular momentum, including new simulations with spin magnitudes up to $S/M^2 = 0.994$. We measure the area and (using approximate Killing vectors) the spin on the individual and common apparent horizons, finding that the inequality $8πS < A$ is satisfied in all cases but is very close to equality on the common apparent horizon at the instant it first appears. We also introduce a gauge-invariant lower bound on the extremality by computing the smallest value that Booth and Fairhurst's extremality parameter can take for any scaling. Using this lower bound, we conclude that the common horizons are at least moderately close to extremal just after they appear. Finally, following Lovelace et al. (2008), we construct quasiequilibrium binary-black-hole initial data with "overspun" marginally trapped surfaces with $8πS > A$ and for which our lower bound on their Booth-Fairhurst extremality exceeds unity. These superextremal surfaces are always surrounded by marginally outer trapped surfaces (i.e., by apparent horizons) with $8πS<A$. The extremality lower bound on the enclosing apparent horizon is always less than unity but can exceed the value for an extremal Kerr black hole. (Abstract abbreviated.)

preprint2014arXiv

Improved methods for simulating nearly extremal binary black holes

Astrophysical black holes could be nearly extremal (that is, rotating nearly as fast as possible); therefore, nearly extremal black holes could be among the binaries that current and future gravitational-wave observatories will detect. Predicting the gravitational waves emitted by merging black holes requires numerical-relativity simulations, but these simulations are especially challenging when one or both holes have mass $m$ and spin $S$ exceeding the Bowen-York limit of $S/m^2=0.93$. We present improved methods that enable us to simulate merging, nearly extremal black holes more robustly and more efficiently. We use these methods to simulate an unequal-mass, precessing binary black hole coalescence, where the larger black hole has $S/m^2=0.99$. We also use these methods to simulate a non-precessing binary black hole coalescence, where both black holes have $S/m^2=0.994$, nearly reaching the Novikov-Thorne upper bound for holes spun up by thin accretion disks. We demonstrate numerical convergence and estimate the numerical errors of the waveforms; we compare numerical waveforms from our simulations with post-Newtonian and effective-one-body waveforms; we compare the evolution of the black-hole masses and spins with analytic predictions; and we explore the effect of increasing spin magnitude on the orbital dynamics (the so-called "orbital hangup" effect).

preprint2013arXiv

A catalog of 174 binary black-hole simulations for gravitational-wave astronomy

This paper presents a publicly available catalog of 174 numerical binary black-hole simulations following up to 35 orbits. The catalog includes 91 precessing binaries, mass ratios up to 8:1, orbital eccentricities from a few percent to $10^{-5}$, black-hole spins up to 98% of the theoretical maximum, and radiated energies up to 11.1% of the initial mass. We establish remarkably good agreement with post-Newtonian precession of orbital and spin directions for two new precessing simulations, and we discuss other applications of this catalog. Formidable challenges remain: e.g., precession complicates the connection of numerical and approximate analytical waveforms, and vast regions of the parameter space remain unexplored.