Pith. sign in

REVIEW 3 major objections 4 minor 45 references

The paper argues that vector dark matter from kinetic-coupling misalignment of an isotropized multi-vector condensate is ruled out by mutually incompatible CMB non-Gaussianity and isocurvature bounds.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-02 02:36 UTC pith:J4RLC3BI

load-bearing objection The strong-mixing half of the no-go argument collapses once the isocurvature bound and Eq. (50) are corrected; the weak-mixing half may be right but the paper's own derivation is invalid. the 3 major comments →

arxiv 2607.14267 v1 pith:J4RLC3BI submitted 2026-07-15 astro-ph.CO gr-qc

Misalignment production of isotropized vector dark matter?

classification astro-ph.CO gr-qc
keywords vector dark mattermisalignment mechanismkinetic couplingisocurvature fluctuationsnon-Gaussianityinflationprimordial perturbationsmulti-vector condensate
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish a no-go result: vector dark matter cannot be produced by the misalignment mechanism of a multi-vector condensate that is kinetically coupled to the inflaton and isotropized to avoid anisotropy constraints. It shows that the same mixing parameter h is bounded from above by curvature non-Gaussianity and from below by isocurvature fluctuations, in both the weak-mixing (h << 1) and strong-mixing (h >> 1) regimes. Specifically, non-Gaussianity requires h < 3 x 10^-4 (weak) or h < 9.2 (strong), while isocurvature requires h > 5.2 x 10^-3 (weak) or h > 45 (strong). Since no value of h satisfies both, the production channel is excluded. A reader should care because this closes off a seemingly natural inflationary route to vector dark matter and points to isocurvature overproduction as the generic obstruction.

Core claim

The paper's central claim is that isotropized vector dark matter from kinetic-coupling misalignment fails against CMB data. The mixing parameter h, defined as the ratio of vector to inflaton kinetic energies, governs all perturbation effects. For h << 1, curvature and bispectrum grow as h^2 N^2 and h^2 N^3 while the entropy spectrum shrinks as 1 - (56/3)h^2 N^2, so non-Gaussianity demands very small h and isocurvature demands larger h. For h >> 1, curvature is exponentially amplified (e^{2.37h}), the local bispectrum decreases with h, and the entropy fluctuation is anti-correlated with power ~ (2/h^2) P_R, again leaving no overlap in allowed h. No value of h survives both bounds.

What carries the argument

The key object is h, the ratio of vector-field kinetic energy to inflaton kinetic energy, which sets the coupling between inflaton and vector perturbations. Weak mixing (h << 1) produces cumulative superhorizon modes with h^2 N^2 and h^2 N^3 growth in the curvature spectrum and bispectrum and a -(56/3)h^2 N^2 correction in the entropy spectrum; strong mixing (h >> 1) produces constant modes with exponential curvature amplification e^{2.37h}, local f_NL = (5/6)(7 - 2h - 6/h), and anti-correlated entropy power about (2/h^2) P_R. The conflict is that non-Gaussianity bounds h from above while isocurvature bounds h from below in both regimes.

Load-bearing premise

The no-go rests on the fluctuation spectra and bispectra quoted from the author's earlier work—especially the weak-mixing entropy formula and the strong-mixing relations—being correct and valid at the h values used; if those formulas are wrong or must be used non-perturbatively, the incompatibility could disappear.

What would settle it

