Strong Gravity

Superradiant instabilities may create clouds of ultralight bosons around rotating black holes, forming so-called "gravitational atoms". The presence of a binary companion can induce transitions between bound states of the cloud ("resonances"), as well as from bound to unbound states ("ionization"). These processes backreact on the binary dynamics and leave characteristic imprints on the emitted GWs, carrying direct information about the mass of the boson and the state of the cloud. Even though some of the resonances can destroy the cloud before the system becomes observable in GWs, their backreaction forces the binary inclination towards attractors, opening up the possibility to infer the existence of the boson from a statistical analysis of a population of binary black holes. We also discuss implications for eccentric orbits and merger rates

Extreme mass-ratio binary black hole systems, known as EMRIs, are expected to radiate low-frequency gravitational waves detectable by planned space-based Laser Interferometer Space Antenna (LISA). We hope to use these systems to probe black hole spacetimes in exquisite detail and make precision measurements of supermassive black hole properties. Accurate models using general relativistic perturbation theory will allow us to unlock the potential of these unique systems. Such models must include post-geodesic corrections, which account for forces driving the smaller black hole away from a geodesic trajectory. When a spinning body orbits a black hole, its spin couples to the curvature of the background spacetime, introducing post-geodesic correction called the spin-curvature force. In this talk, I will present our calculation of EMRI waveforms that include both spin-curvature forces and the leading backreaction due to gravitational radiation. We use a near-identity transformation to eliminate dependence on the orbital phases, allowing for very fast computation of completely generic worldlines of spinning bodies; such efficiency is crucial for LISA data analysis. Finally, I will discuss what aspects still need to be included in future calculations so that we can use EMRIs for a new era of precision gravitational-wave astronomy, addressing outstanding puzzles in astrophysics, cosmology and fundamental theoretical physics.

---

Zoom link https://pitp.zoom.us/j/91917788358?pwd=MWp5OUhxbkRmZDFxWWE4cHR0VlBTUT09

A high specific abundance of millisecond radio pulsars has been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. In this talk, I will first demonstrate how we model millisecond pulsars in GCs using realistic cluster simulations. I will show the importance of electron-capture supernovae for neutron star retention, and how millisecond pulsar formation is greatly enhanced through dynamical interaction processes. I will also present some latest results on isolated millisecond pulsars, which are especially intriguing given the fact that millisecond pulsars are descendants of binary star systems. I will demonstrate the potential formation channels of isolated millisecond pulsars, some of which may also link to the formation of magnetars and the newly discovered fast radio bursts in a GC.

---

Zoom link https://pitp.zoom.us/j/97622593487?pwd=SHNoM1o3T1JjWVROTkJoZ0NWYmdyQT09

We investigate perturbations of the Schwarzschild geometry using a linearization of the Einstein vacuum equations within a Bondi-Sachs, or null cone, formalism. We develop a numerical method to calculate the quasinormal modes, and present results for the case ℓ= 2. The values obtained are different than those of a Schwarzschild black hole, and we interpret them as quasinormal modes of a Schwarzschild white hole.

---

Zoom link https://pitp.zoom.us/j/93648313308?pwd=akZMZXFKejIwdFphOVU0ZjFpeW13dz09

If axion Dark Matter exists, then they can collapse to form self-gravitating exotic compact objects known as axion stars. As mergers of such compact objects can potentially yield detectable gravititational waves which are correlated with a burst of electromagnetic radiation, much effort have been expended to compute these signals. I will discuss both the technical and theoretical challenges of this endeavour, and demonstrate such decays. I will show that axion stars may not be stable to electromagnetic decay, raising the question on whether we should expect to see these objects in the first place.

--- 

Zoom link https://pitp.zoom.us/j/98298572209?pwd=WVFZUFJqQzZQdU8vcjhQRVpzVHZDdz09

