Search NASA⌕ Search

SEARCH · Search NASA

Results for “gravitational waveforms”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Error quantification and comparison of binary neutron star gravitational waveforms from numerical relativity codes

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (fil, spec), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. Finally, we also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors.

gravitational waves↗
DOE OSTI2024

Gravitational waveforms: A tale of two formalisms

We revisit the quantum-amplitude-based derivation of the gravitational waveform emitted by the scattering of two spinless massive bodies at the third order in Newton’s constant, ℎ ∼ 𝐺 + 𝐺 2 + 𝐺 3 (one-loop level), and correspondingly update its comparison with its classically derived multipolar-post-Minkowskian counterpart. A spurious-pole-free reorganization of the one-loop five-point amplitude substantially simplifies the post-Newtonian expansion. We find complete agreement between the two results up to the fifth order in the small velocity expansion after taking into account three subtle aspects of the amplitude derivation: (1) in agreement with [A. Georgoudis et al., J. High Energy Phys. 03 (2024) 089], the term quadratic in the amplitude in the observable-based formalism [D. A. Kosower et al., J. High Energy Phys. 02 (2019) 137] generates a frame rotation by half the classical scattering angle; (2) the dimensional regularization of the infrared divergences of the amplitude introduces an additional (𝑑 − 4)/(𝑑 − 4) finite term; and (3) zero-frequency gravitons are found to contribute additional terms both at order ℎ ∼ 𝐺 1 and at order ℎ ∼ 𝐺 3 when including disconnected diagrams in the observable-based formalism.

Astronomy & Astrophysics↗
DOE OSTI2024

Testing Gravity with Realistic Gravitational Waveforms in Pulsar Timing Arrays

We consider the effects of relaxing the assumption that gravitational waves composing the stochastic gravitational wave background (SGWB) are uncorrelated between frequencies in analyses of the data from Pulsar Timing Arrays (PTAs). While individual monochromatic plane waves are often a good approximation, a background composed of astrophysical sources cannot be monochromatic since an infinite plane wave carries no signal. We consider how relaxing this assumption allows us to extract potential information about modified dispersion relations and other fundamental physics questions, as both the group and phase velocity of waves become relevant. After developing the formalism we carry out simple Gaussian wavepacket examples and then consider more realistic waveforms, such as that from binary inspirals. When the frequency evolves only slowly across the PTA temporal baseline, the monochromatic assumption at an effective mean frequency remains a good approximation and we provide scaling relations that characterize its accuracy.

Hu, Wayne↗
DOE OSTI2023

Numerical relativity higher order gravitational waveforms of eccentric, spinning, nonprecessing binary black hole mergers

We use the open source, community-driven, numerical relativity software, the EINSTEIN TOOLKIT to study the physics of eccentric, spinning, nonprecessing binary black hole mergers with mass-ratios q = {2; 4; 6}, individual dimensionless spin parameters chi 1z = +/- 0.6, chi 2z = +/- 0.3, that include higher order gravitational wave modes l <= 4, except for memory modes. Assuming stellar mass binary black hole mergers that may be detectable by the advanced LIGO detectors, we find that including modes up to l = 4 increases the signal-to-noise of compact binaries between 3.5% to 35%, compared to signals that only include the l = ImI = 2 mode. We use two waveform models, TEOBResumS and SEOBNRE, which incorporate spin and eccentricity corrections in the waveform dynamics, to quantify the orbital eccentricity of our numerical relativity catalog in a gauge-invariant manner through fitting factor calculations. Our findings indicate that the inclusion of higher order wave modes has a measurable effect in the recovery of moderately and highly eccentric black hole mergers, and thus it is essential to develop waveform models and signal processing tools that accurately describe the physics of these astrophysical sources.

Joshi, Abhishek V.↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Adaptive mesh refinement in binary black holes simulations

