Researcher profile

Katerina Chatziioannou

Katerina Chatziioannou contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

17 published item(s)

preprint2022arXiv

Accurate modeling and mitigation of overlapping signals and glitches in gravitational-wave data

The increasing sensitivity of gravitational-wave detectors has brought about an increase in the rate of astrophysical signal detections as well as the rate of "glitches"; transient and non-Gaussian detector noise. Temporal overlap of signals and glitches in the detector presents a challenge for inference analyses that typically assume the presence of only Gaussian detector noise. In this study we perform an extensive exploration of the efficacy of a recently proposed method that models the glitch with sine-Gaussian wavelets while simultaneously modeling the signal with compact-binary waveform templates. We explore a wide range of glitch families and signal morphologies and demonstrate that the joint modeling of glitches and signals (with wavelets and templates respectively) can reliably separate the two. We find that the glitches that most affect parameter estimation are also the glitches that are well modeled by such wavelets due to their compact time-frequency signature. As a further test, we investigate the robustness of this analysis against waveform systematics like those arising from the exclusion of higher-order modes and spin-precession effects. Our analysis provides an estimate of the signal parameters; the glitch waveform to be subtracted from the data; and an assessment of whether some detected excess power consists of a glitch, signal, or both. We analyze the low-significance triggers (191225_215715 and 200114_020818) and find that they are both consistent with glitches overlapping high-mass signals.

preprint2022arXiv

Comparing Bayes factors and hierarchical inference for testing general relativity with gravitational waves

In the context of testing general relativity with gravitational waves, constraints obtained with multiple events are typically combined either through a hierarchical formalism or though a combined multiplicative Bayes factor. We show that the well-known dependence of Bayes factors on the analysis priors in regions of the parameter space without likelihood support can lead to strong confidence in favor of incorrect conclusions when one employs the multiplicative Bayes factor. Bayes factors $\mathcal{O}(1)$ are ambivalent as they depend sensitively on the analysis priors, which are rarely set in a principled way; additionally, combined Bayes factors $>\mathcal{O}(10^3)$ can be obtained in favor of the incorrect conclusion depending on the analysis priors when many $\mathcal{O}(1)$ Bayes factors are multiplied, and specifically when the priors are much wider than the underlying population. The hierarchical analysis that instead infers the ensemble distribution of the individual beyond-general-relativity constraints does not suffer from this problem, and generically converges to favor the correct conclusion. Rather than a naive multiplication, a more reliable Bayes factor can be computed from the hierarchical analysis. We present a number of toy models showing that the practice of multiplying Bayes Factors can lead to incorrect conclusions.

preprint2022arXiv

Equation-of-state insensitive relations after GW170817

The thermodynamic relation between pressure and density (i.e. the equation of state) of cold supranuclear matter is critical in describing neutron stars, yet it remains one of the largest uncertainties in nuclear physics. The extraction of tidal deformabilities from the gravitational waves emitted in the coalescence of neutron star binaries, such as GW170817, is a promising tool to probe this thermodynamic relation. Equation-of-state insensitive relations between symmetric and antisymmetric combinations of individual tidal deformabilities, the so-called "binary Love relations", have proven important to infer the radius of neutron stars, and thus constrain the equation of state, from such gravitational waves. A similar set of relations between the moment of inertia, the tidal deformability, the quadrupole moment, and the compactness of neutron stars, the so-called `I-Love-Q' and `C-Love' relations, allow for future tests of General Relativity in the extreme gravity regime. But even the most insensitive of such relations still presents some degree of equation-of-state variability that could introduce systematic uncertainties in parameter extraction and in model selection. We here reduce this variability by more than 50% by imposing a prior on the allowed set of equations of state, derived from the posteriors generated from the analysis of GW170817. The resulting increase in insensitivity reduces systematic uncertainties in the extraction of the tidal deformability from future gravitational wave observations, although statistical uncertainties currently dominate the error budget, and will continue to do so until the era of Voyager-class detectors.

preprint2022arXiv

Implicit correlations within phenomenological parametric models of the neutron star equation of state

The rapid increase in the number and precision of astrophysical probes of neutron stars in recent years allows for the inference of their equation of state. Observations target different macroscopic properties of neutron stars which vary from star to star, such as mass and radius, but the equation of state allows for a common description of all neutron stars. To connect these observations and infer the properties of dense matter and neutron stars simultaneously, models for the equation of state are introduced. Parametric models rely on carefully engineered functional forms that reproduce a large array of realistic equations of state. Such models benefit from their simplicity but are limited because any finite-parameter model cannot accurately approximate all possible equations of state. Nonparametric models overcome this by increasing model freedom at the cost of increased complexity. In this study, we compare common parametric and nonparametric models, quantify the limitations of the former, and study the impact of modeling on our current understanding of high-density physics. We show that parametric models impose strongly model-dependent, and sometimes opaque, correlations between density scales. Such interdensity correlations result in tighter constraints that are unsupported by data and can lead to biased inference of the equation of state and of individual neutron star properties.