Cosmic-ray-driven instabilities play a crucial role in particle acceleration at shocks and during the propagation of GeV cosmic rays in galaxies and galaxy clusters within the self-confinement picture of CR transport. These instabilities amplify magnetic fields, which, in turn, scatter cosmic rays and thus self-regulate their transport. This leads to a strong coupling between the collisionless cosmic ray population and the thermal background plasma, implying potentially significant dynamic feedback. In this presentation, I discuss a recent discovery of a new cosmic ray-driven instability, referred to as the intermediate-scale instability, which triggers comoving ion-cyclotron electromagnetic waves at sub-ion skin-depth scales. Its growth rate is notably faster compared to the ion gyro scale (streaming) instability, which is commonly assumed to be the dominant instability in the self-confinement picture. Therefore, this new instability could play a vital role in the transport of cosmic rays in galactic and stellar environments. I then explore the implications of this instability for electron acceleration at non-relativistic shocks. Through Particle-in-cell (PIC) simulations, it is demonstrated that the new instability triggers the dominant mechanism for efficient electron acceleration at parallel electron-ion shocks, addressing a persistent issue with electron injection at these shocks. The PIC simulations also reveal that the common practice of using reduced ion-to-electron mass ratios in shock simulations, which artificially suppresses the intermediate instability, not only hinders electron acceleration but also leads to incorrect electron and ion heating in downstream and shock transition areas.

--- 

Zoom link https://pitp.zoom.us/j/91367222746?pwd=REpEdXE3ZGdLeVQ1bnh0NldIQktWQT09

In general relativity, the remnant black hole of a binary black hole merger emits “ringdown” radiation: a superposition of quasinormal modes with characteristic complex frequencies. The detection of more than one of these frequencies would serve as smoking-gun evidence of a black hole and would enable a test of the no-hair theorem. In this talk, I will discuss strategies for identifying quasinormal modes within simulated ringdown waveforms both in linear perturbation theory and full general relativity. I will show that a rich spectrum of modes could exist in the ringdown, including overtones, retrograde modes and nonlinear modes, and that interesting quasi-universal relationships exist between some of their amplitudes. I will also discuss the spectral instability of quasinormal modes and its potential imprints on the ringdown waveform. These results could guide the analysis of real ringdown signals detected by gravitational-wave detectors.

---

Zoom link https://pitp.zoom.us/j/92002625163?pwd=bHNDVVZKYnVNTENLNXJBLzVNKzVSQT09

Anticipating the launch of several next-generation gravitational wave (GW) detectors in the 2030s, we will be able to more precisely measure spacetime ripples from binary black hole (BH) mergers in a larger parameter space. The forthcoming data will require us to develop more accurate predictions of GWs not only in General Relativity (GR) but also in theories beyond GR and diverse astrophysical environments. Black hole perturbation theory is a cornerstone for making these predictions. In recent years, there have been extensive studies of perturbations of BHs in theories beyond GR, but only for non-rotating or slowly rotating BHs. In this talk, I will present a new formalism, based on Teukolsky's seminal work in the 1970s, to study perturbations of BHs with arbitrary spin in beyond-GR theories and in more complicated astrophysical environments. I will first discuss how to derive a modified Teukolsky equation for BHs deforming perturbatively from their counterparts in GR due to beyond-GR or environmental effects and the necessary techniques to evaluate this equation. Subsequently, I will discuss some applications of this formalism. Specifically, I will prescribe utilizing this formalism to investigate the isospectrality breaking of quasinormal modes (QNMs) in beyond-GR theories, compute the QNM frequency shifts in some specific theories, and efficiently extract these shifts from observation data. Furthermore, I will also show how to apply this formalism to study extreme mass-ratio inspirals beyond GR.

--

Zoom link: https://pitp.zoom.us/j/99282316326?pwd=REtBSFUxdlgxUGVwZFFvVEVBVnFTUT09

