Interpretation of Nanohertz Gravitational Waves via Cosmic Strings in Post-BBN Matter and Kination Domination

arXiv:2610.10674 · astro-ph.CO, hep-ph · Submitted 2026-10-07 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera, and with me are Jocelyn and Subrahmanyan, guest researcher.

Jocelyn: Today's paper: "Interpretation of Nanohertz Gravitational Waves via Cosmic Strings in Post-BBN Matter and Kination Domination".

Vera: The gist The work investigates whether a transient period of matter and kination domination in an early Universe can significantly modify the gravitational-wave spectrum from cosmic strings,

Jocelyn: First, who's behind it and why it matters.

Paper summary: Vera: So we’re looking at this paper now called "Interpretation of Nanohertz Gravitational Waves via Cosmic Strings in Post-BBN Matter and Kination Domination". Jocelyn, what's the main idea here?

Jocelyn: Basically, they’re taking the signals we see from pulsar timing arrays, which are these very low frequency gravitational waves, and trying to see if cosmic strings could be what’s making that noise. They focus on a specific type of early universe history between big bang nucleosynthesis and recombination.

Subrahmanyan: That means they are looking at how the matter density in the early universe changed after it was dominated by radiation, specifically exploring an era where matter domination happens before kination, which is when the background energy density redshifts faster than radiation.

Vera: Right, so what’s their central claim about this cosmic string interpretation? What are they trying to prove?

Jocelyn: The thesis of this paper is that a transient period of matter and kination domination in an early universe can significantly modify the gravitational-wave spectrum from cosmic strings, and this modification provides a potential explanation for the signals we observe by pulsar timing arrays.

Subrahmanyan: They argue that this non-standard cosmological history can generate characteristic features in the nanohertz range, which is where these PTA observations are making their signals.

Vera: So why does that matter for us, for the people listening right now? Why should we care about this specific way the universe expanded early on?

Jocelyn: Because the resulting spectral structure changes dramatically depending on whether you have matter domination or kination domination in that timeframe. They introduce characteristic temperatures, TRM and TKR, to describe these transitions.

Subrahmanyan: When they look at the results, they find that this modified expansion history deforms the broad cosmic string spectrum into a triangular shape.

Vera: A triangular shape sounds specific. What does it mean for those frequencies we are listening to right now?

Jocelyn: The segment from fKR to fMK probes loops associated with the kination era, where the background energy density redshifts faster than radiation, which leads to a spectral tilt where GW proportional to f+one <ref:2610.10674#pg2>.

Subrahmanyan: And then there's another segment from f MK to f RM that probes loops from the preceding matter era, which results in a negative spectral tilt towards higher frequencies, where GW proportional to f-one <ref:2610.10674#pg2>.

Vera: That’s quite a contrast between the two regimes. So how do they actually test this idea against real data?

Jocelyn: They use a Bayesian framework to compare competing models, starting with hypotheses H1 through H4 against a baseline hypothesis H0, which is just the standard smooth expansion. They look at the log Bayes factor, B ij, to see how well each model fits the NANOGrav fifteen-year and EPTA DR2 New+ datasets <ref:2610.10674#pg1,the NANOGrav 15-year and EPTA DR2 New+ datasets>.

Paper summary: Subrahmanyan: They specifically compare hypothesis H3 against hypothesis H1 in the NG15 dataset, and they found substantial support for kination over the standard lambda-CDM model, with a log Bayes factor of ln B NG153 = one point two nine plus or minus zero point zero six.

Vera: So you’re saying that for the NG15 data alone, there’s some evidence pointing toward kination being favored over the standard cosmology?

Jocelyn: That's what they found in that comparison, and they also looked at the EPTA DR2 New+ dataset. They found no strong evidence favoring any of these alternative hypotheses over the baseline H0 in either dataset individually.

