Next Article in Journal
Water Temperature and Salinity Measurement Using Frequency Comb
Next Article in Special Issue
The Conundrum of Relaxation Volumes in First-Principles Calculations of Charged Defects in UO2
Previous Article in Journal
An Effective Framework Using Spatial Correlation and Extreme Learning Machine for Moving Cast Shadow Detection
Previous Article in Special Issue
Assessing Relativistic Effects and Electron Correlation in the Actinide Metals Th to Pu
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a CALPHAD Thermodynamic Database for Pu-U-Fe-Ga Alloys

1
Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, CA 94550, USA
2
CEA-Centre de Valduc, 21120 Is sur Tille, France
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2019, 9(23), 5040; https://doi.org/10.3390/app9235040
Submission received: 11 October 2019 / Revised: 7 November 2019 / Accepted: 15 November 2019 / Published: 22 November 2019

Abstract

:
The interaction of actinides and actinide alloys such as the δ-stabilized Pu-Ga alloy with iron is of interest to understand the impurity effects on phase stability. A newly developed and self-consistent CALPHAD thermodynamic database is presented which covers the elements: Pu, U, Fe, Ga across their whole composition and temperature ranges. The phase diagram and thermodynamic properties of plutonium-iron (Pu-Fe) and uranium-iron (U-Fe) systems are successfully reassessed, with emphasis on the actinide rich side. Density functional theory (DFT) calculations are performed to validate the stability of the stoichiometric (Pu,U)6Fe and (Pu,U)Fe2 compounds by computing their formation enthalpies. These data are combined to construct the Pu-U-Fe ternary phase diagram. The thermodynamic assessment of Fe-Ga is presented for the first time and application to the quaternary Pu-U-Fe-Ga system is discussed.

1. Introduction

The interest in actinide elements is due to their complex physics and extends across many nuclear applications; specifically, plutonium and uranium, that are studied by various government entities, largely to focus on energy production. In that context, the phase stability with respect to chemical composition and temperature is paramount to understanding actinide materials behavior under normal, extended/aging and off-normal (accident) conditions. In this work, the focus lies in the effects of iron (a common impurity) on actinide (An) phase stability. In general, Fe has limited solubility in the actinides themselves, but can easily form the stoichiometric composition Pu6Fe and U6Fe, in addition to higher iron containing compounds (PuFe2 and UFe2). Since experiments can be rather tedious and expensive for actinide systems, computational models, such as the CALPHAD method is applied. Emphasis on Pu alloys with Ga acting as δ-stabilizer (fcc phase of plutonium) and U acting as a transmutation product are studied across potential composition and temperature ranges to elucidate phase relations with Fe acting as an impurity. Except for the Fe-Ga system, the binary systems included in this study of the Pu-U-Fe-Ga system have been extensively characterized in experiment (see literature review below) and preliminary thermodynamic models exist in the literature [1,2,3,4,5,6,7]: Pu-U [1], Pu-Fe [2,3], Pu-Ga [1], U-Fe [3,4], U-Ga [5,6], Fe-Ga [7].
As the focus is on small amounts of iron (to model impurity effects) this work aims to improve the currently accepted models of Pu-Fe and U-Fe of [3], especially on the actinide-rich side. Current shortcomings include little to no solubility of iron in the Pu allotropic phases as well as a neglect of the catatectic reaction described in the Pu-Fe section below. The liquidus on the Pu-rich side is also refined. Similarly, the solubility of iron in tetragonal uranium (β-U) is not well described in the currently accepted model and overall improvements in the accuracy of the invariant reactions is achieved. The new CALPHAD assessments are supported by density-functional-theory (DFT) electronic-structure calculations (i.e., enthalpy of formation of compounds). Finally, this paper includes the first attempt at modeling the Fe-Ga and An-Fe-Ga systems (with An = Pu, U).

2. Thermodynamic Modeling

2.1. CALPHAD Method

Computational Thermodynamics in the form of the CALculation of PHAse Diagrams (CALPHAD) method provides self-consistent Gibbs energy functions, for a given phase and/or structure [8], which are used to calculate equilibrium conditions of multicomponent systems. The model parameters are optimized using critically selected thermochemical and constitutive data as input. For the pure elements, the data from the SGTE [9] database is used and is expressed via the following temperature (T in K) dependent function, for element, x in a given phase, ϕ :
G x ϕ 0 = G ϕ H x SER = a + bT + cT ln ( T )   +   dT 2   + eT 3   + fT 1   + g n T n
In the above expression, H x SER is the molar enthalpy of component x in its standard element reference (SER) state at 298.15 K and 105 Pa. In the case of the elements described here, the standard states are α-U (orthorhombic), α-Fe (body centered cubic), α-Ga (orthorhombic) and α-Pu (simple monoclinic), where crystallographic structures are detailed in Table 1. The following phases: α-U, α-Pu, β-Pu, γ-Pu and δ’-Pu are treated as fixed compositions (neglecting solubility) and represented using (Equation (1)). The liquid, bcc (α-Fe, δ-Fe, ε-Pu, γ-U), fcc (γ-Fe, δ-Pu) and β-U are all treated as solution phases and are modelled as substitutional solutions yielding the following molar Gibbs energy expression for a given phase ϕ :
G m ϕ = i x i G i ϕ 0 + RT i x i ln ( x i )   +   G m ϕ phys   +   G m ϕ E
where the first term represents the mechanical mixing of the end-members ( G i ϕ 0 are the Gibbs energies of the pure elements in the structural state ϕ ), the second term represents the contribution due to the ideal entropy of mixing, the third includes physical models such as magnetic contribution and the fourth term ( G m ϕ ex ) represents the excess molar Gibbs energy.
For phases with order-disorder transformations (e.g., bcc_A2/B2), the partitioning model is used, which means the ordering is described as an addition to the Gibbs energy of the disordered phase so that the molar Gibbs energy is described as:
G m bcc = G m dis   +   Δ G m ord
where G m dis is the contribution to the Gibbs energy based on composition expressed by the disordered solution phases formulae (independent of the ordering state of the phase), and Δ G m ord is the additional Gibbs energy due to ordering, i.e., configuration (contribution due to the long-range ordering, and must be zero in the when the phase is disordered):
Δ G m ord = G m ord ( y i ) G m ord ( y i = x i )
where y i is the mole fraction of i per sublattice, also denoted site fraction, and x i the mole fraction. The parameter describing the ordering, G m ord , is first calculated with the original site fractions, y, which describe the ordering. The site fractions are then set equal to the mole fraction, x, which means that each constituent has the same site fraction in all sublattices (i.e., disordered states), and the value of the expression is calculated again. The difference is the contribution to the Gibbs energy due to ordering. If the phase was originally disordered, the two terms are equal and the difference is zero (no driving for ordering).
The excess Gibbs energy term from Equation (2) is modeled using the Redlich-Kister [10] formalism for binary solutions and applied to multi-component systems via the Muggianu extrapolation [11], within the assumption that possible higher-order many-body interaction parameters can be ignored (above binary interactions). The resulting excess Gibbs energy term for a higher order system is written as:
G m ϕ ex = i j > i c i c j v = 0 p L i , j ϕ v ( c i c j ) v
The Redlich-Kister model parameters, L i , j ϕ v , describe the deviation from ideality and are chosen such that they satisfactorily represent all the available thermochemical and phase diagram data of the corresponding binaries. Generally, L i , j ϕ v is expressed as:
L i , j ϕ v = a   +   bT
The excess Gibbs energy can also include, when necessary, ternary interaction parameters (e.g., L Fe , Pu , U ϕ ) for a multicomponent system. For n elements:
G m ϕ ex , tern = i = 1 n 2 j = i   +   1 n 1 k = j   +   1 n x i x j x k L i , j , k ϕ
where:
L i , j , k ϕ = v i L i , j , k ϕ i   + v j L i , j , k ϕ j   + v k L i , j , k ϕ k
with:
v i = x i   +   1 x i x j x k 3
v j = x j   +   1 x i x j x k 3
v k = x k   +   1 x i x j x k 3
In the case of a ternary system (e.g., Fe-Pu-U) where xi + xj + xk = 1, (Equation (8)) reduces to:
L Fe , Pu , U ϕ = x Fe L Fe , Pu , U ϕ Fe   + x Pu L Fe , Pu , U ϕ Pu   + x U L Fe , Pu , U ϕ U
with:
L Fe , Pu , U ϕ i = a   +   bT
Stoichiometric binary and ternary compounds are described by the following function, where x i is the mole fraction of element i in the phase/compound ϕ :
G comp ϕ = a   +   bT   +   i x i G i Φ 0
In this case, a and b are the adjustable parameters used to model the Gibbs energy of binary or ternary stoichiometric compounds.

2.2. Electronic Structure Calculations

