Fact-checked by Grok 2 weeks ago

Pulsar timing array

A pulsar timing array (PTA) is a network of precisely monitored millisecond pulsars—rapidly rotating neutron stars that emit regular radio pulses—used to detect and characterize ultra-low-frequency gravitational waves in the nanohertz regime through high-precision measurements of pulse arrival times.^[1] These arrays typically involve observations of 50 to 100 pulsars distributed across the sky, enabling the identification of correlated timing residuals induced by passing gravitational waves.^[2] The methodology relies on the exceptional rotational stability of millisecond pulsars, which serve as interstellar clocks with timing precision rivaling atomic clocks on Earth.^[2] Gravitational waves, predicted by general relativity, stretch and squeeze spacetime, causing minute delays or advances (on the order of nanoseconds) in the arrival times of pulses from distant pulsars.^[3] By analyzing residuals after accounting for known astrophysical and instrumental effects, PTAs search for spatial correlations among pulsars that follow the distinctive Hellings-Downs curve, a quadrupolar pattern unique to gravitational waves originating from sources like supermassive black hole binaries.^[4] This approach probes frequencies inaccessible to ground-based detectors like LIGO, targeting waves with periods of years to decades.^[1] Major PTA efforts include the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the European Pulsar Timing Array (EPTA), the Parkes Pulsar Timing Array (PPTA), the Indian Pulsar Timing Array (InPTA), and the Chinese Pulsar Timing Array (CPTA), which collaborate under the International Pulsar Timing Array (IPTA) to pool data from over 100 pulsars worldwide.^[2] Observations are conducted using large radio telescopes such as the Green Bank Telescope, Effelsberg, and MeerKAT, with datasets spanning 10–15 years to build sensitivity to faint signals.^[5] These projects not only aim to detect a stochastic gravitational-wave background but also individual sources, while mitigating noise from interstellar medium effects, pulsar intrinsics, and terrestrial interference.^[1] In June 2023, NANOGrav announced compelling evidence for a stochastic gravitational-wave background using its 15-year dataset of 68 pulsars, with the signal exhibiting the expected Hellings-Downs correlation and spectral properties consistent with a cosmic population of supermassive black hole binaries. In January 2025, NANOGrav was awarded the Bruno Rossi Prize for this evidence.^[4] Concurrently, EPTA, PPTA, and CPTA reported similar evidence from their datasets, marking a milestone in multimessenger astronomy and opening avenues to probe galaxy evolution, dark matter, and early-universe cosmology.^[2] Ongoing analyses and future data releases are expected to further characterize the signal and potentially confirm its nature, with efforts continuing into 2025 and beyond.^[3]

Fundamentals

Pulsars and precision timing

Pulsars are rapidly rotating neutron stars, the remnants of massive stars that have exploded as supernovae, characterized by strong magnetic fields typically ranging from $10^{8} to $10^{12} gauss, with millisecond pulsars exhibiting lower fields around $10^{8} to $10^{10} gauss, and rotation periods ranging from milliseconds to seconds.^[6]^[7] These compact objects, with diameters of about 20 kilometers and masses around 1.4 solar masses, emit beams of electromagnetic radiation, primarily in the radio band, from regions near their magnetic poles. When the rotation axis is misaligned with the magnetic axis, the beams sweep across the sky like a lighthouse, producing observed pulses with periods equal to the neutron star's rotation period or a fraction thereof.^[8] The first pulsar was discovered in 1967 by Jocelyn Bell Burnell, a graduate student at the University of Cambridge, while analyzing data from a radio telescope designed to study interplanetary scintillation.^[9] This pulsating radio source, later identified as PSR B1919+21 with a period of 1.337 seconds, was reported in a seminal paper that confirmed pulsars as rotating neutron stars.^[9] In the 1970s and 1980s, as radio telescopes improved in sensitivity and observational techniques advanced, pulsar timing emerged as a key method to study these objects, enabling precise measurements of their rotational properties and environmental effects.^[10] Pulsar timing consists of repeatedly observing the arrival times of radio pulses and fitting a parameterized model to these times-of-arrival (TOAs) to determine the pulsar's intrinsic properties.^[10] The model typically includes the spin period P, its first derivative (spindown rate) \dot{P}, the pulsar's sky position, proper motion, and corrections for relativistic effects, as well as astrophysical influences like orbital motion in binaries.^[10] By comparing observed TOAs to predictions from this model, residuals are computed, revealing any deviations due to unmodeled effects or external perturbations.^[7] Achieving high-precision timing, essential for pulsar timing arrays, relies on selecting stable millisecond pulsars with rotation periods of 1–10 ms, which exhibit low spindown rates and minimal intrinsic noise due to their age and evolutionary history.^[11] Observations are referenced to atomic clocks, such as hydrogen masers, to achieve TOA uncertainties as low as 100 nanoseconds or better over long baselines.^[12] A critical correction is for the interstellar medium's dispersive delay, quantified by the dispersion measure (DM), which is the integrated electron density along the line of sight and varies slowly due to the medium's motion; this is modeled as a frequency-dependent delay \Delta t \propto \nu^{-2}, where \nu is the observing frequency, and fitted iteratively with multi-frequency data.^[10] A representative example is PSR B1937+21, the first millisecond pulsar discovered in 1982 using the Arecibo telescope, with a period of 1.5578 ms and no companion, indicating it was likely spun up by accretion in a prior binary phase.^[13] This pulsar's timing stability rivals that of terrestrial atomic clocks, achieving a frequency stability of at least $6 \times 10^{-14} over intervals longer than 4 months, with pulse TOAs precise to tens of nanoseconds after corrections.^[14] Such exceptional stability underscores the potential of millisecond pulsars as celestial clocks for advanced applications.

Gravitational wave sources for detection

