{"id":"be55fab4-544e-46a4-8ce1-c2e44910629a","arxiv_id":"1908.04675","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulated AQN dark matter crossing Earth produces relativistic axions at around 10^14 eV/cm2/s with about 10% daily and annual modulations and rare 100 to 10,000 times local flash amplifications.","lead":"This paper simulates dark-matter nuggets crashing through the Earth and calculates the axion stream they produce, including daily wobbles and rare bright bursts. It gives axion detectors a concrete new search pattern if this particular dark-matter model is right.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline flux Eq. (2) is normalized by an undefined parameter ε=1 in Table VI; if ε is an efficiency or antimatter fraction below 1, the central numerical predictions scale down linearly.","rationale":"The reader's weakest assumption correctly identifies the geometric cross-section and the 1/3 axion energy fraction as the main physical premises behind the flux normalization. My concern is narrower and more internal: the results table is computed with an explicit parameter ε=1 that is never defined, and the paper's own discussion of antimatter AQNs suggests ε should not be unity. This is the most load-bearing spot because Eq. (2) is the headline claim, and every modulation amplitude, statistical fluctuation, and local-flash rate inherits this normalization. The concern does not overturn the paper's qualitative picture: AQN-induced relativistic axions with strong time modulations remain a concrete, falsifiable prediction. It does mean that the absolute flux scale, and hence the detection prospects, are not yet secured. The recommended verdict stays CONDITIONAL, matching the reader's judgment, because the issue is a missing definition and a plausible rescaling rather than a demonstrated contradiction. I chose 'partial' agreement because the reader's weakest assumption is the same normalization chain, but the reader did not flag the specific undefined ε=1 in Table VI, which is the cleanest way to expose the fragility.","tokens_in":30982,"tokens_out":36558,"duration_ms":382044,"concrete_test":"Locate the definition of ε in Ref. [2], Ref. [47], or the simulation code used for Table VI. If ε denotes an efficiency or the antimatter fraction, recompute Eq. (38) and Table VI with the justified value, e.g., ε = 2/5 for the antimatter fraction or ε = 1/3 for the domain-wall energy fraction. If the resulting flux falls below ~3×10^13 eV cm^-2 s^-1, the headline normalization of Eq. (2) is not established; if the flux remains within a factor of two of 10^14, the central claim is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim, Eq. (2), is a product of the AQN hit rate, the geometric-annihilation mass loss, and the axion conversion efficiency. Table VI, which supplies the absolute flux values, is computed with the caption condition 'ε = 1 unless specified,' but ε is never defined anywhere in the paper. If ε denotes the annihilation efficiency, the axion emission efficiency, or the fraction of AQNs that are antimatter, then every entry in Tables III, IV, and VI, and the normalization of Eq. (2), scale linearly with ε. This is not a minor bookkeeping point: Sec. VI states that only antimatter AQNs annihilate underground, while Eq. (5) and Eq. (35) quote a total AQN hit rate, suggesting that a suppression factor on the order of 0.4 (or smaller) may be missing. Likewise, Eq. (36) already inserts a factor 1/3 for the domain-wall energy, so ε=1 is not that factor. The paper's claim in Sec. VII that the predictions are insensitive to model parameters is unsupported as long as ε is not defined and varied. If the true ε is, say, 0.4, the headline flux drops from ≈10^14 to ≈4×10^13 eV cm^-2 s^-1; if ε is a small efficiency, the drop is larger. Without a definition of ε, the reader cannot determine whether Eq. (2) is an upper limit or a central prediction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies axion production from Axion Quark Nuggets (AQNs) crossing the Earth. After reviewing the AQN model and the domain-wall axion emission mechanism, the authors perform Monte Carlo simulations of AQN trajectories through a five-layer Earth, compute the heat/axion-emission profile and the surface axion flux, and quantify annual and daily modulations, Poisson statistical fluctuations, and rare 'local flash' burst amplifications. The central result, Eq. (2), is an energy flux of approximately 10^14 A(t) eV cm^-2 s^-1 with average axion energy <E_a> ≈ 1.3 m_a and average velocity <v_a> ≈ 0.6c; the time-dependent factor includes annual and daily modulations of order 1-10%, statistical fluctuations of 20-60%, and local-flash amplifications of 10^2-10^4. The paper also argues that gravitational-lensing amplification is negligible under the standard halo model.","tokens_in":31302,"tokens_out":10382,"duration_ms":116460,"significance":"If the normalization and model inputs hold, the paper gives a concrete, falsifiable prediction for a new relativistic-axion channel (v_a ≈ 0.6c) whose time structure is unique to the AQN dark matter model. The simulation machinery is described in enough detail to reproduce the heat-emission profile q(r,θ) and the surface angular distribution P_a(θ), and the internal consistency between Eq. (38) and Tables III and VI is a strength. The work also usefully identifies the daily modulation and local-flash effects as observables that distinguish AQN-induced axions from conventional galactic axions and WIMPs. However, the absolute flux and all derived amplitudes inherit several model inputs from earlier work, and the paper's robustness claims go beyond what is actually varied.","major_comments":[{"comment":"The parameter ε in the Table VI caption ('ε = 1 unless specified') is never defined anywhere in the text. This is load-bearing because every flux and density entry in Table VI, and therefore the normalization of Eq. (2), scales linearly with ε. Moreover, Sec. VI B states that 'only the AQNs made out of antiquarks will be annihilated underground' and mentions a numerical factor of 3/5, but Eq. (38) contains no such factor if Eq. (35) is the total AQN hit rate. The authors must define ε, specify whether Eq. (35) is the total or antiquark-only hit rate, and if necessary insert the antimatter fraction in Eq. (38) and rescale Tables III, IV, and VI. Without this, Eq. (2) cannot be interpreted as a central prediction rather than an upper limit.","section":"Table VI, Eq. (2), Eq. (38)"},{"comment":"The claim that the predictions have 'little to no flexibility' is stronger than the tables show. Table VI spans a factor of roughly 2.8 in the flux (from 3.46×10^13 to 9.67×10^13 eV cm^-2 s^-1), and the daily-modulation amplitude also varies substantially across the models, with <ΔB>/<B> ranging from 14.1% to 33.8%. The sensitivity study in Appendix D varies α and B_min, but it does not vary ε, the 1/3 energy fraction in Eq. (36), or the spectral parameter δ, which are precisely the quantities that set the absolute amplitude. The conclusions should report a range of predictions rather than claiming near-parameter-independence.","section":"Sec. VII and Appendix D"},{"comment":"The conversion of mass loss to axion number uses the factor 1/3 for the domain-wall energy fraction and ⟨E_a⟩ ≃ 1.3 m_a. These are model inputs inherited from Refs. [1, 47], not outputs of the present simulation. The paper should state explicitly that the energy flux and all modulation amplitudes scale linearly with the product of the annihilation efficiency and this energy fraction, and should comment on the theoretical uncertainty in these inputs. This is especially important because the manuscript presents the 1/3 fraction as fixed while later claiming that the results are insensitive to model parameters.","section":"Sec. IV C, Eq. (36)"}],"minor_comments":[{"comment":"The first sentence of Sec. VI B contains a typo: 'obtianed' should be 'obtained'.","section":"Sec. VI B"},{"comment":"The opening sentence of Appendix D contains a typo: 'senstive' should be 'sensitive', and the phrase 'argue that thee main resultss' has repeated letters.","section":"Appendix D"},{"comment":"The caption says 'ϵ = 1 unless specified', but no row in Table VI specifies a different value; if ε is not a physical parameter, the reference to it should be removed or explained.","section":"Table VI caption"},{"comment":"The text and figures use both 63 degrees (Fig. 2) and 60 degrees (Fig. 4 and Sec. III E) for the angle between the DM wind and the ecliptic or celestial equator; the geometry should be stated consistently.","section":"Fig. 2 and Sec. III E"},{"comment":"References [66] and [67] appear to be the same Abramowitz entry duplicated; the bibliography should be checked.","section":"References [66] and [67]"}],"recommendation":"major_revision","confidential_remarks":"The core Monte Carlo computation is plausible and the paper makes a striking, testable prediction, but the undefined ε and the unclear antimatter-fraction treatment are central normalization issues that must be fixed before the numerical claims can be relied upon. I would be willing to review a revised version that addresses these points and softens the parameter-insensitivity claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a model-prediction paper, not a data fit. It computes the axion flux from AQN dark matter passing through the Earth, including annual/daily modulation, Poisson fluctuations, and rare 'local flashes'. The daily modulation mechanism—the AQN cross-section shrinks as it loses mass, so it emits more axions on the wind-facing side of Earth—is new and physically clear. The local-flash estimates (amplification up to 10^2–10^4 for a nearby passage, with event rates in Table IV) are the most interesting output and are concrete enough for broadband axion detectors to take seriously. The gravitational lensing section is a useful service: it shows the earlier large-lensing claims require implausible colinear DM streams.\n\nThe soft spots are real but not fatal. The stress-test note is right: Table VI is computed with ε = 1 and ε is never defined anywhere. If ε is an efficiency or the antimatter fraction, the headline flux Eq. (2) scales linearly down, and the Sec. VII claim that the predictions have 'little to no flexibility' overstates the case. Second, the text says ΔB/B ~ 0.1 in most cases, but Table VI lists 14–34%. One of those is a typo or a misreading, and it matters because the modulation amplitudes are normalized by that quantity. Third, the absolute flux rests on the geometric cross-section and the 1/3 axion energy fraction, both inherited from the authors' earlier work. That is not a fatal circularity—the predictions are new and falsifiable—but the normalization is less robust than the modulation shapes.\n\nWhat the paper does well: the Monte Carlo is internally consistent, the simulated heat profile and Pa(θ) match the analytical daily-modulation estimate (14), the parameter sensitivity across the six model variants is modest (a factor of ~3 in flux, which they understate as 'insensitive'), and the local-flash rate formula (22) is a clean scaling relation. No code or data are shipped, so the 10^8-particle simulation can't be audited, but the equations are transparent enough to reimplement.\n\nWho it's for: experimentalists planning broadband axion searches, and DM model-builders who want concrete AQN signatures. Even a skeptic of AQNs can extract the time-modulation formalism and the lensing correction.\n\nRecommendation: accept for peer review. The required fixes are small: define ε, fix the ΔB/B mismatch, and release the code. The core predictions are testable and the paper is honest about its assumptions, which puts it ahead of many model-prediction papers.","headline":"A genuine AQN model-prediction paper with concrete, testable signatures, but the key table has an undefined ε and the text disagrees with Table VI on ΔB/B; worth refereeing after small fixes.","tokens_in":31863,"tokens_out":5978,"would_cite":false,"duration_ms":61214,"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":"This paper claims that Axion Quark Nugget dark matter produces a relativistic axion flux at Earth with annual and daily modulations and rare burst amplifications large enough for broadband detectors to test.","keywords":["axion quark nuggets","dark matter","relativistic axions","annual modulation","daily modulation","local flashes","broadband axion detection","Monte Carlo simulation"],"falsifier":"A year-long broadband axion search in the $10^{-6}$--$10^{-3}$ eV window sensitive to the predicted $10^{14}$ eV/(cm$^2$ s) flux should see a 10% daily modulation and occasional local flashes; observing no daily modulation while the annual modulation is present, or no 0.3-second flash at the predicted rate of roughly one per five years at amplification $10^3$, would rule out Eq. (2).","tokens_in":30679,"feed_emoji":"📡","tokens_out":8252,"duration_ms":85470,"temperature":0.7,"pith_summary":"The paper claims that the Axion Quark Nugget (AQN) dark matter model, in which macroscopic quark-matter nuggets are stabilized by axion domain walls, generates a flux of relativistic axions when antimatter AQNs annihilate inside Earth. Using Monte Carlo simulations of AQN trajectories through a layered Earth, it predicts a time-dependent energy flux $\\langle E_a\\rangle\\Phi_a(t)\\simeq 10^{14} A(t)$ eV/(cm$^2$ s), where $A(t)$ contains an annual modulation of 1--10%, a daily modulation of order 10%, Poisson statistical fluctuations of 20--60%, and rare local-flash amplifications of $10^2$--$10^4$. The axions have mean speed $\\langle v_a\\rangle\\simeq 0.6c$, which distinguishes them from conventional non-relativistic galactic axions and gives broadband detectors a concrete temporal search strategy. The predicted flux and modulations are stated to be largely insensitive to the model's free parameters, so the time structure is a testable signature of the AQN framework.","feed_headline":"Dark matter nuggets fire relativistic axion flashes at Earth","feed_subtitle":"Simulation predicts 1–10% daily modulations, Poisson noise, and bursts 10,000 times the average flux.","key_machinery":"The mechanism is the axion domain-wall relaxation: when an antimatter AQN annihilates baryonic matter, it loses mass and shrinks, shifting the equilibrium domain-wall configuration and exciting a mode that radiates propagating axions. The quantitative machinery is the mass-loss relation $dm/ds=-\\sigma\\rho$ with geometric cross-section $\\sigma\\simeq\\pi R^2$, combined with the axion number estimate $\\langle N_a\\rangle\\simeq (1/3)\\Delta m c^2/\\langle E_a\\rangle\\simeq \\Delta m c^2/(4m_a)$, and the emission spectrum from the companion calculation giving $\\langle v_a\\rangle\\simeq 0.6c$. These ingredients are integrated along Monte Carlo-generated trajectories through Earth's five-layer density profile, producing the heat-emission profile $q(r,\\theta)$ whose asymmetry generates the daily modulation and whose local concentration generates the local-flash amplifications.","core_discovery":"The central result is Eq. (2): for the AQN model, the time-dependent axion energy flux at Earth's surface is $\\langle E_a\\rangle\\Phi_a(t)\\simeq 10^{14} A(t)$ eV/(cm$^2$ s), with $\\langle E_a\\rangle\\simeq 1.3 m_a c^2$ and $\\langle A(t)\\rangle = 1$ when averaged over long times. The modulation factor $A(t)$ is composed of an annual modulation $\\kappa_a$ of order 1--10%, a daily modulation $\\kappa_d$ of order 10%, Poisson-driven statistical fluctuations of 20--60%, and rare local-flash bursts that amplify the signal by $10^2$--$10^4$ for short durations. The paper derives these numbers by simulating AQN trajectories through a five-layer model of Earth, using a geometric annihilation cross-section and the axion-emission spectrum from the domain-wall relaxation mechanism. It also states that the resulting energy flux and density are independent of axion mass in the window $10^{-6}\\,\\mathrm{eV}\\lesssim m_a\\lesssim 10^{-3}\\,\\mathrm{eV}$, unlike conventional galactic axions.","pith_inferences":["A null search at the predicted baseline flux would compress all three time signatures proportionally, so a single broadband run could bound the product of the annihilation fraction and the axion energy fraction.","The local-flash event-rate scaling $\\propto A^{-3/2}$ is a parameter-free relation; after several flashes, comparing amplitudes with durations would test the geometric-cross-section assumption independently of the absolute flux.","The same simulation machinery could be applied to the neutrino channel produced by the same annihilation events, making the long-standing annual-modulation signal seen by underground detectors a quantitatively testable prediction rather than an independent anomaly.","Two detectors at similar latitude but separated in longitude could separate the globally coherent daily modulation from local, uncorrelated flashes, since the daily component rotates with Earth while bursts do not repeat."],"forward_implications":["Detectors whose observable couples to the axion gradient rather than the axion density would see AQN-induced axions enhanced by roughly $10^3$ relative to cold galactic axions because the signal scales with axion velocity.","The predicted daily modulation of order 10% and the associated north--south spatial asymmetry provide a timing signature that can separate AQN axions from conventional halo axions and from most instrumental backgrounds.","Local flashes give a rare-event test: an amplification of $10^2$ should appear roughly once every two days as a signal lasting about one second, whereas an amplification of $10^4$ is a 0.1-second burst expected about once every five years.","The baseline flux and energy density are independent of axion mass across the stated window, so the prediction does not require fine-tuning of $m_a$.","Gravitational lensing by the Sun and planets does not amplify the AQN-induced axion flux under the standard halo model, so searches should target the modulation and burst signatures rather than lensing enhancements."],"supporting_citations":[{"why":"Introduces the axion-production mechanism in AQNs and the order-of-magnitude density estimate that this paper refines.","marker":"[1]"},{"why":"Supplies the AQN hit rate, baryon-charge distribution, and energy-loss equations that set the simulation input.","marker":"[2]"},{"why":"Provides the axion-emission velocity spectrum and its normalization, giving $\\langle v_a\\rangle\\simeq 0.6c$.","marker":"[47]"},{"why":"Defines the broadband detection strategy whose annual, daily, and local-flash observables this paper computes.","marker":"[48]"},{"why":"Provides the standard-halo annual-modulation framework used for comparison and for the 1--10% amplitude estimate.","marker":"[53, 54]"},{"why":"Supplies the energy-loss formula $dE/ds=-\\sigma\\rho v^2$ used to determine AQN trajectories through Earth.","marker":"[58]"},{"why":"Provides the layered Earth density model used in the trajectory integration.","marker":"[59]"},{"why":"Fits the AQN baryon-charge distribution parameters from solar extreme-ultraviolet observations, defining the model variants.","marker":"[40]"}],"fun_headline_variants":["Axion quark nuggets amplify signals by up to 10,000 times","Relativistic axion bursts from dark matter nuggets","Daily and annual modulations plus bursts in axion flux","AQN model: axion signal bursts 100-10,000 times average"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole flux normalization rests on the assumption that every nucleus in the AQN's geometric path annihilates and that about one third of the liberated mass energy is emitted as axions; if either fraction is smaller, the flux, daily modulation, and local-flash rates all shrink in proportion.","fun_headline_variants_meta":{"raw":{"variants":["Axion quark nuggets amplify signals by up to 10,000 times","Relativistic axion bursts from dark matter nuggets","Daily and annual modulations plus bursts in axion flux","AQN model: axion signal bursts 100-10,000 times average"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000297,"raw_usage":{"total_tokens":1756,"prompt_tokens":1013,"completion_tokens":743,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":629,"completion_tokens_details":{"reasoning_tokens":667}},"tokens_in":629,"tokens_out":743,"duration_ms":7334,"temperature":1.0,"reasoning_tokens":667,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:10:10.826220+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A year-long broadband axion search in the $10^{-6}$--$10^{-3}$ eV window sensitive to the predicted $10^{14}$ eV/(cm$^2$ s) flux should see a 10% daily modulation and occasional local flashes; observing no daily modulation while the annual modulation is present, or no 0.3-second flash at the predicted rate of roughly one per five years at amplification $10^3$, would rule out Eq. (2).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the axion-production mechanism in AQNs and the order-of-magnitude density estimate that this paper refines."},{"cited_title":"local flash","cited_arxiv_id":null,"evidence_quote":"Supplies the AQN hit rate, baryon-charge distribution, and energy-loss equations that set the simulation input."},{"cited_title":"DAMA/LIBRA annual modulation and Axion Quark Nugget Dark Matter Model","cited_arxiv_id":"1909.05320","evidence_quote":"Defines the broadband detection strategy whose annual, daily, and local-flash observables this paper computes."},{"cited_title":"Freese, J","cited_arxiv_id":null,"evidence_quote":"Provides the layered Earth density model used in the trajectory integration."},{"cited_title":"Andriamonje, S","cited_arxiv_id":null,"evidence_quote":"Fits the AQN baryon-charge distribution parameters from solar extreme-ultraviolet observations, defining the model variants."}],"review_version":1}