{"id":"396fdccd-6080-4f5c-be4e-3aa5be7f75fe","arxiv_id":"2508.08249","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A transiently dominant rotating axion sources an induced gravitational wave background whose amplitude and large-scale anisotropy trace the axion's isocurvature fluctuations, with detection prospects for BBO and DECIGO.","lead":"A rotating axion that briefly dominates the early universe can generate a strong, anisotropic gravitational wave background, potentially visible to future detectors. The anisotropy, inherited from the axion's quantum fluctuations, would offer a new way to probe physics beyond the Standard Model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3.32) is the load-bearing step: an O(1) error in its δχ coefficient would propagate quartically into the flat GW amplitude and octically into the anisotropy power, directly moving the benchmark detection curves in Figs. 6–9.","rationale":"I agree with the reader that Eq. (3.32) is the weakest load-bearing step. The reason it is load-bearing rather than merely cosmetic is the high power with which the coefficient enters the flat GW amplitude and its anisotropy power: ΩGW,flat is quartic in δχ/Mpl, and the absolute anisotropy power C_l^GW is octic. Thus a plausible O(1) correction from the factor-2 energy-density approximation or a gradual transition could change detectability conclusions. However, I do not think this forces a rejection: the peaked component from kination-era curvature perturbations, Eq. (3.11), is independent of Eq. (3.32), and the appendix already contains a partial numerical check of the order of magnitude. The appropriate verdict therefore remains CONDITIONAL, as the reader assigned. The proposed check would either resolve the concern or quantify the required revision of the benchmark plots.","tokens_in":49281,"tokens_out":33602,"duration_ms":365363,"concrete_test":"Recompute the late-time asymptotic value of √(2ρS)aξS/(Mpl k^2 δS,i) from the fluid equations in Appendix D.1.1 for the log-potential background, using the exact ρS of Eq. (D.4), for several FS values (0.01, 0.1, 10, 100) and for kηMK = 10^-3, requiring the asymptotic value to be stable to <1% over a decade in η. Compare with Eq. (3.32). Then insert the extracted coefficient into Eq. (3.44) and into C_l^GW for the flat component, and re-plot the benchmark lines of Figs. 6–9. If the coefficient deviates from Eq. (3.32) by more than 20%, the flat-spectrum amplitude and anisotropy detectability claims must be revised; if it matches within 20%, the conditional concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (3.32) converts the primordial axion charge fluctuation δS into the late-time superhorizon phase perturbation δχ. This coefficient is the unique input to the flat part of the induced GW spectrum: Eq. (3.44) gives ΩGW(ηc,k) ≃ 0.175 Aδχ^2, and Aδχ is the square of δχ/Mpl from Eq. (3.32). The flat component is the broad-frequency target and the one used for the anisotropy benchmarks in Fig. 7, since C_l^GW for the flat part scales as ΩGW^2 δGW^2 and therefore as (δχ/Mpl)^8. The derivation of Eq. (3.32) uses the approximate ρS in Eq. (3.22), which is a factor 2 smaller than r^2 θdot^2 in the matter phase, and an instantaneous MK transition. The numerical check in D.1.1 is reported only as an order-of-magnitude consistency; it does not establish the coefficient to the accuracy needed for SNR≈1 detectability. A √2 error in the coefficient changes ΩGW,flat by a factor 4 and the anisotropy power by a factor 16; a factor 1.5 error changes them by roughly 5 and 25, respectively, which is enough to move the benchmark curves in Figs. 6–9 relative to the detector noise curves. The anisotropy ratio itself, Eq. (4.15), is independent of the coefficient, but the absolute anisotropy power is not. This is the most load-bearing quantitative uncertainty in the central claim that both the amplitude and the anisotropy are accessible to future detectors.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies gravitational wave production from a rotating axion field whose energy density temporarily grows during a matter-like phase and then dilutes during a kination phase before returning to radiation domination. The authors show that the transient enhancement of the axion-induced curvature perturbation and the superhorizon growth of axion phase fluctuations can source a strong induced stochastic gravitational wave background, with a flat low-frequency part plus a feature near the matter-to-kination transition. They further show that long-wavelength axion fluctuations produce large-scale anisotropies in this background, and they compare the predicted amplitude and anisotropy with projected sensitivities of LISA, BBO, DECIGO, CE, ET, muAres, and SKA, after imposing constraints from PBH formation, axion radiation, isocurvature, CMB B-modes, and non-Gaussianity.","tokens_in":49670,"tokens_out":15007,"duration_ms":157205,"significance":"If the quantitative predictions are robust, this is a significant and novel contribution: it identifies a realistic rotating-axion scenario in which an induced gravitational wave background carries large isocurvature-type anisotropies, a signature that is absent in the usual adiabatic induced-GW calculations. The paper is strong in its use of standard second-order perturbation theory, in providing analytic formulas that are cross-checked numerically in several appendices, in the transparent separate-universe derivation of the anisotropy coefficients, and in the comprehensive treatment of independent observational constraints. The central mechanism, namely that kination dilution evades large-scale CMB bounds while allowing a strong signal at shorter scales, is well motivated and clearly explained. The main weakness is that the key coefficient governing the flat part of the spectrum, Eq. (3.32), is obtained through approximations whose accuracy is quantified only at the order-of-magnitude level, and the detection claims are sensitive to that coefficient.","major_comments":[{"comment":"","section":"Sec. 3.2.1, Eq. (3.32)"},{"comment":"","section":"Sec. 5.5, Eq. (5.30)"},{"comment":"","section":"Sec. 3.2.2, Eq. (3.44) and Figs. 3-4"}],"minor_comments":[{"comment":"","section":"Eqs. (5.19) and (5.30)"},{"comment":"","section":"Sec. 3.2.1"},{"comment":"","section":"Fig. 7"},{"comment":"","section":"Eq. (3.6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of JHEP and the central idea is novel and interesting. My main reservation is not the validity of the formalism but the quantitative calibration of the coefficient in Eq. (3.32), which controls the flat GW amplitude and the absolute anisotropy power used for detection benchmarks. The authors are transparent about this limitation, and the manuscript has many strengths, including the analytic kernels, numerical cross-checks, and a comprehensive constraints analysis. With a more precise or better-quantified treatment of the delta_chi growth coefficient, and with the non-Gaussianity constraint handled at a level consistent with its use in setting benchmarks, I would be supportive of publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nYou should know: Bodas et al. present a well-motivated scenario where a rotating axion field, temporarily dominating around its matter-to-kination transition, generates induced gravitational waves with large-scale anisotropies inherited from axion isocurvature. The central new results are the explicit anisotropy coefficients—δGW = -8 ζS in rotation dominance, and -4(1+2ΩS)ζS for the flat part—and the argument that the transient MK dynamics evades CMB isocurvature bounds while boosting the signal.\n\nWhat the paper does well: the constraint analysis is thorough and honest. They identify the strongest bounds (non-Gaussianity for long-wavelength modes, PBH overproduction and axion radiation isocurvature for short wavelengths) and show that future detectors like BBO and DECIGO could see both the monopole and low-ℓ anisotropy. The numerical checks in Appendix D and the separate universe derivation of the anisotropy are consistent. They also explicitly state where their approximations break down.\n\nThe soft spot is the coefficient in Eq. (3.32), which converts the primordial charge fluctuation δS into the late-time superhorizon phase perturbation δχ. It uses an approximate energy density that is a factor 2 off in the matter phase and assumes an instantaneous transition. The numerical check in D.1.1 confirms only order-of-magnitude consistency. Because the flat GW amplitude scales as (δχ/Mpl)^2 and the anisotropy power as (δχ/Mpl)^8, a √2 error changes ΩGW by about a factor 4 and Cℓ^GW by about 16, which can shift the benchmark curves in Figs. 6–9 relative to detector noise. The qualitative claim—that the signal could be accessible—survives, but the precise reach is uncertain. A systematic estimate of this uncertainty would make the paper stronger.\n\nMinor issues: the Introduction's claim of \"first example\" of large anisotropies in induced GWBs is overstated, since earlier papers studied anisotropies of scalar-induced GWs (e.g., Refs. [33,35,36]). The phrase should be narrowed to the axion-kination context. The artificial jumps at the KR/MK boundaries are acknowledged and not a problem.\n\nWho should read it: anyone working on induced GWs, axion phenomenology, or early-universe probes. It deserves serious refereeing; the central mechanism is plausible, the exposition is clear, and the honest caveats make it a good foundation despite the O(1) uncertainty in the key coefficient. I'd recommend acceptance after minor revision.","headline":"Solid, honest phenomenology that connects rotating axion dynamics to detectable GW anisotropies, but the main quantitative predictions rest on an O(1) coefficient that is not pinned down.","tokens_in":50200,"tokens_out":2843,"would_cite":true,"duration_ms":31262,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that a rotating axion field that briefly dominates the early universe before a kination transition can source an induced gravitational-wave background strong enough for future detectors, with pronounced large-scale…","keywords":["induced gravitational waves","rotating axion","kination","axion isocurvature","superhorizon growth","gravitational wave anisotropy","primordial black holes","early matter domination"],"falsifier":"A future experiment with sensitivity comparable to DECIGO or BBO that searches the $10^{-3}$-$10^{-2}$ Hz band and finds no stochastic background with the predicted broken spectrum for benchmark parameters ($\\zeta_S(k_s)\\sim10^{-2}$, $T_{\\rm KR}\\sim10^3$ TeV) would rule out the strong-signal claim; alternatively, measuring the large-scale anisotropy and finding the peak coefficient different from $-8\\,\\zeta_S(k_l)$, or finding the flat-branch anisotropy not of the form $-4(1+2\\Omega_{S,MK})\\,\\zeta_S(k_l)$, would directly contradict the mechanism.","tokens_in":49067,"feed_emoji":"🌊","tokens_out":9314,"duration_ms":95418,"temperature":0.7,"pith_summary":"The paper argues that a rotating axion field that temporarily takes over the universe's energy budget—first behaving like matter and then like a kinetic-energy fluid, or kination—can convert its own inflationary fluctuations into a detectable stochastic gravitational-wave background. The central mechanism is a transient enhancement: near the matter-to-kination (MK) transition, the axion's curvature perturbation dominates the total, and its phase fluctuations grow on superhorizon scales. These fluctuations later re-enter the horizon and source induced gravitational waves, both around the MK transition and, through a flat low-frequency branch, during the subsequent radiation era. A distinctive property is that the resulting wave background inherits pronounced large-scale anisotropies from the axion fluctuations—at the level $\\delta_{\\rm GW} = -8\\,\\zeta_S(k_l)$ in the rotation-dominated case and about $4\\,\\zeta_S(k_l)$ when the axion stays subdominant. The transient nature of the axion's dominance is what lets the scenario evade the usual CMB bounds on large perturbations.","feed_headline":"Rotating axions may ring a detectable anisotropic GW background","feed_subtitle":"A short early-universe phase amplifies axion fluctuations and stamps their anisotropy on the wave signal.","key_machinery":"The central object is the rotating complex scalar $S = \\frac{1}{\\sqrt{2}} r e^{i\\theta}$ and its conserved U(1) charge yield $Y_\\theta = \\dot\\theta r^2/s$. The matter-to-kination (MK) transition, where the axion energy density switches from $\\rho_S \\propto a^{-3}$ to $\\rho_S \\propto a^{-6}$, temporarily boosts the axion's contribution to the total curvature perturbation and then dilutes it fast enough to evade CMB constraints. The load-bearing identity is the superhorizon growth of the phase fluctuation, Eq. (3.32), which relates the late-time GW source $\\delta\\chi/M_{\\rm pl}$ to the primordial fluctuation $\\delta_S$ with a coefficient of order one (rotation dominance) or $F_S^{1/2}$ (non-dominance). The induced-GW kernels of Eqs. (3.2), (3.13), and (3.36) carry the second-order metric perturbations that generate the stochastic background, and the separate-universe modulation of the local $Y_\\theta$ produces the anisotropy coefficients in Eq. (4.6).","core_discovery":"This paper shows that a complex scalar field with a nearly quadratic radial potential, set rotating by an Affleck-Dine type kick, produces a detectable induced gravitational-wave background if its fluctuations exceed the adiabatic ones from the inflaton. The conserved U(1) charge converts the initial angular fluctuation into a conserved curvature perturbation $\\zeta_S \\simeq \\delta_S/3$, and the matter-to-kination transition temporarily boosts the axion's share of the total energy-momentum, enhancing the total curvature perturbation for a short time. During this phase, superhorizon modes of the angular field grow to $\\delta\\chi/M_{\\rm pl} \\simeq 2\\sqrt{2/3}\\,\\delta_S$ in the rotation-dominated case and to $\\sqrt{6}\\,F_S^{1/2}\\,\\delta_S$ when the axion never dominates. The paper computes the induced GW power spectrum: modes entering during kination give $\\Omega_{\\rm GW}h^2 \\simeq 2.9\\times10^{-5} A_{\\zeta_S}^2 (\\nu/3.9\\times10^{-5}\\,{\\rm Hz})(T_{\\rm KR}/{\\rm TeV})^{-1}$, while modes entering later during radiation domination give a flat branch with $\\Omega_{\\rm GW}(\\eta_c,k)\\simeq 0.175 A_{\\delta\\chi}^2$. The large-scale anisotropy of the background is derived using the separate-universe picture, yielding coefficients such as Eq. (4.6) for the MK-peak and Eq. (4.15) for the late flat branch. Observational constraints from primordial black hole overproduction, axion sound-wave dark radiation, late kination-to-radiation isocurvature, CMB B-modes, and non-Gaussianity are incorporated; short-wavelength fluctuations can reach $\\sim 10^{-2}$, while long-wavelength fluctuations are capped at $\\lesssim 10^{-3}$, mainly by non-Gaussianity.","pith_inferences":["The same superhorizon phase growth that sources gravitational waves also sources axion sound waves, which behave as dark radiation; cross-correlating the GW anisotropy with the imprinted dark-radiation anisotropy could in principle separate $\\zeta_S(k_l)$ from the background cosmology without relying on the two frequency observables.","If the axion spectrum is blue-tilted enough to suppress $\\zeta_S(k_l)$ while keeping $\\zeta_S(k_s) \\sim 10^{-2}$, the scenario avoids the non-Gaussianity cap and the short early matter era could produce primordial black holes that carry no isocurvature; those black holes would be accompanied by a flat, kination-boosted gravitational-wave signal that tests the matter-era hypothesis directly.","The mechanism should apply generically to any spectator sector that passes through a matter-like phase followed by kination, not only to axions from Affleck-Dine dynamics; a null search for the predicted broken, anisotropic spectrum would then constrain the duration of any such early matter era.","Because the anisotropy coefficient at the peak, $\\delta_{\\rm GW} = 4(-11 + \\Omega_{S,MK} + 24/(2+\\Omega_{S,MK}))\\,\\zeta_S(k_l)$, differs from the flat-branch coefficient $-4(1+2\\Omega_{S,MK})\\,\\zeta_S(k_l)$, measuring the multipole spectrum at two frequencies would provide a consistency check of the model; a mismatch would point to a non-instantaneous transition or a different source."],"forward_implications":["If the axion rotation dominates around the MK transition, the flat part of the GW spectrum below $\\nu_{MK}$ and the peak at $\\nu_{MK}$ are within reach of proposed experiments such as BBO and a SQL-limited DECIGO, with low-multipole anisotropies detectable for $\\zeta_S(k_l) \\sim 6\\times10^{-4}$.","Detecting both the anisotropy of the peak amplitude and the modulation of the peak frequency, $\\delta\\nu_{MK}/\\nu_{MK} = (2\\Omega_{S,MK}-1)\\,\\zeta_S(k_l)$, would allow independent extraction of the long-wavelength axion fluctuation $\\zeta_S(k_l)$ and the axion's energy fraction $\\Omega_{S,MK}$, identifying the rotating-axion origin of the background.","The low-frequency flat branch is sourced by modes that re-enter during the radiation era, long after the axion has diluted, so the signal does not require the axion to dominate at late times; the background nevertheless retains the large-scale anisotropy imprinted during the MK epoch.","In the non-dominant regime, the GW amplitude scales as $\\Omega_{\\rm GW} \\propto F_S^2 \\delta_S^4$, which is less suppressed than the naive curvature-based estimate, and the allowed short-wavelength fluctuation can be as large as $\\sim10^{-2}$ while long-wavelength fluctuations are forced below about $10^{-3}$, mostly by the CMB non-Gaussianity bound.","The combination of amplitude, spectral shape, and anisotropy provides a rare observational window into the quantum fluctuations of a spectator axion field during inflation, because the induced background free-streams from early times to the present."],"supporting_citations":[{"why":"Supplies the rotating-axion charge dynamics and the matter-to-kination background evolution used as the cosmological setting.","marker":"[11]"},{"why":"Provides the axion kination cosmology and its gravitational-wave and CMB probes, forming the background model for the GW calculation.","marker":"[22]"},{"why":"Gives the induced-GW calculation from axion rotations right after the transition to kination, the basis for the kination-era spectrum.","marker":"[30]"},{"why":"Provides the semianalytic induced gravitational-wave kernel formula used in Eq. (3.2).","marker":"[27]"},{"why":"Supplies the induced-GW calculation from isocurvature perturbations used for the case without rotation dominance.","marker":"[69]"},{"why":"Establishes the relation between the axion angular fluctuation and the conserved curvature perturbation $\\zeta_S$ used throughout.","marker":"[63]"},{"why":"Provides the mapping from kination-era GW production to the present-day energy-density parameter, including the $H_{\\rm KR}$ normalization.","marker":"[67]"},{"why":"The multi-component anisotropy calculation in this paper goes beyond the single-component consistency relation derived there.","marker":"[76]"},{"why":"Supplies the Planck $f_{\\rm NL}$ constraint that sets the strongest bound on long-wavelength axion fluctuations.","marker":"[92]"}],"fun_headline_variants":["Rotating axions produce anisotropic gravitational wave background","Axion rotation leaves subtle mark on gravitational wave sky","Anisotropic GWs from brief axion rotation phase","Axion spin-up could make gravitational waves anisotropic","Detectable anisotropic gravitational waves from rotating axions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on the estimate that the axion's angular fluctuation grows on superhorizon scales as in Eq. (3.32), which assumes an instantaneous switch from matter-like to kination behavior and a simplified energy-density formula in the matter phase; a gradual transition or a factor-of-two correction there would shift the wave amplitudes and the anisotropy coefficients.","fun_headline_variants_meta":{"raw":{"variants":["Rotating axions produce anisotropic gravitational wave background","Axion rotation leaves subtle mark on gravitational wave sky","Anisotropic GWs from brief axion rotation phase","Axion spin-up could make gravitational waves anisotropic","Detectable anisotropic gravitational waves from rotating axions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000416,"raw_usage":{"total_tokens":2226,"prompt_tokens":1104,"completion_tokens":1122,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":1047}},"tokens_in":720,"tokens_out":1122,"duration_ms":9312,"temperature":1.0,"reasoning_tokens":1047,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:36:11.472124+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future experiment with sensitivity comparable to DECIGO or BBO that searches the $10^{-3}$-$10^{-2}$ Hz band and finds no stochastic background with the predicted broken spectrum for benchmark parameters ($\\zeta_S(k_s)\\sim10^{-2}$, $T_{\\rm KR}\\sim10^3$ TeV) would rule out the strong-signal claim; alternatively, measuring the large-scale anisotropy and finding the peak coefficient different from $-8\\,\\zeta_S(k_l)$, or finding the flat-branch anisotropy not of the form $-4(1+2\\Omega_{S,MK})\\,\\zeta_S(k_l)$, would directly contradict the mechanism.","supporting_citations":[],"review_version":1}