High Energy Astrophysical Phenomena
Understanding Pulsar Magnetospheres with the SKAO
L. S. Oswald, A. Basu, M. Chakraborty et al. (+21 more)
The SKA
telescopes will bring unparalleled sensitivity across a broad
radio band, a wide field of view across the Southern sky, and the capacity for sub-arraying, all of which make them the ideal instruments for studying the
pulsar magnetosphere. This chapter describes the advances that have been made in
pulsar magnetosphere physics over the last decade, and details how these have been made possible through the advances of modern
radio telescopes, particularly SKA precursors and pathfinders. It explains how the SKA
telescopes would transform the field of
pulsar magnetosphere physics through a combination of large-scale monitoring surveys and in-depth follow-up observations of unique sources and new discoveries. Finally, it describes how the specific observing opportunities available with the AA* and AA4 configurations will achieve the advances necessary to solve the problem of
pulsar radio emission physics in the coming years.
Read the full abstract ▾
All authors L. S. Oswald, A. Basu, M. Chakraborty, B. C. Joshi, N. Lewandowska, K. Liu, M. E. Lower, A. Philippov, X. Song, P. Tarafdar, J. van Leeuwen, A. L. Watts, P. Weltevrede, G. Wright, J. Benáček, A. Beri, S. Cao, P. Esposito, F. Jankowski, J. C. Jiang, A. Karastergiou, K. J. Lee, N. Rea, D. Vohl
Earth and Planetary Astrophysics
Comparing Monte Carlo Models of Impact Alteration of Planetary Atmospheres
MIkayla R. Huffman, David A. Brain, Kelsi N. Singer et al. (+2 more)
One process that affects atmospheric surface pressure is impact bombardment. The
evolution of a
planet's atmosphere under impact bombardment is an open question. We use a Monte Carlo method to evolve a range (0.006 to 92.5 bar) of initial atmospheres at Mars, Earth, and Venus under bombardment of 5x10^6 impactors using seven individual models. Since these seven models are best suited for specific impactor size regimes, we also combine these models into a composite model and compare it to other existing composites. Alterations to the existing models are required to apply to broad initial conditions. If we use each component model for every impactor, starting from present-day atmospheric pressure, we find about two or three orders of magnitude spread in the final atmospheric pressure. Given these differences, we suggest that the use of any one model to determine atmospheric change due to impact bombardment is risky. Most models and starting parameters result in net growth between +0.01 and +100 bar. Our composite model shows that the atmospheres of Venus, Earth, and Mars tend to grow under bombardment, with Earth's atmosphere growing most quickly. For an early Martian (P_0=1 bar) and an early terrestrial (with an initial pressure of P_0=0.25 bar) atmosphere, both tend to grow under bombardment. The results suggested here, where the models are universally applied, suggest that impact bombardment was likely a significant source of volatiles in the early
Solar System. Additional work and careful consideration of how impact events affect the
evolution of planetary atmospheres is needed.
Read the full abstract ▾
All authors MIkayla R. Huffman, David A. Brain, Kelsi N. Singer, Julia C. Johnston, Caitlin A. Caldwell
Solar and Stellar Astrophysics
Coronal Magnetography using Spectropolarimetry with SKA Telescopes
Deepan Patra, Puja Majee, Soham Dey et al. (+6 more)
The
solar coronal
magnetic field drives nearly every aspect of
solar phenomena and activity – from
flares, coronal mass ejections, and
solar wind that governs space weather to the much weaker nanoflares. These
magnetic fields are routinely measured at the visible surface of the Sun, the photosphere. However, detailed and direct measurements of the
magnetic fields in the
solar atmosphere, particularly in the coronal layer, have remained rather limited. Mostly, these are estimated from vector
magnetic field measurements at photospheric heights through different extrapolation models. In the case of the corona, these extrapolations lack observational constraints from the corona, especially during periods of intense activity when
magnetic structures evolve rapidly. Measurements of coronal
magnetic fields from observations, therefore, remain one of the most crucial and unresolved challenges in
solar and space-weather research.
Radio observations of the Sun hold considerable potential in this regard. Observations of diverse emission mechanisms, ranging from
plasma emissions at lower frequencies to thermal Bremsstrahlung and gyro-resonance at higher frequencies, provide multiple avenues to probe the coronal
magnetic fields, unique at
radio wavelengths. SKAO, with its broad frequency coverage (0.05 to 15 GHz), will allow us to probe wide range of coronal layers through unprecedented high-fidelity polarimetric imaging at high temporal, spectral, and spatial resolutions. This chapter details how the coronal
magnetic field measurements can be achieved through spectro-polarimetric imaging of the Sun with the SKAO.
Read the full abstract ▾
All authors Deepan Patra, Puja Majee, Soham Dey, Devojyoti Kansabanik, Divya Oberoi, Anshu Kumari, Surajit Mondal, Ketaki Deshpande, Divya Paliwal
Solar and Stellar Astrophysics
Archival Diagnostics for Potential Background Contaminantsof Project Hephaistos Dyson Sphere Candidates
Tongtian Ren, Michael A. Garrett, Erik Zackrisson et al. (+4 more)
We report on the diagnostic investigation of nine Project Hephaistos Dyson Sphere candidate M-
dwarfs based on archival data. By comparing the
Gaia DR3 positions with epoch 2016.0, propagated to the AllWISE epoch, with the mid-infrared centroids measured from the AllWISE images, together with deep archival optical/near-infrared imaging, we identified significant background contamination in candidates B and C–candidate B coincident with a
radio counterpart showing a power-law-like spectrum with a
radio spectral index alpha = 0.63 +/- 0.11, while candidate C has a near-infrared companion at an offset of 3.75 arcsec. Candidate A provides suggestive evidence through a
radio counterpart with spectral index alpha = 0.40 +/- 0.35, while candidates E, F, H and J show marginal evidence. These systems exhibit either significant astrometric offsets or visible interlopers, indicating that the mid-infrared excess likely arises from line-of-sight contamination by hot,
dust-obscured
galaxies. However, candidates D and I still lack obvious signs of contamination. Dedicated observations are therefore essential to characterise these potential interlopers, eliminate false positives, and ensure that technosignature searches focus on the most robust Dyson Sphere candidates.
Read the full abstract ▾
All authors Tongtian Ren, Michael A. Garrett, Erik Zackrisson, Andreas J. Korn, Andrew P. V. Siemion, Jason T. Wright, Alexis Brandeker
Instrumentation and Methods for Astrophysics
Machine Learning and the SKA for Cosmic Dawn and the Epoch of Reionization
Anshuman Acharya, Michele Bianco, Daniela Breitman et al. (+18 more)
When operational, the SKA will generate unprecedented amounts of data and provide exquisite sensitivity for 21 cm
tomography of
Cosmic Dawn (CD) and the Epoch of Reionization (EoR). With this comes opportunities for new data-driven
algorithms that unlock new methods for instrument modelling, data analysis, theoretical simulation, and inference for understanding the high-
redshift universe. In this chapter, we provide an overview of some machine learning
algorithms that have been proposed for CD and EoR science with the SKA
Read the full abstract ▾
All authors Anshuman Acharya, Michele Bianco, Daniela Breitman, Huaxi Chen, Abhirup Datta, Kangning Diao, Sambit K. Giri, Caroline S. Heneka, Nicholas Kern, Adrian Liu, Yashrajsinh Mahida, Suman Majumdar, Samit Kumar Pal, Shulei Ni, Yannic Pietschke, Davide Piras, Abinash Kumar Shaw, Hayato Shimabukuro, Ce Sui, Anshuman Tripathi, Xiaosheng Zhao
Astrophysics of Galaxies
From Atomic Gas to Star Formation
Erik Rosolowsky, Eric Koch, Sambit Roychowdhury et al. (+9 more)
Since the 2015 SKA Chapters, the advent of panchromatic observations of nearby
galaxies, capable of resolving the interstellar medium (ISM) into individual
star-forming regions and molecular clouds, has revolutionized our understanding of the
star formation cycle in
galaxies. Despite these advances, our understanding of the gas cycle in
galaxies remains limited because we know little about how the cold atomic phase of the ISM, which dominates its mass budget, participates in the cycle. Specifically, the processes regulating the condensation of atomic gas into molecular gas, the balance between warm and cold atomic gas, and their dependence on local environment and
star formation activity have been explored only in the Milky Way and the Magellanic Clouds. The SKA presents a new opportunity to observe the atomic ISM within nearby
galaxies at comparable resolution to other facilities and is the only facility that can make wide-area surveys of atomic gas across a statistical sample of
galaxies at high (<100 pc) spatial resolution. Such observations will directly address the open questions about the
evolution of the
star-forming ISM. We present recent results from SKA-like observations of the Local Group made with SKA pathfinders and the Jansky VLA to illustrate the promise of replicating those observations for hundreds of nearby
galaxies, reaching out to the Virgo cluster.
Read the full abstract ▾
All authors Erik Rosolowsky, Eric Koch, Sambit Roychowdhury, Luca Cortese, Filippo M. Maccagni, Barbara Catinella, Nushkia Chamba, Timothy A. Davis, Cosima Eibensteiner, Eric Murphy, Mamta Pandey-Pommier, Amidou Sorgho
Instrumentation and Methods for Astrophysics
SGN: A python framework for stream-processing pipelines
Yun-Jing Huang, Olivia Godwin, Chad Hanna et al. (+11 more)
We present the Stream Graph Navigator (SGN), a lightweight Python framework for building streaming data applications. In SGN, stream-processing pipelines are built by connecting computational components into directed acyclic graphs that run within an event loop. The time-series extension of the SGN library, SGN-TS, introduces signal processing methods to handle time series data. Together, SGN and SGN-TS provide the foundation for SGNL, a matched-filtering gravitational-wave search pipeline, and are being adopted by multiple projects across the low-latency gravitational-wave data analysis infrastructure as an extensible and maintainable framework for future gravitational-wave observations.
Read the full abstract ▾
All authors Yun-Jing Huang, Olivia Godwin, Chad Hanna, James Kennington, Jameson Rollins, Max Melching, Nathanael E Sovitzky, Aaron Viets, Madeline Wade, Zach Yarbrough, Yu-Kuang Chu, William Wyatt Phillips, Surabhi Sachdev, Rhiannon Udall
Cosmology and Nongalactic Astrophysics
Hubble tension in k-essence: Evidence for robust tension alleviation
Isaac Opio, Didam Duniya, Bishop Mongwane et al. (+1 more)
The Hubble tension has come to stay as a major problem in modern cosmology as it continues to plague the standard cosmological model ($\Lambda$CDM). As one of the viable, self-consistent dark energy theories, k-essence involves nontrivial self-interactions that can modify the background expansion beyond recombination; thereby impacting the sound horizon to last scattering, and hence, the inferred value of the Hubble constant. We examine this tension in two physically motivated k-essence models, dilaton and tachyon, using
datasets from Planck and late-Universe probes including Pantheon+SH0ES,
cosmic chronometer (CC),
Supernova Cosmology Project Union compilation (Union3), Dark Energy Survey Year~5 (DESY5), and Dark Energy Spectroscopic Instrument (DESI) measurements. While $\Lambda$CDM exhibits inconsistent tension inferences, both k-essence models exhibit a substantial tension alleviation that is robust against the inclusion of the independent late-Universe cosmological
datasets, giving consistent tension reduction irrespective of whether the observations are
supernovae (Pantheon+SH0ES, Union3, DESY5) alone or in combination with
cosmic chronometers (CC) and
baryon acoustic oscillation measurements (DESI). The combined late-Universe
dataset leads to only $0.14\sigma$ and $0.69\sigma$ offsets from the Planck prediction in the dilaton and tachyon models, respectively, compared to $5.89\sigma$ tension in $\Lambda$CDM. Both models demonstrate that the inferred tension alleviation is a stable, intrinsic consequence of the underlying k-essence dynamics rather than of model fine tuning: model parameters remain unchanged across
datasets. The results establish that the apparent Hubble tension is not an unavoidable feature of late-Universe cosmology but depends critically on the description of dark energy.
Read the full abstract ▾
All authors Isaac Opio, Didam Duniya, Bishop Mongwane, Hassan Abdalla
Cosmology and Nongalactic Astrophysics
Inflation with nondynamic distortion to leading order in slow roll
Ali Hassan
We study inflation in metric-affine
gravity. We write an action that contains all the second order algebraic distortion terms, and all first order distortion terms with a single covariant derivative, coupled to a scalar field. We include the Einstein–
Hilbert term with nonminimal coupling, a scalar field potential, and impose projective invariance. The distortion equation of motion is algebraic by construction, and the distortion is integrated out analytically. This yields a kinetic term sourced entirely by distortion, with a kinetic coupling function determined by the 13 free coupling constants of the starting action. We compute inflationary observables for three model classes with a monomial distortion coupling. For a monomial potential, the spectral index and tensor-to-scalar ratio depend only on the ratio of the exponents, with the starting coupling constants dropping out entirely; however, this model lies outside the Planck + BK18 $2\sigma$ contours. For a potential of the $\alpha$-attractor form, the observables are governed by a single parameter and approach the Starobinsky predictions as a limit. Including a nonminimal coupling to the
Ricci scalar with a monomial potential can also yield an asymptotically flat effective potential with the same modified Starobinsky observables.
Read the full abstract ▾
Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc
Black Hole Memory Burden and its Signatures in Gravitational Waves from Mergers
Gia Dvali, Michael Zantedeschi and Sebastian Zell
Swift
memory burden (MB) implies that the information stored in a black hole (BH) can modify its classical dynamics when the BH is perturbed. This influences the gravitational waves (GWs) emitted during BH mergers. In this paper, we investigate how the BH
memory load is determined by the features of the collapsing source. We show that the
memory load can vastly exceed the information content of its progenitor. An extreme example is a BH formed in a two-particle collision, which exhibits maximal MB. We then derive bounds for BHs formed through
stellar collapse and examine the impact of swift MB on BH quasinormal modes, quantifying the MB-induced frequency shift of GWs. These findings imply that GW observations probe the fundamental mechanisms of BH information storage as well as their formation history.
Read the full abstract ▾
High Energy Astrophysical Phenomena
Lucas M. Pasquevich, Gustavo E. Romero and Matías M. Reynoso
Multimessenger observations suggest that Seyfert
galaxies are promising sources of high-energy
neutrinos, but their dense inner environments can strongly suppress the emerging very-high-energy gamma-ray emission. Active
galactic nuclei (AGN) undergoing intense
accretion episodes can enter a super-Eddington state, in which the
accretion flow becomes geometrically and optically thick within a critical radius and develops strong
magnetic fields in its innermost region. At the same time, large amounts of matter are expelled from the disk surface in the form of powerful, radiation-driven winds. In this work, we explore a scenario in which the cores of super-Eddington AGN provide suitable conditions for the acceleration of relativistic particles, including
hadrons, via
magnetic reconnection in a magnetically confined region close to the supermassive black hole. The accelerated hadronic component interacts with the intense
photon field of the disk, leading to a copious
neutrino flux peaking at 10-100 TeV that may be detectable with current observatories such as IceCube and KM3NeT, while the surrounding
outflow efficiently absorbs the accompanying gamma-ray emission from the inner core. We also apply the model to the nearby super-Eddington Seyfert 1 NGC 7469 as a representative case with two reported
neutrino events. In this framework, super-Eddington AGN, in particular Seyfert
galaxies undergoing transient intense
accretion episodes, emerge as plausible hidden
neutrino sources, offering a natural explanation for the coexistence of efficient
neutrino production and a strongly attenuated gamma-ray counterpart.
Read the full abstract ▾
High Energy Astrophysical Phenomena
Beyond the Impulsive Approximation: The Dynamics and Radio Emission of Tidal Disruption Jets
Fabio De Colle and Enrico Ramirez-Ruiz
Radio emission from relativistic
tidal disruption events (TDEs) is commonly interpreted using the Blandford–McKee solution for an impulsive relativistic blast wave. Unlike gamma-ray bursts, however, TDE jets are powered by fallback
accretion, injecting energy over weeks to months. We investigate how this sustained energy injection modifies the dynamics and
radio emission of relativistic TDE jets using one- and two-dimensional relativistic hydrodynamic simulations coupled with synchrotron radiation calculations. For a fixed total energy, fallback-powered jets drive slower forward shocks, producing delayed and systematically fainter
radio emission. We show that the impulsive approximation is valid only when the duration of energy injection is less than approximately ten percent of the deceleration time of an equivalent Blandford–McKee blast wave, a condition that is generally not satisfied for relativistic TDEs. As a result, continuously powered
outflows are significantly less beamed than impulsive explosions, leading to a strong suppression of off-axis
radio emission. These results demonstrate that the
radio evolution of relativistic TDEs retains
memory of the central engine and cannot, in general, be described by an impulsive solution. By substantially reducing the expected off-axis
radio emission, sustained fallback-powered jets imply that relativistic jet production in
tidal disruption events may be considerably more common than current
radio observations suggest.
Read the full abstract ▾
Astrophysics of Galaxies
A MIGHTEE robust measurement of the star formation rate-radio correlation
Charlotte L. Jackson, Matt J. Jarvis, Imogen H. Whittam et al. (+9 more)
Determining the relationship between
star-formation rate (SFR) and the
radio luminosity ($L_{1.4}$) is critical if we are to trace the
star-formation history of the Universe
dust-agnostically using current and future
radio facilities. However, until now, such work has relied on potentially biased binary classifications of sources to remove contaminating active
galactic nuclei (AGN). We present a new, statistically-driven methodology for deriving the SFR – $L_{1.4}$ relation, removing the need for problematic cuts. We use a
Bayesian hierarchical mixture model fit to the
radio-detected sources in the deep MIGHTEE COSMOS DR1 catalogue, incorporating the full SFR posterior probability distributions generated by state-of-the-art spectral energy distribution fitting code GRAHSP. This allows us to probabilistically determine a mean SFR – $L_{1.4}$ relation for the SF dominated
galaxies, whilst accounting for changing fractions of SF dominated sources across
redshift,
radio luminosity and
stellar mass ranges. We find that the SFR –
radio luminosity correlation exhibits a significant dependence on
redshift, but a
stellar mass dependence that is weaker than previous studies. Our resultant SFR-
radio correlation is $\log_{10}(\text{SFR}/M_{\odot}\,\text{yr}^{-1}) = 0.790\times(\log _{10}(L_{1.4}/\text{W\,Hz}^{-1})-23) + 1.244 \times(1+z)^{0.122} -0.033 \times (\log_{10}(M_*/M_{\odot})-10)$, with an intrinsic scatter of 0.178 dex. We show that this
redshift evolution could be explained by a moderate
evolution in the
radio spectral index of SF
galaxies. We attribute the lack of observed strong dependence on
stellar mass, compared to recent studies, to the novel statistical approach that does not rely on cuts to remove AGN.
Read the full abstract ▾
All authors Charlotte L. Jackson, Matt J. Jarvis, Imogen H. Whittam, James H. Matthews, Johannes Buchner, Catherine L. Hale, Joel Hamlett, S. Lyla Jung, Mara Salvato, Daniel J. B. Smith, Natalia Stylianou, Rohan G. Varadaraj
High Energy Astrophysical Phenomena
Multiwavelength properties of short gamma ray bursts with extended emission observed by Swift
M.M. Dinatolo, A. Mei, R. Brivio et al. (+11 more)
Short gamma-ray bursts with extended emission (SGRBEEs) are a particular class of long GRBs (LGRBs) which, despite their duration, share several observational features with short GRBs (SGRBs). They are composed by a short, hard initial pulse (IP) followed by a longer and softer extended emission (EE). We investigate whether SGRBEEs originate from the same progenitor as SGRBs despite their duration, representing a peculiar subclass of LGRBs, or if they constitute a distinct population. Given their respective duration, we tested if the IP and the EE share properties with short and long GRBs, respectively. We analysed SGRBEEs from the flux-limited,
redshift-complete SBAT4 sample using
prompt-emission data from Swift/BAT,
Fermi/GBM and Konus-WIND, and X-ray afterglow observations from Swift/XRT. The temporal and spectral properties of the IP and EE components were compared with those of SGRBs from the extended SBAT4 sample and LGRBs from the BAT6 sample. Despite observing a clear spectral
evolution during the
prompt phase of each SGRBEE, IPs and EEs as well as short and long GRBs can not be distinguished by their hardness ratio only. All bursts analysed have a spectral lag consistent with zero. SGRBEEs XRT light curves are consistently more complex than those of SGRBs, requiring the addition of multiple breaks and showing the presence of steep decays and plateaus. Our results indicate that SGRBEEs are not an intermediate class. During
prompt emission, despite common spectral features, IPs and EEs temporally differ from short and long GRBs, respectively. EEs are on average too faint to appear in the Amati plane, but IPs occupy the same parameter space of SGRBs. In the afterglow, SGRBEEs are more luminous than standard GRBs at early-time, suggesting a direct contribution from the EE; at later times, these events behave as standard SGRBs, and both remain systematically less luminous than LGRBs.
Read the full abstract ▾
All authors M.M. Dinatolo, A. Mei, R. Brivio, M. Ferro, M.G. Bernardini, P. D'Avanzo, S. Belova, S. Campana, S. Covino, D. Frederiks, B. Haskell, R. Salvaterra, B. Sbarufatti, A. Tsvetkova
Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph
Gravitational Wave Signatures of Cosmological Stasis: A Unified Spectral Template
Gabriela Barenboim and Anne-Katherine Burns
Proposed in 2022 by Dienes et al., stasis is a dynamical fixed point in the early universe in which the equation of state, $w_s$, is fixed at a constant value. In this study we show that the inflationary gravitational wave imprint of any stasis epoch is captured by a closed-form spectral template controlled by two physical inputs, the equation of state $w_s$ and the stasis duration $\Delta N_\mathrm{stasis}$, that applies uniformly across every microphysical realization. The template presented here yields two independently measurable observables, the spectral tilt of the spectra in the stasis band, $\alpha(w_s)$ and the amplitude step $C^2(w_s)$ at the beginning and end of the stasis band. Eliminating $w_s$ gives a one-parameter consistency curve $C^2 = C^2(\alpha)$ on which the data must lie if the underlying cosmology is any constant-$w$ era. This makes the spectrum falsifiable without knowing $w_s$ in advance: a measured $(\alpha, C^2)$ pair either lands on the curve or rules out the constant-$w$ class. We show that BBO and DECIGO can resolve the perpendicular displacement from the consistency curve to $\sigma_\perp \simeq 1.5\times10^{-5}$ at a tensor-to-scalar ratio, $r = 0.01$, four orders of magnitude below the curve's range in $C^2$ meaning that any off-curve deviation is detectable across a broad range of allowed $r$ values.
Read the full abstract ▾
Astrophysics of Galaxies
Pebbles to Gems: Intermediate-mass black holes in the first star clusters
Benedetta Mestichelli, Manuel Arca Sedda, Marta Volonteri et al. (+9 more)
The rapid assembly of supermassive black holes (SMBHs) observed at $z\gtrsim7$ requires efficient seeding mechanisms in the early Universe. Population III (Pop. III)
star clusters have recently emerged as a promising pathway that may bridge the gap between traditional light- and heavy-seed scenarios by producing intermediate-mass black holes (IMBHs) with masses up to $\sim10^4\,\rm M_{\odot}$. We investigate the properties and number densities of IMBHs forming in Pop. III
star clusters with masses $M_{\rm cl}\sim10^3-4\times10^5\,\rm M_{\odot}$, and hosted in isolated dark matter minihalos, using a suite of direct $N$-body simulations. We adopt cosmologically motivated initial conditions and explore different
stellar evolution prescriptions, binary orbital parameter distributions, and cluster dynamical configurations. By $z\sim19$, the IMBH mass function consistently peaks at $m_{\rm IMBH}\sim200\,\rm M_{\odot}$, with number densities of $n_{\rm IMBH}\sim0.2-5\,\rm cMpc^{-3}$. In sufficiently dense and massive clusters, IMBHs with masses $>10^3\,\rm M_{\odot}$ can already form by $z\sim19$, reaching number densities of $n_{\rm IMBH}\sim10^{-4}-10^{-2}\,\rm cMpc^{-3}$. The most massive IMBHs in our models reach $\sim6200\,\rm M_{\odot}$ through the collapse of very massive
stars assembled by repeated
stellar collisions, a process enhanced in fractal clusters. Lower-mass IMBHs form instead predominantly through single and binary
stellar evolution and binary
stellar mergers. We find that models combining large
stellar radii and tight binaries produce the highest IMBH abundances relative to isolated Pop. III
evolution. Owing to the high retention fraction of IMBHs ($\gtrsim88\%$), massive dense Pop. III
star clusters can act as efficient incubators of both light and heavy SMBH seeds, even if only a fraction of Pop. III
stars formed in such environments.
Read the full abstract ▾
All authors Benedetta Mestichelli, Manuel Arca Sedda, Marta Volonteri, Michela Mapelli, Stefano Torniamenti, Alessandro Lupi, Marica Branchesi, Shingo Hirano, Tomoaki Ishiyama, Ralf S. Klessen, Veronika Lipatova, Boyuan Liu
Cosmology and Nongalactic Astrophysics
Pulsar Timing Sensitivity to Dark Matter Substructure in the Presence of a Stochastic Gravitational-Wave Background
Abhiram Cherukupalli, Vincent S. H. Lee, Kim Berghaus et al. (+1 more)
Pulsar timing arrays (PTAs) can detect dark matter (DM) substructure through the small shifts a transiting DM subhalo imprints on pulse arrival times. Recently found evidence for a
stochastic gravitational-wave background (GWB) acts as red noise and competes with the substructure signal. Here we provide an analytic understanding of how this background degrades PTA sensitivity to DM substructure. We derive the full gauge-invariant proper-time observable induced by a transiting DM subhalo and develop a framework for the expected signal-to-noise ratio in the presence of red noise, accounting for the degeneracy with the
pulsar timing model. From this we obtain simple scaling relations for the reach across the static, dynamic, and
stochastic regimes, and numerically compute the reach for a Square Kilometre Array
benchmark at three representative points of the NANOGrav 15-year posterior. The GWB background suppresses the sensitivity to DM substructure by one to three orders of magnitude compared to forecasts in the presence of only white noise, and the suppression depends on the amplitude and spectral index of the background within a factor of three. The dynamic Shapiro signal suffers the smallest suppression and gives the best sensitivity to DM substructure near $10^{-2}\, M_\odot$. Probing the regime where subhalos make up all or part of the DM remains a challenge even for surveys with more
pulsars and longer observing time. Despite this, PTA measurements remain a competitive probe of DM substructure, and future surveys will increase in sensitivity by up to two orders of magnitude from existing NANOGrav limits.
Read the full abstract ▾
All authors Abhiram Cherukupalli, Vincent S. H. Lee, Kim Berghaus, Kathryn M. Zurek