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

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 La2O2S and γ-Fe from a crystallographic point of view. A rough orientation relationship (OR) was predicted with the minimum values of fr=9.43% and fd=20.72% as follows: [21¯1¯0]La2O2S∥[100]γ-Fe and (0003¯)La2O2S∥(002¯)γ-Fe. The interface energy and bonding characteristics between La2O2S 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 La2O2S/γ-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, La2O2S particles are not an effective heterogeneous nucleation site for the γ-Fe matrix from crystallography and energetics points of view.


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 2 O 2 S and γ-austenite in an Fe-43Ni expansion alloy and reported the lowest δ-value of 5.42% between the (0001) crystal plane of La 2 O 2 S and the (100) crystal plane of γ-austenite. This implied a high heterogeneous nucleation potency of La 2 O 2 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 3 [27], La 2 O 3 [25], Al 2 O 3 [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 2 O 2 S, which is the typical oxysulfide, was selected. The present work aims to clarify the heterogeneous nucleation potency of γ-Fe on La 2 O 2 S particles, combining a crystallographic calculation using the edge-to-edge matching model with the first-principles calculations from an energetics point of view.

Crystallographic Calculations Using the E2EM Model
To investigate the heterogeneous nucleation potency of La 2 O 2 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 2 O 2 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 2 O 2 S.
The γ-Fe phase possesses a face-centered cubic (FCC) lattice structure with four atoms in each unit cell. By contrast, the La 2 O 2 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 = Fm3m. The La 2 O 2 S lattice parameters were [32]: a = 0.4035 nm; c = 0.6914 nm; space group = P3m1.

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 2 O 2 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 2 O 2 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 2 O 2 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 −4 eV/atom, and the residual Hellmann-Feynman force on each relaxed atom reached 0.03 eV/Å [27][28][29].

Nucleation Crystallography of γ-Fe on La 2 O 2 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 2 O 2 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 2 O 2 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 2 O 2 S and γ-Fe could be predicted based on the matching directions and matching planes. The predicted rough OR between La 2 O 2 S and γ-Fe could be expressed as 21     In view of the lattice match of the interface, γ-Fe 002 and La 2 O 2 S 0003 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 002 through the method presented in refs. [25,28,37]. The surface energy can be expressed by [30]: where E slab is the total energy of the slab; E bulk is the total energy of bulk after energy optimization; A is the corresponding surface area; and N slab and N 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 002 surface energy converged to 2.18 J/m 2 when the surface model contained more than five atomic layers. This figure is close to the experimental value 1.95 J/m 2 reported in ref. [39] Therefore, the five-layer γ-Fe 002 surface slab was adopted to ensure the bulk-like interior in the following calculations. For the La 2 O 2 S 0003 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 Laterminated 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. Tables 7-9 show the La 2 O 2 S 0003 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 Tables 7-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.  Tables 7-9 show the La2O2S 0003 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 Tables 7-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.    [8][9] −0.004  (2) [23,[25][26][27][28]42]: where E slab is total energy of the corresponding slab, N i is the number of atomic species i (i = La, O and S) in the La 2 O 2 S slab, µ 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 2 O 2 S slab in equilibrium with those of bulk La 2 O 2 S, the relationship between the chemical potential of La 2 O 2 S bulk µ bulk La 2 O 2 S , its 0 K formation heat (∆H 0 f (La 2 O 2 S)) and chemical potentials of the La bulk, O 2 gas and S bulk (µ bulk La , 1 2 µ gas O 2 and µ bulk S ) are expressed by [43]: where µ 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 µ i (i = La, O and S) at the surface is the corresponding potential in pure La bulk, gas phase of O 2 and pure S bulk: , which is the basis for forming stable slab structure. The range of µ La − µ bulk La was obtained: Based on Equations (3), (5) and (7), the surface energy (Equation (2)) can be rewritten as: where µ bulk La and µ bulk S refers to the single atomic energy in the La and S bulk; µ bulk La 2 O 2 S is equal to the energy of La 2 O 2 S bulk per formula unit (f.u.); and N i (i = La, O and S) has the same meaning as in Equation (2). For La 2 O 2 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 (∆µ S = 0). The range of the ∆µ La value combined with the formula available for La chemical potential (Equations (4)-(6)) can be derived as: The formation heat ∆H 0 f (La 2 O 2 S) was calculated to be −5.12 eV. The range of ∆µ La was, therefore, derived to be: Figure 2 shows La 2 O 2 S 0003 surface energy varying with the variation of the chemical potential of ∆µ La . Clearly, the surface energy of these surfaces was linearly correlated with ∆µ 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 ∆µ La in the range of −0.74 J/m 2~− 6.29 J/m 2 . By contrast, the surface energy of La(O)-terminated, O(La)-terminated and S-terminated surfaces remained constant at 4.79 J/m 2 , 0.68 J/m 2 and 0.06 J/m 2 , respectively. It should be pointed out that the constant value was due to the zero value of (N O − N La ) in Equation (10). Therefore, the La(S)-terminated surface was the most stable one over the entire range of chemical potential.

Structures of La2O2S/γ-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 La2O2S and γ-Fe were constructed. They were denoted by La(O)-terminated, La(S)-terminated, O(La)-terminated, O(O)-terminated and S-terminated La2O2S/γ-Fe interfaces, as shown in Figure 3. A supercell geometry where a five-layer γ-Fe 002 slab was placed under a corresponding, symmetric atomic terminals La2O2S 0003 slab was used to simulate the La2O2S/γ-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 21 1 0 La 2 O 2 S ∥ 100 γ-Fe and 0003 La 2 O 2 S ∥ 002 γ-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 La2O2S and γ-Fe was used. The interface structures were fully re-  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 2 O 2 S and γ-Fe was used. The interface structures were fully relaxed in all dimensions by minimizing the Hellmann-Feynman force on every atom.

Interface Stability
Binding strength of the interface can be predicted by the work of adhesion ( 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 ad is, the stronger the binding force of interface atoms is. The ad of the interface can be calculated by Equation (11) where La 2 O 2 S/γ-Fe is the total energy of interface system; La 2 O 2 S and γ-Fe are the total energy of isolated La2O2S and γ-Fe slabs in the same interface system when one of the slabs is kept and other is replaced with vacuum, respectively; and is the interface area. Generally, a more positive ad means a larger binding strength. Table 10 presents the optimal ad for the relaxed geometries of the five La2O2S/γ-Fe interfaces. The ad of the five La2O2S/γ-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.

Interface Stability
Binding strength of the interface can be predicted by the work of adhesion (W 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 ad is, the stronger the binding force of interface atoms is. The W ad of the interface can be calculated by Equation (11) [23][24][25][26][27][28]30]: where E La 2 O 2 S/γ−Fe is the total energy of interface system; E La 2 O 2 S and E γ−Fe are the total energy of isolated La 2 O 2 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 ad means a larger binding strength. Table 10 presents the optimal W ad for the relaxed geometries of the five La 2 O 2 S/γ-Fe interfaces. The W ad of the five La 2 O 2 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 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 2 O 2 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 2 O 2 S on the γ-Fe substrate. According to the calculation method provided in refs. [23,26,30], the interface energy can be derived as follows: where E La 2 O 2 S/γ−Fe is the total energy of interface system; µ bulk Fe is equal to the energy of γ-Fe bulk per formula unit (f.u.); δ La 2 O 2 S and δ γ−Fe denote the surface energy of La 2 O 2 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 2 O 2 S/γ-Fe interfaces with the increase of ∆µ La . The La 2 O 2 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 eV ≤ ∆µ La ≤ −4.69 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 eV ≤ ∆µ La ≤ 0 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. In addition to the work of adhesion ( 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 La2O2S/γ-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 La2O2S on the γ-Fe substrate. According to the calculation method provided in refs. [23,26,30], the interface energy can be derived as follows: where La 2 O 2 S/γ-Fe is the total energy of interface system; Fe bulk is equal to the energy of γ-Fe bulk per formula unit (f.u.); La 2 O 2 S and γ-Fe denote the surface energy of La2O2S 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 La2O2S/γ-Fe interfaces with the increase of ∆ La . The La2O2S/γ-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.12eV ∆ La 4.69eV, 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.69eV ∆ La 0eV, 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.

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

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]: where ρ int is the total charge density of La 2 O 2 S/γ-Fe interface system; and ρ La 2 O 2 S and ρ γ−Fe are the charge density of corresponding isolated La 2 O 2 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. (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]: where int is the total charge density of La2O2S/γ-Fe interface system; and La 2 O 2 S and γ-Fe are the charge density of corresponding isolated La2O2S 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]: where is the real electron gas density; and ℎ 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 3 . 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 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]: where D r is the real electron gas density; and Dh 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 3 . 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.
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 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. 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.

