Crystallographic Calculations and First-Principles Calculations of Heterogeneous Nucleation Potency of γ-Fe on La2O2S Particles

Abstract

Rare earth (RE) inclusions with high melting points as heterogeneous nucleation in liquid steel have stimulated many recent studies. Evaluating the potency of RE inclusions as heterogeneous nucleation sites of the primary phase is still a challenge. In this work, the edge-to-edge matching (E2EM) model was employed to calculate the atomic matching mismatch and predict the orientation relationship between La₂O₂S and γ-Fe from a crystallographic point of view. A rough orientation relationship (OR) was predicted with the minimum values of ${\mathrm{f}}_{\mathrm{r}}=9.43\%$ and ${\mathrm{f}}_{\mathrm{d}}=20.72\%$ as follows: ${\left[2\overline{1}\overline{1}0\right]}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left[100\right]}_{\mathsf{\gamma }-\mathrm{Fe}}$ and ${\left(000\overline{3}\right)}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left(00\overline{2}\right)}_{\mathsf{\gamma }-\mathrm{Fe}}$. The interface energy and bonding characteristics between La₂O₂S and γ-Fe were calculated on the atomic scale based on a crystallographic study using the first-principles calculation method. The calculations of the interface energy showed that the S-terminated and La(S)-terminated interface structures were more stable. The results of difference charge density, electron localization function (ELF), the Bader charges and the partial density of states (PDOS) study indicated that the La(S)-terminated interface possessed metallic bonds and ionic bonds, and the S-terminated interface exhibited metallic bond and covalent bond characteristics. This work addressed the stability and the characteristics of the La₂O₂S/γ-Fe interface structure from the standpoint of crystallography and energetics, which provides an effective theoretical support to the study the heterogeneous nucleation mechanism. As a result, La₂O₂S particles are not an effective heterogeneous nucleation site for the γ-Fe matrix from crystallography and energetics points of view.

1. Introduction

Previous studies showed that adding the rare earth element La/Ce into steels resulted either in a decrease in the proportion of the columnar grains or the average grain size of the equiaxed grains or both in the case of as-cast austenitic [1,2,3], ferritic [4,5,6] and other engineering steels [7,8]. The most generally accepted effect mechanisms of La/Ce are the heterogeneous nucleation effect of the primary δ-ferrite/γ-austenite on La/Ce-containing inclusions with high melting points, including oxides, oxysulfides and sulfides [1,2,4,5,6,9], or/and the solute effects [2,10] attributed to the solute enrichment or depletion in the liquid phase ahead of the solid-liquid interface. In the past decades, the two-dimensional lattice disregistry model [11] has been widely adopted to evaluate the heterogeneous nucleation potency. A potency evaluation criterion based on the two-dimensional lattice disregistry model was proposed [11]: when the lattice disregistry value (δ-value) is below 6%, the potency is the most effective, and the grains can be significantly refined; when the δ-value is between 6% and 12%, the potency is moderately effective; when the δ-value is above 12%, the potency is poor, and the grains cannot be refined. Based on this model, researchers reported the heterogeneous nucleation potency of La/Ce oxides, oxysulfides and sulfides for the primary δ-ferrite and/or primary γ-austenite during steel solidification, examined using the two-dimensional disregistry model [2,6,12,13]. However, for a given La/Ce-containing inclusion phase and the primary δ-ferrite/γ-austenite, the δ-values they obtained were highly scattered. Relevant studies have been reviewed by the authors [14]. The heterogeneous nucleation potency of La/Ce oxides, oxysulfides and sulfides for the primary δ-ferrite/γ-austenite was inconsistent in relevant studies. The discrepancy might have resulted from different selection standards. Put another way, the selecting criteria of low index planes and directions for δ-value calculation were not defined based on the two-dimensional lattice disregistry model. Yu and co-workers [15] calculated the lattice disregistry between La₂O₂S and γ-austenite in an Fe-43Ni expansion alloy and reported the lowest δ-value of 5.42% between the $\left(0001\right)$ crystal plane of La₂O₂S and the $\left(100\right)$ crystal plane of γ-austenite. This implied a high heterogeneous nucleation potency of La₂O₂S for γ-austenite. However, there was no conclusive experimental evidence to verify it. Robert Tuttle’s work [16,17,18] showed that the addition of lanthanum can refine the solidification structure, but not all lanthanum inclusions act as heterogeneous nucleus of γ-austenite. The heterogeneous nucleation potency of La-containing inclusion phases for the primary γ-austenite needs further investigation.

Recently, the edge-to-edge matching (E2EM) model [19] originally developed by Zhang and Kelly was validated to be a powerful crystallography model in interpreting the heterogeneous nucleation potency in the grain refinement during solidification and predicting effective heterogeneous nucleation agents in Al, Mg and their alloys [20,21,22]. The major advantage of the E2EM model over the two-dimensional lattice disregistry model is that it focuses on actual atom matching rather than lattice matching. It can not only calculate the atomic matching between two phases but also predict the favorable orientation relationships (ORs) between the two phases (see refs. [18,19,20,22] for related details and calculation method). First-principles calculation based on density functional theory is a very effective method for studying alloy interfaces [23,24,25]. The density functional theory (DFT) method is widely applied to investigate the material interface of the magnesium matrix [23,26] and aluminum matrix [24]. Moreover, for the interface between compounds (such as LaAlO₃ [27], La₂O₃ [25], Al₂O₃ [28], TiC [29] and TiN [30]) and the δ-ferrite/γ-austenite matrix, the DFT calculations are also employed to probe the interfacial characteristics and the mechanism of heterogeneous nucleation. In this paper, La₂O₂S, which is the typical oxysulfide, was selected. The present work aims to clarify the heterogeneous nucleation potency of γ-Fe on La₂O₂S particles, combining a crystallographic calculation using the edge-to-edge matching model with the first-principles calculations from an energetics point of view.

2. Calculations

2.1. Crystallographic Calculations Using the E2EM Model

To investigate the heterogeneous nucleation potency of La₂O₂S for γ-austenite, the E2EM model was adopted for crystallographic calculations. The interatomic spacing misfit (f_r) along the close-packed (CP) or nearly CP matching rows, as well as the interplanar spacing mismatch (f_d) of the CP or nearly CP planes between La₂O₂S and γ-austenite, were obtained based on the crystal structures, lattice parameters and atomic positions in the two phases. Table 1 presents the crystal structure, space group, lattice parameters and atomic positions of γ-austenite (γ-Fe) and La₂O₂S.

The γ-Fe phase possesses a face-centered cubic (FCC) lattice structure with four atoms in each unit cell. By contrast, the La₂O₂S phase has a primitive hexagonal structure, and each unit cell contains two lanthanum atoms, two oxygen atoms and one sulfur atom. The γ-Fe lattice parameters were [31]: a = 0.36544 nm; space group = $\mathrm{Fm}\overline{3}m$. The La₂O₂S lattice parameters were [32]: a = 0.4035 nm; c = 0.6914 nm; space group = $\mathrm{P}\overline{3}m1$.

2.2. First-Principles Methodology

Effective heterogeneous nucleation must satisfy the crystallographic match and interface stability between the nucleating agent and the nucleating phase. The plane-wave, pseudo-potential density functional theory (DFT) method implemented in the VASP package [33] was used to investigate the stability of La₂O₂S/γ-Fe interface. The interaction between ion and valence electrons was described by the projector augmented wave (PAW) [34]. The exchange correlation energy between electrons was employed using the local density approximation (LDA) [35] and the generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) functional [36]. Due to the presence of La atoms in the system, strong correlation effects had to be considered. Thus, the LDA+U/GGA+U functions were adopted as the exchange correlation functions. Moreover, the effective value for U was 6 eV [37]. Since Fe atoms are magnetic, the spin-polarization method was employed. The Kohn-Sham equation was solved with self-consistent field (SCF), and the electronic minimization process was conducted to reach the ground state [36].

