Researcher profile

Sascha Husa

Sascha Husa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

preprint2025arXiv

Eccentric or circular? A reanalysis of binary black hole gravitational wave events for orbital eccentricity signatures

We present a reanalysis of 17 gravitational-wave events detected with Advanced LIGO and Advanced Virgo in their first three observing runs, using the new IMRPhenomTEHM model -- a phenomenological time-domain multipolar waveform model for aligned-spin black-hole binaries in elliptical orbits with two eccentric parameters: eccentricity and mean anomaly. We also analyze all events with the underlying quasi-circular model IMRPhenomTHM to study the impact of including eccentricity and compare the eccentric and quasi-circular binary hypotheses. The high computational efficiency of IMRPhenomTEHM enables us to explore the impact of two different eccentricity priors -- uniform and log-uniform -- as well as different sampler and data settings. We find evidence for eccentricity in two publicly available LVK events, GW200129 and GW200208_22, with Bayes factors favoring the eccentric hypothesis over the quasi-circular aligned-spin scenario: $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[1.30^{+0.15}_{-0.15}, 5.14^{+0.15}_{-0.15}\right]$ and $\log_{10}\mathcal{B}_{\mathrm{E/QC}}\in\left[0.49^{+0.08}_{-0.08}, 1.14^{+0.08}_{-0.08}\right]$, respectively. Additionally, the two high-mass events GW190701 and GW190929 exhibit potential eccentric features. For all four events, we conduct further analyses to study the impact of different sampler settings. We also investigate waveform systematics by exploring the support for spin precession using IMRPhenomTPHM and NRSur7dq4, offering new insights into the formation channels of detected binaries. Our results highlight the importance of considering eccentric waveform models in future observing runs, alongside precessing models, as they can help mitigate potential biases in parameter estimation studies. This will be particularly relevant with the expected increase in the diversity of the binary black hole population with new detectors.

preprint2025arXiv

First eccentric inspiral-merger-ringdown analysis of neutron star-black hole mergers

The gravitational wave event GW200105 was the first confident neutron star-black hole (NSBH) merger identified by the LIGO-Virgo-KAGRA collaboration. A recent analysis by Morras et al. with an eccentric precessing waveform model that describes the inspiral phase of the $l=2$ and $m=\{0,\pm 2\}$ modes has identified this event as the first NSBH merger with strong evidence of orbital eccentricity. In this paper we perform the first analysis of this event with an aligned-spin eccentric waveform model that describes the full inspiral, merger, and ringdown, includes subdominant harmonics, and is partially calibrated to numerical relativity simulations. This analysis confirms the results and finds evidence in favor of eccentricity even with a log-uniform prior in eccentricity. We also analyze the NSBH events GW200115 and GW230529, completing the analysis of all NSBHs with IMRPhenomTEHM, and find that these signal are consistent with vanishing eccentricity. Finally, we briefly discuss computational challenges when performing the analysis with time-domain eccentric waveform models.

preprint2025arXiv

Hybrid waveforms for precessing quasi-circular binary systems

The demand for long and accurate gravitational waveforms is increasing as we prepare for the next generation of detectors and seek to improve current waveform models. However, numerical relativity waveforms, while highly accurate, are often too short for these applications due to their high computational cost. Hybrid waveforms, which stitch together gravitational wave signals from different modeling approaches, provide a way to generate complete inspiral-merger-ringdown signals. While hybridization is well-established for aligned-spin systems, precession introduces additional complexities due to gauge ambiguities, frame dependence, or spin dynamics. Here we study the challenges associated with alignment of precessing waveforms and present a systematic approach for constructing hybrid waveforms of precessing quasi-circular systems. Our approach relies on minimal assumptions about the merger waveforms and employs the quadrupole-aligned frame to mitigate mode-mixing. Our method is designed to be robust and broadly applicable, imposing minimal constraints on the input waveforms. This framework expands the applicability of hybridization techniques, facilitating flexible hybrid construction for parameter estimation, model calibration, and gravitational-wave data analysis.