Pulsar timing arrays (PTAs) are primarily sensitive to gravitational waves (GWs) in the nanohertz frequency band, spanning approximately $10^{-9} to $10^{-7} Hz, which corresponds to wavelengths on the order of light-years. This regime bridges the gap between higher-frequency detections by ground-based observatories like LIGO, which probe frequencies from about 10 Hz to 1 kHz, and ultra-low-frequency signals imprinted on the cosmic microwave background at frequencies below $10^{-16} Hz.^[15] The dominant anticipated source of GWs in this band is the stochastic background arising from a cosmological population of supermassive black hole binaries (SMBHBs), formed through mergers of galaxies that host supermassive black holes at their centers. These binaries, typically involving black holes with masses exceeding $10^8 solar masses, emit continuous GWs during their inspiral phase, with orbital periods ranging from years to centuries, producing GW frequencies in the PTA-sensitive range. The characteristic strain amplitude for individual nearby SMBHBs is expected to be around h \sim 10^{-15} at nanohertz frequencies, though most signals will be below direct detection thresholds, contributing instead to an unresolved stochastic background. The energy density of this background is estimated as \Omega_{\rm gw} \sim 10^{-15} at a reference frequency of $10^{-8} Hz, reflecting the integrated emission from billions of such systems across cosmic history.^[15]^[16] In addition to SMBHBs, other astrophysical and cosmological processes can generate detectable GW backgrounds in the nanohertz band. Cosmic strings—topological defects potentially formed in the early universe—produce bursts of GWs from cusps and kinks in their loops, leading to a broadband stochastic signal with a characteristic strain spectrum that scales as h_c \propto f^{-1}. First-order phase transitions in the early universe, such as those associated with electroweak symmetry breaking or QCD transitions, can bubble into existence and collide, sourcing a stochastic background with peak frequencies around nanohertz depending on the transition temperature. Stellar-mass binaries within dense environments like globular clusters may also contribute a subdominant stochastic component through their collective emission, though their signals are weaker and more localized compared to SMBHBs. Finally, a primordial GW background, relic radiation from inflation or other early-universe mechanisms, could manifest in this band if produced on superhorizon scales, offering probes of fundamental physics beyond the Standard Model.^[15]

Theoretical framework

Timing residuals and signal modeling

Timing residuals in pulsar timing arrays are defined as the differences between the observed times of arrival (ToAs) of pulsar pulses and the ToAs predicted by a parameterized timing model. This model fits for the pulsar's sky position, proper motion, spin frequency, spin-down rate, and binary orbital parameters if applicable, effectively removing the deterministic components of the pulse emission and propagation. The resulting residuals, denoted as r(t), encapsulate any unmodeled effects, including instrumental noise, astrophysical perturbations, and potential gravitational wave (GW) signals. These residuals are typically on the order of microseconds or less for millisecond pulsars, enabling sensitivity to nanohertz-frequency GWs.^[17] Gravitational waves induce timing residuals by perturbing the spacetime metric along the line of sight to the pulsar, causing a fractional frequency shift (redshift) z(t) = \Delta \nu / \nu in the received pulses, where \nu is the pulse frequency. The residual is then the time integral r(t) = \int_0^t z(t') \, dt'. For a GW with strain h_{ij}(t, \mathbf{x}), the redshift for a pulsar in direction \hat{\Omega} is z(t, \hat{\Omega}) = \frac{1}{2} \hat{\Omega}^i \hat{\Omega}^j \Delta \bar{h}_{ij}(t, \hat{\Omega}), where \Delta \bar{h}_{ij} is the difference between the transverse-traceless metric perturbation at Earth (\mathbf{x} = 0) and at the pulsar (\mathbf{x} = D \hat{\Omega}, with distance D), projected transverse to the propagation direction. This yields two distinct contributions: the Earth term, which affects the signal from emission to Earth and is coherent across all pulsars, and the pulsar term, which perturbs the signal from Earth to the pulsar at a retarded time t - D(1 - \cos \theta)/c (with \theta the angle between GW propagation and pulsar directions) and is unique to each pulsar, often treated as additional noise due to its lack of correlation. A simplified form for the residual induced by a plus-polarized GW (h_+) is

