# Thermodynamic assessment of the Ni–Te system

Carl-Magnus Arvhult, Christine Gueneau, Stéphane Gossé, Malin Selleby

► To cite this version:

Carl-Magnus Arvhult, Christine Gueneau, Stéphane Gossé, Malin Selleby. Thermodynamic assessment of the Ni–Te system. Journal of Materials Science, 2019, 54 (16), pp.11304-11319. 10.1007/s10853-019-03689-0. hal-03160221

HAL Id: hal-03160221

https://hal.science/hal-03160221

Submitted on 31 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Thermodynamic assessment of the Ni–Te system

Carl-Magnus Arvhult<fnref n="1" />,*, Christine Guéneau<fnref n="2" />, Stéphane Gossé<fnref n="2" />, and Malin Selleby<fnref n="1" />

**Received:** 19 January 2019

**Accepted:** 8 May 2019

**Published online:**
17 May 2019

© The Author(s) 2019

## ABSTRACT

A thermodynamic assessment of the Ni–Te system has been performed using the Calphad method, based on experimental data available in the literature. The proposed description has been developed for use in the modeling of fission-product-induced internal corrosion of stainless steel cladding in Generation IV nuclear reactors. DFT calculations were performed to obtain 0 K properties of solid phases to assist the thermodynamic optimization. The ionic liquid two-sublattice model was used, and most solution phases were modeled using interstitial metal sub-lattices. With a strict number of parameters, the resulting description satisfactorily reproduces all thermodynamic properties and high-temperature phase transitions. The metastable miscibility gap in the Ni-rich liquid that is experimentally suggested is not present in the final description. The δ phase exhibits a metastable order-disorder transition between the Cdl₂ and NiAs types of interstitial nickel distribution. The Cdl₂ prototype is the stable space group at room temperature. Low-temperature ordering phase transitions have been disregarded in this description, since they are not of interest to the application of corrosion in nuclear reactors.

## Introduction

Internal corrosion of stainless steel fuel pins for Generation IV nuclear reactors, induced by the fission products Cs and Te [1–6], might limit the lifetime of the type of reactor, and the corrosion must be modeled in order to predict its potential impact on the fuel pin integrity. Therefore, a thermodynamic database of transition-metal tellurides is underdevelopment to be incorporated into the thermodynamics of

advanced fuels-international database (TAF-ID) [7]. Descriptions of the Fe–Te and Fe–Ni systems are available in the literature [8, 9], and in order to model the Fe–Ni–Te system, Ni–Te has to be assessed and that is the topic of the present paper. DFT calculations were performed in order to get an estimate of the formation enthalpy and lattice stabilities of certain Ni–Te compounds.

<figure><img src="image_1.png" /></figure>

<footnote n="1">¹Dept. of Materials Science and Engineering, Unit of Structures, KTH Royal Institute of Technology, ITM, Brinellävgen 23, 100 44 Stockholm, Sweden</footnote>

<footnote n="2">²DEN-Service de Corrosion et du Comportement des Matériaux dans leur Environnement (SCCME), CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France</footnote>
# Thermodynamic assessment of the Ni–Te system

Carl-Magnus Arvhult<fnref n="1" />,*, Christine Guéneau<fnref n="2" />, Stéphane Gossé<fnref n="2" />, and Malin Selleby<fnref n="1" />

**Received:** 19 January 2019

**Accepted:** 8 May 2019

**Published online:**
17 May 2019

© The Author(s) 2019

## ABSTRACT

A thermodynamic assessment of the Ni–Te system has been performed using the Calphad method, based on experimental data available in the literature. The proposed description has been developed for use in the modeling of fission-product-induced internal corrosion of stainless steel cladding in Generation IV nuclear reactors. DFT calculations were performed to obtain 0 K properties of solid phases to assist the thermodynamic optimization. The ionic liquid two-sublattice model was used, and most solution phases were modeled using interstitial metal sub-lattices. With a strict number of parameters, the resulting description satisfactorily reproduces all thermodynamic properties and high-temperature phase transitions. The metastable miscibility gap in the Ni-rich liquid that is experimentally suggested is not present in the final description. The δ phase exhibits a metastable order-disorder transition between the Cdl₂ and NiAs types of interstitial nickel distribution. The Cdl₂ prototype is the stable space group at room temperature. Low-temperature ordering phase transitions have been disregarded in this description, since they are not of interest to the application of corrosion in nuclear reactors.

## Introduction

Internal corrosion of stainless steel fuel pins for Generation IV nuclear reactors, induced by the fission products Cs and Te [1–6], might limit the lifetime of the type of reactor, and the corrosion must be modeled in order to predict its potential impact on the fuel pin integrity. Therefore, a thermodynamic database of transition-metal tellurides is underdevelopment to be incorporated into the thermodynamics of

advanced fuels-international database (TAF-ID) [7]. Descriptions of the Fe–Te and Fe–Ni systems are available in the literature [8, 9], and in order to model the Fe–Ni–Te system, Ni–Te has to be assessed and that is the topic of the present paper. DFT calculations were performed in order to get an estimate of the formation enthalpy and lattice stabilities of certain Ni–Te compounds.

<figure><img src="image_2.png" /></figure>

<footnote n="1">¹Dept. of Materials Science and Engineering, Unit of Structures, KTH Royal Institute of Technology, ITM, Brinellävgen 23, 100 44 Stockholm, Sweden</footnote>

<footnote n="2">²DEN-Service de Corrosion et du Comportement des Matériaux dans leur Environnement (SCCME), CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France</footnote>
# State of the art on the Ni–Te system

## Phase diagrams and crystallography

The phase relations in the Ni–Te system have been studied by a few authors [10–15]. Lee and Nash [16] assessed the data to construct the most recent phase diagram as seen in Fig. 1. Crystallographic data have been well studied for the β2, γ1 and δ phases (Table 1) [12, 17–19].

The intermediate phases δ, β2 and γ1 consist of almost close-packed lattices of Te-atoms with Ni partially occupying interstitial sites. The β2 phase has at high temperature a Cu2Sb-type structure with additional interstitial Ni partially occupying the Ni(2) octahedral site [12, 17]. The phase has two different ordered structures at the low-temperature solubility limits: The nickel-rich structure is monoclinic, below 140°C [12], and the tellurium-rich one is orthorhombic, below 190–310 °C [20]. The close relation of these structures to the Cu2Sb-type allotrope is well described by Gulay et al. [17], who in a later paper determined the structure of the γ1 phase [19]. Several

authors have verified that the δ phase experiences an order-disorder transition at about 54.8 at.% Te [18] and a maximum in the c-axis, with the Ni-rich side being of the NiAs-type disordered structure and the Te-rich side of CdI2-type ordering of defects into every other interstitial layer [12, 13, 18, 21–23]. Barstad et al. [13] reported that there is no solubility of Te in FCC–Ni based on lattice parameter measurements, and Abakarov et al. [24] published a solubility of 0.022 at.% Ni in Te at room temperature based on the same method.

The β1 and γ2 phases decompose on quenching, hence the lack of data. The β1 phase has an ordered superstructure of doubled axes, β1', on the Ni-rich side between 731 and 793°C [16]. Stevels [20] presented in a thesis high-T XRD data on the β1 phase and its ordered allotrope β1'. With calculations, they reproduced powder patterns of the F̄43m space group with varying metal-site occupancies and got a good description with 3.8 Ni-atoms in 4(c), 1.2 in 4(b) and 1 atom in 16(e). They compared this with the β-Cu2Se phase of very similar lattice parameter of a = 5.69 Å, but that phase has later been re-characterized as

<figure><img src="image_2.png" /><figcaption>**Figure 1** A digital reproduction of The Ni–Te phase diagram assessed by Lee and Nash [16], but with updated phase names. (This figure is reproduced with permission of ASM International. All rights reserved. www.asmin international.org).</figcaption></figure>
# State of the art on the Ni–Te system

## Phase diagrams and crystallography

The phase relations in the Ni–Te system have been studied by a few authors [10–15]. Lee and Nash [16] assessed the data to construct the most recent phase diagram as seen in Fig. 1. Crystallographic data have been well studied for the β2, γ1 and δ phases (Table 1) [12, 17–19].

The intermediate phases δ, β2 and γ1 consist of almost close-packed lattices of Te-atoms with Ni partially occupying interstitial sites. The β2 phase has at high temperature a Cu2Sb-type structure with additional interstitial Ni partially occupying the Ni(2) octahedral site [12, 17]. The phase has two different ordered structures at the low-temperature solubility limits: The nickel-rich structure is monoclinic, below 140°C [12], and the tellurium-rich one is orthorhombic, below 190–310 °C [20]. The close relation of these structures to the Cu2Sb-type allotrope is well described by Gulay et al. [17], who in a later paper determined the structure of the γ1 phase [19]. Several

authors have verified that the δ phase experiences an order-disorder transition at about 54.8 at.% Te [18] and a maximum in the c-axis, with the Ni-rich side being of the NiAs-type disordered structure and the Te-rich side of CdI2-type ordering of defects into every other interstitial layer [12, 13, 18, 21–23]. Barstad et al. [13] reported that there is no solubility of Te in FCC–Ni based on lattice parameter measurements, and Abakarov et al. [24] published a solubility of 0.022 at.% Ni in Te at room temperature based on the same method.

The β1 and γ2 phases decompose on quenching, hence the lack of data. The β1 phase has an ordered superstructure of doubled axes, β1', on the Ni-rich side between 731 and 793°C [16]. Stevels [20] presented in a thesis high-T XRD data on the β1 phase and its ordered allotrope β1'. With calculations, they reproduced powder patterns of the F̄43m space group with varying metal-site occupancies and got a good description with 3.8 Ni-atoms in 4(c), 1.2 in 4(b) and 1 atom in 16(e). They compared this with the β-Cu2Se phase of very similar lattice parameter of a = 5.69 Å, but that phase has later been re-characterized as

<figure><img src="image_3.png" /><figcaption>**Figure 1** A digital reproduction of The Ni–Te phase diagram assessed by Lee and Nash [16], but with updated phase names. (This figure is reproduced with permission of ASM International. All rights reserved. www.asmin international.org).</figcaption></figure>
<table><caption>**Table 1** Crystallographic data on phases of the Ni–Te system available in the literature</caption><thead><tr><th rowspan="2">Phase</th><th rowspan="2">Pearson symbol</th><th rowspan="2">Space group</th><th colspan="4">Lattice parameters [Å]</th><th rowspan="2">Prototype</th><th rowspan="2">References</th></tr><tr><th>a</th><th>b</th><th>c</th><th>β°</th></tr></thead><tbody><tr><td>γ - Ni</td><td>cF4</td><td>Fm<span style="text-decoration: overline">3</span>m</td><td></td><td></td><td></td><td></td><td>Cu</td><td></td></tr><tr><td>β1</td><td>cF28*</td><td>Fm<span style="text-decoration: overline">3</span>m/F4<span style="text-decoration: overline">3</span>m*</td><td>5.71</td><td>—</td><td>—</td><td>—</td><td>Cu<sub>2</sub>Se*</td><td>[20]</td></tr><tr><td>β1'</td><td>N/A</td><td>N/A</td><td>11.44</td><td>—</td><td>—</td><td>—</td><td>N/A</td><td>[20]</td></tr><tr><td>β2 (ht)</td><td>tP6</td><td>P4 / nmm : 2</td><td>3.99</td><td>—</td><td>6.09</td><td>—</td><td>Cu<sub>2</sub>Sb</td><td>[12, 17]</td></tr><tr><td>β2 (hn)</td><td>mP8</td><td>P2<sub>1</sub>/m</td><td>7.54</td><td>3.79</td><td>6.09</td><td>91.16</td><td>N/A</td><td>[17]</td></tr><tr><td>β2 (ln)</td><td>oS10</td><td>Pma2</td><td>7.54</td><td>3.79</td><td>6.06</td><td>—</td><td>N/A</td><td>[17]</td></tr><tr><td rowspan="2">δ</td><td>hP3</td><td>P<span style="text-decoration: overline">3</span>/m1</td><td>3.97</td><td>—</td><td>5.36</td><td>—</td><td>CdI<sub>2</sub></td><td>[18]</td></tr><tr><td>hP4</td><td>P6<sub>3</sub>/mmc</td><td>3.86</td><td>—</td><td>5.26</td><td>—</td><td>NiAs</td><td>[18]</td></tr><tr><td>γ1</td><td>oP19</td><td>Pnma</td><td>12.38</td><td>3.92</td><td>6.88</td><td>—</td><td>Ni<sub>1.10</sub>Se<sub>0.16</sub>Te<sub>0.74</sub></td><td>[19]</td></tr><tr><td>γ2</td><td>N/A</td><td></td><td>N/A</td><td></td><td></td><td></td><td>N/A</td><td></td></tr><tr><td>Te</td><td>hP3</td><td>P3<sub>1</sub>21</td><td>4.466</td><td>—</td><td>5.919</td><td>—</td><td>γ - Se</td><td>[26]</td></tr></tbody></table>

Literature XRD data [20] consistent with a cF lattice, and given the similar lattice parameter it was compared with the β – Cu₂Se phase (a = 5.69 Å [25])

Only unique lattice parameters given, and β° only if not 90° [12, 17–20, 26]

ht, high temperature; ln, low nickel-side; hn, high nickel-side

*Not certain

Fm<span style="text-decoration: overline">3</span>m [25]. The essential difference is that in Fm<span style="text-decoration: overline">3</span>m the number of clusters of 4(c) sites each surrounded by four 16(e) sites is doubled by vertical mirroring.

The structure of the γ2 phase is still unknown.

## Thermodynamic data

Table 2 summarizes thermodynamic quantities on the phases of the Ni–Te system available in the literature. Heat capacity has been measured on the β2 and δ phases by several authors in a cryostat [27], via adiabatic shield calorimetry (AShC) [28, 29], differential scanning calorimetry (DSC) [30] and adiabatic scanning calorimetry (AScC) [31]; they are all consistent, and the β2 phase measurements show several lambda-type transitions. Enthalpy of formation has only been measured by solution calorimetry on the β2 phase by Shukla et al. [32], and supposedly on δ by Predel and Ruge [33], although their composition corresponds to γ1 + δ. Jandl et al. [34] evaluated 0 K formation enthalpies of the δ phase by DFT calculations. Activity data have been collected on most phases of the system via isopiestic measurements [14], Knudsen effusion mass spectrometry (KEMS) [35, 36] and electromotive force (EMF) [21, 37–39]. Maekawa and Yokokawa measured the partial excess

enthalpy of solution of several transition metals in liquid tellurium [40].

## DFT methodology