preprint2022arXiv

Probing neutron stars with the full pre-merger and post-merger gravitational wave signal from binary coalescences

The gravitational wave signal emitted during the coalescence of two neutron stars carries information about the stars' internal structure. During the long inspiral phase the main matter observable is the tidal interaction between the binary components, an effect that can be parametrically modeled with compact-binary solutions to General Relativity. After the binary merger the main observable is frequency modes of the remnant, most commonly giving rise to a short-duration signal accessible only through numerical simulations. The complicated morphology and the decreasing detector sensitivity in the relevant frequencies currently hinder detection of the post-merger signal and motivate separate analyses for the pre-merger and post-merger data. However, planned and ongoing detector improvements could soon put the post-merger signal within reach. In this study we target the whole pre-merger and post-merger signal without an artificial separation at the binary merger. We construct a hybrid analysis that models the inspiral with templates based on analytical calculations and calibrated to numerical relativity and the post-merger signal with a flexible morphology-independent analysis. Applying this analysis to GW170817 we find, as expected, that the post-merger signal remains undetected. We further study simulated signals and find that we can reconstruct the full signal and simultaneously estimate both the pre-merger tidal deformation and the post-merger signal frequency content. Our analysis allows us to study neutron star physics using all the data available and directly test the pre-merger and post-merger signal for consistency thus probing effects such as the onset of the hadron-quark phase transition.

preprint2022arXiv

Uncertainty limits on neutron star radius measurements with gravitational waves

Upcoming observing campaigns with improved detectors will yield numerous detections of gravitational waves from neutron star binary inspirals. Rare loud signals together with numerous signals of moderate strength promise stringent constraints on the properties of neutron star matter, with a projected radius statistical uncertainty of $50-200$m with ${\cal{O}}(2000)$ sources. Given this precision we revisit all analysis assumptions and identify sources of systematic errors, quantify their impact on radius extraction, and discuss their relative importance and ways to mitigate them.

preprint2021arXiv

Modeling compact binary signals and instrumental glitches in gravitational wave data

Transient non-gaussian noise in gravitational wave detectors, commonly referred to as glitches, pose challenges for inference of the astrophysical properties of detected signals when the two are coincident in time. Current analyses aim towards modeling and subtracting the glitches from the data using a flexible, morphology-independent model in terms of sine-gaussian wavelets before the signal source properties are inferred using templates for the compact binary signal. We present a new analysis of gravitational wave data that contain both a signal and glitches by simultaneously modeling the compact binary signal in terms of templates and the instrumental glitches using sine-gaussian wavelets. The model for the glitches is generic and can thus be applied to a wide range of glitch morphologies without any special tuning. The simultaneous modeling of the astrophysical signal with templates allows us to efficiently separate the signal from the glitches, as we demonstrate using simulated signals injected around real O2 glitches in the two LIGO detectors. We show that our new proposed analysis can separate overlapping glitches and signals, estimate the compact binary parameters, and provide ready-to-use glitch-subtracted data for downstream inference analyses.

preprint2021arXiv

The BayesWave analysis pipeline in the era of gravitational wave observations

We describe updates and improvements to the BayesWave gravitational wave transient analysis pipeline, and provide examples of how the algorithm is used to analyze data from ground-based gravitational wave detectors. BayesWave models gravitational wave signals in a morphology-independent manner through a sum of frame functions, such as Morlet-Gabor wavelets or chirplets. BayesWave models the instrument noise using a combination of a parametrized Gaussian noise component and non-stationary and non-Gaussian noise transients. Both the signal model and noise model employ trans-dimensional sampling, with the complexity of the model adapting to the requirements of the data. The flexibility of the algorithm makes it suitable for a variety of analyses, including reconstructing generic unmodeled signals; cross checks against modeled analyses for compact binaries; as well as separating coherent signals from incoherent instrumental noise transients (glitches). The BayesWave model has been extended to account for gravitational wave signals with generic polarization content and the simultaneous presence of signals and glitches in the data. We describe updates in the BayesWave prior distributions, sampling proposals, and burn-in stage that provide significantly improved sampling efficiency. We present standard review checks indicating the robustness and convergence of the BayesWave trans-dimensional sampler.

preprint2020arXiv

A Population-Informed Mass Estimate for Pulsar J0740+6620

