Improving Semi-Analytic Spin Precession with NITs
Michael LaHaye University of Guelph
Search results
3397 - 3408 of 16791 Results
Format results
-
-
Fast Self-Forced Inspirals into a Rotating Black Hole
Philip Lynch University College Dublin
-
A multi-mode time-domain surrogate model for gravitational wave signals from comparable to extreme mass-ratio black hole binaries
Tousif Islam University of Massachusetts Dartmouth
-
Kerr self-force via elliptic PDEs: Numerical methods (part 2)
Thomas Osburn State University of New York (SUNY)
-
Kerr self-force via elliptic PDEs: Background and theory (part 1)
Nami Nishimura State University of New York (SUNY)
-
Conformal numerical method for self force applications in the time domain
Lidia Joana Gomes Da Silva Queen Mary University of London
-
Discontinuous collocation methods and self-force applications
Charalampos Markakis Queen Mary University of London
-
Self force review
Maarten van de Meent Max Planck Institute for Gravitational Physics - Albert Einstein Institute (AEI)
-
-
Time symmetry in operational theories
Lucien Hardy Perimeter Institute for Theoretical Physics
-
Emergent cosmology from tensorial group field theories for quantum gravity
Daniele Oriti Complutense University of Madrid
-
Exotic MOTS, the stability operator and their role in black hole mergers
Ivan Booth Memorial University of Newfoundland
-
Improving Semi-Analytic Spin Precession with NITs
Michael LaHaye University of Guelph
Semi-analytic solutions are useful because they are much faster than full numerical evolutions by virtue of the fact that they do not have to use as many points to achieve similar levels of accuracy. Currently there exists a semi-analytic solution for spin precessing binaries, which is implemented in LIGO and used to generate waveforms for comparison with gravitational waves. This solution comes with a caveat: it was calculated using precession averaged equations and thus has an oscillating error associated with the unaccounted for precession. This error is large enough that it can’t be overlooked for second and third generation detectors. By using near identity transforms (NITs) to reintroduce the effects of precession to the evolution, we can lower the error significantly. When implementing the NIT we found that the phase of oscillations we introduce (which coincides with the precession phase) does not line up with the precession in the numerical evolution, this causes the NIT to be less effective at reducing the error than it could be. To fix this issue we also introduce corrections to the evolution of the precession phase that were previously overlooked. -
Fast Self-Forced Inspirals into a Rotating Black Hole
Philip Lynch University College Dublin
Analysing the data for the upcoming LISA mission will require extreme mass ratio inpsiral (EMRI) waveforms waveforms that are not only accurate but also fast to compute and extensive in the parameter space. To this end, we present a method for rapidly calculating the inspiral trajectory of EMRIs with a spinning primary. We extend the work of van de Meent and Warburton (2018) by applying the technique of near-identity (averaging) transformations (NITs) to the osculating geodesic equations for a rotating (Kerr) black hole, resulting in equations of motion that do not explicitly depend upon the orbital phases. This allows us accurately to calculate the evolving constants of motion, orbital phases and waveform phase to within subradian accuracy, while dramatically reducing computational cost. We have implemented this scheme with an interpolated gravitational self-force model in both the equatorial and the spherical cases as a proof of concept, and present the first inspirals in Kerr spacetime to include all first order self-force effects. -
A multi-mode time-domain surrogate model for gravitational wave signals from comparable to extreme mass-ratio black hole binaries
Tousif Islam University of Massachusetts Dartmouth
We present EMRISur1dq1e6, a reduced-order multi-mode time-domain surrogate model of gravitational waveforms for non-spinning black hole binary systems with comparable- to extreme mass-ratio configurations. This surrogate model is trained on waveform data generated by a point-particle black hole perturbation theory (ppBHPT) framework computed from a high-performance Teukolsky equation solver code. In the comparable mass-ratio regime, the gravitational waveforms generated through ppBHPT agree surprisingly well with those from full numerical relativity after scaling of the ppBHPT’s total mass parameter. This model extends the EMRISur1dq1e4 waveform model, which spans 13,500M in duration and includes modes only up to (l,m)=(5,5). EMRISur1dq1e6, on the other hand, covers mass ratios from 3 to 1,00,000, can generate waveforms of duration up to 350,000M, and includes several spherical harmonic modes up to (l,m)=(10,10). The accuracy of training data is further improved by employing an updated plunge model in the ppBHPT framework. EMRISur1dq1e6 surrogate model has been extended to enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and events of interest to the future observatories such as Einstein Telescope, Voyager and Cosmic Explorer. -
Kerr self-force via elliptic PDEs: Numerical methods (part 2)
Thomas Osburn State University of New York (SUNY)
I will discuss the numerical methods we use to calculate the self-force on a scalar charge orbiting a Kerr black hole. We apply a 2nd-order finite difference scheme on a rectangular grid in the r*-θ plane. By working in the frequency domain and separating the ϕ variable (but not θ) we encounter elliptic PDEs, which present certain numerical challenges. One challenge is that every grid point is coupled to every other grid point so that a simultaneous solution requires solving a large linear system. Another related challenge involves how imperfect boundary conditions can introduce errors that would inevitably pollute the entire domain. We have applied various techniques to overcome these challenges, such as analyzing the boundary behavior to impose more sophisticated boundary conditions with improved accuracy (for a fixed outer boundary position). Various preliminary self-force results will be presented and discussed. -
Kerr self-force via elliptic PDEs: Background and theory (part 1)
Nami Nishimura State University of New York (SUNY)
Our long-term goal is to calculate the Lorenz gauge gravitational self-force for an extreme mass-ratio binary system in Kerr spacetime. Past work in the time-domain has encountered time instabilities for the two lowest modes m=0 and m=1. In order to overcome this problem, we enter the frequency-domain, which introduces elliptic PDEs. To develop an appropriate scheme, we first investigate the scalar self-force in Kerr spacetime by separating the Φ and t variables. To calculate the self-force, we use the effective source method. This presentation discusses the background and theory, which will be followed by a discussion of numerical methods in a later presentation. -
Conformal numerical method for self force applications in the time domain
Lidia Joana Gomes Da Silva Queen Mary University of London
"In 2034 LISA is due to be launched, which will provide the opportunity to extract physics from stellar objects and systems that would not otherwise be possible, among which are EMRIs. Unlike previous sources detected at LIGO, these sources can be simulated using an accurate computation of the gravitational self-force, resulting from the gravitational effects of the compact object orbiting around the massive BH. Whereas the field has seen outstanding progress in the frequency domain, metric reconstruction and self-force calculations are still an open challenge in the time domain. Such computations would not only further corroborate frequency domain calculations/models but also allow for full self-consistent evolution of the orbit under the effect of the self-force . Given we have a priori information about the local structure of the discontinuity at the particle, we will show how we can construct discontinuous spatial and temporal discretizations by operating on discontinuous Lagrange and Hermite interpolation formulae and hence recover higher order accuracy. We will show how this technique in conjunction with well-suited conformal (hyperboloidal slicing) and numerical (discontinuous time symmetric ) methods can provide a relatively simple method of lines numerical recipe approach to the problem. We will show, in particular, how this method can be applied to solve the Regge-Wheeler and Zerilli equations with a moving particle source in the time domain. /************ Note to organizers: if both my talk and my supervisor, Charalampos Markakis, are selected could this please be after his ? Thank you for your consideration *************/" -
Discontinuous collocation methods and self-force applications
Charalampos Markakis Queen Mary University of London
Numerical simulations of extereme mass ratio inspirals face several computational challenges. We present a new approach to evolving partial differential equations occurring in black hole perturbation theory and calculations of the self-force acting on point particles orbiting supermassive black holes. Such equations are distributionally sourced, and standard numerical methods, such as finite-difference or spectral methods, face difficulties associated with approximating discontinuous functions. However, in the self-force problem we typically have access to full a-priori information about the local structure of the discontinuity at the particle. Using this information, we show that high-order accuracy can be recovered by adding to the Lagrange interpolation formula a linear combination of certain jump amplitudes. We construct discontinuous spatial and temporal discretizations by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient method of solving time-dependent partial differential equations, without loss of accuracy near moving singularities or discontinuities. This method is well-suited for the problem of time-domain reconstruction of the metric perturbation via the Teukolsky or Regge-Wheeler-Zerilli formalisms. Parallel implementations on modern CPU and GPU architectures are discussed. -
Self force review
Maarten van de Meent Max Planck Institute for Gravitational Physics - Albert Einstein Institute (AEI)
-
-
Time symmetry in operational theories
Lucien Hardy Perimeter Institute for Theoretical Physics
The standard operational probabilistic framework (within which Quantum Theory can be formulated) is time asymmetric. This is clear because the conditions on allowed operations include a time asymmetric causality condition. This causality condition enforces that future choices do not influence past events. To formulate operational theories in a time symmetric way I modify the basic notion of an operation allowing classical incomes as well as classical outcomes. I provide a new time symmetric causality condition which I call double causality. I apply these ideas to Quantum Theory proving, along the way, a time symmetric version of the Stinespring extension theorem using double causality. I also propose the idea of a conditional frame of reference. We can transform from the time symmetric frame of reference to a forward or a backward frame of reference. This talk is based on arXiv:2104.00071.
-
Emergent cosmology from tensorial group field theories for quantum gravity
Daniele Oriti Complutense University of Madrid
We introduce the basic elements of tensorial group field theories (TGFTs) for quantum gravity, emphasizing how they encode quantum geometry and their relation with canonical loop quantum gravity and spin foam models. Next, we discuss briefly the issue of continuum limit and how it could be understood in this framework.
In the bulk of the talk, we overview the work on the extraction of an effective cosmological dynamics from TGFTs, inspired by the idea of the universe as a quantum gravity condensate. In this context, we emphasize: the need for appropriately coarse-grained states capturing collective dynamics and the role of relational observables and their construction. We discuss what the theory says (so far) about the fate of the big bang singularity at the beginning of our universe and how it suggests a quantum gravity origin for (phantom) dark energy.
-
Exotic MOTS, the stability operator and their role in black hole mergers
Ivan Booth Memorial University of Newfoundland
In the last couple of years it has been demonstrated that black hole spacetimes contain many more marginally outer trapped surfaces (MOTS) than had been previously recognized. For example, there is an infinite family of axially symmetric MOTS even in the Schwarzschild solution, of which the apparent horizon is only the first element. In a recent series of papers (arXiv: 2104.10265, 2104.11343, 2104.11343) we demonstrated that these exotic new MOTS play a key role in black hole mergers and, in fact, are the missing pieces needed to complete the apparent horizon “pair of pants” diagram. This merger turns out to be far richer than that of event horizons and, staying with clothing analogies, is better represented by a rococo ball gown than a pair of pants.
In this talk, I will overview the techniques that we used to find these surfaces, explain the role played by the exotic MOTS in resolving the merger of apparent horizons during a (head-on, non-rotating) black hole merger and also show how the stability operator brings order to what initially appears to be a chaotic melee of (often self-intersecting!) MOTS as they create, annihilate and weave through time. The stability operator was initially introduced by Andersson, Mars and Simon to characterize when a MOTS will smoothly evolve in time, but I will show that it can also be fruitfully understood as the Jacobi operator for MOTS.