Advancing Gravitational Wave Astronomy: Novel Methodologies for Data Analysis and Waveform Modelling of Nanohertz and Millihertz Signals

Die Erforschung von Gravitationswellen hat unsere Sicht auf das Universum revolutioniert. Mit dem bevorstehenden Start von LISA, einem Weltraum-Gravitationswellendetektor, und neuen Berichten über Hinweise auf einen Gravitationswellenhintergrund im Nanohertz-Bereich aus Pulsar Timing Array (PTA)-Experimenten, eröffnen sich neue Möglichkeiten und Herausforderungen. Diese Dissertation entwickelt innovative Datenanalysetechniken und Wellenformmodelle, um Erkenntnisse aus diesen Beobachtungen zu gewinnen. Ein Schwerpunkt liegt auf der Untersuchung von Extreme Mass Ratio Inspirals (EMRIs) durch LISA. Diese Quellen bestehen aus kleinen, kompakten Objekten, die sich um ein zentrales Schwarzes Loch bewegen. Die Wellenformen von EMRIs bieten die Möglichkeit präziser Parametermessungen, sind jedoch aufgrund ihrer langen Signaldauer und harmonischen Komplexität schwer zu berechnen. Wir präsentieren die Implementierung einsatzbereiter EMRI-Wellenformen im Frequenzbereich für Grafikprozessoren (GPUs) und zentrale Recheneinheiten (CPUs). Zudem untersuchen wir das wissenschaftliche Potenzial von EMRIs innerhalb von Akkretionsscheiben, erforschen den Einfluss von Umwelteffekten mittels bayesianischer Methoden und bewerten die Multimessenger-Aussichten dieser Systeme. Im PTA-Bereich entwickeln wir Methoden zur Optimierung der Datenkombinationen für PTA-Analysen und tragen zum European Pulsar Timing Array bei, indem wir alternative Sampling-Pipelines für die Analyse von Gravitationswellenhintergründen und individuellen Quellen implementieren. Mit transdimensionalen Sampling-Methoden suchen wir nach einzelnen supermassiven Schwarzen Löchern und bewerten deren Signifikanz. Diese Dissertation trägt zur Weiterentwicklung der Gravitationswellenastronomie bei, indem sie neue Methoden und Modelle entwickelt, die tiefere Einblicke in die kosmischen Phänomene ermöglichen, die von LISA- und PTA-Beobachtungen erfasst werden.Gravitational wave astronomy has reshaped our understanding of the cosmos. As we look towards the future launch of LISA, a space-based gravitational wave detector, and analyze recent evidence of a nanohertz gravitational wave background from Pulsar Timing Array (PTA) experiments, new opportunities and challenges emerge. This thesis delves into developing novel data analysis techniques and waveform models to extract information from these observations. Focusing on LISA, we delve into Extreme Mass Ratio Inspirals (EMRIs). These sources consist of small compact objects spiralling into massive black holes at the centres of galaxies. Their observations are expected to provide precise parameter measurements for these systems. However, EMRI waveform generation poses challenges due to the long signal duration and large harmonic content. For the first time, we provide a fast implementation of EMRI waveforms in the frequency domain, suitable for both graphics processing units (GPUs) and central processing units (CPUs). In addition, we explore the scientific potential of EMRIs embedded in accretion disks. Employing Bayesian inference, we investigate the measurability of environmental effects and explore these systems' multimessenger prospects. Transitioning to PTA, we develop methods to optimize data combinations for PTA analyses. We present our contributions to the second data release of the European Pulsar Timing Array collaboration, which consists of implementing alternative sampling pipelines for gravitational wave background and individual source analyses. Using trans-dimensional sampling methods, we search for individual supermassive black hole binaries and assess their significance. The burgeoning field of gravitational wave astronomy has the potential to transform our understanding of the Universe. The work in this thesis develops new approaches that will facilitate the delivery of the best possible scientific results from current and future gravitational wave observations

Assessing the impact of transient orbital resonances