R(t) = \frac{1 - \cos \theta}{2} \int_0^t h_+(t' - D(1 - \cos \theta)/c) \, dt',

illustrating the integrated strain projection along the line of sight.^[17] Signal modeling in PTAs distinguishes between deterministic signals from individual sources and stochastic backgrounds. For monochromatic waves, typically from supermassive black hole binaries, the GW is modeled as a quasi-sinusoidal strain h_A(t) = \mathrm{Re} \{ \tilde{h}_A e^{-i 2\pi f t + i \phi} \} (for polarizations A = +, \times), with parameters including amplitude \tilde{h}_A, frequency f, sky location, inclination, and polarization angle; the induced residual becomes r(t, \hat{\Omega}) = \sum_A F^A(\hat{\Omega}, \psi) \tilde{h}_A e^{-i 2\pi f t} [1 - e^{-i 2\pi f \tau_p}] / (i 2\pi f), where F^A are antenna patterns, \psi the polarization angle, and \tau_p = D(1 - \hat{k} \cdot \hat{\Omega})/c the pulsar-term delay. Stochastic backgrounds, arising from an ensemble of unresolved sources, are characterized as isotropic Gaussian processes with a power-law strain spectrum h_c(f) = A (f / f_\mathrm{yr})^{-\gamma}, where A is the characteristic amplitude at reference frequency f_\mathrm{yr} = 1 yr^{-1} (typically $10^{-15} for supermassive black hole binaries) and \gamma = 13/3 for the expected steep spectrum; the residual power spectrum follows S_r(f) \propto f^{-5} h_c^2(f). Parameter estimation employs Fourier-domain decompositions to represent signals and noise as sums of sinusoids, combined with Bayesian inference to sample posteriors on A and \gamma, marginalizing over noise and timing model parameters using Markov chain Monte Carlo methods.^[17] Observed residuals also incorporate noise from multiple sources, modeled separately to isolate GW signals. White noise, uncorrelated between epochs, stems primarily from radiometer noise due to finite observing time and flux density, as well as intrinsic pulse phase jitter, with a flat power spectral density P_w(f) \approx \sigma^2 T / N, where \sigma is the ToA uncertainty (typically 50–200 ns), T the total observation span, and N the number of ToAs. Red noise, exhibiting power-law enhancement at low frequencies (P_r(f) \propto f^{-\beta}, \beta \approx 3–7), arises from the pulsar's intrinsic rotational irregularities (e.g., "timing noise" from crust-superfluid interactions) and interstellar medium effects, including dispersive delays from electron density fluctuations (\Delta \mathrm{DM}) and scattering broadening, the latter being chromatic and modeled via frequency-dependent power laws. These components are fit using Gaussian process regressions or Fourier bases in Bayesian frameworks to mitigate biases in GW searches.^[17]^[18]

Correlation patterns in pulsar arrays

In pulsar timing arrays, the detection of a stochastic gravitational wave background relies on identifying spatial correlations in the timing residuals across multiple pulsars, which exhibit a distinctive quadrupolar pattern known as the Hellings-Downs curve. This curve describes the expected covariance between residuals from pairs of pulsars as a function of their angular separation ζ on the sky, arising from the isotropic and unpolarized nature of the background.^[19] The correlation approaches approximately 1/3 for small angular separations (ζ ≈ 0°) and -1/3 for antipodal pairs (ζ = 180°), with the difference arising from the tensor transverse-traceless nature of gravitational waves, where perturbations can align or oppose depending on pulsar geometry.^[20] It reaches a peak positive correlation at approximately 60° separation, providing a unique signature that scales with the gravitational wave strain amplitude.^[21] The Hellings-Downs curve is mathematically expressed as

\mu(\zeta) = \frac{1}{3} \left[ 1 + \frac{3}{2} (1 - \cos \zeta) \ln \left( \frac{1 - \cos \zeta}{2} \right) - \frac{1}{4} (1 - \cos \zeta) \right],

where the expression emerges from averaging the gravitational wave-induced perturbations over all sky directions and polarizations, assuming a plane-wave approximation valid at nanohertz frequencies.^[17] The derivation originates in general relativity, where gravitational waves passing over Earth induce fractional frequency shifts in pulsar signals, manifesting as integrated residuals in arrival times. For a stochastic background, correlations arise predominantly from the "Earth term," the common perturbation at the observation site, while pulsar-specific terms decorrelate over long baselines; the quadrupolar nature stems from the tensor transverse-traceless gauge of gravitational waves.^[19]^[22] This pattern distinguishes gravitational wave signals from noise sources, as common noise like clock errors produces a monopole (uniform across all pairs), solar system ephemeris errors yield a dipole (depending on pulsar position relative to the ecliptic), and anisotropic backgrounds introduce higher multipoles beyond the smooth quadrupolar profile.^[20] Simulations of pulsar timing data incorporating injected stochastic backgrounds demonstrate that fitting the observed inter-pulsar correlations to the Hellings-Downs curve yields high signal-to-noise ratios (e.g., 3–4σ in recent analyses), serving as a definitive "smoking gun" for gravitational wave detection while rejecting noise-only models with Bayes factors exceeding 10^3.^[21]

Operational PTAs

North American Nanohertz Observatory for Gravitational Waves (NANOGrav)

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) was founded in October 2007 as a collaboration of astronomers from the United States and Canada, initially comprising 17 members and expanding to over 225 scientists across more than 90 institutions by 2023.^[23]^[24] The project leverages pulsar timing arrays to detect low-frequency gravitational waves in the nanohertz regime, primarily using the Green Bank Telescope (GBT) in West Virginia—a 100-meter fully steerable dish—and the Arecibo Observatory in Puerto Rico until its collapse in 2020.^[25]^[26] Following Arecibo's loss, observations shifted to the GBT as the main instrument, supplemented by facilities like the Karl G. Jansky Very Large Array (VLA) and the Canadian Hydrogen Intensity Mapping Experiment (CHIME) for targeted pulsar searches.^[25]^[27] By the early 2020s, NANOGrav routinely monitored approximately 68 millisecond pulsars, selected for their timing stability to form a sensitive galactic-scale interferometer.^[28] NANOGrav's instrumentation emphasizes high-precision observations to minimize noise in timing residuals. The collaboration employs ultra-wideband (UWB) receivers at the GBT, enabling simultaneous coverage from 0.7 to 4.0 GHz for improved sensitivity to faint pulsar signals and reduced interstellar dispersion effects.^[29] Pulsar timing data are processed using the TEMPO2 software package, which fits pulse arrival times to models accounting for astrophysical and instrumental parameters, achieving sub-microsecond precision in many cases.^[30] These advancements, including wideband recording and advanced digital backends like VEGAS at the GBT, have enhanced signal-to-noise ratios, allowing for the integration of data spanning over 15 years from 2004 to 2020.^[29]^[28] A pivotal milestone came with the release of NANOGrav's 15-year data set in June 2023, which provided compelling evidence for a stochastic gravitational-wave background through correlated timing residuals across 67 pulsars, consistent with the predicted Hellings-Downs correlation curve from general relativity.^[21]^[4] This signal, with a characteristic strain amplitude of A = 2.0 \times 10^{-15} at a reference frequency of 3 nHz and a steep power-law spectral index of \gamma = 13/3, suggests origins from a population of supermassive black hole binaries (SMBHBs) in merging galaxies, though alternative sources like cosmic strings remain possible.^[21]^[31] The analysis imposed stringent constraints on the SMBHB population, limiting the merger rate and energy density while ruling out certain models of binary evolution.^[31] Additionally, NANOGrav has contributed to pulsar discoveries and high-precision timing of systems like PSR J1906+0746, a relativistic binary pulsar originally found via Arecibo surveys but now integral to its array for testing strong-field gravity.^[32] Post-2023, the collaboration continues monthly observations with an expanded array, aiming for definitive detection of the background and individual sources.^[33] In January 2025, NANOGrav was awarded the American Physical Society's Bruno Rossi Prize for its pioneering work in pulsar timing array gravitational-wave detection.^[34] Funded primarily by the National Science Foundation (NSF) as a Physics Frontiers Center, NANOGrav received a $17 million grant in 2021 spanning five years to support operations, data analysis, and infrastructure upgrades.^[35]^[36] It integrates with the International Pulsar Timing Array (IPTA) by sharing data and methodologies for joint analyses, enhancing global sensitivity to gravitational waves without supplanting regional efforts.^[37]^[35] This collaborative framework has yielded over 366 publications and 22,800 citations by 2023, underscoring NANOGrav's role in advancing pulsar astronomy and multimessenger astrophysics.^[5]

International Pulsar Timing Array (IPTA)