Total energies at 0 K were calculated via density functional theory (DFT) calculations performed using the Vienna Ab initio Simulations Package (VASP) [41–44] with the projector-augmented wave (PAW) method [45, 46]. The generalized gradient approximation (GGA) [47, 48] was used for the exchange-correlation effects. Cutoff energy and k-points were converged followed by optimization of volume by energy minimization with lattice parameter variation. Where applicable, internal degrees of freedom were optimized by shape relaxation. The pure element phases were compared with experimental data on lattice parameters and Voigt average bulk modulus, K<sub>V</sub>, as well as controlling that the DOS and final magnetic moment are reasonable. The bulk modulus was obtained by displacing every atom in every direction and calculating the elastic tensor, from which the bulk modulus is obtained as K<sub>V</sub> = [C<sub>11</sub> + C<sub>22</sub> + C<sub>33</sub> + 2(C<sub>12</sub> + C<sub>13</sub> + C<sub>23</sub>)]/9. Table 3 shows selected input and converged parameters of pure Ni and Te used in the DFT calculations, as well
<table><caption>**Table 1** Crystallographic data on phases of the Ni–Te system available in the literature</caption><thead><tr><th rowspan="2">Phase</th><th rowspan="2">Pearson symbol</th><th rowspan="2">Space group</th><th colspan="4">Lattice parameters [Å]</th><th rowspan="2">Prototype</th><th rowspan="2">References</th></tr><tr><th>a</th><th>b</th><th>c</th><th>β°</th></tr></thead><tbody><tr><td>γ - Ni</td><td>cF4</td><td>Fm<span style="text-decoration: overline">3</span>m</td><td></td><td></td><td></td><td></td><td>Cu</td><td></td></tr><tr><td>β1</td><td>cF28*</td><td>Fm<span style="text-decoration: overline">3</span>m/F4<span style="text-decoration: overline">3</span>m*</td><td>5.71</td><td>—</td><td>—</td><td>—</td><td>Cu<sub>2</sub>Se*</td><td>[20]</td></tr><tr><td>β1'</td><td>N/A</td><td>N/A</td><td>11.44</td><td>—</td><td>—</td><td>—</td><td>N/A</td><td>[20]</td></tr><tr><td>β2 (ht)</td><td>tP6</td><td>P4 / nmm : 2</td><td>3.99</td><td>—</td><td>6.09</td><td>—</td><td>Cu<sub>2</sub>Sb</td><td>[12, 17]</td></tr><tr><td>β2 (hn)</td><td>mP8</td><td>P2<sub>1</sub>/m</td><td>7.54</td><td>3.79</td><td>6.09</td><td>91.16</td><td>N/A</td><td>[17]</td></tr><tr><td>β2 (ln)</td><td>oS10</td><td>Pma2</td><td>7.54</td><td>3.79</td><td>6.06</td><td>—</td><td>N/A</td><td>[17]</td></tr><tr><td rowspan="2">δ</td><td>hP3</td><td>P<span style="text-decoration: overline">3</span>/m1</td><td>3.97</td><td>—</td><td>5.36</td><td>—</td><td>CdI<sub>2</sub></td><td>[18]</td></tr><tr><td>hP4</td><td>P6<sub>3</sub>/mmc</td><td>3.86</td><td>—</td><td>5.26</td><td>—</td><td>NiAs</td><td>[18]</td></tr><tr><td>γ1</td><td>oP19</td><td>Pnma</td><td>12.38</td><td>3.92</td><td>6.88</td><td>—</td><td>Ni<sub>1.10</sub>Se<sub>0.16</sub>Te<sub>0.74</sub></td><td>[19]</td></tr><tr><td>γ2</td><td>N/A</td><td></td><td>N/A</td><td></td><td></td><td></td><td>N/A</td><td></td></tr><tr><td>Te</td><td>hP3</td><td>P3<sub>1</sub>21</td><td>4.466</td><td>—</td><td>5.919</td><td>—</td><td>γ - Se</td><td>[26]</td></tr></tbody></table>

Literature XRD data [20] consistent with a cF lattice, and given the similar lattice parameter it was compared with the β – Cu₂Se phase (a = 5.69 Å [25])

Only unique lattice parameters given, and β° only if not 90° [12, 17–20, 26]

ht, high temperature; ln, low nickel-side; hn, high nickel-side

*Not certain

Fm<span style="text-decoration: overline">3</span>m [25]. The essential difference is that in Fm<span style="text-decoration: overline">3</span>m the number of clusters of 4(c) sites each surrounded by four 16(e) sites is doubled by vertical mirroring.

The structure of the γ2 phase is still unknown.

## Thermodynamic data

Table 2 summarizes thermodynamic quantities on the phases of the Ni–Te system available in the literature. Heat capacity has been measured on the β2 and δ phases by several authors in a cryostat [27], via adiabatic shield calorimetry (AShC) [28, 29], differential scanning calorimetry (DSC) [30] and adiabatic scanning calorimetry (AScC) [31]; they are all consistent, and the β2 phase measurements show several lambda-type transitions. Enthalpy of formation has only been measured by solution calorimetry on the β2 phase by Shukla et al. [32], and supposedly on δ by Predel and Ruge [33], although their composition corresponds to γ1 + δ. Jandl et al. [34] evaluated 0 K formation enthalpies of the δ phase by DFT calculations. Activity data have been collected on most phases of the system via isopiestic measurements [14], Knudsen effusion mass spectrometry (KEMS) [35, 36] and electromotive force (EMF) [21, 37–39]. Maekawa and Yokokawa measured the partial excess

enthalpy of solution of several transition metals in liquid tellurium [40].

## DFT methodology

Total energies at 0 K were calculated via density functional theory (DFT) calculations performed using the Vienna Ab initio Simulations Package (VASP) [41–44] with the projector-augmented wave (PAW) method [45, 46]. The generalized gradient approximation (GGA) [47, 48] was used for the exchange-correlation effects. Cutoff energy and k-points were converged followed by optimization of volume by energy minimization with lattice parameter variation. Where applicable, internal degrees of freedom were optimized by shape relaxation. The pure element phases were compared with experimental data on lattice parameters and Voigt average bulk modulus, K<sub>V</sub>, as well as controlling that the DOS and final magnetic moment are reasonable. The bulk modulus was obtained by displacing every atom in every direction and calculating the elastic tensor, from which the bulk modulus is obtained as K<sub>V</sub> = [C<sub>11</sub> + C<sub>22</sub> + C<sub>33</sub> + 2(C<sub>12</sub> + C<sub>13</sub> + C<sub>23</sub>)]/9. Table 3 shows selected input and converged parameters of pure Ni and Te used in the DFT calculations, as well
<table><caption>**Table 2** Thermodynamic data on phases of the Ni–Te system available in the literature [14, 21, 27–40]</caption><thead><tr><th>Phase(s)</th><th>Quantity</th><th>Method(s)</th><th>T [K]</th><th>at.% Te</th><th>References</th></tr></thead><tbody><tr><td rowspan="2">Liquid</td><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1170–1207</td><td>46.1–51.8</td><td>[14]</td></tr><tr><td>H<sup>E,∞</sup><sub>Ni</sub></td><td>Calorimeter</td><td>741</td><td>99.5</td><td>[40]</td></tr><tr><td rowspan="4">β2</td><td>C<sub>P</sub></td><td>AShC, DSC</td><td>218–958</td><td>40–41.2</td><td>[28, 30]</td></tr><tr><td>ΔH<sup>298</sup><sub>for</sub></td><td>Sol. Cal.</td><td>298</td><td>40.5</td><td>[32]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>963–1050</td><td>40.5</td><td>[39]</td></tr><tr><td>P<sub>Te2</sub></td><td>KEMS</td><td>981–1039</td><td>40.5</td><td>[35]</td></tr><tr><td rowspan="3">β1</td><td>μ<sub>Ni</sub></td><td>EMF</td><td>1080–1178</td><td>40.5</td><td>[39]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1106–1236</td><td>38.7–42.5</td><td>[14]</td></tr><tr><td>P<sub>Te2</sub></td><td>KEMS</td><td>1036–1150</td><td>39 and 40.5</td><td>[35]</td></tr><tr><td rowspan="4">δ</td><td>C<sub>P</sub></td><td>Cryo, AShC, AScC</td><td>5–930</td><td>52.3–66.7</td><td>[27, 29, 31]</td></tr><tr><td>ΔH<sup>0K</sup><sub>for</sub></td><td>DFT</td><td>0</td><td>50 and 66.67</td><td>[34]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>765–1178</td><td>52.5–66.3</td><td>[14]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>673 and 700</td><td>54.4–66.7</td><td>[21, 37, 38]</td></tr><tr><td rowspan="2">Various 2-phase<br>γ2 + δ</td><td>a<sub>Te</sub></td><td>Isopiestic</td><td>759–1252</td><td>34.9–70.1</td><td>[14]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1084–1126</td><td>51.2–52.4</td><td>[14]</td></tr><tr><td rowspan="2">δ + γ1</td><td>ΔH<sup>298</sup><sub>for</sub></td><td>Sol. Cal.</td><td>298</td><td>50</td><td>[33]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>673 and 700</td><td>43.5–52.6</td><td>[21, 37, 38]</td></tr><tr><td>γ – Ni + β1</td><td>P<sub>Te2</sub></td><td>KEMS</td><td>1020–1190</td><td>25–35</td><td>[36]</td></tr><tr><td>γ – Ni + β2</td><td>P<sub>Te2</sub></td><td>KEMS</td><td>893–993</td><td>25–35</td><td>[36]</td></tr></tbody></table>

as calculated unit cell volumes and bulk moduli $K_V$ compared with experiment [49, 50]. For FCC calculations, the Monkhorst pack k-point mesh was used, for HCP and the γ1 phase Γ-centered grids were used and for β2 the k-points were mapped with the automatic setting. The converged lattice parameters from 0 K DFT are very close to experimental data, as expressed via the volumes in Table 3 while the bulk moduli differ at most 26 %.

The lattice stability of HCP tellurium was evaluated, as well as the energy of formation of the compounds β2 – Ni₂Te, β2 – NiTe and γ1 – Ni₅₂Te₄₀. The energy of a compound NiₙTeₖ at 0 K was evaluated via

$$ \Delta E^{\text{Ni}_n \text{Te}_b} = E^{\text{Ni}_n \text{Te}_b} - a E_{\text{Ni}}^{\text{FCC}} - b E_{\text{Te}}^{\text{A8}} \quad (1) $$

with a and b being the stoichiometric proportions of the constituents of the phase, and the reference energies $E_{\text{Ni}}^{\text{BCC}}$ and $E_{\text{Te}}^{\text{A8}}$ were calculated with the

settings in Table 3. For a certain formation energy, the compound NiₙTeₖ and the reference states were all calculated with the same cutoff energies, selecting the highest necessary $E_{\text{cutoff}}$ of three for convergence. The energy of γ1 – Ni₅₂Te₄₀ was calculated by creating 50 supercells, each constructed of five unit cells stacked in the b-direction, and randomly removing four each of the Ni(2) and Ni(3) interstitials. An average of the total energy was taken; how this average converges with the number of sampled energies was also controlled.

## Calphad methodology

The Calphad method [52] was used to perform a thermodynamic assessment on the Ni–Te system, with solution phases modeled via the compound energy formalism (CEF) [53]. In the CEF, a phase is

<table><caption>**Table 3** DFT model input and converged parameters for nickel and tellurium. [26, 49–51]</caption><caption>Converged</caption><thead><tr><th>Element</th><th>Stable phase</th><th>Potential used</th><th>Valence e<sup>-</sup></th><th>E<sub>cutoff</sub> [eV]</th><th>Converged K-points</th><th>V<sub>rel</sub> [Å] (ref)</th><th>K<sub>V</sub> [GPa]</th><th>K<sub>exp</sub> [GPa]</th></tr></thead><tbody><tr><td>Ni</td><td>FCC_A1</td><td>Ni_sv_GW 2013 (GGA)</td><td>3d<sup>9</sup>4s<sup>1</sup></td><td>570</td><td>20 × 20 × 20</td><td>43.4 (43.8 [51])</td><td>208.3</td><td>185.3 [49]</td></tr><tr><td>Te</td><td>A8 hexagonal</td><td>GGA_GW 2012</td><td>5s<sup>2</sup>5p<sup>4</sup></td><td>330</td><td>19 × 19 × 15</td><td>104.9 (101.8 [26])</td><td>20.3</td><td>27.34 [50]</td></tr></tbody></table>
<table><caption>**Table 2** Thermodynamic data on phases of the Ni–Te system available in the literature [14, 21, 27–40]</caption><thead><tr><th>Phase(s)</th><th>Quantity</th><th>Method(s)</th><th>T [K]</th><th>at.% Te</th><th>References</th></tr></thead><tbody><tr><td rowspan="2">Liquid</td><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1170–1207</td><td>46.1–51.8</td><td>[14]</td></tr><tr><td>H<sup>E,∞</sup><sub>Ni</sub></td><td>Calorimeter</td><td>741</td><td>99.5</td><td>[40]</td></tr><tr><td rowspan="4">β2</td><td>C<sub>P</sub></td><td>AShC, DSC</td><td>218–958</td><td>40–41.2</td><td>[28, 30]</td></tr><tr><td>ΔH<sup>298</sup><sub>for</sub></td><td>Sol. Cal.</td><td>298</td><td>40.5</td><td>[32]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>963–1050</td><td>40.5</td><td>[39]</td></tr><tr><td>P<sub>Te2</sub></td><td>KEMS</td><td>981–1039</td><td>40.5</td><td>[35]</td></tr><tr><td rowspan="3">β1</td><td>μ<sub>Ni</sub></td><td>EMF</td><td>1080–1178</td><td>40.5</td><td>[39]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1106–1236</td><td>38.7–42.5</td><td>[14]</td></tr><tr><td>P<sub>Te2</sub></td><td>KEMS</td><td>1036–1150</td><td>39 and 40.5</td><td>[35]</td></tr><tr><td rowspan="4">δ</td><td>C<sub>P</sub></td><td>Cryo, AShC, AScC</td><td>5–930</td><td>52.3–66.7</td><td>[27, 29, 31]</td></tr><tr><td>ΔH<sup>0K</sup><sub>for</sub></td><td>DFT</td><td>0</td><td>50 and 66.67</td><td>[34]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>765–1178</td><td>52.5–66.3</td><td>[14]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>673 and 700</td><td>54.4–66.7</td><td>[21, 37, 38]</td></tr><tr><td rowspan="2">Various 2-phase<br>γ2 + δ</td><td>a<sub>Te</sub></td><td>Isopiestic</td><td>759–1252</td><td>34.9–70.1</td><td>[14]</td></tr><tr><td>a<sub>Te</sub></td><td>Isopiestic</td><td>1084–1126</td><td>51.2–52.4</td><td>[14]</td></tr><tr><td rowspan="2">δ + γ1</td><td>ΔH<sup>298</sup><sub>for</sub></td><td>Sol. Cal.</td><td>298</td><td>50</td><td>[33]</td></tr><tr><td>μ<sub>Ni</sub></td><td>EMF</td><td>673 and 700</td><td>43.5–52.6</td><td>[21, 37, 38]</td></tr><tr><td>γ – Ni + β1</td><td>P<sub>Te2</sub></td><td>KEMS</td><td>1020–1190</td><td>25–35</td><td>[36]</td></tr><tr><td>γ – Ni + β2</td><td>P<sub>Te2</sub></td><td>KEMS</td><td>893–993</td><td>25–35</td><td>[36]</td></tr></tbody></table>

as calculated unit cell volumes and bulk moduli $K_V$ compared with experiment [49, 50]. For FCC calculations, the Monkhorst pack k-point mesh was used, for HCP and the γ1 phase Γ-centered grids were used and for β2 the k-points were mapped with the automatic setting. The converged lattice parameters from 0 K DFT are very close to experimental data, as expressed via the volumes in Table 3 while the bulk moduli differ at most 26 %.

The lattice stability of HCP tellurium was evaluated, as well as the energy of formation of the compounds β2 – Ni₂Te, β2 – NiTe and γ1 – Ni₅₂Te₄₀. The energy of a compound NiₙTeₖ at 0 K was evaluated via

$$ \Delta E^{\text{Ni}_n \text{Te}_b} = E^{\text{Ni}_n \text{Te}_b} - a E_{\text{Ni}}^{\text{FCC}} - b E_{\text{Te}}^{\text{A8}} \quad (1) $$

