Researcher profile

Vijay Varma

Vijay Varma contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
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

13 published item(s)

preprint2026arXiv

Analysis of GWTC-3 with fully precessing numerical relativity surrogate models

The third Gravitational-Wave Transient Catalog (GWTC-3) contains 90 binary coalescence candidates detected by the LIGO-Virgo-KAGRA Collaboration (LVK). We provide a re-analysis of binary black hole (BBH) events using a recently developed numerical relativity (NR) waveform surrogate model, NRSur7dq4, that includes all $\ell \leq 4$ spin-weighted spherical harmonic modes as well as the complete physical effects of precession. Properties of the remnant black holes' (BH's) mass, spin vector, and kick vector are found using an associated remnant surrogate model NRSur7dq4Remnant. Both NRSur7dq4 and NRSur7dq4Remnant models have errors comparable to numerical relativity simulations and allow for high-accuracy parameter estimates. We restrict our analysis to 47 BBH events that fall within the regime of validity of NRSur7dq4 (mass ratios greater than 1/6 and total masses greater than $60 M_{\odot}$). While for most of these events our results match the LVK analyses that were obtained using the semi-analytical models such as IMRPhenomXPHM and SEOBNRv4PHM, we find that for more than 20\% of events the NRSur7dq4 model recovers noticeably different measurements of black hole properties like the masses and spins, as well as extrinsic properties like the binary inclination and distance. For instance, GW150914_095045 exhibits noticeable differences in spin precession and spin magnitude measurements. Other notable findings include one event (GW191109_010717) that constrains the effective spin $χ_{eff}$ to be negative at a 99.3\% credible level and two events (GW191109_010717 and GW200129_065458) with well-constrained kick velocities. Furthermore, compared to the models used in the LVK analyses, NRSur7dq4 recovers a larger signal-to-noise ratio and/or Bayes factors for several events.

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.

preprint2022arXiv

Evidence of large recoil velocity from a black hole merger signal

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.

preprint2022arXiv

Hints of spin-orbit resonances in the binary black hole population

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, $ϕ_1$ and $ϕ_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $ϕ_1$ and $ϕ_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for $Δϕ= ϕ_1 - ϕ_2 \sim \pm π$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $ϕ_1 \sim -π/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the $ϕ_1$ and $Δϕ$ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available.

preprint2022arXiv

Measuring binary black hole orbital-plane spin orientations

Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $ϕ_1$, $ϕ_2$ and $Δϕ=ϕ_1-ϕ_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed dimensionless time or frequency near the merger. This leads to significantly improved measurements for $ϕ_1$ and $ϕ_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $Δϕ$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$.

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.

preprint2021arXiv

Up-down instability of binary black holes in numerical relativity

Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechanism for the formation of precessing binaries in environments where sources are preferentially formed with (anti) aligned spins. In this paper, we present the first full numerical relativity simulations capturing this instability. These simulations span $\sim 100$ orbits and $\sim 3$-$5$ precession cycles before merger, making them some of the longest numerical relativity simulations to date. Initialized with a small perturbation of $1^{\circ}$-$10^{\circ}$, the instability causes a dramatic growth of the spin misalignments, which can reach $\sim 90^{\circ}$ near merger. We show that this leaves a strong imprint on the subdominant modes of the gravitational wave signal, which can potentially be used to distinguish up-down binaries from other sources. Finally, we show that post-Newtonian and effective-one-body approximants are able to reproduce the unstable dynamics of up-down binaries extracted from numerical relativity.

preprint2020arXiv

A Surrogate Model for Gravitational Wave Signals from Comparable- to Large- Mass-Ratio Black Hole Binaries

Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.

preprint2020arXiv

Extracting the Gravitational Recoil from Black Hole Merger Signals

Gravitational waves carry energy, angular momentum, and linear momentum. In generic binary black hole mergers, the loss of linear momentum imparts a recoil velocity, or a "kick", to the remnant black hole. We exploit recent advances in gravitational waveform and remnant black hole modeling to extract information about the kick from the gravitational wave signal. Kick measurements such as these are astrophysically valuable, enabling independent constraints on the rate of second-generation mergers. Further, we show that kicks must be factored into future ringdown tests of general relativity with third-generation gravitational wave detectors to avoid systematic biases. We find that, although little information can be gained about the kick for existing gravitational wave events, interesting measurements will soon become possible as detectors improve. We show that, once LIGO and Virgo reach their design sensitivities, we will reliably extract the kick velocity for generically precessing binaries--including the so-called superkicks, reaching up to 5000 km/s.

preprint2020arXiv

