{"id":"dab22830-74af-4c5c-8a1c-ffacb0c99e35","arxiv_id":"2608.11740","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Magnetic noise from the softened spin sound mode of a two-component dark-exciton Bose-Einstein condensate is predicted to be detectable by NV center relaxometry, with a distinctive d log d scaling at short distances.","lead":"The paper predicts that a condensate of dark excitons, particles that cannot emit light, produces measurable magnetic noise that a diamond nitrogen-vacancy sensor can pick up. This would give experimenters a new way to spot a type of exciton condensate that is currently invisible to optical probes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted NV relaxation signal is overestimated by μ0² because the dipole-field prefactor μ0 in Eq. (57) is dropped in Eq. (58), so the central observability claim is not quantitatively supported as written.","rationale":"The paper's strongest claim is not merely that a spin-sound mode softens, but that this produces a measurable NV T1 signal with a distinctive d log d fingerprint. The second part depends on the absolute scale of 1/T1. Tracing the derivation, Eq. (57) is an SI dipole formula containing μ0/(4π); the next line (58) gives the Fourier-space coupling as i/2 e^{-kd} k_η, without μ0. This cannot be a harmless convention choice, because γ_NV=2μ_B/ℏ and the susceptibility χ_zz of Eq. (10) are also in SI units; a magnetic field noise spectrum in T²/Hz must carry μ0². Restoring μ0² would suppress the plotted rates by twelve orders of magnitude, making the central observability claim unsupported. I therefore do not base my verdict on the quasi-equilibrium or parameter-transfer assumptions, which are legitimate premises of a proposal paper, nor on the already-noted Eq. (10) sign and d log d factor issues, which affect spectral shape and fingerprint amplitude rather than the overall detection scale. The d log d coefficient error in SM Eq. (65) is real, but it is secondary to the μ0 prefactor. The appropriate verdict is CONDITIONAL: the author should correct the NV coupling prefactor and show that with the stated parameters the resulting 1/T1 is within NV sensitivity; otherwise the paper's experimental proposal cannot be accepted as quantitative.","tokens_in":18775,"tokens_out":28665,"duration_ms":313893,"concrete_test":"Independently re-evaluate the Fourier transform in SM §V by keeping the μ0/(4π) factor from Eq. (57): compute f^η_k = (3 μ0 d/4π) ∫ d²ρ e^{ik·ρ} ρ^η/(ρ²+d²)^{5/2} and check whether it equals i μ0 k_η e^{-kd}/2 rather than i k_η e^{-kd}/2. Then recompute Eq. (13) and Fig. 1 with the corrected μ0² prefactor and compare the resulting 1/T1 rates to the sensitivity of scanning NV relaxometry. If the corrected rates fall below the measurable window, the NV detection proposal is not viable; if the authors intended a unit system with μ0 absorbed, they must state that convention and convert γ_NV accordingly.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Most load-bearing concern: a missing μ0 in the NV–sample coupling. Eq. (57) defines the dipolar field with an explicit μ0/(4π). Fourier transforming to Eq. (58) should give f^η_k = i μ0 e^{-|k|d} k_η / 2, but the text writes f^η_k = i e^{-|k|d} k_η / 2. Substituted into the T1 formula (13), this removes a factor μ0² from the prefactor: γ_NV²/8 should be γ_NV² μ0²/8 if χ_zz is the spin susceptibility of Eq. (10) and B is in tesla. Fig. 1 quotes absolute rates in s^{-1}, so every plotted rate is inflated by μ0² ≈ 1.6×10^{-12}; for the stated parameters the corrected rates are far below NV relaxometry sensitivity, and the headline claim that the mode 'can be measured' loses its numerical support. This is an internal inconsistency at the step from Eq. (57) to Eq. (58), independent of the assumed BEC existence and independent of the d log d coefficient issue.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes detecting a two-dimensional Bose-Einstein condensate of optically dark S_z = ±1 triplet excitons through the magnetic field noise it emits, using the spin relaxation (T1) of a nearby NV center. For a two-component exciton gas with density interaction g and spin-anisotropy interaction g′, the authors derive the Bogoliubov spectrum of the antiferromagnetic (two-component) phase and show that the lower spin-sound mode softens as the Zeeman field approaches B_c = g′n/γ, entering the gigahertz band. They compute the retarded longitudinal spin susceptibility, include Beliaev damping that diverges as 1/√(1−f_z) near the critical point, insert both into the NV T1 formula, and predict both a measurable relaxation rate and, at small NV-sample distances d, a d log d scaling of the low-frequency rate arising from the cubic corrections to the dispersion and damping. Input parameters are taken from prior TMD experiments, with no fitting to the predicted noise; the claim is that NV relaxometry gives a viable, optically-blind detection channel for dark-exciton condensates.","tokens_in":19089,"tokens_out":45530,"duration_ms":450740,"significance":"The mechanism is well chosen and is the strongest part of the paper: the softening of the lower spin-sound branch at B_c derives from a standard two-component Bogoliubov treatment, the estimate B_c ≈ 2.3 T is concrete, and the Beliaev-damping divergence ∝ 1/√(1−f_z) is a physically interesting prediction. The derivation is transparent: the SM contains the full calculation, all inputs (g, g′, n, m, γ, Γ0) come from prior experiments or are stated as estimates, and no parameter is fitted to the predicted noise spectrum. The d log d small-distance scaling is a genuinely falsifiable fingerprint of the cubic corrections. However, the internal prefactor inconsistencies identified below affect precisely the quantities supporting the headline claims: the absolute rates in Fig. 1 and the coefficient quoted for the d log d term are not reliable, and the observability conclusion fails a direct SI-unit check. If the prefactor issues were corrected and the claims re-scoped, the theoretical content could be worth publishing; in the present form the central quantitative conclusion is not supported.","major_comments":[{"comment":"The stray-field coupling in the SM is missing the vacuum permeability. Equation (57) starts from B̂(r) = −(μ0/4π)[m̂/r³ − 3(r·m̂)r/r⁵], but Eq. (58) quotes f^η_k = (3d/4π)∫d²r e^{ik·r} r_η/r⁵ = (i/2)e^{−|k|d} k_η, i.e., the μ0 that is explicit in Eq. (57) drops out of the Fourier transform. The correct result is f^η_k = (iμ0/2)e^{−|k|d} k_η, so the prefactor in main-text Eq. (13) must be γ²_NV μ0²/8 and every absolute rate in Fig. 1 is overestimated by μ0² ≈ 1.6×10⁻¹². For the stated parameters the corrected peak rates are of order 10⁻¹⁴ s⁻¹, far below the sensitivity of NV T1 relaxometry; the headline claim that the softened mode 'can be measured' therefore has no quantitative support. This is an internal inconsistency between Eqs. (57) and (58) of the SM that is independent of whether the assumed dark-exciton BEC exists.","section":"SM §V, Eqs. (57)–(58); main text Eq. (13); Fig. 1"},{"comment":"The retarded susceptibility in Eq. (10) of the main text as typeset is ambiguous and needs to be reconciled with the SM. SM Eq. (23) gives the Matsubara correlator as 2γ² Q²_{η;k} ω_{ηk}/((iω_n)² + ω²_{ηk}), and the Lorentzian in SM Eq. (60) likewise has a numerator proportional to ω_q, so the retarded form should read 2γ² Σ_η Q²_{η;k} ω_{ηk}/((ω+i0⁺)² − ω²_{ηk}). If the intended form in Eq. (10) is Q²_{η;k}/(ω_{ηk}((ω+i0⁺)² − ω²_{ηk})), the spectral weight of the sound-mode resonance is off by 1/ω²_{ηk}, which changes all T1 integrals, including the low-frequency d log d analysis. Please correct Eq. (10) and confirm explicitly that Figs. 1–3 were computed with the Q²ω_{ηk} form.","section":"Main text Eq. (10); SM Eqs. (23) and (60)"},{"comment":"The coefficient of the d log d term in Eq. (65) does not follow from Eq. (63). In the UV regime ω_q = αq³ and Γ_q = Γ_B q³, so the integrand of Eq. (63) is q⁵ · 4ω_qΓ_q/(ω_q² + Γ_q²)² = 4αΓ_B/[(α² + Γ_B²)² q]; carrying out the integral gives I′_UV = 8d · αΓ_B/(α² + Γ_B²)² · log(1/(dk*)). Equation (65) quotes 8d · αΓ_B/(α² + Γ_B²) · log(1/(dk*)), which is off by one power of (α² + Γ_B²). This factor then propagates to Eq. (16) of the main text, so the claimed extraction of α/Γ_B from the amplitude of the d log d term would be quantitatively wrong by the same factor.","section":"SM §VI, Eqs. (63)–(65); main text Eq. (16)"}],"minor_comments":[{"comment":"The temperature entering the reduced rate 1/T1 · sinh(ω/2T) in Figs. 1 and 3 is never stated; please specify T in the figure captions and state explicitly the units of T in Eq. (13), noting the ℏ = k_B = 1 convention.","section":"Figs. 1 and 3; Eq. (13)"},{"comment":"Γ0 is introduced as 'Γ0 = 100 ps the extrinsic damping lifetime'; as a rate this means Γ0 = 1/(100 ps) ≈ 10 GHz, which is comparable to the probe frequencies of Fig. 2. The value used in the numerics should be stated and labeled on the horizontal dotted line of Fig. 2.","section":"Beliaev-damping section; Fig. 2"},{"comment":"Equation (16) is written as a proportionality; once the coefficient of the d log d term is corrected, this should be given as an explicit equality, since the ratio α/Γ_B is the quantity claimed to be extractable.","section":"Main text Eq. (16)"},{"comment":"The quantity Q_{ηk} appearing in Eq. (10) is defined only in the SM; please define it in the main text.","section":"Main text Eq. (10)"},{"comment":"Calling the d log d behavior an 'unconventional UV scaling' overstates the case: a logarithmic short-distance correction of this type is the standard consequence of a 1/q high-momentum tail regulated at q ≈ 1/d; consider rephrasing.","section":"Main text, Results and Conclusion"}],"recommendation":"reject","confidential_remarks":"The missing μ0 in SM Eqs. (58)–(59) and in main-text Eq. (13) is decisive in my assessment: it inflates the headline rates in Fig. 1 by μ0² ≈ 1.6×10⁻¹², and the corrected numbers place the predicted T1 rates orders of magnitude below what NV relaxometry can detect. The algebraic machinery of the paper (Bogoliubov spectrum, softening, Beliaev damping, d log d scaling) is largely sound and could support a revised manuscript with re-scoped claims, but the present version fails an internal SI-unit check that is independent of the assumed existence of the dark-exciton condensate. I also urge the editor to have the authors settle the main-text Eq. (10) versus SM Eqs. (23)/(60) discrepancy explicitly, since if the main-text form is the one used in the numerics, the spectral-weight analysis would need to be redone."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a genuinely new detection scheme for dark exciton condensates—using NV relaxometry to see the softening spin sound mode—and the d log d distance scaling is a nice fingerprint. But there is a load-bearing unit error in the NV coupling: the dipole field in Eq. (57) has the SI μ0/(4π), and the Fourier-transformed kernel in Eq. (58) drops the μ0. I checked the step; the stress-test is right. That makes every absolute rate in Fig. 1 too large by about 1/μ0² ≈ 6×10^11. Putting μ0 back, the predicted 1/T1 values fall far below what NV relaxometry can measure. This is an internal inconsistency in the calculation, not just an uncertainty in the BEC existence or parameter values. The qualitative softening picture survives, but the central observability claim, as written, is not quantitatively supported.\n\nWhat is genuinely good: the choice of S_z = ±1 dark excitons is well motivated by their long lifetime, and the field-tuned softening from AFM to FM brings the mode into a frequency range that, in principle, an NV can access. The derivation of Beliaev damping near the critical field and the resulting d log d scaling as a function of NV–sample distance is a distinctive prediction, and a clean fingerprint if the signal were large enough. The parameter set comes from TMD bilayer experiments, so the numbers are grounded.\n\nOther soft spots are minor by comparison. The Eq. (10) numerator looks like a typo relative to the SM (Q²ω_k vs Q²/ω_k) and should be clarified because it changes the spectral weight and the low-frequency tail. The assumptions of quasi-equilibrium, homogeneous density, and the transferred g, g', n values are strong; if the real condensate is inhomogeneous or short-lived, the softened mode may sit outside the NV window. The d log d coefficient is derived for a low-frequency limit that may be hard to reach in practice.\n\nThis is a paper for the exciton BEC and quantum sensing communities. The idea deserves a serious referee, but the μ0 error must be fixed and corrected rates reported before the observability claim can be taken seriously. I would send it to review, with a clear request to recalculate the absolute rates and discuss whether any parameter regime survives.","headline":"A clever new detection scheme for dark exciton condensates, but a missing μ0 in the NV dipole coupling inflates the predicted rates by orders of magnitude and undermines the central observability claim.","tokens_in":19597,"tokens_out":13039,"would_cite":false,"duration_ms":130483,"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":"An optically invisible exciton condensate can be detected by the magnetic noise it emits; a soft spin mode brings that noise into an NV center's range.","keywords":["dark excitons","Bose-Einstein condensate","NV center magnetometry","magnetic noise","Beliaev damping","T1 relaxometry","spinor BEC","transition metal dichalcogenides"],"falsifier":"Scan an NV center over a TMD bilayer while sweeping the external field from about 1.8 T up to $B_c \\approx 2.3$ T at a distance of 150 nm: the absence of a growing 1–10 GHz resonance as $B$ approaches $B_c$, or the failure of low-frequency small-distance data to follow $I_0 + I_1 d + I' d \\log d$, would falsify the prediction.","tokens_in":2203,"feed_emoji":"🧲","tokens_out":3058,"duration_ms":107418,"temperature":0.7,"pith_summary":"Dark $S_z = \\pm 1$ triplet excitons—bound electron-hole pairs whose optical decay is spin-forbidden—can form long-lived Bose-Einstein condensates that emit no light and therefore evade standard optical detection. This paper argues that such a condensate should still produce a detectable magnetic fingerprint: an external magnetic field along the spin quantization axis tunes the condensate from antiferromagnetic to ferromagnetic spin ordering, and the softer of the two sound modes dips into the gigahertz frequency range as the field approaches the critical value $B_c = g'n/\\gamma$ (about 2.3 T for parameters taken from recent TMD bilayer experiments). In that range, a nitrogen-vacancy (NV) center in diamond measures the stray magnetic-field noise through its spin relaxation rate $T_1^{-1}$. The paper derives the noise spectrum from the longitudinal spin susceptibility, including cubic corrections to the sound dispersion and Beliaev damping, and shows that the small-distance signal scales as $d \\log d$ rather than a plain power law. If correct, NV relaxometry turns an optically invisible condensate into a measurable magnetic object and provides a way to certify its existence.","feed_headline":"Magnetic noise can betray a dark-exciton condensate","feed_subtitle":"A field-driven spin mode softens into the gigahertz band, where NV centers can hear it.","key_machinery":"The engine of the calculation is the retarded longitudinal magnetic susceptibility $\\chi_{zz}(\\omega,k)$, whose poles sit exactly at the Bogoliubov sound-mode energies; this is what turns a condensate's spin dynamics into measurable magnetic noise. It feeds into the NV relaxometry formula $1/T_1 = (\\gamma_{\\mathrm{NV}}^2/8) \\coth(\\omega/2T) \\int d^2k/(2\\pi)^2 e^{-2kd} k^2 \\chi''_{zz}(\\omega,k)$, in which the exponential factor filters the susceptibility at a wavelength set by the NV–sample distance $d$. Around the critical field, the lower sound velocity $c_-$ vanishes, making the $k^3$ term of the dispersion and the Beliaev damping $\\Gamma_B k^3$ (decay of one sound quasiparticle into two) the dominant scales; their competition produces the logarithmic distance dependence. The parameters $g = 5\\,\\mu\\mathrm{eV}\\,\\mu\\mathrm{m}^2$, $g' = 0.2\\,\\mu\\mathrm{eV}\\,\\mu\\mathrm{m}^2$, $n = 2.0\\times10^{11}\\,\\mathrm{cm}^{-2}$, and $m = m_e$ are taken from a recent two-component exciton-condensate experiment in a TMD bilayer, giving $B_c \\approx 2.3$ T.","core_discovery":"The central claim is that the longitudinal spin sound mode of the $S_z = \\pm 1$ dark-exciton BEC softens as the external field approaches $B_c = g'n/\\gamma$, and this softening brings it into the frequency window of NV $T_1$ relaxometry. In the antiferromagnetic phase the two Bogoliubov modes have sound velocities $c_\\pm$ given by $c_\\pm^2 = [n(g+g') \\pm \\sqrt{n^2(g-g')^2 + 4f_z^2 n^2 gg'}]/(2m)$; as $f_z \\to 1$, $c_- \\to 0$ and the $k^3$ correction to the dispersion becomes important. The retarded longitudinal magnetic susceptibility has poles at these sound energies, and the NV relaxation rate $1/T_1 = (\\gamma_{\\mathrm{NV}}^2/8) \\coth(\\omega/2T) \\int d^2k/(2\\pi)^2 e^{-2kd} k^2 \\chi''_{zz}(\\omega,k)$ acts as a momentum filter selecting $k \\approx 1/d$. The paper further claims that Beliaev damping of the soft mode, $\\Gamma_B k^3$ with $\\Gamma_B \\propto (1-f_z)^{-1/2}$, dominates the linewidth near the transition and, together with the cubic dispersion, yields a small-distance scaling proportional to $I_0 + I_1 d + I' d \\log d$. This provides a route to detect exciton BECs through their spin component rather than their optical emission.","pith_inferences":["Inference: if the predicted $d \\log d$ tail were observed at several probe frequencies, one could map the crossover momentum $k_* = c/\\alpha$ as a function of field and thereby measure the cubic dispersion coefficient directly, a quantity the paper currently fixes by theory.","Inference: the same relaxometry could act as an equilibration probe, since a condensate that has not fully thermalized should produce a different frequency distribution of magnetic noise than the quasi-equilibrium spectrum computed here.","Inference: rotating the applied field away from the NV axis should suppress the signal, providing a control experiment that distinguishes exciton spin noise from other magnetic backgrounds."],"forward_implications":["At an NV–sample distance of 150 nm, fields as low as $0.99\\,B_c$ pull the softened sound mode into the 1–10 GHz band, where the predicted $T_1^{-1}$ spectrum becomes observable.","At small separation ($d \\lesssim 30$ nm), the low-frequency noise tail follows $I_0 + I_1 d + I' d \\log d$, with coefficient $I' \\propto \\Gamma_B \\alpha/(\\alpha^2 + \\Gamma_B^2)$, so fitting the tail extracts the ratio $\\alpha/\\Gamma_B$ as a function of the magnetic field.","Near the transition, Beliaev damping overtakes extrinsic damping at high momenta, making the condensate a setting in which to study two-dimensional many-body decay processes beyond ultracold atomic gases.","In the ferromagnetic phase the massive mode does not couple to the longitudinal susceptibility at zero temperature, and the remaining sound mode is too stiff to be detected by $T_1$ relaxometry, so the antiferromagnetic side of the transition carries the observable signature.","The scheme does not require optical emission from the excitons, so it can also certify condensates that may already exist in current bilayer experiments but have produced no optical signature."],"supporting_citations":[{"why":"Supplies the experimental parameter set ($g$, $g'$, $n$, $m$) for a two-component exciton condensate in a TMD bilayer that fixes $B_c \\approx 2.3$ T.","marker":"[23]"},{"why":"Provides the dark-exciton fine structure and lifetime evidence that motivates long-lived $S_z = \\pm 1$ condensates and the gyromagnetic ratio.","marker":"[15]"},{"why":"Supplemental Material containing the detailed derivations of the Bogoliubov modes, the magnetic susceptibility, Beliaev and Landau damping, the NV relaxometry formula, and the $d \\log d$ scaling.","marker":"[24]"},{"why":"Supplies the relaxometry relation between an NV center's $T_1$ rate and the magnetic noise of a two-dimensional magnet.","marker":"[38]"},{"why":"Establishes NV-center magnetometry as the experimental technique capable of measuring such stray-field noise.","marker":"[34]"},{"why":"Experimental observation of Beliaev coupling in a Bose gas, validating the decay channel the paper invokes.","marker":"[45]"},{"why":"Experimental measurement of Beliaev damping of quasiparticles in a condensate, supporting the damping mechanism.","marker":"[46]"},{"why":"Gives the $k^3$ Beliaev-damping scaling in two-dimensional dilute Bose gases that the paper extends to the critical regime.","marker":"[50]"}],"fun_headline_variants":["Dark exciton BEC leaks magnetic noise","Magnetic noise spots dark exciton condensate","NV centers hear dark exciton BEC's spin sound","Field-tuned spin mode exposes dark exciton BEC","Dark exciton BEC betrays itself via stray fields"],"cache_read_input_tokens":21760,"weakest_assumption_plain":"The whole prediction assumes that a homogeneous, long-lived condensate of dark excitons actually forms with the interaction strengths, density, and mass taken from recent bilayer experiments, and that the diamond sensor's spin axis is aligned with the direction the exciton spins point; if any of those fail, the softened mode may sit outside the gigahertz window and the signal would not show up.","fun_headline_variants_meta":{"raw":{"variants":["Dark exciton BEC leaks magnetic noise","Magnetic noise spots dark exciton condensate","NV centers hear dark exciton BEC's spin sound","Field-tuned spin mode exposes dark exciton BEC","Dark exciton BEC betrays itself via stray fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000382,"raw_usage":{"total_tokens":2114,"prompt_tokens":1120,"completion_tokens":994,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":736,"completion_tokens_details":{"reasoning_tokens":918}},"tokens_in":736,"tokens_out":994,"duration_ms":12139,"temperature":1.0,"reasoning_tokens":918,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:32:18.066845+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Scan an NV center over a TMD bilayer while sweeping the external field from about 1.8 T up to $B_c \\approx 2.3$ T at a distance of 150 nm: the absence of a growing 1–10 GHz resonance as $B$ approaches $B_c$, or the failure of low-frequency small-distance data to follow $I_0 + I_1 d + I' d \\log d$, would falsify the prediction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental parameter set ($g$, $g'$, $n$, $m$) for a two-component exciton condensate in a TMD bilayer that fixes $B_c \\approx 2.3$ T."},{"cited_title":"Robert, T","cited_arxiv_id":null,"evidence_quote":"Provides the dark-exciton fine structure and lifetime evidence that motivates long-lived $S_z = \\pm 1$ condensates and the gyromagnetic ratio."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplemental Material containing the detailed derivations of the Bogoliubov modes, the magnetic susceptibility, Beliaev and Landau damping, the NV relaxometry formula, and the $d \\log d$ scaling."},{"cited_title":"Flebus and Y","cited_arxiv_id":null,"evidence_quote":"Supplies the relaxometry relation between an NV center's $T_1$ rate and the magnetic noise of a two-dimensional magnet."},{"cited_title":"Hodby, O","cited_arxiv_id":null,"evidence_quote":"Experimental observation of Beliaev coupling in a Bose gas, validating the decay channel the paper invokes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental measurement of Beliaev damping of quasiparticles in a condensate, supporting the damping mechanism."}],"review_version":1}