Pith. sign in

REVIEW 3 major objections 5 minor 36 references

Magnetic susceptibility of diluted magnetic semiconductors at low carrier densities

T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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…

desk verdict 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. read the letter →

arxiv 2608.05620 v1 pith:MHGXYTNW submitted 2026-08-06 cond-mat.mtrl-sci

Adel Kassaian , Mona Berciu This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 75.50.Pp75.40.Gb75.25.+z
keywords dilutedmagneticsemiconductorssusceptibilityrandomphaseapproximationpositionaldisorderimpuritybandspinwavesGaMnAsmetal-insulatortransition
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Section VI (Conclusions); Section II] 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.
  2. [Section V, Fig. 9 (and Fig. 8)] 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.
  3. [Section IV and Section VI] 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.
minor comments (5)
  1. [Title] The title in the manuscript text contains a typo: "semiconducto rs" should be "semiconductors".
  2. [Figs. 5-9] The horizontal axis label "hω" should be "ħω" to match the notation used in the equations.
  3. [Section II; Figs. 2, 8, 9] 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.
  4. [Fig. 4 caption] The caption does not state the disorder level of the "disordered sample"; please specify whether it is weakly, moderately, or fully disordered.
  5. [Section III, discussion after Eq. (28)] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the susceptibilities are parameter-free RPA outputs from a fixed impurity-band Hamiltonian, and no prediction is fed back into the model.

full rationale

The static longitudinal and transverse dynamic susceptibilities are computed by a fixed RPA/mean-field response formalism from the impurity-band Hamiltonian of Eq. (1). The parameters (J = 15 meV, aB = 8 Å, 1 Ry = 110 meV) are taken from prior literature (Refs. 10,22,23) and are not adjusted to reproduce the computed susceptibilities; no quantity computed here is fed back into the model. The ordered-case χ(q,ω) of Eq. (48) is an explicit RPA expression whose poles reproduce the spin-wave spectrum of Ref. 13, and the disordered-case result is obtained by numerically solving the linear system Eq. (44) for fixed disorder realizations. The two-peak longitudinal susceptibility and the q-independent broad transverse response are outputs, not inputs. The only self-citations (Refs. 10–15) provide the model, mean-field solution, and localization/spin-wave results used as inputs or consistency checks; they are not invoked as a uniqueness theorem and do not define the susceptibility in terms of itself. The Conclusions' universality claim does rely on the auxiliary assumption, stated in Sec. II and drawn from Ref. 11, that on-site disorder, Hubbard U, and longer-range interactions produce only quantitative changes; and Fig. 9's caption acknowledges that disorder averages are 'not yet smooth.' These are robustness and convergence caveats, not circular reductions: within the stated model the calculation is self-contained. No step equates a prediction to a fitted parameter or imports the conclusion via definition.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

No new particles, forces, or entities are introduced. The model uses established Hamiltonians with parameters from prior literature. The only free parameter is the numerical broadening η. The axioms are standard domain assumptions of the impurity-band/RPA approach, plus one ad hoc robustness assumption about neglected terms.

free parameters (1)
  • broadening η = η=0.02 J (ordered); η=0.05 J (disordered)
    Artificial broadening chosen to avoid numerical singularities in χ(q,ω). Its value affects peak widths; the ω→0 limit is checked by varying η (Sec. V). Not fitted to experimental data.
assumptions (5)
  • domain assumption The impurity band model with s-wave impurity states and parameters t(r), J(r) is appropriate for low carrier density (III,Mn)V DMS.
    Invoked in Sec. II; the model is justified qualitatively by comparison with known GaMnAs behavior, but ignores acceptor-level structure and multi-band physics.
  • ad hoc to paper Adding on-site disorder from charged defects, Hubbard U, and longer-range interactions changes results only quantitatively.
    Stated in Sec. II based on Ref. 11; used to justify discarding these terms. Not re-derived here and is load-bearing for the generality claim.
  • domain assumption The self-consistent mean-field ground state is collinear.
    From Ref. 11, cited in Ref. 25; allows the variational ansatz Eq. (3) without spin-flip terms.
  • standard math Allowed transitions in perturbation theory assume all degeneracies are lifted in disordered systems.
    Used in Appendix Eq. (A.4) to write the first-order wavefunction correction without degenerate perturbation theory; generic for disordered systems.
  • domain assumption RPA captures the singular behavior of the static susceptibility near Tc and the spin-wave poles.
    RPA is the standard approximation; no benchmark against exact results or Monte Carlo is provided in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Magnetic susceptibility of diluted magnetic semiconductors at low carrier densities." pith.science (2026). https://pith.science/paper/MHGXYTNW