with a and b being the stoichiometric proportions of the constituents of the phase, and the reference energies $E_{\text{Ni}}^{\text{BCC}}$ and $E_{\text{Te}}^{\text{A8}}$ were calculated with the

settings in Table 3. For a certain formation energy, the compound NiₙTeₖ and the reference states were all calculated with the same cutoff energies, selecting the highest necessary $E_{\text{cutoff}}$ of three for convergence. The energy of γ1 – Ni₅₂Te₄₀ was calculated by creating 50 supercells, each constructed of five unit cells stacked in the b-direction, and randomly removing four each of the Ni(2) and Ni(3) interstitials. An average of the total energy was taken; how this average converges with the number of sampled energies was also controlled.

## Calphad methodology

The Calphad method [52] was used to perform a thermodynamic assessment on the Ni–Te system, with solution phases modeled via the compound energy formalism (CEF) [53]. In the CEF, a phase is

<table><caption>**Table 3** DFT model input and converged parameters for nickel and tellurium. [26, 49–51]</caption><caption>Converged</caption><thead><tr><th>Element</th><th>Stable phase</th><th>Potential used</th><th>Valence e<sup>-</sup></th><th>E<sub>cutoff</sub> [eV]</th><th>Converged K-points</th><th>V<sub>rel</sub> [Å] (ref)</th><th>K<sub>V</sub> [GPa]</th><th>K<sub>exp</sub> [GPa]</th></tr></thead><tbody><tr><td>Ni</td><td>FCC_A1</td><td>Ni_sv_GW 2013 (GGA)</td><td>3d<sup>9</sup>4s<sup>1</sup></td><td>570</td><td>20 × 20 × 20</td><td>43.4 (43.8 [51])</td><td>208.3</td><td>185.3 [49]</td></tr><tr><td>Te</td><td>A8 hexagonal</td><td>GGA_GW 2012</td><td>5s<sup>2</sup>5p<sup>4</sup></td><td>330</td><td>19 × 19 × 15</td><td>104.9 (101.8 [26])</td><td>20.3</td><td>27.34 [50]</td></tr></tbody></table>
<table><caption>**Table 4** Sublattice models used for alloy phases of the Ni–Te system</caption><thead><tr><th>Phase</th><th>Sublattice model</th><th>End-members</th><th>Range at.%<br>Te</th><th>Comment</th></tr></thead><tbody><tr><td>Liquid</td><td>(Ni<sup>+2</sup>)P(Te<sup>-2</sup>, Va<sup>-Q</sup>, Te<sup>0</sup>)Q</td><td>Ni, Ni<sub>2</sub>Te<sub>2</sub>,Te</td><td>0–100</td><td>Ionic liquid</td></tr><tr><td>BCC</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>3</sub></td><td>Ni<sub>1</sub>Va<sub>3</sub>, Te<sub>1</sub>Va<sub>3</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>FCC</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>1</sub></td><td>Ni<sub>1</sub>Va<sub>1</sub>, Te<sub>1</sub>Va<sub>1</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>HCP</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>0.5</sub></td><td>Ni<sub>1</sub>Va<sub>0.5</sub>, Te<sub>1</sub>Va<sub>0.5</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>β1</td><td>(Ni, Te)<sub>2</sub>(Te)<sub>1</sub></td><td>Ni<sub>2</sub>Te<sub>1</sub>, Te<sub>2</sub>Te<sub>1</sub></td><td>33.33–100</td><td></td></tr><tr><td>β2</td><td>(Ni)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>1</sub></td><td>Ni<sub>1</sub>Ni<sub>1</sub>Te<sub>1</sub>, Ni<sub>1</sub>Va<sub>1</sub>Te<sub>1</sub></td><td>33.33–50</td><td></td></tr><tr><td>γ1</td><td>(Ni)<sub>52</sub>(Te)<sub>40</sub></td><td>Ni<sub>52</sub>Te<sub>40</sub></td><td>43.48</td><td></td></tr><tr><td>γ2</td><td>(Ni)<sub>20</sub>(Te)<sub>17</sub></td><td>Ni<sub>20</sub>Te<sub>17</sub></td><td>45.95</td><td></td></tr><tr><td>δ</td><td>(Ni, Va)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>2</sub></td><td>Ni<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>, Ni<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub> Va<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>,<br>Va<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub></td><td>50–100</td><td>Ni<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub> and Va<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>, are<br>equivalent</td></tr></tbody></table>

divided into sublattices and a Gibbs energy surface stretched between end-members, stable or hypothetical (metastable) compounds comprising all possible combinations of the sublattices each filled with one constituent. These are the compounds for which the total energy was computed via DFT; even though some end-members are not stable, in the Calphad method they must have a free energy assigned. If a formation energy has not been evaluated, it is common to use 5 kJ/mol per atom [52].

$$G_m^\alpha = \text{srf} G_m^\alpha + \text{cnf} G_m^\alpha + \text{E} G_m^\alpha + \text{phys} G_m^\alpha \quad (2)$$

$$\text{srf} G_m^\alpha = \sum_{I_0} P_{I_0}(Y) G_{I_0} \quad (3)$$

$$\text{cnf} G_m^\alpha = RT \sum_{s=1}^{n} a_s \sum_{i=1}^{n_s} y_i^{(s)} \ln(y_i^{(s)}) \quad (4)$$

$$\text{E} G_m^\alpha = \sum_{I_1} P_{I_1}(Y) L_{I_1} + \sum_{I_2} P_{I_2}(Y) L_{I_2} + \dots \quad (5)$$

where $\text{srf} G_m^\alpha$ is the reference energy surface contribution from theoretical mechanical mixing of the constituents, with $I_0$ being a zeroth-order constituent array describing an end-member with a constituent in each sublattice, $P_{I_0}$ the product of the fractions of those constituents and $G_{I_0}$ the Gibbs energy of formation of that end-member compound. The configurational entropy of random mixing gives a contribution via $\text{cnf} G_m^\alpha$, the contribution from interaction between constituents on the same sublattice is represented by the excess energy $\text{E} G_m^\alpha$ and other physical contributions, such as magnetism, are usually assigned to $\text{phys} G_m^\alpha$. For a more thorough

explanation of the equations, see relevant literature [52, 53].

## Solution phases

Table 4 summarizes the sublattice models of all intermediate phases used in this assessment of the Ni–Te system, with their respective end-members and composition ranges.

For simplicity, the $\beta1$ phase was modeled with two sublattices as $(Ni, Te)_2(Te)_1$, since the crystal structure is not well known. The 2:1 ratio was chosen to create an end-member close to the Ni-rich phase boundary. The order–disorder transformation into $\beta1'$ at lower temperature and higher Ni content is disregarded in this work. The $\beta2$ phase was modeled as $(Ni)_1(Te)_1(Ni, Va)_1$ according to the tetragonal P4 / nmm space group, compatible to the Fe–Te $\beta$ phase modeled in a previous assessment [8]. The low-temperature ordering into space groups P21 / m (monoclinic) and Pma2 (orthorhombic) at high and low Ni content, respectively, is not modeled here since it is not of particular interest to the application of nuclear reactors, which operate at higher temperatures. The $\delta$ phase was modeled as $(Ni, Va)_1(Ni, Va)_1(Te)_2$ in order to facilitate the second-order transition between the NiAs type disordering of vacancies to the CdI$_2$-type layering of vacancies in every other Ni layer (hence the separation of the (Ni,Va) interstitial sites into two sublattices) as has been experimentally verified [18, 22]. This is how the $\delta$ phase was previously modeled in the Fe–Te system [8], and compatibility
<table><caption>**Table 4** Sublattice models used for alloy phases of the Ni–Te system</caption><thead><tr><th>Phase</th><th>Sublattice model</th><th>End-members</th><th>Range at.%<br>Te</th><th>Comment</th></tr></thead><tbody><tr><td>Liquid</td><td>(Ni<sup>+2</sup>)P(Te<sup>-2</sup>, Va<sup>-Q</sup>, Te<sup>0</sup>)Q</td><td>Ni, Ni<sub>2</sub>Te<sub>2</sub>,Te</td><td>0–100</td><td>Ionic liquid</td></tr><tr><td>BCC</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>3</sub></td><td>Ni<sub>1</sub>Va<sub>3</sub>, Te<sub>1</sub>Va<sub>3</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>FCC</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>1</sub></td><td>Ni<sub>1</sub>Va<sub>1</sub>, Te<sub>1</sub>Va<sub>1</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>HCP</td><td>(Ni, Te)<sub>1</sub>(Va)<sub>0.5</sub></td><td>Ni<sub>1</sub>Va<sub>0.5</sub>, Te<sub>1</sub>Va<sub>0.5</sub></td><td>0–100</td><td>Ideal solution</td></tr><tr><td>β1</td><td>(Ni, Te)<sub>2</sub>(Te)<sub>1</sub></td><td>Ni<sub>2</sub>Te<sub>1</sub>, Te<sub>2</sub>Te<sub>1</sub></td><td>33.33–100</td><td></td></tr><tr><td>β2</td><td>(Ni)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>1</sub></td><td>Ni<sub>1</sub>Ni<sub>1</sub>Te<sub>1</sub>, Ni<sub>1</sub>Va<sub>1</sub>Te<sub>1</sub></td><td>33.33–50</td><td></td></tr><tr><td>γ1</td><td>(Ni)<sub>52</sub>(Te)<sub>40</sub></td><td>Ni<sub>52</sub>Te<sub>40</sub></td><td>43.48</td><td></td></tr><tr><td>γ2</td><td>(Ni)<sub>20</sub>(Te)<sub>17</sub></td><td>Ni<sub>20</sub>Te<sub>17</sub></td><td>45.95</td><td></td></tr><tr><td>δ</td><td>(Ni, Va)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>2</sub></td><td>Ni<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>, Ni<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub> Va<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>,<br>Va<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub></td><td>50–100</td><td>Ni<sub>1</sub>Va<sub>1</sub>Te<sub>2</sub> and Va<sub>1</sub>Ni<sub>1</sub>Te<sub>2</sub>, are<br>equivalent</td></tr></tbody></table>

divided into sublattices and a Gibbs energy surface stretched between end-members, stable or hypothetical (metastable) compounds comprising all possible combinations of the sublattices each filled with one constituent. These are the compounds for which the total energy was computed via DFT; even though some end-members are not stable, in the Calphad method they must have a free energy assigned. If a formation energy has not been evaluated, it is common to use 5 kJ/mol per atom [52].

$$G_m^\alpha = \text{srf} G_m^\alpha + \text{cnf} G_m^\alpha + \text{E} G_m^\alpha + \text{phys} G_m^\alpha \quad (2)$$

$$\text{srf} G_m^\alpha = \sum_{I_0} P_{I_0}(Y)^{\circ} G_{I_0} \quad (3)$$

$$\text{cnf} G_m^\alpha = RT \sum_{s=1}^{n} a_s \sum_{i=1}^{n_s} y_i^{(s)} \ln(y_i^{(s)}) \quad (4)$$

$$\text{E} G_m^\alpha = \sum_{I_1} P_{I_1}(Y)L_{I_1} + \sum_{I_2} P_{I_2}(Y)L_{I_2} + \dots \quad (5)$$

where $\text{srf} G_m^\alpha$ is the reference energy surface contribution from theoretical mechanical mixing of the constituents, with $I_0$ being a zeroth-order constituent array describing an end-member with a constituent in each sublattice, $P_{I_0}$ the product of the fractions of those constituents and $\circ G_{I_0}$ the Gibbs energy of formation of that end-member compound. The configurational entropy of random mixing gives a contribution via $\text{cnf} G_m^\alpha$, the contribution from interaction between constituents on the same sublattice is represented by the excess energy $\text{E} G_m^\alpha$ and other physical contributions, such as magnetism, are usually assigned to $\text{phys} G_m^\alpha$. For a more thorough

explanation of the equations, see relevant literature [52, 53].

## Solution phases

Table 4 summarizes the sublattice models of all intermediate phases used in this assessment of the Ni–Te system, with their respective end-members and composition ranges.

For simplicity, the $\beta1$ phase was modeled with two sublattices as $(Ni, Te)_2(Te)_1$, since the crystal structure is not well known. The 2:1 ratio was chosen to create an end-member close to the Ni-rich phase boundary. The order–disorder transformation into $\beta1'$ at lower temperature and higher Ni content is disregarded in this work. The $\beta2$ phase was modeled as $(Ni)_1(Te)_1(Ni, Va)_1$ according to the tetragonal P4 / nmm space group, compatible to the Fe–Te $\beta$ phase modeled in a previous assessment [8]. The low-temperature ordering into space groups P21 / m (monoclinic) and Pma2 (orthorhombic) at high and low Ni content, respectively, is not modeled here since it is not of particular interest to the application of nuclear reactors, which operate at higher temperatures. The $\delta$ phase was modeled as $(Ni, Va)_1(Ni, Va)_1(Te)_2$ in order to facilitate the second-order transition between the NiAs type disordering of vacancies to the CdI$_2$-type layering of vacancies in every other Ni layer (hence the separation of the (Ni,Va) interstitial sites into two sublattices) as has been experimentally verified [18, 22]. This is how the $\delta$ phase was previously modeled in the Fe–Te system [8], and compatibility
was accommodated for complete exchange of Fe and Ni in the Fe–Ni–Te system.

With experimental Cₚ data available in the literature, Gibbs energy for a generic end-member AᵢBⱼ can be modeled with a power series with temperature as

$$ 
\begin{aligned}
\Delta^{\circ} G_{\mathrm{m}}^{A_{i}B_{j}} &=^{\circ}G_{\mathrm{m}}^{A_{i}B_{j}}-i^{\circ}H_{\mathrm{A}}^{\mathrm{SER}}-j^{\circ}H_{\mathrm{B}}^{\mathrm{SER}} \\
&=a+bT+cT \ln (T)+d T^{2}+e T^{-1}+f T^{3}
\end{aligned}
\quad (6) $$

where H SER m is the standard element reference, i.e., the enthalpy of the stable state at 1 bar of pressure at room temperature. If Cₚ data are not available for a phase, it is more convenient to use the Neumann–Kopp rule (NKR) which assumes a linear combination of the contributions from the Cₚ of pure constituents, usually expressed via the Gibbs energy as:

$$ \Delta^{\circ} G_{\mathrm{m}}^{A_{i}B_{j}}-(i^{\circ} G_{\mathrm{A}}+j^{\circ} G_{\mathrm{B}})=a+bT \quad (7) $$

## Liquid phase

Fitting the eutectic point at around 50 at.% Te together with a high Ni-rich liquidus proved difficult with a substitutional model for the liquid phase. It was therefore decided to use the ionic liquid model as (Ni⁺²)ₚ(Te⁻², Va⁻², Te⁰)ₚ, which is in accordance with published findings showing short-range order (SRO) of ions in the liquid [54]. Here, P and Q, respectively, equal the average charge of the opposite sublattice; hence, Q is always 2 and P varies as P = 2y''_{Te-2} + Qy''_{Va}. This model is compatible with the ionic liquid used in the Fe–Te system [8] and will accommodate addition of oxygen. The neutral Te⁰ species in the anionic sublattice makes it possible to cover the whole composition range; the charged vacancy maintains charge neutrality. This model creates an end-member at 50 at.% Te, making it easier to control the shape of the liquidus separately on both sides of the central eutecticum. This model can be made equivalent to a substitutional solution model with an associate as (Ni,NiTe,Te). Equations 3, 4 and 5 are for the ionic liquid changed into

