{"id":"f26c178c-31d1-49d3-871d-366117415b64","arxiv_id":"1908.05275","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The intact interstellar gas disk of the Milky Way implies charged dark matter must have charge-to-mass ratio below about 10^-13 e/GeV, unless it makes up only a small fraction of dark matter.","lead":"Charged dark matter passing through a spiral galaxy's magnetic field would drag on the galaxy's gas disk and spin it down within a billion years, so the Milky Way's intact disk rules out a huge range of charged dark matter charges and masses. The result targets the same 'millicharged' dark matter that experiments are currently trying to detect.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Assumption 6 (Sec. III.1) posits a pre-formed, stable 5 Gyr ISM disk; qDM near the bound could alter galaxy formation and invalidate this initial condition, so the exclusion is not yet established.","rationale":"The reader's weakest_assumption correctly identifies the formation-timeline premise as the most load-bearing point. My independent reading of Sec. III.1 Assumption 6, Fig. 2, and the Conclusion confirms that the paper's initial condition (a present-day-like rotating disk and ordered magnetic field existing for T_ism ~ 5 Gyr before qDM spin-down) is an assumption rather than a derived result. The same Lorentz coupling responsible for spin-down would act during galaxy formation, so the bound is conditional on qDM not spoiling standard galaxy formation. This is a real concern, but it is explicitly acknowledged by the authors and does not invalidate the analytic derivation of the torque. A conditional verdict with high confidence is therefore appropriate; I would not change the reader's verdict. I also note an internal numerical inconsistency not central to my attack: Eq. (59) implies a central value q/mc^2 < 3.5e-14 e/GeV, while Eq. (65) and the abstract quote 3.5e-13 and 1e-13 e/GeV, respectively. This factor-of-ten discrepancy sits inside the paper's own ad hoc two-order-of-magnitude uncertainty band, so it does not overturn the conclusion, but it should be corrected in any revision.","tokens_in":25814,"tokens_out":9451,"duration_ms":100769,"concrete_test":"Run a cosmological zoom-in magnetohydrodynamic simulation of a Milky Way-mass halo with a charged dark matter component at f_qdm = 1 and q/mc^2 near 3.5e-14 e/GeV, including the qDM Lorentz-force back-reaction on the gas. Check whether the z = 0 disk radius, gas surface density, and ordered field strength match observed spirals; if no normal disk forms, the bound fails, and if a normal disk forms and then shrinks, the bound is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central bound is derived by applying present-day local ISM parameters (Sigma_ism = 17.2 M_sun/pc^2, Phi_eff ~ 1e18 V, rho_dm = 0.008 M_sun/pc^3, T_ism = 5 Gyr) to a disk that is assumed to have existed and remained stable for ~5 Gyr before qDM spin-down acted. This is Assumption 6 in Sec. III.1 and the timeline in Fig. 2. The load-bearing step is the implicit initial condition: a rotation-supported, magnetized ISM disk essentially like today's existed before the qDM torque became important. But the same qDM-B coupling that drives spin-down at q/m ~ 3.5e-14 e/GeV would also have operated during protogalactic collapse, disk assembly, and dynamo growth. If a near-bound qDM population removes angular momentum from infalling gas, prevents formation of a large rotation-supported disk, or suppresses growth of the ordered microgauss field, then a normal spiral galaxy would never develop the initial state from which the spin-down argument starts. In that case the observation that present-day disks are not shrunken would not exclude the qDM parameters; the bound would not follow. The authors explicitly acknowledge this possibility in Sec. III.1 ('It is interesting to explore whether qDM interactions with magnetic fields would actually modify galaxy formation... beyond the scope of this paper') and again in the Conclusion ('we have assumed that the presence of a qDM population at early times does not spoil our conventional understanding of galaxy formation'). Because this assumption is both load-bearing and untested, the exclusion claim runs ahead of the analysis.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an upper limit on the charge-to-mass ratio of a charged dark matter (qDM) population by considering angular momentum transfer between a rotating, magnetized interstellar medium (ISM) disk and a slowly rotating qDM halo. Working in the ballistic regime and using a second-Born perturbative expansion of the Lorentz force, the authors obtain a spin-down timescale for the ISM and show that, for local Milky Way parameters, avoiding order-unity ISM spin-down over T_ism ~ 5 Gyr requires q/(m c^2) below roughly 1e-14 e/GeV for f_qdm = 1, with an ad hoc two-order-of-magnitude uncertainty. The paper also discusses qDM spin-up for small f_qdm, compares the resulting exclusion with freeze-in production and other millicharged-particle bounds, and sketches possible applications to galaxy clusters and the Bullet Cluster.","tokens_in":26214,"tokens_out":8136,"duration_ms":79852,"significance":"If the result is correct, this is a novel and potentially very strong astrophysical bound: it is roughly two orders of magnitude stronger than the simple Larmor-radius coherence argument and would exclude much of the freeze-in millicharged parameter space below about a TeV. The analytic machinery for momentum transfer in the second Born approximation is original and carefully laid out, and the inputs are independently measured quantities rather than quantities fitted to the claimed bound. The paper is also unusually candid: it explicitly lists its assumptions, acknowledges the rough nature of the result, and admits that no observational spin-down test is specified. These strengths make the derivation worth publishing in some form, but the headline claim needs substantial revision before it can stand as a quantitative exclusion.","major_comments":[{"comment":"The central value of the headline bound is internally inconsistent. Eq. (58) with tau_sd = 5 Gyr gives q/(m c^2) ~ 3.5e-14 e/GeV for f_qdm = 1, and Eq. (59) states this 3.5e-14 value. The Conclusion, however, restates the bound as 3.5e-13+/-1 e/GeV (Eq. (65)), and the Abstract says q/e ~ 1e-13+/-1 m c^2/GeV. These differ from Eq. (59) by a factor of ten in the central value. If Eqs. (58)-(59) are correct, the Abstract and Eq. (65) must be corrected; if the weaker value is intended, the derivation leading to Eq. (58) must be revisited. As written, the reader cannot tell which number is the claimed limit.","section":"Eqs. (59), (65), and Abstract"},{"comment":"The bound rests on the timeline in which the ISM disk and its ordered magnetic field already exist in roughly their current form for T_ism ~ 5 Gyr before qDM spin-down acts (Assumption 6 and Fig. 2). But a qDM population with q/m near the bound would have interacted with magnetic fields during protogalactic collapse, disk assembly, and dynamo growth. Those interactions could remove angular momentum from infalling gas, prevent the formation of a large rotation-supported disk, or suppress the ordered microgauss field, so the assumed initial state would never have formed. In that case the observation that present-day disks are not shrunken would not exclude the qDM parameters. The authors acknowledge this in Sec. III.1 and again in the Conclusion, but stating that the question is beyond the scope of the paper leaves the central exclusion unestablished. I request at least a quantitative estimate of the early-time effect, or an explicit reframing of the result as conditional on conventional galaxy formation, with the word 'limit' qualified accordingly.","section":"Sec. III.1, Assumption 6; Fig. 2; Conclusion"},{"comment":"The logarithmic enhancement that enters the spin-down rate comes from the heuristic saturation regulator in Eq. (43), which replaces the divergent <1/|v_z|> average. The authors state that this cut-off could be implemented more rigorously but do not do so. Because this regulator controls the factor 1 - gamma_E - 2 ln(epsilon) in Eq. (49) and hence the magnitude of the dominant spin-down term, its systematic uncertainty propagates directly into the central bound. The ad hoc two-order-of-magnitude uncertainty is broad, but the derivation should still justify the regulator quantitatively or demonstrate explicitly that the bound is insensitive to the regulator's form.","section":"Sec. III.B.2, Eq. (43)"},{"comment":"The exclusion regions in Figs. 3 and 4, and the phrase 'almost certainly ruled out' in the Fig. 3 caption, presuppose that a spun-down ISM would have been noticed, but the paper does not define a specific observable or a statistical criterion. Without an observational spin-down test, Eq. (59) is a theoretical estimate of when spin-down becomes large, not a measurement-driven upper limit. The authors are candid about this in Sec. VII, but the central claim of a new exclusion should either propose a concrete test or clearly label the result as an estimate pending such a test.","section":"Sec. IV.4; Sec. VII"}],"minor_comments":[{"comment":"The sentence 'we set the criteria for no significant spin-down to be tau_sd < T_ism' has the inequality backwards; it should be tau_sd > T_ism.","section":"Sec. II, after Eq. (4)"},{"comment":"The phrase '< vec v / |hat v| >' appears to be a typo; the non-divergent average is presumably < vec v / |v_z| > or similar.","section":"Sec. III.B.2, after Eq. (43)"},{"comment":"The superscript (1) on Delta vec p_{4+/-} in the momentum-transfer-rate equation is inconsistent with the second-Born definition in Eq. (22); the superscript should be (2).","section":"Sec. III.A.5, Eq. (32)"},{"comment":"The statement that the pitch angle is 'hat v independent' should be phrased as psi(hat v) -> psi in the planar limit; as written it is confusing.","section":"Sec. III.B.1, Eq. (40)"},{"comment":"References [32] and [50] appear to be the same paper with the same arXiv number and should be consolidated or cross-referenced.","section":"References [32] and [50]"},{"comment":"Once the central-value inconsistency is fixed, the notation for the uncertainty exponent should be made uniform: Eq. (59) uses -14 with 'minus-plus 1' while Eq. (65) uses -13 with 'plus-minus 1'.","section":"Eqs. (59) and (65)"}],"recommendation":"major_revision","confidential_remarks":"This paper is borderline between major revision and reject. The underlying derivation is careful and the authors are unusually honest about its limitations, so I would not reject outright. However, the headline bound is stated inconsistently by a factor of ten, and the formation-backreaction issue (Assumption 6) undermines the claim that this is an established exclusion. The authors should be asked to correct the numerical inconsistency and either quantify the early-universe/assembly-time effect or explicitly reframe the result as a conditional estimate. If they can do that, the paper would be a valuable contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it is a real contribution: the second-Born momentum-transfer calculation from ordered galactic magnetic fields is new, and the resulting bound on sqrt(f_qdm) q/m of about 1e-14 e/GeV lands squarely on freeze-in millicharge parameter space that is being actively pursued. Second, the central exclusion is not as solid as the abstract suggests, because the whole argument assumes a pre-existing, stable ISM disk and ordered field for ~5 Gyr, and the same qDM-B coupling could have altered galaxy formation and dynamo growth in the first place. The authors know this; they flag it explicitly as Assumption 6 and in the conclusion, but they do not resolve it.\n\nThe derivation itself is careful and mostly persuasive. I traced the logic from the Maxwell-Vlasov setup through the four-particle average, and the cancellation of first-Born terms with survival at second order is clean. The treatment of realistic field geometries via RM maps is a nice touch, and the authors are commendably honest about the lack of an observational spin-down test and the ad hoc two-order-of-magnitude uncertainty. The bound is a genuine first-principles estimate, not a fit.\n\nThe soft spots are real but not fatal. The formation-timeline assumption is the load-bearing one, and the stress-test concern is on target: if near-bound qDM prevents a normal disk from forming, then observing an unsHrunken disk does not exclude anything. That assumption is untested, and the paper's own admission that exploring it is beyond scope is not a substitute for doing so. The saturation regulator in Eq. (43) is heuristic, but the logarithmic divergence means the result is not sensitive to the exact choice, so I regard that as minor. The abstract's wording that the bound 'excludes' freeze-in production runs slightly ahead of the text's careful caveats, but this is a presentation issue, not a mathematical error.\n\nWho is this for? Anyone working on millicharged dark matter, dark photons, or astrophysical signatures of charged relics. It is a useful paper to put on a reading group syllabus, and it deserves a serious referee. My recommendation: send it to review, but the referee should demand a clearer statement that the bound is conditional on standard galaxy formation, or better yet a short simulation section showing that qDM near the bound does not suppress disk assembly. With that revision, the paper would be a solid contribution rather than an overclaimed exclusion.","headline":"A genuinely new and strong bound on charged dark matter from galactic ISM spin-down, honestly caveated, but the exclusion claim rests on a load-bearing assumption about galaxy formation that the paper does not test.","tokens_in":26711,"tokens_out":1917,"would_cite":true,"duration_ms":21458,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that charged dark matter crossing a spiral galaxy's magnetized gas disk would extract most of the disk's angular momentum within about 5 Gyr unless the charge-to-mass ratio is below roughly $3.5\\times10^{-13}$ e/GeV (for…","keywords":["charged dark matter","millicharged particles","galactic magnetic fields","interstellar medium","angular momentum transfer","spin-down","freeze-in production","dark matter constraints"],"falsifier":"Run a magnetohydrodynamic simulation of a Milky-Way-like galaxy from redshift $z\\approx0.4$ to the present with $f_{\\rm qdm}=1$ and $q/(m c^2)=3\\times10^{-13}\\,e/{\\rm GeV}$; if the simulated gas disk retains more than half of its initial angular momentum while the ordered field stays near a microgauss, the paper's central bound would be refuted.","tokens_in":25582,"feed_emoji":"🌀","tokens_out":10265,"duration_ms":95837,"temperature":0.7,"pith_summary":"Most of the dark matter could be neutral, but a small electrically charged component remains viable. This paper argues that such charged dark matter cannot merely orbit a spiral galaxy: each time it crosses the galaxy's gas disk, the disk's embedded magnetic field deflects it, and the reaction drags on the disk's rotation. Over the $\\sim5$ Gyr that the Milky Way's gas disk and its $\\sim\\mu$G ordered field are thought to have persisted, this drag would remove an order-unity fraction of the disk's angular momentum unless the charge-to-mass ratio obeys $\\sqrt{f_{\\rm qdm}}\\, q/(m c^2) \\lesssim 3.5\\times10^{-14}\\, e/{\\rm GeV}$, i.e. $q/(m c^2)\\lesssim 3.5\\times10^{-13}\\, e/{\\rm GeV}$ if all the dark matter is charged. This bound is strong enough to exclude freeze-in production of charged dark matter below about a TeV, making large-scale magnetic fields a new probe of dark matter's non-gravitational interactions.","feed_headline":"Charged dark matter would spin down galaxy disks within 5 Gyr","feed_subtitle":"Milky Way's microgauss field forces charge-to-mass below 3.5e-13 e/GeV for fully charged dark matter.","key_machinery":"The load-bearing object is the line integral of the magnetic field along a qDM trajectory, $\\Phi = \\int \\vec B\\cdot d\\vec\\ell$, a quantity with dimensions of voltage; near the Sun the paper estimates $\\Phi_{\\rm eff}\\approx 3.5\\times10^{17}$ V. The argument runs through the second Born approximation: first-order kicks cancel between opposite charges and opposite directions, and the surviving momentum transfer is proportional to $(q\\Phi/(m c \\sigma))^2$, so the spin-down rate is set by the square of the ratio between the magnetic 'voltage' and the particle's charge-to-mass voltage $m c^2/q$. A four-particle reflection average isolates this surviving second-order term, and a logarithmic cutoff regulates low-inclination orbits that would otherwise be over-deflected. The machinery turns a purely electromagnetic interaction into an angular-momentum drag whose rate depends on $q/m$, the field voltage, the ISM surface density, and the qDM density and velocity dispersion.","core_discovery":"The paper aims to establish that the ordered component of a spiral galaxy's magnetic field acts as a brake on the interstellar medium when charged dark matter passes through it. Because the field is embedded in the highly conducting, co-rotating ISM, the ideal MHD electric field $\\vec E \\approx -\\vec V/c \\times \\vec B$ makes each qDM particle feel a force proportional to $(\\vec v-\\vec V)\\times\\vec B$. For any incoming velocity, opposite charges receive the same kick at leading order, so the first Born contribution to the momentum transfer cancels after averaging over charges and directions; the surviving effect is second order in the field. Averaging over four particles (opposite charges with velocities $\\pm\\vec v$) reduces the answer to two line integrals of $\\vec B$ along the trajectory, and for planar fields the spin-down rate becomes proportional to $(q \\Phi_{\\rm eff} \\sin\\psi/m c \\sigma)^2$, where $\\Phi_{\\rm eff}$ is the line integral of the ordered field, about $3.5\\times10^{17}$ V in the solar neighborhood. Requiring the spin-down time to exceed $T_{\\rm ism}=5$ Gyr yields the limit $\\sqrt{f_{\\rm qdm}}\\, q/(m c^2) \\lesssim 3.5\\times10^{-14}\\, e/{\\rm GeV}$ for $f_{\\rm qdm}\\gtrsim 0.13$.","pith_inferences":["A straightforward observational follow-up the paper does not carry out is to compare the radial sizes of gas disks versus stellar disks across a large sample of spiral galaxies; a systematic excess of small gas disks in galaxies with strong ordered fields would corroborate the spin-down picture.","The same second-order drag mechanism should operate in any magnetized plasma with a relative velocity to a charged halo, so merging galaxy clusters, whose fields can exceed 10 $\\mu$G over tens of kpc, could set even tighter limits; the paper sketches this for the Bullet Cluster but leaves the calculation schematic.","The calculation regulates grazing trajectories with a cutoff rather than evolving the qDM distribution self-consistently; a particle-in-cell or MHD simulation would show whether the drag saturates as slow particles are captured into co-rotation, softening the sharpest part of the excluded region."],"forward_implications":["If all the dark matter is charged, the bound excludes the freeze-in production scenario for masses below roughly a TeV, since that mechanism requires $q/e\\sim10^{-11}$, far above the limit.","For charged-dark-matter fractions below $f_{\\rm qdm}\\approx0.13$ there is not enough inertia in the dark matter to drag the disk; the prediction flips to qDM being spun up into co-rotation, altering its local velocity distribution rather than shrinking the gas disk.","The bound applies not only to directly charged particles but also to neutral dark matter coupled through an ultralight kinetically mixed dark photon, because the ISM plasma gives the photon an effective mass and the dark matter an effective charge.","At unit charge, the constraint reaches masses near $10^{13}$ GeV, cutting into the 'WIMPzilla' regime of supermassive dark matter."],"supporting_citations":[{"why":"Supplies the local dark matter density used to set the qDM flux and mass fraction.","marker":"[13]"},{"why":"Provides the ISM surface density that sets the inertia being spun down.","marker":"[14]"},{"why":"Gives the ordered magnetic field strength and scale height near the Sun used for the effective voltage.","marker":"[15]"},{"why":"MHD simulations motivating the 5 Gyr timescale over which the field is assumed stable.","marker":"[16]"},{"why":"Rotation-measure maps used to infer the local coherent field direction and pitch angle.","marker":"[28]"},{"why":"Defines the freeze-in production cross section whose predicted abundance is excluded below a TeV.","marker":"[11]"},{"why":"Computes freeze-in abundances including plasmon effects, forming the comparison curve.","marker":"[12]"}],"fun_headline_variants":["Galaxy spin-down sets sharpest limit on charged dark matter","Charged dark matter would brake galaxies: new bound on charge","Magnetic fields expose charged dark matter: spin-down limit","Galactic magnetic fields rule out most charged dark matter","Charged dark matter would drain galaxy spin within 5 Gyr"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a rotating gas disk with roughly today's magnetic field already existed for about five billion years before charged dark matter began dragging on it; if charged dark matter would have disrupted galaxy formation or magnetic-field growth earlier, the initial condition from which the spin-down bound is computed never holds.","fun_headline_variants_meta":{"raw":{"variants":["Galaxy spin-down sets sharpest limit on charged dark matter","Charged dark matter would brake galaxies: new bound on charge","Magnetic fields expose charged dark matter: spin-down limit","Galactic magnetic fields rule out most charged dark matter","Charged dark matter would drain galaxy spin within 5 Gyr"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000875,"raw_usage":{"total_tokens":3884,"prompt_tokens":1142,"completion_tokens":2742,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":758,"completion_tokens_details":{"reasoning_tokens":2656}},"tokens_in":758,"tokens_out":2742,"duration_ms":17559,"temperature":1.0,"reasoning_tokens":2656,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:18:38.040291+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a magnetohydrodynamic simulation of a Milky-Way-like galaxy from redshift $z\\approx0.4$ to the present with $f_{\\rm qdm}=1$ and $q/(m c^2)=3\\times10^{-13}\\,e/{\\rm GeV}$; if the simulated gas disk retains more than half of its initial angular momentum while the ordered field stays near a microgauss, the paper's central bound would be refuted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the local dark matter density used to set the qDM flux and mass fraction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the ISM surface density that sets the inertia being spun down."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the ordered magnetic field strength and scale height near the Sun used for the effective voltage."},{"cited_title":"We assume that initially f+ = f− = f0, but f will evolve away from this as the qDM particles travel through the magnetic field","cited_arxiv_id":null,"evidence_quote":"MHD simulations motivating the 5 Gyr timescale over which the field is assumed stable."},{"cited_title":"Here we suppose that these evolutionary effects can be modeled as if the current configuration has been constant over an interval of duration Tism","cited_arxiv_id":null,"evidence_quote":"Rotation-measure maps used to infer the local coherent field direction and pitch angle."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the freeze-in production cross section whose predicted abundance is excluded below a TeV."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Computes freeze-in abundances including plasmon effects, forming the comparison curve."}],"review_version":1}