Gravitational waveforms of binary neutron star inspirals using post-Newtonian tidal splicing

The tidal deformations of neutron stars within an inspiraling compact binary alter the orbital dynamics, imprinting a signature on the gravitational wave signal. Modeling this signal could be done with numerical-relativity simulations, but these are too computationally expensive for many applications. Analytic post-Newtonian treatments are limited by unknown higher-order nontidal terms. This paper further builds upon the "tidal splicing" model in which post-Newtonian tidal terms are "spliced" onto numerical relativity simulations of black-hole binaries. We improve on previous treatments of tidal splicing by including spherical harmonic modes beyond the (2,2) mode, expanding the post-Newtonian expressions for tidal effects to 2.5 order, including dynamical tide corrections, and adding a partial treatment of the spin-tidal dynamics. Furthermore, instead of numerical relativity simulations, we use the spin-aligned binary black hole (BBH) surrogate model "NRHybSur3dq8" to provide the BBH waveforms that are input into the tidal slicing procedure. This allows us to construct spin-aligned, inspiraling TaylorT2 and TaylorT4 splicing waveform models that can be evaluated quickly. These models are tested against existing binary neutron star and black hole-neutron star simulations. We implement the TaylorT2 splicing model as an extension to "NRHybSur3dq8," creating a model that we call "NRHybSur3dq8Tidal."

preprint2020arXiv

Impact of subdominant modes on the interpretation of gravitational-wave signals from heavy binary black hole systems

Over the past year, a handful of new gravitational wave models have been developed to include multiple harmonic modes thereby enabling for the first time fully Bayesian inference studies including higher modes to be performed. Using one recently-developed numerical relativity surrogate model, NRHybSur3dq8, we investigate the importance of higher modes on parameter inference of coalescing massive binary black holes. We focus on examples relevant to the current three-detector network of observatories, with a detector-frame mass set to $120 M_\odot$ and with signal amplitude values that are consistent with plausible candidates for the next few observing runs. We show that for such systems the higher mode content will be important for interpreting coalescing binary black holes, reducing systematic bias, and computing properties of the remnant object. Even for comparable-mass binaries and at low signal amplitude, the omission of higher modes can influence posterior probability distributions. We discuss the impact of our results on source population inference and self-consistency tests of general relativity. Our work can be used to better understand asymmetric binary black hole merger events, such as GW190412. Higher modes are critical for such systems, and their omission usually produces substantial parameter biases.

preprint2019arXiv

Testing the"no-hair" nature of binary black holes using the consistency of multipolar gravitational radiation

Gravitational-wave (GW) observations of binary black holes offer the best probes of the relativistic, strong-field regime of gravity. Gravitational radiation, in the leading order is quadrupolar. However, non-quadrupole (higher-order) modes make appreciable contribution to the radiation from binary black holes with large mass ratios and misaligned spins. The multipolar structure of the radiation is fully determined by the intrinsic parameters (masses and spin angular momenta of the companion black holes) of a binary in quasi-circular orbit. Following our previous work \cite{Dhanpal:2018ufk}, we develop multiple ways of testing the consistency of the observed GW signal with the expected multipolar structure of radiation from binary black holes in general relativity. We call this a "no-hair" test of binary black holes as this is similar to testing the "no-hair" theorem for isolated black holes through mutual consistency of the quasi-normal mode spectrum. We use Bayesian inference on simulated GW signals that are consistent/inconsistent with binary black holes in GR to demonstrate the power of the proposed tests. We also make estimate systematic errors arising as a result of neglecting companion spins.

preprint2019arXiv

The binary black hole explorer: on-the-fly visualizations of precessing binary black holes

Binary black hole mergers are of great interest to the astrophysics community, not least because of their promise to test general relativity in the highly dynamic, strong field regime. Detections of gravitational waves from these sources by LIGO and Virgo have garnered widespread media and public attention. Among these sources, precessing systems (with misaligned black-hole spin/orbital angular momentum) are of particular interest because of the rich dynamics they offer. However, these systems are, in turn, more complex compared to nonprecessing systems, making them harder to model or develop intuition about. Visualizations of numerical simulations of precessing systems provide a means to understand and gain insights about these systems. However, since these simulations are very expensive, they can only be performed at a small number of points in parameter space. We present binaryBHexp, a tool that makes use of surrogate models of numerical simulations to generate on-the-fly interactive visualizations of precessing binary black holes. These visualizations can be generated in a few seconds, and at any point in the 7-dimensional parameter space of the underlying surrogate models. With illustrative examples, we demonstrate how this tool can be used to learn about precessing binary black hole systems.