One of the primary sources for the future space-based gravitational wave detector, the Laser Interferometer Space Antenna, are the inspirals of small compact objects into massive black holes in the centres of galaxies. The gravitational waveforms from such Extreme Mass Ratio Inspiral (EMRI) systems will provide measurements of their parameters with unprecedented precision, but only if the waveforms are accurately modeled. Here we explore the impact of transient orbital resonances which occur when the ratio of radial and polar frequencies is a rational number. We introduce a new Effective Resonance Model, which is an extension of the numerical kludge EMRI waveform approximation to include the effect of resonances, and use it to explore the impact of resonances on EMRI parameter estimation. For one-year inspirals, we find that the few cycle dephasings induced by 3:2 resonances can lead to systematic errors in parameter estimates, that are up to several times the typical measurement precision of the parameters. The bias is greatest for 3:2 resonances that occur closer to the central black hole. By regarding them as unknown model parameters, we further show that observations will be able to constrain the size of the changes in the orbital parameters across the resonance to a relative precision of 10% for a typical one-year EMRI observation with signal-to-noise ratio of 20. Such measurements can be regarded as tests of fundamental physics, by comparing the measured changes to those predicted in general relativity.Comment: 16 pages, 6 figure

Impact of the noise knowledge uncertainty for the science exploitation of cosmological and astrophysical stochastic gravitational wave background with LISA

This paper investigates the impact of a lack of knowledge of the instrumental noise on the characterisation of stochastic gravitational wave backgrounds with the Laser Interferometer Space Antenna (LISA). We focus on constraints on modelled backgrounds that represent the possible backgrounds from the mergers of binary black holes of stellar origin, from primordial black hole generation, from non-standard inflation, and from sound wave production during cosmic fluid phase transitions. We use splines to model generic, slowly varying, uncertainties in the auto and cross-spectral densities of the LISA time delay interferometry channels. We find that allowing for noise knowledge uncertainty in this way leads to one to two orders of magnitude degradation in our ability to constrain stochastic backgrounds, and a corresponding increase in the background energy density required for a confident detection. We also find that to avoid this degradation, the LISA noise would have to be known at the sub-percent level, which is unlikely to be achievable in practice

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI waveforms. Because the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, we also introduce an improved Augmented Analytic Kludge (AAK) model that describes generic Kerr inspirals. Both waveform models can be accelerated using graphics processing unit (GPU) hardware. With the GPU-accelerated waveforms in hand, a variety of studies are performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation. We find relativistic EMRI waveform templates can be generated with fewer harmonic modes ($\sim10-100$) without biasing signal extraction. However, we show for the first time that extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space.Comment: 26 pages, 12 Figures, FastEMRIWaveforms Package: bhptoolkit.org/FastEMRIWaveforms

Probing Accretion Physics with Gravitational Waves

Gravitational-wave observations of extreme mass ratio inspirals (EMRIs) offer the opportunity to probe the environments of active galactic nuclei (AGN) through the torques that accretion disks induce on the binary. Within a Bayesian framework, we study how well such environmental effects can be measured using gravitational wave observations from the Laser Interferometer Space Antenna (LISA). We focus on the torque induced by planetary-type migration on quasicircular inspirals, and use different prescriptions for geometrically thin and radiatively efficient disks. We find that LISA could detect migration for a wide range of disk viscosities and accretion rates, for both $\alpha$ and $\beta$ disk prescriptions. For a typical EMRI with masses $50M_\odot+10^6M_\odot$, we find that LISA could distinguish between migration in $\alpha$ and $\beta$ disks and measure the torque amplitude with $\sim 20\%$ relative precision. Provided an accurate torque model, we also show how to turn gravitational-wave measurements of the torque into constraints on the disk properties. Furthermore, we show that, if an electromagnetic counterpart is identified, the multimessenger observations of the AGN EMRI system will yield direct measurements of the disk viscosity. Finally, we investigate the impact of neglecting environmental effects in the analysis of the gravitational-wave signal, finding 3$\sigma$ biases in the primary mass and spin, and showing that ignoring such effects can lead to false detection of a deviation from general relativity. This work demonstrates the scientific potential of gravitational observations as probes of accretion-disk physics, accessible so far through electromagnetic observations only

Quality over Quantity: Optimizing pulsar timing array analysis for stochastic and continuous gravitational wave signals

The search for gravitational waves using Pulsar Timing Arrays (PTAs) is a computationally expensive complex analysis that involves source-specific noise studies. As more pulsars are added to the arrays, this stage of PTA analysis will become increasingly challenging. Therefore, optimizing the number of included pulsars is crucial to reduce the computational burden of data analysis. Here, we present a suite of methods to rank pulsars for use within the scope of PTA analysis. First, we use the maximization of the signal-to-noise ratio as a proxy to select pulsars. With this method, we target the detection of stochastic and continuous gravitational wave signals. Next, we present a ranking that minimizes the coupling between spatial correlation signatures, namely monopolar, dipolar, and Hellings & Downs correlations. Finally, we also explore how to combine these two methods. We test these approaches against mock data using frequentist and Bayesian hypothesis testing. For equal-noise pulsars, we find that an optimal selection leads to an increase in the log-Bayes factor two times steeper than a random selection for the hypothesis test of a gravitational wave background versus a common uncorrelated red noise process. For the same test but for a realistic EPTA dataset, a subset of 25 pulsars selected out of 40 can provide a log-likelihood ratio that is $89\%$ of the total, implying that an optimally selected subset of pulsars can yield results comparable to those obtained from the whole array. We expect these selection methods to play a crucial role in future PTA data combinations