$$ 
\begin{aligned}
\mathrm{srf} G_{\mathrm{m}}^{\mathrm{liq}} &= y'_{\mathrm{Ni}^{+2}} y''_{\mathrm{Te}^{-2}}^{\circ} G_{\mathrm{Ni}^{+2}:\mathrm{Te}^{-2}} \\
&\quad + Q(y'_{\mathrm{Ni}^{+2}} y''_{\mathrm{Va}^{-2}}^{\circ} G_{\mathrm{Ni}^{+2}:\mathrm{Va}^{-2}} + y''_{\mathrm{Te}^{0}}^{\circ} G_{\mathrm{Te}^{0}})
\end{aligned}
\quad (8) $$

$$ 
\begin{aligned}
\mathrm{cnf} G_{\mathrm{m}}^{\mathrm{liq}} &= \mathrm{QRT}[y''_{\mathrm{Va}^{-2}} \ln(y''_{\mathrm{Va}^{-2}}) + y''_{\mathrm{Te}^{-2}} \ln(y''_{\mathrm{Te}^{-2}}) \\
&\quad + y''_{\mathrm{Te}^{0}} \ln(y''_{\mathrm{Te}^{0}})]
\end{aligned}
\quad (9) $$

$$ 
\begin{aligned}
E G_{\mathrm{m}}^{\mathrm{liq}} = & Q y'_{\mathrm{Ni}^{+2}} y''_{\mathrm{Va}^{-2}} y''_{\mathrm{Te}^{0}} L_{\mathrm{Ni}^{+2}:\mathrm{Va}^{-2},\mathrm{Te}^{0}} \\
& + y'_{\mathrm{Ni}^{+2}} y''_{\mathrm{Va}^{-2}} y''_{\mathrm{Te}^{-2}} L_{\mathrm{Ni}^{+2}:\mathrm{Te}^{-2},\mathrm{Va}^{-2}} \\
& + y'_{\mathrm{Ni}^{+2}} y''_{\mathrm{Te}^{-2}} y''_{\mathrm{Te}^{0}} L_{\mathrm{Ni}^{+2}:\mathrm{Te}^{-2},\mathrm{Te}^{0}}
\end{aligned}
\quad (10) $$

where y'_{Ni+2} = 1 as there is only one cation; it is retained in the equations above for clarity.

## Stoichiometric phases

The stoichiometric γ1 phase was modeled as (Ni)₅₂(Te)₄₀ and the γ2 phase as (Ni)₂₀(Te)₁₇, and the NKR was used for both. The pure Te phase (hexagonal A8) and γ–Ni (FCC), α–Ni (BCC) and ε–Ni (HCP) were taken directly from Dinsdale’s descriptions in the PURE5 database from SGTE [55]. No nickel was added to the A8 phase, but Te was added to the Ni sites of FCC, BCC and HCP. HCP was included in the assessment only to enter the tellurium lattice stability evaluated via DFT.

## Optimization procedure

A thermodynamic assessment was performed using the PARROT module of the Thermo-Calc Software package [56] for parameter optimization. The assessment began by using a substitutional solution model for the liquid, with parameters of appropriate order of magnitude. Thereafter, solid phases were introduced one by one and optimized in the following manner.

### Liquid, gas and terminal phases

The gas phase description was not optimized in this work, but merely extracted from the SGTE SSUB5 database [57–60].

After solid phases had been introduced, liquid model parameters were optimized to fit the congruent points of the β1 and δ phases, as well as a tentative liquidus point representing the “melting effect” found at 1663 K for 15 at.% Te by Klepp and Komarek [15]. There were no further details about that data point in their paper. When those solid phases were approximately in place, the interaction terms in the liquid could be further optimized to fit liquidus and activity data [11, 14, 15], and the liquid Ni₂Te₂ end-member was used to help fit the eutectic point between γ2 and δ.
was accommodated for complete exchange of Fe and Ni in the Fe–Ni–Te system.

With experimental Cₚ data available in the literature, Gibbs energy for a generic end-member AᵢBⱼ can be modeled with a power series with temperature as

$$ \begin{aligned} \Delta^{\circ}G_{m}^{A_{i}B_{j}} &= \Delta^{\circ}G_{m}^{A_{i}B_{j}} - i^{\circ}H_{A}^{\text{SER}} - j^{\circ}H_{B}^{\text{SER}} \\ &= a + bT + cT \ln(T) + dT^{2} + eT^{-1} + fT^{3} \end{aligned} \quad (6) $$

where HSER is the standard element reference, i.e., the enthalpy of the stable state at 1 bar of pressure at room temperature. If Cₚ data are not available for a phase, it is more convenient to use the Neumann–Kopp rule (NKR) which assumes a linear combination of the contributions from the Cₚ of pure constituents, usually expressed via the Gibbs energy as:

$$ \Delta^{\circ}G_{m}^{A_{i}B_{j}} - (i^{\circ}G_{A} + j^{\circ}G_{B}) = a + bT \quad (7) $$

## Liquid phase

Fitting the eutectic point at around 50 at.% Te together with a high Ni-rich liquidus proved difficult with a substitutional model for the liquid phase. It was therefore decided to use the ionic liquid model as (Ni⁺²)ₚ(Te⁻², Va⁻², Te⁰)Q, which is in accordance with published findings showing short-range order (SRO) of ions in the liquid [54]. Here, P and Q, respectively, equal the average charge of the opposite sublattice; hence, Q is always 2 and P varies as P = 2y''_{Te-2} + Qy''_{Va}. This model is compatible with the ionic liquid used in the Fe–Te system [8] and will accommodate addition of oxygen. The neutral Te⁰ species in the anionic sublattice makes it possible to cover the whole composition range; the charged vacancy maintains charge neutrality. This model creates an end-member at 50 at.% Te, making it easier to control the shape of the liquidus separately on both sides of the central eutecticum. This model can be made equivalent to a substitutional solution model with an associate as (Ni,NiTe,Te). Equations 3, 4 and 5 are for the ionic liquid changed into

$$ \begin{aligned} \text{srf} G_{m}^{\text{liq}} &= y'_{\text{Ni}^{+2}} y''_{\text{Te}^{-2}} \,^{\circ}G_{\text{Ni}^{+2}:\text{Te}^{-2}} \\ &\quad + Q(y'_{\text{Ni}^{+2}} y''_{\text{Va}^{-2}} \,^{\circ}G_{\text{Ni}^{+2}:\text{Va}^{-2}} + y''_{\text{Te}^{0}} \,^{\circ}G_{\text{Te}^{0}}) \end{aligned} \quad (8) $$

$$ \begin{aligned} \text{cnf} G_{m}^{\text{liq}} &= QRT[y''_{\text{Va}^{-2}} \ln(y''_{\text{Va}^{-2}}) + y''_{\text{Te}^{-2}} \ln(y''_{\text{Te}^{-2}}) \\ &\quad + y''_{\text{Te}^{0}} \ln(y''_{\text{Te}^{0}})] \end{aligned} \quad (9) $$

$$ \begin{aligned} E G_{m}^{\text{liq}} = & Q y'_{\text{Ni}^{+2}} y''_{\text{Va}^{-2}} y''_{\text{Te}^{0}} L_{\text{Ni}^{+2}:\text{Va}^{-2},\text{Te}^{0}} \\ & + y'_{\text{Ni}^{+2}} y''_{\text{Va}^{-2}} y''_{\text{Te}^{-2}} L_{\text{Ni}^{+2}:\text{Te}^{-2},\text{Va}^{-2}} \\ & + y'_{\text{Ni}^{+2}} y''_{\text{Te}^{-2}} y''_{\text{Te}^{0}} L_{\text{Ni}^{+2}:\text{Te}^{-2},\text{Te}^{0}} \end{aligned} \quad (10) $$

where y'_{Ni+2} = 1 as there is only one cation; it is retained in the equations above for clarity.

## Stoichiometric phases

The stoichiometric γ1 phase was modeled as (Ni)₅₂(Te)₄₀ and the γ2 phase as (Ni)₂₀(Te)₁₇, and the NKR was used for both. The pure Te phase (hexagonal A8) and γ–Ni (FCC), α–Ni (BCC) and ε–Ni (HCP) were taken directly from Dinsdale’s descriptions in the PURE5 database from SGTE [55]. No nickel was added to the A8 phase, but Te was added to the Ni sites of FCC, BCC and HCP. HCP was included in the assessment only to enter the tellurium lattice stability evaluated via DFT.

## Optimization procedure

A thermodynamic assessment was performed using the PARROT module of the Thermo-Calc Software package [56] for parameter optimization. The assessment began by using a substitutional solution model for the liquid, with parameters of appropriate order of magnitude. Thereafter, solid phases were introduced one by one and optimized in the following manner.

### Liquid, gas and terminal phases

The gas phase description was not optimized in this work, but merely extracted from the SGTE SSUB5 database [57–60].

After solid phases had been introduced, liquid model parameters were optimized to fit the congruent points of the β1 and δ phases, as well as a tentative liquidus point representing the “melting effect” found at 1663 K for 15 at.% Te by Klepp and Komarek [15]. There were no further details about that data point in their paper. When those solid phases were approximately in place, the interaction terms in the liquid could be further optimized to fit liquidus and activity data [11, 14, 15], and the liquid Ni₂Te₂ end-member was used to help fit the eutectic point between γ2 and δ.
<table><caption>**Table 5** Relaxed lattice parameters and 0 K formation energies (ΔEf) from DFT calculations of Ni–Te end-members. [17]</caption><thead><tr><th rowspan="2">Phase</th><th rowspan="2">Ecutoff<br>[eV]</th><th rowspan="2">k-points</th><th colspan="3">Relaxed lattice parameters [Å]</th><th rowspan="2">ΔEf<br>[kJ/mol]</th></tr><tr><th>a</th><th>b</th><th>c</th></tr></thead><tbody><tr><td>β2 – Ni<sub>2</sub>Te</td><td>570</td><td>15 × 15 × 10</td><td>3.88</td><td>—</td><td>6.02</td><td>— 7.8</td></tr><tr><td>β2 – NiTe</td><td>570</td><td>15 × 15 × 10</td><td>3.75</td><td>—</td><td>6.03</td><td>— 16.8</td></tr><tr><td>Exp. β2 – Ni [17]</td><td></td><td></td><td>3.78</td><td></td><td>6.06</td><td></td></tr><tr><td>γ1 – Ni<sub>52</sub>Te<sub>40</sub></td><td>430</td><td>3 × 2 × 6</td><td>Volume not relaxed</td><td></td><td></td><td>— 16.7</td></tr><tr><td>HCP_A3 – Te</td><td>370</td><td>31 × 31 × 19</td><td>4.06</td><td>—</td><td>4.39</td><td>17.3</td></tr></tbody></table>

No solubility was modeled in the pure terminal phases, although a small amount of Ni might be soluble in pure Te [24].

## β1 and δ phases

The end-member and interaction parameters of the β1 phase were first manually adjusted in order to get the phase in the approximately correct location in the phase diagram. The parameters were then optimized to fit the $L → FCC + β1$ eutecticum and activity data [14, 35, 36, 39]

All terms except the entropy terms ($bT$ in Eq. 6) of the $δ–NiTe$ and $δ – NiTe_2$ end-members were fitted to the formation enthalpy data evaluated via DFT by Jandl et al. [34] and the heat capacity data of Grønvold [29] and Tsuji and Ishida [31]. The interaction terms were then optimized to fit activity data [14, 21, 38] and the $bT$ terms to fit solidus data [14, 15] and the phase boundaries evaluated by Barstad et al. [13]. The interaction term $^0L_{Ni,Va:Va:Te}^δ = ^0L_{Va:Ni,Va:Te}^δ$, affecting the range $NiTe_2–Te$, was given a constant positive value in order to suppress a higher Te-solubility than about 66.67 at.% Te.

## β2, γ1 and γ2 phases

The end-members of the β2 phase were first manually adjusted to make the phase appear in approximately the correct temperature and composition range. Both end-members, i.e., β2 – Ni<sub>2</sub>Te and β2–NiTe, were modeled using the NKR. The stable composition range of the phase is narrow and lies roughly in the middle between the end-members; the end-members are thus rather far away from the compositions of available heat capacity data [29], which made it difficult to fit the data with end-members described as power series (Eq. 6). Instead, it worked well to model it as a power-series contribution to the interaction parameter $^0L_{Ni:Ni,Va:Te}^{β2}$. The a- and b-terms of the

interaction parameters together with the b-term of β2 – Ni<sub>2</sub>Te could then be optimized to fit activity data [35, 39] and the congruent reaction $β2 ⇒ β1$. The a-terms of the end-members were optimized to fit a sort of compromise between the enthalpy of formation of the end-members evaluated via DFT in this work and the calorimetric point by Shukla et al. [32]. If the latter enthalpy of formation were to be fit well, the phase would be too stable to allow the proper shape of the δ phase.

The a-term of the γ1 phase was first fitted to the DFT enthalpy of formation value evaluated in this work, but was then lowered to fit the phase diagram; with a formation enthalpy of the β2 phase so low, the γ1 phase enthalpy must also be rather low to form the Gibbs energy tangent between the β2 and δ necessary to have them all stable at low T. The b-term was fitted to relevant invariant reactions. The γ2 phase was added last, and its parameters were optimized to fit the phase diagram only, since there are virtually no other data available on the phase.

## DFT Results and discussion

The formation energies of end-members of Ni–Te phases and lattice stability of HCP_A3 tellurium are summarized in Table 5 together with converged parameters and lattice parameters compared with experimental values. For the γ1 phase calculations, full volume relaxation was not performed since it was desired to only affect the energy by removing atoms, without allowing the cell shape or atom movement to compensate. Out of the 50 different cells, the one with the largest total energy had the vacancies clustered close together, and one with the lowest energy had the vacancies spread over the supercell. Figure 2 shows how the average total energy of γ1 – Ni<sub>52</sub>Te<sub>40</sub> converges with the number of
<table><caption>**Table 5** Relaxed lattice parameters and 0 K formation energies (ΔE<sub>f</sub>) from DFT calculations of Ni–Te end-members. [17]</caption><thead><tr><th rowspan="2">Phase</th><th rowspan="2">Ecutoff<br>[eV]</th><th rowspan="2">k-points</th><th colspan="3">Relaxed lattice parameters [Å]</th><th rowspan="2">ΔE<sub>f</sub><br>[kJ/mol]</th></tr><tr><th>a</th><th>b</th><th>c</th></tr></thead><tbody><tr><td>β2 – Ni<sub>2</sub>Te</td><td>570</td><td>15 × 15 × 10</td><td>3.88</td><td>—</td><td>6.02</td><td>— 7.8</td></tr><tr><td>β2 – NiTe</td><td>570</td><td>15 × 15 × 10</td><td>3.75</td><td>—</td><td>6.03</td><td>— 16.8</td></tr><tr><td>Exp. β2 – Ni [17]</td><td></td><td></td><td>3.78</td><td></td><td>6.06</td><td></td></tr><tr><td>γ1 – Ni<sub>52</sub>Te<sub>40</sub></td><td>430</td><td>3 × 2 × 6</td><td>Volume not relaxed</td><td></td><td></td><td>— 16.7</td></tr><tr><td>HCP_A3 – Te</td><td>370</td><td>31 × 31 × 19</td><td>4.06</td><td>—</td><td>4.39</td><td>17.3</td></tr></tbody></table>

No solubility was modeled in the pure terminal phases, although a small amount of Ni might be soluble in pure Te [24].

## β1 and δ phases