Recompute the weak-mixing entropy power spectrum beyond leading order in h N_k. If the coefficient 56/3 is modified or the linear expansion breaks down for h N_k > sqrt(3/56), the isocurvature lower bound h > 5.2 x 10^-3 may evaporate and a window with h < 3 x 10^-4 could open. Alternatively, a future CMB measurement of the anti-correlated isocurvature fraction beta that finds beta < 10^-5 would falsify the strong-mixing branch's h > 45 requirement.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • In the weak-mixing regime, no h can satisfy both the non-Gaussianity bound h < 3 x 10^-4 and the isocurvature bound h > 5.2 x 10^-3, so this production channel is closed.
  • In the strong-mixing regime, the local bispectrum bound h < 9.2 conflicts with the anti-correlated isocurvature bound h > 45, closing that branch as well.
  • As a result, kinetically coupled isotropized vector condensates cannot be the dark matter, unless a fraction of the vector fields acts as a curvaton rather than as dark matter—the loophole the paper notes.
  • In the strong-mixing regime, the case where the vector starts oscillating before reheating is already excluded by the required inflationary scale, independent of the h conflict.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the conflict traces to the sign of the entropy-power correction, so a similar no-go is likely to hold for other inflationary vector-production mechanisms whose entropy spectrum scales with the same mixing parameter; a model-independent proof would be a natural next step.
  • Beyond the paper: a non-perturbative treatment of the weak-mixing entropy spectrum is the most direct way to test the exclusion, since the linear formula is used near its breakdown point h N_k ~ sqrt(3/56).
  • Beyond the paper: if future CMB experiments detect local non-Gaussianity near current limits, the strong-mixing upper bound h < 9.2 becomes a firmer target and the inconsistency with h > 45 would sharpen the no-go.
  • Beyond the paper: a lattice or numerical simulation of the coupled inflaton-vector system during inflation could independently confirm the quoted spectra and either rescue or bury the mechanism.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes dark matter production via the misalignment mechanism of a multi-vector condensate with kinetic coupling during inflation, using an isotropized vector background A_i = A(t) δ_{ai} to avoid the anisotropy problem. It derives the attractor background, the ratio h of vector to inflaton kinetic energy, and the relic abundance of the resulting vector dark matter. The paper then applies CMB constraints on curvature perturbations, non-Gaussianities, and isocurvature fluctuations, separately in the weak-mixing (h << 1) and strong-mixing (h >> 1) regimes. The central claim is that the non-Gaussianity and isocurvature constraints are incompatible in both regimes, so the model is excluded regardless of h. The decisive perturbation spectra and bispectra are quoted from the author's previous papers [27-30], including the unpublished preprint [29], rather than derived in the present manuscript.

Significance. If established, the claimed no-go would be a useful and nontrivial exclusion for isotropized vector dark matter from kinetic-coupling misalignment. The background and relic-abundance parts are clearly organized, and the paper is candid about its assumptions. However, the central no-go is not currently supported: the strong-mixing isocurvature bound contains an internal inconsistency and a misquoted observational limit, and the weak-mixing isocurvature bound is used outside its perturbative validity. The paper therefore needs substantial revision before its main conclusion can be accepted; after correction, the conclusion may even reverse.

major comments (3)
  1. [§IV.C.2, Eqs. (49)-(51)]
  2. [§IV.C.1, Eq. (47)]
  3. [Eqs. (32), (37), (40), (42), (47), (49), (50)]
minor comments (4)
  1. [§II.A, Eqs. (12)-(13)] The sentence 'If we assume that the energy density of vector fields is negligible...' is repeated almost verbatim before Eq. (13). Remove the duplicate.
  2. [Throughout] Typos and wording: 'CONSTRAINS' in the Section IV title; 'scale-invarince' in the Fig. 2 caption; 'bistpecrum' before Eq. (42); 'Guassianities' in Section V; 'So for we have found' in Section V; 't is also expected' after Eq. (46).
  3. [Figs. 2-4] The figures would benefit from direct axis labels and legends describing the shaded regions; the text alone ('orange region', 'gray region') is difficult to follow. Please clarify which constraints correspond to each shaded region.
  4. [Eqs. (17), (31)] The symbol R is used both for the energy-density ratio ρ_A/ρ_φ in Eq. (17) and for the curvature perturbation in Eq. (31) and thereafter. This notational collision is confusing in Section IV and should be fixed (e.g., by using a different symbol for the energy-density ratio).

Circularity Check

0 steps flagged

No significant circularity: the no-go argument applies previously derived perturbation results to a new model and combines them with external CMB bounds.

full rationale