The results from all electronic-structure calculations come from applying density-functional theory (DFT) that took shape in the seminal papers by Hohenberg, Kohn, and Sham [12,13]. This theory is in principle exact, but it relies on a fundamental approximation, namely, the so-called electron exchange and correlation functional. For actinides as well as iron [14] the generalized gradient approximation (GGA) is currently the best approximation and is applied here. Specifically, the assumption for all DFT calculations is GGA either in its original form (PW91) [15] or the similar but simplified; Perdew, Burke and Ernzerhof (PBE) functional form [16].
Before detailing the electronic-structure methods used here, it is necessary to mention that theories which go beyond DFT have been applied to 5f-electron systems including uranium and plutonium metals. Savrasov and Kotliar [17] applied intra-atomic Coulomb correlations with a large Hubbard U of the order of 4–4.5 eV in their DFT + U approach for plutonium and many other works using similar U (see references in [18]). Also, for uranium metal DFT + U has been applied [19] but it has been shown that with carefully executed DFT calculations this empirical parameter is not required for metallic uranium or plutonium (or alloys with iron) [18,20].
For the plutonium system it is required for good accuracy to extend the DFT to include orbital-orbital coupling. This can be done rather straightforwardly by adopting the orbital-polarization scheme [18]. The orbital-orbital interaction enhances the magnitude of the orbital magnetic moment (caused by the spin-orbit interaction) and results in more realistic magnetic and chemical bonding properties. Here all-electron full-potential linear muffin-tin orbital (FPLMTO) calculations include spin-orbit coupling and orbital polarization. The orbital polarization tends to increase the stability of the respective compound relative to the spin-orbit only case.
In terms of the technical details of the calculations, for the most part, the FPLMTO method is applied, which has been described in detail [21]. In this method, no approximations are made for the electron core states that lie deeper in energy than the valence states, unlike the so-called pseudopotential method. The present implementation does not make any assumptions beyond that of the electron exchange and correlation functional. Basis functions, electron densities, and potentials are calculated without any geometrical approximation and these are expanded in spherical harmonics inside non-overlapping (muffin-tin) spheres surrounding each atom and in Fourier series in the region between these muffin-tin spheres. The theory includes all relativistic corrections including spin-orbit coupling for d and f states but not for the p states [18].
The actinide metals have been shown to be very well described with 6s and 6p semi-core states and 7s, 7p, 5f, and 6d valence states [22] and this setup is replicated here.For the calculation of the formation enthalpy one calculates the energies for the An6Fe (An = Pu, U) compound and subtract the energies of the constituents in their ground-state phase (α-An and α-Fe). Ideally, if computationally feasible, these structures need to be optimized (adjusted or relaxed) so that the calculations produce the lowest total energy of the phase. For the FPLMTO calculations this has been done by calculating numerical forces from small atomic displacements [23]. For An6Fe, α-Pu, α-U, α-Fe, a total of 128, 54, 256, and 2000 k points are included, respectively in the electronic-structure calculations. Each energy eigenvalue is broadened with a Gaussian having a width of 20 mRy.
When calculating the formation enthalpy, it is important to know precisely what the crystal structure is. For the studied compounds they are known (Table 2), except for Pu6Fe where detailed experimental information is lacking. The prototype structure is the same as U6Fe, but the exact crystallographic parameters are not known. Therefore, it is necessary to calculate these parameters, or the energy will be too high for the compound. Hence, we have relaxed the crystal structure (i.e., found the structural parameters that correspond to the lowest total energy) for Pu6Fe. The obtained parameters are significantly different than those for U6Fe. For completeness, structural relaxation is performed for all FPLMTO calculations, but we find that for U6Fe the optimized structure is very close to the experimentally reported one.
The calculations referred to as exact muffin-tin orbital method (EMTO), are performed using the Green’s-function technique based on the improved screened Korringa-Kohn-Rostoker method, where the one-electron potential is represented by optimized overlapping muffin-tin (OOMT) potential spheres [24,25]. Inside the potential spheres the potential is spherically symmetric, and it is constant between the spheres. The radius of the potential spheres, the spherical potential inside these spheres, and the constant value in the interstitial region are determined by minimizing (i) the deviation between the exact and overlapping potentials, and (ii) the errors caused by the overlap between the spheres. Within the EMTO formalism, the one-electron states are calculated exactly for the OOMT potentials. As an output of the EMTO calculations, one can determine self-consistent Green’s function of the system and the complete, non-spherically symmetric charge density. Finally, the total energy is calculated using the full charge-density technique [26]. The valence states are treated as the 7s, 6p, 6d, and 5f states for U and Pu and 4s and 3d states for Fe. The corresponding Kohn-Sham orbitals are expanded in terms of spdf exact muffin-tin orbitals, i.e., an orbital momentum cutoff is applied, where lmax = 3. The EMTO orbitals, in turn, consist of the spdf partial waves (solutions of the radial Schrödinger equation for the spherical OOMT potential wells) and the spdf screened spherical waves (solutions of the Helmholtz equation for the OOMT muffin-tin zero potential). The completeness of the muffin-tin basis was discussed in detail in Ref. [25] and it was shown that for metals crystallizing in close-packed lattices lmax = 3 (spdf orbitals) leads to the well converged charge density and total energy. For the electron exchange and correlation energy functional, GGA (PBE96) is considered [27]. Integration over the Brillouin zone is performed using 25 × 25 × 25 k-points grid generated according to the Monkhorst-Pack scheme [28]. The moments of the density of states, needed for the kinetic energy and valence charge density, are calculated by integrating the Green’s function over a complex energy contour (with 1.9–2.4 Ry diameter) using a Gaussian integration technique with 40 points on a semi-circle enclosing the occupied states. When treating compositional disorder, the EMTO method is combined with the coherent potential approximation (CPA) [24]. For instance, the EMTO-CPA formalism has been successfully used to describe thermodynamic properties of the metallic nuclear fuels, including: U-Zr, U-Mo, Pu-Zr, Pu-Mo, Pu-U, Pu-Np, Pu-Am, U-Ti, U-Nb, Np-Mo, Pu-Mo, for fast breeder reactors [29,30,31,32,33,34,35,36].
Relaxing (or optimizing for lowest energy) the crystal structures are not currently feasible and for this reason all calculations applying the EMTO assumes the experimental structure. For Pu6Fe the experimental structure is not known in detail and in this case the parameters for U6Fe have also been applied for Pu6Fe.

3. Literature Review

As the thermodynamic assessments and their applications to the Pu-U [1,37], Pu-Ga [1,38,39,40,41], U-Ga [5,6] and Pu-U-Ga [1,42] systems are already available in the literature, the present study will only focus on the remaining systems (i.e., Pu-Fe, U-Fe, Fe-Ga, Pu-U-Fe, Pu-Fe-Ga, and U-Fe-Ga).

3.1. The Plutonium-Iron Phase Diagram

The Pu-Fe system has been experimentally investigated by Konobeevsky [43], Avivi [44], Mardon [45] and Ofte [46] using thermal analysis, dilatometry, metallography, x-ray diffraction (XRD), microhardness and viscosity measurements. The full phase diagram was studied by [43] and [45] whereas the authors in [44] focused on the PuFe2-Fe region and [46] investigated the liquidus range on the Pu-rich side up to 13 at. % Fe. Formation enthalpies of the intermetallic PuFe2 and Pu6Fe phases were measured/estimated by [47,48,49] and [49] respectively. Elemental plutonium exists in 6 allotropic phases, including: simple-monoclinic (α-Pu), body-centered monoclinic (β-Pu), face-centered orthorhombic (γ-Pu), face-centered cubic (fcc or δ-Pu), body-centered tetragonal (δ’-Pu) and the body-centered cubic (bcc or ε-Pu). The lower-temperature Pu-allotropes (α, β and γ) exhibit limited (negligible) solubility of iron [50], whereas the higher-temperature phases (except δ’-Pu) were found to have between 0.6–1.5 at. % Fe and 2–2.5 at. % Fe for δ-Pu and ε-Pu respectively, with Mardon [45] predicting higher solubilities than Elliott [51]. An interesting attribute of the Pu-rich side of the phase diagram includes a catatectic reaction (430 °C) between ε-Pu and δ-Pu; upon cooling of ε-Pu, a reappearance of the liquid phase (partial melting) occurs and is in equilibrium with the δ-Pu phase. The eutectic reaction: [Liquid ⟷ Pu6Fe + δ-Pu] was found at 90 at. % Pu by [45], whereas other authors located the eutectic around 89–89.5 and 91.5 at. % Pu by [44] and [46] respectively. Pu6Fe decomposes peritectically into PuFe2 and liquid phases at 430 °C according to [43] or 428 °C according to [45]. Thermal analysis measurements of Pu6Fe by [52] indicate an onset melting temperature at 411.5 °C with a melting enthalpy of 6.77 kJ/mol. The PuFe2 intermetallic compound consists of 3 distinct polymorphs, with transition temperatures of 771 °C and 1020 °C for α-PuFe2→β-PuFe2 and β-PuFe2→γ-PuFe2 respectively; the compound melts congruently at 1240 °C. A eutectic reaction [Liquid ⟷ PuFe2 + γ-Fe] was identified by [45] at 18 at. % Pu and 1165 °C, while other studies [43,44] placed the same reaction at 19 at. % Pu and 1180 °C [43]. Elemental iron transforms from bcc (α-Fe) to fcc (γ-Fe) at 912 °C and back to bcc (δ-Fe) at 1394 °C and melts at 1538 °C. Very low solubility of Pu is reported in α-Fe, at less than 0.02 at. % Pu [53] and a eutectoid reaction [γ-Fe ⟷ α-Fe + Fe2Pu] at 907 °C is observed by [44]. Higher solubility of Pu in the fcc (γ-Fe) phase was reported by [44] with an average of 0.7 at. % Pu. The peritectic transformation [L + δ-Fe ⟷ γ-Fe] is suggested by [54] due to the presumed lower solubility of Pu in the high temperature bcc (δ-Fe) phase. Based on these experimental investigations, a complete thermodynamic assessment of the Pu-Fe phase diagram has been reported by [2,3]. In these assessments, the Pu-rich side did not include the solubility of Fe in the δ-Pu and ε-Pu, which are pertinent to this study. Additionally, the intermetallic compounds were modelled using the liquid phase as their reference state, which is not consistent with the most recent and generally accepted CALPHAD-formulation for the standard refence state.

3.2. The Uranium-Iron Phase Diagram

The U-Fe system has been the source of many experimental investigations, with many authors performing thermal analysis including [55,56,57,58,59,60,61,62], some of whom also used dilatometry, metallography [63], XRD, microprobe analysis and calorimetric measurements [64]. The authors of [57,58,63] focused on the U-rich side of the phase diagram, where the solubility of iron in γ-U (bcc structure) was found to be 1.3 at. % Fe on average by [55,56,57,63]. The other two phases of uranium have no solubility reported for the orthorhombic (α-U) phase and limited information for the tetragonal (β-U) phase with solubility at 0.37 at. % Fe or less [58]. The allotropic phase transitions occur at 668 °C and 776 °C from α-U to β-U and β-U to γ-U, respectively. A eutectoid reaction is reported around 0.1 at. % Fe between 661 °C and 675 °C for [β-U ⟷ α-U + U6Fe]. A similar reaction was reported for [γ-U ⟷ β-U + U6Fe] between 0.5 and 0.8 at. %. Fe, with an average reported temperature of 762 °C [54,58,63]. The compound U6Fe forms from a peritectic decomposition of: [Liquid + γ-U ⟷ U6Fe] with authors reporting temperatures between 795 °C and 829 °C [55,56,59,60,64]. The low melting eutectic of [Liquid ⟷ U6Fe + UFe2] occurs at 725 °C at 66 at. % U according to [60]. The compound UFe2 was found to melt congruently between 1228 °C and 1235 °C [55,56,60,64] and includes a eutectic reaction [Liquid ⟷ UFe2 + γ-Fe] at 1080 °C and 17 at. % uranium. The solubility of U in α-Fe and γ -Fe are limited according to [55,56], since the δ-Fe phase is analogous to that of α-Fe, it is assumed that the solubility is negligible, although quantitative data is not available. The currently accepted model by Kurata [3] for the U-Fe phase diagram includes similar discrepancies as the Pu-Fe system; there is little to no Fe solubility considered within the U-phases and the intermetallic compounds are not consistent with the most recent and generally accepted CALPHAD-formulation for the standard refence state.

