Daniel Joseph Oliver

NANOGrav Physics Frontiers Center Postdoctoral Fellow at Oregon State University

The Gravitational Wave Spectrum

My research spans two complementary regions of the gravitational wave spectrum: the nanohertz band where pulsar timing arrays like NANOGrav operate, and the millihertz band that the upcoming space-based observatory LISA will probe. Click on a detector below to jump to the related research section.

Gravitational wave spectrum showing detector sensitivity curves for NANOGrav, LISA, and aLIGO
Sensitivity curves for current and proposed gravitational wave detectors. Adapted from GWPlotter (Jeremy Baier branch); original from GWPlotter (Main branch).

Pulsar Timing Arrays

Pulsar timing arrays (PTAs) detect gravitational waves at nanohertz frequencies (10-10–10-7 Hz) by monitoring tiny variations in the arrival times of pulses from millisecond pulsars. The NANOGrav PTA consisted of 68 MSPs at the time of the 15 yr data release (G. Agazie et al. 2023), all of which were chosen for their pulse shape stability, rotational stability, and/or brightness. (Jacobson-Bell et al. 2026)

Diagram of a millisecond pulsar showing its rotation axis, magnetic field axis, radio beams, and light cylinder
Diagram of a millisecond pulsar. Figure credit: H. Thankful Cromartie.

Anisotropic Sensitivity

With strong evidence now in hand for a stochastic gravitational wave background, most likely sourced by inspiraling supermassive black hole binaries, one of the next big questions is how that background is distributed across the sky. A non-uniform population of binaries should imprint anisotropy on the background, and detecting it would be a strong hint of the underlying astrophysics.

I have recently been developing an analytical framework for the sensitivity of a PTA to anisotropic backgrounds, pulsar-by-pulsar, direction-by-direction. The formalism is described in full in Oliver et al. (2026b), and is released as an extension to the hasasia sensitivity package.

Radiometer effective sensitivity skymap S_eff for the 16-year SPTA
Per-pixel radiometer effective sensitivity Seffrad at the peak sensitivity frequency (f ≈ 8.54 nHz) for the 16-year SPTA. Lighter regions are more sensitive; sensitivity is best where pulsars (stars) are clustered. Figure 6a from Oliver et al. (2026b).

Customized Chromatic Noise Models

The frequency-dependent (chromatic) delays imprinted on pulse arrival times by the interstellar and interplanetary medium can be up to two orders of magnitude louder than the gravitational wave background itself. I work on building customized chromatic noise models for each pulsar in NANOGrav’s datasets, testing the significance of competing noise prescriptions on a per-pulsar basis.

After applying our chromatic models, we observe significant impacts on the inference of achromatic noise in 19 out of 67 pulsars, finding in several cases that a previously significant achromatic noise process can be partially or entirely described as chromatic. (Larsen, Baier, Oliver, et al. 2026)

We find increased evidence for the Hellings-Downs (HD) correlation signature of the stochastic gravitational wave background (GWB), with a Bayes factor of 1571 ± 14 for HD-correlations over a common uncorrelated red-noise process using a power-law model with 14 Fourier modes. We find this ∼8× increase in Bayes factor from Agazie et al. (2023a) is a result of improved noise mitigation. (Agarwal et al. 2026)

Bar chart of the favored customized chromatic noise model for each NANOGrav 15-yr pulsar, grouped by basis type
The final models after the model selection process are represented in this bar chart. The pulsars are split into three columns corresponding to the favored basis type. Figure 3 from Larsen, Baier, Oliver, et al. (2026).
Customized chromatic noise model results, NANOGrav 12.5-yr dataset
Figure 10 from Hazboun, Simon, Larsen, Baier, Oliver, et al. (2026): comparison of customized chromatic noise modeling (blue) against the standard noise approach at 5 (green) and 30 (red) frequencies. The chromatic models bring the recovered gravitational wave background results more in line with the common process.

Laser Interferometer Space Antenna

LISA, the planned space-based gravitational wave detector, will be sensitive between roughly 10-5 and 10-1 Hz.

LISA sensitivity curve with the primary LISA source classes in the frequency-amplitude plane
Illustration of the primary LISA source classes in the gravitational wave (GW) frequency-amplitude plane. Figure 2.2 from the LISA Definition Study Report (Red Book; Colpi et al. 2024).

Extreme Mass Ratio Inspirals

A prime LISA source is the extreme mass ratio inspiral (EMRI): a stellar-mass compact object (~10 M☉) spiraling into a massive black hole (105–108 M☉) at the center of a galaxy.

The Gravitational Wave Peep

EMRI orbits typically start out highly eccentric and only emit detectable radiation in brief bursts as the compact object swings past periapsis. Each one of these bursts is a gravitational wave “peep”, firmly in the LISA band even when the long-period orbit itself is not. A single peep is unlikely to be individually resolvable, but a cosmological population of them out to redshift z = 3 contributes to LISA’s signal confusion noise. My dissertation and follow-up work focuses on modeling that background and quantifying how much it can obscure otherwise-detectable LISA sources.

Full EMRI waveform from capture to merger
A fully modeled EMRI waveform from capture to merger using a numerical kludge code (a = 0.8 M, p = 119.999916 M, e = 0.9999981). Repeated bursts at periapsis give rise to the gravitational wave peep. From Oliver et al. (2024).
Stochastic background from gravitational wave peeps
Stochastic background from gravitational wave peeps under three different population assumptions, plotted against the LISA A and E channel sensitivity curves. From Oliver et al. (2026a).

Post-Mission EMRIs

A stellar-mass compact object can radiate in the millihertz band for years to decades before it plunges into its massive black hole, so a substantial fraction of the EMRIs that LISA can detect will still be inspiraling when the mission ends. We find that post-mission EMRIs are roughly one in four of the LISA-detectable population at ρ > 20 over the four-year mission, and roughly one in six over the ten-year extended mission. (Oliver et al. 2026c)

Four-year spectrogram of an EMRI that plunges at the end of the window
(a) Merges within 4 years
Four-year spectrogram of the same EMRI starting 20 years before plunge
(b) Merges at 20 years
Four-year spectrograms of the strain for one EMRI (with parameters M = 106 M☉, m = 10 M☉, a = 0.9, ef = 0.05, at redshift z = 1), but starting at two different times. Clearly, both four-year segments are basically in-band for LISA, but are qualitatively different. From Oliver et al. (2026c).