The end-member and interaction parameters of the β1 phase were first manually adjusted in order to get the phase in the approximately correct location in the phase diagram. The parameters were then optimized to fit the $L → FCC + β1$ eutecticum and activity data [14, 35, 36, 39]

All terms except the entropy terms (bT in Eq. 6) of the δ–NiTe and δ – NiTe<sub>2</sub> end-members were fitted to the formation enthalpy data evaluated via DFT by Jandl et al. [34] and the heat capacity data of Grønvold [29] and Tsuji and Ishida [31]. The interaction terms were then optimized to fit activity data [14, 21, 38] and the bT terms to fit solidus data [14, 15] and the phase boundaries evaluated by Barstad et al. [13]. The interaction term $^0L_{Ni,Va:Va:Te}^δ = ^0L_{Va:Ni,Va:Te}^δ$ affecting the range NiTe<sub>2</sub>–Te, was given a constant positive value in order to suppress a higher Te-solubility than about 66.67 at.% Te.

## β2, γ1 and γ2 phases

The end-members of the β2 phase were first manually adjusted to make the phase appear in approximately the correct temperature and composition range. Both end-members, i.e., β2 – Ni<sub>2</sub>Te and β2–NiTe, were modeled using the NKR. The stable composition range of the phase is narrow and lies roughly in the middle between the end-members; the end-members are thus rather far away from the compositions of available heat capacity data [29], which made it difficult to fit the data with end-members described as power series (Eq. 6). Instead, it worked well to model it as a power-series contribution to the interaction parameter $^0L_{Ni:Ni,Va:Te}^{β2}$. The a- and b-terms of the

interaction parameters together with the b-term of β2 – Ni<sub>2</sub>Te could then be optimized to fit activity data [35, 39] and the congruent reaction $β2 ⇒ β1$. The a-terms of the end-members were optimized to fit a sort of compromise between the enthalpy of formation of the end-members evaluated via DFT in this work and the calorimetric point by Shukla et al. [32]. If the latter enthalpy of formation were to be fit well, the phase would be too stable to allow the proper shape of the δ phase.

The a-term of the γ1 phase was first fitted to the DFT enthalpy of formation value evaluated in this work, but was then lowered to fit the phase diagram; with a formation enthalpy of the β2 phase so low, the γ1 phase enthalpy must also be rather low to form the Gibbs energy tangent between the β2 and δ necessary to have them all stable at low T. The b-term was fitted to relevant invariant reactions. The γ2 phase was added last, and its parameters were optimized to fit the phase diagram only, since there are virtually no other data available on the phase.

## DFT Results and discussion

The formation energies of end-members of Ni–Te phases and lattice stability of HCP_A3 tellurium are summarized in Table 5 together with converged parameters and lattice parameters compared with experimental values. For the γ1 phase calculations, full volume relaxation was not performed since it was desired to only affect the energy by removing atoms, without allowing the cell shape or atom movement to compensate. Out of the 50 different cells, the one with the largest total energy had the vacancies clustered close together, and one with the lowest energy had the vacancies spread over the supercell. Figure 2 shows how the average total energy of γ1 – Ni<sub>52</sub>Te<sub>40</sub> converges with the number of
<figure><img src="image_3.png" /><figcaption>**Figure 2** Convergence of average total energy of γ1 – Ni₅₂Te₄₀ with number of data samples out of the total of 50 different cells.</figcaption></figure>

sampled results. The last two points differ about 50 J/mol-atom.

## Calphad modeling results and discussion

### Thermodynamic properties

The formation enthalpy by Jandl et al. is very close to the calorimetric data point of a NiTe sample by Predel and Ruge [33], whose analyzed composition corresponds to a two-phase equilibrium of $δ + γ1$. Figure 3 shows the calculated enthalpy of formation, separately calculated, for all phases in the system. The enthalpy of the β2 phase closer fits DFT values than the experimental data point, while the $δ$ phase has a lower enthalpy than predicted by DFT at the NiTe end-member [34]. Therefore, the γ1 phase has a lower enthalpy than predicted by DFT in order to lie on a line in Gibbs energy between β2 and $δ$.

Figure 4 shows the calculated heat capacity of β2 compared with experimental data of 40 at.% Te and 41.1 at.% Te [28]. The lambda peaks were not modeled; it is a satisfying fit to the pseudo-linear portions of the datasets. The heat capacity of the $δ$ phase fits experimental data well [27, 29, 31], as shown in Figs. 5 and 6. Figure 6 shows for 57.1 at.% Te a bump in $C_p$ at around 350 K; this was not modeled but the general slope of the data is satisfactory.

<figure><img src="image_4.png" /><figcaption>**Figure 3** Calculated enthalpy of formation at 298 K for all phases of the Ni–Te system [32, 34].</figcaption></figure>

Figures 7 and 8 show the negative logarithm of the calculated thermodynamic activity of the Ni–Te system compared with the isopiestic data by Ettenberg [14] at 600 and 875 °C, respectively. It was difficult to obtain a good fit together with other activity data (as seen in Figs. 9, 10, 11 and 12), and the greater deviation from experiments in Fig. 8 is due to the Te-rich solidus of the $δ$ phase not matching the experimental phase diagram data [11, 15] at that temperature, as
<figure><img src="image_4.png" /><figcaption>**Figure 2** Convergence of average total energy of γ1 – Ni₅₂Te₄₀ with number of data samples out of the total of 50 different cells.</figcaption></figure>

sampled results. The last two points differ about 50 J/mol-atom.

## Calphad modeling results and discussion

### Thermodynamic properties

The formation enthalpy by Jandl et al. is very close to the calorimetric data point of a NiTe sample by Predel and Ruge [33], whose analyzed composition corresponds to a two-phase equilibrium of $δ + γ1$. Figure 3 shows the calculated enthalpy of formation, separately calculated, for all phases in the system. The enthalpy of the β2 phase closer fits DFT values than the experimental data point, while the $δ$ phase has a lower enthalpy than predicted by DFT at the NiTe end-member [34]. Therefore, the γ1 phase has a lower enthalpy than predicted by DFT in order to lie on a line in Gibbs energy between β2 and $δ$.

Figure 4 shows the calculated heat capacity of β2 compared with experimental data of 40 at.% Te and 41.1 at.% Te [28]. The lambda peaks were not modeled; it is a satisfying fit to the pseudo-linear portions of the datasets. The heat capacity of the $δ$ phase fits experimental data well [27, 29, 31], as shown in Figs. 5 and 6. Figure 6 shows for 57.1 at.% Te a bump in $C_p$ at around 350 K; this was not modeled but the general slope of the data is satisfactory.

<figure><img src="image_5.png" /><figcaption>**Figure 3** Calculated enthalpy of formation at 298 K for all phases of the Ni–Te system [32, 34].</figcaption></figure>

Figures 7 and 8 show the negative logarithm of the calculated thermodynamic activity of the Ni–Te system compared with the isopiestic data by Ettenberg [14] at 600 and 875 °C, respectively. It was difficult to obtain a good fit together with other activity data (as seen in Figs. 9, 10, 11 and 12), and the greater deviation from experiments in Fig. 8 is due to the Te-rich solidus of the $δ$ phase not matching the experimental phase diagram data [11, 15] at that temperature, as
<figure><img src="image_5.png" /><figcaption>**Figure 4** Calculated heat capacity of the β2 phase compared with experimental data [28].</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5** Calculated C<sub>P</sub> of the δ phase at low temperatures compared with experimental data [27].</figcaption></figure>

will be presented in "Phase diagram" section. Note that the chemical potential in Fig. 10 seems to deviate much from experimental data, but that the scale of the maximum deviation is only about 6%; it was difficult to improve the fit beyond this.

The description sufficiently reproduces the chemical potential data derived from EMF measurements by Carbonara and Hoch [21] and Geiderikh et al. [38], as seen in Figs. 11 and 12. The three upper points in the latter dataset (Fig. 12) lie in the γ2 + δ region, and the discrepancy is probably due to the γ2 description.

<figure><img src="image_7.png" /><figcaption>**Figure 6** Calculated C<sub>P</sub> of the δ phase at high temperatures compared with experimental data [29, 31].</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7** Negative logarithm of the tellurium activity of the Ni–Te system at 600 °C [14], with liquid as reference state.</figcaption></figure>

Phase diagram

The calculated phase diagram of the Ni–Te system is shown in Fig. 13. The final description is a compromise of fitting thermodynamic data and the phase diagram as best as reasonably achievable, and it is apparent from the figure that some aspects have been sacrificed. As seen in the zoomed-in Fig. 14 the β1 and β2, congruent and solvus lines do not fit perfectly, whereas the invariant reactions are well reproduced. The β1 liquidus on the Te-rich side is slightly higher than measured by Klepp and Komarek [15] while the solvus fits well.
<figure><img src="image_6.png" /><figcaption>**Figure 4** Calculated heat capacity of the β2 phase compared with experimental data [28].</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5** Calculated C<sub>P</sub> of the δ phase at low temperatures compared with experimental data [27].</figcaption></figure>

will be presented in "Phase diagram" section. Note that the chemical potential in Fig. 10 seems to deviate much from experimental data, but that the scale of the maximum deviation is only about 6%; it was difficult to improve the fit beyond this.

The description sufficiently reproduces the chemical potential data derived from EMF measurements by Carbonara and Hoch [21] and Geiderikh et al. [38], as seen in Figs. 11 and 12. The three upper points in the latter dataset (Fig. 12) lie in the γ2 + δ region, and the discrepancy is probably due to the γ2 description.

<figure><img src="image_8.png" /><figcaption>**Figure 6** Calculated C<sub>P</sub> of the δ phase at high temperatures compared with experimental data [29, 31].</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7** Negative logarithm of the tellurium activity of the Ni–Te system at 600 °C [14], with liquid as reference state.</figcaption></figure>

## Phase diagram

The calculated phase diagram of the Ni–Te system is shown in Fig. 13. The final description is a compromise of fitting thermodynamic data and the phase diagram as best as reasonably achievable, and it is apparent from the figure that some aspects have been sacrificed. As seen in the zoomed-in Fig. 14 the β1 and β2, congruent and solvus lines do not fit perfectly, whereas the invariant reactions are well reproduced. The β1 liquidus on the Te-rich side is slightly higher than measured by Klepp and Komarek [15] while the solvus fits well.
<figure><img src="image_9.png" /><figcaption>**Figure 8** Negative logarithm of the tellurium activity of the Ni–Te system at 875 °C [14], with liquid as reference state.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 9** Chemical potential of Ni and Te for the β1 and β2 phases, with temperature, compared with partial pressure and EMF data [35, 39].</figcaption></figure>

Furthermore, the δ phase was rather difficult to optimize. Prioritizing a decent fit of the solubility limits by Barstad et al. [13], the congruent point and δ + L liquidus together with activity data resulted in a rather large discrepancy in the Te-rich solidus at high temperature from the data by Westrum and Machol [11]. The δ–NiTe end-member is stable from 0 to 38 K (Fig. 13), an artifact of the enthalpies of formation that does not affect the application to nuclear reactors since the description is not valid below room temperature. The phase diagram largely resembles the experimental ones in the literature,

<figure><img src="image_11.png" /><figcaption>**Figure 10** Chemical potential of Te in the two-phase regions of FCC and respective β1 and β2 phases, with temperature, compared with partial pressure data [36].</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11** Chemical potential of Ni in δ phase at ca 400 °C compared with EMF data by *Carbonara* and Hoch [21].</figcaption></figure>

showing a wide solubility of δ at room temperature, as well as estimating that β2, γ1 are both stable there.

As a result of the optimization, a transition between NiAs-type disorder of interstitials and CdI₂ order appears in the δ phase. This transition was calculated for the metastable δ phase diagram and overlies the zoomed-in Ni–Te diagram in Fig. 15. The position of this line could be optimized to fit the data in the literature [18, 22], but doing so resulted in a poor fit of the phase boundaries; therefore, it was not prioritized with the application in mind. Figure 16 shows the site fractions of all constituents in the δ
<figure><img src="image_10.png" /><figcaption>**Figure 8** Negative logarithm of the tellurium activity of the Ni–Te system at 875 °C [14], with liquid as reference state.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9** Chemical potential of Ni and Te for the β1 and β2 phases, with temperature, compared with partial pressure and EMF data [35, 39].</figcaption></figure>

Furthermore, the δ phase was rather difficult to optimize. Prioritizing a decent fit of the solubility limits by Barstad et al. [13], the congruent point and δ + L liquidus together with activity data resulted in a rather large discrepancy in the Te-rich solidus at high temperature from the data by Westrum and Machol [11]. The δ–NiTe end-member is stable from 0 to 38 K (Fig. 13), an artifact of the enthalpies of formation that does not affect the application to nuclear reactors since the description is not valid below room temperature. The phase diagram largely resembles the experimental ones in the literature,

<figure><img src="image_12.png" /><figcaption>**Figure 10** Chemical potential of Te in the two-phase regions of FCC and respective β1 and β2 phases, with temperature, compared with partial pressure data [36].</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 11** Chemical potential of Ni in δ phase at ca 400 °C compared with EMF data by *Carbonara* and Hoch [21].</figcaption></figure>

showing a wide solubility of δ at room temperature, as well as estimating that β2, γ1 are both stable there.

As a result of the optimization, a transition between NiAs-type disorder of interstitials and $CdI_2$ order appears in the δ phase. This transition was calculated for the metastable δ phase diagram and overlies the zoomed-in Ni–Te diagram in Fig. 15. The position of this line could be optimized to fit the data in the literature [18, 22], but doing so resulted in a poor fit of the phase boundaries; therefore, it was not prioritized with the application in mind. Figure 16 shows the site fractions of all constituents in the δ
<figure><img src="image_13.png" /><figcaption>**Figure 12** Chemical potential of Ni in δ phase at 700 K compared with EMF data of Geiderikh et al. [38].</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 13** Calculated phase diagram of the Ni–Te system compared with experimental data [11, 13–15].</figcaption></figure>

phase from 50 to 66.67 at.% Te at 962 K. It is seen that the interstitials are fully disordered (NiAs-like) up to 52.06 at.% Te, i.e., outside the equilibrium solubility limit of Ni in the phase at 52.17 at.% Te, above which nickel rapidly shifts to the second sublattice; therefore, the completely disordered region only exists mainly in the metastable delta phase region, although there is a very narrow stable disordered interval right by the leftmost δ corner at the invariant $γ2 ↔ γ1 + δ$. At 54.3 at.% Te, the phase is very close to $CdI_2$ type ordering with $y''_{Ni} = 0.99$.

The experimental phase diagram seems to feature a metastable liquid miscibility gap evident by the kink

<figure><img src="image_15.png" /><figcaption>**Figure 14** Phase diagram in Fig. 13 zoomed-in on the high-temperature phase relations, compared with experimental data [11, 13–15].</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 15** Zoomed-in region of Ni–Te phase diagram showing the order–disorder transition in the δ phase. Dashed line described by the condition $y'_{Ni} - y''_{Ni} = 0.001$.</figcaption></figure>

in curvature of the Ni-rich liquidus (see Fig 1). Miscibility gaps in transition-metal telluride liquids are not surprising. It is acceptable that the liquidus in this region is about 100 K lower than the “melting effect” noted in the literature [15], since this data point might not be very accurate at such high temperatures with the experimental method used (above the stability limit of silica capsules). It is, however, possible that there is a stable miscibility gap in the liquid and that the reported melting effect [15] is the invariant
<figure><img src="image_14.png" /><figcaption>**Figure 12** Chemical potential of Ni in δ phase at 700 K compared with EMF data of Geiderikh et al. [38].</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 13** Calculated phase diagram of the Ni–Te system compared with experimental data [11, 13–15].</figcaption></figure>

