BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Jordan Moxon

Jordan Moxon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEhep-thastro-ph.COphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
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

14 published item(s)

preprint2024arXiv

Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling

We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to $5\%$) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counter-rotating and co-rotating quasinormal modes as black hole spin increases.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Numerical relativity surrogate model with memory effects and post-Newtonian hybridization

Numerical relativity simulations provide the most precise templates for the gravitational waves produced by binary black hole mergers. However, many of these simulations use an incomplete waveform extraction technique -- extrapolation -- that fails to capture important physics, such as gravitational memory effects. Cauchy-characteristic evolution (CCE), by contrast, is a much more physically accurate extraction procedure that fully evolves Einstein's equations to future null infinity and accurately captures the expected physics. In this work, we present a new surrogate model, NRHybSur3dq8$\_$CCE, built from CCE waveforms that have been mapped to the post-Newtonian (PN) BMS frame and then hybridized with PN and effective one-body (EOB) waveforms. This model is trained on 102 waveforms with mass ratios $q\leq8$ and aligned spins $χ_{1z}, \, χ_{2z} \in \left[-0.8, 0.8\right]$. The model spans the entire LIGO-Virgo-KAGRA (LVK) frequency band (with $f_{\text{low}}=20\text{Hz}$) for total masses $M\gtrsim2.25M_{\odot}$ and includes the $\ell\leq4$ and $(\ell,m)=(5,5)$ spin-weight $-2$ spherical harmonic modes, but not the $(3,1)$, $(4,2)$ or $(4,1)$ modes. We find that NRHybSur3dq8$\_$CCE can accurately reproduce the training waveforms with mismatches $\lesssim2\times10^{-4}$ for total masses $2.25M_{\odot}\leq M\leq300M_{\odot}$ and can, for a modest degree of extrapolation, capably model outside of its training region. Most importantly, unlike previous waveform models, the new surrogate model successfully captures memory effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code

Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple days to solve these problems at high resolution and when matter is involved, because they are either hard to parallelize or require a large amount of computational resources. Here we present a new solver for linear and nonlinear elliptic problems that is designed to scale with resolution and to parallelize on computing clusters. To achieve this we employ a discontinuous Galerkin discretization, an iterative multigrid-Schwarz preconditioned Newton-Krylov algorithm, and a task-based parallelism paradigm. To accelerate convergence of the elliptic solver we have developed novel subdomain-preconditioning techniques. We find that our multigrid-Schwarz preconditioned elliptic solves achieve iteration counts that are independent of resolution, and our task-based parallel programs scale over 200 million degrees of freedom to at least a few thousand cores. Our new code solves a classic initial data problem for binary black holes faster than the spectral code SpEC when distributed to only eight cores, and in a fraction of the time on more cores. It is publicly accessible in the next-generation SpECTRE numerical relativity code. Our results pave the way for highly parallel elliptic solves in numerical relativity and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fixing the BMS Frame of Numerical Relativity Waveforms

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE's Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during postprocessing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincaré charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincaré transformations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational-wave echoes from numerical-relativity waveforms via space-time construction near merging compact objects

We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbations of the final Kerr black hole, and a nonlinear region. At late times, the boundary between the two regions is an infalling shell. The BHP region contains late-time gravitational-waves emitted toward the future horizon, as well as those emitted toward future null infinity. In this region, by imposing no-ingoing wave conditions at past null infinity, and matching out-going waves at future null infinity with waveforms computed from numerical relativity, we can obtain waves that travel toward the future horizon. In particular, the Newman-Penrose $ψ_0$ associated with the in-going wave on the horizon is related to tidal deformations measured by fiducial observers floating above the horizon. We further determine the boundary of the BHP region on the future horizon by imposing that $ψ_0$ inside the BHP region can be faithfully represented by quasi-normal modes. Using a physically-motivated way to impose boundary conditions near the horizon, and applying the so-called Boltzmann reflectivity, we compute the quasi-normal modes of non-rotating ECOs, as well as gravitational-wave echoes. We also investigate the detectability of these echoes in current and future detectors, and prospects for parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High Precision Ringdown Modeling: Multimode Fits and BMS Frames

Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), an extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim10^5$ when using multimode fitting as opposed to only fitting the $(2,\pm2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms$-$which exhibit the memory effect$-$to a BMS frame known as the super rest frame, there is an improvement of $\sim10^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim4$ better than using the publicly-available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by NR and comparing our results to previous QNM studies.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adding Gravitational Memory to Waveform Catalogs using BMS Balance Laws

Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration's catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Lastly, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package $\texttt{sxs}$, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Computation of Displacement and Spin Gravitational Memory in Numerical Relativity

We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS Collaboration's $\texttt{SpEC}$ code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using $\texttt{SpECTRE}$'s version of a Cauchy-characteristic extraction. We find that we can accurately resolve the strain's traditional $m=0$ memory modes and some of the $m\not=0$ oscillatory memory modes that have previously only been theorized. We also perform a separate calculation of the memory using equations for the Bondi-Metzner-Sachs charges as well as the energy and angular momentum fluxes at asymptotic infinity. Our new calculation uses only the gravitational wave strain and two of the Weyl scalars at infinity. Also, this computation shows that the memory modes can be understood as a combination of a memory signal throughout the binary's inspiral and merger phases, and a quasinormal mode signal near the ringdown phase. Additionally, we find that the magnetic memory, up to numerical error, is indeed zero as previously conjectured. Lastly, we find that signal-to-noise ratios of memory for LIGO, the Einstein Telescope (ET), and the Laser Interferometer Space Antenna (LISA) with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or Minimal Waveform models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extending Gravitational Wave Extraction Using Weyl Characteristic Fields

We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in postprocessing. These methods have been used to produce $Ψ_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes (SXS) waveform catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy characteristic extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dissipation in extreme-mass ratio binaries with a spinning secondary

We present the gravitational-wave flux balance law in an extreme mass-ratio binary with a spinning secondary. This law relates the flux of energy (angular momentum) radiated to null infinity and through the event horizon to the local change in the secondary's orbital energy (angular momentum) for generic (non-resonant) bound orbits in Kerr spacetime. As an explicit example we compute these quantities for a spin-aligned body moving on a circular orbit around a Schwarzschild black hole. We perform this calculation both analytically, via a high-order post-Newtonian expansion, and numerically in two different gauges. Using these results we demonstrate explicitly that our new balance law holds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-accuracy waveforms for black hole-neutron star systems with spinning black holes

The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-precessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to non-spinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy Characteristic Extraction. The new waveforms have $<0.1\,{\rm rad}$ phase errors during inspiral, rising to $\sim (0.2-0.4)\,{\rm rad}$ errors at merger, and $\lesssim 1\%$ error in their amplitude. We compute the faithfulness of recent analytical models to these numerical results, and find that models specifically designed for BHNS binaries perform well ($F>0.99$) for binaries seen face-on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical relativity simulation of GW150914 beyond general relativity

We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $\gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spectral Cauchy-characteristic extraction of the gravitational wave news function

We present an improved spectral algorithm for Cauchy-characteristic extraction and characteristic evolution of gravitational waves in numerical relativity. The new algorithms improve spectral convergence both at the poles of the spherical-polar grid and at future null infinity, as well as increase the temporal resolution of the code. The key to the success of these algorithms is a new set of high-accuracy tests, which we present here. We demonstrate the accuracy of the code and compare with the existing PittNull implementation.

0 citations0 reviewsNo code linkedOpen access