The paper's central exclusion claim is not circular by construction. The background dynamics (Sec. II) and relic abundance formulas (Sec. III) are derived within the paper from the stated action. The incompatibility conclusion in Secs. IV–V is obtained by combining two kinds of inputs: external CMB limits from Planck/BICEP (refs [43,44]) and perturbation spectra/bispectra quoted from the author's earlier papers [27–30] (Eqs. 32, 37, 40, 42, 47, 49). These quoted formulas are load-bearing, but they do not assume the conclusion. They are prior computations of the same model's fluctuations under stated slow-roll and mixing assumptions, not fits to the present no-go result; they are parameter-free in the sense that no parameter is adjusted to make the incompatibility appear. Applying one's own previously derived formulas to a new phenomenological question is not circular. The weak-mixing isocurvature spectrum Eq. (47) is also cited from independent work [37]. Possible internal inconsistencies in the strong-mixing isocurvature branch (Eq. (49) vs Eq. (50), and whether h>45 follows from beta<1e-4) are arithmetic or validity objections, not circularity: even if the strong-mixing branch is miscomputed, no step defines its conclusion into its inputs. No fitted parameter is relabeled as a prediction; h and H_inf are free parameters scanned against constraints. Therefore no specific reduction of a predicted quantity to an input quantity can be exhibited, and the appropriate finding is no significant circularity.

Axiom & Free-Parameter Ledger

3 free parameters · 6 axioms · 0 invented entities

The paper introduces no new particles or forces; the condensate is a configuration of known vector fields. The free parameters are the mixing strength h, the inflationary Hubble scale H_inf, and the benchmark reheating temperature T_reh. The main fragility is the unstated reliance on the author's earlier perturbation results, which are not independently derived here and are used inconsistently in places (Eq. 47 negativity; Eq. 49 vs Eq. 50).

free parameters (3)
  • h (mixing strength)
    Free parameter controlling perturbation mixing; related to the kinetic-coupling constant c via h ≈ sqrt((c-1)/2) on the attractor. The paper scans h and derives bounds from CMB constraints.
  • H_inf (inflationary Hubble scale)
    Treated as a free parameter when converting R to vector mass (Eqs. 28, 30, 34, 38); constrained only by the tensor-to-scalar ratio r < 0.036.
  • T_reh (reheating temperature) = 10^12 GeV (fixed for plots)
    Fixed to 10^12 GeV for quantitative constraints and figures; the paper states 'We also fixed T_reh = 10^12 GeV'. Not fitted, but an arbitrary benchmark.
axioms (6)
  • domain assumption The isotropized configuration A_i^(a) = A(t) δ_ai with more than two vector fields is a dynamical attractor.
    Adopted from ref. [24]; motivates replacing anisotropic single-field vector hair with an isotropic multi-vector hair.
  • domain assumption Vector bosons are light during inflation, m_A/f << H, so the mass term is neglected and af^2 ∂_t A = p_A is constant.
    Stated before Eq. (11); required for the attractor solution.
  • domain assumption The back-reaction attractor (f ∝ a^{-2}, ϕ˙ ≈ -V_φ/(3cH), h constant) is reached before CMB modes exit, i.e., N_back > 60.
    Section IV preamble; all perturbation results assume this attractor.
  • domain assumption Instantaneous reheating with f(ϕ)→1 at the end of inflation (f_e = 1), so the vector energy density at a_e is the initial condition for the post-inflationary scaling.
    Section III: 'we assume that at a_e we have f_e = 1 immediately'.
  • ad hoc to paper The fluctuation spectra and bispectra quoted from Refs. [27-30] (Eqs. 32, 37, 40, 42, 43, 47, 49) are correct and apply to this model; the weak-mixing series remains valid at the h values used.
    These formulas are not derived in this paper. Eq. (47) becomes negative at h N_k ≈ 0.23, so the linear correction is used beyond its regime of validity for the isocurvature bound.
  • standard math Leading-order slow-roll approximation and the Planck/BICEP observational bounds are adopted as stated.
    Section IV; standard practice for inflationary model constraints.

pith-pipeline@v1.3.0-alltime-deepseek · 12599 in / 20605 out tokens · 207409 ms · 2026-08-02T02:36:34.273925+00:00 · methodology