In order to verify the methodology proposed in this work, the bulk properties of La₂O₂S and γ-Fe were investigated, and the results were compared with published experimental data. To estimate the necessary layer numbers of bulk-like slabs N used in the subsequent interface systems, the surface properties of slabs with different layer numbers were studied. The interface properties of different atomic terminals, including adhesion work, interface energy and interface bonding properties, were calculated to further investigate the interface stability.

In the structural relaxation, the plane-wave cutoff E_cut was 450 eV. In the γ-Fe and La₂O₂S bulk phase optimization, the Brillouin zone integral was divided by 9 × 9 × 9 and 10 × 10 × 5 Monkhorst-Pack k-point grids (centered on the gamma point), respectively, which was proven to be sufficient in our convergence test. The k-points sampling grids of the surface and interface were generated in the Monkhorst-Pack method of 0.04 grid density by VASPKIT [38]. A vacuum gap of 15 Å was inserted to eliminate periodical influence in the Z direction. According to the previous DFT calculations regarding other condensed matter interfaces [24,28,39], a vacuum depth of 15 Å is thick enough to generate accurate results. The energy changes in the structural optimization process converged to 1.0 × 10⁻⁴ eV/atom, and the residual Hellmann–Feynman force on each relaxed atom reached 0.03 eV/Å [27,28,29].

3. Results and Discussion

3.1. Nucleation Crystallography of γ-Fe on La₂O₂S Particles

Based on crystal structures, lattice parameters and atom positions, the E2EM calculations of the CP or nearly CP rows and the CP or nearly CP planes of γ-Fe and La₂O₂S were obtained, as shown in Table 2. It can be seen that the smaller the atomic spacing, the denser the arrangement; the larger the face spacing, the denser the arrangement. The CP atomic rows were straight, pseudo-straight or zigzag, denoted by SS, PS and ZZ, respectively. Then, based on the identified CP rows and CP planes, the values of f_r and f_d between La₂O₂S and γ-Fe were calculated, as listed in Table 3. It can be seen that f_r was 9.43%, and the minimum f_d was 4.58%.

Because the matching rows must be on the corresponding matching planes, a favorable orientation relationship (OR) between La₂O₂S and γ-Fe could be predicted based on the matching directions and matching planes. The predicted rough OR between La₂O₂S and γ-Fe could be expressed as ${\left[2\overline{1}\overline{1}0\right]}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left[100\right]}_{\mathsf{\gamma }-\mathrm{Fe}}$ and ${\left(000\overline{3}\right)}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left(00\overline{2}\right)}_{\mathsf{\gamma }-\mathrm{Fe}}$, and the values of f_r and f_d were 9.43% and 20.72%, respectively. Based on this orientation relationship, the properties of the heterogeneous interface between La₂O₂S and γ-Fe were calculated by the first principles. The effectiveness of La₂O₂S as the heterogeneous nucleation site of γ-Fe was discussed from the point of view of energy.

3.2. Bulk and Surface Properties

3.2.1. Bulk Calculations

Table 4 and Table 5 present the optimized γ-Fe and La₂O₂S lattice constants (a and c), bulk module (B₀) and phase equilibrium volume per atom (V₀) by LDA and GGA-PBE function. For γ-Fe bulk, the calculated constants by GGA-PBE function were a = b = c = 3.636 Å, 0.5% below the experimental values [31] a = b = c = 3.654 Å. For La₂O₂S phase, the calculated lattice constants were a = b = 4.134 Å, and c = 7.073 Å. They are roughly consistent with the experimental values [32] a = b = 4.035 Å and c = 6.914 Å. The results of GGA-PBE functional also agreed well with the reported calculations in Table 4 and Table 5. This ensured the reliability of the presented calculations.

3.2.2. Surface Convergence Test

In view of the lattice match of the interface, γ-Fe$\left(00\overline{2}\right)$ and La₂O₂S$\left(000\overline{3}\right)$ surfaces were selected to construct their interface. The surface model should be thick enough to ensure the bulk-like character interiors so as to obtain accurate surface properties. It was, therefore, necessary to perform a convergence test on the slabs to construct the interface calculations. The 5-layer, 7-layer, 9-layer and 11-layer slabs were respectively modeled for the convergence test of γ-Fe$\left(00\overline{2}\right)$ through the method presented in refs. [25,28,37]. The surface energy can be expressed by [30]:

(1)$E_{\mathrm{surf}}=\left(E_{\mathrm{slab}}-\left(N_{\mathrm{slab}}/N_{\mathrm{bulk}}\right)E_{\mathrm{bulk}}\right)/2A,$

where $E_{\mathrm{slab}}$ is the total energy of the slab; $E_{\mathrm{bulk}}$ is the total energy of bulk after energy optimization; $A$ is the corresponding surface area; and $N_{\mathrm{slab}}$ and $N_{\mathrm{bulk}}$ are the number of atoms in the slab and bulk structure, respectively. Table 6 gives the surface energy of γ-Fe stoichiometric slab corresponding to the number of atom layers. Obviously, the γ-Fe$\left(00\overline{2}\right)$ surface energy converged to 2.18 J/m² when the surface model contained more than five atomic layers. This figure is close to the experimental value 1.95 J/m² reported in ref. [39] Therefore, the five-layer γ-Fe$\left(00\overline{2}\right)$ surface slab was adopted to ensure the bulk-like interior in the following calculations.

For the La₂O₂S$\left(000\overline{3}\right)$ slab, there were three different kinds of surface terminations (La-, O- and S-terminated). Both La- and O-terminated surfaces had two nonequivalent surfaces with different sub-layer atoms. The sub-layer atoms were O or S for the La-terminated surface and La or O for O-terminated. Symmetric slabs were adopted in the surface convergence tests to eliminate the spurious dipole effects. Therefore, there were five possible types of slab surface, as shown in Figure 1.

Table 7, Table 8 and Table 9 show the La₂O₂S$\left(000\overline{3}\right)$ surface relaxations of the different atomic layers with different terminals before and after the change of the atomic layer spacing accounted for the percentage of the interlayer spacing of the bulk. It can be seen that the effects of atomic relaxations were mainly localized within the top three atomic layers. With the increasing atomic layers, the slab interiors became gradually bulk-like in nature. In the case of the 14-layer, 13-layer, 15-layer, 12-layer and 11-layer structure, respectively, the La(O)-terminated, La(S)-terminated, O(La)-terminated, O(O)-terminated and S-terminated surface interlayer relaxations were well converged. Therefore, the corresponding slabs were applied to simulate their interface structure. The results in Table 7, Table 8 and Table 9 demonstrate that the outermost and second interlayer relaxations for S-terminated surface were smallest among the five types of surface, followed by the La(S)-terminated surface. It was observed that the outermost interlayer distances of S-terminated and La(O)-terminated surfaces with eleven and thirteen atomic layers were reduced/enlarged by 1.037% and 5.663%, and the second interlayer distances were enlarged/reduced by 2.080% and 2.121%, respectively. Thus, the S-terminated and La(O)-terminated surfaces tended to be more stable, with no surface reconstruction.

3.2.3. Surface Stability of La₂O₂S$\left(000\overline{3}\right)$

The atomic species of La₂O₂S surface have significant impact on the interfacial stability. Based on the surface convergence test result, the stability of different terminal surfaces was investigated. For surface energy of La₂O₂S, characterized by non-stoichiometry, La, O and S chemical potentials needed to be considered comprehensively to accurately determine the stability in the experimental environment. The surface energy of La₂O₂S was calculated by Equation (2) [23,25,26,27,28,42]:

(2)${\delta }_{{\text{La}}_{\text{2}}{\text{O}}_{\text{2}}\text{S}\left(000\overline{3}\right)}=\left[E_{\text{slab}}-N_{\text{La}}{\mu }_{\text{La}}-N_{\text{O}}{\mu }_{\text{O}}-N_{\text{S}}{\mu }_{\text{S}}+PV-TS\right]/2A,$

where $E_{\mathrm{slab}}$ is total energy of the corresponding slab, $N_{\mathrm{i}}$ is the number of atomic species i (i = La, O and S) in the La₂O₂S slab, ${\mu }_{\mathrm{i}}$ is chemical potentials of atomic species i (i = La, O and S) at ground state and $A$ is surface area. The calculations were carried out at 0 K and condensed-matter states, so the PV and TS terms can be neglected [28].

Combining the chemical potentials of elements La, O and S of the La₂O₂S slab in equilibrium with those of bulk La₂O₂S, the relationship between the chemical potential of La₂O₂S bulk ${\mu }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}^{\mathrm{bulk}}$, its 0 K formation heat ($\mathsf{\Delta }{H}_f^0\left({\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}\right)$) and chemical potentials of the La bulk, O₂ gas and S bulk (${\mu }_{\mathrm{La}}^{\mathrm{bulk}}$, $\frac{1}{2}{\mu }_{O_2}^{\mathrm{gas}}$ and ${\mu }_{\mathrm{S}}^{\mathrm{bulk}}$) are expressed by [43]:

(3)${\mu }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}^{\mathrm{bulk}}=2{\mu }_{\mathrm{La}}+2{\mu }_{\mathrm{O}}+{\mu }_{\mathrm{S}},$

(4)${\mu }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}^{\mathrm{bulk}}=2{\mu }_{\mathrm{La}}^{\mathrm{bulk}}+{\mu }_{O_2}^{\mathrm{gas}}+{\mu }_{\mathrm{S}}^{\mathrm{bulk}}+\mathsf{\Delta }{H}_f^0\left({\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}\right),$

where ${\mu }_{\mathrm{i}}$ is chemical potentials of atomic species i (i = La, O and S), the actual value of which is located in the surface within a certain range and is determined by experimental conditions. The upper bound of ${\mu }_{\mathrm{i}}$ (i = La, O and S) at the surface is the corresponding potential in pure La bulk, gas phase of O₂ and pure S bulk: ${\mu }_{\mathrm{La}}\le {\mu }_{\mathrm{La}}^{\mathrm{bulk}}$, ${\mu }_{\mathrm{O}}\le \frac{1}{2}{\mu }_{O_2}^{\mathrm{gas}}$ and ${\mu }_{\mathrm{S}}\le {\mu }_{\mathrm{S}}^{\mathrm{bulk}}$, which is the basis for forming stable slab structure. The range of ${\mu }_{\mathrm{La}}-{\mu }_{\mathrm{La}}^{\mathrm{bulk}}$ was obtained:

(5)$\mathsf{\Delta }{\mu }_{\mathrm{La}}={\mu }_{\mathrm{La}}-{\mu }_{\mathrm{La}}^{\mathrm{bulk}}\le 0,$

(6)$\mathsf{\Delta }{\mu }_{\mathrm{O}}={\mu }_{\mathrm{O}}-\frac{1}{2}{\mu }_{O_2}^{\mathrm{gas}}\le 0,$

(7)$\mathsf{\Delta }{\mu }_{\mathrm{S}}={\mu }_{\mathrm{S}}-{\mu }_{\mathrm{S}}^{\mathrm{bulk}}\le 0,$

Based on Equations (3), (5) and (7), the surface energy (Equation (2)) can be rewritten as:

(8)${\delta }_{{\text{La}}_{\text{2}}{\text{O}}_{\text{2}}\text{S}\left(000\overline{3}\right)}=\left[\begin{array}{c}{E}_{\text{slab}}-\frac{N_{\text{O}}}{2}{\mu }_{{\text{La}}_{\text{2}}{\text{O}}_{\text{2}}\text{S}}^{\text{bulk}}+\left(N_{\text{O}}-N_{\text{La}}\right){\mu }_{\text{La}}^{\text{bulk}}+\left(N_{\text{O}}-N_{\text{La}}\right)\mathsf{\Delta }{\mu }_{\text{La}}\\ +\left(\frac{N_{\text{O}}}{2}-N_{\text{S}}\right){\mu }_{\text{S}}^{\text{bulk}}+\left(\frac{N_{\text{O}}}{2}-N_{\text{S}}\right)\mathsf{\Delta }{\mu }_{\text{S}}\end{array}\right]/2A,$

where ${\mu }_{\mathrm{La}}^{\mathrm{bulk}}$ and ${\mu }_{\mathrm{S}}^{\mathrm{bulk}}$ refers to the single atomic energy in the La and S bulk; ${\mu }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}^{\mathrm{bulk}}$ is equal to the energy of La₂O₂S bulk per formula unit (f.u.); and $N_{\mathrm{i}}$ (i = La, O and S) has the same meaning as in Equation (2). For La₂O₂S bulk as a heterogeneous nucleus, the actual solidification process was in a S-rich environment, so the S chemical potential was approximated to zero ($\mathsf{\Delta }{\mu }_{\mathrm{S}}=0$). The range of the $\mathsf{\Delta }{\mu }_{\mathrm{La}}$ value combined with the formula available for La chemical potential (Equations (4)–(6)) can be derived as:

(9)$\mathsf{\Delta }{H}_f^0\left({\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}\right)\le \mathsf{\Delta }{\mu }_{\mathrm{La}}\le 0,$

The formation heat $\mathsf{\Delta }{H}_f^0\left({\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}\right)$ was calculated to be −5.12 eV. The range of $\mathsf{\Delta }{\mu }_{\mathrm{La}}$ was, therefore, derived to be:

(10)$-5.12 \mathrm{eV}\le \mathsf{\Delta }{\mu }_{\mathrm{La}}\le 0 \mathrm{eV},$

Figure 2 shows La₂O₂S$\left(000\overline{3}\right)$ surface energy varying with the variation of the chemical potential of $\mathsf{\Delta }{\mu }_{\mathrm{La}}$. Clearly, the surface energy of these surfaces was linearly correlated with $\mathsf{\Delta }{\mu }_{\mathrm{La}}$. Generally, the lower the surface energy is, the more stable the surface structure is. In the whole range of La chemical potential, the surface energy of the La(S)-terminated surface was the lowest, and the values decreased with the increases of $\mathsf{\Delta }{\mu }_{\mathrm{La}}$ in the range of −0.74 J/m²~−6.29 J/m². By contrast, the surface energy of La(O)-terminated, O(La)-terminated and S-terminated surfaces remained constant at 4.79 J/m², 0.68 J/m² and 0.06 J/m², respectively. It should be pointed out that the constant value was due to the zero value of ($N_{\mathrm{O}}-N_{\mathrm{La}}$) in Equation (10). Therefore, the La(S)-terminated surface was the most stable one over the entire range of chemical potential.

3.3. Interface Properties between La₂O₂S and γ-Fe

3.3.1. Structures of La₂O₂S/γ-Fe Interface