The International Pulsar Timing Array (IPTA) was formed in 2010 to foster international collaboration among regional pulsar timing efforts, merging datasets from the European Pulsar Timing Array (EPTA), NANOGrav, Parkes Pulsar Timing Array (PPTA), Indian Pulsar Timing Array (InPTA), and other groups to create a unified global observatory. This consortium collectively times approximately 100 millisecond pulsars, leveraging observations from radio telescopes worldwide to achieve greater sensitivity for detecting nanohertz gravitational waves than any single array could provide alone.^[38] IPTA data handling relies on specialized joint pipelines that integrate heterogeneous datasets, addressing discrepancies arising from different telescopes, observing frequencies, and solar system ephemerides used in pulsar position modeling.^[39] These pipelines employ techniques such as fitting for inter-observatory time jumps and noise parameter adjustments to produce coherent timing residuals across the array, enabling robust analyses of correlated signals.^[40] Among its key outputs, the IPTA established the first joint upper limits on an isotropic stochastic gravitational wave background in 2016, using early combined data from 49 pulsars to constrain the amplitude at nanohertz frequencies to below levels expected from supermassive black hole binaries. In 2023, analyses of updated datasets revealed strong evidence for a common-spectrum process across multiple pulsars, confirming the presence of a nanohertz stochastic gravitational wave background with a significance approaching 3σ in the joint IPTA framework. The IPTA Development Program promotes the identification and precise timing of new millisecond pulsars to expand the array's coverage and improve sky localization capabilities for gravitational wave sources. For gravitational wave searches, the collaboration utilizes the ENTERPRISE software package, a Bayesian framework that models timing residuals including red noise, multipath propagation, and Hellings-Downs correlations to infer signal properties. IPTA governance features annual science meetings and student workshops to coordinate research and training, with open data policies implemented following the 2023 background evidence, making pulsar timing datasets publicly available to accelerate community-wide analyses and verification.^[38] The IPTA held its 2025 science meeting to advance collaborative efforts in pulsar timing and gravitational wave research.^[41]

Proposed and emerging PTAs

European Pulsar Timing Array (EPTA)

The European Pulsar Timing Array (EPTA) was established in 2008 through a collaboration of major European radio astronomy institutes, aiming to detect nanohertz-frequency gravitational waves via precise timing of millisecond pulsars.^[42] The consortium coordinates observations using five large radio telescopes: the 100-m Effelsberg Radio Telescope in Germany, the 76-m Lovell Telescope at Jodrell Bank Observatory in the United Kingdom, the 94-m equivalent Nançay Radio Telescope in France, the 64-m Sardinia Radio Telescope in Italy, and the Westerbork Synthesis Radio Telescope (a 94-m equivalent array) in the Netherlands.^[43] These facilities enable the monitoring of an array of approximately 50 millisecond pulsars, with the second data release (DR2) in 2023 providing high-precision timing data for 25 such pulsars spanning up to 24.7 years.^[44] A key milestone for the EPTA was the 2023 release of its 24.7-year dataset, which, through analysis of timing residuals, confirmed evidence for a stochastic gravitational wave background at nanohertz frequencies, marking a significant step toward direct detection of supermassive black hole binaries and other cosmic sources.^[45] The collaboration prioritizes high-cadence observing campaigns, often conducted via the Large European Array for Pulsars (LEAP), which coherently combines signals from the five telescopes to achieve sub-microsecond timing precision on select targets.^[46] The EPTA emphasizes international coordination in data sharing and analysis protocols, contributing uniquely to global efforts by developing customized extensions to the TEMPO2 pulsar timing software for handling noise models and gravitational wave searches.^[43] Notable achievements include stringent upper limits on cosmic string networks, with Bayesian analyses of DR2 data yielding log₁₀(Gμ) < -10.5 at 95% confidence for certain loop distribution models, constraining the dimensionless string tension Gμ to below approximately 3 × 10^-11.^[47] The millisecond pulsar PSR J1713+0747 stands out as a flagship stable source in EPTA observations, prized for its low noise and long-term timing stability that enables sensitive probes of gravitational wave signals at higher nanohertz frequencies.^[48] As a foundational member of the International Pulsar Timing Array (IPTA), the EPTA played a central role in the 2023 joint collaboration that combined datasets from multiple regional arrays to robustly characterize the stochastic gravitational wave background. In January 2025, the EPTA received the Royal Astronomical Society Group Award in Astronomy and an ERC Advanced Grant for its pioneering work in gravitational wave astronomy.^[49]

Parkes Pulsar Timing Array (PPTA) and future extensions

The Parkes Pulsar Timing Array (PPTA) is a long-term observational program based at the Parkes Observatory in Australia, where regular high-precision timing measurements of millisecond pulsars began in early 2005.^[50] The project targets over 20 millisecond pulsars selected for their stability and low noise, with the third data release in 2023 encompassing observations of 32 millisecond pulsars spanning up to 18 years, enabling detailed analyses of timing residuals.^[51] These data have been instrumental in refining pulsar timing models and probing astrophysical phenomena. Key contributions from the PPTA include establishing some of the earliest stringent upper limits on the stochastic gravitational wave background in the nanohertz regime, based on analyses of timing data from 2013 onward, which constrained the energy density parameter to below 10^{-15} relative to the cosmic critical density at frequencies around 10^{-8} Hz.^[52] Additionally, the PPTA has advanced pulsar timing precision through observations of PSR J0437-4715, achieving sub-nanosecond residual errors that highlight the limits of interstellar medium effects and instrumental noise in high signal-to-noise observations. Looking toward future extensions, the PPTA is poised to integrate with the Square Kilometre Array (SKA), which is expected to expand pulsar timing capabilities by discovering and timing over 1,000 millisecond pulsars by the 2030s, dramatically increasing array density and sensitivity. The SKA's enhanced flux sensitivity, approximately 10 times greater than current facilities like Parkes, will enable timing precisions below 100 nanoseconds for faint pulsars, boosting detection prospects for gravitational waves.^[53] Complementing this, the Chinese Pulsar Timing Array (CPTA), utilizing the Five-hundred-meter Aperture Spherical Telescope (FAST), has initiated plans for a dedicated network since 2019, with its first data release (DR1) in 2025 from observations of 57 millisecond pulsars, demonstrating its potential to contribute to global efforts.^[54] The PPTA played a pivotal role in the 2023 International Pulsar Timing Array (IPTA) detection of a nanohertz gravitational wave background, providing key southern sky data in the second data release that corroborated the common-spectrum signal across multiple collaborations. However, the PPTA's focus on southern hemisphere pulsars introduces a sky coverage bias, limiting angular resolution for certain gravitational wave sources; ongoing global integration through the IPTA and SKA aims to achieve uniform coverage by incorporating northern and equatorial arrays.

Observations and analyses

Data collection and processing

