{"id":"07114617-b07c-49bf-a165-a20f7f538556","arxiv_id":"2506.13237","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Charged cosmic rays spiraling in galactic magnetic fields produce a stochastic gravitational wave background peaking at a few millihertz, though at undetectably low amplitude.","lead":"Cosmic rays moving through the Milky Way's magnetic field should emit a faint hum of gravitational waves, and this paper calculates its spectrum for the first time. The signal peaks near millihertz frequencies, in the range LISA will listen to, but is far too weak to ever detect.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The monochromatic-emission assumption and the ζ≈5 interpolation are coupled: the mHz peak location depends on the assumed spectrum shape, but this can be tested concretely against the published spectral formula.","rationale":"The reader identified the same load-bearing concern: the phenomenological interpolation and the monochromatic emission assumption are the weakest links in the argument. The stress-test confirms this is the right target because the headline result—the mHz peak in Ω(f)—is a convolution of the CR spectrum with a per-particle power that is both approximated (ζ≈5) and assumed monochromatic. The paper's own caveat about alternative spectral viewpoints (Refs. [26,27]) strengthens this. No internal contradiction was found in the derivation; the issue is a robustness/uncertainty gap rather than a demonstrated error. Given the authors' pedagogical goal and their explicit 'approximation' framing, the conditional verdict is appropriate: the qualitative claim is plausible, but the quantitative peak and shape would benefit from a full-spectrum computation. I do not see grounds for rejection since the authors are transparent about the approximations and the detectability claim is already extremely modest.","tokens_in":9515,"tokens_out":1481,"duration_ms":14125,"concrete_test":"Take the spectral formula of Eq. (16) (from Ref. [16]) for dP/dω, and recompute the galactic Ω(f) without the monochromatic assumption: instead of P(f)δ(ν−f_GW(E)), use the full spectral density dP/dν(E,ν) and integrate over the CR energy distribution. If the resulting peak shifts by less than a factor of 2 and the amplitude changes by less than an order of magnitude, the reader's conditional verdict stands; if the peak moves out of the mHz band or the shape changes qualitatively, the headline claim needs revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central prediction is the location and shape of the Ω(f) peak (≈5 mHz for electrons, ≈2 mHz for protons). Two assumptions control this: (i) that all GW power is emitted at f_GW = 2f_B = qBc²/(πE), because the spectrum is monotonically decreasing and dominated by the fundamental mode; and (ii) that the phenomenological power P ≈ ζ G c m² β⁶ γ⁴/R² with ζ≈5 (Eq. 19) is accurate across the intermediate regime. These are not independent: Eq. (22) and the Ω(f) convolution use Eq. (19) as a monochromatic power at each frequency, so if the true spectrum is broad (even if monotonically decreasing), or if Eq. (19) mis-estimates the power at the energies that dominate Ω(f), the peak position and amplitude shift. The paper itself flags that the spectrum result of Refs. [26,27] represents an alternative view, so the monochromatic assumption is acknowledged as a modeling choice rather than a derived result. The weakest point is therefore the unquantified error on the peak frequency and Ω amplitude from combining a phenomenological power law with a single-frequency replacement of a continuous spectrum.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript investigates gravitational radiation emitted by cosmic-ray electrons and protons spiraling in a uniform magnetic field. It derives the non-relativistic and ultra-relativistic power formulas for a charged particle, proposes a phenomenological interpolation P ≈ ζ G c m² β⁶ γ⁴ / R² with ζ ≈ 5, replaces the gravitational-wave spectrum by a monochromatic line at f_GW = qBc²/(πE), and convolves the result with galactic and extra-galactic cosmic-ray spectra. The resulting stochastic background Ω(f) is found to peak at ≈5 mHz for electrons and ≈2 mHz for protons, within the LISA band but with amplitudes far below any foreseeable sensitivity. The paper also provides a pedagogical comparison with electromagnetic synchrotron radiation and emphasizes that the signal is a guaranteed but unobservable background.","tokens_in":9719,"tokens_out":30368,"duration_ms":306162,"significance":"If the monochromatic-spectrum approximation is valid, the paper identifies a new, physically guaranteed stochastic gravitational-wave background from known astrophysical populations, with a non-trivial spectral shape that happens to peak in the LISA band. The manuscript is honest about its limitations (constant B, uniform CR density, isotropic emission, alternative spectral views) and is clearly written. It provides explicit analytic formulas that can be checked and reused, and it correctly stresses that measurable predictions are not the only purpose of a calculation. The main claim, however, hinges on a delta-function replacement of the actual gravitational-wave spectrum (Eq. (16)), and the paper does not quantify the error introduced by that replacement; this should be addressed before the peak frequencies can be considered robust.","major_comments":[{"comment":"The peak frequency and shape of Ω(f), which are the paper's central claims, are obtained by assuming that all gravitational power of a particle of energy E is emitted at the single frequency f_GW = qBc²/(πE). The justification given is that the spectrum of Eq. (16) is \"monotonically decreasing ... and thus dominated by the fundamental mode\". This inference is not demonstrated quantitatively. For the electrons that dominate the peak (E ≈ 3.4 GeV, γ ≈ 6.7×10³), the characteristic frequency in Eq. (16) is ω_c = γ³ω_B ≈ 3×10¹¹ ω_B, so \"dominated by the fundamental\" must mean an enormous suppression of the high-frequency tail; a spectrum that is merely monotonically decreasing could still have most of its power at frequencies far above the fundamental. I ask the authors to quantify, using Eq. (16) or an equivalent published result, the fraction of the total power radiated within a logarithmic interval around the fundamental for the energies that determine the peak, and to state how f_peak and the shape of Ω(f) change when the exact spectral distribution is used instead of the delta-function approximation. The existence of an alternative viewpoint (Refs. [26,27]) makes this check essential. Without it, the headline \"maximum around a few mHz\" is not yet established.","section":"Section III, Eqs. (16)-(22); Section VI, Eq. (28)"}],"minor_comments":[{"comment":"The sentence \"f_GW = 4πωB\" is dimensionally inconsistent and contradicts the later definition f_GW = qBc²/(πE). Since ω_B is an angular frequency, the correct GW frequency is f_GW = ω_B/π (equivalently ω_GW = 2ω_B).","section":"Section II"},{"comment":"The Jacobian relation as written, d f = qBc²/(2π f²) dE, is not consistent with f_GW = qBc²/(πE); differentiating that mapping gives dE/df = qBc²/(π f²). The factor 2 and the direction of the relation should be corrected, and the text should clarify whether f in this equation denotes the GW frequency or the orbital frequency f_B = f_GW/2.","section":"Section V, Eq. (27)"},{"comment":"Eq. (16) uses Φ₂(y) before it is defined in Eq. (17), and the brackets \"3x−1/3\" appear garbled. Please reorder the equations and typeset the Airy-function terms unambiguously.","section":"Section III, Eqs. (16)-(17)"},{"comment":"The paper should specify the value of R_max used for the galactic component, and it should state explicitly how the flux dP/dSdf of Eq. (28) is converted into the energy density ρ used in Eq. (29), including the appropriate geometric factor (e.g., 1/c versus 4π/c for an isotropic background).","section":"Section VI, Eq. (28) and Eq. (29)"},{"comment":"The proton parametrization is said to differ from that in Ref. [23], but no reference or derivation is given for the modified form; please provide one.","section":"Section V, Eqs. (25)-(26)"},{"comment":"There are several typographical errors, including \"inflationnary\" (Introduction), \"synchroton\" (Conclusion), and \"AKNOWLEDGEMENTS\" (Section IX).","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the physics is interesting. The make-or-break point is the monochromatic approximation: if the authors can demonstrate from Eq. (16) that the fundamental mode indeed carries the dominant power, or at least quantify how much the peak shifts when the full spectrum is used, the paper should be publishable. The factor-2 Jacobian inconsistency in Eq. (27) is a concrete error that must be fixed, even though it does not change the peak frequency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a clean completeness result, not a discovery. The genuinely new part is the first computation of the stochastic gravitational-wave background from cosmic-ray electrons and protons gyrating in the Galactic magnetic field, with Ω(f) peaking near 5 mHz for electrons and 2 mHz for protons. The amplitude is far below any conceivable detector, and the authors say so up front. I agree with the reader: the qualitative conclusion is very likely correct, and the calculation is worth taking seriously.\n\nWhat the paper does well: it carefully separates the non-relativistic and ultra-relativistic limits of gravitational synchrotron power, points out the non-monotonic artifact in Eq. (15), and gives a simple all-regime interpolation, Eq. (19), that matches both limits. The observation that the strain is independent of B for fixed speed is counterintuitive and correct. Rewriting the interstellar CR spectra as functions of GW frequency is the right move, and the convolution leading to Ω(f) is laid out clearly. The citation pattern is honest: the difficult relativistic calculation is credited to [13], the spectrum to [16], and they explicitly correct the γ-dependence used in the storage-ring paper [15]. The authors also flag the limitations of a constant 6 μG field, isotropic emission, and the alternative spectral viewpoint in Refs [26,27].\n\nThe soft spots are real but manageable. The monochromatic assumption — all power at f_GW = 2f_B — is doing much of the work in fixing the peak location. The spectrum from Eq. (16) is monotonically decreasing, so the assumption is plausible, but it is an assumption. Because the paper itself flags an alternative view, a referee should ask for a quantitative check: use the published spectral formula, convolve it with the CR spectrum, and show the peak does not move much. I suspect it won't move by orders of magnitude, but the paper does not demonstrate that. The phenomenological ζ ≈ 5 is a multiplicative constant, so it does not affect the peak frequency, but it sets the amplitude and no uncertainty is attached. There are also minor notation slips (Eq. (27), and the f_GW = 4πω_B line), and no code or tabulated data. None of this undermines the central result, because the amplitude is unobservable and the shape is the point.\n\nWho is this for? People who care about the completeness of the SGWB source catalogue and readers who enjoy seeing known physics applied cleanly. It deserves peer review, not a desk rejection, with a request for the spectrum-shape check and a notation cleanup. I would not cite it in my own work, but I would not object to it appearing.","headline":"A clean completeness calculation of a new SGWB source — tiny amplitude, plausible mHz peak — that should be refereed, with one requested spectrum-shape check.","tokens_in":10279,"tokens_out":5675,"would_cite":false,"duration_ms":55872,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.30.-w"],"model":"deepseek-v4-flash","headline":"Cosmic rays spiraling in galactic magnetic fields produce a gravitational-wave background whose energy density peaks at a few millihertz, far too weak to detect but with a shape set by the cosmic-ray spectrum.","keywords":["gravitational waves","stochastic background","cosmic rays","synchrotron radiation","gravitational synchrotron radiation","LISA","millihertz"],"falsifier":"Recompute the background using the full gravitational synchrotron spectrum of Eq. (16) instead of the monochromatic approximation, and compare the resulting Ω(f) peak with the claimed 5 mHz for electrons and 2 mHz for protons; if the peak moves by more than a factor of a few, the central claim fails. A second falsifier would be a future detector with sensitivity around $10^{-47}$ J $m^{-2}$ in the mHz band, which could either find the predicted background or rule it out.","tokens_in":9293,"feed_emoji":"🌌","tokens_out":5856,"duration_ms":59604,"temperature":0.7,"pith_summary":"This paper argues that cosmic rays, the charged particles moving through the Galactic magnetic field, do not only emit synchrotron photons but also gravitational radiation, and that their collective emission forms a stochastic gravitational-wave background. Using a phenomenological formula for the radiated power that interpolates between the non-relativistic and ultra-relativistic regimes, the authors compute the normalized energy density Ω(f) for both electrons and protons. They find that Ω(f) peaks around 5 mHz for electrons and 2 mHz for protons, inside the LISA band, even though the amplitude is far below any foreseeable sensitivity. The shape of the spectrum is non-trivial and encodes the interstellar cosmic-ray spectrum, rewritten as a function of gyration frequency. The paper also points out that for hypothetical ultra-heavy charged particles, gravitational radiation could dominate over electromagnetic synchrotron radiation.","feed_headline":"Cosmic-ray background of gravitational waves peaks near 5 mHz","feed_subtitle":"Even though far too weak to detect, the predicted peak falls in the LISA band and encodes the galactic cosmic-ray spectrum.","key_machinery":"The key machinery is the phenomenological power formula P ≈ ζ G c $m^{2}$ $β^{6}$ $γ^{4}$ / $R^{2}$ with ζ ≈ 5, which matches both the non-relativistic result (scaling as $β^{6}$) and the ultra-relativistic result (scaling as $γ^{4}$), combined with the assumption that each particle radiates all its gravitational power at the single frequency f_GW = $qBc^{2}$/(πE), twice the orbital frequency. This frequency mapping, together with the measured interstellar cosmic-ray spectra (rewritten as functions of frequency), allows a direct convolution to produce the background Ω(f).","core_discovery":"The central claim is that the stochastic gravitational-wave background from Galactic cosmic rays has a computable, non-trivial spectral shape with a maximum around a few millihertz, about 5 mHz for electrons and 2 mHz for protons, arising from the convolution of the single-particle gravitational radiation power with the cosmic-ray number density per unit frequency. The amplitude is extremely small, with power flux in the range $10^{-47}$ J $m^{-2}$ for electrons and $10^{-50}$ J $m^{-2}$ for protons, so it is not a measurable signal; the interest lies in the fact that a well-established physical process produces a background with a distinctive peak in a band that future space-based interferometers will explore.","pith_inferences":["The monochromatic approximation, which places all radiation at twice the orbital frequency, could be replaced by the full spectrum of Eq. (16); if the true spectrum is broad, the peak of Ω(f) might be slightly wider or shifted, so the exact peak frequency is tied to that approximation.","The same formalism could be applied to other galaxies; if their cosmic-ray spectra and magnetic fields differ, the extra-galactic background might have a different effective peak, and the 'typical galaxy' assumption could be relaxed.","The ratio of gravitational to electromagnetic power for a charged particle in a magnetic field could, in principle, be used to search for new ultra-heavy charged particles, since gravitational-wave emission would dominate for masses above roughly 10^-10 kg.","Although undetectable today, this background is a guaranteed foreground for any future detector with sensitivity better than about 10^-47 J m^-2 in the mHz band; computing it removes a potential confusion source for such instruments."],"forward_implications":["The Galactic cosmic-ray background has a peak in the mHz range, overlapping the most sensitive band of LISA, but with an amplitude many orders of magnitude below the instrument's reach.","The shape of Ω(f) is determined by the convolution of the cosmic-ray spectrum with the radiated power, so future precise cosmic-ray measurements could in principle be cross-checked against the predicted gravitational-wave spectrum.","The extra-galactic component, though extremely small, acquires a redshift-smeared cutoff and is essentially isotropic; its shape near the cutoff differs from the Galactic one.","For ultra-heavy charged particles with mass above about 6.8×10^-10 kg, gravitational synchrotron power would exceed electromagnetic synchrotron power at relativistic speeds.","The strain amplitude of a single particle depends only on its speed, not on the magnetic field strength, a counterintuitive feature of the emission."],"supporting_citations":[{"why":"Supplies the rigorous derivation of the ultra-relativistic gravitational radiation power from a charged particle in a magnetic field.","marker":"[13]"},{"why":"Gives the spectrum of gravitational synchrotron radiation, shown to be monotonically decreasing and dominated by the fundamental mode.","marker":"[16]"},{"why":"Provides the textbook quadrupole formula and the non-relativistic expression for the radiated power used in the low-energy limit.","marker":"[11]"},{"why":"Provides the interstellar electron and proton flux parameterizations used to build the volume density per frequency.","marker":"[22]"},{"why":"Supplies the proton parametrization that the paper adjusts to better fit interstellar fluxes in the relevant energy range.","marker":"[23]"},{"why":"Gives the Galactic magnetic field value of about 6 µG used in the calculation.","marker":"[24]"},{"why":"Defines the LISA band that the predicted peak frequency of the background falls into.","marker":"[25]"}],"fun_headline_variants":["Cosmic-ray gravitational wave background peaks near 5 mHz","Gravitational waves from cosmic rays peak at millihertz","Cosmic-ray background: gravitational wave peak at 5 mHz","Weak cosmic-ray gravitational wave background has mHz peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predicted peak frequency and shape of the background depend on the approximate formula for the radiated power that blends the slow and fast limits with a single coefficient, and on the simplification that each particle radiates all its gravitational power at exactly twice its orbital frequency.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic-ray gravitational wave background peaks near 5 mHz","Gravitational waves from cosmic rays peak at millihertz","Cosmic-ray background: gravitational wave peak at 5 mHz","Weak cosmic-ray gravitational wave background has mHz peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000492,"raw_usage":{"total_tokens":2317,"prompt_tokens":741,"completion_tokens":1576,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":357,"completion_tokens_details":{"reasoning_tokens":1507}},"tokens_in":357,"tokens_out":1576,"duration_ms":11497,"temperature":1.0,"reasoning_tokens":1507,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:37:38.457171+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the background using the full gravitational synchrotron spectrum of Eq. (16) instead of the monochromatic approximation, and compare the resulting Ω(f) peak with the claimed 5 mHz for electrons and 2 mHz for protons; if the peak moves by more than a factor of a few, the central claim fails. A second falsifier would be a future detector with sensitivity around $10^{-47}$ J $m^{-2}$ in the mHz band, which could either find the predicted background or rule it out.","supporting_citations":[{"cited_title":"Pustovoit and M","cited_arxiv_id":null,"evidence_quote":"Supplies the rigorous derivation of the ultra-relativistic gravitational radiation power from a charged particle in a magnetic field."},{"cited_title":"Davis, R","cited_arxiv_id":null,"evidence_quote":"Gives the spectrum of gravitational synchrotron radiation, shown to be monotonically decreasing and dominated by the fundamental mode."},{"cited_title":"Maggiore,Gravitational Waves","cited_arxiv_id":null,"evidence_quote":"Provides the textbook quadrupole formula and the non-relativistic expression for the radiated power used in the low-energy limit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the interstellar electron and proton flux parameterizations used to build the volume density per frequency."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the proton parametrization that the paper adjusts to better fit interstellar fluxes in the relevant energy range."},{"cited_title":"Non-parametric determination of H and He interstellar fluxes from cosmic-ray data","cited_arxiv_id":"1511.08650","evidence_quote":"Gives the Galactic magnetic field value of about 6 µG used in the calculation."}],"review_version":1}