preprint2025arXiv

Time-domain phenomenological multipolar waveforms for aligned-spin binary black holes in elliptical orbits

We introduce IMRPhenomTEHM, a new phenomenological time-domain model for eccentric aligned-spin binary black holes. Building upon the accurate quasi-circular IMRPhenomTHM model, IMRPhenomTEHM integrates the eccentric post-Newtonian (PN) dynamics and introduces eccentric corrections into the waveform multipoles up to 3PN, including spin effects. The model incorporates the dominant (2, $\pm$2) spherical harmonic mode, as well as the subdominant modes (2, $\pm$1), (3, $\pm$3), (4, $\pm$4), and (5, $\pm$5), assuming the binary has circularized by the time of merger. This approach ensures a smooth transition to the non-eccentric limit, providing an accurate quasi-circular limit against the IMRPhenomTHM model. When comparing against 28 public eccentric numerical relativity simulations from the Simulating eXtreme Spacetimes catalog, IMRPhenomTEHM achieves lower than 2% unfaithfulness, confirming its accurate description without calibration to numerical relativity eccentric datasets. IMRPhenomTEHM provides a reliable description of the evolution of eccentric black hole binaries with aligned spins and eccentricities lower than $e=0.4$ at a frequency of 10 Hz, making it suitable for upcoming gravitational-wave observing runs. We validate the model's accuracy through parameter estimation studies, recovering injected parameters within 90% credible intervals for three numerical relativity eccentric simulations and reanalyzing GW150914 and GW190521, obtaining results consistent with the literature.

preprint2021arXiv

A detailed analysis of GW190521 with phenomenological waveform models

In this paper we present an extensive analysis of the GW190521 gravitational wave event with the current (fourth) generation of phenomenological waveform models for binary black hole coalescences. GW190521 stands out from other events since only a few wave cycles are observable. This leads to a number of challenges, one being that such short signals are prone to not resolve approximate waveform degeneracies, which may result in multi-modal posterior distributions. The family of waveform models we use includes a new fast time-domain model IMRPhenomTPHM, which allows us extensive tests of different priors and robustness with respect to variations in the waveform model, including the content of spherical harmonic modes. We clarify some issues raised in a recent paper [Nitz&Capano], associated with possible support for a high-mass ratio source, but confirm their finding of a multi-modal posterior distribution, albeit with important differences in the statistical significance of the peaks. In particular, we find that the support for both masses being outside the PISN mass-gap, and the support for an intermediate mass ratio binary are drastically reduced with respect to what Nitz&Capano found. We also provide updated probabilities for associating GW190521 to the potential electromagnetic counterpart from ZTF.

preprint2021arXiv

New twists in compact binary waveform modelling: a fast time domain model for precession

We present IMRPhenomTPHM, a phenomenological model for the gravitational wave signals emitted by the coalescence of quasi-circular precessing binary black holes systems. The model is based on the "twisting up" approximation, which maps non-precessing signals to precessing ones in terms of a time dependent rotation described by three Euler angles, and which has been utilized in several frequency domain waveform models that have become standard tools in gravitational wave data analysis. Our model is however constructed in the time domain, which allows several improvements over the frequency domain models: we do not use the stationary phase approximation, we employ a simple approximation for the precessing Euler angles for the ringdown signal, and we implement a new method for computing the Euler angles through the evolution of the spin dynamics of the system, which is more accurate and also computationally efficient.

preprint2020arXiv

A first survey of spinning eccentric black hole mergers: Numerical relativity simulations, hybrid waveforms, and parameter estimation