phase from 50 to 66.67 at.% Te at 962 K. It is seen that the interstitials are fully disordered (NiAs-like) up to 52.06 at.% Te, i.e., outside the equilibrium solubility limit of Ni in the phase at 52.17 at.% Te, above which nickel rapidly shifts to the second sublattice; therefore, the completely disordered region only exists mainly in the metastable delta phase region, although there is a very narrow stable disordered interval right by the leftmost δ corner at the invariant $γ2 ↔ γ1 + δ$. At 54.3 at.% Te, the phase is very close to $CdI_2$ type ordering with $y''_{Ni} = 0.99$.

The experimental phase diagram seems to feature a metastable liquid miscibility gap evident by the kink

<figure><img src="image_16.png" /><figcaption>**Figure 14** Phase diagram in Fig. 13 zoomed-in on the high-temperature phase relations, compared with experimental data [11, 13–15].</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 15** Zoomed-in region of Ni–Te phase diagram showing the order–disorder transition in the δ phase. Dashed line described by the condition $y'_{Ni} - y''_{Ni} = 0.001$.</figcaption></figure>

in curvature of the Ni-rich liquidus (see Fig 1). Miscibility gaps in transition-metal telluride liquids are not surprising. It is acceptable that the liquidus in this region is about 100 K lower than the “melting effect” noted in the literature [15], since this data point might not be very accurate at such high temperatures with the experimental method used (above the stability limit of silica capsules). It is, however, possible that there is a stable miscibility gap in the liquid and that the reported melting effect [15] is the invariant
<figure><img src="image_17.png" /><figcaption>**Figure 16** Variation of site fractions of all constituents with composition in the δ phase, at 962 K. $y'''_{Te} = 1$ and is omitted.</figcaption></figure>

Liquid#1 ↔ Liquid#2 + FCC reaction, which would make it similar to our assessment of the Fe–Te system [8]. Until whether there are miscibility gaps in these systems or not is experimentally proven, it is merely a matter of taste and convenience if one chooses to model them. The Ni–Te liquidus can be described rather well without even a metastable miscibility gap in the liquid, while the Fe–Te liquidus cannot.

## Conclusions and future work

A thermodynamic assessment of the Ni–Te system has been performed, supported by DFT calculations. The thermodynamic description is deemed good for the conditions present in the application to nuclear reactors. The liquid was modeled using the ionic two-sublattice model. There is an order–disorder transition in the metastable δ phase but its exact position was not optimized. The description does not model all known ordered superstructures of the Ni–Te alloys; for such applications where that is desired, the

description should be modified to include those phases, e.g., by modifying the sublattice models accordingly, introducing ordering parameters, or modeling them as separate phases. An assessment of the Fe–Ni–Te system is in progress.

## Acknowledgements

The authors are grateful to the Swedish Research Council (Vetenskapsrådet) for funding the SAFARI project. Carl-Magnus thanks his colleagues at CEA and KTH, as well as Drs Bonnie Lindahl and Nathalie Dupin for invaluable advice on modeling, Andrei Ruban for advice about DFT calculations. C.-M. further extends thanks to Mats Kronberg at the National Supercomputer Centre, Linköping, for keeping a backup of the DFT data when it was believed to be erased.

## Compliance with ethical standards

**Conflicts of interest** The authors declare that they have no conflict of interest.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Appendix 1: thermodynamic tables

See Tables 6 and 7.
<figure><img src="image_18.png" /><figcaption>**Figure 16** Variation of site fractions of all constituents with composition in the δ phase, at 962 K. $y'''_{Te} = 1$ and is omitted.</figcaption></figure>

Liquid#1 ↔ Liquid#2 + FCC reaction, which would make it similar to our assessment of the Fe–Te system [8]. Until whether there are miscibility gaps in these systems or not is experimentally proven, it is merely a matter of taste and convenience if one chooses to model them. The Ni–Te liquidus can be described rather well without even a metastable miscibility gap in the liquid, while the Fe–Te liquidus cannot.

## Conclusions and future work

A thermodynamic assessment of the Ni–Te system has been performed, supported by DFT calculations. The thermodynamic description is deemed good for the conditions present in the application to nuclear reactors. The liquid was modeled using the ionic two-sublattice model. There is an order–disorder transition in the metastable δ phase but its exact position was not optimized. The description does not model all known ordered superstructures of the Ni–Te alloys; for such applications where that is desired, the

description should be modified to include those phases, e.g., by modifying the sublattice models accordingly, introducing ordering parameters, or modeling them as separate phases. An assessment of the Fe–Ni–Te system is in progress.

## Acknowledgements

The authors are grateful to the Swedish Research Council (Vetenskapsrådet) for funding the SAFARI project. Carl-Magnus thanks his colleagues at CEA and KTH, as well as Drs Bonnie Lindahl and Nathalie Dupin for invaluable advice on modeling, Andrei Ruban for advice about DFT calculations. C.-M. further extends thanks to Mats Kronberg at the National Supercomputer Centre, Linköping, for keeping a backup of the DFT data when it was believed to be erased.

## Compliance with ethical standards

**Conflicts of interest** The authors declare that they have no conflict of interest.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Appendix 1: thermodynamic tables

See Tables 6 and 7.
<table><caption>**Table 6** Thermodynamic functions of the Ni–Te system produced by this work. *: 3 times lattice stability of FCC-Te calculated in an assessment of Fe–Te [8]</caption><thead><tr><th colspan="2">Liquid: (Ni<sup>+2</sup>)P(Te<sup>-2</sup>, Va<sup>-Q</sup>, Te<sup>0</sup>)Q</th></tr></thead><tbody><tr><td>°G<sub>Ni<sup>+2</sup>,Va<sup>-Q</sup></sub> = GLIQNI</td><td>0L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Va<sup>-Q</sup></sub> = -64653 + 50.504 · T</td></tr><tr><td>°G<sub>Te<sup>0</sup></sub> = GLIQTE</td><td>0L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Te<sup>0</sup></sub> = -13203</td></tr><tr><td>°G<sub>Ni<sup>+2</sup>,Te<sup>-2</sup></sub> = -96269 + 4.125 · T + 2 · GLIQNI + 2 · GLIQTE</td><td>1L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Va<sup>-Q</sup></sub> = -16.231 · T, 1L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Te<sup>0</sup></sub> = -2.918 · T</td></tr><tr><td>β1: (Ni, Te)<sub>2</sub>(Te)<sub>1</sub></td><td>Wyckoff: Model not based on sites</td></tr><tr><td>°G<sub>Ni<sub>2</sub>Te</sub> = -34833 - 22.432 · T + 2 · GHSERNI + GHSERTE</td><td>0L<sub>Ni,Te:Te</sub> = -230574</td></tr><tr><td>°G<sub>Te<sub>2</sub>Te</sub> = 3 · GHSERTE + 107100*</td><td>1L<sub>Ni,Te:Te</sub> = -55942</td></tr><tr><td>β2: (Ni)<sub>1</sub>(Te)<sub>1</sub>(Ni, Va)<sub>1</sub></td><td>Wyckoff: (1a)(1b)(2d)</td></tr><tr><td>°G<sub>NiTeNi</sub> = -25019 - 25.706 · T + 2 · GHSERNI + GHSERTE</td><td>°G<sub>NiTeVa</sub> = -26042 - 10.813 · T + GHSERNI + GHSERTE</td></tr><tr><td>0L<sub>Ni:Te:Ni,Va</sub> = -121972 + 7.208 · T + 16.1631 · T ln T - 0.037774 · T<sup>2</sup></td><td>1L<sub>Ni:Te:Ni,Va</sub> = 19.998 · T</td></tr><tr><td>δ: (Ni, Va)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>2</sub></td><td>Wyckoff: (1a)(1b)(2d)</td></tr><tr><td>°G<sub>NiNiTe<sub>2</sub></sub> = -120222 + 465.516 · T - 93.3232 · T ln T - 0.019815 · T<sup>2</sup> + 106000 · T<sup>-1</sup></td><td></td></tr><tr><td>°G<sub>NiVaTe<sub>2</sub></sub> = -92171 + 350.832 · T - 71.0051 · T ln T - 0.011639 · T<sup>2</sup> + 95000 · T<sup>-1</sup> + (7.1 × 10<sup>-8</sup>) · T<sup>3</sup></td><td></td></tr><tr><td>°G<sub>VaNiTe<sub>2</sub></sub> = -92171 + 350.832 · T - 71.0051 · T ln T - 0.011639 · T<sup>2</sup> + 95000 · T<sup>-1</sup> + (7.1 × 10<sup>-8</sup>) · T<sup>3</sup></td><td></td></tr><tr><td>°G<sub>VaVaTe<sub>2</sub></sub> = 13547 + 2 · GHSERTE</td><td>0L<sub>Ni,Va:Va:Te</sub> = 0 L<sub>Va:Ni,Va:Te</sub> = 20000</td></tr><tr><td>0L<sub>Ni,Va:Ni:Te</sub> = 0 L<sub>Ni:Ni,Va:Te</sub> = 3897 - 9.541 · T</td><td>1L<sub>Ni,Va:Ni:Te</sub> = 1 L<sub>Ni:Ni,Va:Te</sub> = 5 · T</td></tr><tr><td>γ1: (Ni)<sub>52</sub>(Te)<sub>40</sub></td><td>Wyckoff: (4c)(4c)</td></tr><tr><td>°G<sub>Ni<sub>2</sub>Te<sub>40</sub></sub> = -2175569 - 128.644 · T + 52 · GHSERNI + 40 · GHSERTE</td><td></td></tr><tr><td>γ2: (Ni)<sub>20</sub>(Te)<sub>17</sub></td><td>Wyckoff: N/A</td></tr><tr><td>°G<sub>Ni<sub>20</sub>Te<sub>17</sub></sub> = -747938 - 202.534 · T + 20 · GHSERNI + 17 · GHSERTE</td><td></td></tr><tr><td>HCP_A3: (Ni, Te)<sub>1</sub>(Va)<sub>0.5</sub></td><td>Wyckoff: (2c)(2a)</td></tr><tr><td>°G<sub>TeVa<sub>0.5</sub></sub> = 17300 + GHSERTE</td><td></td></tr><tr><td>FCC_A1: (Ni, Te)<sub>1</sub>(Va)<sub>1</sub></td><td>Wyckoff: (4a)(4b)</td></tr><tr><td>°G<sub>TeVa</sub> = 35700 + GHSERTE [8]</td><td></td></tr></tbody></table>

<table><caption>**Table 7** Invariant equilibria above 300 K in the Ni–Te system, calculated from the thermodynamic description in this work compared with the assessment of Lee and Nash [16]</caption><thead><tr><th>Reaction</th><th>at.% Te phase 1</th><th>at.% Te phase 2</th><th>at.% Te phase 3</th><th>T [K]</th><th>Ref. T [K]</th><th>Reaction type</th></tr></thead><tbody><tr><td>Liq ↔ β1</td><td>38.14</td><td>38.14</td><td></td><td>1295.17</td><td>1294.7</td><td>Congruent</td></tr><tr><td>Liq ↔ FCC + β1</td><td>32.77</td><td>0.007</td><td>37.10</td><td>1271.0</td><td>1277.7</td><td>Eutectic</td></tr><tr><td>β1 ↔ β2</td><td>40.77</td><td>40.77</td><td></td><td>1057.8</td><td>1063.2</td><td>Congruent</td></tr><tr><td>β1 ↔ FCC + β2</td><td>37.44</td><td>0.001</td><td>39.45</td><td>1004.6</td><td>1004.7</td><td>Eutectoid</td></tr><tr><td>β1 + Liq ↔ γ2</td><td>43.37</td><td>48.87</td><td>45.95</td><td>1153.4</td><td>1153.2</td><td>Peritectic</td></tr><tr><td>β1 ↔ β2 + γ2</td><td>41.82</td><td>41.15</td><td>45.95</td><td>1052.8</td><td>1048.2</td><td>Eutectoid</td></tr><tr><td>β2 + γ2 ↔ γ1</td><td>41.20</td><td>45.95</td><td>43.48</td><td>1012.0</td><td>1015.7</td><td>Peritectoid</td></tr><tr><td>Liq ↔ γ2 + δ</td><td>49.50</td><td>45.95</td><td>52.86</td><td>1147.9</td><td>1146.2</td><td>Eutectic</td></tr><tr><td>γ2 ↔ γ1 + δ</td><td>45.95</td><td>43.48</td><td>52.03</td><td>962.2</td><td>963.2</td><td>Eutectoid</td></tr><tr><td>Liq ↔ δ</td><td>56.87</td><td>56.87</td><td></td><td>1172.8</td><td>1173.7</td><td>Congruent</td></tr><tr><td>Liq ↔ δ + Te - A8</td><td>99.32</td><td>66.71</td><td>100</td><td>720.9</td><td>721.7</td><td>Eutectic</td></tr></tbody></table>
<table><caption>**Table 6** Thermodynamic functions of the Ni–Te system produced by this work. *: 3 times lattice stability of FCC-Te calculated in an assessment of Fe–Te [8]</caption><thead><tr><th colspan="2">Liquid: (Ni<sup>+2</sup>)P(Te<sup>-2</sup>, Va<sup>-Q</sup>, Te<sup>0</sup>)Q</th></tr></thead><tbody><tr><td>°G<sub>Ni<sup>+2</sup>,Va<sup>-Q</sup></sub> = GLIQNI</td><td>0L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Va<sup>-Q</sup></sub> = -64653 + 50.504 · T</td></tr><tr><td>°G<sub>Te<sup>0</sup></sub> = GLIQTE</td><td>0L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Te<sup>0</sup></sub> = -13203</td></tr><tr><td>°G<sub>Ni<sup>+2</sup>,Te<sup>-2</sup></sub> = -96269 + 4.125 · T + 2 · GLIQNI + 2 · GLIQTE</td><td>1L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Va<sup>-Q</sup></sub> = -16.231 · T, 1L<sub>Ni<sup>+2</sup>,Te<sup>-2</sup>,Te<sup>0</sup></sub> = -2.918 · T</td></tr><tr><td>β1: (Ni, Te)<sub>2</sub>(Te)<sub>1</sub></td><td>Wyckoff: Model not based on sites</td></tr><tr><td>°G<sub>Ni<sub>2</sub>Te</sub> = -34833 - 22.432 · T + 2 · GHSERNI + GHSERTE</td><td>0L<sub>Ni,Te:Te</sub> = -230574</td></tr><tr><td>°G<sub>Te<sub>2</sub>Te</sub> = 3 · GHSERTE + 107100*</td><td>1L<sub>Ni,Te:Te</sub> = -55942</td></tr><tr><td>β2: (Ni)<sub>1</sub>(Te)<sub>1</sub>(Ni, Va)<sub>1</sub></td><td>Wyckoff: (1a)(1b)(2d)</td></tr><tr><td>°G<sub>NiTeNi</sub> = -25019 - 25.706 · T + 2 · GHSERNI + GHSERTE</td><td>°G<sub>NiTeVa</sub> = -26042 - 10.813 · T + GHSERNI + GHSERTE</td></tr><tr><td>0L<sub>Ni:Te:Ni,Va</sub> = -121972 + 7.208 · T + 16.1631 · T ln T - 0.037774 · T<sup>2</sup></td><td>1L<sub>Ni:Te:Ni,Va</sub> = 19.998 · T</td></tr><tr><td>δ: (Ni, Va)<sub>1</sub>(Ni, Va)<sub>1</sub>(Te)<sub>2</sub></td><td>Wyckoff: (1a)(1b)(2d)</td></tr><tr><td>°G<sub>NiNiTe<sub>2</sub></sub> = -120222 + 465.516 · T - 93.3232 · T ln T - 0.019815 · T<sup>2</sup> + 106000 · T<sup>-1</sup></td><td></td></tr><tr><td>°G<sub>NiVaTe<sub>2</sub></sub> = -92171 + 350.832 · T - 71.0051 · T ln T - 0.011639 · T<sup>2</sup> + 95000 · T<sup>-1</sup> + (7.1 × 10<sup>-8</sup>) · T<sup>3</sup></td><td></td></tr><tr><td>°G<sub>VaNiTe<sub>2</sub></sub> = -92171 + 350.832 · T - 71.0051 · T ln T - 0.011639 · T<sup>2</sup> + 95000 · T<sup>-1</sup> + (7.1 × 10<sup>-8</sup>) · T<sup>3</sup></td><td></td></tr><tr><td>°G<sub>VaVaTe<sub>2</sub></sub> = 13547 + 2 · GHSERTE</td><td>0L<sub>Ni,Va:Va:Te</sub> = 0 L<sub>Va:Ni,Va:Te</sub> = 20000</td></tr><tr><td>0L<sub>Ni,Va:Ni:Te</sub> = 0 L<sub>Ni:Ni,Va:Te</sub> = 3897 - 9.541 · T</td><td>1L<sub>Ni,Va:Ni:Te</sub> = 1 L<sub>Ni:Ni,Va:Te</sub> = 5 · T</td></tr><tr><td>γ1: (Ni)<sub>52</sub>(Te)<sub>40</sub></td><td>Wyckoff: (4c)(4c)</td></tr><tr><td>°G<sub>Ni<sub>2</sub>Te<sub>40</sub></sub> = -2175569 - 128.644 · T + 52 · GHSERNI + 40 · GHSERTE</td><td></td></tr><tr><td>γ2: (Ni)<sub>20</sub>(Te)<sub>17</sub></td><td>Wyckoff: N/A</td></tr><tr><td>°G<sub>Ni<sub>20</sub>Te<sub>17</sub></sub> = -747938 - 202.534 · T + 20 · GHSERNI + 17 · GHSERTE</td><td></td></tr><tr><td>HCP_A3: (Ni, Te)<sub>1</sub>(Va)<sub>0.5</sub></td><td>Wyckoff: (2c)(2a)</td></tr><tr><td>°G<sub>TeVa<sub>0.5</sub></sub> = 17300 + GHSERTE</td><td></td></tr><tr><td>FCC_A1: (Ni, Te)<sub>1</sub>(Va)<sub>1</sub></td><td>Wyckoff: (4a)(4b)</td></tr><tr><td>°G<sub>TeVa</sub> = 35700 + GHSERTE [8]</td><td></td></tr></tbody></table>