Subrahmanyan: However, Jocelyn mentioned something interesting about how kination affects the spectrum generally; it actually enhances the SGWB spectrum toward the high-frequency end, resulting in a more pronounced blue tilt when compared to its counterpart in lambda-CDM.

Vera: So even though they didn't find strong evidence for kination over standard cosmology in either dataset alone, it still changes how we interpret the signal. What’s the catch there?

Jocelyn: The paper points out that while both datasets mildly favor kination over lambda-CDM within the cosmic string framework, there are overlapping pulsars between NG15 and EPTA DR2 New+.

Subrahmanyan: That overlap is a complication because combining those two likelihoods would double count information associated with the common pulsars, which could potentially bias the model constraints.

Vera: So what’s the practical implication for us listening? What does this mean for our understanding of cosmic string signals?

Jocelyn: The authors suggest that a statistically consistent joint analysis, one that properly treats those correlations between overlapping datasets, would be necessary to get the tightest constraints on model parameters and model evidences.

Subrahmanyan: For future work, they point out that the next step is for future studies to clarify the nature of this SGWB and what it implies about the early expansion history of our universe.

Vera: So in conclusion, what’s the big picture takeaway from this paper on "Interpretation of Nanohertz Gravitational Waves via Cosmic Strings in Post-BBN Matter and Kination Domination"?

Jocelyn: The authors conclude that while they found no strong evidence favoring any single alternative hypothesis over the baseline H0 in either dataset, the kination scenario considered enhances the SGWB spectrum toward higher frequencies, which gives it a better fit compared to its counterpart in lambda-CDM.

Subrahmanyan: They also note that within the cosmic string framework, with or without contributions from SMBHBs, both datasets mildly favor kination over lambda-CDM.

Vera: And they highlight that the inferred cosmic-string-tension posteriors remain largely consistent between the two cosmologies they analyzed.

Jocelyn: But again, because of those overlapping pulsars in the NG15 and EPTA DR2 New+ datasets, combining them needs careful statistical treatment to avoid biasing the evidence.

Vera: And as a final thought for our listeners, what’s on the horizon for this research?

Subrahmanyan: The next step involves future work to clarify what this stochastic gravitational wave background actually is and what it tells us about how the universe expanded in its earliest moments.

Conclusion: Vera: So we've been looking at this paper, "Interpretation of Nanohertz Gravitational Waves via Cosmic Strings in Post-BBN Matter and Kination Domination," and now we're wrapping up.

Jocelyn: It basically boils down to whether a specific early universe history—one with matter domination followed by kination—can actually explain the gravitational wave signals we see from pulsar timing arrays.

Subrahmanyan: The authors are showing how this transient period of matter and kination domination modifies the spectrum of cosmic strings in a way that might match what we observe across different frequency bands.

Vera: So, what’s the main conclusion here about these cosmic string signals? Are they actually pointing toward this early universe history?

Jocelyn: They find that while there isn't one single smoking gun, the kination scenario does create a spectral tilt that makes it look better than the standard lambda-CDM model in certain comparisons.

Subrahmanyan: The paper uses Bayesian analysis to compare different cosmological models against real data from pulsar timing arrays like NG15 and EPTA DR2 New+.

Vera: But what's the catch with that comparison, Jocelyn? Because they mention those overlapping pulsars.

Jocelyn: Yeah, they point out that combining those two datasets would double count information from the same pulsars if you don't use a specific statistical treatment for those correlations.

Subrahmanyan: That’s a real issue for getting tight constraints on the model parameters because of how those datasets overlap across different sky regions.

Vera: So, what does this mean for us listening right now? What's the big picture implication of finding kination might be favored in these models?

Jocelyn: It suggests that even if we can’t definitively rule out standard cosmology yet, the kination scenario is a plausible way to generate those specific nanohertz gravitational wave shapes we are hunting for.

Subrahmanyan: It opens up a new avenue for understanding the physics of the very early universe expansion, which is something we've been trying to constrain with these observations.