While Einstein's theory of gravity is formulated in a smooth setting, the singularity theorems of Hawking and Penrose describe many physical situations in which this smoothness must eventually breakdown. In positive-definite signature, there is a highly successful theory of metric and metric-measure geometry which includes Riemannian manifolds as a special case, but permits the extraction of nonsmooth limits under curvature and dimension bounds analogous to the energy conditions in relativity: here sectional curvature is reformulated through triangle comparison, while Ricci curvature is reformulated using entropic convexity along geodesics of probability measures.This lecture explores recent progress in the development of an analogous theory in Lorentzian signature, whose ultimate goal is to provide a nonsmooth theory of gravity. We highlight how the null energy condition of Penrose admits a nonsmooth formulation as a variable lower bound on timelike Ricci curvature. We discuss an example showing the null energy condition does not survive the same sort of limits that uniform bounds on the timelike Ricci curvature respect. 

---

Zoom link: https://pitp.zoom.us/j/92777399316?pwd=QmNaM1J0elgyRGdmaGxCTTQxT2Y4QT09

A stochastic gravitational wave background (SGWB) originates from numerous faint gravitational wave (GW) signals arising from coalescing compact binary objects. Based on the current estimated merger rate, the SGWB signal is expected to originate from non-overlapping GW waveforms where the chirping nature of individual events is expected to be preserved. In this talk, we present a novel technique, Spectrogram Correlated Stacking (or SpeCs), which goes beyond the usual cross-correlation (and to higher frequencies) by exploiting the higher-order statistics in the time-frequency domain. This method would account for the chirping nature of the individual events that comprise SGWB and enable us to extract more information from the signal due to its intrinsic non-gaussianity. We show that SpeCs improve the signal-to-noise for the detection of SGWB by a factor close to 8, compared to standard optimal cross-correlation methods which are tuned to measure only the power spectrum of the signal. SpeCs can probe beyond the power spectrum and its application to the GW data available from the current and next-generation detectors would speed up the SGWB discovery.

--

Zoom link: https://pitp.zoom.us/j/91002244803?pwd=a0dnMjZEYTEwSHBCVGRSeHB2Y2pJdz09

Two major approaches are used when numerically solving the Einstein field equations. The first one is to use spatial Cauchy slices and treat the system as a standard Cauchy initial value problem. Cauchy-characteristic evolution (CCE) serves as the second approach, which evolves spacetime based on null hypersurfaces. The Cauchy formulation is suitable for the strong field region but is computationally expensive to extend to the wave zone, whereas the Characteristic approach is fast in the wave zone but fails near the binary system where the null surfaces are ill-defined. By combining those two techniques — simulating the inner region with Cauchy evolution and the outer region with CCE, Cauchy-Characteristic matching (CCM) enables us to take advantage of both methods. In this talk, I present our recent implementation of CCM based on a numerical relativity code SpECTRE. I also discuss how CCM improves the accuracy of Cauchy boundary conditions — a benefit that allows us to evolve less of the wave zone in the Cauchy code without losing precision.

---

Zoom link https://pitp.zoom.us/j/98246275227?pwd=QWtmUDNkMlF6bXROLzBoYXVVTGpldz09

Merging neutron stars have long been speculated to produce Gamma-Ray Bursts (GRBs). While this has been the standard assumption for short GRBs for decades, this was only officially confirmed with the still relatively recent detection of gravitational wave counterpart GRB 170817A. In this presentation, I will discuss how we model short GRBs and their afterglows produced by the relativistic jets launched during the merger. A recent development of interest includes the apparent confirmation of some 'long' GRBs to also follow from mergers rather than collapsing massive stars, which has implication for the prospects of future multi-messenger detections. Joint modelling of GRB / afterglow data and GW signals makes it possible to further constrain the shared model parameters between the two, most notably the orientation of the system.

Zoom link: https://pitp.zoom.us/j/95551884959?pwd=bUVGVzhXK3pLb3MxMGV4ejZzUmRYdz09