@misc{pith2026260805620,
  author       = {Pith},
  title        = {Pith review of: Magnetic susceptibility of diluted magnetic semiconductors at low carrier densities},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MHGXYTNW}},
  note         = {Machine review of arXiv:2608.05620}
}
read the original abstract

We calculate the static longitudinal and the transverse dynamic magnetic susceptibilities of (III,Mn)V diluted magnetic semiconductors, using the random phase approximation, for a simple impurity band model appropriate for the low charge carrier concentration regime. The magnetic susceptibilities are shown to depend sensitively on the amount of positional disorder of the Mn impurities. The results we obtain are consistent with previous studies of the spin wave spectrum and of the spatially inhomogeneous ferromagnetic state of these materials.

Figures

Figures reproduced from arXiv: 2608.05620 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4: Left: total susceptibility of a disordered sample, n [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figures from the paper (4 more)
Figure 6
Figure 6. Figure 6: FIG. 6: Spin-wave energies ¯h [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p009_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    H. Ohno, J. Magn. Magn. Mat. 200 , 110 (1999)

    work page 1999

  2. [2]

    K. W. Edmonds, P. Bogusawski, K. Y. Wang, R. P. Campion, S. N. Novikov, N. R. S. Farley, B. L. Gallagher, C. T. Foxon, M. Sawicki, T. Dietl, M. Buongiorno Nardelli and J. Bernholc, Phys. Rev. Lett. 92 , 037201 (2004)

    work page 2004

  3. [3]

    A. M. Nazmul, S. Sugahara, and M. Tanaka, Phys. Rev. B 67 , 241308(R) (2003)

    work page 2003

  4. [4]

    Beschoten, P.A

    B. Beschoten, P.A. Crowell, I. Malajovich, D.D. Awschalom, F. Matsukura, A. Shen, and H. Ohno, Phys. Rev. Lett. 83 , 3073 (1999)

    work page 1999

  5. [5]

    K. M. Yu, W. Walukiewicz, T. Wojtowicz, I. Kuryliszyn, X. Liu, Y. Sasaki and J. K. Furdyna, Phys. Rev. B 65 , 201303(R) (2002)

    work page 2002

  6. [6]

    Zar\'and and B

    G. Zar\'and and B. Janko, Phys. Rev. Lett. 89 , 047201 (2002)

    work page 2002

  7. [7]

    G. A. Fiete, G. Zar\'and and K. Damle, Phys. Rev. Lett. 91 , 097202 (2003)

    work page 2003

  8. [8]

    Brey and G

    L. Brey and G. Gomez-Santos, Phys. Rev. B 68 , 115206 (2003)

    work page 2003

Show all 36 references
  1. [9]

    Zhou, M.P

    C. Zhou, M.P. Kennett, X. Wan, M. Berciu and R.N. Bhatt, Phys. Rev. B 69 , 144419 (2004)

    2004

  2. [10]

    Berciu and R

    M. Berciu and R. N. Bhatt, Phys. Rev. Lett. 87 , 107203 (2000)

    2000

  3. [11]

    Berciu and R

    M. Berciu and R. N. Bhatt, Phys. Rev. B 69 , 045202 (2004)

    2004

  4. [12]

    M. P. Kennett, M. Berciu and R. N. Bhatt, Phys. Rev. B 66 , 045207 (2002)

    2002

  5. [13]

    Berciu and R

    M. Berciu and R. N. Bhatt, Phys. Rev. B 66 , 085207 (2002)

    2002

  6. [14]

    Berciu and R

    M. Berciu and R. N. Bhatt, Phys. Rev. Lett. 90 , 029702 (2003)

    2003

  7. [15]

    M. P. Kennett, Mona Berciu and R. N. Bhatt, Phys. Rev. B 65 , 115308 (2002)

    2002

  8. [16]

    Theodoropoulos, A

    N. Theodoropoulos, A. F. Hebard, M. E. Overberg, C. R. Abernathy, S. J. Pearton, S. N. G. Chu and R. G. Wilson, Appl. Phys. Lett. 78 , 3475 (2001)

    2001

  9. [17]

    Y. D. Park, A. T. Hanbicki, S. C. Erwin, C. S. Hellberg, J. M. Sullivan, J. E. Mattson, T. F. Ambrose, A. Wilson, G. Spanos and B. T. Jonker, Science 295 , 651 (2002)

    2002

  10. [18]

    onig, H.-H. Lin and A. H. MacDonald, Phys. Rev. Lett. 84 , 5628 (2000); J. Schliemann, J. K\

    J. K\"onig, H.-H. Lin and A. H. MacDonald, Phys. Rev. Lett. 84 , 5628 (2000); J. Schliemann, J. K\"onig, H.-H. Lin and A. H. MacDonald, Appl. Phys. Lett. 78 , 1550 (2001)

    2000

  11. [19]

    Chattopadhyay, S

    A. Chattopadhyay, S. Das Sarma and A. J. Millis, Phys. Rev. Lett. 87 , 222702 (2001)

    2001

  12. [20]

    Dietl, H

    T. Dietl, H. Ohno and T. Matsukura, Phys. Rev. B 63 , 195205 (2001)

    2001

  13. [21]

    C. Timm, F. von Oppen, and F. H\"ofling, Phys. Rev. B 69 , 115202 (2004); Y. Qi and S. Zhang, Phys. Rev. B 67 , 052407 (2003)

    2004

  14. [22]

    R. N. Bhatt, Phys. Rev. B 24 , 3630 (1981); Phys. Rev. B 26 , 1082 (1982)

    1981

  15. [23]

    Bhattacharjee and C.B

    A.K. Bhattacharjee and C.B. \'a la Guillaume, Solid State Commun. 113 , 17 (2000)

    2000

  16. [24]

    C. Timm, F. Sch\"afer, and F. von Oppen, Phys. Rev. Lett. 89 , 137201 (2002)

    2002

  17. [25]

    However, a previous study of this model HF showed that the self-consistent mean-field ground state is always collinear, and thus the variational guess of Eq

    The most general variational form would allow spin flips _ i,j, ^ h^B_ i , j c^ _ i c_ j and non-collinearity _ i H ^B_ i S _i in K . However, a previous study of this model HF showed that the self-consistent mean-field ground state is always collinear, and thus the variationa...

  18. [26]

    Blaizot and G

    J.-P. Blaizot and G. Ripka, Quantum theory of finite systems (MIT Press, Cambridge, Mass., 1986)

    1986

  19. [27]

    Furdyna, J

    J. Furdyna, J. Appl. Phys. 64 , R29 (1988); S. Lee, M. Dobrowolska, J. K. Furdyna and L. R. Ram-Mohan, Phys. Rev. B 61 , 2120 (2000)

    1988

  20. [28]

    M. A. Zudov, J. Kono, Y. H. Matsuda, T. Ikaida, N. Miura, H. Munekata, G. D. Sanders, Y. Sun and C. J. Stanton, Phys. Rev. B 66 , 161307(R) (2002)

    2002

  21. [29]

    Berciu and B

    M. Berciu and B. Janko, Phys. Rev. Lett. 90 , 246804 (2003)

    2003

  22. [30]

    Berciu, T

    M. Berciu, T. Rappoport and B. Janko (unpublished)

  23. [31]

    Schliemann, J

    J. Schliemann, J. K\"onig and A.H. MacDonald, Phys. Rev. B 64 , 165201 (2001)

    2001

  24. [32]

    Kaminski and S

    A. Kaminski and S. Das Sarma, Phys. Rev. Lett. 88 , 247202 (2002); Phys. Rev. B 68 , 235210 (2003)

    2002

  25. [33]

    M. Mayr, G. Alvarez and E. Dagotto, Phys. Rev. B 65 , 241202 (2002)

    2002

  26. [34]

    Thesis, (University of British Columbia, 2004)

    Adel Kassaian, ``Magnetic susceptibility of diluted magnetic semiconductors'' , M.Sc. Thesis, (University of British Columbia, 2004)

    2004

  27. [35]

    Alvarez and E

    G. Alvarez and E. Dagotto, Phys. Rev. B 68 , 045202 (2003)

    2003

  28. [36]

    D. J. Priour, Jr., E. H. Hwang, and S. Das Sarma, Phys. Rev. Lett. 92 , 117201 (2004); S. Das Sarma, E. H. Hwang, and A. Kaminski, Phys. Rev. B 67 , 155201 (2003)

    2004

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.