{"id":"d2d15821-8834-4f12-a62c-bec7e77b0b3d","arxiv_id":"2411.10542","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Random nuclear spin orientations break momentum conservation, so axion absorption in crystals excites phonons across the whole Brillouin zone, yielding a broadband detection rate proportional to the phonon density of states.","lead":"The paper finds that axion dark matter in the 1 to 100 meV mass range can excite single quantized vibrations (phonons) in crystals with randomly oriented nuclear spins, producing a broad, continuous signal. This could give a new detection channel for a mass range that is currently very hard to probe experimentally.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (31) assumes degenerate spin configurations (Eq. (9), footnote 4); in molecular crystals H2, D2, H2O, and D2O exchange symmetry couples spin to rotation with meV-scale ortho-para gaps, so the low-mass projected rates and the broadband coverage claim are not established.","rationale":"The reader's weakest_assumption identifies the same load-bearing point: the derivation of Eq. (31) requires the crystal state to factorise into independent nuclear-spin and phonon degrees of freedom with all spin configurations degenerate. That assumption is legitimate for atomic crystals with isolated nuclei and negligible magnetic fields, where spin-projections are degenerate and hyperfine/dipolar splittings are far below the meV range. It is not legitimate for molecular crystals of identical nuclei, where exchange symmetry entangles the nuclear-spin state with the rotational state. The paper itself flags the degeneracy assumption in footnote 4, and the full text uses the phonon DOS of solid parahydrogen for H2, which is a spin-zero species; that strengthens the concern beyond a generic caveat. My read does not move the verdict: the proposed broadband mechanism for atomic crystals appears plausible and the derivation up to Eq. (31) is clean for those targets, so conditional acceptance remains appropriate. The fix is specific, namely to redo the molecular-target rates with spin-rotation states and energy gaps, and to verify which mass windows survive.","tokens_in":23120,"tokens_out":13303,"duration_ms":148231,"concrete_test":"Recompute the solid H2 rate of Eq. (31) using the exact two-proton molecular basis rather than the factorised ansatz of Eq. (9): take initial/final states |I=0,J=0> and |I=1,J=1> with Delta = E(J=1)-E(J=0) ~ 15 meV, use the para-H2 phonon DOS from [52], and impose energy conservation delta(ma - omega_k - Delta) for spin-changing transitions. Compare with the H2 curve in Fig. 2 (left panel); if the rate for ma below 15 meV is zero or suppressed by more than an order of magnitude, the low-mass H2/D2/D2O contours and the stated 1-100 meV coverage need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central rate formula (31) follows from replacing the final-spin sum by the identity after Eq. (25), which is valid only if every spin configuration reached by the axion is degenerate with the initial configuration (the assumption asserted in footnote 4 after Eq. (25)). This fails for the molecular crystals that dominate the projected rates. In H2 and D2, the two identical nuclei in each molecule are exchanged by molecular rotation, so the nuclear-spin and rotational Hilbert spaces are entangled: para-H2 has I=0 and is inert, while ortho-H2 has I=1, J=1, with E(J=1)-E(J=0) around 15 meV (about 7.5 meV for D2). A one-proton spin flip changes the total molecular spin, forcing a rotational transition, so energy conservation should be delta(ma - omega_k - Delta_rot), not delta(ma - omega_k). The same obstruction applies to H2O and D2O through ortho/para water and librational modes. There is also a material-specific inconsistency: the phonon DOS used for solid H2 is measured on parahydrogen [52], which has I=0, yet Eq. (28) and Table I assign xi_H = 0.75 GeV^-1 to every H site, treating spin-carrying ortho molecules as if they formed the same crystal. Consequently the low-mass portions of the H2/D2/H2O/D2O contours in Fig. 2 are not supported by the calculation as written, and the claim that 2-3 materials cover the 1-100 meV range is not established. Atomic crystals such as Be, Li2O, Al2O3, and GaAs do not have this exchange-symmetry obstruction, so the broadband mechanism itself is not called into question; the flaw is in the molecular realisations that give the largest rates.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new phonon-based direct-detection channel for axion dark matter in the meV to 100 meV mass range. Starting from the non-relativistic axion-nucleon interaction, the authors derive the rate for incoherent single-phonon emission in crystals with randomly oriented nuclear spins. Their central result, Eq. (31), states that the rate per unit mass is proportional to the axion mass, the axion-proton coupling squared, and the atom-projected phonon density of states, with no momentum-conservation restriction. They compute this rate for H2, D2, Al2O3, GaAs, H2O, D2O, Be, and Li2O, compare it with the previously studied coherent absorption, estimate the subleading contribution from the axion-induced nuclear electric dipole moment, and discuss the experimental prospects for TESSERACT-style phonon sensors and the relevant backgrounds. The headline claim is that only two or three target materials suffice to cover the 1-100 meV axion mass range without frequency scanning.","tokens_in":23487,"tokens_out":10885,"duration_ms":120847,"significance":"If Eq. (31) is valid for a given material, the paper makes an important methodological contribution: it shows that spin disorder in an unpolarised crystal removes the momentum-matching condition, so the axion absorption spectrum becomes proportional to the phonon density of states. This turns a narrowband search into a broadband one and gives falsifiable predictions for Be, Li2O, Al2O3, and GaAs, which do not involve any fitting parameters and depend only on the axion couplings and independently measured or computed phonon spectra. The coherent limit is checked against the existing literature, and the paper is careful about nuclear form factors and isotope averaging. However, the most promising projected rates and the '2-3 materials cover the full range' claim rely on molecular crystals (H2, D2, H2O, D2O), where a load-bearing spin-degeneracy assumption fails because identical-nucleus exchange symmetry couples nuclear spin to molecular rotation. The mechanism as rigorously proven is therefore restricted to atomic crystals, and the headline coverage claim is not yet established.","major_comments":[{"comment":"The derivation of Eq. (31) replaces the final-spin sum in Eq. (25) by the identity using the assertion in footnote 4 that all spin configurations have degenerate energies. This is the load-bearing assumption of the incoherent rate. It is valid for atomic crystals in the absence of magnetic fields, but it fails for the molecular crystals H2, D2, H2O, and D2O, where exchange symmetry of identical nuclei entangles nuclear spin with molecular rotation. For H2, para-H2 has total nuclear spin I=0 and even rotational J, while ortho-H2 has I=1 and odd J, with E(J=1)-E(J=0) about 15 meV; D2 has a comparable gap of several meV. A one-proton spin flip in a para molecule changes the total molecular spin and forces a rotational transition, so the energy-conserving delta function should involve delta(ma - omega_k - Delta_rot), not delta(ma - omega_k). The footnote's observation that the axion flips only a single spin limits which final states can contribute, but it does not imply that those states are degenerate with the initial state. Equations (25)-(31) and the low-mass rates for these molecular targets are therefore not established as written.","section":"Sec. II D, Eqs. (25)-(31), footnote 4"},{"comment":"The material-specific input for solid H2 is internally inconsistent even apart from the general degeneracy issue. The phonon density of states is taken from Ref. [52], which is a measurement on solid parahydrogen. Parahydrogen molecules have total nuclear spin I=0 and do not supply the independent spin-1/2 degree of freedom per proton that Eq. (28) and Table I assume when assigning xi_H = 0.75 GeV^-1 to every H site. A single-proton spin flip in such a molecule is an off-diagonal transition to an ortho-like state, not a reorientation within the degenerate manifold used in Eq. (27). The same objection applies to D2 and to the hydrogen/deuterium nuclei in H2O and D2O. Consequently, the low-mass parts of the H2, D2, H2O, and D2O contours in Fig. 2 are not supported by the calculation as written, and the claim that 2-3 materials cover the 1-100 meV range is not established. The rates for Be, Li2O, Al2O3, and GaAs do not have this defect because their nuclei are not subject to identical-nucleus exchange within the crystal unit.","section":"Table I, Eq. (28), and Fig. 2"}],"minor_comments":[{"comment":"The statement that random spin orientations 'break translation symmetry' is imprecise: the ensemble-averaged crystal is translation invariant, and what removes the momentum-conservation delta function is the suppression of the off-diagonal interference term in Eq. (19) by uncorrelated spin orientations. Consider rewording to avoid the impression that the ensemble itself is a disordered static configuration.","section":"Abstract and Sec. II D"},{"comment":"Section IV states that the aluminium pair-breaking threshold of 7.2 meV sets a theoretical lower limit on the detectable axion mass, yet Fig. 2 shows contours extending down to 1 meV. Please clarify which readout or sensor is assumed for masses below 7.2 meV, or restrict the plotted range accordingly.","section":"Sec. IV and Fig. 2"},{"comment":"There are several typographical issues: 'Brioullin zone' in the Introduction, 'gann' in Appendix B, and the notation in Eq. (34) for the spin operators appears garbled. A careful proofread would be helpful.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the atomic-crystal part of the calculation appears sound. The main issue is whether the authors can either justify the spin-degeneracy assumption for molecular crystals, which currently underpins the most promising projected rates, or substantially revise the claims about H2, D2, H2O, and D2O. I see no problem with novelty or citation practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about this paper. First, it identifies a genuinely new detection channel: in unpolarised crystals, the axion sees a target whose translation symmetry is broken by randomly oriented nuclear spins, so phonon momentum need not match the axion momentum and the absorption rate is proportional to the phonon density of states. That is real physics and a nice extension of the coherent phonon absorption of Mitridate et al. Second, the advertised rates for the best molecular targets—solid H2, D2, H2O, D2O—are not actually established by the calculation as written.\n\nThe derivation for atomic crystals is clean and standard FGR mechanics. The paper also deserves credit for including nuclear form factors, which the prior coherent-phonon paper omitted, and for being honest about the sensitivity: at the benchmark couplings the QCD axion rates are a few events per 10 kg-year, comparable to or weaker than existing astrophysical bounds. The paper does not oversell the discovery potential.\n\nThe soft spot is the spin-degeneracy assumption in Eq. (9) and the paragraph after Eq. (25). The derivation replaces the final-spin sum with the identity, which requires all spin configurations reached by the axion to be degenerate with the initial one. That holds for isolated atomic nuclei in non-magnetic crystals. It fails for the molecular crystals, where exchange symmetry of identical hydrogen/deuterium nuclei couples nuclear spin to molecular rotation. A spin flip changes the total molecular spin, so it must be accompanied by a rotational transition; energy conservation should be delta(m_a - omega_k - Delta_rot), not delta(m_a - omega_k). The ortho-para splittings are in the meV range—15 meV for H2, 7.5 for D2—right in the claimed sensitivity window. There is also a material-specific inconsistency: the phonon DOS used for solid H2 was measured on parahydrogen, which has I=0, yet Eq. (28) assigns a spin-dependent coupling to every H site. So the low-mass contours for the four molecular targets in Fig. 2 are not supported as written.\n\nNone of this undermines the central physical mechanism. The atomic crystals (Be, Li2O, Al2O3, GaAs) don't have the ortho-para obstruction, and the claim that random spin disorder produces broadband absorption is credible. The fix is clear: treat the molecular crystals with the spin-rotation level structure built in, or restrict the projected sensitivity to atomic targets.\n\nThis paper is for the TESSERACT-style single-phonon direct detection community and people working on axion-nucleon couplings. I would send it to a serious referee: the idea is new, the core derivation is sound, and the molecular crystal gap is fixable rather than fatal. My own confidence in the molecular rates would stay low until the rotational level structure is included, but this is not a desk reject.","headline":"New broadband axion-phonon mechanism is real for atomic crystals, but the molecular targets H2/D2/H2O/D2O need a fix for ortho-para spin-rotation coupling before their projected rates can be trusted.","tokens_in":24032,"tokens_out":2864,"would_cite":true,"duration_ms":25154,"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":"The paper claims that axion dark matter from 1 to 100 meV can excite phonons in crystals with randomly oriented nuclear spins, breaking translation symmetry so the rate tracks the phonon density of states without frequency scanning.","keywords":["axion dark matter","broadband phonon detection","incoherent absorption","nuclear spin","phonon density of states","meV-scale axions","crystal targets","QCD axion"],"falsifier":"Measure the phonon density of states of solid H$_2$ and D$_2$ at meV energies and compare a cryogenic calorimeter's axion-absorption spectrum with the prediction $R \\propto \\sum_j \\xi_j D_j(m_a)$: a rate that drops steeply below the ortho--para rotational splitting, rather than tracing $D(m_a)$, would falsify the degenerate-spin assumption.","tokens_in":22923,"feed_emoji":"📡","tokens_out":10017,"duration_ms":82992,"temperature":0.7,"pith_summary":"This paper claims that axion dark matter in the 1--100 meV mass range can be detected through incoherent phonon excitation in ordinary crystals, without magnetic fields or frequency scanning. The key mechanism is that randomly oriented nuclear spins break the translational symmetry of the crystal as seen by the axion, so momentum conservation no longer ties the phonon momentum to the tiny axion momentum; any phonon whose energy equals $m_a$ can be created. The absorption rate per unit mass then equals $R = 2\\pi \\rho_a m_a \\, g_p^2 \\sum_j \\xi_j D_j(m_a) / \\sum_j m_j$, proportional to the atom-projected phonon density of states $D_j(\\omega)$ rather than to a narrow resonance. The paper evaluates this rate for eight solids, finds that light nuclei with unpaired spins (hydrogen, deuterium, lithium, beryllium) are the best targets, and estimates QCD-axion rates of a few events per 10 kg-year exposure. If correct, this turns the difficult meV-to-100 meV axion window into a broadband search that needs only two or three target materials.","feed_headline":"Random nuclear spins make axion absorption broadband","feed_subtitle":"The rate follows the phonon density of states, so a few crystals cover 1-100 meV axions.","key_machinery":"The load-bearing object is the atom-projected phonon density of states, $D_j(\\omega) = (1/3N)\\sum_{\\nu k} |\\epsilon^j_{\\nu k}|^2 \\delta(\\omega - \\omega_{\\nu k})$, together with the lattice-averaged spin coupling $\\xi_j = (1/N)\\sum_\\ell \\lambda_{\\ell j}^2 J_{\\ell j}(J_{\\ell j}+1)/m_{\\ell j}$. The argument splits the absorption rate into coherent and incoherent terms: the coherent terms carry the phase $e^{-i\\mathbf{k}\\cdot(\\ell-\\ell')}$ that enforces $\\mathbf{k}\\approx 0$ in a polarised crystal, while the incoherent term has no such phase after the random-spin ensemble average. Because the final-spin sum becomes a resolution of the identity, only $D_j(\\omega)$ and $\\xi_j$ control the rate, and since $D_j(\\omega)$ has support across a wide range of phonon energies, the detection scheme is broadband.","core_discovery":"The discovery is that axion absorption in an unpolarised crystal is broadband because the random spin configuration acts as disorder that fully breaks translation symmetry in the phonon effective theory. In a polarised crystal the coherent terms in the rate sum to a momentum-conservation condition $\\delta_{\\mathbf{k},0}$, so only zero-momentum optical phonons can absorb the axion and the response is a narrow Breit-Wigner peak. In an unpolarised crystal the ensemble average over spin configurations kills those interference terms, leaving the diagonal, incoherent term with no momentum-restricting phase; the sum over final spin states collapses to the identity because all spin configurations are degenerate in energy. The result is Eq. (31), $R = 2\\pi \\rho_a m_a \\, g_p^2 \\sum_j \\xi_j D_j(m_a)/\\sum_j m_j$, where $\\xi_j$ is the lattice-averaged spin coupling of the $j$-th nucleus and $D_j(\\omega)$ is its atom-projected phonon density of states. The same mechanism operates for the axion-induced time-dependent nuclear electric dipole moment, and the paper shows that this second channel is subdominant in standard QCD axion models.","pith_inferences":["Beyond the paper: if the degenerate-spin assumption fails in molecular solids through ortho--para exchange couplings, the promised low-mass coverage for H$_2$, D$_2$, H$_2$O and D$_2$O would shift, and the rate formula would need a spin-flip energy threshold.","Beyond the paper: the same 'disorder makes absorption broadband' principle suggests that other forms of quenched disorder, such as isotope mixtures, impurities, or amorphous targets, could generate broadband phonon absorption for other dark matter candidates.","Beyond the paper: a testable extension is to measure the partial phonon density of states of solid H$_2$ and D$_2$ with inelastic neutron scattering and use it to predict the full spectral shape that a future single-phonon calorimeter should see.","Beyond the paper: for axion-like particles with suppressed photon coupling but large nucleon coupling, this channel may be the only terrestrial probe in the 1--100 meV window."],"forward_implications":["A single unpolarised crystal target can search the full energy range where its phonon density of states has support, with no frequency scanning or applied magnetic field.","Two or three light-nucleus targets (hydrogen- or deuterium-bearing solids) cover essentially the whole 1--100 meV axion mass window.","The previously studied coherent, narrow-band optical-phonon channel is not needed for sensitivity; the incoherent channel gives a comparable or larger rate in unpolarised crystals.","For standard QCD axion models, the spin coupling channel dominates over the induced nuclear electric dipole channel, so the rate is controlled by $\\xi_j D_j(m_a)$.","Reaching the predicted QCD-axion rate of a few events per 10 kg-year requires a corresponding reduction of low-energy backgrounds in single-phonon calorimeters."],"supporting_citations":[{"why":"Supplies the coherent single-phonon absorption calculation that this paper recomputes and contrasts with the new incoherent channel.","marker":"[44]"},{"why":"Provides the relations between the QCD axion couplings $g_p,g_n$ and the UV Wilson coefficients, used to present rates model-independently.","marker":"[47]"},{"why":"Gives the non-relativistic reduction of the axion-nucleon Hamiltonian that produces the spin-momentum operator used in the rate calculation.","marker":"[48]"},{"why":"Supplies the phonon density of states for H$_2$O and D$_2$O used in the rate calculation.","marker":"[51]"},{"why":"Supplies the phonon density of states of solid parahydrogen.","marker":"[52]"},{"why":"Supplies the phonon density of states of solid deuterium.","marker":"[53]"},{"why":"Provides the odd-group-model nuclear spin form factors used to match axion-nucleon to axion-nucleus couplings.","marker":"[88]"},{"why":"Provides ab initio structure factors for spin-dependent scattering used for the $^{27}$Al form factor.","marker":"[90]"}],"fun_headline_variants":["Axion absorption goes broadband with random spins","Random spins broaden axion-phonon absorption","Spin disorder unlocks broadband axion detection","Axion phonons get broadband via random nuclear spins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that flipping a nuclear spin costs no energy beyond the phonon it creates, an assumption that molecular crystals with identical nuclei, where exchange symmetry couples spin to molecular rotation, may violate.","fun_headline_variants_meta":{"raw":{"variants":["Axion absorption goes broadband with random spins","Random spins broaden axion-phonon absorption","Spin disorder unlocks broadband axion detection","Axion phonons get broadband via random nuclear spins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1535,"prompt_tokens":1008,"completion_tokens":527,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":624,"completion_tokens_details":{"reasoning_tokens":470}},"tokens_in":624,"tokens_out":527,"duration_ms":4910,"temperature":1.0,"reasoning_tokens":470,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:36:10.285619+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the phonon density of states of solid H$_2$ and D$_2$ at meV energies and compare a cryogenic calorimeter's axion-absorption spectrum with the prediction $R \\propto \\sum_j \\xi_j D_j(m_a)$: a rate that drops steeply below the ortho--para rotational splitting, rather than tracing $D(m_a)$, would falsify the degenerate-spin assumption.","supporting_citations":[{"cited_title":"Colognesi, M","cited_arxiv_id":null,"evidence_quote":"Supplies the phonon density of states of solid parahydrogen."},{"cited_title":"Engel and P","cited_arxiv_id":null,"evidence_quote":"Provides the odd-group-model nuclear spin form factors used to match axion-nucleon to axion-nucleus couplings."}],"review_version":1}