BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Prayush Kumar

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.CO

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

17 published item(s)

preprint2026arXiv

Rapid inference of gravitational-wave signals in the time domain using a heterodyned likelihood

Parameter estimation of gravitational wave signals is computationally intensive and typically requires millions of likelihood evaluations to construct posterior probability distributions. This computational cost increases significantly in the time domain, which requires non-diagonal covariance matrices to compute the likelihood. Consequently, parameter estimation of long-duration gravitational wave signals, such as binary neutron star mergers, becomes computationally infeasible in time domain. In this work, we detail a framework for the heterodyned likelihood that enables rapid inference in the time domain. Our method is applicable to signals with arbitrary mode content, and leverages the smoothness of the ratio of complex-valued waveform modes, approximating the ratio as a linear function within appropriately chosen time bins. This allows downsampling of the waveform modes and a reformulation of the likelihood, such that it depends only on the bin edges. We demonstrate that this likelihood recovers posteriors that are indistinguishable from those obtained using the standard likelihood in the time domain. We also observe dramatic improvement in speed - for a 128 seconds-long gravitational wave signal, our method is at least $\sim 400$ times faster than the standard time-domain analysis, reducing the wall clock time to just a few hours. We also demonstrate the reliability and unbiasedness of the likelihood using percentile-percentile tests for binary black hole and binary neutron star injections. We use the Gohberg-Semencul representation of the inverse of Toeplitz covariance matrix to accelerate matrix-vector products, which has potential applications even in non heterodyned time-domain inference.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network.

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.

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
preprint2020arXiv

Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Accuracy of binary black hole waveform models for aligned-spin binaries

Coalescing binary black holes are among the primary science targets for second generation ground-based gravitational wave (GW) detectors. Reliable GW models are central to detection of such systems and subsequent parameter estimation. This paper performs a comprehensive analysis of the accuracy of recent waveform models for binary black holes with aligned spins, utilizing a new set of $84$ high-accuracy numerical relativity simulations. Our analysis covers comparable mass binaries ($1\le m_1/m_2\le 3$), and samples independently both black hole spins up to dimensionless spin-magnitude of $0.9$ for equal-mass binaries and $0.85$ for unequal mass binaries. Furthermore, we focus on the high-mass regime (total mass $\gtrsim 50M_\odot$). The two most recent waveform models considered (PhenomD and SEOBNRv2) both perform very well for signal detection, losing less than 0.5\% of the recoverable signal-to-noise ratio $ρ$, except that SEOBNRv2's efficiency drops mildly for both black hole spins aligned with large magnitude. For parameter estimation, modeling inaccuracies of SEOBNRv2 are found to be smaller than systematic uncertainties for moderately strong GW events up to roughly $ρ\lesssim 15$. PhenomD's modeling errors are found to be smaller than SEOBNRv2's, and are generally irrelevant for $ρ\lesssim 20$. Both models' accuracy deteriorates with increased mass-ratio, and when at least one black hole spin is large and aligned. The SEOBNRv2 model shows a pronounced disagreement with the numerical relativity simulation in the merger phase, for unequal masses and simultaneously both black hole spins very large and aligned. Two older waveform models (PhenomC and SEOBNRv1) are found to be distinctly less accurate than the more recent PhenomD and SEOBNRv2 models. Finally, we quantify the bias expected from all GW models during parameter estimation for recovery of binary's masses and spins.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the accuracy and precision of numerical waveforms: Effect of waveform extraction methodology

We present a new set of 95 numerical relativity simulations of non-precessing binary black holes (BBHs). The simulations sample comprehensively both black-hole spins up to spin magnitude of 0.9, and cover mass ratios 1 to 3. The simulations cover on average 24 inspiral orbits, plus merger and ringdown, with low initial orbital eccentricities $e<10^{-4}$. A subset of the simulations extends the coverage of non-spinning BBHs up to mass ratio $q=10$. Gravitational waveforms at asymptotic infinity are computed with two independent techniques, extrapolation, and Cauchy characteristic extraction. An error analysis based on noise-weighted inner products is performed. We find that numerical truncation error, error due to gravitational wave extraction, and errors due to the finite length of the numerical waveforms are of similar magnitude, with gravitational wave extraction errors somewhat dominating at noise-weighted mismatches of $\sim 3\times 10^{-4}$. This set of waveforms will serve to validate and improve aligned-spin waveform models for gravitational wave science.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accuracy and precision of gravitational-wave models of inspiraling neutron star -- black hole binaries with spin: comparison with numerical relativity in the low-frequency regime

Coalescing binaries of neutron stars (NS) and black holes (BH) are one of the most important sources of gravitational waves for the upcoming network of ground based detectors. Detection and extraction of astrophysical information from gravitational-wave signals requires accurate waveform models. The Effective-One-Body and other phenomenological models interpolate between analytic results and $10-30$ orbit numerical relativity (NR) merger simulations. In this paper we study the accuracy of these models using new NR simulations that span $36-88$ orbits, with mass-ratios and black hole spins $(q,χ_{BH}) = (7, \pm 0.4), (7, \pm 0.6)$, and $(5, -0.9)$. We find that: (i) the recently published SEOBNRv1 and SEOBNRv2 models of the Effective-One-Body family disagree with each other (mismatches of a few percent) for black hole spins $\geq 0.5$ or $\leq -0.3$, with waveform mismatch accumulating during early inspiral; (ii) comparison with numerical waveforms indicate that this disagreement is due to phasing errors of SEOBNRv1, with SEOBNRv2 in good agreement with all of our simulations; (iii) Phenomenological waveforms disagree with SEOBNRv2 over most of the NSBH binary parameter space; (iv) comparison with NR waveforms shows that most of the model's dephasing accumulates near the frequency interval where it switches to a phenomenological phasing prescription; and finally (v) both SEOBNR and post-Newtonian (PN) models are effectual for NSBH systems, but PN waveforms will give a significant bias in parameter recovery. Our results suggest that future gravitational-wave detection searches and parameter estimation efforts targeted at NSBH systems with $q\lesssim 7$ and $χ_\mathrm{BH} \approx [-0.9, +0.6]$ will benefit from using SEOBNRv2 templates. For larger black hole spins and/or binary mass-ratios, we recommend the models be further investigated as suitable NR simulations become available.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accurate and efficient waveforms for compact binaries on eccentric orbits