Data collection for pulsar timing arrays (PTAs) involves regular, high-precision observations of millisecond pulsars using large radio telescopes, typically conducted on monthly to bi-weekly cadences to achieve the necessary temporal resolution for detecting nanohertz gravitational waves.^[55] Observations are performed at radio frequencies between 1 and 2 GHz, with some arrays extending to lower (e.g., 327–430 MHz) or higher bands (up to 3–4 GHz) for multi-frequency coverage to mitigate interstellar effects.^[43]^[51] Integration times per pulsar session range from 10 to 80 minutes, depending on the telescope and backend, aiming for signal-to-noise ratios exceeding 100 to ensure timing precision.^[43]^[55] For example, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) uses the Green Bank Telescope and Very Large Array with up to 30-minute integrations, while the Parkes Pulsar Timing Array (PPTA) employs 64-minute sessions at the Parkes telescope.^[55]^[51] The data processing pipeline begins with raw voltage or filterbank data reduction using specialized software to generate pulse profiles and times of arrival (TOAs). Radio frequency interference (RFI) is excised through automated methods such as median filtering, spectral kurtosis detection, or tools like pfits_zapUWL and MeerGuard, followed by manual inspection to remove contaminated segments.^[43]^[51] Polarization calibration is applied using noise diode injections at the start of observations, with flux calibration via primary standards like PKS B1937-15 for wideband systems.^[51]^[55] Folding searches are conducted with DSPSR to align pulses using predicted spin-down models, producing dedispersed profiles via coherent dedispersion backends like ROACH-based systems. TOAs are then extracted using PSRCHIVE's pat routine or similar, fitting templates to observed profiles.^[43] Timing residuals are fitted with Tempo2 or PINT software, incorporating pulsar parameters, orbital ephemerides, and Earth orientation corrections to minimize phase inconsistencies.^[56] Bayesian extensions like TEMPONEST handle parameter estimation for complex models. Noise mitigation is essential to isolate potential gravitational wave signals from instrumental and astrophysical effects. Profile variations due to interstellar scattering or intrinsic pulsar changes are modeled through template matching and multi-frequency observations, enabling dispersion measure (DM) estimation by fitting frequency-dependent delays with a cold plasma model.^[43] Clock corrections account for telescope instabilities and jumps between observing sessions, often parameterized in Tempo2.^[51] White noise (e.g., radiometer and template-fitting errors) is parameterized with EFAC and EQUAD factors, while red noise from spin irregularities or DM variations is modeled as power-law spectra using tools like ENTERPRISE.^[55] Typical PTA datasets span over 15 years, yielding 10^4 to 10^5 residuals per pulsar with error bars around 100 ns, as seen in NANOGrav's 15-year set with 676,000 narrowband TOAs across 68 pulsars.^[55]^[51] Quality control emphasizes pulsar selection for long-term stability and geometric distribution. Pulsars are chosen based on timing stability better than 10^{-15} in spin frequency derivative, low scatter (≤1 μs in initial tests), and sky positions ensuring quadrupolar correlation sensitivity, typically 20–70 per array.^[55]^[43] For instance, the European Pulsar Timing Array (EPTA) selected 25 from 42 candidates by evaluating noise properties and gravitational wave detectability via coupling matrices.^[43] Ongoing monitoring flags pulsars with excessive red noise or profile evolution for exclusion or special modeling.^[51]

Key results and detections

In June 2023, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the European Pulsar Timing Array (EPTA), the Parkes Pulsar Timing Array (PPTA), and the International Pulsar Timing Array (IPTA) announced compelling evidence for a stochastic gravitational-wave background (GWB) at nanohertz frequencies, consistent with predictions from a population of supermassive black hole binaries (SMBHBs). These results, derived from pulsar timing datasets spanning 15 years for NANOGrav and comparable durations for the others, showed correlated timing residuals across multiple pulsars matching the expected quadrupolar spatial pattern. The significance levels ranged from approximately 2σ for PPTA to over 3σ for NANOGrav and EPTA, with Bayesian analyses favoring a GWB model over noise-only explanations by factors exceeding 10^3 in some cases.^[21]^[45]^[57] Analyses from the individual collaborations yielded similar characteristic strain amplitudes around (2.0–2.5) × 10^{-15} at a reference frequency of f = 3 × 10^{-8} Hz (≈1 yr^{-1}), with a power-law spectral index near γ = 13/3, as expected for an SMBHB-origin GWB; for example, NANOGrav reported A = 2.4 × 10^{-15}.^[57] No individual continuous gravitational-wave sources were detected, but stringent upper limits were established, such as sky-averaged 95% confidence limits on the strain h_{95} < 9.1 × 10^{-15} at 10 nHz from EPTA data. Additionally, these observations constrained alternative GWB sources, ruling out cosmic string models with tension parameter Gμ > 10^{-10} and setting upper bounds around Gμ ≲ 10^{-11} in favored scenarios. A fit of the observed inter-pulsar correlations to the Hellings-Downs curve confirmed the quadrupolar signature expected from a tensor-mode gravitational-wave background, with a p-value of approximately 0.03 favoring the gravitational-wave hypothesis over uncorrelated noise.^[57]^[58]^[59] Following the 2023 announcements, refined analyses on the 15-year datasets have sharpened the GWB spectral characterization, confirming the power-law form while tightening constraints on the amplitude. These updates reveal hints of anisotropy in the GWB, potentially arising from source clustering or shot noise in the SMBHB population, though current data provide only upper limits on anisotropic power at the 10-20% level relative to the isotropic component. In January 2025, the EPTA received the Royal Astronomical Society Group Achievement Award for its contributions to gravitational wave astronomy. Preparations for the IPTA third data release, incorporating ongoing observations, are underway as of 2025. No definitive individual sources or deviations from general relativity have emerged, underscoring the stochastic nature of the detected signal.^[60]^[61]

Applications and future prospects

Scientific implications