We analyze a new numerical relativity dataset of spinning but nonprecessing binary black holes on eccentric orbits, with eccentricities from approximately $0.1$ to $0.5$, with dimensionless spins up to $0.75$ included at mass ratios $q=m_1/m_2 = (1, 2)$, and further nonspinning binaries at mass ratios $q = (1.5, 3, 4)$. A comparison of the final mass and spin of these simulations with noneccentric data extends previous results in the literature on circularization of eccentric binaries to the spinning case. For the $(l,m)=(2,2)$ spin-weighted spherical harmonic mode we construct eccentric hybrid waveforms that connect the numerical relativity data to a post-Newtonian description for the inspiral, and we discuss the limitations in the current knowledge about post-Newtonian theory which complicate the generation of eccentric hybrid waveforms. We also perform a Bayesian parameter estimation study, quantifying the parameter biases introduced when using three different quasicircular waveform models to estimate the parameters of highly eccentric binary systems. We find that the used aligned-spin quasicircular model including higher order modes produces lower bias in certain parameters than the nonprecessing quasicircular model without higher order modes and the quasicircular precessing model.

preprint2020arXiv

Accelerating the evaluation of inspiral-merger-ringdown waveforms with adapted grids

This paper presents an algorithm to accelerate the evaluation of inspiral-merger-ringdown waveform models for gravitational wave data analysis. While the idea can also be applied in the time domain, here we focus on the frequency domain, which is most typically used to reduced computational cost in gravitational wave data analysis. Our work extends the idea of multibanding, which has been developed to accelerate frequency domain waveforms, to include the merger and ringdown and spherical harmonics beyond the dominant quadrupole spherical harmonic. The original method is based on a heuristic algorithm based on the inspiral to de-refine the equi-spaced frequency grid used for data analysis where a coarser grid is sufficient for accurate evaluation of a waveform model. Here we use a different criterion, based on the local interpolation error, which is more flexible and can easily be adapted to general waveforms, if their phenomenology is understood. We discuss our implementation in the LIGO Algorithm Library for the PhenomXHM frequency domain model, and report the acceleration in different parts of the parameter space of compact binary systems.

preprint2020arXiv

Impact of eccentricity on the gravitational wave searches for binary black holes: High mass case

The possible formation of stellar-mass binary black holes through dynamical interactions in dense stellar environments predicts the existence of binaries with non-negligible eccentricity in the frequency band of ground-based gravitational wave detectors; the detection of binary black hole mergers with measurable orbital eccentricity would validate the existence of this formation channel. Waveform templates currently used in the matched-filter gravitational-wave searches of LIGO-Virgo data neglect effects of eccentricity which is expected to reduce their efficiency to detect eccentric binary black holes. Meanwhile, the sensitivity of coherent unmodeled gravitational-wave searches (with minimal assumptions about the signal model) have been shown to be largely unaffected by the presence of even sizable orbital eccentricity. In this paper, we compare the performance of two state-of-the-art search algorithms recently used by LIGO and Virgo to search for binary black holes in the second Observing Run (O2), quantifying their search sensitivity by injecting numerical-relativity simulations of inspiral-merger-ringdown eccentric waveforms into O2 LIGO data. Our results show that the matched-filter search PyCBC performs better than the unmodeled search cWB for the high chirp mass ($>20 M_{\odot}$) and low eccentricity region ($e_{30 Hz} < 0.3$) of parameter space. For moderate eccentricities and low chirp mass, on the other hand, the unmodeled search is more sensitive than the modeled search.

preprint2020arXiv

Setting the cornerstone for the IMRPhenomX family of models for gravitational waves from compact binaries: The dominant harmonic for non-precessing quasi-circular black holes

In this paper we present IMRPhenomXAS, a thorough overhaul of the IMRPhenomD [1,2] waveform model, which describes the dominant $l=2, \:| m | = 2$ spherical harmonic mode of non-precessing coalescing black holes in terms of piecewise closed form expressions in the frequency domain. Improvements include in particular the accurate treatment of unequal spin effects, and the inclusion of extreme mass ratio waveforms. IMRPhenomD has previously been extended to approximately include spin precession [3] and subdominant spherical harmonics [4], and with its extensions it has become a standard tool in gravitational wave parameter estimation. Improved extensions of IMRPhenomXAS are discussed in companion papers [5,6].