Compact binaries that emit gravitational waves in the sensitivity band of ground-based detectors can have non-negligible eccentricities just prior to merger, depending on the formation scenario. We develop a purely analytic, frequency-domain model for gravitational waves emitted by compact binaries on orbits with small eccentricity, which reduces to the quasi-circular post-Newtonian approximant TaylorF2 at zero eccentricity and to the post-circular approximation of Yunes et al. (2009) at small eccentricity. Our model uses a spectral approximation to the (post-Newtonian) Kepler problem to model the orbital phase as a function of frequency, accounting for eccentricity effects up to ${\cal{O}}(e^8)$ at each post-Newtonian order. Our approach accurately reproduces an alternative time-domain eccentric waveform model for eccentricities $e\in [0, 0.4]$ and binaries with total mass less than 12 solar masses. As an application, we evaluate the signal amplitude that eccentric binaries produce in different networks of existing and forthcoming gravitational waves detectors. Assuming a population of eccentric systems containing black holes and neutron stars that is uniformly distributed in co-moving volume, we estimate that second generation detectors like Advanced LIGO could detect approximately 0.1-10 events per year out to redshift $z\sim 0.2$, while an array of Einstein Telescope detectors could detect hundreds of events per year to redshift $z \sim 2.3$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-forced evolutions of an implicit rotating source: A natural framework to model comparable and intermediate mass-ratio systems from inspiral through ringdown

We develop a waveform model to describe the inspiral, merger and ringdown of binary systems with comparable and intermediate mass-ratios. This model incorporates first-order conservative self-force corrections to the energy and angular momentum, which are valid in the strong-field regime [1]. We model the radiative part of the self-force by deriving second-order radiative corrections to the energy flux. These corrections are obtained by minimizing the phase discrepancy between our self-force model and the effective one body model [2, 3] for a variety of mass-ratios. We show that our model performs substantially better than post-Newtonian approximants currently used to model neutron star-black hole mergers from early inspiral to the innermost stable circular orbit. In order to match the late inspiral evolution onto the plunge regime, we extend the 'transition phase' developed by Ori and Thorne [4] by including finite mass-ratio corrections and modelling the orbital phase evolution using an implicit rotating source [5]. We explicitly show that the implicit rotating source approach provides a natural transition from late-time radiation to ringdown that is equivalent to ringdown waveform modelling based on a sum of quasinormal modes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Template Banks for Binary black hole searches with Numerical Relativity waveforms

Gravitational waves (GW) from coalescing stellar-mass black hole binaries (BBH) are expected to be detected by the Advanced Laser Interferometer Gravitational-wave Observatory and Advanced Virgo. Detection searches operate by matched-filtering the detector data using a bank of waveform templates. Traditionally, template banks for BBH are constructed from intermediary analytical waveform models which are calibrated against numerical relativity simulations and which can be aluated for any choice of BBH parameters. This paper explores an alternative to the traditional approach, namely the construction of template banks directly from numerical BBH simulations. Using non-spinning BBH systems as an example, we demonstrate which regions of the mass-parameter plane can be covered with existing numerical BBH waveforms. We estimate the required number and required length of BBH simulations to cover the entire non-spinning BBH parameter plane up to mass-ratio 10, thus illustrating that our approach can be used to guide parameter placement of future numerical simulations. We derive error bounds which are independent of analytical waveform models; therefore, our formalism can be used to independently test the accuracy of such waveform models. The resulting template banks are suitable for advanced LIGO searches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accurate modeling of intermediate-mass-ratio inspirals: Exploring the form of the self-force in the intermediate-mass-ratio regime

In this paper we develop a waveform model that accurately reproduces the dynamical evolution of intermediate-mass-ratio inspirals, as predicted by the effective-one-body (EOB) model introduced in [1], and which enables us to shed some light on the form of the self-force for events with mass-ratio 1:6, 1:10 and 1:100. To complement this study, we make use of self-force results in the extreme-mass-ratio regime, and of predictions of the EOB model introduced in [1], to derive a prescription for the shift of the orbital frequency at the innermost stable circular orbit which consistently captures predictions from the extreme, intermediate and comparable mass-ratio regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Search for the PN coefficients for the Energy flux through Gravitational Waves from Black-Hole Binaries using Markov Chain Monte Carlo

In this work, the focus is on the improvement of the existing post-Newtonian approximation for the gravitational flux from Super Massive Black Hole Binaries. In order to improve the existing templates for LISA, we need more accurate post-Newtonian expansions for the gravitational flux. Stochastic search techniques like the Markov Chain Monte Carlo (MCMC) have been used extensively for searching for sky parameters etc. The idea is to combine the two and approach the problem of finding post-Newtonian coefficients using MCMC. It has been shown that matching against a 5.5PN signal, with noise, the last coefficient can be found by MCMC very easily and displays fast convergence. Also the space for higher dimensional searches are explored.

0 citations0 reviewsNo code linkedOpen access