pith. sign in

arxiv: 1907.09343 · v1 · pith:S6AJSLVMnew · submitted 2019-07-22 · ❄️ cond-mat.mtrl-sci

GdFe₂ Laves phase intermetallic system under pressure: an ab-initio study

Pith reviewed 2026-05-24 18:02 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords GdFe2Curie temperaturehydrostatic pressureab initioHeisenberg modelexchange interactionsLaves phase
0
0 comments X

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.

The authors calculate the pressure dependence of the Curie temperature in GdFe2 using first-principles DFT/LDA methods to obtain exchange integrals, then insert those values into a mean-field Heisenberg model that treats multiple magnetic sublattices separately. They find that the iron sublattice alone reaches only about 75% of the measured Tc while the gadolinium sublattice contributes almost nothing, so the iron-gadolinium cross terms close the remaining gap. The calculated Tc increases with hydrostatic pressure up to 70 kbar in line with experiment, showing that the inter-sublattice exchanges are essential for a quantitative account.

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

These are editorial extensions of the paper, not claims the author makes directly.

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

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We provide point-by-point responses below.

read point-by-point responses
  1. 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

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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of LDA exchange integrals and the validity of the mean-field Heisenberg solution for a multi-sublattice system; no free parameters are explicitly fitted to Tc in the abstract, but LDA itself contains standard approximations.

axioms (2)
  • domain assumption Mean-field approximation to the Heisenberg model yields a reliable estimate of Tc for this multi-sublattice intermetallic.
    Invoked to convert exchange integrals into Tc; stated in the abstract as the method used.
  • domain assumption DFT/LDA produces sufficiently accurate local moments and exchange integrals for the pressure range studied.
    Basis for all numerical inputs; no validation against higher-level methods shown in abstract.

pith-pipeline@v0.9.0 · 5706 in / 1273 out tokens · 38771 ms · 2026-05-24T18:02:25.314716+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    Wallace, ed., Rare Earth Intermetallics (Aca- demic Press, 1973)

    W . Wallace, ed., Rare Earth Intermetallics (Aca- demic Press, 1973)

  2. [2]

    Buschow, Reports on Progress in Physics 54, 1123 (1991)

    K.H.J. Buschow, Reports on Progress in Physics 54, 1123 (1991)

  3. [3]

    Lukoyanov, E.E

    A.V . Lukoyanov, E.E. Kokorina, M.V . Medvedev, I.A. Nekrasov, Phys. Rev. B 80, 104409 (2009)

  4. [4]

    Igoshev, E.E

    P .A. Igoshev, E.E. Kokorina, M.V . Medvedev, I.A. Nekrasov, Journal of Experimental and Theoretical Physics 127, 342 (2018)

  5. [5]

    Tesluk, ed., Metallic compounds with Laves phase (Moskow: Nauka, 1969)

    M.Y . Tesluk, ed., Metallic compounds with Laves phase (Moskow: Nauka, 1969)

  6. [6]

    C. D. Creagh, S. A yling, Journal of Materials Science 13, 113 (1978)

  7. [7]

    Anisimov, F

    V .I. Anisimov, F. Aryasetiawan, A.I. Lichtenstein, Journal of Physics: Condensed Matter 9, 767 (1997)

  8. [8]

    Andersen, Phys

    O.K. Andersen, Phys. Rev. B 12, 3060 (1975)

  9. [9]

    Gunnarsson, O

    O. Gunnarsson, O. Jepsen, O.K. Andersen, Phys. Rev. B 27, 7144 (1983)

  10. [10]

    Andersen, O

    O.K. Andersen, O. Jepsen, Phys. Rev. Lett. 53, 2571 (1984)

  11. [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)

  12. [12]

    Palasyuk, H

    T. Palasyuk, H. Figiel, M. Tkacz, Journal of Alloys and Compounds 375, 62 (2004)

  13. [13]

    Elalfy, R.M

    G.M. Elalfy, R.M. Shabara, S.H. Aly, S. Y ehia, Com- putational Condensed Matter 5, 24 (2015)

  14. [14]

    Brouha, K

    M. Brouha, K. Buschow, Journal of Applied Physics 44, 1813 (1973)

  15. [15]

    Néel, Ann

    L. Néel, Ann. de Phys. 3, 137 (1948)

  16. [16]

    Brouers, A

    F. Brouers, A. V edyayev, Solid State Communica- tions 9, 1521 (1971)