Fast and Fourier: Extreme Mass Ratio Inspiral Waveforms in the Frequency Domain

Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the first ready-to-use Schwarzschild eccentric EMRI waveform implementation in the frequency domain for use with either graphics processing units (GPUs) or central processing units (CPUs). We present the overall waveform implementation and test the accuracy and performance of the frequency domain waveforms against the time domain implementation. On GPUs, the frequency domain waveform takes in median $0.044$ seconds to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses $\geq 2 \times 10^6 \,{\rm M_\odot}$ and initial eccentricities $e_0 > 0.2$. On CPUs, the median waveform evaluation time is $5$ seconds, and it is five times faster in the frequency domain than in the time domain. Using a sparser frequency array can further speed up the waveform generation, reaching up to $0.3$ seconds. This enables us to perform, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs. Future EMRI models which encompass wider source characteristics (particularly black hole spin and generic orbit geometries) will require significantly more harmonics. Frequency-domain models will be essential analysis tools for these astrophysically realistic and important signals.Comment: 23 pages, 6 figure

Workshop on Gravitational-Wave Astrophysics for Early Career Scientists

Gravitational-wave science is rapidly growing in maturity as a research area; in May 2021 the next generation of gravitational-wave scientists gathered together to create a vision of the future of the field.Non peer reviewe

Waveform Modelling for the Laser Interferometer Space Antenna

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.Comment: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from communit

Cosmology with the Laser Interferometer Space Antenna