<table><caption>**Table 7** Invariant equilibria above 300 K in the Ni–Te system, calculated from the thermodynamic description in this work compared with the assessment of Lee and Nash [16]</caption><thead><tr><th>Reaction</th><th>at.% Te phase 1</th><th>at.% Te phase 2</th><th>at.% Te phase 3</th><th>T [K]</th><th>Ref. T [K]</th><th>Reaction type</th></tr></thead><tbody><tr><td>Liq ↔ β1</td><td>38.14</td><td>38.14</td><td></td><td>1295.17</td><td>1294.7</td><td>Congruent</td></tr><tr><td>Liq ↔ FCC + β1</td><td>32.77</td><td>0.007</td><td>37.10</td><td>1271.0</td><td>1277.7</td><td>Eutectic</td></tr><tr><td>β1 ↔ β2</td><td>40.77</td><td>40.77</td><td></td><td>1057.8</td><td>1063.2</td><td>Congruent</td></tr><tr><td>β1 ↔ FCC + β2</td><td>37.44</td><td>0.001</td><td>39.45</td><td>1004.6</td><td>1004.7</td><td>Eutectoid</td></tr><tr><td>β1 + Liq ↔ γ2</td><td>43.37</td><td>48.87</td><td>45.95</td><td>1153.4</td><td>1153.2</td><td>Peritectic</td></tr><tr><td>β1 ↔ β2 + γ2</td><td>41.82</td><td>41.15</td><td>45.95</td><td>1052.8</td><td>1048.2</td><td>Eutectoid</td></tr><tr><td>β2 + γ2 ↔ γ1</td><td>41.20</td><td>45.95</td><td>43.48</td><td>1012.0</td><td>1015.7</td><td>Peritectoid</td></tr><tr><td>Liq ↔ γ2 + δ</td><td>49.50</td><td>45.95</td><td>52.86</td><td>1147.9</td><td>1146.2</td><td>Eutectic</td></tr><tr><td>γ2 ↔ γ1 + δ</td><td>45.95</td><td>43.48</td><td>52.03</td><td>962.2</td><td>963.2</td><td>Eutectoid</td></tr><tr><td>Liq ↔ δ</td><td>56.87</td><td>56.87</td><td></td><td>1172.8</td><td>1173.7</td><td>Congruent</td></tr><tr><td>Liq ↔ δ + Te - A8</td><td>99.32</td><td>66.71</td><td>100</td><td>720.9</td><td>721.7</td><td>Eutectic</td></tr></tbody></table>
References

[1] Adamson MG, Aitken EA, Lindemer TB (1985) Chemical thermodynamics of Cs and Te fission product interactions in irradiated LMFBR mixed-oxide fuel pins. J Nuclear Mater 130:375–392. https://doi.org/10.1016/0022-3115(85)90325-3

[2] Adamson MG, Aitken EA (1985) Cs and Te fission product-induced attack and embrittlement of stainless steel cladding in oxide fuel pins. J Nuclear Mater 132:160–166

[3] Pulham RJ, Richards MW (1990a) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys Part I—the corrosion by tellurium. J Nuclear Mater 171:319–326. https://doi.org/10.1016/0022-3115(90)90378-Z

[4] Pulham RJ, Richards MW (1990b) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part II—the corrosion by caesium-oxygen mixtures and oxygen. J Nuclear Mater 172:47–53. https://doi.org/10.1016/0022-3115(90)90008-B

[5] Pulham RJ, Richards MW (1990) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part III—the effect of oxygen potential on the corrosion by caesium-tellurium mixtures. J Nuclear Mater 172:206–219. https://doi.org/10.1016/0022-3115(90)90439-T

[6] Pulham RJ, Richards MW (1990d) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part IV—the corrosion of ferritic steels. J Nuclear Mater 172:304–313. https://doi.org/10.1016/0022-3115(90)90285-U

[7] OECD-NEA (2017) NEA Nuclear Science Committee: thermodynamics of advanced fuels—International Database (TAF-ID). https://www.oecd-nea.org/science/taf-id/. Accessed 03 Dec 2018

[8] Arvhult CM, Guéneau C, Gossé S, Selleby M (2018) Thermodynamic assessment of the Fe–Te system. Part II: thermodynamic modelling. J Alloys Compd 767:883–893. https://doi.org/10.1016/j.jallcom.2018.07.051

[9] Cacciamani G, Dinsdale A, Palumbo M, Pasturel A (2010) The Fe–Ni system: thermodynamic modelling assisted by atomistic calculations. Intermetallics 18(6):1148–1162. http://doi.org/10.1016/j.intermet.2010.02.026

[10] Uchida E, Kondoh H (1956) Magnetic properties of nickel telluride. J Phys Soc Jpn 11(1):21–27. https://doi.org/10.1143/JPSJ.11.21

[11] Westrum EF, Machol RE (1958) Thermodynamics of non-stoichiometric nickel tellurides. II. Dissociation pressures and phase relations of tellurium-rich compositions. J Chem Phys 29(4):824–828. https://doi.org/10.1063/1.1744597

[12] Kok RB, Wiegers GA, Jellinek F (1965) The system nickel-tellurium. 1. Structure and some superstructures of the Ni3+-qTe2 phase. Recueil 84:1585–1588

[13] Barstad J, Grønvold F, Røst E, Vestersjø E (1966) On the tellurides of Nickel. Acta Chem Scand 20:2865–2879

[14] Ettenberg M, Komarek KL, Miller E (1970) Thermodynamic properties of Nickel–Tellurium alloys. J Solid State Chem 1:583–592

[15] Klepp KO, Komarek KL (1972) Übergangsmetall-Chalkogensysteme, 3. Mitt.: das system nickel-tellur. Monatshefte fur Chemie - Chemical Monthly Chem 103:934–946

[16] Lee SY, Nash P (1990) Ni–Te (Nickel–Tellurium). In: Massalski TB, Okamoto H, Subramanian PR, Kacprzak L (eds) Binary alloy phase diagrams, vol 3, 2nd edn. ASM International, pp 2869–2872

[17] Gulay LD, Olekseyuk ID (2004) Crystal structures of the compounds Ni3Te2, Ni3-δTe2(δ = 0.12) and Ni1.29Te. J Alloys Compd 376(1–2):131–138. https://doi.org/10.1016/j.jallcom.2003.12.022

[18] Norén L, Ting V, Withers RL, van Tendeloo G (2001) An electron and X-Ray diffraction investigation of Ni1+xTe2 and Ni1+xSe2 CdI2/NiAs type solid solution phases. J Solid State Chem 161:266–273. https://doi.org/10.1006/jssc.2001.9309

[19] Gulay LD, Daszkiewicz M, Pietraszko A (2007) Evidence of a centre of symmetry: redetermination of Ni 2.60Te2 from single-crystal data. Acta Crystallogr Sect E Struct Rep Online 63(11):i188. https://doi.org/10.1107/S1600536807050568

[20] Stevels ALN (1969) Phase transitions in Nickel and Copper Selenides and Tellurides. Ph.D. thesis, University of Groningen

[21] Carbonara RS, Hoch M (1972) Thermodynamics and Structure of Ni1-xTe. Monatshefte fur Chemie Chemical Monthly Chem 715(103):695–715

[22] Coffin P, Jacobson AJ, Fender BEF (1974) The observation of an order-disorder transition in the Ni–Te system by neutron diffraction. J Phys C Solid State Phys 7:2781–2790

[23] Bensch W, Heid W (1996) Anionic polymeric bonds in nickel ditelluride: crystal structure, and experimental and theoretical band structure. J Solid State Chem 121(13):87–94. https://doi.org/10.1006/jssc.1996.0013

[24] Abakarov SA, Bagduev GB, Dazhaev PS (1970) Solubility of Various Elements in Tellurium. Izvestiya Akademii Nauk SSSR, Neorganicheskie Materialy 6(6):1169–1170

[25] Yamamoto K, Kashida S (1991) X-ray study of the average structures of Cu2Se and Cu1.8S in the room temperature and the high temperature phases. J Solid State Chem. https://doi.org/10.1016/0022-4596(91)90289-T
References

[1] Adamson MG, Aitken EA, Lindemer TB (1985) Chemical thermodynamics of Cs and Te fission product interactions in irradiated LMFBR mixed-oxide fuel pins. J Nuclear Mater 130:375–392. https://doi.org/10.1016/0022-3115(85)90325-3

[2] Adamson MG, Aitken EA (1985) Cs and Te fission product-induced attack and embrittlement of stainless steel cladding in oxide fuel pins. J Nuclear Mater 132:160–166

[3] Pulham RJ, Richards MW (1990a) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys Part I—the corrosion by tellurium. J Nuclear Mater 171:319–326. https://doi.org/10.1016/0022-3115(90)90378-Z

[4] Pulham RJ, Richards MW (1990b) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part II—the corrosion by caesium-oxygen mixtures and oxygen. J Nuclear Mater 172:47–53. https://doi.org/10.1016/0022-3115(90)90008-B

[5] Pulham RJ, Richards MW (1990) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part III—the effect of oxygen potential on the corrosion by caesium-tellurium mixtures. J Nuclear Mater 172:206–219. https://doi.org/10.1016/0022-3115(90)90439-T

[6] Pulham RJ, Richards MW (1990d) Chemical reactions of caesium, tellurium and oxygen with fast breeder reactor cladding alloys. Part IV—the corrosion of ferritic steels. J Nuclear Mater 172:304–313. https://doi.org/10.1016/0022-3115(90)90285-U

[7] OECD-NEA (2017) NEA Nuclear Science Committee: thermodynamics of advanced fuels—International Database (TAF-ID). https://www.oecd-nea.org/science/taf-id/. Accessed 03 Dec 2018

[8] Arvhult CM, Guéneau C, Gossé S, Selleby M (2018) Thermodynamic assessment of the Fe–Te system. Part II: thermodynamic modelling. J Alloys Compd 767:883–893. https://doi.org/10.1016/j.jallcom.2018.07.051

[9] Cacciamani G, Dinsdale A, Palumbo M, Pasturel A (2010) The Fe–Ni system: thermodynamic modelling assisted by atomistic calculations. Intermetallics 18(6):1148–1162. http://doi.org/10.1016/j.intermet.2010.02.026

[10] Uchida E, Kondoh H (1956) Magnetic properties of nickel telluride. J Phys Soc Jpn 11(1):21–27. https://doi.org/10.1143/JPSJ.11.21

[11] Westrum EF, Machol RE (1958) Thermodynamics of non-stoichiometric nickel tellurides. II. Dissociation pressures and phase relations of tellurium-rich compositions. J Chem Phys 29(4):824–828. https://doi.org/10.1063/1.1744597

[12] Kok RB, Wiegers GA, Jellinek F (1965) The system nickel-tellurium. 1. Structure and some superstructures of the Ni3+-qTe2 phase. Recueil 84:1585–1588

[13] Barstad J, Grønvold F, Røst E, Vestersjø E (1966) On the tellurides of Nickel. Acta Chem Scand 20:2865–2879

[14] Ettenberg M, Komarek KL, Miller E (1970) Thermodynamic properties of Nickel–Tellurium alloys. J Solid State Chem 1:583–592

[15] Klepp KO, Komarek KL (1972) Übergangsmetall-Chalkogensysteme, 3. Mitt.: das system nickel-tellur. Monatshefte für Chemie - Chemical Monthly Chem 103:934–946

[16] Lee SY, Nash P (1990) Ni–Te (Nickel–Tellurium). In: Massalski TB, Okamoto H, Subramanian PR, Kacprzak L (eds) Binary alloy phase diagrams, vol 3, 2nd edn. ASM International, pp 2869–2872

[17] Gulay LD, Olekseyuk ID (2004) Crystal structures of the compounds Ni3Te2, Ni3-δTe2(δ = 0.12) and Ni1.29Te. J Alloys Compd 376(1–2):131–138. https://doi.org/10.1016/j.jallcom.2003.12.022

[18] Norén L, Ting V, Withers RL, van Tendeloo G (2001) An electron and X-Ray diffraction investigation of Ni1+xTe2 and Ni1+xSe2 CdI2/NiAs type solid solution phases. J Solid State Chem 161:266–273. https://doi.org/10.1006/jssc.2001.9309

[19] Gulay LD, Daszkiewicz M, Pietraszko A (2007) Evidence of a centre of symmetry: redetermination of Ni 2.60Te2 from single-crystal data. Acta Crystallogr Sect E Struct Rep Online 63(11):i188. https://doi.org/10.1107/S1600536807050568

[20] Stevels ALN (1969) Phase transitions in Nickel and Copper Selenides and Tellurides. Ph.D. thesis, University of Groningen

[21] Carbonara RS, Hoch M (1972) Thermodynamics and Structure of Ni1-xTe. Monatshefte für Chemie Chemical Monthly Chem 715(103):695–715