3.3. The Iron-Gallium Phase Diagram

The Fe-Ga system was largely investigated by [65,66,67,68,69,70,71,72] with various techniques including thermal analysis, dilatometry, XRD, light microscope and Mössbauer spectroscopy to cover the full composition range of the system and elucidate phase relations, especially the complex ordering of the iron-rich BCC (including D03) phases. The liquidus and solidus curves were measured by [65], who deduced the general shape of the phase diagram, but did not take into account the above-mentioned ordering on the iron rich side. Elemental gallium melts at 30 °C and does not exhibit any solubility of iron within its structure. Two Ga-rich phases were found to form peritectically; FeGa3 phase decomposes at 824 °C into the monoclinic Fe3Ga4 and liquid phases with a composition of 81.5 and 85 at. % Ga, according to [72] and [65] respectively. The second peritectic reaction occurs at 906 °C [L + α’(B2) ⟷ Fe3Ga4]. Both phases were found to have a small composition ranges by [65], which was confirmed by [67] for the Fe3Ga4 phase, but not the FeGa3 [67,72]. Additionally, Köster et al. [66,67] identified the Fe6Ga5 phase, which was not considered by [65], that forms via a peritectoid reaction [α’ + Fe3Ga4 ⟷ β-Fe6Ga5] at 800 °C and exists in its α-Fe6Ga5 polymorphic form below 778 °C. The Fe3Ga phase is found to be stable until 600 °C with a composition varying from 25–29 at. % Ga according to [65,71]. Köster et al. identified this phase as having two different polymorphs; the α-Fe3Ga which exhibits a composition range between 26.1 and 29.5 at. % Ga until 605 and 619 °C respectively, above these temperatures the β-Fe3Ga phase is stable, with a slightly narrower composition range and transforms at 681 °C and 82.5 at. % Fe into the α’’ (B2′) phase. Gallium solubility in BCC iron extends to a maximum composition of 47.5 at. % gallium and is accompanied by various ordered structures, including the α’ (B2) CsCl-type phase with composition range from 31.5 at. % Ga up to the aforementioned maximum. The existence of this phase was verified by [68] and [71] via quenching from the liquid. Köster et al. [66] found that the α’ (B2) phase crystallizes directly from the melt at 1037 °C and 35.5 at. % Ga, and the phase is stable down to 950 °C, below which the α’’ (B2′) exists until 650 °C where it undergoes a eutectoid decomposition to α’’ (D03) and α-Fe (A2) and is coupled to the magnetic transition temperature of α-Fe. Köster et al. [66,67] also identified an ordered phase α’’’ that is involved in a eutectoid reaction: [α’’’ ⟷ α-Fe + α-Fe3Ga]. The authors however were skeptical of the α’’’ phase being an equilibrium structure, since most of their samples in this phase region also contained traces of β-Fe3Ga and even found some samples containing two-phases: α-Fe + β-Fe3Ga, making the latter structure the probable equilibrium phase. Köster et al. [66,67] indicated that the phase field between 21–26 at. % Ga and 580–680 °C is a likely metastable domain as reactions proceed extremely slowly and with difficulty in reproducing the data. Finally, there is limited solubility of gallium in γ-Fe (FCC), with a characteristic γ-loop in which the maximum solubility is 2.8 at. % Ga at 1140 °C.

3.4. The Plutonium-Uranium-Iron Phase Diagram

The Pu-U-Fe system was experimentally investigated by [2,3,73,74,75,76]. The first two papers consider a metallic fuel, either U-Pu-Zr [73] or U-Pu [74] with iron as a diffusion couple at 650 °C and 670 °C respectively. The experimental techniques include electroprobe microanalysis, EDX and SEM. It was found that the higher Pu containing fuels lead to increased liquid formation, since the liquid phase field extends into the (Pu,U)6Fe domain from the low-melting Pu6Fe compound. The solubility limits were found to be 16 at. % Pu at 650 °C and estimated (based on models) to be 12 at. % Pu at 670 °C in the U6Fe compound. Nakamura [75] and Kurata [2] prepared ternary Pu-U-Fe samples of varying composition by arc-melting, with [75] using DTA to elucidate ternary phase transitions from the prepared samples and [2] annealing the samples at 650 °C and characterizing them using XRD, SEM-EDS and chemical analysis. They were therefore able to establish a thermodynamic database derived from the information of these studies. This was later amended by Kurata et al. [3] to incorporate minor actinides within a CALPHAD database.

3.5. The Plutonium-Iron-Gallium Phase Diagram

No experimental data exists for this ternary system.

3.6. The Uranium-Iron-Gallium Phase Diagram

There have been six experimentally reported intermetallic compounds in the U-Fe-Ga phase space. They are UFeGa [77], UFeGa5, [78,79,80,81] UFe6Ga6 [82], UFe5Ga7 [83], U2FeGa8 [84] and U4FeGa12 [85]. The crystal structures of these phases have been investigated and low temperature experiments (up to 25 °C), including heat capacity, magnetic susceptibility and electrical resistivity have been performed [77,78,79,80,81,82,83,84,85]. There is no known information on phase stability above room temperature and a thermodynamic assessment with the given information is not possible.

4. Crystal Structure

The crystallographic information for the unary, binary and ternary systems described above are listed in Table 1, Table 2 and Table 3 respectively, the elemental crystal structures are also derived in the text. Information on structures related to the remaining binary systems are detailed in [1] for Pu-U, U-Ga and Pu-Ga. There are no specific compounds known for the ternary systems investigated.

5. Results and Discussion

The re-assessed parameters for the Pu-Fe, U-Fe and Pu-U-Fe system are listed in Table 4 along with the newly assessed Fe-Ga parameters. Since no thermodynamic information is available on the Pu-Fe-Ga and U-Fe-Ga ternary systems, the extrapolation from binary parameters were deemed sufficient. Table 5 compares the invariant reactions for each of the binaries with experiment and available models in the literature, Table 6 compares the formation enthalpies found in experimental and theoretical literature to the values calculated from the electronic structure methods described in this work and those obtained with the newly assessed CALPHAD model.

5.1. Pu-Fe Assessment

The assessment of the Pu-Fe phase diagram was performed across the full composition and temperature ranges. The calculated phase diagram for the Pu-Fe system is indicated in Figure 1 and compared to available experimental data from Mardon [45] and Ofte [46].
In Table 5, the current assessment is also compared to that of Kurata [3] and considerable improvements are to be noted. For the iron-rich side, the congruent melting temperature of the PuFe2 compound is improved and the higher solubility limit of Pu in γ-Fe is realized. Note that most of the amendments were made of the Pu-rich side (Figure 2), as that is the domain where iron-impurities are of interest. For the low-temperature allotropes (α, β and γ) of plutonium, iron solubility was not considered as there is limited (or negligible) experimental evidence. An overall improved fit of the liquidus and solidus data on the Pu-rich side was achieved. Specifically, the solubility of iron in both δ-Pu and ε-Pu were improved (from 0.4 at. % to 0.9 at. %, and 1.9 at. % to 2.7 at. %, respectively) and in better agreement with literature (1.3 at. % Fe and 2.4 at. % Fe, respectively). The invariant reaction temperatures were also ameliorated across the phase diagram, compared to the previous assessment, especially in the Pu-Pu6Fe phase domain. The stoichiometric compounds were re-assessed to include the generally used standard reference states (α-Pu and α-Fe) compared to Kurata et al. [3], who used the liquid phase as their reference state. Table 6 indicates that the formation enthalpy of PuFe2 calculated from the thermodynamic database in this work is in good agreement with the solution calorimetry measurements of [49] but is much less stable than the experimental value reported in [48] using e.m.f. techniques and the EMTO calculations presented in this work. Initial attempts to use the EMTO calculated values to optimize the Pu-Fe system included a more stable formation enthalpy for PuFe2 were unsuccessful.
Since the relaxed structure calculations using FPLMTO resulted in a less stable formation enthalpy, it is regarded as a more suitable starting point for fitting the PuFe2 intermetallic. In order to fit both liquidus data of [45] and include an accurate description of the PuFe2 compounds congruent melting, the formation enthalpy was lowered and found to be less stable than that of the EMTO calculations, the estimated values (based on the analogous UFe2) of [47] and the experimental work of [48]. The value for Pu6Fe calculated from CALPHAD lies in between that calculated via DFT and the experimentally available values, where the EMTO value (with an uncertainty of 1.3 kJ/mol) was used as a starting point for the parameters during optimization. The enthalpy of melting (6.77 kJ/mol) derived by [52] at 411.5 °C for the Pu6Fe compound agrees well with the peritectic decomposition temperature of 423.5 °C and reaction enthalpy of 6.51 kJ/mol calculated from this work.

5.2. U-Fe Assessment

The calculated phase diagram resulting from the present CALPHAD assessment is depicted in Figure 3. Overall improvements including the solubility of iron into the uranium lattice (β-U and γ-U) is provided and is indicated in the invariant reactions described in Table 5.
The improved assessment includes better agreement for the congruent melting of UFe2 compared to Kurata [3] and Chatain [4], as well as some adjustments to the eutectic temperatures. Similarly, the eutectoid reaction [γ-U ⟷ β-U + U6Fe] was improved with respect to U-composition. The peritectic reaction (reaction 20, Table 5) temperature was updated to include the most recent measurements by [64], which was also considered in the works by Chatain et al. [4], but neglected in the most recent assessment by Kurata [3]. This CALPHAD assessment is in good agreement for the formation enthalpies for both compounds (UFe2 and U6Fe) when compared to literature and the values calculated by DFT from this work, which were used initial points for the fit.

