{"id":"0bbe766b-156d-4ae6-8540-de9617fea2e6","arxiv_id":"2608.05620","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Positional disorder of manganese impurities strongly changes both the static and dynamic magnetic susceptibility of diluted magnetic semiconductors at low carrier densities.","lead":"This paper calculates how the magnetic response of a diluted magnetic semiconductor changes when the magnetic impurities are placed randomly rather than in an ordered lattice. It predicts that even slight randomness gives a qualitatively different, broad magnetic response that experiments could detect.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universality claim rests on untested robustness of χ(q,ω) to disordered on-site potentials, Hubbard U, and multi-band structure.","rationale":"The reader's weakest assumption correctly identifies the robustness of the simplified model as the most load-bearing gap: the paper's central claim is a universality statement, yet the numerical evidence is generated only for a minimal impurity-band model. My reading of the manuscript confirms that Section II justifies neglecting charged-defect disorder, Hubbard U, and multi-band structure by citing Ref. 11, which studied magnetization and Tc, not the susceptibility quantities central to this paper. The Conclusions extend the claim to 'all x below and near the MIT' without a direct test. I also note the manuscript's own admission that the disorder-averaged curves are 'not yet smooth,' which compounds the uncertainty. A concrete computational extension, as proposed, would settle whether the predicted signatures are robust or an artifact of the simplified model. The paper is otherwise internally coherent: the RPA derivation is standard, the ordered-case limits reproduce known results, and the q-independence mechanism (localized spin-wave modes excited regardless of q) is physically reasonable. Thus the reader's CONDITIONAL verdict is appropriate; no change is needed, provided the stated conditions — robustness checks and/or stronger caveats — are addressed.","tokens_in":17391,"tokens_out":14755,"duration_ms":166385,"concrete_test":"Add to the impurity-band Hamiltonian (1) a random on-site energy term Σ_i ε_i n_i with ε_i drawn from a Gaussian of width W ~ 0.5–1 eV (modeling charged compensation centers) and a Hubbard U ~ 1 eV; recompute the disorder-averaged χ(q,ω) and static χ(T) at x=0.00924, p=0.1, using the same RPA framework and at least 100 disorder realizations with error bars. If the broad, roughly q-independent spectrum (for q≥q2 in Fig. 8) and the low-T free-spin plus high-T mean-field peaks survive, the universality claim is supported; if the q-dependence or peak structure is qualitatively altered, the conclusions must be restricted to the minimal model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result is the universality claim in the Conclusions: 'The only ingredient necessary for this dramatic change in the shape of χ(q,ω) is the existence of some charge carrier localized states.' Within the specific impurity-band model, the RPA calculation supports a broad, q-independent response for moderate-to-large q. However, the leap to real (III,Mn)V DMS assumes that the neglected terms listed in Section II (random on-site energies from charged compensation centers, a Hubbard U, longer-range interactions, and the multi-band acceptor structure with spin-orbit coupling) lead only to quantitative changes. This assumption is validated in Ref. 11 only for magnetization curves and Tc, not for the dynamic transverse susceptibility or the two-peak longitudinal response computed here. It is conceivable that a strong random on-site potential changes the nature of localized states (e.g., producing a different energy landscape) or that spin-orbit coupling creates anisotropic/non-collinear local moments, altering the RPA spin-wave spectrum and the q-independence. Thus the central claim's generality is underdetermined by the presented evidence, and the manuscript itself acknowledges the disorder averages are 'not yet smooth' (Fig. 9 caption), so even the model-specific result lacks converged statistics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper calculates the static longitudinal and the transverse dynamic magnetic susceptibility of an impurity-band model for (III,Mn)V diluted magnetic semiconductors with positional disorder, using self-consistent mean-field theory and the RPA. In ordered configurations the static susceptibility has a single peak at the mean-field Tc and the dynamic transverse susceptibility has momentum-resolved spin-wave poles. In disordered configurations the static susceptibility develops an additional low-temperature peak attributed to weakly coupled Mn spins, and the dynamic susceptibility becomes broad and approximately q-independent over the spin-wave spectrum. The authors conclude that this behavior requires only the existence of localized charge carriers and predict it for all x below and near the metal-insulator transition.","tokens_in":17591,"tokens_out":9487,"duration_ms":99168,"significance":"If the central claim survives closer scrutiny, the paper would provide concrete, experimentally testable signatures of disorder-induced localization in DMSs: a broad, roughly q-independent transverse dynamic susceptibility observable in neutron scattering, and a two-peak static longitudinal susceptibility. The derivation has several internal checks: the ordered-case formulas reduce to the Pauli susceptibility, the RPA denominator yields a mean-field Tc consistent with known scaling, and the poles match the spin-wave spectrum of Ref. 13. The static disordered case is reduced to a closed system of linear equations. The paper also gives a useful caution about the conventional fluctuation formula for the susceptibility in mean-field Monte Carlo analyses. The main limitation is that the headline universality claim goes beyond the evidence presented for the specific model.","major_comments":[{"comment":"The statement in the Conclusions that \"the only ingredient necessary for this dramatic change in the shape of χ(q,ω) is the existence of some charge carrier localized states\" is a leap from the model calculation. Section II lists several neglected terms — random on-site energies from charged compensation centers, a Hubbard U, longer-range interactions, and the multi-band acceptor structure — and asserts, citing Ref. 11, that they lead only to quantitative changes. Ref. 11, however, validated that claim for magnetization curves and Tc, not for the dynamic transverse susceptibility or the two-peak longitudinal response. Because the prediction is made for all x below and near the MIT, where charged compensation centers are always present, this is load-bearing. Please either extend the RPA calculation to include a random on-site potential (the authors note the formalism can be straightforwardly generalized) and verify that the q-independent peak and the two-peak χ(T) survive, or explicitly restrict the conclusion to the positional-disorder-only model.","section":"Section VI (Conclusions); Section II"},{"comment":"The evidence for the central q-independence claim rests on disorder averages that the caption itself describes as \"not yet smooth\" (15 to 40 realizations, Nd = 125 or 216). The inset of Fig. 9 compares Nd=125 and Nd=216 using different q-vectors and invokes the q-independence claim to interpret the comparison; it therefore does not independently test q-independence or finite-size effects. Please provide convergence data with respect to the number of realizations and system size, ideally with error bars, and present a size comparison at the same q.","section":"Section V, Fig. 9 (and Fig. 8)"},{"comment":"The RPA derivation starts from a collinear mean-field state (footnote 25) and neglects spin-orbit coupling, which the Introduction notes may be significant. The Conclusions assert that adding anisotropies due to strain or spin-orbit coupling \"can only lead to quantitative changes,\" but no calculation or estimate supports this for the dynamic susceptibility. Non-collinear or anisotropic local moments could introduce gaps or alter the q-dependence of the spin-wave response. Please provide a concrete estimate (e.g., the spin-orbit splitting relative to J and the bandwidth) or a calculation for a simple anisotropic term, or soften the claim to a conjecture.","section":"Section IV and Section VI"}],"minor_comments":[{"comment":"The title in the manuscript text contains a typo: \"semiconducto rs\" should be \"semiconductors\".","section":"Title"},{"comment":"The horizontal axis label \"hω\" should be \"ħω\" to match the notation used in the equations.","section":"Figs. 5-9"},{"comment":"The paper uses \"weak,\" \"moderate,\" and \"full\" disorder without defining them; please add a sentence or a precise reference to Ref. 10 so the reader can interpret the disorder levels in the figures.","section":"Section II; Figs. 2, 8, 9"},{"comment":"The caption does not state the disorder level of the \"disordered sample\"; please specify whether it is weakly, moderately, or fully disordered.","section":"Fig. 4 caption"},{"comment":"The claim that the \"conventional\" fluctuation formula \"gives very wrong results for T ∼ T_c\" would be more convincing if a direct numerical comparison between that formula and the full solution were shown in a figure.","section":"Section III, discussion after Eq. (28)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-executed model calculation with good internal checks, but the published conclusion is broader than the evidence. I would ask the authors to either add the straightforward extension to on-site disorder or temper the universality claim, and to address the disorder-averaging convergence. No concerns about novelty or attribution; the authors are careful to reference their prior work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Solid, careful RPA calculation with a real result: positional disorder of Mn turns the transverse dynamic susceptibility from a sharp spin-wave peak into a broad, roughly q-independent response, because disorder-localized carriers couple the field to all spin-wave modes. The ordered-case limits check out (Pauli susceptibility, Goldstone mode, spin-wave spectrum from their earlier work), and the static longitudinal susceptibility develops a two-peak structure (free-spin peak from weakly coupled spins, mean-field peak at Tc). The calculation is parameter-free in the sense that J, aB, Ry are taken from literature, not fitted to these susceptibilities. The authors also make a genuinely useful methodological point: the conventional formula chi = beta<S.S> evaluated with the mean-field density matrix fails near Tc because it drops hole-polarization and cross terms.\n\nSoft spots: (1) The Conclusions claim that localized carrier states are 'the only ingredient necessary' and that the effect survives 'at all x below and near the MIT.' Within the impurity-band model, the claim is supported. But the model omits on-site disorder from charged compensation centers, Hubbard U, longer-range repulsion, and the multi-band acceptor structure with spin-orbit coupling. The paper asserts, citing Ref. 11, that these change little. Ref. 11 checked magnetization curves and Tc, not the dynamic susceptibility or the two-peak longitudinal response. So the generality claim is an extrapolation, not a demonstrated result. This matters because the paper's own falsifiable neutron-scattering prediction is framed in terms of real samples. (2) The numerics are real but limited: Fig. 9's caption admits the disorder averaging is 'not yet smooth', and the central q-independence plot (Fig. 8) uses only 15 configurations with no error bars. The qualitative effect is convincing, but the claim of q-independence for weak disorder deserves more averaging. Minor: some finite-size checks are good (inset of Fig. 9), others are single-configuration anecdotes (Figs. 3 and 7).\n\nWho this is for: anyone working on DMS magnetism, disordered magnets, or RPA response of inhomogeneous systems. It is a legitimate subfield advance, internally consistent, and it makes a testable prediction. It deserves a serious referee. I would send it to review with a request that the authors either soften the universality language or test the robustness (for example, add random on-site energies and see if the q-independence survives), and add more disorder realizations with error bars to the key figures. My verdict: accept after these revisions; the core result holds within the model.","headline":"A careful RPA calculation that makes a sharp falsifiable prediction about disorder-induced q-independent susceptibility, wrapped in a universality claim broader than the model supports.","tokens_in":18119,"tokens_out":2441,"would_cite":true,"duration_ms":27357,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.50.Pp","75.40.Gb","75.25.+z"],"model":"deepseek-v4-flash","headline":"The paper argues that even very little positional disorder of Mn impurities turns the transverse dynamic magnetic susceptibility of diluted magnetic semiconductors into a broad, roughly wavevector-independent response, because localized…","keywords":["diluted magnetic semiconductors","magnetic susceptibility","random phase approximation","positional disorder","impurity band","spin waves","GaMnAs","metal-insulator transition"],"falsifier":"A neutron-scattering measurement of the transverse dynamic susceptibility of a Ga$_{1-x}$Mn$_x$As film with $x$ near the metal-insulator transition (for example, $x \\approx 0.03$–$0.05$) that resolves sharp, $\\mathbf{q}$-dependent spin-wave peaks at finite wavevectors, instead of a broad response spanning the spin-wave spectrum at every $\\mathbf{q}$, would contradict the paper's central claim.","tokens_in":17142,"feed_emoji":"🧲","tokens_out":6159,"duration_ms":68789,"temperature":0.7,"pith_summary":"The paper argues that positional disorder of Mn impurities—not just carrier density or exchange strength—controls the magnetic response of diluted magnetic semiconductors at low carrier densities. Using a random-phase-approximation treatment of an impurity-band model, it shows that even weak disorder produces localized charge-carrier states, which in turn localize spin-wave modes; the result is a transverse dynamic susceptibility $\\chi(\\mathbf{q},\\omega)$ that is broad in frequency and roughly independent of wavevector $\\mathbf{q}$, instead of the sharp $\\mathbf{q}$-dependent spin-wave peaks of an ordered sample. It also predicts that the static longitudinal susceptibility of disordered samples has two temperature peaks, one from weakly coupled 'free' Mn spins and one marking local ferromagnetic cluster formation. This matters because magnetic properties, unlike transport, depend on all occupied states, so disorder should affect GaMnAs and related materials even on the metallic side of the metal-insulator transition, and the prediction can be checked by neutron scattering.","feed_headline":"Weak Mn disorder makes magnetic response momentum-blind","feed_subtitle":"Predicted broad, q-independent spin response below and near the metal-insulator transition, testable by neutron scattering.","key_machinery":"The central object is a variational mean-field plus generalized random-phase-approximation (RPA) treatment of a kinetic-exchange impurity-band Hamiltonian $H = \\sum_{i,j,\\sigma} t_{ij} c^\\dagger_{i\\sigma} c_{j\\sigma} + \\sum_{i,j} J_{ij} \\mathbf{S}_i \\cdot \\mathbf{s}_j$, in which Mn spins ($S=5/2$) occupy disordered positions and interact with holes through an exponentially decaying exchange $J_{ij}$. The machinery that carries the argument is the set of linear response equations for the transverse spin deviations $\\delta S_+(i,\\omega)$ at finite temperature, derived from an action principle; these show that a $\\mathbf{q}$-dependent external field couples to every spin-wave mode once translational invariance is broken. For the static longitudinal susceptibility, the analogous site-resolved system separates the response into strongly coupled clusters and weakly coupled 'free' spins. The mechanism is that positional disorder localizes some charge carriers, which localizes spin-wave modes, which removes momentum conservation and makes $\\chi(\\mathbf{q},\\omega)$ roughly $\\mathbf{q}$-independent.","core_discovery":"On its own terms, the paper establishes that the only ingredient needed to change the shape of the transverse dynamic susceptibility from a Lorentzian centered at a well-defined spin-wave frequency to a broad, roughly $\\mathbf{q}$-independent peak extending over the entire spin-wave spectrum is the existence of some charge-carrier localized states; on general grounds such states exist at all Mn concentrations below and near the metal-insulator transition. In the same impurity-band model, the static longitudinal susceptibility of disordered samples develops a two-peak temperature structure—a low-temperature $1/T$-like peak from weakly coupled Mn spins far from the hole-rich regions and a higher-temperature mean-field peak marking the polarization of strongly coupled clusters—with the higher peak broadening and shifting as disorder increases. The paper further claims that conventional fluctuation-based formulas for the susceptibility fail near $T_c$ because they omit the hole-mediated channel, and that these results are consistent with earlier studies of the inhomogeneous ferromagnetic state and spin-wave spectrum of these materials.","pith_inferences":["If the $\\mathbf{q}$-independence is generic, a clean experiment would compare a disordered alloy with an ordered or digitally doped heterostructure of the same composition: the ordered one should show dispersive spin-wave peaks, the disordered one a broad continuum, isolating disorder as the cause.","The two-peak longitudinal susceptibility implies that simple magnetization or susceptibility measurements may misidentify the transition temperature; combined ac-susceptibility and neutron or muon-spin-rotation experiments could separate $T^*$ from the ordering temperature.","The same mechanism—localized carriers from positional disorder inducing localized spin-wave modes and a $\\mathbf{q}$-independent response—should apply to other disordered impurity-band ferromagnets, such as GaMnN or GeMn, if their ferromagnetism is carrier-mediated, giving a testable family of materials.","The single-realization narrow peaks suggest mesoscopic cluster-level structure, so small-sample local probes could observe individual cluster polarization temperatures."],"forward_implications":["Neutron scattering on GaMnAs samples below and near the metal-insulator transition should observe a broad, roughly wavevector-independent dynamic response spanning the whole spin-wave spectrum, not sharp $\\mathbf{q}$-dependent spin-wave peaks.","Measured static susceptibility of disordered samples should show two peaks in temperature: a low-temperature free-spin peak and a higher-temperature peak at the mean-field cluster-polarization temperature $T^*$, with the higher peak broadening as disorder increases.","The higher-temperature susceptibility peak marks local ferromagnetic cluster formation, not the true long-range-ordering transition, so Monte Carlo studies that identify $T_c$ from the susceptibility should use the full susceptibility including hole-Mn correlations.","Disorder raises the mean-field $T_c$ relative to the ordered case, and individual realizations show multiple narrow peaks, one per strongly coupled cluster, implying sample-to-sample variation.","Because only localized carriers are needed, the qualitative predictions should survive adding on-site disorder, Hubbard $U$, longer-range interactions, strain, or spin-orbit anisotropy, with only quantitative changes."],"supporting_citations":[{"why":"Supplies the disordered impurity-band model and the weak, moderate, and full disorder definitions used throughout the paper.","marker":"[10]"},{"why":"Provides the self-consistent mean-field solution, the localized-carrier phenomenology, and the argument that additional terms are only quantitative.","marker":"[11]"},{"why":"Provides the spin-wave spectrum showing localized modes at low and high energies, the basis for the q-independent susceptibility interpretation.","marker":"[13]"},{"why":"The Monte Carlo study whose approximation $\\chi \\approx \\chi_{Mn}$ is shown to be questionable near $T_c$.","marker":"[12]"},{"why":"Supplies the hopping amplitude formula that defines the impurity band with a mobility edge.","marker":"[22]"},{"why":"Supplies the material parameters ($J=15$ meV, 1 Ry $\\approx 110$ meV, $a_B \\approx 8$ Å) used in the numerics.","marker":"[23]"},{"why":"Establishes the experimental context, including compensation and the metal-insulator transition near $x \\approx 0.03$.","marker":"[1]"}],"fun_headline_variants":["Mn disorder flattens spin response across momenta","Spin response goes momentum-blind with Mn disorder","Disorder turns spin-wave peak into broad q-independent hump","Two-peak susceptibility emerges from Mn disorder","Hole-mediated channel breaks conventional susceptibility near Tc"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction that disorder qualitatively reshapes the susceptibility rests on the premise that the complications left out of the model—charged compensation centers, electron-electron repulsion, multi-band acceptor structure, and strong spin-orbit coupling—change the results only quantitatively and never eliminate the localized charge-carrier states.","fun_headline_variants_meta":{"raw":{"variants":["Mn disorder flattens spin response across momenta","Spin response goes momentum-blind with Mn disorder","Disorder turns spin-wave peak into broad q-independent hump","Two-peak susceptibility emerges from Mn disorder","Hole-mediated channel breaks conventional susceptibility near Tc"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000647,"raw_usage":{"total_tokens":2897,"prompt_tokens":798,"completion_tokens":2099,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":414,"completion_tokens_details":{"reasoning_tokens":2026}},"tokens_in":414,"tokens_out":2099,"duration_ms":20732,"temperature":1.0,"reasoning_tokens":2026,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T05:17:49.146332+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A neutron-scattering measurement of the transverse dynamic susceptibility of a Ga$_{1-x}$Mn$_x$As film with $x$ near the metal-insulator transition (for example, $x \\approx 0.03$–$0.05$) that resolves sharp, $\\mathbf{q}$-dependent spin-wave peaks at finite wavevectors, instead of a broad response spanning the spin-wave spectrum at every $\\mathbf{q}$, would contradict the paper's central claim.","supporting_citations":[{"cited_title":"Berciu and R","cited_arxiv_id":null,"evidence_quote":"Supplies the disordered impurity-band model and the weak, moderate, and full disorder definitions used throughout the paper."},{"cited_title":"Berciu and R","cited_arxiv_id":null,"evidence_quote":"Provides the self-consistent mean-field solution, the localized-carrier phenomenology, and the argument that additional terms are only quantitative."},{"cited_title":"Berciu and R","cited_arxiv_id":null,"evidence_quote":"Provides the spin-wave spectrum showing localized modes at low and high energies, the basis for the q-independent susceptibility interpretation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The Monte Carlo study whose approximation $\\chi \\approx \\chi_{Mn}$ is shown to be questionable near $T_c$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the hopping amplitude formula that defines the impurity band with a mobility edge."},{"cited_title":"Bhattacharjee and C.B","cited_arxiv_id":null,"evidence_quote":"Supplies the material parameters ($J=15$ meV, 1 Ry $\\approx 110$ meV, $a_B \\approx 8$ Å) used in the numerics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the experimental context, including compensation and the metal-insulator transition near $x \\approx 0.03$."}],"review_version":1}