Based on the surface convergence test results and the orientation relationship predicted by edge-to-edge matching crystallography theory, five interface models between La₂O₂S and γ-Fe were constructed. They were denoted by La(O)-terminated, La(S)-terminated, O(La)-terminated, O(O)-terminated and S-terminated La₂O₂S/γ-Fe interfaces, as shown in Figure 3. A supercell geometry where a five-layer γ-Fe$\left(00\overline{2}\right)$ slab was placed under a corresponding, symmetric atomic terminals La₂O₂S$\left(000\overline{3}\right)$ slab was used to simulate the La₂O₂S/γ-Fe interface. For interfacing the close-packed planes, the slabs were rotated about an axis normal to the interface so as to align the close-packed directions, resulting in an orientation relationship of ${\left[2\overline{1}\overline{1}0\right]}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left[100\right]}_{\mathsf{\gamma }-\mathrm{Fe}}$ and ${\left(000\overline{3}\right)}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left(00\overline{2}\right)}_{\mathsf{\gamma }-\mathrm{Fe}}$. To accommodate the periodic boundary conditions inherent in a supercell calculation, the coherent interface approximation [44] was invoked, and the average of the lattice parameters for La₂O₂S and γ-Fe was used. The interface structures were fully relaxed in all dimensions by minimizing the Hellmann-Feynman force on every atom.

3.3.2. Interface Stability

Binding strength of the interface can be predicted by the work of adhesion ($W_{\mathrm{ad}}$), which is defined as the reversible work required to separate the interface between A and B phase to form two free surfaces. The larger the $W_{\mathrm{ad}}$ is, the stronger the binding force of interface atoms is. The $W_{\mathrm{ad}}$ of the interface can be calculated by Equation (11) [23,24,25,26,27,28,30]:

(11)$W_{\mathrm{ad}}=\left(E_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}+E_{\mathsf{\gamma }-\mathrm{Fe}}-E_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}/\mathsf{\gamma }-\mathrm{Fe}}\right)/A,$

where $E_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}/\mathsf{\gamma }-\mathrm{Fe}}$ is the total energy of interface system; $E_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ and $E_{\mathsf{\gamma }-\mathrm{Fe}}$ are the total energy of isolated La₂O₂S and γ-Fe slabs in the same interface system when one of the slabs is kept and other is replaced with vacuum, respectively; and $A$ is the interface area. Generally, a more positive $W_{\mathrm{ad}}$ means a larger binding strength.

Table 10 presents the optimal $W_{\mathrm{ad}}$ for the relaxed geometries of the five La₂O₂S/γ-Fe interfaces. The $W_{\mathrm{ad}}$ of the five La₂O₂S/γ-Fe interfaces in descending order was: S-terminated > La(S)-terminated > La(O)-terminated > O(O)-terminated > O(La)-terminated. This indicates that the S-terminated interface had the largest binging strength, followed by the La(S)-terminated interface. This partly resulted from the bonding type between the interface atoms. In addition, the S-terminated and La(S)-terminated surface structures were more stable than others with a smaller surface energy. These demonstrate that it is more likely that the S-terminated and La(S)-terminated interfaces are maintained.

In addition to the work of adhesion ($W_{\mathrm{ad}}$), interface energy is another major indicator for measuring interface stability. In general, the lower the interface energy becomes, the more stable the interface structure is. The interface energy of La₂O₂S/γ-Fe was calculated from the viewpoint of thermodynamics, and the interface stability of the different terminated interface configurations was quantitatively analyzed. The results can provide theoretical data for the heterogeneous nucleation mechanism of La₂O₂S on the γ-Fe substrate. According to the calculation method provided in refs. [23,26,30], the interface energy can be derived as follows:

(12)${\gamma }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}/\mathsf{\gamma }-\mathrm{Fe}}=\left[\begin{array}{c}{E}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}/\mathsf{\gamma }-\mathrm{Fe}}-\frac{N_{\mathrm{O}}}{2}{\mu }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}^{\mathrm{bulk}}+\left(N_{\mathrm{O}}-N_{\mathrm{La}}\right){\mu }_{\mathrm{La}}^{\mathrm{bulk}}+\left(N_{\mathrm{O}}-N_{\mathrm{La}}\right)\mathsf{\Delta }{\mu }_{\mathrm{La}}\\ +\left(\frac{N_{\mathrm{O}}}{2}-N_{\mathrm{S}}\right){\mu }_{\mathrm{S}}^{\mathrm{bulk}}-N_{\mathrm{Fe}}{\mu }_{\mathrm{Fe}}^{\mathrm{bulk}}\end{array}\right]/A-{\delta }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}-{\delta }_{\mathsf{\gamma }-\mathrm{Fe}},$

where $E_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}/\mathsf{\gamma }-\mathrm{Fe}}$ is the total energy of interface system; ${\mu }_{\mathrm{Fe}}^{\mathrm{bulk}}$ is equal to the energy of γ-Fe bulk per formula unit (f.u.); ${\delta }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ and ${\delta }_{\mathsf{\gamma }-\mathrm{Fe}}$ denote the surface energy of La₂O₂S and γ-Fe free surfaces in the same interface system; and the other parameters are consistent with those in Equation (8).

Figure 4 shows the interface energy of La₂O₂S/γ-Fe interfaces with the increase of $\mathsf{\Delta }{\mu }_{\mathrm{La}}$. The La₂O₂S/γ-Fe structures of S-terminated and La(S)-terminated interfaces, which have smaller interface energy, were more stable than others. It, therefore, can be concluded that S-terminated and La(S)-terminated interfaces are most likely to form. Within the range of $-5.12 \mathrm{eV}\le \mathsf{\Delta }{\mu }_{\mathrm{La}}\le -4.69 \mathrm{eV}$, the interface energy of the S-terminated interface structure was smaller than that of the La(S)-terminated interface structure. This indicates that the S-terminated interface had a more stable structure. In the case of $-4.69 \mathrm{eV}\le \mathsf{\Delta }{\mu }_{\mathrm{La}}\le 0 \mathrm{eV}$, the interface energy of S-terminated interface structure was larger than that of La(S)-terminated interface structure. Thus, the La(S)-terminated interface structure was more stable.

3.3.3. Electronic Structure and Bonding

The bonding strength of the interface is closely related to the bonding type of the interface atoms, as well as the electronic structure of the interface. To better understand interface atomic bonding, the difference charge density, the electron localization function (ELF), the Bader charges and the partial density of states (PDOS) were employed. The electronic structure and bonding of the more stable S-terminated and La(S)-terminated interface structures were compared to explore the interface stability.

The difference charge density can be used to analyze the transfer of interatomic charge. The charge accumulation and depletion regions were ascertained. The charge density difference was calculated by the equation below [28,30]:

(13)$\mathsf{\Delta }\rho ={\rho }_{\mathrm{int}}-{\rho }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}-{\rho }_{\mathsf{\gamma }-\mathrm{Fe}},$

where ${\rho }_{\mathrm{int}}$ is the total charge density of La₂O₂S/γ-Fe interface system; and ${\rho }_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ and ${\rho }_{\mathsf{\gamma }-\mathrm{Fe}}$ are the charge density of corresponding isolated La₂O₂S slab and γ-Fe slab in the same interface structure, respectively.

Figure 5 shows the charge density differences around La(S)-terminated and S-terminated interface structures. There was an obvious charge transfer or redistribution in or near the interface region. By comparing the three-dimensional charge density difference of both interface structures, it was found that the electron accumulation was mainly concentrated below the La and S atoms, while the electron depletion was mostly on the upward side of the Fe atoms. The interface electron transfer was mainly contributed from the Fe–La and Fe–S atom pairs. Therefore, the interfacial bonds mainly came from the electronic interactions of the Fe–La and Fe–S atom pairs.

The localized distribution characteristics of electrons can be characterized by electron localization function (ELF), which can also confirm the bond types. The ELF was calculated by Equation (14) [45]:

(14)$ELF=1/\left[1+{\left(D_r/D{h}_r\right)}^2\right],$