254 pags:, 44 figs.The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational-wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational-wave observations by LISA to probe the universe.This work is partly supported by: A.G. Leventis Foundation; Academy of Finland Grants 328958 and 345070; Alexander S. Onassis Foundation, Scholarship ID: FZO 059-1/2018-2019; Amaldi Research Center funded by the MIUR program “Dipartimento di Eccellenza” (CUP: B81I18001170001); ASI Grants No. 2016-24-H.0 and No. 2016-24-H.1-2018; Atracción de Talento Grant 2019-T1/TIC-15784; Atracción de Talento contract no. 2019-T1/TIC-13177 granted by the Comunidad de Madrid; Ayuda ‘Beatriz Galindo Senior’ by the Spanish ‘Ministerio de Universidades’, Grant BG20/00228; Basque Government Grant (IT-979-16); Belgian Francqui Foundation; Centre national d’Etudes spatiales; Ben Gurion University Kreitman Fellowship, and the Israel Academy of Sciences and Humanities (IASH) & Council for Higher Education (CHE) Excellence Fellowship Program for International Postdoctoral Researchers; Centro de Excelencia Severo Ochoa Program SEV-2016-0597; CERCA program of the Generalitat de Catalunya; Cluster of Excellence “Precision Physics, Fundamental Interactions, and Structure of Matter” (PRISMA? EXC 2118/1); Comunidad de Madrid, Contrato de Atracción de Talento 2017-T1/TIC-5520; Czech Science Foundation GAČR, Grant No. 21-16583M; Delta ITP consortium; Department of Energy under Grant No. DE-SC0008541, DE-SC0009919 and DESC0019195; Deutsche Forschungsgemeinschaft (DFG), Project ID 438947057; Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy - EXC 2121 Quantum Universe - 390833306; European Structural and Investment Funds and the Czech Ministry of Education, Youth and Sports (Project CoGraDS - CZ.02.1.01/0.0/0.0/15 003/0000437); European Union’s H2020 ERC Consolidator Grant “GRavity from Astrophysical to Microscopic Scales” (Grant No. GRAMS-815673); European Union’s H2020 ERC, Starting Grant Agreement No. DarkGRA-757480; European Union’s Horizon 2020 programme under the Marie Sklodowska-Curie Grant Agreement 860881 (ITN HIDDeN); European Union’s Horizon 2020 Research and Innovation Programme Grant No. 796961, “AxiBAU” (K.S.); European Union’s Horizon 2020 Research Council grant 724659 MassiveCosmo ERC-2016-COG; FCT through national funds (PTDC/FIS-PAR/31938/2017) and through project “BEYLA – BEYond LAmbda” with Ref. Number PTDC/FIS-AST/0054/2021; FEDER-Fundo Europeu de Desenvolvimento Regional through COMPETE2020 - Programa Operacional Competitividade e Internacionalização (POCI-01-0145- FEDER-031938) and research Grants UIDB/04434/2020 and UIDP/04434/2020; Fondation CFM pour la Recherche in France; Foundation for Education and European Culture in Greece; French ANR project MMUniverse (ANR-19-CE31-0020); FRIA Grant No.1.E.070.19F of the Belgian Fund for Research, F.R. S.-FNRS Fundação para a Ciência e a Tecnologia (FCT) through Contract No. DL 57/2016/CP1364/ CT0001; Fundação para a Ciência e a Tecnologia (FCT) through Grants UIDB/04434/2020, UIDP/04434/ 2020, PTDC/FIS-OUT/29048/2017, CERN/FIS-PAR/0037/2019 and “CosmoTests – Cosmological tests of gravity theories beyond General Relativity” CEECIND/00017/2018; Generalitat Valenciana Grant PROMETEO/2021/083; Grant No. 758792, project GEODESI; Government of Canada through the Department of Innovation, Science and Economic Development and Province of Ontario through the Ministry of Colleges and Universities; Grants-in-Aid for JSPS Overseas Research Fellow (No. 201960698); I?D Grant PID2020-118159GB-C41 of the Spanish Ministry of Science and Innovation; INFN iniziativa specifica TEONGRAV; Israel Science Foundation (Grant No. 2562/20); Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Nos. 20H01899 and 20H05853; IFT Centro de Excelencia Severo Ochoa Grant SEV-2; Kavli Foundation and its founder Fred Kavli; Minerva Foundation; Ministerio de Ciencia e Innovacion Grant PID2020-113644GB-I00; NASA Grant 80NSSC19K0318; NASA Hubble Fellowship grants No. HST-HF2-51452.001-A awarded by the Space Telescope Science Institute with NASA contract NAS5-26555; Netherlands Organisation for Science and Research (NWO) Grant Number 680-91-119; new faculty seed start-up grant of the Indian Institute of Science, Bangalore, the Core Research Grant CRG/2018/002200 of the Science and Engineering; NSF Grants PHY-1820675, PHY-2006645 and PHY-2011997; Polish National Science Center Grant 2018/31/D/ ST2/02048; Polish National Agency for Academic Exchange within the Polish Returns Programme under Agreement PPN/PPO/2020/1/00013/U/00001; Pró-Reitoria de Pesquisa of Universidade Federal de Minas Gerais (UFMG) under Grant No. 28359; Ramón y Cajal Fellowship contract RYC-2017-23493; Research Project PGC2018-094773-B-C32 [MINECO-FEDER]; Research Project PGC2018-094773-B-C32 [MINECO-FEDER]; ROMFORSK Grant Project. No. 302640; Royal Society Grant URF/R1/180009 and ERC StG 949572: SHADE; Shota Rustaveli National Science Foundation (SRNSF) of Georgia (Grant FR/18-1462); Simons Foundation/SFARI 560536; SNSF Ambizione grant; SNSF professorship Grant (No. 170547); Spanish MINECO’s “Centro de Excelencia Severo Ochoa” Programme Grants SEV-2016- 0597 and PID2019-110058GB-C22; Spanish Ministry MCIU/AEI/FEDER Grant (PGC2018-094626-BC21); Spanish Ministry of Science and Innovation (PID2020-115845GB-I00/AEI/10.13039/ 501100011033); Spanish Proyectos de I?D via Grant PGC2018-096646-A-I00; STFC Consolidated Grant ST/T000732/1; STFC Consolidated Grants ST/P000762/1 and ST/T000791/1; STFC Grant ST/ S000550/1; STFC Grant ST/T000813/1; STFC Grants ST/P000762/1 and ST/T000791/1; STFC under the research Grant ST/P000258/1; Swiss National Science Foundation (SNSF), project The Non-Gaussian Universe and Cosmological Symmetries, Project Number: 200020-178787; Swiss National Science Foundation Professorship Grants No. 170547 and No. 191957; SwissMap National Center for Competence in Research; “The Dark Universe: A Synergic Multi-messenger Approach” Number 2017X7X85K under the MIUR program PRIN 2017; UK Space Agency; UKSA Flagship Project, Euclid.Peer reviewe