Galactic double neutron star systems have a tight mass distribution around $\sim 1.35 M_{\odot}$, but the mass distribution of all known pulsars is broader. Here we reconstruct the Alsing, et al. (2018) bimodal mass distribution of pulsars observed in binary systems, incorporating data from observations of J0740+6620 which were not available at the time of that work. Because J0740+6620 is an outlier in the mass distribution with non-negligible uncertainty in its mass measurement, its mass receives a large correction from the population, becoming $m_{J0740+6620} = 2.03^{+0.10}_{-0.08} \, M_\odot$ (median and 68\% CI). Stochastic samples from our population model, including population-informed pulsar mass estimates, are available at https://github.com/farr/AlsingNSMassReplication .

preprint2020arXiv

Detection and parameter estimation of binary neutron star merger remnants

Detection and parameter estimation of binary neutron star merger remnants can shed light on the physics of hot matter at supranuclear densities. Here we develop a fast, simple model that can generate gravitational waveforms, and show it can be used for both detection and parameter estimation of post-merger remnants. The model consists of three exponentially-damped sinusoids with a linear frequency-drift term. The median fitting factors between the model waveforms and numerical-relativity simulations exceed 0.90. We detect remnants at a post-merger signal-to-noise ratio of $\ge 7$ using a Bayes-factor detection statistic with a threshold of 3000. We can constrain the primary post-merger frequency to $\pm_{1.2}^{1.4}\%$ at post-merger signal-to-noise ratios of 15 with an increase in precision to $\pm_{0.2}^{0.3}\%$ for post-merger signal-to-noise ratios of 50. The tidal coupling constant can be constrained to $\pm^{9}_{12}\%$ at post-merger signal-to-noise ratios of 15, and $\pm 5\%$ at post-merger signal-to-noise ratios of 50 using a hierarchical inference model.

preprint2020arXiv

Distinguishing binary neutron star from neutron star-black hole mergers with gravitational waves

The gravitational-wave signal from the merger of two neutron stars cannot be easily differentiated from the signal produced by a comparable-mass mixed binary of a neutron star and a black hole. Indeed, both binary types can account for the gravitational-wave signal GW170817 even if its electromagnetic counterpart emission is taken into account. We propose a method that {requires neither information from the post-inspiral phase of the binary nor an electromagnetic counterpart} to identify mixed binaries of neutron stars merging with low-mass black holes using gravitational-waves alone. This method is based on the fact that certain neutron star properties that can be measured with gravitational-waves are common or similar for all neutron stars. For example all neutron stars share the same equation of state and if the latter is hadronic, neutron stars have similar radii. If a mixed binary is misidentified as a neutron star binary, the inferred neutron star properties will be misestimated and appear as outliers in a population of low-mass binaries. We show that as few as $\sim 5$ low-mass events will allow for the identification of the type of one event at the $80\%$ confidence level. We model the population of low-mass binaries with a hierarchical mixture model and show that we can constrain the existence of mixed binaries or measure their abundance relative to neutron star binaries to $\sim 0.1$ at the $68\%$ credible level with 100 events.

preprint2020arXiv

Including higher order multipoles in gravitational-wave models for precessing binary black holes

Estimates of the source parameters of gravitational-wave (GW) events produced by compact binary mergers rely on theoretical models for the GW signal. We present the first frequency-domain model for inspiral, merger and ringdown of the GW signal from precessing binary-black-hole systems that also includes multipoles beyond the leading-order quadrupole. Our model, {\tt PhenomPv3HM}, is a combination of the higher-multipole non-precessing model {\tt PhenomHM} and the spin-precessing model {\tt PhenomPv3} that includes two-spin precession via a dynamical rotation of the GW multipoles. We validate the new model by comparing to a large set of precessing numerical-relativity simulations and find excellent agreement across the majority of the parameter space they cover. For mass ratios $<5$ the mismatch improves, on average, from $\sim6\%$ to $\sim 2\%$ compared to {\tt PhenomPv3} when we include higher multipoles in the model. However, we find mismatches $\sim8\%$ for the mass-ratio $6$ and highly spinning simulation. As a first application of the new model we have analysed the binary black hole event GW170729. We find larger values for the primary black hole mass of $58.25^{+11.73}_{-12.53} \, M_\odot$ (90\% credible interval). The lower limit ($\sim 46 \, M_\odot$) is comparable to the proposed maximum black hole mass predicted by different stellar evolution models due to the pulsation pair-instability supernova (PPISN) mechanism. If we assume that the primary \ac{BH} in GW170729 formed through a PPISN then out of the four PPISN models we considered only the model of Woosley (2017) is consistent with our mass measurements at the 90\% level.

preprint2020arXiv

Inference of the neutron star equation of state from cosmological distances

Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of signal will be redshifted towards the detector most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the post-merger frequency, are better measured from the distant weak sources from $z\sim 1$.