Abstract We discuss refinement criteria for the Berger–Rigoutsos (block-based) refinement algorithm in our numerical relativity code GR-Athena++ in the context of binary black hole (BBH) merger simulations. We compare three different strategies: the ‘box-in-box’ approach, the ‘sphere-in-sphere’ approach and a local criterion for refinement based on the estimation of truncation error of the finite difference scheme. We extract and compare gravitational waveforms using the three different mesh refinement methods and compare their accuracy against a calibration waveform and demonstrate that the sphere-in-sphere approach provides the best strategy overall when considering computational cost and the waveform accuracy. Ultimately, we demonstrate the capability of each mesh refinement method in accurately simulating gravitational waves from BBH systems—a crucial aspect for their application in next-generation detectors. We quantify the mismatch achievable with the different strategies by extrapolating the gravitational wave mismatch to higher resolution.

Astronomy & Astrophysics↗
DOE OSTI2022

Orbital Decay of Double White Dwarfs: Beyond Gravitational-wave Radiation Effects

The traditional description of the orbital evolution of compact-object binaries, like double white dwarfs (DWDs), assumes that the system is driven only by gravitational-wave (GW) radiation. However, the high magnetic fields with intensities of up to gigagausses measured in WDs alert a potential role of the electromagnetic (EM) emission in the evolution of DWDs. We evaluate the orbital dynamics of DWDs under the effects of GW radiation, tidal synchronization, and EM emission by a unipolar inductor generated by the magnetic primary and the relative motion of the nonmagnetic secondary. We show that the EM emission can affect the orbital dynamics for magnetic fields larger than megagausses. We applied the model to two known DWDs, SDSS J0651+2844 and ZTF J1539+5027, for which the GW radiation alone does not fully account for the measured orbital decay rate. We obtain upper limits to the primary's magnetic field strength, over which the EM emission causes an orbital decay faster than observed. The contribution of tidal locking and the EM emission is comparable, and together they can contribute up to 20% to the measured orbital decay rate. We show that the gravitational waveform for a DWD modeled as purely driven by GWs and including tidal interactions and EM emission can have large relative dephasing detectable in the mHz regime of frequencies relevant for space-based detectors like LISA. Therefore, including physics besides GW radiation in the waveform templates is essential to calibrate the GW detectors using known sources, e.g., ZTF J1539+5027, and to infer binary parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Supertranslations from Scattering Amplitudes

On shell methods have found a new application to local observables such as asymptotic radiation fields and gravitational waveforms. While these observables are invariant under small gauge transformations, they are known to depend on a choice of asymptotic gauge; in gravity on asymptotically Minkowski spacetimes, this is a choice of Bondi-Van der Burg-Metzner-Sachs frame. In this Letter, we provide a method for capturing these supertranslations, to all orders in perturbation theory, using the on shell framework of scattering amplitudes.

Classical black holes↗
DOE OSTI2024

Leading-order gravitational radiation to all spin orders

Starting with on-shell amplitudes compatible with the scattering of Kerr black holes up to Compton-amplitude contact terms, we produce the gravitational waveform and memory effect including spin at their leading post-Minkowskian orders to all orders in the spins of both scattering objects. For the memory effect, we present results at next-to-leading order as well, finding a closed form for all spin orders when the spins are anti-aligned and equal in magnitude. Considering instead generically oriented spins, we produce the next-to-leading-order memory to sixth order in spin. Compton-amplitude contact terms up to sixth order in spin are included throughout our analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational Waves from Neutrino Emission Asymmetries in Core-collapse Supernovae

We present a broadband spectrum of gravitational waves (GWs) from core-collapse supernovae (CCSNe) sourced by neutrino emission asymmetries for a series of full 3D simulations. The associated GW strain probes the long-term secular evolution of CCSNe and small-scale turbulent activity and provides insight into the geometry of the explosion. For nonexploding models, both the neutrino luminosity and the neutrino gravitational waveform will encode information about the spiral SASI. The neutrino memory will be detectable for a wide range of progenitor masses for a galactic event. Our results can be used to guide near-future decihertz and long-baseline GW detection programs, including aLIGO, the Einstein Telescope, and DECIGO.

Astronomy & Astrophysics↗
DOE OSTI2023

New Covariant Density Functionals of Nuclear Matter for Compact Star Simulations

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of the not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gravitational wave propagation beyond general relativity: waveform distortions and echoes