Vera: So, what's the next step for this research team? Where does this paper leave us heading?

Jocelyn: The authors suggest that future work needs to focus on getting a better understanding of exactly what this stochastic gravitational wave background looks like and what it really tells us about the expansion history.

Raymond T. Co, Siu Cheung Lam, Taegyu Lee

Indiana University

astro-ph.CO, hep-ph

Submitted: 2026-10-07

Updated: 2026-10-07

Comments: 28 pages, 10 figures, and 8 tables

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 92/100

The gist: The gist The work investigates whether a transient period of matter and kination domination in an early Universe can significantly modify the gravitational-wave spectrum from cosmic strings,

Key concepts

Cosmic String Gravitational Waves: Cosmic strings are hypothetical one-dimensional topological defects in the early Universe. They are expected to produce a specific stochastic gravitational wave background (SGWB) in the nanohertz frequency range, which is sensitive to how the Universe expanded during its formation.
Matter-Kination Domination Sequence: This describes a non-standard cosmological history between big bang nucleosynthesis and recombination where an early matter-dominated era transitions into a kination phase. Kination occurs when the background energy density redshifts faster than radiation, significantly altering the gravitational wave spectrum produced by cosmic strings.
Spectral Tilt: The spectral tilt refers to how the strength of gravitational waves changes across different frequencies. The paper shows that matter-kination domination deforms the broad cosmic string spectrum into a triangular shape, leading to specific frequency dependencies ($\Omega_{GW} \propto f+1$ or $\Omega_{GW} \propto f-1$) depending on whether the loop formation occurred during the matter or kination era.

Terminology

Abundance — 16× in this paper · explained in 1 episode
Bayesian Inference — 4× in this paper · explained in 5 episodes
European Pulsar Timing Array — 3× in this paper
Phase Transitions — 3× in this paper · explained in 2 episodes
pulsar timing arrays — 3× in this paper · explained in 2 episodes
Symmetry Breaking — 2× in this paper · explained in 4 episodes
Anisotropy — 1× in this paper · explained in 4 episodes
Bayesian Framework — 1× in this paper · explained in 3 episodes
Benchmark — 1× in this paper · explained in 2 episodes
Core — 1× in this paper · explained in 2 episodes
Cosmic Linear Anisotropy Solving System — 1× in this paper
Gravitational Lensing — 1× in this paper · explained in 2 episodes
Gravitational Wave Emission — 1× in this paper · explained in 2 episodes
Matter Power Spectrum — 1× in this paper · explained in 2 episodes
Reheating — 1× in this paper · explained in 2 episodes

Summary

The gist The work investigates whether a transient period of matter and kination domination in an early Universe can significantly modify the gravitational-wave spectrum from cosmic strings, providing a potential explanation for signals observed by pulsar timing arrays.

Cosmic String Gravitational Waves and Cosmological History

Cosmic strings are considered relic processes that can produce a stochastic gravitational wave background (SGWB) in the nanohertz regime, and this signal is sensitive to the background expansion history of the UniverseThe resulting SGWB depends not only on string parameters such as the tension, emission spectrum, and intercommutation probability, but also heavily on the background expansion history of the UniverseThis makes explicit that omegaGW(f) is determined jointly by the loop formation history, the loop decay law, and the background expansion history

Matter-Kination Domination Sequence

The paper focuses on a non-standard cosmological history between big bang nucleosynthesis and recombination featuring an early matter-dominated era followed by a kination phaseThis late matter–kination sequence can modify the resulting spectral structure in the nanohertz rangeThe characteristic temperatures TRM and TKR are introduced to describe these transitions, with TMK automatically determined by Eq. (2.19)

Impact on Gravitational Wave Spectra

