Doping-dependent orbital magnetism in Chromium pnictides
Pith reviewed 2026-05-21 15:44 UTC · model grok-4.3
The pith
Electron doping in LaCrAsO switches the system between commensurate antiferromagnetic states tied to localized d3z2-r2 electrons and incommensurate states driven by itinerant dxy electrons.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using the matrix random-phase approximation on a tight-binding or DFT-derived band structure, the work shows that electron doping in LaCrAsO first stabilizes a commensurate antiferromagnetic state, then a stripe antiferromagnetic state; at higher doping both reappear with incommensurate magnetic vectors. The commensurate instabilities are carried by more localized electrons in the Cr d3z2-r2 orbital, whereas the incommensurate instabilities are carried by the dxy orbital whose stronger overlap favors itinerant-electron magnetism.
What carries the argument
Matrix random-phase approximation applied to the multi-orbital band structure, which computes the doping-dependent magnetic susceptibility and identifies the leading orbital-resolved instabilities.
If this is right
- At low electron doping the ground state is a commensurate antiferromagnet with opposite spins on the two Cr sublattices.
- Intermediate doping stabilizes a stripe antiferromagnetic phase.
- Still higher doping restores both phases but with incommensurate ordering vectors.
- The switch from commensurate to incommensurate order tracks a change from localized d3z2-r2 to itinerant dxy orbital character.
Where Pith is reading between the lines
- The same orbital-selectivity mechanism may operate in other chromium pnictides when carriers are added by substitution or pressure.
- If the dxy-driven incommensurate state can be stabilized near a superconducting dome, it would suggest a route to test whether itinerant magnetism competes with or promotes pairing.
- The sequence implies that modest further doping or strain could be used to toggle between localized and itinerant magnetism in a single material.
Load-bearing premise
The chosen tight-binding or DFT band structure together with a fixed set of interaction parameters is enough to determine the sequence of magnetic instabilities without vertex corrections or exact treatment of competing orders.
What would settle it
Neutron scattering measurements that track the evolution of magnetic ordering vectors with electron doping, combined with orbital-sensitive spectroscopy that confirms whether d3z2-r2 or dxy character dominates at each doping level.
Figures
read the original abstract
We present results for the phase diagram of the parent compound LaCrAsO under electron doping using the matrix random-phase approximation. At low doping levels, the system stabilizes an antiferromagnetic state in which different Cr sublattices carry opposite spins, consistent with experimental observations. As the doping concentration increases, a stripe-type antiferromagnetic phase becomes favored. At even higher doping, the system repeats the two former magnetic states, but with incommensurate magnetic ordering vectors. The commensurate magnetic phases are associated with more localized electrons in the Cr $d_{3z^2-r^2}$ orbital, whereas the incommensurate phases are linked to the $d_{xy}$ orbital, whose stronger overlap favors itinerant-electron magnetism.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript applies the matrix random-phase approximation to a multi-orbital model of electron-doped LaCrAsO to map the evolution of magnetic instabilities with doping. It reports an antiferromagnetic state with opposite spins on Cr sublattices at low doping (consistent with experiment), a transition to stripe-type antiferromagnetism at intermediate doping, and repetition of both phases with incommensurate ordering vectors at higher doping. The commensurate instabilities are attributed to more localized electrons in the Cr d_{3z^2-r^2} orbital, while the incommensurate ones are linked to the d_{xy} orbital whose stronger overlap favors itinerant magnetism.
Significance. If the orbital-resolved assignment of instabilities survives scrutiny, the work would provide a concrete microscopic picture of how orbital character selects between localized and itinerant magnetism in Cr pnictides. The matrix RPA approach permits direct extraction of orbital weights in the leading eigenvectors, which is a useful technical feature for multi-orbital systems when the underlying band structure and interactions are adequately justified.
major comments (2)
- [Methods] Methods section (band-structure and interaction parameters): the central claim that d_{3z^2-r^2} drives commensurate order while d_{xy} drives incommensurate order rests on the orbital character of the RPA susceptibility eigenvectors. No explicit values for the interaction parameters (U, J, etc.), no source or convergence details for the tight-binding/DFT band structure, and no benchmark against experimental ordering vectors or other observables are provided. These inputs directly control which orbital dominates the leading instability at each doping, so their justification is load-bearing for the reported sequence.
- [Results] Results on doping evolution: the transition from commensurate to incommensurate order is obtained from the poles of the RPA susceptibility. The manuscript does not examine whether vertex corrections or self-energy renormalizations (known to be relevant for intermediate-correlation Cr d states) would shift the crossover doping or alter the orbital content of the leading eigenvector. A direct comparison with a method that includes these effects would be required to confirm that the reported orbital association is robust.
minor comments (2)
- [Abstract] The abstract states that the low-doping AF state is 'consistent with experimental observations' but does not cite the specific experimental ordering vector or reference; adding this citation would strengthen the comparison.
- [Figure captions] Figure captions for the phase diagram should explicitly label the doping intervals corresponding to each magnetic phase to improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We respond to each major comment below, indicating the revisions we will implement where appropriate.
read point-by-point responses
-
Referee: [Methods] Methods section (band-structure and interaction parameters): the central claim that d_{3z^2-r^2} drives commensurate order while d_{xy} drives incommensurate order rests on the orbital character of the RPA susceptibility eigenvectors. No explicit values for the interaction parameters (U, J, etc.), no source or convergence details for the tight-binding/DFT band structure, and no benchmark against experimental ordering vectors or other observables are provided. These inputs directly control which orbital dominates the leading instability at each doping, so their justification is load-bearing for the reported sequence.
Authors: We agree that explicit documentation of the model parameters and their provenance is necessary to substantiate the orbital assignments. In the revised version we will add a dedicated Methods subsection that lists the interaction parameters (U, J_H, U', J') together with their origin in constrained-RPA estimates for Cr pnictides, specifies the DFT functional and k-mesh used to generate the tight-binding Hamiltonian, and reports convergence checks. We will also include a direct comparison of the calculated zero-doping ordering vector with the experimental neutron-scattering result for LaCrAsO. revision: yes
-
Referee: [Results] Results on doping evolution: the transition from commensurate to incommensurate order is obtained from the poles of the RPA susceptibility. The manuscript does not examine whether vertex corrections or self-energy renormalizations (known to be relevant for intermediate-correlation Cr d states) would shift the crossover doping or alter the orbital content of the leading eigenvector. A direct comparison with a method that includes these effects would be required to confirm that the reported orbital association is robust.
Authors: We acknowledge that the RPA neglects vertex corrections and self-energy renormalizations, which can be quantitatively important for Cr d states. Within the present methodological scope, however, the matrix RPA already provides orbital-resolved eigenvectors that directly link the leading instability to a given orbital character; this diagnostic is unavailable in many other approximations. We will add a concise limitations paragraph in the Discussion that states the RPA approximation explicitly and notes that a systematic study including vertex corrections would require a different framework (e.g., fRG or DMFT+RPA) and is left for future work. No change to the reported RPA phase sequence is planned. revision: partial
Circularity Check
RPA-derived orbital associations are direct outputs of susceptibility eigenvectors, not definitional reductions
full rationale
The derivation proceeds by constructing a multi-orbital tight-binding or DFT band structure, adding Hubbard interactions, and computing the matrix RPA susceptibility to locate doping-dependent poles that signal magnetic instabilities. The orbital assignments (commensurate phases linked to d_{3z^2-r^2} localization, incommensurate to d_{xy} itinerancy) are read off from the orbital-resolved eigenvectors at those poles. This is a standard forward calculation whose outputs are not equivalent to the input band parameters or interaction strengths by construction. No self-citation chain, fitted parameter renamed as prediction, or ansatz smuggled via prior work is required for the central claim. The paper remains self-contained against external benchmarks such as experimental ordering vectors and doping trends.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Matrix random-phase approximation accurately identifies the leading magnetic instability as a function of doping.
Reference graph
Works this paper leans on
-
[1]
P. A. Lee, N. Nagaosa, X.-G. Wen, Doping a Mott insulator: Physics of high-temperature superconductivity, Rev. Mod. Phys. 78 (2006) 17.doi: 10.1103/RevModPhys.78.17
-
[2]
N. P. Armitage, P. Fournier, R. L. Greene, Progress and perspectives on electron-doped cuprates, Rev. Mod. Phys. 82 (2010) 2421.doi:10.1103/ RevModPhys.82.2421
work page 2010
- [3]
-
[4]
Note that in this case the cuprates are being hole-doped and the Fe-SC are being electron-doped. 14 □π 0 π a) dx2□y2 Leading λ =0.000243 b) dxy Leading λ =0.000308 □π 0 π□π 0 π c) dxy Leading λ =0.000366 □π 0 π d) gxy(x2□y2) Leading λ =0.000329 □0.10 □0.05 0.0 0.05 0.10 □0.10 □0.05 0.0 0.05 0.10 □0.10 □0.05 0.0 0.05 0.10 □0.10 □0.05 0.0 0.05 0.10 Figure A...
-
[5]
J. M. Pizarro, M. J. Calderón, J. Liu, M. C. Muñoz, E. Bascones, Strong correlations and the search for high-Tc superconductivity in chromium pnic- tides and chalcogenides, Phys. Rev. B 95 (2017) 075115.doi:10.1103/ PhysRevB.95.075115
work page 2017
-
[6]
M. Edelmann, G. Sangiovanni, M. Capone, L. de’ Medici, Chromium analogs of iron-based superconductors, Phys. Rev. B 95 (2017) 205118. doi:10.1103/PhysRevB.95.205118
-
[7]
W.-S. Wang, M. Gao, Y . Yang, Y .-Y . Xiang, Q.-H. Wang, Possible supercon- ductivity in the electron-doped chromium pnictide LaOCrAs, Phys. Rev. B 95 (2017) 144507.doi:10.1103/PhysRevB.95.144507
-
[8]
Note that, as Cu-doping in NaFe 1−xCuxAs has been shown to result in elec- tron doping of the FeAs planes [36], changingxfrom∼0.5 to 0, as done in Ref. [37], is an example of hole-doping a Fe-SC Mott-insulating phase and obtaining a SC phase with aT max c ∼11.5 K atx=0.019 [36], which is not as high as the above reasoning would suggest. However, one may a...
-
[9]
Y . Singh, A. Ellern, D. C. Johnston, Physical Review B 79 (9) (2009) 094519.doi:10.1103/PhysRevB.79.094519, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevB.79.094519
-
[10]
Y . Singh, M. A. Green, Q. Huang, A. Kreyssig, R. J. McQueeney, D. C. Johnston, A. I. Goldman, Physical Review B 80 (10) (2009) 100403.doi: 10.1103/PhysRevB.80.100403, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevB.80.100403
-
[11]
B. S. Jacobs, A. Pandey, bamn 2p2: Highest magnetic ordering tem- perature 122-pnictide compound, Phys. Rev. Mater. 7 (2023) 044410. doi:10.1103/PhysRevMaterials.7.044410. URLhttps://link.aps.org/doi/10.1103/PhysRevMaterials.7. 044410
-
[12]
J.-K. Bao, H. Jiang, Y .-L. Sun, W.-H. Jiao, C.-Y . Shen, H.-J. Guo, Y . Chen, C.-M. Feng, H.-Q. Yuan, Z.-A. Xu, G.-H. Cao, R. Sasaki, T. Tanaka, K. Matsubayashi, Y . Uwatoko, Weakly ferromagnetic metallic state in heav- ily doped ba 1−xkxmn2as2, Phys. Rev. B 85 (2012) 144523.doi:10.1103/ 16 PhysRevB.85.144523. URLhttp://link.aps.org/doi/10.1103/PhysRevB....
-
[13]
A. Pandey, R. S. Dhaka, J. Lamsal, Y . Lee, V . K. Anand, A. Kreyssig, T. W. Heitmann, R. J. McQueeney, A. I. Goldman, B. N. Harmon, A. Kamin- ski, D. C. Johnston, Physical Review Letters 108 (8) (2012) 087005. doi:10.1103/PhysRevLett.108.087005, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevLett.108. 087005
-
[14]
N. S. Sangeetha, V . Smetana, A.-V . Mudring, D. C. Johnston, Antiferromag- netism in semiconducting srmn2sb2 and bamn2sb2 single crystals, Phys. Rev. B 97 (2018) 014402.doi:10.1103/PhysRevB.97.014402. URLhttps://link.aps.org/doi/10.1103/PhysRevB.97.014402
-
[15]
A. T. Satya, A. Mani, A. Arulraj, N. V . C. Shekar, K. Vinod, C. S. Sundar, A. Bharathi, Physical Review B 84 (18) (2011) 180515.doi:10.1103/ PhysRevB.84.180515, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevB.84.180515
-
[16]
A. Pandey, V . K. Anand, D. C. Johnston, Physical Review B 84 (1) (2011) 014405.doi:10.1103/PhysRevB.84.014405, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevB.84.014405
-
[17]
J. Lamsal, G. S. Tucker, T. W. Heitmann, A. Kreyssig, A. Jesche, A. Pandey, W. Tian, R. J. McQueeney, D. C. Johnston, A. I. Goldman, Physical Review B 87 (14) (2013) 144418.doi:10.1103/PhysRevB.87.144418, [link]. URLhttp://link.aps.org/doi/10.1103/PhysRevB.87.144418
-
[18]
D. J. Singh, A. S. Sefat, M. A. McGuire, B. C. Sales, D. Mandrus, L. H. Van- Bebber, V . Keppens, Itinerant antiferromagnetism in bacras2: Experimental characterization and electronic structure, Physical Review B 79 (9) (2009) 094429.doi:10.1103/PhysRevB.79.094429. URLhttp://link.aps.org/doi/10.1103/PhysRevB.79.094429
-
[19]
Z. Ning, P. Das, Y . Lee, N. S. Sangeetha, D. L. Abernathy, D. C. John- ston, R. J. McQueeney, D. Vaknin, L. Ke, Spin dynamics in the itinerant antiferromagnet srcr2as2, Phys. Rev. B 111 (2025) 134438.doi:10.1103/ PhysRevB.111.134438. URLhttps://link.aps.org/doi/10.1103/PhysRevB.111.134438 17
-
[20]
U. B. Paramanik, R. Prasad, C. Geibel, Z. Hossain, Itinerant and local- moment magnetism in eucr2 as2 single crystals, Physical Review B - Con- densed Matter and Materials Physics 89 (14) (2014) 1–8.doi:10.1103/ PhysRevB.89.144423
work page 2014
-
[21]
S. Nandi, Y . Xiao, N. Qureshi, U. B. Paramanik, W. T. Jin, Y . Su, B. Oulad- diaf, Z. Hossain, T. Brückel, Magnetic structures of the eu and cr moments ineucr 2as2: Neutron diffraction study, Phys. Rev. B 94 (2016) 094411. doi:10.1103/PhysRevB.94.094411. URLhttp://link.aps.org/doi/10.1103/PhysRevB.94.094411
-
[22]
K. Kubo, Pairing symmetry in a two-orbital hubbard model on a square lattice, Phys. Rev. B 75 (2007) 224509.doi:10.1103/PhysRevB.75. 224509. URLhttps://link.aps.org/doi/10.1103/PhysRevB.75.224509
- [23]
-
[24]
R. M. Fernandes, A. V . Chubukov, Low-energy microscopic models for iron-based superconductors: a review, Rep. Prog. Phys. 80 (2016) 014503. doi:10.1088/1361-6633/80/1/014503. URLhttps://iopscience.iop.org/article/10.1088/ 1361-6633/80/1/014503/meta?casa_token=XjO8PnhuG4oAAAAA: gofMFXJjs9D7SwV90rTULD2loWli2wwOS8hchbMWXOqwjPYXp1iJtesOJMSPzmsFPBZ0AlRiuO4ENr...
-
[25]
M. Altmeyer, D. Guterding, P. J. Hirschfeld, T. A. Maier, R. Valentí, D. J. Scalapino, Role of vertex corrections in the matrix formulation of the ran- dom phase approximation for the multiorbital Hubbard model, Phy. Rev. B 94 (2016) 214515.doi:10.1103/PhysRevB.94.214515. URLhttps://doi.org/10.1103/PhysRevB.94.214515
-
[26]
T. Ozaki, Continued fraction representation of the fermi-dirac function for large-scale electronic structure calculations, Phys. Rev. B 75 (2007) 035123. doi:10.1103/PhysRevB.75.035123. URLhttps://link.aps.org/doi/10.1103/PhysRevB.75.035123 18
-
[27]
L. V . Boehnke, Susceptibilities in materials with multiple strongly correlated orbitals, Ph.D. thesis, Universität Hamburg, dissertation, oral defence: 2014- 10-28 (2015). URLhttps://ediss.sub.uni-hamburg.de/handle/ediss/6373
work page 2014
-
[28]
M. R. Cantarino, K. R. Pakuszewski, B. Salzmann, P. H. A. Moya, W. R. d. S. Neto, G. S. Freitas, P. G. Pagliuso, W. H. Brito, C. Monney, C. Adriano, F. A. Garcia, Incoherent electronic band states in mn-substituted bafe 2as2, Phys. Rev. B 108 (2023) 245124.doi:10.1103/PhysRevB.108.245124. URLhttps://link.aps.org/doi/10.1103/PhysRevB.108.245124
-
[29]
Heavy operators and hydrodynamic tails
M. dos Reis Cantarino, K. R. Pakuszewski, B. Salzmann, P. H. A. Moya, W. R. da Silva Neto, G. S. Freitas, P. G. Pagliuso, C. Adriano, W. H. Brito, F. A. Garcia, Hole doping and electronic correlations in Cr substituted BaFe2As2, SciPost Phys. 17 (2024) 141.doi:10.21468/SciPostPhys. 17.5.141. URLhttps://scipost.org/10.21468/SciPostPhys.17.5.141
-
[30]
X. Chen, J. Qi, D. Shi, Strain-engineering of magnetic coupling in two- dimensional magnetic semiconductor crsite 3: Competition of direct ex- change interaction and superexchange interaction, Phys. Lett. A 379 (2015) 60–63.doi:10.1016/j.physleta.2014.10.042
-
[31]
Y . Lou, D. Wu, Y . Liao, C. Fang, Y . Pu, J. Xie, A in-plane biaxial strain tunable electronic structures and magnetic properties of fe2c monolayer, J. Magn. Magn. Mater. 563 (2022) 169959.doi:10.1016/j.jmmm.2022. 169959
-
[32]
E. Bascones, B. Valenzuela, M. J. Calderón, Orbital differentiation and the role of orbital ordering in the magnetic state of fe superconductors, Phys. Rev. B 86 (2012) 174508.doi:10.1103/PhysRevB.86.174508
-
[33]
J. Nayak, K. Filsinger, G. H. Fecher, S. Chadov, J. Minár, E. D. L. Rienks, B. Büchner, S. P. Parkin, J. Fink, C. Felser, Observation of a remarkable reduction of correlation effects in bacr2as2 by arpes, Proc. Natl. Acad. Sci. U.S.A. 114 (47) (2017) 12425–12429.doi:10.1073/pnas.1702234114
-
[34]
P. Richard, A. van Roekeghem, B. Q. Lv, T. Qian, T. K. Kim, M. Hoesch, J.-P. Hu, A. S. Sefat, S. Biermann, H. Ding, IsBaCr 2As2 symmetrical to 19 BaFe2As2 with respect to half 3dshell filling?, Phys. Rev. B 95 (2017) 184516.doi:10.1103/PhysRevB.95.184516
-
[35]
A. G. de Figueiredo, M. R. Cantarino, W. R. da Silva Neto, K. R. Pakuszewski, R. Grossi, D. S. Christovam, J. C. Souza, M. M. Piva, G. S. Freitas, P. G. Pagliuso, C. Adriano, F. A. Garcia, Orbital localization and the role of the fe and as 4porbitals in bafe 2as2 probed by xanes, Phys. Rev. B 105 (2022) 045130.doi:10.1103/PhysRevB.105.045130. URLhttps://l...
-
[36]
A. F. Wang, J. J. Lin, P. Cheng, G. J. Ye, F. Chen, J. Q. Ma, X. F. Lu, B. Lei, X. G. Luo, X. H. Chen, Phase diagram and physical properties of nafe1−xcuxas single crystals, Phys. Rev. B 88 (2013) 094516.doi:10.1103/ PhysRevB.88.094516
work page 2013
-
[37]
Y . Song, Z. Yamani, C. Cao, Y . Li, C. Zhang, J. S. Chen, Q. Huang, H. Wu, J. Tao, Y . Zhu, W. Tian, S. Chi, H. Cao, Y .-B. Huang, M. Dantz, T. Schmitt, R. Yu, A. H. Nevidomskyy, E. Morosan, Q. Si, P. Dai, A Mott insulator continuously connected to iron pnictide superconductors, Nat. Commun. 7 (2016).doi:10.1038/ncomms13879. 20
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.