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.

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called ``Erebor,'' designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog $\sim12000$ Galactic Binaries ($\sim8000$ as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

We construct a one-parameter family of lattice models starting from a two-dimensional rational conformal field theory on a torus with a regular lattice of holes, each of which is equipped with a conformal boundary condition. The lattice model is obtained by cutting the surface into triangles with clipped-off edges using open channel factorisation. The parameter is given by the hole radius. At finite radius, high energy states are suppressed and the model is effectively finite. In the zero-radius limit, it recovers the CFT amplitude exactly. In the touching hole limit, one obtains a topological field theory. If one chooses a special conformal boundary condition which we call "cloaking boundary condition", then for each value of the radius the fusion category of topological line defects of the CFT is contained in the lattice model. The fact that the full topological symmetry of the initial CFT is realised exactly is a key feature of our lattice models. We provide an explicit recursive procedure to evaluate the interaction vertex on arbitrary states. As an example, we study the lattice model obtained from the Ising CFT on a torus with one hole, decomposed into two lattice cells. We numerically compare the truncated lattice model to the CFT expression obtained from expanding the boundary state in terms of the hole radius and we find good agreement at intermediate values of the radius.

We investigate Refined Algebraic Quantization (RAQ) with group averaging in a constrained Hamiltonian system with unreduced phase space T^*R^4 and gauge group SL(2,R). The reduced phase space M is connected and contains four mutually disconnected `regular' sectors with topology R x S^1, but these sectors are connected to each other through an exceptional set where M is not a manifold and where M has non-Hausdorff topology. The RAQ physical Hilbert space H_{phys} decomposes as H_{phys} = (direct sum of) H_i, where the four subspaces H_i naturally correspond to the four regular sectors of M. The RAQ observable algebra A_{obs}, represented on H_{phys}, contains natural subalgebras represented on each H_i. The group averaging takes place in the oscillator representation of SL(2,R) on L^2(R^{2,2}), and ensuring convergence requires a subtle choice for the test state space: the classical analogue of this choice is to excise from M the exceptional set while nevertheless retaining information about the connections between the regular sectors. A quantum theory with the Hilbert space H_{phys} and a finitely-generated observable subalgebra of A_{obs} is recovered through both Ashtekar's Algebraic Quantization and Isham's group theoretic quantization.

The future space based gravitational wave detector LISA (Laser Interferometer Space Antenna) will observe millions of Galactic binaries constantly present in the data stream. A small fraction of this population (of the order of several thousand) will be individually resolved. One of the challenging tasks from the data analysis point of view will be to estimate the parameters of resolvable galactic binaries while disentangling them from each other and from other gravitational wave sources present in the data. This problem is quite often referred to as a global fit in the field of LISA data analysis. A Bayesian framework is often used to infer the parameters of the sources and their number. The efficiency of the sampling techniques strongly depends on the proposals, especially in the multi-dimensional parameter space. In this paper we demonstrate how we can use neural density estimators, and in particular Normalising flows, in order to build proposals which significantly improve the convergence of sampling. We also demonstrate how these methods could help in building priors based on physical models and provide an alternative way to represent the catalogue of identified gravitational wave sources.

The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts which have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop "Challenges and opportunities of high-frequency gravitational wave detection" held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative.

A detailed description of how black holes grow in full, non-linear general relativity is presented. The starting point is the notion of dynamical horizons. Expressions of fluxes of energy and angular momentum carried by gravitational waves across these horizons are obtained. Fluxes are local and the energy flux is positive. Change in the horizon area is related to these fluxes. A notion of angular momentum and energy is associated with cross-sections of the horizon and balance equations, analogous to those obtained by Bondi and Sachs at null infinity, are derived. These in turn lead to generalizations of the first and second laws of black hole mechanics. The relation between dynamical horizons and their asymptotic states --the isolated horizons-- is discussed briefly. The framework has potential applications to numerical, mathematical, astrophysical and quantum general relativity.

We review aspects of loop quantum gravity in a pedagogical manner, with the aim of enabling a precise but critical assessment of its achievements so far. We emphasise that the off-shell (`strong') closure of the constraint algebra is a crucial test of quantum space-time covariance, and thereby of the consistency, of the theory. Special attention is paid to the appearance of a large number of ambiguities, in particular in the formulation of the Hamiltonian constraint. Developing suitable approximation methods to establish a connection with classical gravity on the one hand, and with the physics of elementary particles on the other, remains a major challenge.

One of the most exciting prospects for the LISA gravitational wave observatory is the detection of gravitational radiation from the inspiral of a compact object into a supermassive black hole. The large inspiral parameter space and low amplitude of the signal makes detection of these sources computationally challenging. We outline here a first cut data analysis scheme that assumes realistic computational resources. In the context of this scheme, we estimate the signal-to-noise ratio that a source requires to pass our thresholds and be detected. Combining this with an estimate of the population of sources in the Universe, we estimate the number of inspiral events that LISA could detect. The preliminary results are very encouraging -- with the baseline design, LISA can see inspirals out to a redshift z=1 and should detect over a thousand events during the mission lifetime.

An astrophysical population of supermassive black hole binaries is thought to be the strongest source of gravitational waves in the frequency range covered by Pulsar Timing Arrays (PTAs). A potential cause for concern is that the standard cross-correlation method used in PTA data analysis assumes that the signals are isotropically distributed and Gaussian random, while the signals from a black hole population are likely to be anisotropic and deterministic. Here we argue that while the conventional analysis is not optimal, it is not hopeless either, as the standard Hellings-Downs correlation curve turns out to hold for point sources, and the small effective number of signal samples blurs the distinction between Gaussian and deterministic signals. Possible improvements to the standard cross-correlation analysis that account for the anisotropy of the signal are discussed.

We present the first search for gravitational waves from sub-solar mass compact-binary mergers which allows for non-negligible orbital eccentricity. Sub-solar mass black holes are a signature of primordial origin black holes, which may be a component of dark matter. To produce binary coalescences, primordial black holes may form close binaries either in the early universe or more recently through dynamical interactions. A signature of dynamical formation would be the observation of non-circularized orbits. We search for black hole mergers where the primary mass is $0.1-7 M_{\odot}$ and the secondary mass is $0.1-1 M_{\odot}$. We allow for eccentricity up to $\sim0.3$ at a dominant-mode gravitational-wave frequency of 10 Hz for binaries with at least one component with mass $>0.5 M_{\odot}$. We find no convincing candidates in the public LIGO data. The two most promising candidates have a false alarm rate of 1 per 3 and 4 years, respectively, which combined is only a $\sim 2.4\sigma$ deviation from the expected Poisson rate. Given the marginal statistical significance, we place upper limits on the rate of sub-solar mass mergers under the assumption of a null observation and compare how these limits may inform the possible dark matter contribution.

We present an exhaustive analysis of the numerical evolution of the Einstein-Klein-Gordon equations for the case of a real scalar field endowed with a quadratic self-interaction potential. The self-gravitating equilibrium configurations are called oscillatons and are close relatives of boson stars, their complex counterparts. Unlike boson stars, for which the oscillations of the two components of the complex scalar field are such that the spacetime geometry remains static, oscillatons give rise to a geometry that is time-dependent and oscillatory in nature. However, they can still be classified into stable (S-branch) and unstable (U-branch) cases. We have found that S-oscillatons are indeed stable configurations under small perturbations and typically migrate to other S-profiles when perturbed strongly. On the other hand, U-oscillatons are intrinsically unstable: they migrate to the S-branch if their mass is decreased and collapse to black holes if their mass is increased even by a small amount. The S-oscillatons can also be made to collapse to black holes if enough mass is added to them, but such collapse can be efficiently prevented by the gravitational cooling mechanism in the case of diluted oscillatons.

The self field of a charged particle has a component that diverges at the particle. We use both coordinate and covariant approaches to compute an expansion of this singular field for generic geodesic orbits in Schwarzschild spacetime for scalar, electromagnetic and graviational cases. We check agreement of both approaches and give, as an application, the calculation of previously unknown regularisation parameters. In this so-called "mode-sum regularization" approach, each mode of the field is finite, while their sum diverges. The sum may be rendered finite and convergent by the subtraction of "regularization parameters". Higher order parameters lead to faster convergence in the mode-sum. As a second example application, we compute high order expressions for the effective source approach to self-force calculations.

The dynamical evolution of self-gravitating scalar field configurations in numerical relativity is studied. The previous analysis on ground state boson stars of non-interacting fields is extended to excited states and to fields with self couplings. Self couplings can significantly change the physical dimensions of boson stars, making them much more astrophysically interesting (e.g., having mass of order 0.1 solar mass). The stable ($S$) and unstable ($U$) branches of equilibrium configurations of boson stars of self-interacting fields are studied; their behavior under perturbations and their quasi-normal oscillation frequencies are determined and compared to the non-interacting case. Excited states of boson stars with and without self-couplings are studied and compared. Excited states also have equilibrium configurations with $S$ and $U$ branch structures; both branches are intrinsically unstable under a generic perturbation but have very different instability time scales. We carried out a detailed study of the instability time scales of these configurations. It is found that highly excited states spontaneously decay through a cascade of intermediate states similar to atomic transitions.

The one-dimensional Heisenberg XXX spin chain appears in a special limit of the AdS/CFT integrable system. We review various ways of proving its integrability, and discuss the associated methods of solution. In particular, we outline the coordinate and the algebraic Bethe ansatz, giving reference to literature suitable for learning these techniques. Finally we speculate which of the methods might lift to the exact solution of the AdS/CFT system, and sketch a promising method for constructing the Baxter Q-operator of the XXX chain. It allows to find the spectrum of the model using certain algebraic techniques, while entirely avoiding Bethe's ansatz.

We examine the generic theory of a partially massless (PM) spin-two field interacting with gravity in four dimensions from a bottom-up perspective. By analyzing the most general form of the Lagrangian, we first show that if such a theory exists, its de Sitter background must admit either so(1, 5) or so(2, 4) global symmetry depending on the relative sign of the kinetic terms: the former for a positive sign the latter for a negative sign. Further analysis reveals that the coupling constant of the PM cubic self-interaction must be fixed with a purely imaginary number in the case of a positive sign. We conclude that there cannot exist a unitary theory of a PM spin-two field coupled to Einstein gravity with a perturbatively local Lagrangian. In the case of a negative sign we recover conformal gravity. As a special case of our analysis, it is shown that the PM limit of massive gravity also lacks the PM gauge symmetry.

We construct explicit consistent Kaluza-Klein reductions of type IIB supergravity on HK_4 x S^1, where HK_4 is an arbitrary four-dimensional hyper-Kaehler manifold, and on SE5, an arbitrary five-dimensional Sasaki-Einstein manifold. In the former case we obtain the bosonic action of D=5 N=4 (ungauged) supergravity coupled to two vector multiplets. For the SE_5 case we extend a known reduction, which leads to minimal D=5 N=2 gauged supergravity, to also include a multiplet of massive fields, containing the breathing mode of the SE_5. We show that the resulting D=5 action is also consistent with N=4 gauged supergravity coupled to two vector multiplets. This theory has a supersymmetric AdS_5 vacuum, which uplifts to the class of supersymmetric AdS_5 x SE_5 solutions, that spontaneously breaks N=4 to N=2, and also a non-supersymmetric AdS_5 vacuum which uplifts to a class of solutions first found by Romans.

The Laser Interferometer Space Antenna (LISA) will detect and characterize $\sim10^4$ Galactic Binaries consisting predominantly of two White Dwarfs (WD). An interesting prospect within this population is a third object--another WD star, a Circumbinary Exoplanet (CBP), or a Brown Dwarf (BD)--in orbit about the inner WD pair. We present the first fully Bayesian detection and posterior analysis of sub-stellar objects with LISA, focusing on the characterization of CBPs. We used an optimistic astrophysically motivated catalogue of these CBP third-body sources, including their orbital eccentricity around the inner binary for the first time. We examined Bayesian Evidence computations for detectability, as well as the effects on the posterior distributions for both the inner binary parameters and the third-body parameters. We find that the posterior behavior bifurcates based on whether the third-body period is above or below half the observation time. Additionally, we find that undetectable third-body sources can bias the inner binary parameters whether or not the correct template is used. We used the information retrieved from the study of the CBP population to make an initial conservative prediction for the number of detectable BD systems in the original catalogue. We end with commentary on the predicted qualitative effects on LISA global fitting and Galactic Binary population analysis. The procedure used in this work is generic and can be directly applied to other astrophysical effects expected within the Galactic Binary population.

Using the pure spinor formalism in part I [1] we compute the complete tree-level amplitude of N massless open strings and find a striking simple and compact form in terms of minimal building blocks: the full N-point amplitude is expressed by a sum over (N-3)! Yang-Mills partial subamplitudes each multiplying a multiple Gaussian hypergeometric function. While the former capture the space-time kinematics of the amplitude the latter encode the string effects. This result disguises a lot of structure linking aspects of gauge amplitudes as color and kinematics with properties of generalized Euler integrals. In this part II the structure of the multiple hypergeometric functions is analyzed in detail: their relations to monodromy equations, their minimal basis structure, and methods to determine their poles and transcendentality properties are proposed. Finally, a Groebner basis analysis provides independent sets of rational functions in the Euler integrals.

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.