where $D_r$ is the real electron gas density; and $D{h}_r$ is the uniform electron gas density. The values of ELF are between 0 and 1. The upper limit ELF = 1 corresponds to perfect localization; ELF = 1/2 corresponds to electron gas-like pair probability. While the value ELF = 0 means that electron might be entirely delocalized (or that there are no electrons).

Figure 6 shows the ELF images of the La(S)-terminated and S-terminated interfaces when the isosurface level value was 0.1 eV/Bohr³. Chemical bonds were formed between the different types of interface atoms, and the charge was redistributed after the interface was formed. As shown in Figure 6a,b, the electronic spatial distribution in the atomic layer near the interface differed greatly from that in the innermost atomic layer. Moreover, the extranuclear electron orbitals of the La and S atoms had obvious interactions with that of the Fe atoms at the interface. It can be observed from Figure 6c that there were regions with ELF ≈ 0.3 at the interface. In this case, there was low electron localization, and an ionic bond was formed between the interfaces. Combined with Bader charge analysis, the electrons of La atoms at the interface were transferred to the Fe side, and a metal bond was formed between the interfaces. As shown in Figure 6d, there were green regions with ELF ≈ 0.5 between the S atoms and Fe atoms at the interface. Because of the strong electronegativity of S atoms, the electron loss was mostly concentrated on the Fe atoms at the S-terminated interface, while the electron accumulation was mainly on the S atoms, which contributed to the formation of covalent bonds at the S-terminated interface.

To further investigate the bonding characteristics, the partial density of states (PDOS) was also calculated, as shown in Figure 7. The dashed line indicates the position of the Fermi level. In addition, the DOS near both interfaces was not zero at the Fermi energy, indicating certain metallic characteristics. Obviously, the DOS of the La and Fe atoms at the first layer of the La(S)-terminated interface peaked at Fermi level. The metal characteristics of the La(S)-terminated interface were contributed by La atoms and Fe atoms, while those of the S-terminal interface were made by the Fe atomic contributions. As can be seen from Figure 7a, there was a weak electronic interaction between Fe and La atoms at the La(S)-terminated interface. By comparing the PDOS of different layers of Fe and La atoms, it was found that the d-orbital electrons of the subsurface and bulk La atoms had a more instinct localized feature than the interface La atoms. This suggests that the inner La atoms were more covalent. Moreover, several overlapped peaks appeared between the third layer of La atoms and the second layer of S atoms near the interface, indicating the existence of covalent bonds near the interface. It can be observed from Figure 7b that there was strong electron interaction between Fe atoms and S atoms at the S-terminated interface. In addition, there were some overlapped peaks between Fe-d orbital and the S-p orbital, forming a covalent bond between them at the S-terminated interface. Some overlapped peaks also appeared between the second layer of La atoms and the third layer of O atoms near the interface, which contributed to the formation of covalent bonds near the interface.

3.3.4. Heterogeneous Nucleation Potency of γ-Fe on La₂O₂S Particles

To ensure sufficient atomic matching to retain a coherent or semi-coherent interface between the two phases, empirical criteria of the f_r and f_d values were both less than 10 pct in an as-cast light alloy based on the E2EM model [18,19,20,22]. According to the crystallographic calculation results in Table 3 and the predicted rough OR between La₂O₂S and γ-Fe, the minimum values of f_r and f_d in construction of the OR were 9.43% and 20.72%, respectively. This implies that, from a crystallographic point view, La₂O₂S is not an effective nucleant for the nucleation of γ-Fe.

Although all of the first-principles calculations were performed at 0 K, the results proved the experimental phenomena of solid-phase systems and solid–liquid interfaces at high temperatures [46,47]. To make La₂O₂S the effective heterogeneous nucleus of γ-Fe during the solidification, the solid–solid interface energy of La₂O₂S/γ-Fe should be between 0 J/m² and 0.204 J/m² from the viewpoint of thermodynamics [37], for which the theoretical value of liquid–solid interface energy is 0.241 J/m² [48], as shown in Figure 4. As can be seen from Figure 4, the interface energy of La₂O₂S/γ-Fe could not meet the interface energy requirements as a heterogeneous nucleation interface. It did not comply with the thermodynamic conditions that make La₂O₂S particles effective as potent heterogeneous nucleation substrates for the γ-Fe matrix. In addition, factors, such as the size and distribution of La₂O₂S particles, also impact whether La₂O₂S can be an effective heterogeneous nucleus of γ-Fe during solidification. Future research will focus on the impact of such factors. The orientation relationship of crystallography will be verified by both electron backscatter diffraction and transmission electron microscopy.

4. Conclusions

• With the crystallographic calculations using the E2EM model, a rough HCP/FCC orientation relationship (OR) was predicted between La₂O₂S particles and the γ-Fe matrix, the values of the interatomic spacing misfit and the interplanar spacing mismatch of which were 9.43% and 20.72%, respectively. This OR was ${\left[2\overline{1}\overline{1}0\right]}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left[100\right]}_{\mathsf{\gamma }-\mathrm{Fe}}$ and ${\left(000\overline{3}\right)}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$∥${\left(00\overline{2}\right)}_{\mathsf{\gamma }-\mathrm{Fe}}$.

• For γ-Fe/La₂O₂S$\left(000\overline{3}\right)$ interfaces, the La(S)-terminated and S-terminated interfaces were more stable than others with smaller interface energy. In the range of $-5.12 \mathrm{eV}\le \mathsf{\Delta }{\mu }_{\mathrm{La}}\le -4.69 \mathrm{eV}$, the S-terminated interface structure was the most stable. Nevertheless, for the range of $\mathsf{\Delta }{\mu }_{\mathrm{La}}$ from −4.69 eV to Fermi energy level, the La(S)-terminated interface structure had the best stability. The S-terminated interface had the largest $W_{\mathrm{ad}}$, which was 4.403 J/m², followed by the La(S)-terminated interface, which was 3.372 J/m².

• According to the bonding analysis of the electronic structure and charge transfer on the two different interfaces, the interface interaction of the S-terminated interface exhibited metallic and ionic bonds, whereas the La(S)-terminated interface exhibited metallic and covalent bonds.

• Crystallographic calculations indicated that the interatomic spacing misfit and the interplanar spacing mismatch were large enough, compared with other reported efficient nucleating systems, leading to weak potency of La₂O₂S as a nucleation site for γ-Fe. The first-principles calculations revealed that the interfacial energy of La₂O₂S/γ-Fe cannot meet the interface energy requirements as a heterogeneous nucleation interface. Thus, La₂O₂S particles are noneffective heterogeneous nucleation substrates for a γ-Fe matrix from crystallography and thermodynamics viewpoints.

Author Contributions

Conceptualization, Y.Z. and Y.J.; methodology, Y.Z., Y.J. and J.Q.; investigation, Y.Z., H.X. and Y.J.; resources, Y.J., H.R. and J.Q.; data curation, Y.Z. and Y.J.; writing—original draft preparation, Y.Z.; writing—review and editing, Y.Z., Y.J., Y.L. and H.R.; visualization, Y.Z. and Y.J.; supervision, Y.J. and H.R.; project administration, Y.J. and H.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China, grant number 51761034.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (grant no. 51761034).

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Slab configuration: (a) γ-Fe[Forumla omitted. See PDF.] surface and La2O2S[Forumla omitted. See PDF.] surfaces terminated with: (b) La(O)-terminated; (c) La(S)-terminated; (d) O(La)-terminated; (e) O(O)-terminated; (f) S-terminated. Coordinate dimension a, b and c represent x-axis, y-axis and z-axis respectively.

Figure 2. Surface energies for the La2O2S[Forumla omitted. See PDF.] surfaces as [Forumla omitted. See PDF.]. The corresponding number of atomic layers adopted in the surface configuration are 14, 13, 15, 12 and 11 for La(O)-terminated, La(S)-terminated, O(La)-terminated, O(O)-terminated and S-terminated surfaces, respectively.