5.3. Fe-Ga Assessment

This is the first published CALPHAD assessment of the Fe-Ga system (Figure 4), based on a preliminary internal report by coauthors [7]. Since ordering of the iron-rich side is not pertinent to the scope of this investigation and given the uncertainties reported in experimental studies [65,66,67], it is therefore not considered at this time. However, the overall shape of the phase diagram is reproduced in the present work, relevant invariant reactions and stoichiometric compounds are included, and improvements to the preliminary work of Turchi et al. [7] have been made so that some of the invariant reactions are better represented. This includes adjustments to the peritectic reaction of Fe3Ga4, where the reaction temperature is improved by 30 °C compared to [7] and the composition is in better agreement with experiments by [65,66,67]. Similarly, reaction temperatures and composition ranges are improved for reactions 3 and 5 in Table 5, however the best fit of experimental data resulted in a eutectic instead of a peritectic reaction for reaction 3. Major modifications to the liquidus and solidus ranges result in a more acceptable fit compared to the preliminary assessment. Since the α-Fe (BCC_A2) phase is considered without any possible ordering, some disagreement on the Fe-rich side is to be expected. While the literature review (presented above) includes adequate descriptions of the ordering of α’-Fe, α’’-Fe and α’’’-Fe (Fe3Ga composition dependence), it is to be noted that experiments proved difficulty in achieving equilibrium conditions in this phase space, which supports the current decision to exclude ordering on Fe-rich side. Similarly, the focus is on small amounts of Fe-impurities, therefore, only stoichiometric compositions were presented for each of the compounds and parameters for any ordered iron phase were neglected. This provided some difficulty when optimizing this region of the phase diagram and further investigation would be useful. The liquid phase is modeled using three interaction parameters, two of them temperature dependent. Even with multiple parameters, the liquidus curve on the Ga-rich side is difficult to reproduce due to the sharp decrease close to 100 at.% Ga. As additional terms or parameters for the liquid phase do not lead to significant improvement, the authors believe that including ordering for the BCC phase is necessary to have a better representation on the liquidus and solidus line across the phase diagram. Otherwise, this work shows rather acceptable improvements over the preliminary assessment of [7] and is a satisfactory representation of the Fe-Ga phase equilibria as seen experimentally.

5.4. Pu-U-Fe Assessment

The binary Pu-Fe and U-Fe reassessed parameters were used to extrapolate to the ternary Pu-U-Fe system, where some ternary interaction parameters were necessary to improve agreement with experiments. The most recent work of Kurata et al. [3] included an assessment of both the Pu-Fe and U-Fe binary, but did not cover the Pu-U-Fe ternary system. A prior assessment by the same authors [2] includes ternary Pu-U-Fe parameters, however their binary descriptions are not consistent with their most up to date publications [3]. Therefore, a new consistent assessment based on the newly assessed binary parameters presented in this study was necessary to properly fit the experimental data presented by [2]. The calculated phase diagram for the available isotherm at 650 °C is depicted in Figure 5 and compared with experimental data points reported by Kurata et al. [2]. The current assessment shows improvement over the binary extrapolation and includes more detailed description of the An-rich corners, which was previously lacking. Some enhancements in the tie-lines compared previous work is achieved on the U-rich corner, however the fit of the Liquid ⟷ (Pu,U)Fe2 equilibrium could be improved by providing more experimental data on the Pu-rich side. The adjustment of further ternary parameters proved excessive without much improvement compared to the model of [2]. The same phase diagram is also calculated at room temperature, with the Pu-rich corner shown for clarity (Figure 6). Figure 6b shows that Pu-rich alloys include the (Pu,U)6Fe intermetallic as well as low temperature ζ-phase (from the Pu-U system), with both of these phases consuming the α-Pu matrix as a function of increased U content, whereas Fe leads to the formation of a second intermetallic (Pu,U)Fe2.

5.5. Extrapolation to An-Fe-Ga Systems (An = Pu, U)

The extrapolation to the Pu-Fe-Ga and U-Fe-Ga systems are indicated at 1000 °C and room temperature in Figure 7, Figure 8 and Figure 9. In the Pu-Fe-Ga system, there are no additional known ternary compounds so that the predicted phase diagrams are constructed via extrapolation of the pertinent binary systems. While there are some known intermetallics in the U-Fe-Ga system, there is no thermodynamic information about their stability above room temperature as all these compounds have been studied at T < 25 °C for magnetic purposes as indicated in the literature review. They are therefore not considered in this work, however, future experiments on these compounds would be useful for thermodynamic assessments. The Pu-Fe-Ga phase diagram is plotted at 1000 °C in Figure 7; phases from the binary systems which are stable at high temperature phases such as α-Fe (BCC), γ-Fe (FCC) and the PuFe2 intermetallic, and Pu-Ga compounds are among the first to form from the liquid phase. The U-analogous diagram in Figure 8 shows similar phase formation with the addition of γ-U (BCC) being stable for the uranium rich portion. More importantly, the room temperature phase diagrams show the formation of both Pu6Fe and U6Fe intermetallics for the An-rich region of the Pu-Fe-Ga and U-Fe-Ga systems (Figure 9a,b), respectively, which are accompanied by An-Ga phases within the actinide matrix. Closer interpretation is given to the phase stability and formation with respect to temperature for two different alloy compositions in the following section.

5.6. Prediction of Quaternary Pu-U-Fe-Ga System