We study the cosmological propagation of gravitational waves (GWs) beyond general relativity (GR) across homogeneous and isotropic backgrounds. We consider scenarios in which GWs interact with an additional tensor field and use a parametrized phenomenological approach that generically describes their coupled equations of motion. We analyze four distinct classes of derivative and non-derivative interactions: mass, friction, velocity, and chiral. We apply the WKB formalism to account for the cosmological evolution and obtain analytical solutions to these equations. We corroborate these results by analyzing numerically the propagation of a toy GW signal. We then proceed to use the analytical results to study the modified propagation of realistic GWs from merging compact binaries, assuming that the GW signal emitted is the same as in GR. We generically find that tensor interactions lead to copies of the originally emitted GW signal, each one with its own possibly modified dispersion relation. These copies can travel coherently and interfere with each other leading to a scrambled GW signal, or propagate decoherently and lead to echoes arriving at different times at the observer that could be misidentified as independent GW events. Depending on the type of tensor interaction, the detected GW signal may exhibit amplitude and phase distortions with respect to a GW waveform in GR, as well as birefringence effects. Furthermore, we discuss observational probes of these tensor interactions with both individual GW events, as well as population studies for both ground- and space-based detectors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Probing compactified extra dimensions with gravitational waves

We study the effect of compact extra dimensions on the gravitational wave luminosity and waveform. We consider a toy model, with a compactified fifth dimension, and matter confined on a brane. We work in the context of five-dimensional (5D) general relativity, though we do make connections with the corresponding Kaluza-Klein effective 4D theory. We show that the luminosity of gravitational waves emitted in 5D gravity by a binary with the same characteristics (same masses and separation distance) as a 4D binary is 20.8% less relative to the 4D case, to leading post-Newtonian order. The phase of the gravitational waveform differs by 26% relative to the 4D case, to leading post-Newtonian order. Such a correction arises mainly due to the coupling between matter and dilaton field in the effective 4D picture and agrees with previous calculations when we set black holes’ scalar charges to be those computed from the Kaluza-Klein reduction. The above corrections to the waveform and the luminosity are inconsistent with the gravitational-wave and binary pulsar observations, and they thus effectively rule out the possibility of such a simple compactified higher dimensions scenario. We also comment on how our results change if there are several compactified extra dimensions and show that the discrepancy with 4D general relativity only increases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

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. Here, 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.1rad phase errors during inspiral, rising to ~(0.2-0.4)rad errors at merger, and ≲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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes ($\ell$, |m|) = {(2,2), (2,1), (3,3), (3,2), (4,4)}, and mode mixing effects in the $\ell$ = 3, |m| = 2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Statistical and systematic uncertainties in extracting the source properties of neutron star-black hole binaries with gravitational waves

Gravitational waves emitted by neutron star black hole mergers encode key properties of neutron stars - such as their size, maximum mass and spins - and black holes. However, the presence of matter and the high mass ratio makes generating long and accurate waveforms from these systems hard to do with numerical relativity, and not much is known about systematic uncertainties due to waveform modeling. In this work, we simulate gravitational waves from neutron star black hole mergers by hybridizing numerical relativity waveforms produced with the SpEC code with a recent numerical relativity surrogate NRHybSur3dq8Tidal. These signals are analyzed using a range of available waveform families, and statistical and systematic errors are reported. We find that at a network signal-to-noise ratio (SNR) of 30, statistical uncertainties are usually larger than systematic offsets, while at an SNR of 70 the two become comparable. The individual black hole and neutron star masses, as well as the mass ratios, are typically measured very precisely, though not always accurately at high SNR. At a SNR of 30 the neutron star tidal deformability can only be bound from above, while for louder sources it can be measured and constrained away from zero. All neutron stars in our simulations are non-spinning, but in no case we can constrain the neutron star spin to be smaller than ~0.4 (90% credible interval). Waveform families whose late inspiral has been tuned specifically for neutron star black hole signals typically yield the most accurate characterization of the source parameters. Their measurements are in tension with those obtained using waveform families tuned against binary neutron stars, even for mass ratios that could be relevant for both binary neutron stars and neutron star black holes mergers.

79 ASTRONOMY AND ASTROPHYSICS↗