0 comments
read the original abstract

We present dark matter production by the misalignment mechanism of a multi-vector condensate through kinetic coupling during inflation. We impose isotropized background vector fields to release the model from the stringent constraint of anisotropy. However, it turns out that the constraints imposed by non-Gaussianity and isocurvature fluctuations are incompatible with each other, regardless of whether the fluctuations are in the weak-mixing or strong-mixing regime.

Figures

Figures reproduced from arXiv: 2607.14267 by Chong-Bin Chen.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the evolution of [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Constraints on the curvature fluctuation for [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Constraints on the curvature fluctuation for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Constraints on the curvature fluctuation for [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 40 linked inside Pith

  1. [1]

    P. W. Graham, J. Mardon, and S. Rajendran, Vector Dark Matter from Inflationary Fluctuations, Phys. Rev. D93, 103520 (2016), arXiv:1504.02102 [hep-ph]

  2. [2]

    orthogo- nal

    +O H 2 m2 A .(24) Therefore, in the deep coherent oscillation regime, the vector field behaves as non-relativistic matter with ⟨ρA⟩ ∝a−3, and can thus serve as a viable dark mat- ter candidate. 4 III. RELIC ABUNDANCE OF DM After inflation ends, the universe begins reheating. We assume that reheating ends at timea reh, which is after the end of inflation a...

  3. [3]

    ¨Ozsoy and G

    O. ¨Ozsoy and G. Tasinato, Vector dark matter, inflation, and non-minimal couplings with gravity, JCAP06, 003, arXiv:2310.03862 [astro-ph.CO]

  4. [4]

    Bastero-Gil, J

    M. Bastero-Gil, J. Santiago, L. Ubaldi, and R. Vega- Morales, Vector dark matter production at the end of inflation, JCAP04, 015, arXiv:1810.07208 [hep-ph]

  5. [5]

    Bastero-Gil, J

    M. Bastero-Gil, J. Santiago, R. Vega-Morales, and L. Ubaldi, Dark photon dark matter from a rolling in- flaton, JCAP02(02), 015, arXiv:2103.12145 [hep-ph]

  6. [6]

    R. T. Co, A. Pierce, Z. Zhang, and Y. Zhao, Dark Photon Dark Matter Produced by Axion Oscillations, Phys. Rev. D99, 075002 (2019), arXiv:1810.07196 [hep-ph]

  7. [7]

    J. A. Dror, K. Harigaya, and V. Narayan, Parametric Resonance Production of Ultralight Vector Dark Matter, Phys. Rev. D99, 035036 (2019), arXiv:1810.07195 [hep- ph]

  8. [8]

    Zhang, P

    H.-Y. Zhang, P. Arias, A. Cheek, E. D. Schiappacasse, L. Visinelli, and L. Roszkowski, Dark photon dark mat- ter from flattened axion potentials, JHEP10, 142, arXiv:2507.20484 [hep-ph]

  9. [9]

    Adshead, K

    P. Adshead, K. D. Lozanov, and Z. J. Weiner, Dark pho- ton dark matter from an oscillating dilaton, Phys. Rev. D107, 083519 (2023), arXiv:2301.07718 [hep-ph]

  10. [10]

    A. J. Long and L.-T. Wang, Dark Photon Dark Matter from a Network of Cosmic Strings, Phys. Rev. D99, 063529 (2019), arXiv:1901.03312 [hep-ph]

  11. [11]

    Kitajima and K

    N. Kitajima and K. Nakayama, Dark photon dark matter from cosmic strings and gravitational wave background, JHEP08, 068, arXiv:2212.13573 [hep-ph]

  12. [12]

    Nakai, R

    Y. Nakai, R. Namba, and Z. Wang, Light Dark Photon Dark Matter from Inflation, JHEP12, 170, arXiv:2004.10743 [hep-ph]

  13. [13]

    Firouzjahi, M

    H. Firouzjahi, M. A. Gorji, S. Mukohyama, and B. Sale- hian, Dark photon dark matter from charged inflaton, JHEP06, 050, arXiv:2011.06324 [hep-ph]

  14. [14]

    Salehian, M

    B. Salehian, M. A. Gorji, H. Firouzjahi, and S. Muko- hyama, Vector dark matter production from inflation with symmetry breaking, Phys. Rev. D103, 063526 (2021), arXiv:2010.04491 [hep-ph]

  15. [15]

    Nakai, R

    Y. Nakai, R. Namba, and I. Obata, Peaky produc- tion of light dark photon dark matter, JCAP08, 032, arXiv:2212.11516 [hep-ph]

  16. [16]

    A. E. Nelson and J. Scholtz, Dark Light, Dark Matter and the Misalignment Mechanism, Phys. Rev. D84, 103501 (2011), arXiv:1105.2812 [hep-ph]

  17. [17]

    Nakayama, Vector Coherent Oscillation Dark Matter, JCAP10, 019, arXiv:1907.06243 [hep-ph]

    K. Nakayama, Vector Coherent Oscillation Dark Matter, JCAP10, 019, arXiv:1907.06243 [hep-ph]

  18. [18]

    Nakayama, Constraint on Vector Coherent Oscilla- tion Dark Matter with Kinetic Function, JCAP08, 033, arXiv:2004.10036 [hep-ph]

    K. Nakayama, Constraint on Vector Coherent Oscilla- tion Dark Matter with Kinetic Function, JCAP08, 033, arXiv:2004.10036 [hep-ph]

  19. [19]

    Kitajima and K

    N. Kitajima and K. Nakayama, Viable vector coherent oscillation dark matter, JCAP07, 014, arXiv:2303.04287 [hep-ph]

  20. [20]

    Kaneta, H.-S

    K. Kaneta, H.-S. Lee, J. Lee, and J. Yi, Misalignment mechanism for a mass-varying vector boson, JCAP09, 017, arXiv:2306.01291 [astro-ph.CO]

  21. [21]

    Fujita, K

    T. Fujita, K. Murai, K. Nakayama, and W. Yin, Mis- alignment production of vector boson dark matter from axion-SU(2) inflation, JCAP04, 007, arXiv:2312.06889 [hep-ph]

  22. [22]

    M. a. Watanabe, S. Kanno, and J. Soda, Inflationary Universe with Anisotropic Hair, Phys. Rev. Lett.102, 191302 (2009), arXiv:0902.2833 [hep-th]

  23. [23]

    M. a. Watanabe, S. Kanno, and J. Soda, The Nature of Primordial Fluctuations from Anisotropic Inflation, Prog. Theor. Phys.123, 1041 (2010), arXiv:1003.0056 [astro-ph.CO]

  24. [24]

    Kanno, J

    S. Kanno, J. Soda, and M. a. Watanabe, Anisotropic Power-law Inflation, JCAP12, 024, arXiv:1010.5307 [hep-th]

  25. [25]

    Yamamoto, M

    K. Yamamoto, M. a. Watanabe, and J. Soda, Inflation with Multi-Vector-Hair: The Fate of Anisotropy, Class. Quant. Grav.29, 145008 (2012), arXiv:1201.5309 [hep- th]

  26. [26]

    Maleknejad, M

    A. Maleknejad, M. M. Sheikh-Jabbari, and J. Soda, 10 Gauge Fields and Inflation, Phys. Rept.528, 161 (2013), arXiv:1212.2921 [hep-th]

  27. [27]

    Chen and J

    C.-B. Chen and J. Soda, Anisotropic hyperbolic inflation, JCAP09, 026, arXiv:2106.04813 [hep-th]

  28. [28]

    Chen and J

    C.-B. Chen and J. Soda, Geometric structure of multi- form-field isotropic inflation and primordial fluctuations, JCAP05(05), 029, arXiv:2201.03160 [hep-th]

  29. [29]

    Chen, Inflation with vector fields revisited: heavy entropy perturbations and primordial black holes, JCAP 11, 063, arXiv:2312.06105 [astro-ph.CO]

    C.-B. Chen, Inflation with vector fields revisited: heavy entropy perturbations and primordial black holes, JCAP 11, 063, arXiv:2312.06105 [astro-ph.CO]

  30. [30]

    Chen, Inflation with vector fields revisited: non- Gaussianities, (2026), arXiv:2605.28752 [hep-th]

    C.-B. Chen, Inflation with vector fields revisited: non- Gaussianities, (2026), arXiv:2605.28752 [hep-th]

  31. [31]

    Chen, B.-X

    C.-B. Chen, B.-X. An, and F.-W. Shu, Primordial black holes with anisotropic hair, Phys. Rev. D112, 10 (2025), arXiv:2507.16807 [astro-ph.CO]

  32. [32]

    M. C. Bento, O. Bertolami, P. V. Moniz, J. M. Mourao, and P. M. Sa, On the cosmology of massive vector fields with SO(3) global symmetry, Class. Quant. Grav.10, 285 (1993), arXiv:gr-qc/9302034

  33. [33]

    Golovnev, V

    A. Golovnev, V. Mukhanov, and V. Vanchurin, Vector Inflation, JCAP06, 009, arXiv:0802.2068 [astro-ph]

  34. [34]

    Maleknejad and M

    A. Maleknejad and M. M. Sheikh-Jabbari, Non-Abelian Gauge Field Inflation, Phys. Rev. D84, 043515 (2011), arXiv:1102.1932 [hep-ph]

  35. [35]

    Maleknejad and M

    A. Maleknejad and M. M. Sheikh-Jabbari, Gauge-flation: Inflation From Non-Abelian Gauge Fields, Phys. Lett. B 723, 224 (2013), arXiv:1102.1513 [hep-ph]

  36. [36]

    Demozzi, V

    V. Demozzi, V. Mukhanov, and H. Rubinstein, Magnetic fields from inflation?, JCAP08, 025, arXiv:0907.1030 [astro-ph.CO]

  37. [37]

    Firouzjahi, M

    H. Firouzjahi, M. A. Gorji, S. A. Hosseini Mansoori, A. Karami, and T. Rostami, Charged Vector Inflation, Phys. Rev. D100, 043530 (2019), arXiv:1812.07464 [hep- th]

  38. [38]

    M. A. Gorji, S. A. Hosseini Mansoori, and H. Firouz- jahi, Inflation with multiple vector fields and non- Gaussianities, JCAP11, 041, arXiv:2008.08195 [astro- ph.CO]

  39. [39]

    Ratra, Expressions for linearized perturbations in a massive scalar field dominated cosmological model, Phys

    B. Ratra, Expressions for linearized perturbations in a massive scalar field dominated cosmological model, Phys. Rev. D44, 352 (1991)

  40. [40]

    Hwang, Roles of a coherent scalar field on the evolu- tion of cosmic structures, Phys

    J.-c. Hwang, Roles of a coherent scalar field on the evolu- tion of cosmic structures, Phys. Lett. B401, 241 (1997), arXiv:astro-ph/9610042

  41. [41]

    Hwang and H

    J.-c. Hwang and H. Noh, Axion as a Cold Dark Matter candidate, Phys. Lett. B680, 1 (2009), arXiv:0902.4738 [astro-ph.CO]

  42. [42]

    Kofman, A

    L. Kofman, A. D. Linde, and A. A. Starobinsky, Towards the theory of reheating after inflation, Phys. Rev. D56, 3258 (1997), arXiv:hep-ph/9704452

  43. [43]

    Braden, L

    J. Braden, L. Kofman, and N. Barnaby, Reheating the Universe After Multi-Field Inflation, JCAP07, 016, arXiv:1005.2196 [hep-th]

  44. [44]

    Akramiet al.(Planck), Planck 2018 results

    Y. Akramiet al.(Planck), Planck 2018 results. X. Con- straints on inflation, Astron. Astrophys.641, A10 (2020), arXiv:1807.06211 [astro-ph.CO]

  45. [45]

    P. A. R. Adeet al.(BICEP, Keck), Improved Con- straints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season, Phys. Rev. Lett.127, 151301 (2021), arXiv:2110.00483 [astro-ph.CO]