The modified expansion history deforms the broad cosmic string spectrum into a triangular shape, with asymptotic behavior described by Eq. (3.13)The segment fKR ≲ f ≲ fMK probes loops associated with the kination era, where the background energy density redshifts faster than radiation, leading to omegaGW ∝ f+1 The segment fMK ≲ f ≲ fRM probes loops associated with the preceding matter era, resulting in a negative spectral tilt towards higher frequencies and omegaGW ∝ f-1

Bayesian Analysis and Model Comparison

The analysis employs a Bayesian framework to compare competing models, defined by hypotheses H1 through H4, against a baseline SMBHB hypothesis H0The log Bayes factor ln BDij is used to compare two competing models, Mi and Mj with the dataset DThe decomposition of the log evidence using Eqs. (4.5) and (4.6) enables contributions from model fit and the effective model complexity to be assessed separately across different models The comparison between H3 and H1 in the NG15 dataset indicates substantial support for kination over ΛCDM, yielding ln BNG153 = 1.29 ± 0.06, or equivalently BNG153 = 3.63+0.22−0.21

Conclusion on Observational Constraints

The analysis finds no strong evidence favoring any of the alternative hypotheses over the SMBHB-only baseline H0, in either dataset The kination scenario considered in this work enhances the SGWB spectrum toward the high-frequency end, resulting in a more pronounced blue tilt and thereby providing a better fit compared to its counterpart in ΛCDM In contrast to NG15, the EPTA DR2 New+ data slightly favor a flatter SGWB spectrum omegaGWh2 Nonetheless, within the cosmic string framework, with or without contributions from SMBHBs, both datasets mildly favor kination over ΛCDM The inferred cosmic-string-tension posteriors remain largely consistent between the two cosmologies However, the two PTA datasets contain overlapping pulsars; combining the two likelihoods will double count information associated with the common pulsars, potentially biasing the model constraints and model evidences A statistically consistent joint analysis will therefore require dedicated statistical treatments of the correlations between overlapping datasets The forthcoming International Pulsar Timing Array Third Data Release (IPTA DR3) [89] will combine observations from multiple PTA collaborations and will provide a more comprehensive joint dataset for characterizing the SGWB

How it works

The cosmic string gravitational wave spectrum is calculated using Eq. (3.1), where the loop formation efficiency depends on the background equation of state at tk, commonly approximated by Ceff(tk) ≃ 5.4 for RD, 0.39 for MD, and KD
The characteristic frequencies fKR, fMK, and fRM are determined by the transition temperatures TRM, TMK (determined by Eq. (2.19)), and TKR
The Bayesian evidence is decomposed into posterior-weighted log likelihood and KL divergence terms to assess model support

The paper concludes that a joint analysis of the NG15 and EPTA DR2 NEW+ datasets could, in principle, improve the statistical sensitivity and provide tighter constraints on model parameters and model evidences The next step is for future work to clarify the nature of the SGWB and its implication for the early expansion history

--- Page 30 ---

[89] IPTA collaboration, From a Rumble to a Roar: The International Pulsar Timing Array’s Third Data Release, J. Phys. Conf. Ser. 3177 (2026) 012065.

--- Page 31 ---

[9] K.N. Ananda, C. Clarkson and D. Wands, The Cosmological gravitational wave background from primordial density perturbations, Phys Rev D 75 (2007) 123518 [gr-qc/0612013].

--- Page 32 ---

[81] T.N. Collaboration, Kde representations of the gravitational wave background free spectra present in the nanograv 15-year dataset, June, 2023. 10.5281/zenodo.877 [81].

--- Page 33 ---

[73] D. Foreman-Mackey, D.W. Hogg, D. Lang and J. Goodman, emcee: The MCMC Hammer, Publ Astron Soc Pac 125 (2013) 306 [73].

--- Page 34 ---

[84] A. Mitridate, D. Wright, R. von Eckardstein, T. Schröder, J. Nay, K. Olum et al., PTArcade, 2306.16377 [84].

--- Page 35 ---

[51] J. Buchner, UltraNest - a robust, general purpose Bayesian inference engine, J Open Source Softw. 6 (2021) 3001 [51].

