TbFe₂D₄.2 Deuteride: Structure & Transitions

• TbFe₂D₄.2 deuteride is an intermetallic hydride featuring a rare earth–transition metal framework stabilized by deuterium in interstitial sites.
• Neutron and synchrotron diffraction reveal a reversible order–disorder transition from a monoclinic to cubic phase between 320–380 K alongside multipeak thermal desorption from 400–550 K.
• Precise D-atom occupancy and systematic phase evolution, including monoclinic distortions and tetragonal superstructures, make this compound a model for studying hydrogen storage and phase diagram behavior.

TbFe$_{2}$D$_{4.2}$ deuteride is an intermetallic hydride featuring a rare earth–transition metal backbone (TbFe$_2$) stabilized by deuterium occupancy at interstitial sites. Its intricate phase behavior and structural transitions, particularly under variable thermodynamic conditions, place this system among the model compounds for studying order–disorder transformations and polymorphism in Laves-phase deuterides. The compound exhibits magnetic, crystallographic, and hydrogen–deuterium storage properties whose detailed understanding is relevant for both fundamental solid-state physics and materials engineering.

1. Room-Temperature Crystal Structure

At 300 K, TbFe$_{2}$D$_{4.2}$ crystallizes in a monoclinic structure (space group $Pc$, No. 7), derived from the MgCu$_2$-type (C15) cubic parent phase by long-range ordering of D atoms. The refined lattice parameters (as determined by neutron powder diffraction) are:

• $a = 5.5171(2)$ Å
• $b = 11.5061(5)$ Å
• $c = 9.4382(3)$ Å
• $_{4.2}$0
• $_{4.2}$1 Å$_{4.2}$2 per unit cell (2 formula units, $_{4.2}$3 Å$_{4.2}$4/f.u.)

The structure features 4 distinct Tb sites (each fully occupied on Wyckoff 2a), 8 unique Fe sites (on 4c, all fully occupied), and 18 tetrahedral interstitial D sites with varying occupancies:

• 15 of the [Tb$_{4.2}$5Fe$_{4.2}$6] type (occupancy $_{4.2}$7 to $_{4.2}$8)
• 3 of the [TbFe$_{4.2}$9] type (occupancy $_2$0 to $_2$1) The total D content is $_2$2 atoms per formula unit. D-atom occupancies and atomic positions are precisely determined, confirming significant long-range ordering driven by deuterium arrangement.

2. Order–Disorder Transition (320–380 K)

Upon heating, TbFe$_2$3D$_2$4 displays a two-step reversible order–disorder (OD) transition, probed by high-resolution in-situ XRD and NPD alongside DSC. Calorimetry (5 K/min) reveals exothermic transitions at $_2$5 K and $_2$6 K (and endothermic on cooling at $_2$7 K and $_2$8 K), each contributing $_2$9–$_{2}$0 J g$_{2}$1 to the transition enthalpy ($_{2}$2–$_{2}$3 J g$_{2}$4 or $_{2}$5–$_{2}$6 kJ mol$_{2}$7 D), indicating two first-order phase transformations.

The high-temperature disordered phase adopts the cubic MgCu$_{2}$8-type structure ($_{2}$9, No. 227), with D atoms occupying [Tb$_{4.2}$0Fe$_{4.2}$1] tetrahedral sites in a statistically disordered fashion—superstructure reflections due to D order disappear entirely. Two-phase coexistence (monoclinic + cubic) spans $_{4.2}$2–$_{4.2}$3 K (from NPD), and above $_{4.2}$4 K, the structure is fully cubic and disordered.

3. Thermal Desorption and High-Temperature Phase Evolution

Between $_{4.2}$5 and $_{4.2}$6 K, thermal desorption proceeds via a multipeak sequence, reflecting successive phase transitions between cubic deuterides of varying D content. DSC registers broad exothermic features with principal maxima at $_{4.2}$7 K, $_{4.2}$8 K, and $_{4.2}$9 K, and an integrated desorption heat of $Pc$0 J g$Pc$1. In-situ NPD links these calorimetric events to desorption-rate maxima at $Pc$2 K, $Pc$3 K, $Pc$4 K, and $Pc$5 K.

These transitions are interpreted as a succession of distinct cubic phases (labeled $Pc$6, with $Pc$7 as the nearly D-free $Pc$8-phase of TbFe$Pc$9). Two-phase plateaus appear between each transition (e.g., $_2$0 between $_2$1–$_2$2 K, $_2$3 at $_2$4–$_2$5 K), signaling first-order character and two-phase coexistence. No activation energies are reported; enthalpic data derive from DSC integration.

