Research

Full list of publications can be found on Google Scholar, arXiv, ORCID and Inspire.

Gravitational waves from Gyroscopes inspiraling into gargantuan black holes

The study of spinning bodies/gyroscopes moving in curved spacetime is important for understanding binary black hole systems with extreme mass ratios. As the smaller black hole orbits the larger one, the system emits gravitational waves. The backreaction due to this gravitational-wave emission leads to the inspiral of the small body into the black hole. The spin of the small black hole modifies the waveform associated with the inspiral and must be included in our models in order to have sufficiently accurate templates for low frequency gravitational wave detectors such as LISA.
In our paper arXiv:2305.08919, we compute generic inspirals of spinning bodies and their associated waveforms (omitting several key aspects). We lay out a framework which we hope to serve as a foundation for including secondary spin effects in large-mass-ratio waveform models. Click for more details.

Extreme mass-ratio binary black hole systems 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 the spacetime in exquisite detail and make precision measurements of the larger black hole’s properties. Accurate models using general relativistic perturbation theory will allow us to realize the potential of these large mass-ratio systems. Such models must include post-geodesic corrections, including the backreaction due to gravitational-wave emission that leads to the inspiral of the small body into the black hole. When a spinning body orbits a black hole, its spin couples to the curvature of the background spacetime. This introduces a second post-geodesic correction called the spin-curvature force. In our paper arXiv:2305.08919, we calculate spinning-body inspirals and associated waveforms that include both spin-curvature forces and the leading gravitational wave backreaction.  Aspects of the self force have been neglected, and these must be included in future work. We build a framework using an osculating geodesic formulation combined with a near-identity transformation to eliminate dependence on the orbital phases, allowing for very fast computation of completely generic worldlines.

Precisely computing orbits of spinning bodies around black holes

Binary black hole systems with extreme mass ratios are expected to be important sources for low frequency gravitational wave detectors, such as LISA. In such systems, the smaller black hole can be treated as a test body moving in the background spacetime of the larger black hole. At lowest order, the small black hole moves along geodesics of the background spacetime. However, astrophysical black holes will in general be spinning! Therefore, we need to include a post-geodesic effect called spin-curvature coupling, which is the coupling of the small body’s spin to the larger black hole’s spacetime curvature.
In our papers arXiv:2201.13334 and arXiv:2201.13335, we present a frequency-domain approach for precisely characterizing orbits of spinning bodies in the presence of spin-curvature coupling. Click for more details.

Post-geodesic effects such as spin-curvature coupling must be included in models of extreme mass-ratio binary black hole systems in order to build sufficiently accurate templates for LISA. Exploiting the fact that in the large mass-ratio limit spinning-body orbits are close to geodesics, we develop a frequency-domain formulation of the motion which can be solved precisely. We examine a range of orbits with this formulation. We investigate orbits which are eccentric and nearly equatorial but for which the small body’s spin is arbitrarily oriented (in DOI:10.1103/PhysRevD.105.124040, arXiv:2201.13334) and we also discuss generic orbits with general small-body spin orientation (in DOI:10.1103/PhysRevD.105.124041, arXiv:2201.13335). We characterize the behavior of these orbits and show how the small body’s spin shifts the frequencies which affect orbital motion. These frequency shifts change accumulated phases which are direct gravitational-wave observables, illustrating the importance of precisely characterizing these quantities for gravitational-wave observations.

Predicting pulsar glitches using a state-dependent Poisson process

On the left, a magnetized, rotating neutron star, called a pulsar, emits regular pulses of electromagnetic radiation. The pulses can be observed when a beam of emission from the neutron star’s magnetic poles sweeps over the Earth, like a lighthouse. Therefore, these pulses allow us to deduce the rotational frequency of the pulsar, which typically decreases over time. Occasionally a “glitch” interrupts this decrease; a sketch of a canonical pulsar glitch is shown on the right.
In our paper, arXiv:1910.05503, we use a state-dependent Poisson process to model pulsar glitches, enabling us to predict the epoch of the next glitch for three different pulsars. Click for more details.

Glitches in some pulsars display power-law size and exponential waiting time distributions. These statistics are consistent with a state-dependent Poisson process, where the glitch rate is an increasing function of a global stress variable, which in this case is the angular velocity lag between the pulsar’s crust and the superfluid in its interior. In our paper (DOI:10.3847/1538-4357/ab44c3, arXiv:1910.05503), we estimate the parameters for this model for three pulsars (PSR J1740−3015, PSR J0534+2200, and PSR J0631+1036) and predict the epochs for each of their next glitches according to our model.

A Neutron star super-mixture of Interlinked vortex and flux tube arrays

On the left, we show the stratified interior of a neutron star. In the outer core (shown in purple), there are superfluid neutrons and superconducting protons, which are threaded by arrays of quantum vortices (red shading) and magnetic flux tubes (dark blue shading) respectively. Broadly, the superfluid vortices tend to align with the rotation axis and the flux tubes with the magnetic axis; misalignment between the two axes leads to frustration in the system and ‘glassy’ behaviour. The vortex array tangles as a consequence of coupling to the misaligned flux tube array, as shown in the inset.
In our papers, arXiv:1709.02254 and arXiv:1712.02938 we investigate the tangling of interlinked superfluid vortex and superconducting flux tube arrays. Click for more details.

The outer core of a neutron star contains two interpenetrating fluids: superfluid neutrons and superconducting protons. The protons and rigid crust corotate, while the angular velocity of the neutrons is determined by the number and disposition of the superfluid vortices, each of which carries a quantum of circulation. The proton superconductor is type II, implying that the magnetic field is concentrated into flux tubes, each carrying a magnetic flux quantum. We investigate the complex microscopic interaction between neutron vortices and proton flux tubes in detail. We investigate the way in which the the vortex array rearranges and deforms geometrically in equilibrium (in DOI:10.1093/mnras/stx2301, arXiv:1709.02254) and under far-from-equilibrium conditions (in DOI:10.1093/mnras/stx3197, arXiv:1712.02938). For the idealized model presented in our papers, we find that an initially rectilinear vortex array bends macroscopically and tangles microscopically in certain regimes (forming a ‘vortex crystal’), challenging the conventional picture of the outer core of the neutron star.