preprint2020arXiv

Reconstructing gravitational wave signals from binary black hole mergers with minimal assumptions

We present a systematic comparison of the binary black hole (BBH) signal waveform reconstructed by two independent and complementary approaches used in LIGO and Virgo source inference: a template-based analysis, and a morphology-independent analysis. We apply the two approaches to real events and to two sets of simulated observations made by adding simulated BBH signals to LIGO and Virgo detector noise. The first set is representative of the 10 BBH events in the first Gravitational Wave Transient Catalog (GWTC-1). The second set is constructed from a population of BBH systems with total mass and signal strength in the ranges that ground based detectors are typically sensitive. We find that the reconstruction quality of the GWTC-1 events is consistent with the results of both sets of simulated signals. We also demonstrate a simulated case where the presence of a mismodelled effect in the observed signal, namely higher order modes, can be identified through the morphology-independent analysis. This study is relevant for currently progressing and future observational runs by LIGO and Virgo.

preprint2020arXiv

Studying strong phase transitions in neutron stars with gravitational waves

The composition of neutron stars at the extreme densities reached in their cores is currently unknown. Besides nuclear matter of normal neutrons and protons, the cores of neutron stars might harbor exotic matter such as deconfined quarks. In this paper we study strong hadron-quark phase transitions in the context of gravitational wave observations of inspiraling neutron stars. We consider upcoming detections of neutron star coalescences and model the neutron star equations of state with phase transitions through the Constant-Speed-of-Sound parametrization. We use the fact that neutron star binaries with one or more hadron-quark hybrid stars can exhibit qualitatively different tidal properties than binaries with hadronic stars of the same mass, and hierarchically model the masses and tidal properties of simulated populations of binary neutron star inspiral signals. We explore the parameter space of phase transitions and discuss under which conditions future observations of binary neutron star inspirals can identify this effect and constrain its properties, in particular the threshold density at which the transition happens and the strength of the transition. We find that if the detected population of binary neutron stars contains both hadronic and hybrid stars, the onset mass and strength of a sufficiently strong phase transition can be constrained with 50-100 detections. If the detected neutron stars are exclusively hadronic or hybrid, then it is possible to place lower or upper limits on the transition density and strength.

preprint2020arXiv

You Can&#39;t Always Get What You Want: The Impact of Prior Assumptions on Interpreting GW190412

GW190412 is the first observation of a black hole binary with definitively unequal masses. GW190412&#39;s mass asymmetry, along with the measured positive effective inspiral spin, allowed for inference of a component black hole spin: the primary black hole in the system was found to have a dimensionless spin magnitude between 0.17 and 0.59 (90% credible range). We investigate how the choice of priors for the spin magnitudes and tilts of the component black holes affect the robustness of parameter estimates for GW190412, and report Bayes factors across a suite of prior assumptions. Depending on the waveform family used to describe the signal, we find either marginal to moderate (2:1-6:1) or strong ($\gtrsim$ 20:1) support for the primary black hole being spinning compared to cases where only the secondary is allowed to have spin. We show how these choices influence parameter estimates, and find the asymmetric masses and positive effective inspiral spin of GW190412 to be qualitatively, but not quantitatively, robust to prior assumptions. Our results highlight the importance of both considering astrophysically motivated or population-based priors in interpreting observations and considering their relative support from the data.

preprint2019arXiv

A morphology-independent search for gravitational wave echoes in data from the first and second observing runs of Advanced LIGO and Advanced Virgo

Gravitational wave echoes have been proposed as a smoking-gun signature of exotic compact objects with near-horizon structure. Recently there have been observational claims that echoes are indeed present in stretches of data from Advanced LIGO and Advanced Virgo immediately following gravitational wave signals from presumed binary black hole mergers, as well as a binary neutron star merger. In this paper we deploy a morphology-independent search algorithm for echoes introduced in Tsang et al., Phys. Rev. D 98, 024023 (2018), which (a) is able to accurately reconstruct a possible echoes signal with minimal assumptions about their morphology, and (b) computes Bayesian evidences for the hypotheses that the data contain a signal, an instrumental glitch, or just stationary, Gaussian noise. Here we apply this analysis method to all the significant events in the first Gravitational Wave Transient Catalog (GWTC-1), which comprises the signals from binary black hole and binary neutron star coalescences found during the first and second observing runs of Advanced LIGO and Advanced Virgo. In all cases, the ratios of evidences for signal versus noise and signal versus glitch do not rise above their respective &#34;background distributions&#34; obtained from detector noise, the smallest $p$-value being 3% (for event GW170823). Hence we find no statistically significant evidence for echoes in GWTC-1.