preprint2020arXiv

Simple procedures to reduce eccentricity of binary black hole simulations

We present simple procedures to construct quasi-circular initial data for numerical evolutions of binary black hole spacetimes. Our method consists of using Post-Newtonian theory in three ways: first to provide an initial guess for the initial momenta at 3.5PN order that implies low residual eccentricity, second to measure the resulting eccentricity, and third to calculate corrections to the momenta or initial separation which further reduce the eccentricity. Regarding the initial guess, we compare numerical evolutions in post-Newtonian theory to the post-circular and post-post-circular analytical approximations to quasi-circular data. We discuss a robust fitting procedure to measure eccentricity from numerical simulations using the orbital frequency, and derive from the quasi-Keplerian parametrization at 1PN oder the correction factors for the tangential and radial momentum components required to achieve reduce the measured eccentricity to zero. We first test our procedure integrating PN equations of motion at $3.5$PN where low eccentric initial data is easily obtained, and then apply our method to sets of binary black hole numerical relativity simulations with different mass ratios (q=1,2,...,8), spin configurations and separations. Our set of simulations contains non-spinning, spin-aligned and precessing simulations. We observe that the iterative procedure produces low eccentric simulations with eccentricities of the order (10^{-4}) with only one iteration. The simplicity of the procedure allows to obtain low eccentric NR simulations easily and saving computational resources. Moreover, the analytical PN formulas derived in this paper will be useful to generate eccentric hybrid waveforms.

preprint2020arXiv

Towards the routine use of subdominant harmonics in gravitational-wave inference: re-analysis of GW190412 with generation X waveform models

We re-analyse the gravitational-wave event GW190412 with state-of-the-art phenomenological waveform models. This event, which has been associated with a black hole merger, is interesting due to the significant contribution from subdominant harmonics. We use both frequency-domain and time-domain waveform models. The PhenomX waveform models constitute the fourth generation of frequency-domain phenomenological waveforms for black hole binary coalescence; they have more recently been complemented by the time-domain PhenomT models, which open up new strategies to model precession and eccentricity, and to perform tests of general relativity with the phenomenological waveforms approach. Both PhenomX and PhenomT have been constructed with similar techniques and accuracy goals, and due to their computational efficiency this &#34;generation X&#34; model family allows the routine use of subdominant spherical harmonics in Bayesian inference. We show the good agreement between these and other state-of-the-art waveform models for GW190412, and discuss the improvements over the previous generation of phenomenological waveform models. We also discuss practical aspects of Bayesian inference such as run convergence, variations of sampling parameters, and computational cost.

preprint2020arXiv

Validity of common modelling approximations for precessing binary black holes with higher-order modes

The current paradigm for constructing waveforms from precessing compact binaries is to first construct a waveform in a non-inertial, co-precessing binary source frame followed by a time-dependent rotation to map back to the physical, inertial frame. A key insight in the construction of these models is that the co-precessing waveform can be effectively mapped to some equivalent aligned spin waveform. Secondly, the time-dependent rotation implicitly introduces $m$-mode mixing, necessitating an accurate description of higher-order modes in the co-precessing frame. We assess the efficacy of this modelling strategy in the strong field regime using Numerical Relativity simulations. We find that this framework allows for the highly accurate construction of $(2,\pm 2)$ modes in our data set, while for higher order modes, especially the $(2,|1|), (3,|2|)$ and $(4,|3|)$ modes, we find rather large mismatches. We also investigate a variant of the approximate map between co-precessing and aligned spin waveforms, where we only identify the slowly varying part of the time domain co-precessing waveforms with the aligned-spin one, but find no significant improvement. Our results indicate that the simple paradigm to construct precessing waveforms does not provide an accurate description of higher order modes in the strong-field regime, and demonstrate the necessity for modelling mode asymmetries and mode-mixing to significantly improve the description of precessing higher order modes.