The evidence for a nanohertz gravitational-wave background by pulsar timing arrays provides crucial insights into the population of supermassive black hole binaries (SMBHBs), which are expected to form during the mergers of massive galaxies and dominate the observed signal through their incoherent superposition. These binaries, typically with total masses exceeding $10^8 M_\odot, offer a direct probe of galaxy merger rates across cosmic history, revealing that mergers occur efficiently on timescales of approximately 2.8 Gyr, with higher rates in dense environments that facilitate black hole pairing. The background implies a population of around $10^8 to $10^9 such binaries across the observable universe, contributing to the stochastic signal at frequencies of 2–30 nHz, and highlights black hole growth mechanisms where binaries exhibit masses slightly larger than previously anticipated from quasar luminosity functions, suggesting significant accretion post-merger.^[62]^[63]^[64] In cosmology, future resolutions of individual SMBHBs by pulsar timing arrays could serve as standard sirens, enabling model-independent measurements of the Hubble constant H_0 with precision comparable to current distance-ladder methods, potentially resolving the tension between early- and late-universe estimates by leveraging luminosity distances derived solely from gravitational-wave signals.^[65] Pulsar timing arrays complement space-based detectors like LISA in multi-messenger astronomy, as PTAs probe the low-frequency inspiral phase of massive SMBHBs ($10^9–$10^{10} M_\odot) while LISA targets higher-frequency mergers of lighter systems ($10^6–$10^8 M_\odot), with PTA constraints on merger rates predicting up to thousands of LISA detections over its mission lifetime and enabling joint searches for electromagnetic counterparts in active galactic nuclei. Pulsar searches could further identify counterparts to resolved sources, linking gravitational-wave events to host galaxy emissions.^[66]^[64] The 2023 NANOGrav evidence specifically implies the presence of approximately $10^3 nearby merging SMBHBs within a few hundred megaparsecs, consistent with the observed background amplitude and offering tests of general relativity in the strong-field regime through deviations in the gravitational-wave spectrum, such as broken power-law shapes that could indicate post-Newtonian corrections at the -2PN or -3PN level.^[62]^[67] As of 2025, ongoing analyses of extended datasets continue to support this evidence, with recent publications comparing results across PTAs reinforcing the signal's consistency without yet achieving definitive confirmation.^[68] Interdisciplinary applications extend to fundamental physics, where the background constrains variations in fundamental constants, such as \dot{G}/G \lesssim 10^{-5} yr^{-1} over redshifts z = 0.1–1.^[67]

Technological challenges and advancements

One of the primary technological challenges in pulsar timing arrays (PTAs) is the limited number of millisecond pulsars available for monitoring, currently around 100 across global efforts like the International Pulsar Timing Array (IPTA). This constraint restricts the array's ability to form a sufficiently dense network for robust gravitational wave detection, as fewer pulsars reduce the statistical power to distinguish correlated signals from noise. Additionally, the distribution of these pulsars, primarily concentrated within the Milky Way's disk, introduces sky coverage biases that unevenly sample the celestial sphere and limit sensitivity to sources in certain directions.^[69]^[70] At low frequencies, red noise from the interstellar medium (ISM), including turbulence and dispersion measure variations, dominates the timing residuals and masks potential gravitational wave signals. These effects arise from plasma scattering and multipath propagation through ionized ISM structures, producing chromatic noise that correlates with observing frequency and complicates data analysis. Such red noise processes, alongside intrinsic pulsar spin irregularities, often exceed white noise contributions on timescales of years to decades, posing significant hurdles for achieving the required precision in pulse arrival times.^[71]^[72] Advancements in receiver technology have addressed some of these issues, notably through the adoption of wideband systems that span broader frequency ranges to mitigate ISM-induced noise. For instance, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) implemented wideband timing in its 12.5-year dataset released in 2020, enabling simultaneous modeling of dispersive effects across 0.8–2.0 GHz and improving timing precision for 47 millisecond pulsars. Artificial intelligence and machine learning techniques are increasingly employed for noise subtraction, such as modeling nonstationary red noise components to enhance signal-to-noise ratios in PTA datasets. Efforts to improve clock stability, including the use of optical fiber links for synchronizing terrestrial time standards at observatories, further reduce instrumental errors in long-term monitoring.^[30]^[73] Looking ahead, the Square Kilometre Array (SKA) Phase 1, anticipated to begin operations around 2027, promises to expand PTA capabilities by discovering 700–900 new millisecond pulsars, potentially increasing the array size by an order of magnitude and enhancing sky coverage. Cryogenic receivers, which lower system noise temperatures, offer up to a 50% boost in sensitivity by improving signal detection thresholds for faint pulsars. Current PTA noise floors, characterized by root-mean-square timing residuals around 100 ns, limit sensitivity; future goals aim for sub-10 ns precision with arrays of 1000 pulsars by 2040 to reach gravitational wave strain sensitivities below 10^{-15}. Mitigation strategies include developing international pulsar catalogs for standardized data sharing and applying machine learning to profile fitting, which dynamically adjusts pulse templates to account for profile evolution and reduce fitting errors.^[74]

References