4. Structural Sequence and Phase Diagram Versus Deuterium Content

Systematic ex-situ synchrotron XRD has been performed on partially desorbed TbFe$_2$6D$_2$7 samples spanning $_2$8 to $_2$9. The sequence of phases at 300 K as a function of D content is:

$a = 5.5171(2)$0 (D/f.u.)	Observed Phase(s)	Symmetry/Structure
$a = 5.5171(2)$1 2.12	Single cubic	$a = 5.5171(2)$2
$a = 5.5171(2)$3	Tetragonal superstructure (weak)	$a = 5.5171(2)$4
2.3–3.0	Two cubic phases (“plateau” region)	$a = 5.5171(2)$5 + $a = 5.5171(2)$6
3.45–3.76	Monoclinic $a = 5.5171(2)$7	$a = 5.5171(2)$8
$a = 5.5171(2)$9	Mixture: $b = 11.5061(5)$0, $b = 11.5061(5)$1, cubic	$b = 11.5061(5)$2 (at $b = 11.5061(5)$3)
$b = 11.5061(5)$4	Single cubic (high-P)	$b = 11.5061(5)$5

Two-phase regions always separate adjacent structures, in agreement with observed plateaus in pressure-composition isotherms (PCI). For $b = 11.5061(5)$6, a tetragonal (I$b = 11.5061(5)$7) superstructure reminiscent of YFe$b = 11.5061(5)$8D$b = 11.5061(5)$9 emerges. For $c = 9.4382(3)$0, a monoclinic $c = 9.4382(3)$1 variant (labeled $c = 9.4382(3)$2) dominates, and at higher $c = 9.4382(3)$3 ($c = 9.4382(3)$4), mixtures of monoclinic ($c = 9.4382(3)$5, $c = 9.4382(3)$6) and cubic phases prevail. At $c = 9.4382(3)$7 (high-pressure synthesis), fully cubic symmetry is restored.

The cell volume dependence on D content obeys an empirical relationship analogously established for Y–Fe$c = 9.4382(3)$8D$c = 9.4382(3)$9 systems: $_{4.2}$00 where $_{4.2}$01 Å$_{4.2}$02/f.u., $_{4.2}$03 Å$_{4.2}$04/D, $_{4.2}$05 Å$_{4.2}$06/D$_{4.2}$07.

5. Reassessment of “Rhombohedral” Phases

Earlier XRD studies such as Berthier et al. (1985) and Mushnikov et al. (1997) reported rhombohedral phases at $_{4.2}$08–$_{4.2}$09 based on observed distortion metrics. High-resolution synchrotron and neutron diffraction now reveal that these phases are in fact monoclinic $_{4.2}$10 distortions arising from symmetry lowering of the parent $_{4.2}$11 cubic cell, rather than true rhombohedral $_{4.2}$12. The ideal rhombohedral parameters, for a hypothetical cubic $_{4.2}$13 Å, are $_{4.2}$14 Å and $_{4.2}$15 Å. In contrast, in the monoclinic $_{4.2}$16 structure, $_{4.2}$17 Å, $_{4.2}$18 Å, $_{4.2}$19 Å, $_{4.2}$20, with $_{4.2}$21 significantly less than the $_{4.2}$22 of the undistorted rhombohedral limit.

Thus, regions formerly assigned as “rhombohedral” in the literature for $_{4.2}$23–$_{4.2}$24 and $_{4.2}$25 are more precisely described as monoclinic $_{4.2}$26 distortions. This clarification aligns TbFe$_{4.2}$27D$_{4.2}$28 phase assignments with the symmetries resolved by contemporary diffraction experiments.

6. Summary of Phase Evolution and Significance

The phase diagram of TbFe$_{4.2}$29D$_{4.2}$30 at 300 K is summarized as: cubic ($_{4.2}$31) for $_{4.2}$32, then tetragonal ($_{4.2}$33) at $_{4.2}$34, followed by two-phase cubic plateaus for $_{4.2}$35, monoclinic $_{4.2}$36 ($_{4.2}$37) for $_{4.2}$38, further monoclinic distortions ($_{4.2}$39 at $_{4.2}$40), and re-entrant cubic symmetry for $_{4.2}$41. All transitions occur via two-phase regions, indicative of first-order transitions and plateau formation in PCIs. The overall sequence and phase symmetry-lowering closely parallel observations in YFe$_{4.2}$42D$_{4.2}$43, providing a unified framework for rare-earth Laves-phase deuterides and clarifying ambiguous rhombohedral assignments in previous work (Paul-Boncour et al., 18 Dec 2025).