GdFe₂ Laves phase intermetallic system under pressure: an ab-initio study
Pith reviewed 2026-05-24 18:02 UTC · model grok-4.3
The pith
Fe-Gd exchange pairs supply the missing 25% of experimental Curie temperature in GdFe2.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In GdFe2, mean-field solution of the Heisenberg model with DFT/LDA exchange integrals shows the Fe sublattice supplies about 75% of Tc^exp, the Gd sublattice supplies very little, and the Fe-Gd exchange pairs supply the missing 25%. The same calculation reproduces the experimental rise of Tc under hydrostatic pressure corresponding to 2% volume compression.
What carries the argument
Mean-field solution of the multi-sublattice Heisenberg model fed by DFT/LDA exchange integrals and local moments
If this is right
- Tc increases with pressure in agreement with experiment.
- Fe sublattice alone accounts for roughly 75% of Tc^exp.
- Gd sublattice contributes negligibly.
- Fe-Gd pairs close the remaining 25% gap.
Where Pith is reading between the lines
- The same decomposition into sublattice and cross terms may be needed in other rare-earth transition-metal Laves phases where transition-metal-only estimates fall short.
- Because the pressure range studied is only 2% volume change, the relative importance of Fe-Gd pairs is likely to persist at ambient conditions in related compounds.
- Applying the same mean-field decomposition to systems with different rare-earth or transition-metal partners could test whether the 25% cross-term fraction is general.
Load-bearing premise
The mean-field Heisenberg model with DFT/LDA exchange integrals reproduces the experimental Curie temperature once contributions from Fe, Gd, and Fe-Gd pairs are all included.
What would settle it
A measured Curie temperature under pressure that deviates from the calculated increase, or an independent evaluation of the Fe-Gd exchange integral showing its contribution differs from 25%.
read the original abstract
Here we perform $ab-initio$ study of Curie temperature $T_C$ under hydrostatic pressure for intermetallic compound GdFe$_2$. To calculate $T_C$ for GdFe$_2$ we applied mean-field solution of the Heisenberg model for several magnetic sublattices with DFT/LDA calculated values of necessary exchange interaction integrals and local magnetic moments. To compare with available experimental data pressure values were taken from zero up to about 70 Kbar. It corresponds to 2\% compression of the volume of the unit cell. In agreement with experimental data $T_C$ grows under pressure. It was shown that Fe ions magnetic sublattice alone provides only about 75\% of the experimental Curie temperature $T_C^{exp}$. Gd sublattice is found to give very weak contribution to the $T_C^{exp}$. Here we show that the missing 25 \% of $T_C^{exp}$ comes from Fe-Gd exchange pairs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports an ab-initio DFT/LDA study of the GdFe2 Laves phase, computing exchange integrals and local moments, then solving the multi-sublattice Heisenberg model in mean-field approximation to obtain the pressure dependence of Tc up to ~70 kbar (2% volume compression). It finds Tc increasing with pressure in agreement with experiment and decomposes the contributions to conclude that the Fe sublattice supplies ~75% of experimental Tc, the Gd sublattice contributes negligibly, and Fe-Gd pairs account for the remaining 25%.
Significance. If the numerical decomposition is reliable, the work supplies a concrete attribution of inter-sublattice exchange contributions to Tc in an RFe2 compound under pressure, which could inform models of magnetovolume effects in related intermetallics.
major comments (2)
- [Abstract] Abstract: the claim that Fe-Gd pairs supply the missing 25% of T_C^exp is obtained from the largest eigenvalue of the full J matrix in the multi-sublattice mean-field solution. This decomposition is directly sensitive to the accuracy of the LDA-computed J_Fe-Gd integrals; the manuscript provides no test of LDA versus LDA+U (or other 4f treatments) for the localized Gd 4f states, leaving open the possibility that delocalization errors in LDA systematically misestimate the inter-sublattice terms and thereby falsify the 75%/25% split.
- [Heisenberg model solution] The section describing the Heisenberg model solution: no validation or error estimate is given for the mean-field approximation itself, nor are details supplied on how the exchange integrals were extracted (e.g., Liechtenstein formula parameters) or on their sensitivity; without these, the precise numerical attribution to Fe-Gd pairs cannot be assessed as robust.
minor comments (1)
- The pressure range is stated as 'up to about 70 Kbar' corresponding to 2% compression; a table or explicit equation linking the two would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We provide point-by-point responses below.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that Fe-Gd pairs supply the missing 25% of T_C^exp is obtained from the largest eigenvalue of the full J matrix in the multi-sublattice mean-field solution. This decomposition is directly sensitive to the accuracy of the LDA-computed J_Fe-Gd integrals; the manuscript provides no test of LDA versus LDA+U (or other 4f treatments) for the localized Gd 4f states, leaving open the possibility that delocalization errors in LDA systematically misestimate the inter-sublattice terms and thereby falsify the 75%/25% split.
Authors: We agree that the 75%/25% attribution depends on the LDA J integrals. LDA is the standard choice in the literature for RFe2 compounds and yields pressure trends consistent with experiment. We will add a paragraph discussing the choice of LDA for Gd 4f states and the possible influence of LDA+U, while noting that a full comparative study lies beyond the present scope. revision: partial
-
Referee: [Heisenberg model solution] The section describing the Heisenberg model solution: no validation or error estimate is given for the mean-field approximation itself, nor are details supplied on how the exchange integrals were extracted (e.g., Liechtenstein formula parameters) or on their sensitivity; without these, the precise numerical attribution to Fe-Gd pairs cannot be assessed as robust.
Authors: We will expand the methods section to specify the Liechtenstein formula implementation details (cutoff radii, k-mesh, etc.) and include a short discussion of MFA limitations together with a note that the pressure dependence of Tc matches experiment. These additions will allow readers to assess the robustness of the decomposition. revision: yes
Circularity Check
No circularity: Tc and pair contributions obtained from independent DFT J integrals via mean-field Heisenberg model
full rationale
The derivation computes exchange integrals from DFT/LDA (Liechtenstein formula implied), inserts them into the multi-sublattice mean-field Heisenberg model, and decomposes Tc into Fe-Fe (~75% of Tc^exp), Gd-Gd (weak), and Fe-Gd (~25%) contributions by comparing partial versus full calculations. These steps use first-principles inputs without fitting parameters to the target Tc^exp or renaming fitted quantities as predictions. No self-citation is load-bearing, no ansatz is smuggled, and the chain does not reduce to its own outputs by construction. The paper remains self-contained against external DFT and mean-field benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Mean-field approximation to the Heisenberg model yields a reliable estimate of Tc for this multi-sublattice intermetallic.
- domain assumption DFT/LDA produces sufficiently accurate local moments and exchange integrals for the pressure range studied.
Reference graph
Works this paper leans on
-
[1]
Wallace, ed., Rare Earth Intermetallics (Aca- demic Press, 1973)
W . Wallace, ed., Rare Earth Intermetallics (Aca- demic Press, 1973)
work page 1973
-
[2]
Buschow, Reports on Progress in Physics 54, 1123 (1991)
K.H.J. Buschow, Reports on Progress in Physics 54, 1123 (1991)
work page 1991
-
[3]
A.V . Lukoyanov, E.E. Kokorina, M.V . Medvedev, I.A. Nekrasov, Phys. Rev. B 80, 104409 (2009)
work page 2009
-
[4]
P .A. Igoshev, E.E. Kokorina, M.V . Medvedev, I.A. Nekrasov, Journal of Experimental and Theoretical Physics 127, 342 (2018)
work page 2018
-
[5]
Tesluk, ed., Metallic compounds with Laves phase (Moskow: Nauka, 1969)
M.Y . Tesluk, ed., Metallic compounds with Laves phase (Moskow: Nauka, 1969)
work page 1969
-
[6]
C. D. Creagh, S. A yling, Journal of Materials Science 13, 113 (1978)
work page 1978
-
[7]
V .I. Anisimov, F. Aryasetiawan, A.I. Lichtenstein, Journal of Physics: Condensed Matter 9, 767 (1997)
work page 1997
- [8]
- [9]
- [10]
-
[11]
Murnaghan, Proceedings of the National Academy of Sciences 30, 244 (1944)
F.D. Murnaghan, Proceedings of the National Academy of Sciences 30, 244 (1944)
work page 1944
-
[12]
T. Palasyuk, H. Figiel, M. Tkacz, Journal of Alloys and Compounds 375, 62 (2004)
work page 2004
-
[13]
G.M. Elalfy, R.M. Shabara, S.H. Aly, S. Y ehia, Com- putational Condensed Matter 5, 24 (2015)
work page 2015
- [14]
- [15]
- [16]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.