Figure 3. Supercell structures before relaxation: (a) La(O)-terminated La2O2S/γ-Fe interface; (b) La(S)-terminated La2O2S/γ-Fe interface; (c) O(La)-terminated La2O2S/γ-Fe interface; (d) O(O)-terminated La2O2S/γ-Fe interface; (e) S-terminated La2O2S/γ-Fe interface. Coordinate dimension a, b and c represent x-axis, y-axis and z-axis respectively.

Figure 4. Interface energies for the La2O2S/γ-Fe interfaces with different terminated interfaces as [Forumla omitted. See PDF.].

Figure 5. Three-dimensional charge density difference: (a) the La(S)-terminated La2O2S/γ-Fe interface; (b) S-terminated La2O2S/γ-Fe interface. The yellow regions indicate charge accumulation, and the blue regions represent charge loss.

Figure 6. ELF isosurface and internal cross-section along plane [Forumla omitted. See PDF.] of γ-Fe/La2O2S interface: (a) ELF isosurface of the La(S)-terminated La2O2S/γ-Fe interface; (b) ELF isosurface of S-terminated La2O2S/γ-Fe interface; (c) ELF cross-section of S-terminated La2O2S/γ-Fe interface; (d) ELF cross-section of S-terminated La2O2S/γ-Fe interface (isosurface level = 0.1 eV/Bohr3). The numbers marked in the section denote the Bader charges (unit: e) of interfacial atoms. Positive and negative denote gain and loss of electrons respectively. Coordinate dimension a, b and c represent x-axis, y-axis and z-axis respectively.

Figure 7. The partial density of states different atomic layer near the interface: (a) the La(S)-terminated La2O2S/γ- Fe interface; (b) S-terminated La2O2S/γ-Fe interface (the vertical dashed line indicates the position for Fermi level).

Table 1

Crystal structure, space group, lattice parameters and atomic positions of γ-Fe and La₂O₂S.

Phase	Crystal Structure	Space Group	Lattice Parameters (nm)	Atomic Positions
				Atom	Mult.	Wyck.	x	y	z
γ-Fe	FCC	$\mathrm{Fm}\overline{3}m$	a = 0.3654	Fe	4	a	0	0	0
La2O2S	Hexagonal	$\mathrm{P}\overline{3}m1$	a = 0.4035c = 0.6914	La	2	d	0.3333	0.6667	0.279
				O	2	d	0.3333	0.6667	0.629
				S	1	a	0	0	0

Table 2

CP rows and CP planes of γ-Fe and La₂O₂S.

Phase	CP Rows	Interatomic Spacing (nm)	Type	CP Planes	Interplanes Spacing (nm)
γ-Fe	$⟨2\overline{1}1⟩$	0.224	ZZ	$\left\{111\right\}$	0.211
	$⟨10\overline{1}⟩$	0.258	SS	$\left\{020\right\}$	0.183
	$⟨100⟩$	0.365	SS	$\left\{022\right\}$	0.129
				$\left\{131\right\}$	0.110
La2O2S	$⟨0001⟩$	0.229	PS	$\left\{01\overline{1}1\right\}$	0.312
	$⟨\overline{1}101⟩$	0.325	PS	$\left\{01\overline{1}2\right\}$	0.246
	$⟨2\overline{1}\overline{1}0⟩$	0.398	SS	$\left\{0003\right\}$	0.230
				$\left\{11\overline{2}0\right\}$	0.202

Table 3

Interatomic spacing misfit (f_r) along CP or nearly CP matching rows and interplanar spacing mismatch (f_d) of CP or nearly CP matching planes between γ-Fe and La₂O₂S.