--- Page 36 ---

[72] Y. Xu, K. Takahashi, S. Goriely, M. Arnould, M. Ohta and H. Utsunomiya, NACRE II: an update of the NACRE compilation of charged-particle-induced thermonuclear reaction rates for nuclei with mass number A < 16, Nucl Phys A 918 (2013) 61 [72].

--- Page 37 ---

[40] Y. Cui, M. Lewicki, D.E. Morrissey and J.D. Wells, Cosmic Archaeology with Gravitational Waves from Cosmic Strings, Phys Rev D 97 (2018) 123505 [40].

--- Page 38 ---

[41] S. Blasi, V. Brdar and K. Schmitz, Fingerprint of low-scale leptogenesis in the primordial gravitational-wave spectrum, Phys Rev Res 2 (2020) 043321 [41].

--- Page 39 ---

[5] T. Vachaspati and A. Vilenkin, Gravitational Radiation from Cosmic Strings, Phys Rev D 31 (1985) 3052 [5].

--- Page 40 ---

[6] R.R. Caldwell and B. Allen, Cosmological constraints on cosmic string gravitational radiation, Phys Rev D 45 (1992) 3447 [6].

--- Page 41 ---

[7] C.J. Hogan, Gravitational radiation from cosmological phase transitions, Mon Not Roy Astron Soc 218 (1986) 629 [7].

--- Page 42 ---

[8] M. Kamionkowski, A. Kosowsky and M.S. Turner, Gravitational radiation from first order phase transitions, Phys Rev D 49 (1994) 2837 [8].

--- Page 43 ---

[11] R.T. Co, D. Dunsky, N. Fernandez, A. Ghalsasi, L.J. Hall, K. Harigaya et al., Gravitational wave and CMB probes of axion kination, JHEP 09 (2022) 116 [11].

--- Page 44 ---

[57] E. Aver, D.A. Berg, A.S. Hirschauer, K.A. Olive, R.W. Pogge, N.S.J.

Improvements for AI systems

This paper provides a framework for Bayesian analysis comparing different SGWB sources (SMBHBs vs. cosmic strings) under various cosmological models (ΛCDM vs. kination). The improved AI system can perform rigorous model comparison by calculating Bayes factors, as detailed in Section 4, allowing it to quantify the extent to which such a cosmological history is favored by the data and determine if a model offers a better description of current PTA data than the standard cosmology.
The system can be trained to distinguish between spectral features imprinted by different expansion histories. Specifically, it can analyze how the modified expansion history can significantly modify the cosmic-string-induced gravitational wave spectrum, generating characteristic features in the nanohertz range, allowing it to identify which cosmological sequence (e.g., RD → MD → KD → RD) better matches observed spectral breaks like those inferred from NG15 or EPTA data.
The system can be used to constrain fundamental string parameters and cosmological history simultaneously by incorporating observational priors. It can implement the methodology described in Section 5.1, where BBN and CMB observations constrain TRM and TKR, allowing it to derive constraints on the nonstandard cosmology parameters TRM and TKR based on light-element abundances and CMB anisotropy data before applying them to the PTA signal.
The AI can perform joint likelihood analyses across multiple datasets, such as NG15 and EPTA DR2 New+, while accounting for observational overlap. It can address the issue that the two PTA datasets contain overlapping pulsars; combining the two likelihoods will double count information associated with the common pulsars, by implementing dedicated statistical treatments of the correlations between overlapping datasets to yield a more comprehensive constraint.
The system can generate physically motivated predictions for gravitational wave spectra under non-standard expansion histories. It can calculate how a sequence like RD → MD → KD → RD imprints characteristic breaks and changes of slope in the gravitational wave spectrum, such as the triangular shape described in Eq. (3.13), allowing it to predict spectral enhancements or deviations relative to standard cosmology for specific PTA frequency ranges.