[1]
[1707.01615] Gravitational wave research using pulsar timing arrays
Jul 6, 2017 · A pulsar timing array (PTA) refers to a program of regular, high-precision timing observations of a widely distributed array of millisecond pulsars.
[2]
Pulsar Timing Arrays - Albert Einstein Institute
Pulsar Timing Arrays are searching for low-frequency gravitational waves by regularly observing many millisecond pulsars and analyzing the arrival times of ...
[3]
Evidence for a Gravitational-Wave Background - NANOGrav
NANOGrav found evidence for a gravitational-wave background by tracking pulsars, showing correlated timing deviations consistent with general relativity, and ...
[4]
The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave ...
Jun 29, 2023 · We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15 yr pulsar timing data set ...Abstract · Introduction · Checks and Validation · Discussion
[5]
NANOGrav
NANOGrav is an international collaboration exploring the low-frequency gravitational wave universe through radio pulsar timing.People · Low-frequency gravitational... · NANOGrav Collaboration · DataMissing: 2023 | Show results with:2023
[6]
Determining the rotation direction in pulsars | Nature Communications
Jul 19, 2019 · Pulsars are strongly magnetised rotating neutron stars. Because of rotation, pulsars emit two intense radiation beams. For a distant ...
[7]
Quasi-periodic sub-pulse structure as a unifying feature for radio ...
Nov 23, 2023 · Magnetars are highly magnetized rotating neutron stars that are predominantly observed as high-energy sources. Six of this class of neutron ...
[8]
Observation of a Rapidly Pulsating Radio Source - Nature
Observation of a Rapidly Pulsating Radio Source ... Unusual signals from pulsating radio sources have been recorded at the Mullard Radio Astronomy Observatory.
[9]
tempo2, a new pulsar timing package – II. The timing model and ...
tempo2 is a new software package for the analysis of pulsar pulse times of arrival. In this paper, we describe in detail the timing model used by tempo2.Missing: tutorial | Show results with:tutorial
[10]
Chapter 6 Pulsars
Pulsar timing is the regular monitoring of the rotation of the neutron star by tracking (nearly exactly) the arrival times of the radio pulses. The key point to ...
[11]
Milky Way Accelerometry via Millisecond Pulsar Timing
Apr 7, 2021 · The temporal stability of millisecond pulsars is remarkable, rivaling even some terrestrial atomic clocks at long timescales.
[12]
Pulsar Timing and Its Application for Navigation and Gravitational ...
Jan 17, 2018 · Pulsars are natural cosmic clocks. On long timescales they rival the precision of terrestrial atomic clocks.Pulsar Timing And Its... · 1.2 Pulsar Timing: Pulsars... · 3 Gravitational Wave...
[13]
A millisecond pulsar - Nature
Dec 16, 1982 · A fast pulsar, 1937 + 214, with a period of 1.558 ms in the constellation Vulpecula only a few degrees from the direction to the original pulsar, 1919+21.
[14]
Millisecond Pulsar PSR 1937+21: A Highly Stable Clock | Science
The stable rotation and sharp radio pulses of PSR 1937+21 make this pulsar a clock whose long-term frequency stability approaches and may exceed that of the ...
[15]
[2307.02376] Gravitational wave sources for Pulsar Timing Arrays
Potential gravitational wave sources include supermassive black hole binary mergers, first-order phase transitions, cosmic strings, domain walls, and large- ...
[16]
A Practical Theorem on Gravitational Wave Backgrounds - arXiv
The theorem relates the gravitational wave background spectrum to the source's energy spectrum and remnant density, allowing computation of the background's ...
[17]
[PDF] Principles of Gravitational-Wave Detection with Pulsar Timing Arrays
Dec 14, 2021 · This paper aims to review the science behind PTA, the basis of detection, the main scientific goals, and the current results. 2. Pulsar Timing ...
[18]
Identifying and mitigating noise sources in precision pulsar timing ...
Nov 6, 2020 · Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency ...INTRODUCTION · SIGNAL MODELS · RESULTS · CONCLUSIONMissing: seminal | Show results with:seminal
[19]
https://ui.adsabs.harvard.edu/abs/1983ApJ...265L..39H/abstract
[20]
Answers to frequently asked questions about the pulsar timing array ...
Aug 10, 2023 · Abstract page for arXiv paper 2308.05847: Answers to frequently asked questions about the pulsar timing array Hellings and Downs curve.
[21]
The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave ...
Jun 28, 2023 · We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set ...
[22]
Understanding the gravitational-wave Hellings and Downs curve for ...
Dec 3, 2014 · In this paper, we give a pedagogical discussion of the Hellings and Downs curve for pulsar timing arrays, considering simpler analogous scenarios involving ...Missing: original | Show results with:original
[23]
NANOGrav Collaboration
NANOGrav, or the North American Nanohertz Observatory for Gravitational Waves, was founded in October 2007 and has since grown to over 225 members at over 90 ...Missing: history formation
[24]
Scientists Find Evidence for Slow-Rolling Sea of Gravitational Waves
Jun 28, 2023 · NANOGrav was founded in October 2007 and has grown to more than 190 members at more than 70 institutions. In 2015, it was designated a ...
[25]
Telescopes - NANOGrav
The Arecibo Telescope was one of the primary instruments in pulsar searching, leading to the discovery of many pulsars, including highly stable millisecond ...
[26]
Gravitational Wave Search Provides Insights into Galaxy Evolution ...
Apr 5, 2016 · NANOGrav is currently monitoring 54 pulsars, using the National Science Foundation's Green Bank Telescope in West Virginia and Arecibo Radio ...
[27]
North American Nanohertz Observatory for Gravitational Waves
Aug 2, 2019 · NANOGrav observes a pulsar timing array (PTA) of dozens of millisecond pulsars with approximately monthly cadence using the Robert C. Byrd Green ...
[28]
The NANOGrav 15 yr Data Set: Observations and Timing of 68 ...
We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15 yr data set of the North American Nanohertz Observatory for ...Missing: curve SMBHB J1906- 0746
[29]
Glossary - NANOGrav
Ultra-wideband (UWB) Receiver. This next-generation frontend will be able to receive 0.7–4.0 GHz simultaneously. Term. VEGAS. Versatile GBT Astronomical ...
[30]
Wideband Timing of 47 Millisecond Pulsars - IOPscience
Dec 21, 2020 · Ultrawideband receivers and simultaneous multiband observations are becoming norms in pulsar timing, and along with the anticipated increase ...Missing: ultra- | Show results with:ultra-
[31]
Astrophysical Interpretation of a Gravitational-Wave Background ...
In this study, we produce simulations of SMBH binary populations that contain billions of sources, and compare their predicted GW signatures with NANOGrav's ...<|control11|><|separator|>
[32]
Arecibo Pulsar Survey using ALFA. 2. The young, highly relativistic ...
We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected ...Missing: NANOGrav 2023 SMBHB
[33]
NANOGrav's 15-Year Data Release
Jun 28, 2023 · NANOGrav's 15-Year Data Release. Pub: 28 Jun 2023. Image. Artist Rendering of a Pulsar Timing Array with a Gravitational Wave Background.
[34]
NSF Funds NANOGrav Physics Frontiers Center
Jun 21, 2021 · NANOGrav is also a member of the International Pulsar Timing Array (IPTA) collaboration ... “West Virginia University Advances Collaborative ...
[35]
NSF funds NANOGrav Physics Frontiers Center
Jun 25, 2021 · NSF renewed funding for the North American Nanohertz Observatory for Gravitational Waves (NSF NANOGrav) with a $17 million grant over five ...Missing: IPTA | Show results with:IPTA
[36]
Partner Institutions - NANOGrav
IPTA members agree to share data and strategies and to work together collaboratively with the goal of detecting gravitational waves more quickly. Through the ...
[37]
International Pulsar Timing Array
### Summary of IPTA from http://ipta4gw.org/
[38]
International Pulsar Timing Array: second data release
ABSTRACT. In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three.Missing: outputs | Show results with:outputs
[39]
https://academic.oup.com/mnras/article/490/4/4666/5586597
[40]
https://ui.adsabs.harvard.edu/abs/2016MNRAS.458.1267V/abstract
[41]
The second data release from the European Pulsar Timing Array
The EPTA uses data from six European radio telescopes: the Effelsberg 100 m radio telescope (EFF) in Germany, the 76 m Lovell Telescope at Jodrell Bank ...<|control11|><|separator|>
[42]
The second data release from the European Pulsar Timing Array I ...
Jun 28, 2023 · The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe.
[43]
The second data release from the European Pulsar Timing Array
Such a common-spectrum process has already been observed in pulsar timing data. We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset ...
[44]
LEAP: the Large European Array for Pulsars - Oxford Academic
The five telescopes presently included in LEAP are the Effelsberg Telescope, the Lovell Telescope at Jodrell Bank, the Nançay Radio Telescope, the Sardinia ...
[45]
Practical approaches to analyzing PTA data: Cosmic strings with six ...
Jun 21, 2023 · We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars.
[46]
High-precision timing of 42 millisecond pulsars with the European ...
The EPTA also operates the Large European Array for Pulsars (LEAP), where data from the EPTA telescopes are coherently combined to form a tied-array telescope ...
[47]
The Parkes Pulsar Timing Array Project
Jan 24, 2013 · The EPTA, NANOGrav, and the PPTA are collaborating to form the International Pulsar Timing Array (IPTA) (Hobbs et al. Reference Hobbs2010b).<|separator|>
[48]
The Parkes Pulsar Timing Array project: second data release
Jun 5, 2020 · We describe 14 yr of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of ...Missing: CPTA | Show results with:CPTA
[49]
[1307.2629] The Parkes Pulsar Timing Array - arXiv
Jul 9, 2013 · The Parkes Pulsar Timing Array (PPTA) aims to detect gravitational waves, improve planetary ephemeris, and develop a pulsar-based time scale.Missing: future SKA
[50]
[PDF] Forecasting the sensitivity of Pulsar Timing Arrays to gravitational ...
Apr 3, 2024 · The SKA will provide high-precision pulsar timing measurement with uncer- tainties below ∼100ns [100], making it roughly 10 times better ...
[51]
The Chinese Pulsar Timing Array Data Release I
The Chinese Pulsar Timing Array (CPTA) has collected observations from 57 millisecond pulsars using the Five-hundred-meter Aperture Spherical Radio Telescope ( ...
[52]
The NANOGrav 15 yr Data Set: Observations and Timing of 68 ...
Jun 29, 2023 · We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15 yr data set of the North American Nanohertz ...
[53]
The Parkes Pulsar Timing Array third data release
Jul 19, 2023 · The EPTA, InPTA, MPTA, NANOGrav, and PPTA data sets will be combined to form the third IPTA data release, which will be the most sensitive ...
[54]
PINT: A Modern Software Package for Pulsar Timing - IOPscience
Apr 14, 2021 · In this paper, we describe the design, use, and validation of PINT, and we compare timing results between it and Tempo and Tempo2. Export ...
[55]
Comparing Recent Pulsar Timing Array Results on the Nanohertz ...
These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which ...
[56]
Searching for continuous Gravitational Waves in the second data ...
Mar 19, 2023 · We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. Subjects ...
[57]
The second data release from the European Pulsar Timing Array
A cosmic string origin would allow narrowing down the string tension to values of −11 ≲ log10Gμ ≲ −9.5, depending on the specific distribution of loops ...
[58]
[2406.16031] Source anisotropies and pulsar timing arrays - arXiv
Pulsar timing arrays (PTA) hunt for gravitational waves (GW) by searching for the correlations that GWs induce in the time-of-arrival residuals ...Missing: 2025 refined spectra
[59]
The NANOGrav 15 yr Data Set: Constraints on Supermassive Black ...
The NANOGrav 15 yr data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source ...
[60]
Implications of the pulsar timing array detections for massive black ...
Jul 23, 2023 · PTA data suggests massive black hole binaries merge efficiently, with masses larger than expected. LISA may detect at least a dozen, up to ...
[61]
Pulsar Timing Array Detections of Supermassive Binary Black Holes
Sep 26, 2023 · Pulsar timing arrays (PTAs) are anticipated to detect the stochastic gravitational-wave background (GWB) from supermassive binary black holes (BBHs)
[62]
Achieving precision cosmology with gravitational-wave bright sirens ...
Jan 3, 2022 · The GW signals from individual SMBHBs can be employed as standard sirens to measure absolute cosmic distances, and further provide constraints ...
[63]
Implications of pulsar timing array observations for LISA detections ...
May 10, 2023 · Pulsar timing arrays (PTAs) and the Laser Interferometer Space Antenna (LISA) will open complementary observational windows on massive black-hole binaries ( ...
[64]
Novel tests of gravity using nano-Hertz stochastic gravitational-wave background signals
### Summary: Testing General Relativity (GR) in the Strong-Field Limit Using PTA Data for Supermassive Black Hole Binaries (SMBHBs)
[65]
International Pulsar Timing Array
The goal of the IPTA is to detect and characterize the low-frequency gravitational wave universe through timing a global array of approximately 100 millisecond ...
[66]
Extending gravitational wave burst searches with pulsar timing arrays
More pulsars in the array can therefore increase the sensitivity for certain sky areas. It should be noted that the sky coverage may be biased by selection ...
[67]
[PDF] Pulsar Timing Arrays - Inspire HEP
Jul 18, 2018 · An array of millisecond pulsars distributed across the sky can be used to detect low-frequency gravitational waves by searching for a ...
[68]
Regularizing the pulsar timing array likelihood: A path toward ...
Pulsar timing arrays (PTAs) are sensitive to gravitational wave backgrounds (GWBs) in the nano-Hertz frequency band ( 10 - 9 – 10 - 7 Hz ), a regime where the ...<|control11|><|separator|>
[69]
Machine learning for nanohertz gravitational wave detection and ...
Oct 20, 2020 · Machine learning for nanohertz gravitational wave detection and parameter estimation with pulsar timing array ... noise can be simulated ...
[70]
[PDF] Pulsar Timing Arrays and the SKA
Oct 5, 2021 · Pulsar surveys with. SKA1 will reveal a large new population of millisecond pulsars (MSPs): we expect a total of 700-900 new MSPs with SKA-MID ...