Matching Row Pairs	Misfit (fr) (%)	Matching Type	Matching Plane Pairs	Mismatch (fd)(%)
${⟨2\overline{1}\overline{1}0⟩}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {⟨100⟩}_{\mathsf{\gamma }-\mathrm{Fe}}$	9.43	SS-SS	${\left\{11\overline{2}0\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{111\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	4.58
			${\left\{0003\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{111\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	8.45
			${\left\{11\overline{2}0\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{020\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	9.43
			${\left\{01\overline{1}2\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{111\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	14.15
			${\left\{0003\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{020\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	20.72
			${\left\{01\overline{1}2\right\}}_{{\mathrm{La}}_2{\mathrm{O}}_2\mathrm{S}}$ $\parallel {\left\{020\right\}}_{\mathsf{\gamma }-\mathrm{Fe}}$	25.65

Note: ∥— Matching.

Table 4

Calculated lattice constants, bulk module and equilibrium volume of the bulk γ-Fe compared with other references and experimental data.

Phase		Methods	a (nm)	B0 (GPa)	V0 (Å3/Cell)
γ-Fe	This work	LDA	0.3372	344	38.3
		GGA-PBE	0.3636	152	48.1
	Others	LDA [39]	0.3382	320	38.7
		GGA-PBE [25]	0.344	-	40.9
		GGA-PW91 [39]	0.3472	240	41.9
	Exp. [31]		0.3654	-	48.8

Table 5

Calculated lattice constants, bulk module and equilibrium volume per atom of the bulk La2O2S compared with other references and experimental data.

Phase		Methods	a (nm)	c (nm)	c/a	B0 (GPa)	V0 (Å3/Cell)
La2O2S	This work	LDA	0.4065	0.6939	1.707	111	99.32
		GGA-PBE	0.4134	0.7073	1.711	102	104.67
	Others	LDA [40]	0.3982	0.6826	1.714	128	93.73
		GGA-PBE [40]	0.4055	0.6944	1.712	115	98.91
		GGA-PBE [41]	0.406	0.695	1.712	-	-
	Exp. [32]		0.4035	0.6914	1.713	-	97.49

Table 6

Surface energy of γ-Fe$\left(00\overline{2}\right)$.

Layer	Surface Energy (J/m2)
	5	7	9	11
γ-Fe$\left(00\overline{2}\right)$	2.18	2.18	2.19	2.20

Table 7

Change of the interlayer distance of the La-terminated surface as a percentage of the respective bulk spacing.

Slab Thickness (n)
d-Layer	La(O)-Terminated Surface				La(S)-Terminated Surface
	9	14	19	24	8	13	18	23
d1-2	11.677	10.600	12.295	12.164	4.500	5.663	6.178	5.534
d2-3	−2.300	−1.508	−2.361	−2.443	−2.387	−2.121	−2.748	−1.732
d3-4		−0.475	−0.860	−0.774		0.184	−0.084	0.331
d4-5		−0.137	0.370	0.077		−0.773	−0.208	−0.788
d5-6			−0.454	−0.507			0.356	−0.323
d6-7				1.534			−0.079	−0.866
d7-8				0.466				−0.620

Table 8

Change of the interlayer distance of the O-terminated surface as a percentage of the respective bulk spacing.

Slab Thickness (n)
d-Layer	O(La)-Terminated Surface				O(O)-Terminated Surface
	5	10	15	20	7	12	17	22
d1-2	8.076	8.532	8.116	8.426	−7.891	−9.082	−8.793	−9.164
d2-3	−0.899	−0.386	−0.410	−0.306	17.474	19.441	18.994	19.533
d3-4		−0.880	−0.576	−0.497		−1.594	−1.289	−1.489
d4-5			−0.221	−0.595		−0.192	−0.178	0.049
d5-6			−0.199	−0.039			−0.882	−0.280
d6-7				−0.043			−0.256	−0.286
d7-8				−0.010				0.424
d8-9								−0.004

Table 9

Change of the interlayer distance of the S-terminated surface as a percentage of the respective bulk spacing.

Slab Thickness(n)
d-Layer	S-Terminated Surface
	6	11	16	21
d1-2	0.243	−1.037	−1.269	−0.929
d2-3	−1.093	−2.080	−1.917	−2.001
d3-4		−0.060	0.112	−0.089
d4-5			0.268	0.335
d5-6			−0.039	−0.768
d6-7				−0.466

Table 10

Interface ideal adhesion energy $W_{\mathrm{ad}}$ for the five γ-Fe/La₂O₂S interface systems.

Interface	${\mathit{W}}_{\mathbf{ad}} (\mathbf{J}/{\mathbf{m}}^2)$
La(O)-terminated	2.887
La(S)-terminated	3.372
O(La)-terminated	2.132
O(O)-terminated	2.584
S-terminated	4.403

References

1. Bartlett, L.N.; Avila, B.R. Grain Refinement in Lightweight Advanced High-Strength Steel Castings. Int. J. Metalcast.; 2016; 10, pp. 401-420. [DOI: https://dx.doi.org/10.1007/s40962-016-0048-0]

2. Li, G.; Lan, P.; Zhang, J.Q.; Wu, G.X. Refinement of the Solidification Structure of Austenitic Fe-Mn-C-Al TWIP Steel. Metall. Mater. Trans. B; 2020; 51, pp. 452-466. [DOI: https://dx.doi.org/10.1007/s11663-020-01773-4]

3. Suito, H.; Ohta, H.; Morioka, S. Refinement of Solidification Microstructure and Austenite Grain by Fine Inclusion Particles. ISIJ Int.; 2006; 46, pp. 840-846. [DOI: https://dx.doi.org/10.2355/isijinternational.46.840]

4. Ji, Y.P.; Li, Y.M.; Zhang, M.X.; Qu, W.; Zhao, T.X.; Ren, H.P. Grain Refinement Mechanism of the δ-Ferrite in Steels Through Cerium Addition. Metall. Mater. Trans. A; 2020; 51, pp. 1707-1718. [DOI: https://dx.doi.org/10.1007/s11661-020-05645-z]

5. Wurlitzer, M.; Tonn, B. Influence of cerium on the microstructure and mechanical properties of an as-cast and Q&T steel. Int. J. Cast Met. Res.; 2019; 32, pp. 171-180.

6. Luo, S.; Shen, Z.S.; Yu, Z.M.; Wang, W.L.; Zhu, M.Y. Effect of Ce Addition on Inclusions and Grain Structure in Gear Steel 20CrNiMo. Steel Res. Int.; 2020; 92, 2000394. [DOI: https://dx.doi.org/10.1002/srin.202000394]

7. Wang, L.M.; Lin, Q.; Yue, L.J.; Liu, L.; Guo, F.; Wang, F.M. Study of application of rare earth elements in advanced low alloy steels. J. Alloys Compd.; 2008; 451, pp. 534-537. [DOI: https://dx.doi.org/10.1016/j.jallcom.2007.04.234]

8. Gao, J.Z.; Fu, P.X.; Liu, H.W.; Li, D.Z. Effects of Rare Earth on the Microstructure and Impact Toughness of H13 Steel. Metals; 2015; 5, pp. 383-394. [DOI: https://dx.doi.org/10.3390/met5010383]

9. Guo, M.X.; Suito, H. Influence of Dissolved Cerium and Primary Inclusion Particles of Ce₂O₃ and CeS on Solidification Behavior of Fe-0.20 mass%C–0.02 mass%P Alloy. ISIJ Int.; 1999; 39, pp. 722-729. [DOI: https://dx.doi.org/10.2355/isijinternational.39.722]

10. Li, H.; McLean, A.; Rutter, J.W.; Sommerville, I.D. Influence of rare earth metals on the nucleation and solidification behavior of iron and 1045 steel. Metall. Trans. B; 1988; 19, pp. 383-395. [DOI: https://dx.doi.org/10.1007/BF02657736]

11. Bramfitt, B.L. The effect of carbide and nitride additions on the heterogeneous nucleation behavior of liquid iron. Metall. Trans.; 1970; 1, pp. 1987-1995. [DOI: https://dx.doi.org/10.1007/BF02642799]

12. Chen, L.H.; Li, C.R.; Li, Z.Y. Effect of rare earth Ce on martensite structure of Fe-Mn damping alloy. Mater. Res. Express; 2020; 6, 1265i4. [DOI: https://dx.doi.org/10.1088/2053-1591/ab6d2d]

13. Gao, S.; Wang, M.; Guo, J.L.; Wang, W.; Zhi, J.G.; Bao, Y.P. Characterization Transformation of Inclusions Using Rare Earth Ce Treatment on Al-Killed Titanium Alloyed Interstitial Free Steel. Steel Res. Int.; 2019; 90, 1900194. [DOI: https://dx.doi.org/10.1002/srin.201900194]

14. Ji, Y.P.; Zhang, M.X.; Ren, H.P. Roles of Lanthanum and Cerium in Grain Refinement of Steels during Solidification. Metals; 2018; 8, 884. [DOI: https://dx.doi.org/10.3390/met8110884]

15. Yu, Y.C.; Zhang, S.H.; Li, H.; Wang, S.B. Effect of Rare Earth Lanthanum on the Solidification Sturcture and Hot Ductility of Fe-43Ni Expansion Alloy. High Temp. Mater. Proc.; 2018; 37, pp. 261-269. [DOI: https://dx.doi.org/10.1515/htmp-2016-0157]

16. Tuttle, R.; Kapadia, H.A. Thermal Analysis of Rare Earth Grain Refined 4130. Int. J. Metalcast.; 2019; 13, pp. 273-285. [DOI: https://dx.doi.org/10.1007/s40962-018-0274-8]

17. Tuttle, R. Study of Rare Earth Additions Effect on the Solidification and Properties of 4130 Steel. J. Mater. Eng. Perform.; 2019; 28, pp. 6720-6727. [DOI: https://dx.doi.org/10.1007/s11665-019-04410-7]

18. Tuttle, R. Effect of Rare Earth Master Alloys on 4130. Int. J. Metalcast.; 2022; 16, pp. 35-44. [DOI: https://dx.doi.org/10.1007/s40962-021-00584-w]

19. Zhang, M.X.; Kelly, P.M. Edge-to-edge matching model for predicting orientation relationships and habit planes—the improvements. Scr. Mater.; 2005; 52, pp. 963-968. [DOI: https://dx.doi.org/10.1016/j.scriptamat.2005.01.040]

20. Zhang, M.X.; Kelly, P.M.; Easton, M.A.; Taylor, J.A. Crystallographic study of grain refinement in aluminum alloys using the edge-to-edge matching model. Acta Mater.; 2005; 53, pp. 1427-1438. [DOI: https://dx.doi.org/10.1016/j.actamat.2004.11.037]

21. Qiu, D.; Zhang, M.X.; Kelly, P.M. Crystallography of the heterogeneous nucleation of Mg grains on Al₂Y nucleation particles in an Mg-10 wt.% Y alloy. Scr. Mater.; 2009; 61, pp. 312-315. [DOI: https://dx.doi.org/10.1016/j.scriptamat.2009.04.011]

22. Liu, Z.L.; Qiu, D.; Wang, F.; Taylor, J.A.; Zhang, M.X. The grain refining mechanism of cast zinc through silver inoculation. Acta Mater.; 2014; 79, pp. 315-326. [DOI: https://dx.doi.org/10.1016/j.actamat.2014.07.026]

23. Wang, J.H.; Li, Y.S.; Xu, R. First-principles calculations on electronic structure and interfacial stability of Mg/NbB₂ heterogeneous nucleation interface. Surf. Sci.; 2020; 691, 121487. [DOI: https://dx.doi.org/10.1016/j.susc.2019.121487]

24. Li, J.; Cui, Y.M.; Zhang, M.; Zhao, J.Y.; Luo, X. Thermodynamic evidence of α-Al heterogeneous nucleation on Al₂MgC₂ and the interfacial bonding mechanism: A first-principles study. J. Solid State Chem.; 2020; 288, 121431. [DOI: https://dx.doi.org/10.1016/j.jssc.2020.121431]

25. Jiao, X.Y.; Fu, W.T.; Shao, W.; Zhu, X.W.; Zhou, Y.F.; Xing, X.L.; Wang, Z.R.; Yang, Q.X. First-principles calculation on γ-Fe/La₂O₃ interface properties and austenite refinement mechanism by La₂O₃. Mater. Chem. Phys.; 2021; 259, 124194. [DOI: https://dx.doi.org/10.1016/j.matchemphys.2020.124194]

26. Liu, R.; Yin, X.M.; Feng, K.X.; Xu, R. First-principles calculations on Mg/TiB₂ interfaces. Comput. Mater. Sci.; 2018; 149, pp. 373-378. [DOI: https://dx.doi.org/10.1016/j.commatsci.2018.03.045]

27. Yang, J.; Hou, X.R.; Zhang, P.F.; Zhou, Y.F.; Xing, X.L.; Ren, X.J.; Yang, Q.X. First-principles calculations on LaAlO₃ as the heterogeneous nucleus of austenite. Comput. Theor. Chem.; 2014; 1029, pp. 48-56. [DOI: https://dx.doi.org/10.1016/j.comptc.2013.12.023]

28. Li, R.W.; Chen, Q.C.; Ouyang, L.; Ding, Y.L. Adhesion strength and bonding mechanism of γ-Fe (111)/α-Al₂O₃ (0001) interfaces with different terminations. J. Alloys Compd.; 2021; 870, 159529. [DOI: https://dx.doi.org/10.1016/j.jallcom.2021.159529]

29. Yang, J.; Zhang, P.F.; Zhou, Y.F.; Guo, J.; Ren, X.J.; Yang, Y.L.; Yang, Q.X. First-principles study on ferrite/TiC heterogeneous nucleation interface. J. Alloys Compd.; 2013; 556, pp. 160-166. [DOI: https://dx.doi.org/10.1016/j.jallcom.2012.12.099]

30. Yuan, X.H.; Xiao, Y.Y.; Wang, G.C.; Zhang, L.F. TiN inducing ferrite nucleation based on the bcc-Fe/TiN interfaces formation at atomic scale by first-principles calculation. Comput. Mater. Sci.; 2021; 197, 110570. [DOI: https://dx.doi.org/10.1016/j.commatsci.2021.110570]

31. Straumanis, M.E.; Kim, D.C. Lattice Constants, Thermal Expansion Coefficients, Densities and Perfection of Structure of Pure Iron and of Iron Loaded with Hydrogen. Int. J. Mater. Res.; 1965; 233, pp. 1519-1525. [DOI: https://dx.doi.org/10.1515/ijmr-1969-600404]

32. Mikami, M.; Nakamura, S. Electronic structure of rare-earth sesquioxides and oxysulfides. J. Alloys Compd.; 2006; 408–412, pp. 687-692. [DOI: https://dx.doi.org/10.1016/j.jallcom.2005.01.068]

33. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B; 1996; 54, 11169. [DOI: https://dx.doi.org/10.1103/PhysRevB.54.11169] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9984901]

34. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B; 1999; 59, 1758. [DOI: https://dx.doi.org/10.1103/PhysRevB.59.1758]

35. Perdew, J.P.; Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B; 1981; 23, 5048. [DOI: https://dx.doi.org/10.1103/PhysRevB.23.5048]

36. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.; 1997; 78, pp. 3865-3868. [DOI: https://dx.doi.org/10.1103/PhysRevLett.77.3865] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10062328]

37. Yang, J.; Huang, J.H.; Fan, D.Y.; Chen, S.H.; Zhao, X.K. Comparative investigation on RE(La,Ce)AlO₃(100)/γ-Fe(100) interfaces: A first-principles calculation. Appl. Surf. Sci.; 2016; 384, pp. 207-216. [DOI: https://dx.doi.org/10.1016/j.apsusc.2016.05.024]

38. Wang, V.; Xu, N.; Liu, J.C.; Tang, G.; Geng, W.T. VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code. Comput. Phys. Commun.; 2021; 267, 108033. [DOI: https://dx.doi.org/10.1016/j.cpc.2021.108033]

39. Yu, J.; Lin, X.; Wang, J.J.; Chen, J.; Huang, W.D. First-principles study of the relaxation and energy of bcc-Fe, fcc-Fe and AISI-304 stainless steel surfaces. Appl. Surf. Sci.; 2009; 255, pp. 9032-9039. [DOI: https://dx.doi.org/10.1016/j.apsusc.2009.06.087]

40. Alemany, P.; Moreira, I.P.R.; Castillo, R.; Llanos, J. Electronic, structural, and optical properties of host materials for inorganic phosphors. J. Alloys Compd.; 2012; 513, pp. 630-640. [DOI: https://dx.doi.org/10.1016/j.jallcom.2011.11.036]

41. Polfus, J.M.; Norby, T.; Bredesen, R. Protons in Oxysulfides, Oxysulfates, and Sulfides: A First-Principles Study of La₂O₂S, La₂O₂SO₄, SrZrS₃, and BaZrS₃. J. Phys. Chem. C; 2015; 119, pp. 23875-23882. [DOI: https://dx.doi.org/10.1021/acs.jpcc.5b08278]

42. Reuter, K.; Scheffler, M. Composition, structure, and stability of RuO₂(110) as a function of oxygen pressure. Phys. Rev. B; 2007; 78, 049901. [DOI: https://dx.doi.org/10.1103/PhysRevB.75.049901]

43. Wei, C.; Ren, Q.Q.; Fan, J.L.; Gong, H.R. Cohesion properties of W/La₂O₃ interfaces from first principles calculation. J. Nucl. Mater.; 2015; 466, pp. 234-238. [DOI: https://dx.doi.org/10.1016/j.jnucmat.2015.08.013]

44. Zhao, X.B.; Yuan, X.M.; Liu, S.; Zhao, C.C.; Wang, C.X.; Zhou, Y.F.; Yang, Q.X. Investigation on WC/LaAlO₃ heterogeneous nucleation interface by first-principles. J. Alloys Compd.; 2017; 695, pp. 1753-1762. [DOI: https://dx.doi.org/10.1016/j.jallcom.2016.11.005]

45. Becke, A.D.; Edgecombe, K.E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys.; 1990; 92, 5397. [DOI: https://dx.doi.org/10.1063/1.458517]

46. Zhang, W.; Smith, J.R.; Evans, A.G. The connection between ab initio calculations and interface adhesion measurements on metal/oxide systems: Ni/Al₂O₃ and Cu/Al₂O₃. Acta Mater.; 2002; 50, pp. 3803-3816. [DOI: https://dx.doi.org/10.1016/S1359-6454(02)00177-5]

47. Han, Y.F.; Bai, Y.B.; Wang, J.; Shu, D.; Sun, B.D. First-principles calculations on Al/AlB₂ interfaces. Appl. Surf. Sci.; 2011; 257, pp. 7831-7836. [DOI: https://dx.doi.org/10.1016/j.apsusc.2011.04.038]

48. Jiang, Q.; Lu, H.M. Size dependent interface energy and its applications. Surf. Sci. Rep.; 2008; 63, pp. 427-464. [DOI: https://dx.doi.org/10.1016/j.surfrep.2008.07.001]