[22] Coffin P, Jacobson AJ, Fender BEF (1974) The observation of an order-disorder transition in the Ni–Te system by neutron diffraction. J Phys C Solid State Phys 7:2781–2790

[23] Bensch W, Heid W (1996) Anionic polymeric bonds in nickel ditelluride: crystal structure, and experimental and theoretical band structure. J Solid State Chem 121(13):87–94. https://doi.org/10.1006/jssc.1996.0013

[24] Abakarov SA, Bagduev GB, Dazhaev PS (1970) Solubility of Various Elements in Tellurium. Izvestiya Akademii Nauk SSSR, Neorganicheskie Materialy 6(6):1169–1170

[25] Yamamoto K, Kashida S (1991) X-ray study of the average structures of Cu2Se and Cu1.8S in the room temperature and the high temperature phases. J Solid State Chem. https://doi.org/10.1016/0022-4596(91)90289-T
[26] Adenis C, Langer V, Lindqvist O (1989) Reinvestigation of the structure of tellurium. Acta Crystallogr Sect C Cryst Struct Commun 45(6):941–942. https://doi.org/10.1107/S0108270188014453

[27] Westrum EF, Chou JC, Machol RE (1958) Thermodynamics of nonstoichiometric Nickel tellurides. 1. Heat Capacity and Thermodynamic Functions of the delta phase from 5 to 350 K. J Chem Phys 28(3):497–503

[28] Grønvold F, Kveseth NJ, Sveen A (1972) The Ni3+-xTe2 phases. Heat capacities in the range 298 to 900 k and transition behavior. J Chem Thermodyn 4:337–348

[29] Grønvold F (1973) Heat capacities and thermodynamic properties of Cr3Te4and Ni3Te4from 298 to 950 K. Structural order-disorder transitions. J Chem Thermodyn 5(4):545–551. https://doi.org/10.1016/S0021-9614(73)80102-8

[30] Mills KC (1974) Heat capacity of nickel and cobalt tellurides. J Chem Soc Faraday Trans 1 Phys Chem Condens Phases 70(7):2224–2231. https://doi.org/10.1039/f19747002224

[31] Tsuji T, Ishida K (1995) heat capacities of Cr5Te8 and Ni5Te8 from 325 to 930 K. Thermochim Acta 253:11–18

[32] Shukla NK, Prasad R, Roy KN, Sood DD (1990) Standard molar enthalpies of formation at the temperature 298.15 K of iron telluride (FeTe0.9) and of nickel telluride (Ni0.595Te0.405). J Chem Thermodyn 22(9):899–903. http://doi.org/10.1016/0021-9614(90)90178-S

[33] Predel B, Ruge H (1972) Bildungsenthalpien und bindungsverhältnisse in einigen intermetallischen verbindungen vom NiAs-Typ. Thermochim Acta 3(5):411–418. https://doi.org/10.1016/0040-6031(72)87055-2

[34] Jandl I, Ipser H, Richter KW (2015) Thermodynamic modelling of the general NiAs-type structure: a study of first principle energies of formation for binary Ni-containing B8 compounds. Calphad Comput Coupling Phase Diagr Thermochem 50:174–181. https://doi.org/10.1016/j.calphad.2015.06.006

[35] Prasad R, IYer VS, Venugopal V, Sundaresh V, Singh Z, Sood DD (1987) Partial and integral molar thermodynamic properties of NixTe1-x(s, x=0.595 to 0.630) alloys. J Chem Thermodyn 19:891–896

[36] Viswanathan R, Baba MS, Raj DDA, Balasubramanian R, Saha B, Mathews CK (1987) Vaporisation thermodynamics of the nickel-rich phases in the Ni–Te binary system–a high-temperature mass-spectrometric study. J Nuclear Mater 149:302–311

[37] Kutsenok IB, Geiderikh VA, Valuev IA (1980) Use of the instantaneous measuring of EMF method for the study of

thermodynamic properties of solid alloys. Vestnik Moskovskogo Universiteta Seriya 2 Khimiya 21(6):554–558

[38] Geiderikh VA, Sheveleva SN, Kutsenok IB, Krivovsheya NS (1980) Thermodynamic properties of nickel tellurides. Zhurnal Fizicheskoi Khimii 54(4):1068–1071

[39] Prasad R, Iyer VS, Venugopal V, Singh Z, Sood DD (1987) Thermodynamic Study of Ni3Te2 by E.M.F. measurements. J Chem Thermodyn 19:613–616

[40] Maekawa T, Yokokawa T (1975) The enthalpies of binary mixtures of chalcogen elements with various metals. J Chem Thermodyn 7:505–506

[41] Kresse G, Hafner J (1993) Ab initio molecular dynamics for liquid metals. Phys Rev B. https://doi.org/10.1103/PhysRevB.47.558

[42] Kresse G, Hafner J (1994) Ab initio molecular-dynamics simulation of the liquid-metalamorphous- semiconductor transition in germanium. Phys Rev B. https://doi.org/10.1103/PhysRevB.49.14251

[43] Kresse G, Furthmüller J (1996a) Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys Rev B. https://doi.org/10.1103/PhysRevB.54.11169

[44] Kresse G, Furthmüller J (1996b) Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci. https://doi.org/10.1016/0927-0256(96),00008-0

[45] Blöchl PE (1994) Projector augmented-wave method. Phys Rev B. 50:17953–17979. https://doi.org/10.1103/PhysRevB.50.17953

[46] Joubert D (1999) From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B Condens Matter Mater Phys. https://doi.org/10.1103/PhysRevB.59.1758

[47] Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made simple. Phys Rev Lett. https://doi.org/10.1103/PhysRevLett.78.1396

[48] Perdew JP, Burke K, Ernzerhof M (1997) Erratum: generalized gradient approximation made simple. Phys Rev Lett. https://doi.org/10.1103/PhysRevLett.78.1396

[49] Ledbetter HM, Reed RP (1973) Elastic properties of metals and alloys. I. Iron, nickel, and iron-nickel alloys. J Phys Chem Ref Data. https://doi.org/10.1063/1.3253127

[50] SpringerMaterials (1998) Tellurium (Te) elastic moduli. Non-tetrahedrally bonded elements and binary compounds I, vol 41. Springer, Berlin, pp 1–2. https://doi.org/10.1007/10681727_1279

[51] Van Ingen RP, Fastenau RH, Mittemeijer EJ (1994) Laser ablation deposition of Cu–Ni and Ag–Ni films: nonconservation of alloy composition and film microstructure. J Appl Phys 76(3):1871–1883. https://doi.org/10.1063/1.357711

<figure><img src="image_18.png" /></figure>
[26] Adenis C, Langer V, Lindqvist O (1989) Reinvestigation of the structure of tellurium. Acta Crystallogr Sect C Cryst Struct Commun 45(6):941–942. https://doi.org/10.1107/S0108270188014453

[27] Westrum EF, Chou JC, Machol RE (1958) Thermodynamics of nonstoichiometric Nickel tellurides. 1. Heat Capacity and Thermodynamic Functions of the delta phase from 5 to 350 K. J Chem Phys 28(3):497–503

[28] Grønvold F, Kveseth NJ, Sveen A (1972) The Ni3+-xTe2 phases. Heat capacities in the range 298 to 900 k and transition behavior. J Chem Thermodyn 4:337–348

[29] Grønvold F (1973) Heat capacities and thermodynamic properties of Cr3Te4and Ni3Te4from 298 to 950 K. Structural order-disorder transitions. J Chem Thermodyn 5(4):545–551. https://doi.org/10.1016/S0021-9614(73)80102-8

[30] Mills KC (1974) Heat capacity of nickel and cobalt tellurides. J Chem Soc Faraday Trans 1 Phys Chem Condens Phases 70(7):2224–2231. https://doi.org/10.1039/f19747002224

[31] Tsuji T, Ishida K (1995) heat capacities of Cr5Te8 and Ni5Te8 from 325 to 930 K. Thermochim Acta 253:11–18

[32] Shukla NK, Prasad R, Roy KN, Sood DD (1990) Standard molar enthalpies of formation at the temperature 298.15 K of iron telluride (FeTe0.9) and of nickel telluride (Ni0.595Te0.405). J Chem Thermodyn 22(9):899–903. http://doi.org/10.1016/0021-9614(90)90178-S

[33] Predel B, Ruge H (1972) Bildungsenthalpien und bindungsverhältnisse in einigen intermetallischen verbindungen vom NiAs-Typ. Thermochim Acta 3(5):411–418. https://doi.org/10.1016/0040-6031(72)87055-2

[34] Jandl I, Ipser H, Richter KW (2015) Thermodynamic modelling of the general NiAs-type structure: a study of first principle energies of formation for binary Ni-containing B8 compounds. Calphad Comput Coupling Phase Diagr Thermochem 50:174–181. https://doi.org/10.1016/j.calphad.2015.06.006

[35] Prasad R, IYer VS, Venugopal V, Sundaresh V, Singh Z, Sood DD (1987) Partial and integral molar thermodynamic properties of NixTe1-x(s, x=0.595 to 0.630) alloys. J Chem Thermodyn 19:891–896

[36] Viswanathan R, Baba MS, Raj DDA, Balasubramanian R, Saha B, Mathews CK (1987) Vaporisation thermodynamics of the nickel-rich phases in the Ni–Te binary system–a high-temperature mass-spectrometric study. J Nuclear Mater 149:302–311

[37] Kutsenok IB, Geiderikh VA, Valuev IA (1980) Use of the instantaneous measuring of EMF method for the study of

thermodynamic properties of solid alloys. Vestnik Moskovskogo Universiteta Seriya 2 Khimiya 21(6):554–558

[38] Geiderikh VA, Sheveleva SN, Kutsenok IB, Krivovsheya NS (1980) Thermodynamic properties of nickel tellurides. Zhurnal Fizicheskoi Khimii 54(4):1068–1071

[39] Prasad R, Iyer VS, Venugopal V, Singh Z, Sood DD (1987) Thermodynamic Study of Ni3Te2 by E.M.F. measurements. J Chem Thermodyn 19:613–616

[40] Maekawa T, Yokokawa T (1975) The enthalpies of binary mixtures of chalcogen elements with various metals. J Chem Thermodyn 7:505–506

[41] Kresse G, Hafner J (1993) Ab initio molecular dynamics for liquid metals. Phys Rev B. https://doi.org/10.1103/PhysRevB.47.558

[42] Kresse G, Hafner J (1994) Ab initio molecular-dynamics simulation of the liquid-metalamorphous- semiconductor transition in germanium. Phys Rev B. https://doi.org/10.1103/PhysRevB.49.14251

[43] Kresse G, Furthmüller J (1996a) Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys Rev B. https://doi.org/10.1103/PhysRevB.54.11169

[44] Kresse G, Furthmüller J (1996b) Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput Mater Sci. https://doi.org/10.1016/0927-0256(96),00008-0

[45] Blöchl PE (1994) Projector augmented-wave method. Phys Rev B. 50:17953–17979. https://doi.org/10.1103/PhysRevB.50.17953

[46] Joubert D (1999) From ultrasoft pseudopotentials to the projector augmented-wave method. Phys Rev B Condens Matter Mater Phys. https://doi.org/10.1103/PhysRevB.59.1758

[47] Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made simple. Phys Rev Lett. https://doi.org/10.1103/PhysRevLett.78.1396

[48] Perdew JP, Burke K, Ernzerhof M (1997) Erratum: generalized gradient approximation made simple. Phys Rev Lett. https://doi.org/10.1103/PhysRevLett.78.1396

[49] Ledbetter HM, Reed RP (1973) Elastic properties of metals and alloys, I. Iron, nickel, and iron-nickel alloys. J Phys Chem Ref Data. https://doi.org/10.1063/1.3253127

[50] SpringerMaterials (1998) Tellurium (Te) elastic moduli. Non-tetrahedrally bonded elements and binary compounds I, vol 41. Springer, Berlin, pp 1–2. https://doi.org/10.1007/10681727_1279

[51] Van Ingen RP, Fastenau RH, Mittemeijer EJ (1994) Laser ablation deposition of Cu–Ni and Ag–Ni films: nonconservation of alloy composition and film microstructure. J Appl Phys 76(3):1871–1883. https://doi.org/10.1063/1.357711
[52] Lukas HL, Fries SG, Sundman B (eds) (2007) Computational thermodynamics: the Calphad method. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511804137

[53] Hillert M (2001) The compound energy formalism. J Alloys Compd 320(2):161–176. https://doi.org/10.1016/S0925-8388(00)01481-X

[54] Nguyen VT, Gay M, Enderby JE, Newport RJ, Howe RA (1982) The structure and electrical properties of liquid semiconductors: I. The structure of liquid NiTe2 and NiTe. J Phys Chem 15:4627–4634

[55] Dinsdale AT (1991) SGTE data for pure elements. Calphad 15(4):317–425

[56] Andersson JO, Helander T, Höglund L, Shi P, Sundman B (2002) Thermo-Calc & DICTRA, computational tools for materials science. Calphad Comput Coupl Phase Diagr Thermochem 26:273–312. https://doi.org/10.1016/S0364-5916(02)00037-8

[57] SGTE 2019 SGTE Substance database: SGTE—Scientific group thermodata Europe. https://www.sgte.net/en/neu. Accessed 15 Jan 2019

[58] Davydov AV, Rand MH, Argent BB (1995) Review of Heat Capacity Data for Tellurium. Calphad Comput Coupl Phase Diagr Thermochem 19(3):375–387

[59] Cheynet B, Chevalier PY, Fischer E (2002) Thermosuite. Calphad 26(2):167–174. https://doi.org/10.1016/S0364-5916(02)00033-0

[60] Belov G, Iorish V, Yungman VS (1999) IVTANTHERMO for Windows-database on thermodynamic properties and related software. Calphad 23(2):173–180. https://doi.org/10.1016/S0364-5916(99)00023-1

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
[52] Lukas HL, Fries SG, Sundman B (eds) (2007) Computational thermodynamics: the Calphad method. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511804137

[53] Hillert M (2001) The compound energy formalism. J Alloys Compd 320(2):161–176. https://doi.org/10.1016/S0925-8388(00)01481-X

[54] Nguyen VT, Gay M, Enderby JE, Newport RJ, Howe RA (1982) The structure and electrical properties of liquid semiconductors: I. The structure of liquid NiTe2 and NiTe. J Phys Chem 15:4627–4634

[55] Dinsdale AT (1991) SGTE data for pure elements. Calphad 15(4):317–425

[56] Andersson JO, Helander T, Höglund L, Shi P, Sundman B (2002) Thermo-Calc & DICTRA, computational tools for materials science. Calphad Comput Coupl Phase Diagr Thermochem 26:273–312. https://doi.org/10.1016/S0364-5916(02)00037-8

[57] SGTE 2019 SGTE Substance database: SGTE—Scientific group thermodata Europe. https://www.sgte.net/en/neu. Accessed 15 Jan 2019

[58] Davydov AV, Rand MH, Argent BB (1995) Review of Heat Capacity Data for Tellurium. Calphad Comput Coupl Phase Diagr Thermochem 19(3):375–387

[59] Cheynet B, Chevalier PY, Fischer E (2002) Thermosuite. Calphad 26(2):167–174. https://doi.org/10.1016/S0364-5916(02)00033-0

[60] Belov G, Iorish V, Yungman VS (1999) IVTANTHERMO for Windows-database on thermodynamic properties and related software. Calphad 23(2):173–180. https://doi.org/10.1016/S0364-5916(99)00023-1

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.