-
Implications for Galactic Electron Density Structure from Pulsar Sightlines Intersecting HII Regions
Authors:
S. K. Ocker,
L. D. Anderson,
J. Lazio,
J. M. Cordes,
V. Ravi
Abstract:
Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures; DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against HII region catalogs, we sh…
▽ More
Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures; DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against HII region catalogs, we show that the majority of pulsars with $\rm DM > 600$ pc cm$^{-3}$ and scattering delays $τ(1\ {\rm GHz}) > 10$ ms lie behind HII regions, and that HII region intersections may be relevant to as much as a third of the observed pulsar population. Accounting for HII regions resolves apparent discrepancies where Galactic electron density models place high-DM pulsars beyond the Galactic disk. By comparing emission measures (EMs) inferred from recombination line observations to pulsar DMs, we show that HII regions can contribute tens to hundreds of pc cm$^{-3}$ in electron column density along a pulsar LOS. We find that nearly all pulsars with significant excess (and deficit) scattering from the mean $τ$-DM relation are spatially coincident with known discrete ionized gas structures, including HII regions. Accounting for HII regions is critical to the interpretation of radio dispersion and scattering measurements as electron density tracers, both in the Milky Way and in other galaxies.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
A Heavily Scattered Fast Radio Burst Is Viewed Through Multiple Galaxy Halos
Authors:
Jakob T. Faber,
Vikram Ravi,
Stella Koch Ocker,
Myles B. Sherman,
Kritti Sharma,
Liam Connor,
Casey Law,
Nikita Kosogorov,
Gregg Hallinan,
Charlie Harnach,
Greg Hellbourg,
Rick Hobbs,
David Hodge,
Mark Hodges,
James W. Lamb,
Paul Rasmussen,
Jean J. Somalwar,
Sander Weinreb,
David P. Woody
Abstract:
We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the F…
▽ More
We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the FRB with a Milky Way-like host galaxy at $z_{\mathrm{host}} = 0.554$ of stellar mass $\mathrm{log}_{10}(M_{\star, \mathrm{host}}) = 10.20^{+0.04}_{-0.03} ~M_\odot$. We identify two intervening galaxy halos at redshifts $z_{\mathrm{igh1}} = 0.492$ and $z_{\mathrm{igh2}} = 0.438$, with low impact parameters, $b_{\mathrm{igh1}} = 43.0_{-11.3}^{+11.3}$ kpc and $b_{\mathrm{igh2}} = 36.1_{-11.3}^{+11.3}$ kpc, and intermediate stellar masses, $\mathrm{log}_{10}(M_{\star, \mathrm{igh1}}) = 10.01^{+0.02}_{-0.02} ~M_\odot$ and $\mathrm{log}_{10}(M_{\star, \mathrm{igh2}}) = 10.60^{+0.02}_{-0.02} ~M_\odot$. The presence of two such galaxies suggests that the sightline is significantly overcrowded compared to the median sightline to this redshift, as inferred from the halo mass function. We perform a detailed analysis of the sightline toward FRB 20221219A, constructing both DM and scattering budgets. Our results suggest that, unlike most well-localized sources, the host galaxy does not dominate the observed scattering. Instead, we posit that an intersection with a single partially ionized cloudlet in the circumgalactic medium of an intervening galaxy could account for the substantial scattering in FRB 20221219A and remain in agreement with typical electron densities inferred for extra-planar dense cloud-like structures in the Galactic and extragalactic halos (e.g., high-velocity clouds).
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
A Cyclic Spectroscopy Scintillation Study of PSR B1937+21 I. Demonstration of Improved Scintillometry
Authors:
Jacob E. Turner,
Timothy Dolch,
James M. Cordes,
Stella K. Ocker,
Daniel R. Stinebring,
Shami Chatterjee,
Maura A. McLaughlin,
Victoria E. Catlett,
Cody Jessup,
Nathaniel Jones,
Christopher Scheithauer
Abstract:
We use cyclic spectroscopy to perform high frequency-resolution analyses of multi-hour baseband Arecibo observations of the millisecond pulsar PSR B1937+21. This technique allows for the examination of scintillation features in far greater detail than is otherwise possible under most pulsar timing array observing setups. We measure scintillation bandwidths and timescales in each of eight subbands…
▽ More
We use cyclic spectroscopy to perform high frequency-resolution analyses of multi-hour baseband Arecibo observations of the millisecond pulsar PSR B1937+21. This technique allows for the examination of scintillation features in far greater detail than is otherwise possible under most pulsar timing array observing setups. We measure scintillation bandwidths and timescales in each of eight subbands across a 200 MHz observing band in each observation. Through these measurements we obtain intra-epoch estimates of the frequency scalings for scintillation bandwidth and timescale.Thanks to our high frequency resolution and the narrow scintles of this pulsar, we resolve scintillation arcs in the secondary spectra due to the increased Nyquist limit, which would not have been resolved at the same observing frequency with a traditional filterbank spectrum using NANOGrav's current time and frequency resolutions, and the frequency-dependent evolution of scintillation arc features within individual observations. We observe the dimming of prominent arc features at higher frequencies, possibly due to a combination of decreasing flux density and the frequency dependence of the plasma refractive index of the interstellar medium. We also find agreement with arc curvature frequency dependence predicted by Stinebring et al. (2001) in some epochs. Thanks to the frequency resolution improvement provided by cyclic spectroscopy, these results show strong promise for future such analyses with millisecond pulsars, particularly for pulsar timing arrays, where such techniques can allow for detailed studies of the interstellar medium in highly scattered pulsars without sacrificing the timing resolution that is crucial to their gravitational wave detection efforts.
△ Less
Submitted 21 June, 2024; v1 submitted 21 April, 2024;
originally announced April 2024.
-
The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background
Authors:
Gabriella Agazie,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Lucas Brown,
Sarah Burke-Spolaor,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Megan E. DeCesar,
Paul B. Demorest,
Heling Deng,
Timothy Dolch,
Elizabeth C. Ferrara,
William Fiore,
Emmanuel Fonseca,
Gabriel E. Freedman,
Nate Garver-Daniels
, et al. (58 additional authors not shown)
Abstract:
The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt…
▽ More
The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Searching for Magnetar Binaries Disrupted by Core-Collapse Supernovae
Authors:
Myles B. Sherman,
Vikram Ravi,
Kareem El-Badry,
Kritti Sharma,
Stella Koch Ocker,
Nikita Kosogorov,
Liam Connor,
Jakob T. Faber
Abstract:
Core-collapse Supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that $\sim45\%$ of magnetars should be nearby high-velocity stars. Here we conduct a multi-wavelength search for unbound stars, magnetar bi…
▽ More
Core-collapse Supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that $\sim45\%$ of magnetars should be nearby high-velocity stars. Here we conduct a multi-wavelength search for unbound stars, magnetar binaries, and SNR shells using public optical ($uvgrizy-$bands), infrared ($J-$, $H-$, $K-$, and $K_s-$bands), and radio ($888$ MHz, $1.4$ GHz, and $3$ GHz) catalogs. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalog we associate with SGR J1822.3-1606, and an X-ray pulsar we associate with 3XMM J185246.6+003317. Using a Markov-Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to $5\lesssim f_u \lesssim 24\%$, which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require $31\lesssim f_{nc} \lesssim 66\%$ of magnetars to descend from such channels. Our results support a high fraction ($48\lesssim f_m \lesssim 86\%$) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars.
△ Less
Submitted 15 May, 2024; v1 submitted 7 April, 2024;
originally announced April 2024.
-
NE2001p: A Native Python Implementation of the NE2001 Galactic Electron Density Model
Authors:
S. K. Ocker,
J. M. Cordes
Abstract:
The Galactic electron density model NE2001 describes the multicomponent ionized structure of the Milky Way interstellar medium. NE2001 forward models the dispersion and scattering of compact radio sources, including pulsars, fast radio bursts, AGNs, and masers, and the model is routinely used to predict the distances of radio sources lacking independent distance measures. Here we present the open-…
▽ More
The Galactic electron density model NE2001 describes the multicomponent ionized structure of the Milky Way interstellar medium. NE2001 forward models the dispersion and scattering of compact radio sources, including pulsars, fast radio bursts, AGNs, and masers, and the model is routinely used to predict the distances of radio sources lacking independent distance measures. Here we present the open-source package NE2001p, a fully Python implementation of NE2001. The model parameters are identical to NE2001 but the computational architecture is optimized for Python, yielding small (<1%) numerical differences between NE2001p and the Fortran code. NE2001p can be used on the command-line and through Python scripts available on PyPI. Future package releases will include modular extensions aimed at providing short-term improvements to model accuracy, including a modified thick disk scale height and additional clumps and voids. This implementation of NE2001 is a springboard to a next-generation Galactic electron density model now in development.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Jeremy Baier,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Dallas DeGan,
Paul B. Demorest
, et al. (74 additional authors not shown)
Abstract:
Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl…
▽ More
Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.
△ Less
Submitted 18 October, 2023;
originally announced October 2023.
-
The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B
Authors:
Gabriella Agazie,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Harsha Blumer,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Belinda D. Cheeseboro,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Megan E. DeCesar,
Paul B. Demorest,
Lankeswar Dey,
Timothy Dolch,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara
, et al. (63 additional authors not shown)
Abstract:
The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin…
▽ More
The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.
△ Less
Submitted 15 January, 2024; v1 submitted 29 September, 2023;
originally announced September 2023.
-
Pulsar Scintillation through Thick and Thin: Bow Shocks, Bubbles, and the Broader Interstellar Medium
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
D. R. Stinebring,
T. Dolch,
V. Pelgrims,
J. W. McKee,
C. Giannakopoulos,
D. J. Reardon
Abstract:
Observations of pulsar scintillation are among the few astrophysical probes of very small-scale ($\lesssim$ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially over-pressurized plasma in local ISM inhomogeneities, such as supernova remn…
▽ More
Observations of pulsar scintillation are among the few astrophysical probes of very small-scale ($\lesssim$ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially over-pressurized plasma in local ISM inhomogeneities, such as supernova remnants, HII regions, and bow shocks. Here we present a survey of eight pulsars conducted at the Five-hundred-meter Aperture Spherical Telescope (FAST), revealing a diverse range of scintillation arc characteristics at high sensitivity. These observations reveal more arcs than measured previously for our sample. At least nine arcs are observed toward B1929$+$10 at screen distances spanning $\sim 90\%$ of the pulsar's $361$ pc path-length to the observer. Four arcs are observed toward B0355$+$54, with one arc yielding a screen distance as close as $\sim10^5$ au ($<1$ pc) from either the pulsar or the observer. Several pulsars show highly truncated, low-curvature arcs that may be attributable to scattering near the pulsar. The scattering screen constraints are synthesized with continuum maps of the local ISM and other well-characterized pulsar scintillation arcs, yielding a three-dimensional view of the scattering media in context.
△ Less
Submitted 27 November, 2023; v1 submitted 24 September, 2023;
originally announced September 2023.
-
How to Detect an Astrophysical Nanohertz Gravitational-Wave Background
Authors:
Bence Bécsy,
Neil J. Cornish,
Patrick M. Meyers,
Luke Zoltan Kelley,
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Tyler Cohen,
James M. Cordes,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch
, et al. (71 additional authors not shown)
Abstract:
Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam…
▽ More
Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.
△ Less
Submitted 1 December, 2023; v1 submitted 8 September, 2023;
originally announced September 2023.
-
Comparing recent PTA results on the nanohertz stochastic gravitational wave background
Authors:
The International Pulsar Timing Array Collaboration,
G. Agazie,
J. Antoniadis,
A. Anumarlapudi,
A. M. Archibald,
P. Arumugam,
S. Arumugam,
Z. Arzoumanian,
J. Askew,
S. Babak,
M. Bagchi,
M. Bailes,
A. -S. Bak Nielsen,
P. T. Baker,
C. G. Bassa,
A. Bathula,
B. Bécsy,
A. Berthereau,
N. D. R. Bhat,
L. Blecha,
M. Bonetti,
E. Bortolas,
A. Brazier,
P. R. Brook,
M. Burgay
, et al. (220 additional authors not shown)
Abstract:
The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA…
▽ More
The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. 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 will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory
Authors:
Gabriella Agazie,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Harsha Blumer,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Megan E. DeCesar,
Dallas DeGan,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler,
Justin A. Ellis
, et al. (65 additional authors not shown)
Abstract:
We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign…
▽ More
We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline
Authors:
Aaron D. Johnson,
Patrick M. Meyers,
Paul T. Baker,
Neil J. Cornish,
Jeffrey S. Hazboun,
Tyson B. Littenberg,
Joseph D. Romano,
Stephen R. Taylor,
Michele Vallisneri,
Sarah J. Vigeland,
Ken D. Olum,
Xavier Siemens,
Justin A. Ellis,
Rutger van Haasteren,
Sophie Hourihane,
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Bence Bécsy,
J. Andrew Casey-Clyde
, et al. (71 additional authors not shown)
Abstract:
This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P…
▽ More
This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data.
△ Less
Submitted 7 July, 2023; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan DeCesar,
Paul B. Demorest,
Matthew C. Digman,
Timothy Dolch,
Brendan Drachler
, et al. (74 additional authors not shown)
Abstract:
Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav…
▽ More
Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler,
Elizabeth C. Ferrara,
William Fiore
, et al. (68 additional authors not shown)
Abstract:
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these…
▽ More
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Alexander Bonilla,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Curt J. Cutler
, et al. (89 additional authors not shown)
Abstract:
The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul…
▽ More
The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.
△ Less
Submitted 18 July, 2023; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Search for Signals from New Physics
Authors:
Adeela Afzal,
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Jose Juan Blanco-Pillado,
Laura Blecha,
Kimberly K. Boddy,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie
, et al. (98 additional authors not shown)
Abstract:
The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string…
▽ More
The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. Decesar,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler,
Elizabeth C. Ferrara,
William Fiore,
Emmanuel Fonseca
, et al. (66 additional authors not shown)
Abstract:
Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated…
▽ More
Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars
Authors:
Gabriella Agazie,
Md Faisal Alam,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Laura Blecha,
Victoria Bonidie,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Bence Bécsy,
Christopher Chapman,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler
, et al. (75 additional authors not shown)
Abstract:
We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival…
▽ More
We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Becsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Curt J. Cutler,
Megan E. DeCesar
, et al. (89 additional authors not shown)
Abstract:
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr…
▽ More
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
Authors:
Zaven Arzoumanian,
Paul T. Baker,
Laura Blecha,
Harsha Blumer,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Bence Bécsy,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Siyuan Chen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler,
Justin A. Ellis,
E. C. Ferrara,
William Fiore,
Emmanuel Fonseca,
Gabriel E. Freedman
, et al. (53 additional authors not shown)
Abstract:
Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s…
▽ More
Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B.
△ Less
Submitted 6 June, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
Scattering variability detected from the circumsource medium of FRB 20190520B
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
D. Li,
C. H. Niu,
J. W. McKee,
C. J. Law,
R. Anna-Thomas
Abstract:
Fast radio bursts (FRBs) are millisecond-timescale radio transients, the origins of which are predominantly extragalactic and likely involve highly magnetized compact objects. FRBs undergo multipath propagation, or scattering, from electron density fluctuations on sub-parsec scales in ionized gas along the line-of-sight. Scattering observations have located plasma structures within FRB host galaxi…
▽ More
Fast radio bursts (FRBs) are millisecond-timescale radio transients, the origins of which are predominantly extragalactic and likely involve highly magnetized compact objects. FRBs undergo multipath propagation, or scattering, from electron density fluctuations on sub-parsec scales in ionized gas along the line-of-sight. Scattering observations have located plasma structures within FRB host galaxies, probed Galactic and extragalactic turbulence, and constrained FRB redshifts. Scattering also inhibits FRB detection and biases the observed FRB population. We report the detection of scattering times from the repeating FRB 20190520B that vary by up to a factor of two or more on minutes to days-long timescales. In one notable case, the scattering time varied from $7.9\pm0.4$ ms to less than 3.1 ms ($95\%$ confidence) over 2.9 minutes at 1.45 GHz. The scattering times appear to be uncorrelated between bursts or with dispersion and rotation measure variations. Scattering variations are attributable to dynamic, inhomogeneous plasma in the circumsource medium, and analogous variations have been observed from the Crab pulsar. Under such circumstances, the frequency dependence of scattering can deviate from the typical power-law used to measure scattering. Similar variations may therefore be detectable from other FRBs, even those with inconspicuous scattering, providing a unique probe of small-scale processes within FRB environments.
△ Less
Submitted 29 November, 2022; v1 submitted 4 October, 2022;
originally announced October 2022.
-
Heliosphere Meets Interstellar Medium, in a Galactic Context
Authors:
Stella Koch Ocker,
James Cordes,
Shami Chatterjee,
Jeffrey Hazboun,
Timothy Dolch,
Daniel Stinebring,
Dustin Madison,
Stephen White,
Gregory Taylor,
Natalia Lewandowska,
Michael Lam
Abstract:
The physical conditions within our heliosphere are driven by the Sun's motion through an evolving interstellar environment that remains largely unexplored. The next generation of outer heliosphere and interstellar explorers will answer fundamental questions about the heliosphere's relationship with the very local interstellar medium (VLISM) by diving deeper into the Sun's interstellar surroundings…
▽ More
The physical conditions within our heliosphere are driven by the Sun's motion through an evolving interstellar environment that remains largely unexplored. The next generation of outer heliosphere and interstellar explorers will answer fundamental questions about the heliosphere's relationship with the very local interstellar medium (VLISM) by diving deeper into the Sun's interstellar surroundings. The impact of these future missions will be vastly enhanced by concurrent, interdisciplinary studies that examine the direct connections between conditions within the heliosphere, the heliosphere's immediate interstellar environment, and the larger-scale Galactic ISM. Comparisons of the heliosphere and VLISM to their analogs across the Galaxy will constrain the global processes shaping both stellar astrospheres and their sustained impact on the ISM.
△ Less
Submitted 24 August, 2022;
originally announced August 2022.
-
A Scintillation Arc Survey of 22 Pulsars with Low to Moderate Dispersion Measures
Authors:
Dan R. Stinebring,
Barney J. Rickett,
Anthony H. Minter,
Alex S. Hill,
Adam P. Jussila,
Lele Mathis,
Maura A. McLaughlin,
Stella Koch Ocker,
Scott M. Ransom
Abstract:
Context: By providing information about the location of scattering material along the line of sight (LoS) to pulsars, scintillation arcs are a powerful tool for exploring the distribution of ionized material in the interstellar medium. Here, we present observations that probe the ionized ISM on scales of $\sim$~0.001 -- 30~au.
Aims: We have surveyed pulsars for scintillation arcs in a relatively…
▽ More
Context: By providing information about the location of scattering material along the line of sight (LoS) to pulsars, scintillation arcs are a powerful tool for exploring the distribution of ionized material in the interstellar medium. Here, we present observations that probe the ionized ISM on scales of $\sim$~0.001 -- 30~au.
Aims: We have surveyed pulsars for scintillation arcs in a relatively unbiased sample with DM < 100 pc cm-3. We present multi-frequency observations of 22 low to moderate DM pulsars. Many of the 54 observations were also observed at another frequency within a few days.
Methods: For all observations we present dynamic spectra, autocorrelation functions, and secondary spectra. We analyze these data products to obtain scintillation bandwidths, pulse broadening times, and arc curvatures.
Results: We detect definite or probable scintillation arcs in 19 of the 22 pulsars and 34 of the 54 observations, showing that scintillation arcs are a prevalent phenomenon. The arcs are better defined in low DM pulsars. We show that well-defined arcs do not directly imply anisotropy of scattering. Only the presence of reverse arclets and a deep valley along the delay axis, which occurs in about 20\% of the pulsars in the sample, indicates substantial anisotropy of scattering.
Conclusions: The survey demonstrates substantial patchiness of the ionized ISM on both au size scales transverse to the line of sight and on $\sim$~100~pc scales along it. We see little evidence for distributed scattering along most lines of sight in the survey.
△ Less
Submitted 18 July, 2022;
originally announced July 2022.
-
Radio Scattering Horizons for Galactic and Extragalactic Transients
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
M. R. Gorsuch
Abstract:
Radio wave scattering can cause severe reductions in detection sensitivity for surveys of Galactic and extragalactic fast ($\sim$ms duration) transients. While Galactic sources like pulsars undergo scattering in the Milky Way interstellar medium (ISM), extragalactic fast radio bursts (FRBs) can also experience scattering in their host galaxies and other galaxies intervening their lines-of-sight. W…
▽ More
Radio wave scattering can cause severe reductions in detection sensitivity for surveys of Galactic and extragalactic fast ($\sim$ms duration) transients. While Galactic sources like pulsars undergo scattering in the Milky Way interstellar medium (ISM), extragalactic fast radio bursts (FRBs) can also experience scattering in their host galaxies and other galaxies intervening their lines-of-sight. We assess Galactic and extragalactic scattering horizons for fast radio transients using a combination of NE2001 to model the dispersion measure (DM) and scattering time ($τ$) contributed by the Galactic disk, and independently constructed electron density models for the Galactic halo and other galaxies' ISMs and halos that account for different galaxy morphologies, masses, densities, and strengths of turbulence. For source redshifts $0.5\leq z_{\rm s}\leq1$, an all-sky, isotropic FRB population has simulated values of $τ(1\rm~GHz)$ ranging from $\sim1~μ$s to $\sim2$ ms ($90\%$ confidence, observer frame) that are dominated by host galaxies, although $τ$ can be $\gg2$ ms at low Galactic latitudes. A population at $z_{\rm s}=5$ has $0.01\lesssimτ\lesssim300$ ms at 1 GHz ($90\%$ confidence), dominated by intervening galaxies. About $20\%$ of these high-redshift FRBs are predicted to have $τ>5$ ms at 1 GHz (observer frame), and $\gtrsim40\%$ of FRBs between $z_{\rm s}\sim0.5-5$ have $τ\gtrsim1$ ms for $ν\leq 800$ MHz. Our scattering predictions may be conservative if scattering from circumsource environments is significant, which is possible under specific conditions. The percentage of FRBs selected against from scattering could also be substantially larger than we predict if circumgalactic turbulence causes more small-scale ($\ll1$ au) density fluctuations than observed from nearby halos.
△ Less
Submitted 8 June, 2022; v1 submitted 30 March, 2022;
originally announced March 2022.
-
The Large Dispersion and Scattering of FRB 20190520B are Dominated by the Host Galaxy
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
C. -H. Niu,
D. Li,
J. W. McKee,
C. J. Law,
C. -W. Tsai,
R. Anna-Thomas,
J. -M. Yao,
M. Cruces
Abstract:
The repeating FRB 20190520B is localized to a galaxy at $z=0.241$, much closer than expected given its dispersion measure $\rm DM=1205\pm4\ pc\ cm^{-3}$. Here we assess implications of the large DM and scattering observed from FRB 20190520B for the host galaxy's plasma properties. A sample of 75 bursts detected with the Five-hundred-meter Aperture Spherical radio Telescope shows scattering on two…
▽ More
The repeating FRB 20190520B is localized to a galaxy at $z=0.241$, much closer than expected given its dispersion measure $\rm DM=1205\pm4\ pc\ cm^{-3}$. Here we assess implications of the large DM and scattering observed from FRB 20190520B for the host galaxy's plasma properties. A sample of 75 bursts detected with the Five-hundred-meter Aperture Spherical radio Telescope shows scattering on two scales: a mean temporal delay $τ(1.41\ {\rm GHz})=10.9\pm1.5$ ms, which is attributed to the host galaxy, and a mean scintillation bandwidth $ν_{\rm d}(1.41\ {\rm GHz})=0.21\pm0.01$ MHz, which is attributed to the Milky Way. Balmer line measurements for the host imply an H$α$ emission measure (galaxy frame) $\rm EM_s=620$ pc cm$^{-6} \times (T/10^4 {\rm K})^{0.9}$, implying $\rm DM_{\rm Hα}$ of order the value inferred from the FRB DM budget, $\rm DM_h=1121^{+89}_{-138}$ pc cm$^{-3}$ for plasma temperatures greater than the typical value $10^4$ K. Combining $τ$ and $\rm DM_h$ yields a nominal constraint on the scattering amplification from the host galaxy $\tilde{F} G=1.5^{+0.8}_{-0.3}$ (pc$^2$ km)$^{-1/3}$, where $\tilde{F}$ describes turbulent density fluctuations and $G$ represents the geometric leverage to scattering that depends on the location of the scattering material. For a two-screen scattering geometry where $τ$ arises from the host galaxy and $Δν_{\rm d}$ from the Milky Way, the implied distance between the FRB source and dominant scattering material is $\lesssim100$ pc. The host galaxy scattering and DM contributions support a novel technique for estimating FRB redshifts using the $τ-\rm DM$ relation, and are consistent with previous findings that scattering of localized FRBs is largely dominated by plasma within host galaxies and the Milky Way.
△ Less
Submitted 9 May, 2022; v1 submitted 27 February, 2022;
originally announced February 2022.
-
Magnetic field reversal in the turbulent environment around a repeating fast radio burst
Authors:
Reshma Anna-Thomas,
Liam Connor,
Shi Dai,
Yi Feng,
Sarah Burke-Spolaor,
Paz Beniamini,
Yuan-Pei Yang,
Yongkun Zhang,
Kshitij Aggarwal,
Casey J. Law,
Di Li,
Chenhui Niu,
Shami Chatterjee,
Marilyn Cruces,
Ran Duan,
Miroslav D. Filipovi,
George Hobbs,
Ryan S. Lynch,
Chenchen Miao,
Jiarui Niu,
Stella K. Ocker,
Chao-Wei Tsai,
Pei Wang,
Mengyao Xue,
Jumei Yao
, et al. (5 additional authors not shown)
Abstract:
Fast radio bursts (FRBs) are brief, intense flashes of radio waves from unidentified extragalactic sources. Polarized FRBs originate in highly magnetized environments. We report observations of the repeating FRB 20190520B spanning seventeen months , which show its amount of Faraday rotation is highly variable and twice changes its sign. The FRB also depolarizes below radio frequencies around 1 to…
▽ More
Fast radio bursts (FRBs) are brief, intense flashes of radio waves from unidentified extragalactic sources. Polarized FRBs originate in highly magnetized environments. We report observations of the repeating FRB 20190520B spanning seventeen months , which show its amount of Faraday rotation is highly variable and twice changes its sign. The FRB also depolarizes below radio frequencies around 1 to 3 GHz. We interpret these properties as due to change in the parallel component of the integrated magnetic field along the line-of-sight, including reversals. This could result from propagation through a turbulent, magnetized screen of plasma located between $10^{-5}$ to 100 parsecs of the FRB source. This is consistent with the bursts passing through the stellar wind of a binary companion of the FRB source.
△ Less
Submitted 12 May, 2023; v1 submitted 22 February, 2022;
originally announced February 2022.
-
Redshift Estimation and Constraints on Intergalactic and Interstellar Media from Dispersion and Scattering of Fast Radio Bursts
Authors:
James M. Cordes,
Stella Koch Ocker,
Shami Chatterjee
Abstract:
A sample of 14 FRBs with measured redshifts and scattering times is used to assess contributions to dispersion and scattering from the intergalactic medium (IGM), galaxy halos, and the disks of host galaxies. The IGM and galaxy halos contribute significantly to dispersion measures but evidently not to scattering, which is then dominated by host galaxies. This enables usage of scattering times for…
▽ More
A sample of 14 FRBs with measured redshifts and scattering times is used to assess contributions to dispersion and scattering from the intergalactic medium (IGM), galaxy halos, and the disks of host galaxies. The IGM and galaxy halos contribute significantly to dispersion measures but evidently not to scattering, which is then dominated by host galaxies. This enables usage of scattering times for estimating DM contributions from host galaxies and also for a combined scattering-dispersion redshift estimator. Redshift estimation is calibrated using scattering of Galactic pulsars after taking into account different scattering geometries for Galactic and intergalactic lines of sight. The DM-only estimator has a bias $\sim 0.1$ and RMS error $\sim 0.15$ in the redshift estimate for an assumed ad-hoc value of 50 pc cm$^{-3}$ for the host galaxy's DM contribution. The combined redshift estimator shows less bias by a factor of four to ten and a 20 to 40% smaller RMS error. We find that values for the baryonic fraction of the ionized IGM $f_{\rm igm} \sim 0.85 \pm 0.05$ optimize redshift estimation using dispersion and scattering. Our study suggests that two of the 14 candidate galaxy associations (FRB 20190523 and FRB~20190611B) should be reconsidered.
△ Less
Submitted 20 April, 2022; v1 submitted 2 August, 2021;
originally announced August 2021.
-
An In Situ Study of Turbulence Near Stellar Bow Shocks
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
T. Dolch
Abstract:
Stellar bow shocks are observed in a variety of interstellar environments and are shaped by the conditions of gas in the interstellar medium (ISM). In situ measurements of turbulent density fluctuations near stellar bow shocks are only achievable with a few observational probes, including H$α$ emitting bow shocks and the Voyager Interstellar Mission (VIM). In this paper, we examine density variati…
▽ More
Stellar bow shocks are observed in a variety of interstellar environments and are shaped by the conditions of gas in the interstellar medium (ISM). In situ measurements of turbulent density fluctuations near stellar bow shocks are only achievable with a few observational probes, including H$α$ emitting bow shocks and the Voyager Interstellar Mission (VIM). In this paper, we examine density variations around the Guitar Nebula, an H$α$ bow shock associated with PSR B2224$+$65, in tandem with density variations probed by VIM near the boundary of the solar wind and ISM. High-resolution Hubble Space Telescope observations of the Guitar Nebula taken between 1994 and 2006 trace density variations over scales from 100s to 1000s of au, while VIM density measurements made with the Voyager 1 Plasma Wave System constrain variations from 1000s of meters to 10s of au. The power spectrum of density fluctuations constrains the amplitude of the turbulence wavenumber spectrum near the Guitar Nebula to ${\rm log}_{10}C_{\rm n}^2 = -0.8\pm0.2$ m$^{-20/3}$ and for the very local ISM probed by Voyager ${\rm log}_{10}C_{\rm n}^2 = -1.57\pm0.02$ m$^{-20/3}$. Spectral amplitudes obtained from multi-epoch observations of four other H$α$ bow shocks also show significant enhancements in $C_{\rm n}^2$ from values that are considered typical for the diffuse, warm ionized medium, suggesting that density fluctuations near these bow shocks may be amplified by shock interactions with the surrounding medium, or by selection effects that favor H$α$ emission from bow shocks embedded in denser media.
△ Less
Submitted 28 September, 2021; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Persistent Plasma Waves in Interstellar Space Detected by Voyager 1
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
D. A. Gurnett,
W. S. Kurth,
S. R. Spangler
Abstract:
In 2012, Voyager 1 became the first in situ probe of the very local interstellar medium. The Voyager 1 Plasma Wave System has given point estimates of the plasma density spanning about 30 astronomical units (au) of interstellar space, revealing a large-scale density gradient and compressive turbulence outside the heliopause. Previous studies of the plasma density relied exclusively on the detectio…
▽ More
In 2012, Voyager 1 became the first in situ probe of the very local interstellar medium. The Voyager 1 Plasma Wave System has given point estimates of the plasma density spanning about 30 astronomical units (au) of interstellar space, revealing a large-scale density gradient and compressive turbulence outside the heliopause. Previous studies of the plasma density relied exclusively on the detection of discrete plasma oscillation events that are triggered ahead of shocks propagating outwards from the Sun, which were used to infer the plasma frequency and hence density. We present the detection of a class of very weak, narrowband plasma wave emission in the Voyager 1 Plasma Wave System data that persists from 2017 onwards and enables a steadily sampled measurement of the interstellar plasma density over about 10 au with an average sampling time of 3 days, or 0.03 au. We find au-scale density fluctuations that trace turbulence in the interstellar medium between episodes of previously detected plasma oscillations. Possible mechanisms for the narrowband emission include thermally excited plasma oscillations and quasi-thermal noise, and could be clarified by new findings from Voyager or a future interstellar mission. The persistence of the emission suggests that Voyager 1 may be able to continue tracking the interstellar plasma density in the absence of shock-generated plasma oscillation events.
△ Less
Submitted 13 December, 2021; v1 submitted 9 May, 2021;
originally announced May 2021.
-
Constraining Galaxy Haloes from the Dispersion and Scattering of Fast Radio Bursts and Pulsars
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee
Abstract:
Fast radio bursts (FRBs) can be scattered by ionized gas in their local environments, host galaxies, intervening galaxies along their lines-of-sight, the intergalactic medium, and the Milky Way. The relative contributions of these different media depend on their geometric configuration and the internal properties of the gas. When these relative contributions are well understood, FRB scattering is…
▽ More
Fast radio bursts (FRBs) can be scattered by ionized gas in their local environments, host galaxies, intervening galaxies along their lines-of-sight, the intergalactic medium, and the Milky Way. The relative contributions of these different media depend on their geometric configuration and the internal properties of the gas. When these relative contributions are well understood, FRB scattering is a powerful probe of density fluctuations along the line-of-sight. The precise scattering measurements for FRB 121102 and FRB 180916 allow us to place an upper limit on the amount of scattering contributed by the Milky Way halo to these FRBs. The scattering time $τ\propto(\tilde{F} \times {\rm DM}^2) A_τ$, where ${\rm DM}$ is the dispersion measure, $\tilde{F}$ quantifies electron density variations with $\tilde{F}=0$ for a smooth medium, and the dimensionless constant $A_τ$ quantifies the difference between the mean scattering delay and the $1/e$ scattering time typically measured. A likelihood analysis of the observed scattering and halo DM constraints finds that $\tilde{F}$ is at least an order of magnitude smaller in the halo than in the Galactic disk. The maximum pulse broadening from the halo is $τ\lesssim12$ $μ$s at 1 GHz. We compare our analysis of the Milky Way halo with other galaxy haloes by placing limits on the scattering contributions from haloes intersecting the lines-of-sight to FRB 181112 and FRB 191108. Our results are consistent with haloes making negligible or very small contributions to the scattering times of these FRBs.
△ Less
Submitted 9 May, 2022; v1 submitted 12 January, 2021;
originally announced January 2021.
-
Electron Density Structure of the Local Galactic Disk
Authors:
Stella Koch Ocker,
James M. Cordes,
Shami Chatterjee
Abstract:
Pulsar dispersion measures (DMs) have been used to model the electron density of the interstellar medium (ISM) in the Galactic disk as a plane-parallel medium, despite significant scatter in the DM-distance distribution and strong evidence for inhomogeneities in the ISM. We use a sample of pulsars with independent distance measurements to evaluate a model of the local ISM in the thick disk of the…
▽ More
Pulsar dispersion measures (DMs) have been used to model the electron density of the interstellar medium (ISM) in the Galactic disk as a plane-parallel medium, despite significant scatter in the DM-distance distribution and strong evidence for inhomogeneities in the ISM. We use a sample of pulsars with independent distance measurements to evaluate a model of the local ISM in the thick disk of the Galaxy that incorporates turbulent fluctuations, clumps, and voids in the electron density. The latter two components are required because $\sim 1/3$ of the lines of sight are discrepant from a strictly plane parallel model. A likelihood analysis for smooth components of the model yields a scale height $z_0=1.57^{+0.15}_{-0.14}$ kpc and a mid-plane density $n_0=0.015 \pm 0.001$ cm$^{-3}$. The scatter in the DM-distance distribution is dominated by clumps and voids but receives significant contributions from a broad spectrum of density fluctuations, such as a Kolmogorov spectrum. The model is used to identify lines of sight with outlier values of DM. Three of these pulsars, J1614$-$2230, J1623$-$0908, and J1643$-$1224, lie behind known HII regions, and the electron density model is combined with H$α$ intensity data to constrain the filling factors and other substructure properties of the HII regions (Sh 2-7 and Sh 2-27). Several pulsars also exhibit enhanced DM fluctuations that are likely caused by their lines of sight intersecting the superbubble GSH 238+00+09.
△ Less
Submitted 2 June, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
The Frequency Dependence of Scintillation Arc Thickness in Pulsar B1133+16
Authors:
Dan R. Stinebring,
Barney J. Rickett,
Stella Koch Ocker
Abstract:
Scintillation arcs have become a powerful tool for exploring scattering in the ionized interstellar medium. There is accumulating evidence that the scattering from many pulsars is extremely anisotropic resulting in highly elongated, linear brightness functions. We present a three-frequency (327~MHz, 432~MHz, 1450~MHz) Arecibo study of scintillation arcs from one nearby, bright, high-velocity pulsa…
▽ More
Scintillation arcs have become a powerful tool for exploring scattering in the ionized interstellar medium. There is accumulating evidence that the scattering from many pulsars is extremely anisotropic resulting in highly elongated, linear brightness functions. We present a three-frequency (327~MHz, 432~MHz, 1450~MHz) Arecibo study of scintillation arcs from one nearby, bright, high-velocity pulsar, PSR~B1133+16. We show that a one-dimensional (1D), linear brightness function is in good agreement with the data at all three observing frequencies. We use two methods to explore the broadening of the 1D brightness function $B(θ)$ as a function of frequency: 1) crosscuts of the forward arc at constant delay and 2) a 1D modeling of $B(θ)$ using a comparison between model and observed secondary spectrum as a goodness-of-fit metric. Both methods show that the half-power width of $B(θ)$ deviates from the expected dependence $\propto ν^{-a}$, where $ν$ is the observing frequency . Our estimates of $a$ have moderately large uncertainties but imply $a \lesssim1.8$, and so are inconsistent with the expected $a = 2.0$ for plasma refraction or $a = 2.2$ for Kolmogorov turbulence. In addition the shape of $B(θ)$ cuts off more steeply than predicted for Kolmogorov turbulence. Ultimately, we conclude that the underlying physics of the broadening mechanism remains unexplained. Our results place the scattering screen at a distance that is broadly consistent with an origin at the boundary of the Local Bubble.
△ Less
Submitted 11 November, 2018;
originally announced November 2018.