Heterogeneous Nucleation Potency of γ-Fe on La2O2S Particles
To ensure sufficient atomic matching to retain a coherent or semi-coherent interface between the two phases, empirical criteria of the fr and fd 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 La2O2S and γ-Fe, the minimum values of fr and fd in construction of the OR were 9.43% and 20.72%, respectively. This implies that, from a crystallographic point view, La2O2S 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 La2O2S the effective heterogeneous nucleus of γ-Fe during the solidification, the solid-solid interface energy of La2O2S/γ-Fe should be between 0 J/m 2 and 0.204 J/m 2 from the viewpoint of thermodynamics [37], for which the theoretical value of liquid-solid interface energy is 0.241 J/m 2 [48], as shown in Figure 4. As can be seen from Figure 4, the interface energy of La2O2S/γ-Fe could not meet the interface energy requirements as a heterogeneous nucleation interface. It did not comply with the thermodynamic conditions that make La2O2S particles effective as potent heterogeneous nucleation substrates for the γ-Fe matrix. In addition, factors, such as the size and distribution of La2O2S particles, also impact whether La2O2S 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.

Conclusions
1. With the crystallographic calculations using the E2EM model, a rough HCP/FCC orientation relationship (OR) was predicted between La2O2S 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 21 1 0 La 2 O 2 S ∥ 100 γ-Fe and 0003 La 2 O 2 S ∥ 002 γ-Fe . 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 2 O 2 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 2 O 2 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 2 O 2 S the effective heterogeneous nucleus of γ-Fe during the solidification, the solid-solid interface energy of La 2 O 2 S/γ-Fe should be between 0 J/m 2 and 0.204 J/m 2 from the viewpoint of thermodynamics [37], for which the theoretical value of liquid-solid interface energy is 0.241 J/m 2 [48], as shown in Figure 4. As can be seen from Figure 4, the interface energy of La 2 O 2 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 2 O 2 S particles effective as potent heterogeneous nucleation substrates for the γ-Fe matrix. In addition, factors, such as the size and distribution of La 2 O 2 S particles, also impact whether La 2 O 2 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.

1.
With the crystallographic calculations using the E2EM model, a rough HCP/FCC orientation relationship (OR) was predicted between La 2 O 2 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 21 10 La 2 O 2 S [100] γ−Fe and 0003 La 2 O 2 S 002 γ−Fe .

2.
For γ-Fe/La 2 O 2 S 0003 interfaces, the La(S)-terminated and S-terminated interfaces were more stable than others with smaller interface energy. In the range of −5.12 eV ≤ ∆µ La ≤ −4.69 eV, the S-terminated interface structure was the most stable. Nevertheless, for the range of ∆µ 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 ad , which was 4.403 J/m 2 , followed by the La(S)-terminated interface, which was 3.372 J/m 2 .

3.
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. 4.
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 2 O 2 S as a nucleation site for γ-Fe. The first-principles calculations revealed that the interfacial energy of La 2 O 2 S/γ-Fe cannot meet the interface energy requirements as a heterogeneous nucleation interface. Thus, La 2 O 2 S particles are noneffective heterogeneous nucleation substrates for a γ-Fe matrix from crystallography and thermodynamics viewpoints. Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.