Two property diagrams are given to indicate the phase behavior of the quaternary system as function of temperature: Pu89-U5-Ga5-Fe1 and Pu92.9-U5-Ga2-Fe0.1 (at. %). These calculations assume that Ga stabilizes the δ-Pu phase down to room temperature (metastable state stabilized by the extremely slow eutectoid decomposition of δ-Pu into α-Pu and Pu3Ga [94], therefore low temperature allotropes of Pu are not considered for the Ga concentration considered here (2 and 5 at. %). All of the Ga is therefore contained within the FCC phase of δ-Pu. In Figure 10 (Pu89-U5-Ga5-Fe1), the predominant phases at room temperature are δ-Pu, the ζ-phase, which exists in the binary Pu-U phase space, and the (Pu,U)6Fe intermetallic. The ζ-phase exists at low temperatures and exhibits increasing uranium solubility with respect to temperature, extending up to almost 70 at. % uranium at temperatures up to 588 °C [1] in the Pu-U system. It is this phase that accommodates most of the U from the alloy composition in Figure 10; a small remaining amount of uranium is contained in the (Pu,U)6Fe intermetallic, where the uranium solubility in the intermetallic slightly increases over a small temperature. As the latter phase disappears, the Fe is accommodated within another lower An-containing intermetallic (Pu,U)Fe2, with this phase accounting for slightly more U than Pu. Once (Pu,U)Fe2 is saturated at its stoichiometric ratio, the remaining Pu released from (Pu,U)6Fe is reincorporated into the Pu-matrix as clearly indicated by the increase of δ-Pu phase fraction in Figure 10. Finally, with the disappearance of the low temperature ζ-phase, the Pu-U binary η-phase forms at higher temperatures, accommodating a larger amount of Pu leading to a sharp decrease in δ-Pu. At 462 °C the first appearance of liquid occurs, indicative of the upper stability limit of (Pu,U)Fe2, which is accompanied by a decrease in the η-phase, resulting in an iron and uranium rich liquid. At 491 °C, the Pu-U binary η-phase transitions into a BCC solid solution of ε-Pu and γ-U, while δ-Pu phase continues to transform into BCC (ε-Pu) and eventually melts at 666 °C. It is interesting to note that even minor additions (1 at. %) of iron induces a drastic decrease of the incipient melting point of the alloy by 195 °C, compared to Pu90-U5-Ga5. The evolution of the elemental distribution per phase is plotted as a function of temperature in Figure 11 for the Pu89-U5-Ga5-Fe1 alloy composition.
Figure 12 simulates an iron impurity of about 1000 ppm with the atomic composition of Pu92.9U5-Ga2-Fe0.1. As with the previous property diagram, the three stable phases at room temperature are δ-Pu, the ζ-phase and (Pu,U)6Fe. With small (trace) amounts of Fe, only the An-rich intermetallic is formed. As the stability limit of (Pu,U)6Fe is reached as a function of temperature, the Fe becomes soluble in δ-Pu which increases with temperature. The uranium in the alloy is either accommodated within the ζ-phase at low temperature and η-phase at higher temperature. The formation of η-phase is able to include higher Pu content and therefore depletes δ-Pu until the transition to ε-Pu around 483 °C where the BCC lattice also accommodates γ-U and Fe until its melting point.

6. Conclusions

From the extensive literature review presented here, the Pu-Fe and U-Fe systems have been reassessed to include more details on the solubility of Fe in the actinide rich composition ranges. Density functional theory calculations of the formation enthalpy were performed to use for the CALPHAD optimizations for both intermetallic phases (An6Fe and AnFe2). Pu-Fe and U-Fe assessments were improved to include a description of impurity level iron in the respective actinides. These binary assessments were used to reassess the Pu-U-Fe system, which is in agreement with available experimental data. A first CALPHAD assessment of the Fe-Ga phase diagram is presented, where the iron rich ordered phases are neglected. Finally, property diagrams were used to simulate alloys with varying Fe content to show the impact of impurities on Pu-U-Ga alloys. Experimental investigations of the Pu-Fe-Ga and U-Fe-Ga systems would be useful for future assessments.

Author Contributions

Conceptualization, E.E.M., P.E.A.T., B.O., S.A.S. and A.P.; Funding acquisition, J.L.B., S.A.S. and A.P.; Investigation, E.E.M., A.L. and P.S.; Project administration, P.E.A.T., J.L.B., S.A.S. and A.P.; Supervision, A.P.; Validation, E.E.M.; Writing—original draft, E.E.M., A.L. and P.S.; Writing—review & editing, E.E.M., P.E.A.T., B.O., J.L.B., S.A.S. and A.P.

Funding

This research received no external funding.

Acknowledgments

This work was performed under the auspices of the U.S Department of Energy by Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344. This work was done as part of the international agreement on cooperation between DOE-NNSA and CEA-DAM in fundamental science supporting stockpile stewardship.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Perron, A.; Turchi, P.E.A.; Landa, A.; Söderlind, P.; Ravat, B.; Oudot, B.; Delaunay, F.; Kurata, M. Thermodynamic re-assessment of the Pu-U system and its application to the ternary Pu-U-Ga system. J. Nucl. Mater. 2014, 454, 81–95. [Google Scholar] [CrossRef]
  2. Kurata, M.; Nakamura, K.; Ogata, T. Thermodynamic evaluation of the quaternary U-Pu-Zr-Fe system-assessment of cladding temperature limits of metallic fuel in a fast reactor. J. Nucl. Mater. 2001, 294, 123–129. [Google Scholar] [CrossRef]
  3. Kurata, M. Thermodynamic database on U-Pu-Zr-Np-Am-Fe alloy system II- Evaluation of Np, Am and Fe containing systems. In Proceedings of the IOP Conference Series: Material Science and Engineering, San Francisco, CA, USA, 12–17 July 2009. [Google Scholar] [CrossRef]
  4. Chatain, S.; Guéneau, C.; Labroche, D. Thermodynamic assessment of the Fe-U binary system. J. Phase Equilib. 2003, 24, 122–131. [Google Scholar] [CrossRef]
  5. Wang, J.; Liu, Z.J.; Wang, C.P. Thermodynamic calculation of phase equilibria of the U-Ga and U-W systems. J. Nucl. Mater. 2008, 380, 105–110. [Google Scholar] [CrossRef]
  6. Moussa, G.; Berche, A.; Barbosa, J.; Pasturel, M.; Stepnik, B.; Tougait, O. Experimental investigation of the phase equilibria and thermodynamic assessment in the U-Ga and U-Al-Ga systems. J. Nucl. Mater. 2018, 499, 361–371. [Google Scholar] [CrossRef]
  7. Turchi, P.E.A.; Kaufman, L.; Liu, Z.-K.; Zhou, S. Thermodynamics and Kinetics of Phase Transformations in Plutonium Alloys-Part 1; UCRL-TR-206658; Lawrence Livermore National Lab.: Livermore, CA, USA, 2004. [Google Scholar]
  8. Lukas, H.L.; Fries, S.G.; Sundman, B. Computational Thermodynamics, The Calphad Method; Cambridge University Press: New York, NY, USA, 2007. [Google Scholar]
  9. Dinsdale, A.T. SGTE data for pure elements. Calphad 1991, 15, 317–425. [Google Scholar] [CrossRef]
  10. Redlich, O.; Kister, A.T. Algebraic representation of thermodynamic properties and the classification of solutions. Ind. Eng. Chem. 1948, 40, 345–348. [Google Scholar] [CrossRef]
  11. Muggianu, Y.; Gambino, M.; Bros, J. Enthalpies of formation of liquid alloys bismuth-gallium-tin at 723 K—Choice of an analytical representation of integral and partial thermodynamic functions of mixing for this ternary system. J. Chim. Phys. Phys-Chim. Biol. 1975, 72, 83–88. [Google Scholar] [CrossRef]
  12. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864–B871. [Google Scholar] [CrossRef]
  13. Kohn, W.; Sham, L. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1338. [Google Scholar] [CrossRef]
  14. Söderlind, P.; Gonis, A. Assessing a solids-biased density-gradient functional for actinide metals. Phys. Rev. B 2010, 82, 033102. [Google Scholar] [CrossRef]
  15. Perdew, J.P. Electronic Structures of Solids; Ziesche, P., Eschrig, H., Eds.; Springer: Berlin, Germany, 1991; pp. 11–20. [Google Scholar]
  16. Perdew, J.P.; Ruzsinszky, A.; Gabor, G.O.; Vydrov, O.A.; Scuseria, G.E.; Constantin, L.A.; Zhou, X.; Kieron, B. Restoring the Density-Gradient Expansion for Exchange in Solids and Surfaces. Phys. Rev. Lett. 2008, 100, 136406. [Google Scholar] [CrossRef] [PubMed]
  17. Savrasov, S.Y.; Kotliar, G. Ground-State Theory of δ-Pu. Phys. Rev. Lett. 2000, 84, 3670–3673. [Google Scholar] [CrossRef] [PubMed]
  18. Söderlind, P.; Landa, A.; Sadigh, B. Density-functional theory for plutonium. Adv. Phys. 2019, 68, 1–47. [Google Scholar] [CrossRef]
  19. Xie, W.; Xiong, W.; Marianetti, C.A.; Morgan, D. Correlation and relativistic effects in U metal and U-Zr alloy: Validation of ab initio approaches. Phys. Rev. B 2013, 88, 235128. [Google Scholar] [CrossRef]
  20. Söderlind, P.; Landa, A.; Turchi, P.E.A. Comment on “Correlation and relativistic effects in U metal and U-Zr alloy: Validation of ab initio approaches”. Phys. Rev. B 2014, 90, 157101. [Google Scholar] [CrossRef]
  21. Wills, J.M.; Alouani, M.; Andersson, P.; Delin, A.; Eriksson, O.; Grechnev, O. Full Potential Electronic Structure Method; Springer: Berlin, Germany, 2010. [Google Scholar]
  22. Söderlind, P. Theory of the crystal structures of cerium and the light actinides. Adv. Phys. 1998, 47, 959–998. [Google Scholar] [CrossRef]
  23. Söderlind, P.; Grabowski, B.; Yang, L.; Landa, A.; Björkman, T.; Souvatzis, P.; Eriksson, O. High-temperature phonon stabilization of γ-uranium from relativistic first-principles theory. Phys. Rev. B 2012, 85, 060301. [Google Scholar] [CrossRef]
  24. Vitos, L.; Abrikosov, I.A.; Johansson, B. Anisotropic Lattice Distortions in Random Alloys from First-Principals Theory. Phys. Rev. Lett. 2001, 87, 156401. [Google Scholar] [CrossRef]
  25. Vitos, L. Computational Quantum Mechanics for Materials Engineers: The EMTO Method and Application; Springer: London, UK, 2007. [Google Scholar]
  26. Kollar, J.; Vitos, L.; Skriver, H.L. Electronic Structure and Physical Properties in Solids: The Uses of the LMTO Method; Springer: Berlin, Germany, 2000; pp. 85–113. [Google Scholar]
  27. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  28. Monkhorst, H.P.J. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  29. Landa, A.; Söderlind, P.; Turchi, P.E.A. Density-functional study of the U-Zr system. J. Alloy Compd. 2009, 478, 103–110. [Google Scholar] [CrossRef]
  30. Landa, A.; Söderlind, P.; Turchi, P.E.A.; Vitos, L.; Ruban, A. Density-functional study of Zr-based actinide alloys: 2. U-Pu-Zr system. J. Nucl. Mater. 2009, 393, 141–145. [Google Scholar] [CrossRef]
  31. Landa, A.; Söderlind, P.; Turchi, P.E.A.; Vitos, L.; Ruban, A. Density functional study of Zr-based actinide alloys. J. Nucl. Mater. 2009, 385, 68–71. [Google Scholar] [CrossRef]
  32. Landa, A.; Söderlind, P.; Turchi, P.E.A.; Vitos, L.; Ruban, A. Density-functional study of bcc Pu-U, Pu-Np, Pu-Am, Pu-Cm alloys. J. Nucl. Mater. 2009, 408, 61–66. [Google Scholar] [CrossRef]
  33. Bajaj, S.; Garay, A.; Landa, A.; Söderlind, P.; Turchi, P.E.A.; Arróyave, R. Thermodynamic study of the Np-Zr system. J. Nucl. Mater. 2011, 409, 1–8. [Google Scholar] [CrossRef]
  34. Landa, A.; Söderlind, P.; Turchi, P.E.A. Density-functional study of bcc U-Mo, Np-Mo, Pu-Mo, and Am-Mo alloys. J. Nucl. Mater. 2013, 434, 31–37. [Google Scholar] [CrossRef]
  35. Turchi, P.E.A.; Landa, A.; Söderlind, P.A. Thermodynamic assessment of the Am-Pu system with input from ab initio. J. Nucl. Mater. 2011, 418, 165–173. [Google Scholar] [CrossRef]
  36. Bajaj, S.; Söderlind, P.; Turchi, P.E.A.; Arróyave, R. The U-Ti system: Strengths and weaknesses of the CALPHAD method. J. Nucl. Mat. 2011, 419, 177–185. [Google Scholar] [CrossRef]
  37. Perron, A.; Turchi, P.E.A.; Landa, A.; Söderlind, P.; Ravat, B.; Oudot, B.; Delaunay, F. The Pu-U-Am system: An ab initio informed CALPHAD thermodynamic study. J. Nucl. Mater. 2015, 458, 425–441. [Google Scholar] [CrossRef] [Green Version]
  38. Turchi, P.E.A.; Kaufman, L.; Zhou, S.; Liu, Z.-K. Thermostatics and kinetics of transformations in Pu-based alloys. J. Alloy Compd. 2007, 444, 28–35. [Google Scholar] [CrossRef] [Green Version]
  39. Perron, A.; Ravat, B.; Oudot, B.; Lalire, F.; Mouturat, K.; Delaunay, F. Phase transformations in Pu-Ga alloy: Synergy between simulations and experiments to elucidate direct and indirect reversion competition. Acta Mater. 2013, 61, 7109–7120. [Google Scholar] [CrossRef]
  40. Ravat, B.; Oudot, B.; Perron, A.; Lalire, F.; Delaunay, F. Phase transformations in PuGa 1 at. % alloy: Study of whole reversion process following martensitic transformation. J. Alloy Compd. 2013, 580, 298–309. [Google Scholar] [CrossRef]
  41. Ravat, B.; Lalire, F.; Oudot, B.; Appolaire, B.; Aeby-Gautier, E.; Panisiot, J.; Delaunay, F. Phase transformations in PuGa 1 at. % alloy: Influence of stress on δ→α’ martensitic transformation at low temperatures. Materialia 2019, 6, 100304. [Google Scholar] [CrossRef]
  42. Perron, A.; Turchi, P.E.A.; Landa, A.; Oudot, B.; Ravat, B.; Delaunay, F. Thermodynamic assessments and inter-relationships between systems involving Al, Am, Ga, Pu and U. J. Nucl. Mater. 2016, 482, 187–200. [Google Scholar] [CrossRef] [Green Version]
  43. Konobeevsky, S.T. Phase Diagram of Some Plutonium Systems, Conference on the Peaceful Use of Atomic Energy, Moscow, Soviet Union. 1955.
  44. Avivi, E. Studies of Plutonium-Iron and Uranium-Plutonium-Iron Alloys. Doctorate Thesis/Dissertation, CEA-Fontenay-aux-Roses, Université de Paris, Paris, France, 1964. [Google Scholar]
  45. Mardon, P.G.; Haines, H.R.; Pearce, J.H.; Waldron, M.B. The Plutonium-Iron system. J. Inst. Met. 1957, 86, 166–171. [Google Scholar]
  46. Ofte, D.; Wittenberg, L.J. Viscosity-Composition Relationships in Molten Plutonium-Iron Alloys. Trans. ASM 1964, 57, 916–924. [Google Scholar]
  47. Chiotti, P.; Akhachinskij, V.V.; Ansara, I.; Rand, M.H. The Chemical Thermodynamics of Actinide Elements and Compounds; International Atomic Energy Agency: Vienna, Austria, 1981; pp. 114–117, 228–230. [Google Scholar]
  48. Campbell, G. Plutonium and Other Actinides; Blanck, H., Lindner, R., Eds.; North Holland: Amsterdam, The Netherlands, 1976. [Google Scholar]
  49. Akhachinskij, V.V.; Koputin, L.M.; Ivanov, I.; Podol’skaya, N.S. Thermodynamics of Nuclear Materials. In Proceedings of the Symposium on Thermodynamics of Nuclear Materials, Vienna, Austria, 21–25 May 1962. [Google Scholar]
  50. Schonfeld, F.W. Plutonium Phase Diagrams Studied at Los Alamos; Coffinberry, A.S., Miner, W.N., Eds.; University of Chicago Press: Chicago, IL, USA, 1961. [Google Scholar]
  51. Elliott, R.O.; Larson, A.C. The Metal Plutonium; Coffinberry, A.S., Miner, W.N., Eds.; University of Chicago Press: Chicago, IL, USA, 1961; pp. 265–280. [Google Scholar]
  52. Schwartz, D.S.; Tobash, P.H.; Richmond, S. Thermal Analysis of Pu6Fe Synthesized from Hydride Precursor; Material Research Society: San Francisco, CA, USA, 2014. [Google Scholar]
  53. Moreau, G.; Calais, D. Solubility of Plutonium in Iron. J. Nucl. Mater. 1967, 24, 121–124. [Google Scholar] [CrossRef]
  54. Okamoto, H. Phase Diagrams of Binary Actinide Alloys; Kassner, M.E., Peterson, D.E., Eds.; ASM International: Materials Park, OH, USA, 1995; pp. 314–317. [Google Scholar]
  55. Gordon, P.; Kaufmann, A.R. Uranium-Aluminum and Uranium-Iron. Trans. Metall. AIME 1950, 188, 182–194. [Google Scholar]
  56. Grogan, J.D. The Uranium-Iron System. J. Inst. Met. 1950, 77, 571–576. [Google Scholar]
  57. Bellot, J.; Blanchon, A.; Chazot, R.; Dosiere, P.; Henry, J.-M.; Colas, M. Uranium-Iron Equilibrium Diagram for Dilute Iron Concentration Range. C. R. 1958, 246, 3063–3065. [Google Scholar]
  58. Straatman, J.A.; Neumann, N.F. Equilibrium Structures in the High Uranium-Iron Alloy System; Topical Report. No. MCW-1487 MCW-1487; Mallinckrodt Chemical Works: Welding Spring, Missouri, USA, 1964. [Google Scholar]
  59. Michaud, G.G. A study of the iron-rich portion of the Fe-U phase diagram. Can. Metall. Q. 1966, 5, 355–365. [Google Scholar] [CrossRef]
  60. Chapman, L.R.; Holcombe, C.F. Revision of the Uranium-Iron Phase Diagram. J. Nucl. Mater. 1984, 126, 323–326. [Google Scholar] [CrossRef]
  61. Leibowitz, L.; Blomquist, R.A. Thermodynamics and Phase Equilibria of the Iron-Uranium System. J. Nucl. Mater. 1991, 184, 47–52. [Google Scholar] [CrossRef]
  62. Gardie, P.; Bordier, G.; Poupeau, J.J.; Le Ny, J. Thermodynamic Activity Measurements of U-Fe and U-Ga Alloys by Mass Spectrometry. J. Nucl. Mater. 1992, 189, 85–96. [Google Scholar] [CrossRef]
  63. Swindells, N. The Solubility of Iron in Solid Uranium Between 0.003 wt. % and 0.3 wt. % Iron. J. Nucl. Mater. 1966, 18, 261–271. [Google Scholar] [CrossRef]
  64. Labroche, D. Contribution a l’etude Thermodynamique du Systeme Ternaire U-Fe-O. Ph.D. Thesis/Dissertation, Institut National Polytechnique de Grenoble, Grenoble, France, 2000. [Google Scholar]
  65. Dasarathy, C.; Hume Rothery, W. The system iron-gallium. Proc. R. Soc. Lond. Ser. A 1965, 286, 141–157. [Google Scholar]
  66. Köster, W.; Gödecke, T. Über den Aufbau des Systems Eisen-Gallium zwischen 10 und 50 At.-% Ga und dessen Abhägigkeit von der Wärmebehandlung, I. Das Diagramm der raumzentrierten Phasen. Z. Metallk. 1977, 68, 582–589. [Google Scholar]
  67. Köster, W.; Gödecke, T. Über den Aufbau des Systems Eisen-Gallium z zwischen 10 und 50 At.-% Ga und dessen Abhägigkeit von der Wärmebehandlung II. Das Gleichgewichtsdiagramm. Z. Metallk. 1977, 68, 661–666. [Google Scholar]
  68. Luo, H.L. Lattice Parameters of Iron-Rich Iron-Gallium Alloys. Trans. Metall. AIME 1967, 239, 119–120. [Google Scholar]
  69. Malaman, B.; Philippe, M.J.; Roques, B.; Courtois, A.; Protas, J. Structures cristallines des phases Fe6Ge5 et Fe5Ga5. Acta Crystallogr. B 1974, 30, 2081–2087. [Google Scholar] [CrossRef] [Green Version]
  70. Meissner, H.G.; Schubert, K. Constitution of Some Systems Homologous and Quasihomologous to T5-Ga, II. The Systems Chromium-Gallium, Manganese-Gallium, and Iron-Gallium and some Notes on the Systems Vanadium-Antimony and Vanadium-Arsenic. Z. Metallk. 1965, 56, 523–530. [Google Scholar]
  71. Schubert, K.; Bhan, S.; Burkhart, W.; Gohle, R.; Meissner, H.G.; Poetschke, M.; Stolz, E. Structural Data on Metallic Phases. Naturwissenschaften 1960, 47, 303. [Google Scholar] [CrossRef]
  72. Okamoto, H. Binary Alloy Phase Diagrams, 2nd ed.; Massalski, T.B., Okamoto, H., Subramnian, P.F., Kacprzak, L., Eds.; ASM International: Materials Park, OH, USA, 1993. [Google Scholar]
  73. Ogata, T.; Nakamura, K.; Kurata, M.; Yokoo, T.; Mignanelli, M.A. Reactions between U-Pu-Zr Alloys and Fe at 923 K. J. Nucl. Sci. Technol. 2000, 37, 244–252. [Google Scholar] [CrossRef]
  74. Nakamura, K.; Ogata, T.; Kurata, M.; Yokoo, T.; Mignanelli, M.A. Reactions of Uranium-Plutonium Alloys with Iron. J. Nucl. Sci. Technol. 2001, 38, 112–119. [Google Scholar] [CrossRef]
  75. Nakamura, K.; Kurata, M.; Ogata, T.; Yokoo, T.; Mignanelli, M.A. Phase Relations in the Fe-Pu-U Ternary System. J. Phase Equilib. 2001, 22, 259–264. [Google Scholar] [CrossRef]
  76. Nakamura, K.; Ogata, T.; Kurata, M. Analysis of metal fuel/cladding metallurgical interaction during off-normal transient events with phase diagram of the U-Pu-Zr-Fe system. J. Phys. Chem. Solids 2005, 66, 643–646. [Google Scholar] [CrossRef]
  77. Dwight, A.E.; Mueller, M.H.; Conner, R.A.; Downey, J.W., Jr.; Knott, H. Ternary Compounds with the Fe2P type Structure. Trans. Metall. AIME 1968, 242, 2075–2080. [Google Scholar]
  78. Grun, Y.N.; Rogl, P.; Hiebl, K. Structural Chemistry and Magnetic Behaviour of Ternary Uranium Gallides U{Fe,Co,Ni,Ru,Rh,Pd,Os,Ir,Pt}Ga5. J. Less-Common Met. 1986, 121, 497–505. [Google Scholar] [CrossRef]
  79. Tokiwa, Y.; Maehira, T.; Ikeda, S.; Haga, Y.; Yamamoto, E.; Nakamura, A.; Onuki, Y.; Higuchi, M.; Hasegawa, A. Magnetic and Fermi Surface Properties of UFeGa5. J. Phys. Soc. Jpn. 2001, 70, 2982–2988. [Google Scholar] [CrossRef]
  80. Ikeda, S.; Tokiwa, Y.; Okubo, T.; Yamada, M.; Matsuda, T.D.; Inada, Y.; Settai, R.; Yamamoto, E.; Haga, Y.; Onuki, Y. Magnetic and Fermi surface properties of UTGa5 (T:Fe, Co and Pt). Physica B 2003, 329, 610–611. [Google Scholar] [CrossRef]
  81. Moreno, N.O.; Bauer, E.D.; Sarrao, J.L.; Hundley, M.F.; Thompson, J.D.; Fisk, Z. Thermodynamic and transport properties of single-crystalline UMGa5 (M = Fe, Co, Ni, Ru, Rh). Phys. Rev. B 2005, 72, 035119. [Google Scholar] [CrossRef] [Green Version]
  82. Gonçalves, A.P.; Werenborgh, J.C.; Sério, S.; Paixão, J.A.; Godinho, M.; Almeida, M. Structural and magnetic properties of UFe6Ga6. Intermetallics 2006, 13, 530–536. [Google Scholar] [CrossRef] [Green Version]
  83. Henriques, M.A.; Mora, P.; Cruz, M.M.; Noël, H.; Tougait, O.; Tran, V.H.; Gonçalves, A.P. Crystal structure and magnetic properties of UFe5Ga7. J. Nucl. Mater. 2009, 389, 160–163. [Google Scholar] [CrossRef]
  84. Ikeda, S.; Okubo, T.; Tokiwa, Y.; Kaneko, K.; Matsuda, T.D.; Yamamoto, E.; Haga, Y.; Onuki, Y. Magnetic properties of U2RhGa8 and U2FeGa8. J. Phys. Condens. Matter 2003, 15, S2015–S2018. [Google Scholar] [CrossRef]
  85. Jardin, R.; Colineau, E.; Griveau, J.-C.; Boulet, P.; Wastin, F.; Rebizant, J. A new family of heavy-fermion compounds. J. Alloys Compd. 2007, 432, 39–44. [Google Scholar] [CrossRef]
  86. Couderc, J.J.; Bras, J.; Fagot, F. Precipitation dans le system fer-gallium au voisinage de 25 at. % Ga. Phys. Status Solidi A 1977, 41, 595–605. [Google Scholar] [CrossRef]
  87. Nishino, Y.; Matsuo, M.; Asano, S.; Kawamiya, N. Stability of the D03 phase in (FeM1-xMx)3Ga (M = 3d transition metals). Scr. Metall. Mater. 1991, 25, 2291–2296. [Google Scholar] [CrossRef]
  88. Suzuki, T.; Oya, Y.; Ochiai, S. The Mechanical Behavior of Nonstoichiometric Compounds Ni3Si, Ni3Ge and Fe3Ga. Metall. Trans. A 1984, 15, 173–181. [Google Scholar] [CrossRef]
  89. Philippe, M.J.; Malaman, B.; Roques, B.; Courtois, A.; Protas, J. Structures cristallines des phases Fe3Ga4 et Cr3Ga4. Acta Crystallogr. B 1975, 31, 477–482. [Google Scholar] [CrossRef]
  90. Hausserman, U.; Bostrom, M.; Viklund, P.; Rapp, O.; Bjornangen, T. FeGa3 and RuGa3: Semiconducting intermetallic Compounds. J. Solid State Chem. 2002, 165, 94–99. [Google Scholar] [CrossRef]
  91. Lebech, B.; Wullf, M.; Lander, G.H.; Spirlet, J.C.; Delapalme, A. Neutron diffraction studies of the crystalline and magnetic properties of UFe2. J. Phys. Condens. Mater. 1989, 1, 10229–10248. [Google Scholar] [CrossRef]
  92. Kimball, C.W.; Vaishnava, P.P.; Jorgensen, F.Y. Phonon anomalies and local atomic displacements in the exchange-enhanced superconductor U6Fe. Phys. Rev. B 1985, 32, 4419–4425. [Google Scholar] [CrossRef] [PubMed]
  93. Ivanov, M.I.; Podol’skaya, N.S. Heats of formation of U6Fe and UFe2 Translated from. Soviet. J. Atomic. Energy. 1962, 13, 572–575. [Google Scholar]
  94. Hecker, S.S.; Timofeeva, L.F. A Tale of Two Diagrams; Los Alamos Science: Los Alamos, NM, USA, 2000. [Google Scholar]
Figure 1. Calculated Pu-Fe phase diagram based on the present CALPHAD assessment and compared with experimental data [45,46].
Figure 1. Calculated Pu-Fe phase diagram based on the present CALPHAD assessment and compared with experimental data [45,46].
Applsci 09 05040 g001
Figure 2. Calculated Pu-rich side of the Pu-Fe phase diagram based on the present CALPHAD assessment.
Figure 2. Calculated Pu-rich side of the Pu-Fe phase diagram based on the present CALPHAD assessment.
Applsci 09 05040 g002
Figure 3. U-Fe phase diagram based on the present CALPHAD assessment compared with selected experimental data [56,58,59,60,61,62,64].
Figure 3. U-Fe phase diagram based on the present CALPHAD assessment compared with selected experimental data [56,58,59,60,61,62,64].
Applsci 09 05040 g003
Figure 4. Calculated Fe-Ga phase diagram based on the present CALPHAD assessment and compared with experimental data [65,66,67].
Figure 4. Calculated Fe-Ga phase diagram based on the present CALPHAD assessment and compared with experimental data [65,66,67].
Applsci 09 05040 g004
Figure 5. Calculated Pu-U-Fe phase diagram at 650 °C based on the present CALPHAD assessment compared to experiment [75].
Figure 5. Calculated Pu-U-Fe phase diagram at 650 °C based on the present CALPHAD assessment compared to experiment [75].
Applsci 09 05040 g005
Figure 6. Calculated (a) Pu-U-Fe ternary section at room temperature with (b) a magnification of the Pu-rich part based on the present CALPHAD assessment.
Figure 6. Calculated (a) Pu-U-Fe ternary section at room temperature with (b) a magnification of the Pu-rich part based on the present CALPHAD assessment.
Applsci 09 05040 g006aApplsci 09 05040 g006b
Figure 7. Pu-Fe-Ga ternary section at 1000 °C based on the present CALPHAD assessment.
Figure 7. Pu-Fe-Ga ternary section at 1000 °C based on the present CALPHAD assessment.
Applsci 09 05040 g007
Figure 8. Calculated U-Fe-Ga ternary section at 1000 °C based on the present CALPHAD assessment.
Figure 8. Calculated U-Fe-Ga ternary section at 1000 °C based on the present CALPHAD assessment.
Applsci 09 05040 g008
Figure 9. Calculated (a) Pu-Fe-Ga and (b) U-Fe-Ga ternary sections at room temperature based on the present CALPHAD assessments.
Figure 9. Calculated (a) Pu-Fe-Ga and (b) U-Fe-Ga ternary sections at room temperature based on the present CALPHAD assessments.
Applsci 09 05040 g009
Figure 10. Property diagram of Pu89-U5-Ga5-Fe1 (at. %) calculated using the present CALPHAD database.
Figure 10. Property diagram of Pu89-U5-Ga5-Fe1 (at. %) calculated using the present CALPHAD database.
Applsci 09 05040 g010
Figure 11. Distribution of elemental components within a given phase for the Pu89-U5-Ga5-Fe1 for (a) plutonium, (b) uranium, (c) gallium and (d) iron.
Figure 11. Distribution of elemental components within a given phase for the Pu89-U5-Ga5-Fe1 for (a) plutonium, (b) uranium, (c) gallium and (d) iron.
Applsci 09 05040 g011
Figure 12. Property diagram of Pu92.9U5-Ga2-Fe0.1 (at. %) calculated using the present CALPHAD database.
Figure 12. Property diagram of Pu92.9U5-Ga2-Fe0.1 (at. %) calculated using the present CALPHAD database.
Applsci 09 05040 g012
Table 1. Crystal structure of unary systems.
Table 1. Crystal structure of unary systems.
PhaseCalphad DesignationPearson SymbolSpace GroupStrukturberichtPrototypeTransition Temp. (°C)
α-FeBCC_A2cF4Im 3 ¯ mA2W912
γ-FeFCC_A1cI2Fm 3 ¯ mA1Cu1394
δ-FeBCC_A2cF4Im 3 ¯ mA2W1538
α-GaOrthorhombic_A11oC8CmcaA11Ga30
α-PuAlpha_PumP16P21/m-α-Pu124
β-PuBeta_PumC34C2/m-β-Pu215
γ-PuGamma_PuoF8Fddd-γ-Pu320
δ-PuFCC_A1cF4Fm 3 ¯ mA1Cu463
δ-PuTetragonal_A6tI2I4/mmmA6In483
ε-PuBCC_A2cI2Im 3 ¯ mA2W640
α-UOrthorhombic_A20oC4CmcmA20α-U668
β-UTetragonal_UtP30P42/mmmAbβ-U776
γ-UBCC_A2cI2Im 3 ¯ mA2W1135
Table 2. Crystal structure of binary systems.
Table 2. Crystal structure of binary systems.
CompoundCompositionPearson SymbolSpace GroupStrukturberichtPrototypeRefs.
α’Fe0.75Ga0.25cP2Pm 3 ¯ mB2CsCl[86]
α’’Fe3GacF16Fm 3 ¯ m--[65]
α’’’Fe3GacF16Fm 3 ¯ mD03BiF3[87]
α-Fe3GaFe2.8Ga1.2cP4Pm 3 ¯ mL12AuCu3[88]
β-Fe3GaFe3GahP8P63/mmcD019Ni3Sn[87]
α-Fe6Ga5Fe6Ga5mC44C2/m-Al8Cl5[69]
β-Fe6Ga5Fe6Ga5hR26R 3 ¯ m--
Fe3Ga4Fe3Ga4mC42C2/m--[89]
FeGa3FeGa3tP16P 4 ¯ n2-CoGa3[90]
α-(Pu,U)Fe2PuFe2/UFe2cF24Fd 3 ¯ mC15Cu2Mg[45,91]
β-PuFe2PuFe2hP24P63/mmcC36MgNi2[7]
γ-PuFe2PuFe2c ** - [7]
(Pu,U)6FePu6Fe/U6Fetl28I 4 ¯ /mcmD2cMnU6[45,92]
Table 3. Crystal Structure of U-Fe-Ga intermetallic compounds.
Table 3. Crystal Structure of U-Fe-Ga intermetallic compounds.
CompoundPearson SymbolSpace GroupStrukturberichtPrototypeRefs.
UFeGahP9P62mC22Fe2P[77]
UFeGa5tP7P4/mmm-HoCoGa5[78,81]
UFe6Ga6tI26I4/mmmD2bThMn12[82]
UFe5Ga7tI26I4/mmmD2bThMn12[83]
U2FeGa8tP11P4/mmm-Ho2CoGa8[84]
U4FeGa12cI34Im 3 ¯ m-Y4PdGa12[85]
Table 4. Parameters for Pu-Fe, U-Fe, Pu-U-Fe and Fe-Ga systems.
Table 4. Parameters for Pu-Fe, U-Fe, Pu-U-Fe and Fe-Ga systems.
Phase ParameterKurata [2,3] **This Work
Liquid0LFe,Pu−35,332 + 27.530 × T−23,000 + 2.1 × T
1LFe,Pu−8149.0380
2LFe,Pu−4933.02680
0LFe,U−46,128−0.13459 × T −30,613.93−22.81 × T
1LFe,U−11,776−57,241.80 + 41.34 × T
2LFe,U9258.5−1988.06 + 7.3308 × T
0LFe,Pu,U10,000 **−14,000
1LFe,Pu,UN/A−9500
2LFe,Pu,UN/A−10,000
0LFe,GaN/A−86,500 + 18 × T
1LFe,GaN/A−15,363 + 3.5 × T
2LFe,GaN/A−13,000
BCC (ε-Pu, γ-U)0LFe,Pu13,0007295
1LFe,Pu850012,150
0LFe,U1204.553,000
1LFe,U066,000
0LFe,GaN/A−104,669 + 26.3 × T
1LFe,GaN/A8000−19 × T
FCC (δ-Pu, γ-Fe)0LFe,Pu18,00011,250
1LFe,Pu30007250
0LFe,U−3595.318,142.90
0LFe,GaN/A−107,800 + 28 × T
1LFe,GaN/A19,800-24 × T
TETRAGONAL_U (β-U)0LFe,U30,000−8400
(Pu,U)Fe2GFe:Pu−61902 + 26.18 × T + GPULIQ + 2 × GFELIQ−15850 + 0.53 × T + 0.333 × GHSERPU + 0.667 × GHSERFE
GFe:U−106,537 + 33.251 + GULIQ + 2 × GFELIQ−21,061.115-0.281944 × T + 0.333 × GHSERUU + 0.667 × GHSERFE
0LFe:Pu,U0−5000
(Pu,U)6FeGFe:Pu−91,210 + 90.6 × T + 6 × GPULIQ + GLIQFE−17,850-31.2 × T + 6 × GHSERPU + GHSERFE
GFe:U−149,660 + 88.270 × T + 6 × GPULIQ + GLIQFE−49,520.86−3.709854 × T + 6 × GHSERUU + GHSERFE
LFe:Pu,U04500
Fe3GaGFe:GaN/A−25,875 + 3.7 × T + 0.75 × GHSERFE + 0.25 × GHSERGA
Fe6Ga5GFe:GaN/A−32,594 + 3.3 × T + 0.546 × GHSERFE + 0.454 × GHSERGA
Fe3Ga4GFe:GaN/A−33,545 + 3.0 × T + 0.429 × GHSERFE + 0.571 × GHSERGA
FeGa3GFe:GaN/A−27,275 + 0.1 × T + 0.25 × GHSERFE + 0.75 × GHSERGA
** Ternary parameter of Kurata [3].
Table 5. Invariant reactions for Fe-Ga, Pu-Fe and U-Fe systems compared with experiments and models.
Table 5. Invariant reactions for Fe-Ga, Pu-Fe and U-Fe systems compared with experiments and models.
Reaction TypeReactionComposition at. % FeTemp (°C/K)Refs.
1 Congruentα’’⟷β-Fe3Ga72.572.572.5680/953[65]
757575703/976This work
2 PeritecticL + α’⟷Fe3Ga438.252/52.543906/1173[64,65,66]
39.749.342.9911/1178This work
3 PeritecticL + Fe3Ga4⟷FeGa318.54225824/1097[65,66]
Eutectic 28.442.925825/1097This work
4 PeritecticL + Fe3Ga4⟷α-Ga0.125-34/307[65,66]
0.125-34/303This work
5 Peritectoidα’ + Fe3Ga4⟷β-Fe6Ga55843.555800/1073[65,66]
54.742.954.6800/1073This work
6 Peritectoidβ-Fe6Ga5 + Fe3Ga4⟷α-Fe6Ga554.543.554.5778/1051[65,66]
----This work
7 Peritectoidβ-Fe3Ga4 + Fe6Ga5⟷α-Fe3Ga7156.570.8619/892[65,66]
----This work
8 Eutectoidβ-Fe6Ga5⟷α’ + α-Fe6Ga555.55955.5770/1043[65,66]
----This work
9 Eutectoidα’⟷β-Fe3Ga + Fe6Ga566.57156.5625/898[65,66]
67.17556.5605/883This work
10 Eutectoidβ-Fe3Ga⟷α’’ + α-Fe3Ga747573.8605/878[65,66]
----This work
11 Eutectoidα’’’⟷α-Fe + α-Fe3Ga76.679.473.7588/861[65,66]
----This work
12 CongruentL⟷PuFe266.6766.6766.671240/1513[45]
66.6766.6766.671240/1513This work
66.6766.6766.671258/1531[3]
13 PeritecticL + PuFe2⟷Pu6Fe11.566.6714.29428/701[45]
11.366.6714.29423.5/696.5This work
13.766.6714.29429.2/702.2[3]
14 CatatecticL + δ-Fe⟷γ-Fe9498.81001400/1673[44]
92.999.61001402/1675This work
92.998.61001404/1677[3]
15 CatatecticL + δ-Pu⟷ε-Pu8.51.32.4430/703[45]
8.70.92.7429.5/702.5This work
7.60.41.9444.5/716.5[3]
16 Eutectoidγ-Fe⟷α-Fe + PuFe29910066.67907/1180[46]
99.710066.67905/1178This work
98.910066.67894/1166[3]
17 EutecticL⟷γ-Fe + PuFe28298.866.71165/1438[45]
81.298.966.71164/1437This work
81.29766.71135/1408[3]
18 EutecticL⟷δ-Pu + Pu6Fe9.50.514.29413/686[46]
9.350.914.29416.6/689.6This work
8.90.414.29420/693[3]
19 CongruentL⟷UFe266.6766.6766.671230/1503[55]
66.6766.6766.671230/1503This work
66.6766.6766.671236/1509[4]
66.6766.6766.671227/1500[3]
20 PeritecticL + γ-U⟷U6Fe151.514.39832/1102[64]
15.251.3214.29832/1102This work
15.061.3514.29834/1104[4]
19.091.1414.291075.5[3]
21 Eutectoidγ-U⟷β-U + U6Fe0.80.3714.29762/1035[60]
0.850.3714.29764/1037This work
N/A0.4614.29766/1039[4]
0.70014.29763/1036[3]
22 Eutectoidβ-U⟷α-U + U6Fe0.1750.0514.29669/942[60]
0.172014.29664.4/937.4This work
0.16014.29664.8/937.8[4]
23 EutecticL⟷γ-Fe + UFe283066.671080/1353[55,56,60]
82.7066.671079/1352This work
82.55066.671078/1351[4]
83066.671080/1353[3]
24 EutecticL⟷UFe2 + U6Fe3466.6714.29723/996[60]
33.8566.6714.29723.4/996.4This work
33.4966.6714.29719.4/992.5[4]
31.1666.6714.29725.1/998.1[3]
Table 6. ΔfH (kJ/mol) of Pu-Fe and U-Fe compounds compared to DFT calculations and experiment.
Table 6. ΔfH (kJ/mol) of Pu-Fe and U-Fe compounds compared to DFT calculations and experiment.
Method/Refs.Pu6FeU6FePuFe2UFe2
CALPHAD (This work)−1.2−5.8−9.8−15
Experiment [64]/CALPHAD [4]N/A−6.8/−7.3N/A−16.3/−15.8
CALPHAD [3]−4.4−7.8−5.1−5.2
Experiment [47]N/A−2.33−13.1−11.30
Experiment [93]N/A−2.33 ± 0.72N/A−10.74
Experiment [49]−1.97 ± 0.65N/A−9.06 ± 0.56N/A
Experiment [48]N/AN/A−16.6N/A
SR-EMTO (This work)−0.50−6.17−16.25−16.0
FPLMTO (This work)0−2.5−11.6−15.5

Share and Cite

MDPI and ACS Style

Moore, E.E.; Turchi, P.E.A.; Landa, A.; Söderlind, P.; Oudot, B.; Belof, J.L.; Stout, S.A.; Perron, A. Development of a CALPHAD Thermodynamic Database for Pu-U-Fe-Ga Alloys. Appl. Sci. 2019, 9, 5040. https://doi.org/10.3390/app9235040

AMA Style

Moore EE, Turchi PEA, Landa A, Söderlind P, Oudot B, Belof JL, Stout SA, Perron A. Development of a CALPHAD Thermodynamic Database for Pu-U-Fe-Ga Alloys. Applied Sciences. 2019; 9(23):5040. https://doi.org/10.3390/app9235040

Chicago/Turabian Style

Moore, Emily E., Patrice E.A. Turchi, Alexander Landa, Per Söderlind, Benoit Oudot, Jonathan L. Belof, Stephen A. Stout, and Aurélien Perron. 2019. "Development of a CALPHAD Thermodynamic Database for Pu-U-Fe-Ga Alloys" Applied Sciences 9, no. 23: 5040. https://doi.org/10.3390/app9235040

APA Style

Moore, E. E., Turchi, P. E. A., Landa, A., Söderlind, P., Oudot, B., Belof, J. L., Stout, S. A., & Perron, A. (2019). Development of a CALPHAD Thermodynamic Database for Pu-U-Fe-Ga Alloys. Applied Sciences, 9(23), 5040. https://doi.org/10.3390/app9235040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop