Titles and Abstracts for Virtual Infinity Workshop 2024

Titles and abstracts complementing the schedule of the Virtual Infinity Workshop 2024 workshop that takes place online, July 1—2, 2024.

Pedro Baptista

Instituto Superior Técnico Lisbon

Particle creation with numerical hyperboloidal evolution

In quantum field theory particle pairs can be spontaneously excited out of the quantum vacuum by a dynamical spacetime or background. Well-known examples include Hawking radiation by collapsing stars, or particle creation by accelerated mirrors. Although this topic is well-understood from a theoretical viewpoint, only few explicit examples can be calculated in full closed form using analytical methods. In this work we explore this phenomenon using numerical techniques instead, and simulate the scattering of a massless scalar field from past null infinity to future null infinity with a set of effective potentials. As radiation (revealing the creation of particles) is only unambiguously defined at null infinity, this evolution naturally resorts to ingoing and outgoing hyperboloidal slices – spacelike slices that reach past and future null infinity, respectively. By comparing the field’s spectra there, we aim to ascertain the appearance of new frequency modes, which indicates the creation of particles. These numerical experiments are executed in spherical symmetry on a Minkowski spacetime background with a stationary, radially-dependent potential, as well as with a time-dependent potential that serves as a toy model and mimics the dynamics of the non-stationary regime of gravitational collapse. This presentation will report on the progress achieved so far.

Patrick Bourg

Radboud University

Quadratic quasi-normal mode dependence on linear mode parity

Quasi-normal modes (QNMs) uniquely describe the gravitational-wave ringdown of post-merger black holes. While the linear QNM regime has been extensively studied, recent work has highlighted the importance of second-perturbative-order, quadratic QNMs (QQNMs) arising from the nonlinear coupling of linear QNMs. Previous attempts to quantify the magnitude of these QQNMs have shown discrepant results. Using a new hyperboloidal framework, we show how the discrepancy can be resolved by showing that the QQNM/QNM ratio is a function not only of the black hole parameters, but also of the ratio between even- and odd-parity linear QNMs. In other words, the ratio QQNM/QNM depends on what created the ringing black hole, but the ratio of even to odd parity is the only dependence on the black hole pregenitor.

Carla Cederbaum

University of Tübingen

Geometric asymptotic foliations in initial data sets

Using asymptotic foliations as a tool to gain insights into the asymptotic behavior of initial data sets and in particular their charges has a long history in GR. The most well-known such foliations arise as coordinate spheres in suitable asymptotic coordinate charts. From a geometric and from a physical perspective, it is preferable to use asymptotic foliations of initial data sets which are themselves determined geometrically, e.g. as coordinate spheres in geometrically defined coordinate charts such as (spatially) harmonic coordinates or directly determined by a geometric condition such as Huisken—Yau’s constant mean curvature (CMC) foliations and the spacetime constant mean curvature (STCMC) foliation introduced by Sakovich and the speaker. The STCMC foliation is closely related to the center of mass charge of asymptotically Euclidean initial data sets.

In this talk, we will briefly present the state of the art of the STCMC foliation in asymptotically Euclidean initial data sets, in particular its relation to the center of mass, and then sketch open questions about the existence, uniqueness, and relevance of STCMC foliations in asymptotically hyperboloidal initial data sets.

Brad Cownden

University College Dublin

The Pseudospectrum of Black Holes in AdS

In this talk, based on the work in JHEP 05 (2024) 202, we explore the stability of the spectrum of black hole quasinormal modes (QNMs) in the context of planar Reissner-Nordström black holes in 5-dimensional Anti-de Sitter (AdS) spacetime. Our analysis encompasses the entire spectrum of gravitoelectric perturbations, employing ingoing Eddington-Finkelstein coordinates to derive a generalized eigenvalue system for these QNMs. While this approach produces a numerically convergent pseudospectrum, there are lingering issues regarding the invariance of the pseudospectrum under an overall rescaling of the eigenvalue problem. We will discuss how this issue poses challenges for the robustness of spectral stability analysis in physically relevant systems and how the requirement of invariance of the eigenvalue problem can limit the definition of the pseudospectrum.

Marina De Amicis

Niels Bohr Institute

Inspiral-inherited ringdown tails

Recently, studies on numerical evolutions of eccentric binary inspirals found a several orders of magnitude enhancement of the post-ringdown tail amplitude. This characteristic might render the tail a phenomenon of observational interest, opening the way to experimental verification of this general relativistic prediction in the near future.

Considering a source term describing an infalling test-particle in generic non-circular orbits, driven by post-Newtonian radiation reaction, I derive an integral expression over the system’s entire history, showing how the post-ringdown tail is inherited from the non-circular inspiral in a non-local fashion. Beyond its excellent agreement with numerical evolutions, the model explains the tail amplification with the progenitors’ binary eccentricity.

I will prove the tail to be a superposition of many power-laws, with each term’s excitation coefficient depending on the specific inspiral history. A single power law is recovered only in the limit of asymptotically late times, consistent with Price’s results and the classical soft-graviton theorem. I will conclude by discussing future directions, including the non-linear extension to comparable masses.

João Dinis Álvares

Instituto Superior Técnico Lisbon

Electrically charged numerical simulations on hyperboloidal slices

Gravitational wave radiation is only unambiguously defined at future null infinity: the location in spacetime where light rays arrive and where global properties of spacetimes can be measured. Within the context of numerical relativity we set up simulations reaching future null infinity by using hyperboloidal slices, as opposed to traditional Cauchy slices that reach spacelike infinity. Extending previous work in spherical symmetry, the Einstein-Maxwell-Klein-Gordon system is evolved on hyperboloidal slices, allowing to model gravity coupled to electromagnetism and a complex massless scalar field. This allows us to simulate the evolution of a charged scalar field and/or a Reissner-Nördstrom (electrically charged) black hole, where this last scenario serves as a useful toy model for a rotating (Kerr) black hole. We will report on current progress on these charged evolutions, where we retrieve their emitted signals at future null infinity as they would be seen by detectors on Earth.

Shalabh Gautam


Spherical Constraint Equations on Hyperboloidal Slices in Generalized Harmonic Gauge Formulation

Being a free evolution scheme, the only way to assure that constraints associated with the Einstein Field Equations in generalized harmonic gauge (GHG) are satisfied for all times is that the former are satisfied in the initial data (ID) up to their first-time derivative. This boils down to solving the system of Hamiltonian, momentum and GHG constraint equations in the ID. However, spherical polar coordinates and compactification make the equations singular at the origin and future null infinity. Nonlinearities introduce further complications. In this talk, I propose a method to regularize these equations and a choice of variables that make these equations linear in spherical symmetry and give the choice of variables a physical interpretation. Promising numerical results are obtained, offering a hope to extend this method in full 3d.

Benjamin Leather

Max Planck Institute for Gravitational Physics

Lorenz-Gauge Metric Perturbations with hyperboloidal slicing and spectral methods

Gravitational self-force theory is the leading approach for modelling gravitational wave emission from small mass-ratio compact binaries. This method perturbatively expands the metric of the binary in powers of the mass ratio. Traditional frequency domain approaches employ the variation of parameters method and compute the perturbation on constant time slices of the spacetime with numerical boundary conditions supplied at finite radius from series expansions of the asymptotic behaviour. This approach has been very successful but the boundary conditions calculations are tedious and the approach is not well suited to non-compact sources where homogeneous solutions must be computed at all radii. The latter, in particular, is important for emerging second-order in the mass ratio calculations. More recently, a new methodology has emerged based on approach where the spacetime is foliated by horizon-penetrating hyperboloidal slices - see Phys. Rev. D 105, 104033. Further compactifying the coordinates along these slices allows for simple treatment of the boundary conditions. In this talk, I shall present one of the first applications of these methods to gravitational perturbations in the Lorenz gauge. We find this method works efficiently for self-force calculations at first-order in the mass ratio and has some distinct advantages over the traditional approach. Finally, I shall discuss the importance of these methods for emerging second-order in the mass ratio calculations.

Philippe G. LeFloch

Sorbonne University

Optimal shielding for Einstein gravity

I will review recent advances about the construction of asymptotically Euclidean initial data sets for Einstein’s field equations. In collaboration with Bruno Le Floch (Sorbonne), I proved that solutions to the Einstein constraints can be glued together along conical domains. The solutions may have arbitrarily low decay at infinity, while enjoying (super-)harmonic estimates within possibly narrow (and nested) cones at infinity. Our results generalize earlier work by Alessandro Carlotto and Rick Schoen (sub-harmonic localization) and P. LeFloch and Nguyen (super-harmonic asymptotic localization). This presentation will be based on https://arxiv.org/abs/2312.17706 and https://arxiv.org/abs/2402.17598 Blog: philippelefloch.org

Michele Lenzi

Institute of Space Sciences (ICE-CSIC)

Integrable (hidden) structures in black hole perturbation theory

Perturbation theory of vacuum spherically symmetric spacetimes is a crucial tool for understanding the dynamics of black hole (BH) perturbations as well as BH scattering phenomena and the ringdown signal of binary BHs. Since the pioneering work of Regge and Wheeler it is known that the equations for the perturbations can be decoupled in terms of (gauge-invariant) master functions that satisfy 1 + 1 wave equations. However, the full landscape of master equations was recently found, clarifying that Einstein equations actually allows for an infinite set of them. These findings paved the way to the introduction of some new hidden symmetries governing the dynamics of perturbed non-rotating BHs: Darboux covariance and the infinite hierarchy of Korteweg-de Vries (KdV) deformations of the master equation. This generates a novel connection with integrable systems that relates the physical description of the perturbed BH, such as the greybody factors, to the KdV conserved quantities. The hyperboloidal approach gives a different point of view on these symmetries, promoting them to “bulk symmetries”. We explore this interplay between BH perturbation theory, integrable structures and hyperboloidal approach.

Christian Peterson Bórquez

CENTRA, Instituto Superior Técnico

Hyperboloidal evolution in spherical symmetry using the generalized harmonic formulation

In this talk I will discuss the use of the dual-foliation formalism to numerically evolve Einstein equations in the generalized harmonic formulation minimally coupled to a massless scalar field, restricted to spherical symmetry. I will show results for different types of initial data, with or without a black hole, all of which converge at the expected rate with increasing resolution. I will end by mentioning some undergoing ideas on how to deal with formally singular terms appearing when one wants to get O(1) variables all the way to future null infinity.

Anna Sancassani

University of Tübingen

Asymptotically corresponding initial data sets and associated geometric invariants

Motivated by the work of Marcus Khuri and Ye Sle Cha [1], we study the correspondence between geometric invariants of asymptotically AdS- hyperbolic and asymptotically hyperboloidal initial data sets. This is obtained from a detailed investigation of a correspondence first high- lighted by Piotr Chruściel and Paul Tod [2], which expresses the relationship between CMC initial data sets for the Einstein Equation under varying cosmological constants. This work is part of my ongoing PhD project and allows us to generalize the results cited above.

[1] Cha, Y.S., Khuri, M. Transformations of asymptotically AdS hyperbolic initial data and associated geometric inequalities. Gen Relativ Gravit 50, 3 (2018).

[2] Chruściel, Piotr T., and Paul Tod. ”An angular momentum bound at null infinity.” Advances in Theoretical and Mathematical Physics 13.5 (2009): 1317-1334.

Subhodeep Sarkar

Indian Institute of Technology Madras

The pseudospectrum of asymptotically de Sitter black holes

We will discuss the stability of the quasinormal mode spectrum of asymptotically de Sitter black holes. To capture the notion of spectral instability, we shall introduce and explain the idea behind the pseudospectrum of non-self-adjoint operators. Our findings will be based on a study of the spectrum and the pseudospectrum of asymptotically de Sitter black holes carried out numerically using spectral methods in a hyperboloidal foliation. Our analysis demonstrates how small deformations or perturbations to the scattering potential may cause the quasinormal modes to deviate from their unperturbed values. We shall compare our results with those obtained for asymptotically flat and anti-de Sitter black holes. We shall also briefly highlight recent results and open questions regarding the pseudospectrum of black holes.

Alex Vañó Viñuales

Instituto Superior Técnico Lisbon

Towards relating Cauchy and hyperboloidal initial data

A practical goal of non-linear hyperboloidal evolutions in numerical relativity is to extract (gravitational wave) signals at future null infinity. Assessment of the results will include comparisons with state-of-the-art methods to obtain signals there, currently Cauchy-characteristic evolution. The natural requirement is that the physical system under consideration be the same in both approaches. However, there is no direct way to relate initial data on a hyperboloidal slice to that on a Cauchy one, or vice versa, other than (numerical) evolution including future null infinity. I will illustrate the problem and sketch possible ways to overcome it.

Manas Vishal

UMass Dartmouth

(Towards) Fast and Accurate Simulation Technique for EMRIs in Space Telescopes

Laser Interferometer Space Antenna (LISA), a future space-borne gravitational wave detector, is primarily sensitive to Extreme Mass Ratio Inspirals (EMRIs), which are black hole mergers with a mass ratio greater than 100,000. We require an extremely precise template wave bank for the matched filtering procedure. In this talk, I will present a Discontinuous Galerkin (DG) technique for gravitational wave simulations of EMRI systems. The Teukolsky equation, which governs the behavior of EMRIs, is first reduced to a set of coupled 1+1D wave equations with a delta source term that functions as the secondary black hole in an EMRI system. Our DG approach, in contrast to previous numerical schemes, can precisely incorporate the smaller black hole’s point particle behavior in the form of a delta function. Compared to previous methods, our efficient method produces a very accurate waveform due to the scheme’s spectral (super)convergence property and exact treatment of the Dirac delta function that models the smaller black hole. Our time domain solver now includes hyperboloidal layers, which allow us to retrieve the solution at null infinity. We verify our computation by calculating Kerr and Schwarzschild energy fluxes from circular orbits at null infinity and Price tail power laws.

Network Meeting

The Infinity on a Gridshell 2023 meeting at the Niels Bohr Institute in Copenhagen marked the foundation of our Hyperboloidal Research Network. In this talk, we will present a brief update on our activities and achievements over this first year (monthly Virtual Infinity Seminars, Hyperbolic Times newsletter, YouTube channel, and topical collection). We plan to discuss networking options and strategic developments for the future, mainly focusing on fund applications for staff mobility for the network.