Next Article in Journal
Significance or Scatter—Statistical Evaluation of Rapid Chloride Migration Test Results for Sprayable Cement Mortar
Previous Article in Journal
Strength of Partially Encased Steel-Concrete Composite Column for Modular Building Structures
 
 
Correction published on 14 July 2023, see Materials 2023, 16(14), 4997.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Principle Surface Analysis of YF3 and Isostructural HoF3

Institute for Chemistry and Biochemistry, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
*
Author to whom correspondence should be addressed.
Materials 2022, 15(17), 6048; https://doi.org/10.3390/ma15176048
Submission received: 3 August 2022 / Revised: 22 August 2022 / Accepted: 24 August 2022 / Published: 1 September 2022 / Corrected: 14 July 2023
(This article belongs to the Section Materials Chemistry)

Abstract

:
The trifluorides of the two high field strength elements yttrium and holmium are studied by periodic density functional theory. As a lanthanide, holmium also belongs to the group of rare earth elements (REE). Due to their equivalent geochemical behavior, both elements form a geochemical twin pair and consequently, yttrium is generally associated with the REE as REE+Y. Interestingly, it has been found that DFT/DFT+U describe bulk HoF3 best, when the 4f-electrons are excluded from the valence region. An extensive surface stability analysis of YF3 (PBE) and HoF3 (PBE+U d /3 eV/4f-in-core) using two-dimensional surface models (slabs) is performed. All seven low-lying Miller indices surfaces are considered with all possible stoichiometric or substoichiometric terminations with a maximal fluorine-deficit of two. This leads to a scope of 24 terminations per compound. The resulting Wulff plots consists of seven surfaces with 5–26% abundance for YF3 and six surfaces with 6–34% for HoF3. The stoichiometric (010) surface is dominating in both compounds. However, subtle differences have been found between these two geochemical twins.

Graphical Abstract

1. Introduction

Yttrium and holmium form a geochemical twin pair. The term emphasizes their identical geochemical behavior caused by the equal ratio of charge to radius in their only stable oxidation state +III. According to their small ionic radii of 1.075 Å (Y) and 1.072 Å (Ho) in nine-fold coordination [1] and their high oxidation state, both belong to the interesting group of high field strength elements (HFSE). As a lanthanide, holmium also belongs to the rare earth elements (REE). Due to their twin character, yttrium is also often associated with that group [2,3,4].
As fluorides, both metals can be used for different specific applications. The wide-band-gap material YF3 has very good properties for laser applications [5,6,7,8]. Doped with trivalent REE cations, YF3 is also applicable as an optical filter in 157-nm photolithography [9]. Another emerging field of application is solid-state fluoride batteries, resulting from the very high conductivity of fluoride anions [10,11,12,13,14]. HoF3 is interesting for magnetic high-field applications as, e.g., a contrast agent, due to the very high magnetic moment of holmium [15,16]. Moreover, YF3 and HoF3 are important precursors for the synthesis of the respective pure metallic compounds [17,18]. In nature, YF3 is found within the mineral waimirite-(Y), which contains high concentrations of other REE [8]. Fluoride plays a significant role in accumulating HFSE and REE within hydrothermal fluids, as these cations do not form such stable complexes with chloride [2,19,20,21]. Interestingly, those fluoride-rich hydrothermal fluids produce ores with a non-chondritic excess of yttrium over holmium. It is suggested that one underlying reason for the twin separation is their different affinity to fluorine, which was found in dissolving experiments of YF3 and HoF3 in diluted hydrofluoric acid [20]. To lay one foundation for future quantum chemical studies on the different fluorine-affinity of yttrium and holmium, we started with an investigation of the respective trifluorides and their surfaces.
In accordance with their twin behavior, solid YF3 and HoF3 occur in the same crystal structure type of β -YF3 (space group Pnma, fully occupied Wyckoff positions (Y) 4c, (F) 4c, (F) 8d (Figure 1)) [18,22,23]. This is their only stable phase up to 1343 K (HoF3) or 1350 K (YF3), which is well beyond the temperature regime of hydrothermal fluids of typically 323–873 K [7,18,24,25]. The ionic radii of the middle and late lanthanides Sm(III)–Lu(III) differ by only ≤6 pm compared to Y(III) [1]. Consequently, all their trifluorides crystallize as well in β -YF3 [18,22,23,26,27,28,29,30,31]. The same low-temperature phase is also found for the two actinides Bk(III) and Cf(III) [32,33,34]. Due to the often observed analogy of actinides to lanthanides, the same crystal structure is assumed for the, so far, experimentally unknown heavier actinide trifluorides of Es(III)–Lr(III) [35]. Interestingly, the known orthorhombic low-temperature phase of plutonium trihydride is also reasoned to be an exotic example of a β -YF3 structure [36]. Outside the f-block, the β -YF3-structure is experimentally known for bismuth trifluoride [23,37,38] and predicted as an accessible meta-stable phase for the trichlorides of Y(III) and Bi(III) [38].
To the best knowledge of the authors, no first principle surface stability analysis of any compound within the whole structure type is available in the literature. The only surface calculation of any compound of β -YF3 structure was published in 2013 by Ye et al. [39] on two selected surfaces of DyF3 (001) and (101). They only calculated the surfaces matching their experimentally obtained nano-plates. However, these and other experiments on this class of compounds clearly demonstrate that the obtained surface structures are very dependent on the experimental conditions, especially on the utilized nature and geometry of the substrate, as well as on the solvent and fluoride concentration [14,16,40,41,42,43,44,45]. The present work analyzes the inherent quantum chemical stability of all of the seven low Miller indices ( h k l ) surfaces, namely (001), (010), (100), (011), (101), (110) and (111). A previous study on another metal trifluoride, AlF 3 revealed stoichiometric or substoichiometric surfaces with a small fluorine-deficit as the most stable terminations [46]. Additionally, a substoichiometric fluorine content has also been found for YbF3 thin films made from ion assisted deposition [47]. Consequently, this study includes all possible stoichiometric terminations and those with a small-to-moderate fluorine-deficit of 1–2 fluorine atoms per surface unit cell. This results in a scope of 24 terminations. The obtained surface energy results are combined with the geometry of the surface cut by a Wulff analysis to examine the expected surface abundance [48,49].

2. Methodology

2.1. Computational Details

All calculations were performed in the Vienna Ab Initio Simulation Package (VASP, version 5.4.4) [50] on the supercomputer cluster HLRN in Berlin and Göttingen, Germany. using periodic density functional theory (DFT) with a generalised gradient approximation (GGA). As an exchange–correlation functional, the one of Perdew–Burke–Ernzerhof (PBE) is applied [51]. The inner shell electrons were described by the projector augmented wave (PAW) method [52,53]. The outer shell electrons were expanded in plane waves.
For converged YF3 total bulk energies, the VASP potential files F_h (“hard”, 7 electrons) and Y_sv (11 electrons) were applied together with a 9 × 9 × 9 Monkhorst–Pack grid. In accordance with the F_h potential file, a kinetic energy cut-off of 772.6 eV was used. For HoF3, both available Ho potential files Ho_3 (9 electrons, 4f-in-core) and Ho (21 electrons, 4f-in-valence) were evaluated with respective grid sizes of 7 × 7 × 7 and 3 × 3 × 3 . On holmium, the Hubbard-type correction in the simple Dudarev formalism was applied [54]. In a test series of 1–10 eV in 1 eV steps with U d (with Ho_3) and U f (with Ho), PBE+U d with 3 eV agreed best with the crystal structure and the presumed electronic structure (Table 1 and Figures S1 and S2). As an electronic structure reference, bulk HoF3 was also calculated with the Heyd–Scuseria–Ernzerhof hybrid functional (HSE06) [55].
For electron smearing, tests on several bulk and slabs structures of both trifluorides were performed, comparing Gaussian smearing with the tetrahedron method with Blöchl correction [56]. No energy difference within the applied self-consistent field (SCF) convergence criteria could be found. We therefore used Gaussian smearing on our insulating trifluorides.
Apart from the trifluorides, molecular fluorine, as well as metallic yttrium and holmium, were also considered. The first was calculated in a cubic box of 25 Å length. For the latter two, Gaussian smearing could not be applied. A convergence test with 1st and 2nd-order Methfessel–Paxton smearing with widths of 0.05–0.35 eV yielded 2nd-order Methfessel–Paxton smearing with a width of 0.10 eV (Y) or 0.15 eV (Ho) as the best combination to minimize the difference between total energy and free energy.
Each bulk structure started from the respective, experimental crystal structure (YF3 [23], HoF3 [22], Y [57], Ho [58]) and was fully relaxed in atomic positions, lattice constants and volume. The accurate precision setting was applied. As convergence criteria, 0.01 meV per unit cell was used for SCF total energies and 0.1 meV per unit cell for the difference in total energy between two ionic steps. Final total energies, density of states (DOS) and Bader charges were performed with an SCF criteria of 0.001 meV. All DOS plots and Bader charges, as well as all HoF3 data, were calculated with allowed spin polarization. To aid SCF convergence, an additional support grid (.ADDGRID.) and/or a reduced minimal mixing parameter for Kerker’s initial approximation [59] (AMIN) of <0.01 were applied on most slabs.
Symmetric slabs were built from the relaxed bulk structure with the Python package pymatgen [60,61]. The vacuum height perpendicular to the surface was tested for one stoichiometric termination of YF3 (001). The converged value of 25 Å was applied for all slabs. For slab calculations, only one k-point was used perpendicular to the surface. For the other two directions, we applied the same k-point grid size as in bulk. The complete slabs were relaxed in atomic positions.
DOS plots and band structures were generated with pymatgen. Wulff plots were constructed with the WulffPack Python package [62]. Atomic structures were visualized in VESTA [63].

2.2. Choice of Electronic Structure Method

The effect of dispersion was tested by applying Grimme’s dispersion correction with Becke–Johnson damping (D3(BJ)) [64]. From PBE to PBE+D3(BJ), the lattice constants changed only by 1.9–4.5 pm or 0.3–1.0% during the full optimization of atomic positions, lattice constants and volume of YF3. Due to this small deviation, we neglected dispersion correction for our highly ionic systems.
For HoF3, a test series was performed to decide whether to treat the 4f-electrons inside the core or at the valence level. Hubbard-type Coulomb parameters of 1–10 eV were scanned for the 4f-in-core with U d acting on Ho-d orbitals, as well as for 4f-in-valence with U f acting on Ho-f orbitals. It should be noted that the Ho-5d orbitals mainly constituted the broad conduction band in both approaches. Yet, they also hybridized in the valence band mainly constructed by F-2p (Figure S3). The PBE+U benchmark plots for unit cell parameters and band gaps are given in the SI with further discussion (Figures S1 and S2). All HoF3 (PBE+U d /4f-in-core) band structures resembled the YF3 (PBE) one and produced comparable F-2p to Ho-5d or Y-4d charge transfer band gaps of 7–8 eV (Figure S3). By adding exact exchange via HSE06/4f-in-core, these bands were further separated to 11 eV. Whereas, HSE06/4f-in-valence predicted an Ho-4f to Ho-4f transition of 8 eV. In contrast, PBE/4f-in-valence was not able to separate the partially filled 4f 10 into un-/occupied bands. Instead, it placed the Fermi-level ( E F ) inside the 4f band, predicting a pseudo-metal. When introducing the additional Coulomb potential of 1–10 eV onto the 4f in PBE+U f , this 4f–4f gap was tuneable from 1 eV to a maximum of 6 eV. At U f 5 eV, the nature of the band gap changed to a charge transfer of F-2p to Ho-4f. Unfortunately, no measured band gap exists in the literature for HoF3. Therefore, it was not possible to pin-point the true band gap, nor to evaluate the correct nature of that transition. Nevertheless, based on a purely empirical model derived from other lanthanide compounds, HoF3 is expected to have a band gap of ca. 9 eV [65]. This empirically estimated band gap, as well as the calculated HSE06 reference, were best reproduced without including the 4f-electrons explicitly.
Another quantity upon which to judge the applied electronic structure method was the Bader charges obtained by applying the atoms in molecules (AIM) population analysis [66,67,68,69,70]. For both bulk materials of YF3 and HoF3, all tested methods predicted a metal charge of 2.2–2.4 e and fluorine charge of −(0.7–0.8) e. The Bader charges of all applied methods with 4f-in-core or valence agreed well with each other and thus suggested that including 4f explicitly was not necessary for HoF3.
Furthermore, all methods used with 4f-in-valence predicted a high-spin bulk unit cell with all four holmium aligned resulting in an electronic magnetic moment of 16 µ B . This ferromagnetic result was obtained even when starting from anti-ferromagnetic spin arrangements. According to the experimentally known magnetic structures, the physically correct spin arrangement is anti-ferromagnetic below 0.53 K or paramagnetic above [71].
To summarize, not including the 4f-electrons explicitly provided the best electronic structure results. The differences between simple PBE and PBE+U d were minor. When considering the unit cell parameters given in Table 1, PBE+U d /3 eV/4f-in-core performed best with deviations of as little as 0.1–0.8%.

2.3. Surface Energy

The surface formation energy ( E surf ) is generally calculated from the total energy of the 2D-periodic slab ( E n ), the energy of the 3D-periodic bulk unit cell ( E bulk ) and the surface area of the slab (A):
E surf bd = E n n E bulk 2 A .
n is the slab thickness measured in unit cells. We label this bulk-derived surface energy E surf bd . Equation (1) is used for all YF3 surface energies. In this work, we also considered surfaces with a substoichiometric amount of fluorine. For these, the fluorine potential µ F for each missing fluorine was added to the numerator of Equation (1). µ F was obtained from E bulk , the bulk energy per atom of the pure metal of yttrium or holmium (µ M ), as well as the number of metal ( n M = 4 ) and fluorine ( n F = 12 ) atoms within the bulk MF3 unit cell:
μ F = E bulk n M μ M n F .
Yet, as pointed out by Boettger, this bulk-derived surface energy ( E surf bd ) can lead to diverging E surf with respect to n [72]. This can be avoided by using slab-derived (sd) energies only:
E surf sd = E n n ( E n E n 1 ) 2 A .
E bulk is then replaced by the difference of E n to the total energy of the next smaller slab ( E n 1 ). For HoF3, we indeed observed linearly diverging E surf bd when applying Equation (1), despite system sizes of up to 7 UC or Ho28F84. Depending on the ( h k l ), this stoichiometry corresponds to 12, 24 or 26 HoF3-layers. Likely, this is a result of the allowed spin-polarization with Hubbard-type correction and atomic relaxation of the whole slab. It can be seen in Table S3, that this linear divergence only appears after relaxation in E surf , opt bd . The unrelaxed surface energies E surf , SP bd show no divergence. In YF3, no divergent E surf , opt bd are observed. Here, no Hubbard-type correction is applied and the atomic relaxation is performed without spin polarization. A comparison of slab convergence by both equations is given in Tables S2 and S3. Due to the divergence issue, all HoF3 surface energies given within this paper are slab-derived using Equation (3), which nicely converge. As each E surf sd is derived from two slabs differing by one unit cell in size, at least three slabs are needed to determine convergence. Whereas for E surf bd , these are just two. Due to the observed convergence of E surf bd in YF3, only two slab thicknesses are modeled for many terminations. Therefore, the convergence of the respective E surf sd cannot be evaluated. As a consequence, we used the converged E surf bd for YF3 to compare with the converged HoF3  E surf sd . All YF3 bulk-derived surface energies converged within 0.03 J m 2 at slab thickness of about 5–5.5 UC or 10–22 YF3-layers (Table S2). The HoF3 slab-derived surface energies of 14 terminations, including all of the most stable ones per Miller indices, converged to 0.01 J m 2 or less within a slab thickness of about 6–6.5 UC or 12–26 HoF3-layers (Table S3). Some of the higher energy terminations converged only to 0.02–0.04 J m 2 at that thickness, while four high energy terminations did not converge even to 0.1 J m 2 . Fortunately, it is clear from their surface energies that even within the present uncertainty, those high energy terminations do not compete with the lowest energy ones. The slab thickness convergence for HoF3 is visualized by error bars in Figure S5.

3. Results and Discussion

3.1. Bulk Properties

For YF3, the PBE relaxed lattice constants, given in Table 1, agree very well with both experimental values, which are underestimated by as little as 0.5–1.5% [22,23]. The resulting unit cell volume is underestimated by 2.6–2.7%, which is still in good agreement for a GGA functional. The best performing HoF3 method against the only available experimental unit cell data and the calculated HSE06 band gaps was found to be PBE+U d /3 eV/4f-in-core. The resulting unit cell parameters deviate by as little as 0.1–0.8%.
The respective F–M bond length on the PBE (YF3) and PBE+U d level (HoF3) are R F - Y = { 2.28; 2.29; 2.46} Å and R F - Ho = { 2.30; 2.32; 2.45} Å. These agree perfectly with the measured interatomic distances of 2.3–2.6 Å [22,23].
Before we come to the surfaces, we evaluate possible energetic differences between the two geochemical twins as bulk materials. We calculate the electronic contribution to the formation enthalpies ( Δ H f ) according to:
2 M ( s ) + 3 F 2 ( g ) 2 Δ H f 2 MF 3 ( s )
The electronic energies are taken from the bulk metals in hcp (P63/mmc) structure and the bulk trifluorides, as well as molecular fluorine. For YF3, we obtained an electronic contribution of 1591.1 kJ mol 1 versus 1587.3 kJ mol 1 for HoF3. Thus, judged by the electronic energies only, both trifluorides are equally strong bound with a very small favor of 3.7 kJ mol 1 or 0.2% for YF3 over HoF3.

3.2. Surface Energies

The surface energies of all calculated terminations are given in Table 2. The given metal surface coordination number (CN surf ) is determined with a bond length cut-off of R F - M 2.60 Å. Table 2 also includes the nominal net surface charge ( q surf ) caused by substoichiometric fluoride. Finally, the last column includes the surface abundance for each respective most stable termination predicted by Wulff construction (% surf ).
The two terminations, (110)-1 and -2 greatly illustrate the importance of atomic relaxation of the surface prior analysis. Before relaxation, nothing but the very surface layer differs within each ( h k l ) cut. As both terminations are stoichiometric, they are also identical in composition. However, for both trifluorides, the unrelaxed (110)-1 surface is by 0.6 J m 2 more stable than the one of (110)-2 (see E surf unrel . in Table 2). When allowed to relax in atomic positions, the {5,9,8} surface coordinations of (110)-2 rearrange into {6,8,8} (Table S1 and Figure S4). Hence, the surface energy reduces by as much as 1.41 J m 2 for YF3 or 1.18 J m 2 for HoF3. In contrast, termination (110)-1 already starts at a higher surface coordination of {6,9,8}, before it also rearranges into {6,8,8}. According to the lesser degree of rearrangement, its surface energy only reduces by 0.79 J m 2 for YF3 or 0.60 J m 2 for HoF3. Both rearranged terminations are structurally equivalent.
The argumentation in CN surf cannot only be applied to explain the high E surf unrel . of some terminations, but is also partially applicable to the relaxed E surf . Within all YF3 ( h k l ) subsets, except those of (101) and (111), the respective smallest CN surf value correlates with E surf . Thus, the smaller the smallest coordination polyhedron, the less stable the surface and the higher its surface energy. For example, a 4-fold coordination, as present in many surfaces with a fluorine-deficit of two, is only found for the highest E surf within the ( h k l ) subset. Yet, this correlation holds only for the very minimal value within a set of CN surf . No correlation can be found for the remaining, higher CN surf values of the same surface. Therefore, these cannot explain the energetic order of two terminations showing the same smallest CN surf value (as seen e.g., in YF3 (100)-1 and -3 in Table 2). For HoF3, this correlation of surface coordination and stability has two more exceptions. Here, the stability of the least and second-least stable (100) and (001) terminations flip compared to YF3, without any change in CN surf . Prior to surface relaxation, all coordination polyhedrons of YF3 and HoF3 are identical as they share the same crystal bulk structure. After relaxation, this is still true for twenty terminations (Table S1). Only four rearranged terminations differ slightly in surface coordination between YF3 and HoF3. All of these four terminations belong to the less stable surfaces within the respective ( h k l ). All most or second-most stable terminations are identical in surface coordination between YF3 and HoF3. The most stable surface structure termination for each of the seven Miller indices is shown in Figure 2.
As shown in Figure 3, the obtained E surf are similar in magnitude and, within most Miller indices, the order of terminations is equal between YF3 and HoF3. Within convergence, this is also true for the two stoichiometric terminations of (110) and (101), which are very similar in surface energy. For (100) and (001), the least and second-least stable terminations switch their order between YF3 and HoF3. Here, HoF3 prefers the surface with a nominal surface net charge of + 2 over the stoichiometric one. The only difference in termination order between the two compounds, which also affects the most stable surface, is found in (111). For YF3, the most stable surface is (111)-2, which shows a surface coordination of CN surf = { 6 , 5 , 8 , 8 } . Whereas, HoF3 prefers (111)-3 with CN surf = { 6 , 7 , 6 , 9 } (Figure 2). However, both of these terminations are equal in constitution with a fluorine-deficit of 1 per surface.
Even though the order within one ( h k l ) is largely the same between YF3 and HoF3, the order between the different ( h k l ) does change. For YF3, the overall two most stable surfaces are (010)-1 and (001)-2, which are both stoichiometric and give a surface energy of 0.58 J m 2 . This is closely followed by the stoichiometric surface (011)-2. Medium stable surfaces are found for (101)-3 and (111)-2, which are both substoichiometric surfaces missing a single fluorine. The two least stable surfaces (110)-1/-2 and (100)-2 prefer a stoichiometric termination again.
(010)-1 is also the overall most stable surface for HoF3, but the remaining surfaces differ in order. (100)-2, which is the most unstable ( h k l ) in YF3, is the second-most stable one in HoF3. The moderately stable surfaces (001)-2, (011)-2 and (101)-3 have equivalent surface energies within the slab thickness convergence of 0.01 J m 2 . (111)-3 is the second least stable surface. The least stable surface is stoichiometric (110)-1/-2, which is not even part of the Wulff plot (Figure 4).
From the corresponding energies of the respective most stable surfaces shown in Figure 2, the Wulff plots are constructed. A Wulff plot visualizes the thermodynamically most stable crystal shape at quantum chemical conditions of 0 K and vacuum. To test the dependence of surface ratios on the slab thickness convergence, an estimation on the maximum possible error is given in the SI (Table S4).
The largest surface area of over one quarter in YF3 or one third in HoF3 is formed by (010). That HoF3 prefers (010) even stronger is a consequence of its surface energy being 0.11 J m 2 more stable than any other. This contrasts YF3, for which (001) has the same surface energy as (010), as well as a very closely (0.03 J m 2 ) following (011). Nonetheless, the geometric interdependence of surfaces cause a much smaller abundance of only 10% for (001) versus more than double for (011) in YF3. The third most abundant surface in YF3 is (101) with 20%, which is one of the two obtained substoichiometric surfaces. As these same three surfaces (001), (011) and (101) are only medium stable in HoF3, they also only constitute 6%, 13% and 14% of the overall surface. The second substoichiometric surface present in both Wulff plots is (111), which forms an area of 10% in YF3 and 7% in HoF3. Thus, almost one third of the YF3 crystal is made from terminations with a nominal positive net charge of +1. Whereas for HoF3, these are only about a fifth. The two least stable surfaces of YF3 are (100) and (110), which constitute about 7% and 5%. In HoF3, the latter is to such an extend energetically unstable, that it is completely excluded from the Wulff plot. (100), on the other hand, turns to be the second most stable and second most abundant surface in HoF3. It constructs one quarter of the total surface.
The overall scale of most stable surface energies per Miller indices is comparable between YF3 and HoF3. The respective ranges are 0.58–1.03 J m 2 for YF3 and 0.47–1.00 J m 2 for HoF3. However, the resulting average by the Wulff plot is 16% higher for YF3 with Ø E surf = 0.70 J m 2 than for HoF3 with Ø E surf = 0.59 J m 2 . This means, that forming surfaces from the bulk crystal involves a higher thermodynamic barrier in YF3 than in HoF3. This is interesting, as there is no significant difference within the formation enthalpies of the bulk. The difference in average surface energy also hints, that the thermodynamic barrier of crystal nucleation is also higher in YF3 than in HoF3. Though, to accurately predict the nucleation, the nature of the respective precursors and the media needs to be considered.

3.2.1. Bader Charges

Figure 5 shows the partial charges obtained by Bader analysis for all 24 thickness-converged slabs of YF3 or HoF3. As these slabs are built from up to 104 atoms, a large number of very similar charges are obtained. To ease comparison, the Bader charges given in Figure 5 are rounded to 0.1 e. The Bader charges of the central slab atoms reproduce the bulk values with 2.4 e for yttrium, 2.3 e for holmium, as well as 0.8 e for fluorine in both compounds. Only HoF3 (011)-3 shows a marginally increased central slab value of 2.4 e for holmium. In general, the Miller indices do not seem to affect the Bader charges within the analyzed accuracy. The highly ionic partial charge on fluorine does not significantly change for any slab with an overall range of −(0.7–1.0) e. Moreover, all stoichiometric surfaces have practically the same metal partial charges with 2.3–2.4 e. On the contrary, substoichiometric surfaces with a fluorine-deficit of one do all have at least one metal center charged less at the surface. These might be as low as 1.7 e, as found for the third termination of (100) in both trifluorides. Subsequently, surfaces with a fluorine-deficit of two contain even less ionic metal centers at the surface. The least charged metal center is observed again in (100) with only 1.2–1.3 e by the forth termination. The Bader charges suggest that all substoichiometric surfaces missing two fluorine have at at least a single surface metal center in an oxidation state of +II.

3.2.2. Surface Band Gaps

Investigating the electronic properties, we found that the band gap, total DOS and projected DOS of atoms central within the slab converged already at the smallest slab thickness. All band gaps are plotted in Figure S6. For the most stable termination of each ( h k l ), the DOS near the Fermi level are also shown in Figure S7. The surface DOS narrow the band gap within all terminations. For stoichiometric slabs, the direct band gap is reduced from the bulk value of ca. 8 eV to 4–7 eV. Thus, all stoichiometric surfaces remain fully within the insulating regime. In contrast, for substoichiometric terminations, the direct band gap collapses to 0–1 eV, predicting a pseudo-metallic or narrow-band-gap surface. For HoF3 (101) and (111), only one spin direction shows a nearly metallic character. The other stays insulating (5–6 eV). It should be noted that this pseudo-metallic or narrow-band-gap electronic structure at the substoichiometric surfaces might be strongly effected by the chosen neutral 2D-periodic model. However, as this paper is focusing on the relative stability of surfaces, we are not investigating the nature of the band gaps, nor the observed surface magnetism or spin-asymmetric band gaps of some HoF3 surfaces further.

4. Conclusions

The aim of this study was to obtain the relative surface stabilities, in order to find the most abundant surfaces of the two REE trifluorides, YF3 and HoF3 according to their inherent quantum chemical stability. While YF3 can be treated on the DFT level, the 4f-electrons of HoF3 required an extensive electronic structure benchmark evaluating DFT, DFT+U and hybrid DFT against the crystal unit cell, band gaps and Bader charges. On the DFT or DFT+U level, our results show that including the 4f-electrons explicitly within the plane wave expansion worsens the geometrical and band gap results, while the Bader charges stay unaffected. Considering also the experimentally not observed high-spin preference of the 4f-electrons, as well as the strongly increased computational demand, we treated HoF3 by a 4f-in-core DFT+U approach, in which the Hubbard-type correction is applied on the Ho-5d orbitals, which mix into the valence band mainly constructed by F-2p.
From the relaxed bulk, surface models were created for any of the seven low-lying Miller indices. Our analysis included all possible stoichiometric terminations, as well as those showing a small to moderate fluorine-deficit. The surfaces were quantified by Bader charges, band gaps and DOS. From the resulting scope of 24 surfaces, we constructed the first Wulff plots for the whole class of β -YF3-structured compounds.
We found that, within each Miller indices, both trifluorides prefer the same termination with the exception of (111), in which different surface coordinations are favored.
Comparing the different Miller indices, both compounds clearly show stoichiometric (010) as the most stable surface. The preference of the other surfaces, though, varies between the two. The greatest difference is found for (100), which is the second-most stable surface for HoF3, but the second-least stable one for YF3. On average, the surface energy predicted by the Wulff plot is higher for YF3 than for HoF3. This suggests a higher thermodynamical barrier for the formation of YF3 surfaces from the bulk.
In total, one third of the predicted equilibrium crystal shape of YF3 is made from the substoichiometric terminations of (101) and (111) missing a single fluorine per surface. In HoF3, these only constitute a fifth. In the search for the underlying reason between the different fluorine affinity of the two compounds, this different availability of substoichiometric surfaces is an interesting finding. However, to evaluate possible effects, further studies are needed that actually model binding interactions with these surfaces. These should also apply more elaborate binding analysis tools than simple population analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/1996-1944/15/17/6048/s1, Figure S1: Calculated HoF3 direct band gaps with HSE06 (blue), PBE+U (green) and pure PBE (red) applied on 4f-in-core (full markers) or 4f-in-valence (crosses); HSE06/4f-in-valence is not relaxed but done on-top of the crystal structure; the area between the two HSE06 values is highlighted in blue; Figure S2: Calculated HoF3 relaxed unit cell parameters with HSE06 (blue), PBE+U (green) and pure PBE (red) applied on 4f-in-core (full markers) or 4f-in-valence (crosses) compared to the experimental values (horizontal line); Figure S3: Bulk band structure, total DOS (tDOS: gray) and DOS projected onto the metal d band (blue) or fluorine 2p band (green): (a) YF3 (PBE) and (b) HoF3 (PBE+U d /3 eV/4f-in-core); Table S1: Comparison of unrelaxed versus relaxed (or rearranged) slabs in metal coordination number at the surface (CN surf ), as well as in metal centers of the non-surface layers (CN non - surf ) as determined with a bond distance cut-off of 2.6 Å; Figure S4: Effect of surface rearrangement on the stoichiometric surface terminations of (110)-1 (left), (110)-2 (middle) and (101)-2 (right). Atomic positions are shown before (gray) and after relaxation (M: blue, F: green). For the latter, all polyhedra are shown but the one from the initially lowest surface coordination number (CN surf unrel ). Given are the surface energies in J m 2 of the unrelaxed surfaces ( E surf unrel ) for YF3 (first) and HoF3 (second); Table S2: YF3 (PBE) bulk-derived ( E surf bd ) and slab-derived ( E surf sd ) surface energies without (SP) and with atomic position relaxation (OPT); all energies in J m 2 ; the E surf , opt bd values are used within the main paper; Table S3: HoF3 (PBE+U d /3 eV/4f-in-core) bulk-derived ( E surf bd ) and slab-derived ( E surf sd ) surface energies without (SP) and with atomic position relaxation (OPT); all energies in J m 2 ; all magnetic moments in µ B ; the E surf , opt sd values are used within the main paper; Figure S5: Relaxed slab-derived surface energies of HoF3 (PBE+U d /3 eV/4f-in-core). The uncertainty due to slab thickness convergence is given by error bars on each termination; Table S4: Effect of maximal error accumulation due to the convergence in slab thickness of maximal ±0.03 J m 2 for YF3 and ±0.01 J m 2 for HoF3 onto Wulff construction; i denotes the initial value of average surface energy ( Ø E surf ) or surface abundance (% surf ) given by the Wulff plots in the main paper Figure 4; Figure S6: YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core) band gaps of surfaces compared with the respective bulk value (gray). Minimal band gaps, direct or indirect are given by solid bars. In the case, the minimal band gap was found to be indirect, also the direct band gap is given by a transparent bar. For HoF3 (101) and (111) with +1 nominal charges, the band gaps are not spin-symmetric and both direct transitions are given; Figure S7: DOS comparison between the most stable surfaces ordered by their abundance (in %): YF3 (left, PBE), HoF3 (right, PBE+U d /3 eV/4f-in-core), total DOS (tDOS: gray, downscaled to the bulk tDOS) and projected DOS of a single surface atom (Y, Ho: blue; F: green). Substoichiometric slabs with a fluorine-deficit of 1 per surface are framed in green. The top row gives the bulk tDOS with projected DOS of a single bulk atom as reference.

Author Contributions

Conceptualization, B.P.; Formal analysis, J.A. and N.L.; Funding acquisition, B.P.; Investigation, J.A.; Resources, J.A. and B.P.; Supervision, B.P.; Validation, J.A.; Visualization, J.A.; Writing - original draft, J.A.; Writing-review & editing, N.L. and B.P. All authors have read and agreed to the published version of the manuscript.

Funding

The project was funded by the Freie Universität Berlin and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—Project-ID 387284271—CRC 1349—fluorine-specific interactions. funded by the Freie Universität Berlin and the German Science Foundation (DFG) via CRC 1349—Fluorine-Specific Interactions.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

See supplementary materials. Further data can be requested from the authors.

Acknowledgments

The authors thank the North-German Supercomputing Alliance (Norddeutscher Verbund zur Förderung des Hoch-und Höchstleistungsrechnens HLRN) and the Zentraleinrichtung für Datenverarbeitung (ZEDAT) at the Freie Universität Berlin for computational resources, the German Science Foundation (DFG) for funding within the CRC 1349—Fluorine-Specific Interactions, as well as Anselm Loges for scientific discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shannon, R.D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Cryst. A 1976, 32, 751–767. [Google Scholar] [CrossRef]
  2. Bau, M.; Dulski, P. Comparative study of yttrium and rare-earth element behaviours in fluorine-rich hydrothermal fluids. Contrib. Mineral. Petrol. 1995, 119, 213–223. [Google Scholar] [CrossRef]
  3. Minuzzi, O.R.R.; Bastos Neto, A.C.; Formoso, M.L.L.; Andrade, S.; Janasi, V.A.; Flores, J.A. Rare earth element and yttrium geochemistry applied to the genetic study of cryolite ore at the Pitinga Mine (Amazon, Brazil). An. Acad. Bras. Cienc. 2008, 80, 719–733. [Google Scholar] [CrossRef]
  4. Goodenough, K.M.; Schilling, J.; Jonsson, E.; Kalvig, P.; Charles, N.; Tuduri, J.; Deady, E.A.; Sadeghi, M.; Schiellerup, H.; Müller, A.; et al. Europe’s rare earth element resource potential: An overview of REE metallogenetic provinces and their geodynamic setting. Ore Geol. Rev. 2016, 72, 838–856. [Google Scholar] [CrossRef]
  5. Kaminski, A. Laser Crystals: Their Physics and Properties, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 1990. [Google Scholar]
  6. Kollia, Z.; Sarantopoulou, E.; Cefalas, A.C.; Nicolaides, C.A.; Naumov, A.K.; Semashko, V.V.; Abdulsabirov, R.Y.; Korableva, S.L.; Dubinskii, M.A. Vacuum-ultraviolet interconfigurational 4f3 → 4f25d absorption and emission studies of the Nd3+ ion in KYF, YF, and YLF crystal hosts. J. Opt. Soc. Am. B 1995, 12, 782–785. [Google Scholar] [CrossRef]
  7. Lage, M.M.; Righi, A.; Matinaga, F.M.; Gesland, J.Y.; Moreira, R.L. Raman-spectroscopic study of lanthanide trifluorides with the - YF3 structure. J. Phys. Condens. Matter 2004, 16, 3207–3218. [Google Scholar] [CrossRef]
  8. Atencio, D.; Bastos Neto, A.C.; Pereira, V.P.; Ferron, J.T.M.M.; Hoshino, M.; Moriyama, T.; Watanabe, Y.; Miyawaki, R.; Coutinho, J.M.V.; Andrade, M.B.; et al. Waimirite-(Y), orthorhombic YF3, a new mineral from the Pitinga mine, Presidente Figueiredo, Amazonas, Brazil and from Jabal Tawlah, Saudi Arabia: Description and crystal structure. Mineral. Mag. 2015, 79, 767–780. [Google Scholar] [CrossRef]
  9. Sarantopoulou, E.; Kollia, Z.; Cefalas, A.C. YF3:Nd3+, Pr3+, Gd3+ wide band gap crystals as optical materials for 157-nm photolithography. Opt. Mater. 2001, 18, 23–26. [Google Scholar] [CrossRef]
  10. O’Keeffe, M. Ionic Conductivity of Yttrium Fluoride and Lutetium Fluoride. Science 1973, 180, 1276–1277. [Google Scholar] [CrossRef]
  11. Pogorenko, Y.V.; Pshenychnyi, R.M.; Lutsyk, V.I.; Omel’chuk, A.O. Transport Properties of Aliovalent Substitution Solid Solutions of the System (1-x)PbF2-xYF3-SnF2. IOP Conf. Ser. Mater. Sci. Eng. 2017, 175, 012039. [Google Scholar] [CrossRef]
  12. Fedorov, P.P.; Sorokin, N.I. Stabilization of the α-YF3 structure type by isomorphous substitutions. Inorg. Mater. 2017, 53, 1307–1311. [Google Scholar] [CrossRef]
  13. Cui, X.; Hu, T.; Wang, J.; Zhang, J.; Zhong, X.; Chen, Y.; Li, X.; Yang, J.; Gao, C. Ionic transportation and dielectric properties of YF3:Eu3+ nanocrystals. Nanomaterials 2018, 8, 995. [Google Scholar] [CrossRef] [PubMed]
  14. Dai, X.; Komatsu, Y.; Shimizu, R.; Hitosugi, T. Diffusion of F atoms from fluoride substrates promotes the epitaxial growth of metal fluorides. Appl. Phys. Express 2020, 13, 85507. [Google Scholar] [CrossRef]
  15. Hoard, R.; Mance, S.; Leber, R.; Dalder, E.; Chaplin, M.; Blair, K.; Nelson, D.; Dyke, D.V. Field enhancement of a 12.5-T magnet using holmium poles. IEEE Trans. Magn. 1985, 21, 448–450. [Google Scholar] [CrossRef]
  16. González-Mancebo, D.; Becerro, A.I.; Rojas, T.C.; García-Martín, M.L.; de la Fuente, J.M.; Ocaña, M. HoF3 and DyF3 Nanoparticles as Contrast Agents for High-Field Magnetic Resonance Imaging. Part. Part. Syst. Charact. 2017, 34, 1700116. [Google Scholar] [CrossRef]
  17. Rudzitis, E.; Feder, H.M.; Hubbard, W.N. Fluorine bomb calorimetry. XI. The enthalpy of formation of yttrium trifluoride. J. Phys. Chem. 1965, 69, 2305–2307. [Google Scholar] [CrossRef]
  18. Spedding, F.H.; Henderson, D.C. High-temperature heat contents and related thermodynamic functions of seven trifluorides of the rare earths: Y, La, Pr, Nd, Gd, Ho, and Lu. J. Chem. Phys. 1971, 54, 2476–2483. [Google Scholar] [CrossRef]
  19. Loges, A.; Migdisov, A.A.; Wagner, T.; Williams-Jones, A.E.; Markl, G. Fluoride complexation of hafnium under hydrothermal conditions. In Proceedings of the Goldschmidt Conference, Montréal, QC, Canada, 24–29 June 2012; p. 2030. [Google Scholar]
  20. Loges, A.; Migdisov, A.A.; Wagner, T.; Williams-Jones, A.E.; Markl, G. An experimental study of the aqueous solubility and speciation of Y(III) fluoride at temperatures up to 250C. Geochim. Cosmochim. Acta 2013, 123, 403–415. [Google Scholar] [CrossRef]
  21. Bau, M. Controls on the fractionation of isovalent trace elements in magmatic and aqueous systems: Evidence from Y/Ho, Zr/Hf, and lanthanide tetrad effect. Contrib. Mineral. Petrol. 1996, 123, 323–333. [Google Scholar] [CrossRef]
  22. Zalkin, A.; Templeton, D.H. The Crystal Structures of YF3 and Related Compounds. J. Am. Chem. Soc. 1953, 75, 2453–2458. [Google Scholar] [CrossRef]
  23. Cheetham, A.K.; Norman, N. The Structures of Yttrium and Bismuth Trifluorides by Neutron Diffraction. Acta Chem. Scand. 1974, 28, 55–60. [Google Scholar] [CrossRef]
  24. Khairulin, R.A.; Stankus, S.V.; Lyapunov, K.M. The Thermal Properties of Holmium Trifluoride at High Temperatures. High Temp. 2000, 38, 149–151. [Google Scholar] [CrossRef]
  25. Loges, A.; Freie Universität Berlin, Berlin, Germany. Private communication, 2020.
  26. Vologzhanina, A.V.; Pushkin, D.V.; Serezhkin, V.N. Coordination polyhedra LnFn (Ln = La-Lu) in crystal structures. Russ. J. Inorg. Chem. 2006, 51, 747–758. [Google Scholar] [CrossRef]
  27. Rotereau, K.; Daniel, P.; Desert, A.; Gesland, J.Y. The high-temperature phase transition in samarium fluoride: Structural and vibrational investigation. J. Phys. Condens. Matter 1998, 10, 1431–1446. [Google Scholar] [CrossRef]
  28. Zinchenko, V.; Efryushina, N.; Eryomin, O.; Markiv, V.; Belyavina, N.; Mozkova, O.; Zakharenko, M. Synthesis, structure and optical properties of EuF3 film-forming material. J. Alloys Compd. 2002, 347, L1–L3. [Google Scholar] [CrossRef]
  29. Piotrowski, M.; Murasik, A. Exchange and Dipolar Interactions in TbF3. Phys. Status Solidi A 1985, 89, 571–580. [Google Scholar] [CrossRef]
  30. Piotrowski, M.; Ptasiewicz-ba̧k, H.; Murasik, A. The crystal structures of HoF3 and TbF3 by neutron diffraction. Phys. Status Solidi A 1979, 55, K163–K166. [Google Scholar] [CrossRef]
  31. Bukvetskii, B.; Garashina, L. Crystal-chemical investigation of the orthorhombic trifluorides of samarium, holmium, and ytterbium. Russ. J. Coord. Chem. 1977, 3, 1024–1029. [Google Scholar]
  32. Peterson, J.R.; Cunningham, B.B. Crystal structures and lattice parameters of the compounds of berkelium—IV berkelium trifluoride. J. Inorg. Nucl. Chem. 1968, 30, 1775–1784. [Google Scholar] [CrossRef]
  33. Cunningham, B.; Ehrlich, P. USAEC Report UCRL-20426; Technical Report; USAEC: Houston, TX, USA, 1970. [Google Scholar]
  34. Stevenson, J.N.; Peterson, J.R. The trigonal and orthorhombic crystal structures of CfF3 and their temperature relationship. J. Inorg. Nucl. Chem. 1973, 35, 3481–3486. [Google Scholar] [CrossRef]
  35. Morss, L.R.; Edelstein, N.M.; Fuger, J. (Eds.) The Chemistry of the Actinide and Transactinide Elements, 3rd ed.; Springer International Publishing: Dordrecht, The Netherlands, 2008; p. 1468. [Google Scholar]
  36. Haschke, J.M.; Hodges, A.E.; Lucas, R.L. Equilibrium and structural properties of the PuH system. J. Less-Common Met. 1987, 133, 155–166. [Google Scholar] [CrossRef]
  37. Greis, O.; Martinez-Ripoll, M. Darstellung, Temperaturverhalten und Kristallstruktur von BiF3. Z. Anorg. Allg. Chem. 1977, 436, 105–112. [Google Scholar] [CrossRef]
  38. Deng, Z.; Wei, F.; Wu, Y.; Seshadri, R.; Cheetham, A.K.; Canepa, P. Understanding the Structural and Electronic Properties of Bismuth Trihalides and Related Compounds. Inorg. Chem. 2020, 59, 3377–3386. [Google Scholar] [CrossRef]
  39. Ye, X.; Chen, J.; Engel, M.; Millan, J.A.; Li, W.; Qi, L.; Xing, G.; Collins, J.E.; Kagan, C.R.; Li, J.; et al. Competition of shape and interaction patchiness for self-assembling nanoplates. Nat. Chem. 2013, 5, 466–473. [Google Scholar] [CrossRef]
  40. Griffiths, C.L.; Macdonald, J.E.; Williams, R.H. The growth and characterisation of epitaxial insulating HoF3 layers on silicon. Appl. Surf. Sci. 1992, 56–58, 782–788. [Google Scholar] [CrossRef]
  41. Guo, F.; Li, H.; Zhang, Z.; Meng, S.; Li, D. Synthesis of mesoporous YF3 nanoflowers via solvent extraction route. Mater. Sci. Eng. B 2009, 163, 134–137. [Google Scholar] [CrossRef]
  42. Wen, C.; Sun, L.; Yan, J.; Liu, Y.; Song, J.; Zhang, Y.; Lian, H.; Kang, Z. Mesoporous rare earth fluoride nanocrystals and their photoluminescence properties. J. Colloid Interface Sci. 2011, 357, 116–120. [Google Scholar] [CrossRef]
  43. Tian, Y.; Chen, B.; Li, X.; Zhang, J.; Tian, B.; Sun, J.; Cheng, L.; Zhong, H.; Zhong, H.; Hua, R. Solvothermal synthesis and tunable luminescence of Tb3+, Eu3+ codoped YF3 nano- and micro-crystals with uniform morphologies. J. Solid State Chem. 2012, 196, 187–196. [Google Scholar] [CrossRef]
  44. Chen, J.x.; Wang, X.p.; Wang, L.J.; Yang, X.w.; Yang, Y. White electroluminescence of diamond—HoF3- diamond composite film. J. Lumin. 2020, 224, 117310. [Google Scholar] [CrossRef]
  45. Nanda, S.S.; Nayak, P.; Goutam, U.K.; Dash, S. Influence of Eu3+ on the Structure and Photophysical Properties in (Y,Gd)F3 Nanophosphors. J. Fluoresc. 2021, 31, 129–139. [Google Scholar] [CrossRef]
  46. Pandharkar, R.; Becker, C.; Budau, J.; Kaawar, Z.; Paulus, B. A Computational Study of AlF3 and ACF Surfaces. Inorganics 2018, 6, 124. [Google Scholar] [CrossRef]
  47. Schnellbügel, A.; Anton, R. On background subtraction for quantitative analysis of X-ray photoelectron spectra of rare earth fluorides. Surf. Sci. 2001, 492, 305–314. [Google Scholar] [CrossRef]
  48. Wulff, G. XXV. Zur Frage der Geschwindigkeit des Wachstums und der Auflösung der Krystallflächen. Z. Kristallogr.-Cryst. Mater. 1901, 34, 449–530. [Google Scholar] [CrossRef]
  49. Laue, M.V. Der Wulffsche Satz für die Gleidigewichtsform von Kristallen. Z. Kristallogr.-Cryst. Mater. 1943, 105, 124–133. [Google Scholar] [CrossRef]
  50. Kresse, G.; Hafner, J. Norm-conserving and ultrasoft pseudopotentials for first-row and transition elements. J. Phys. Condens. Matter 1994, 6, 8245–8257. [Google Scholar] [CrossRef]
  51. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  52. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  53. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  54. Dudarev, S.L.; Botton, G.A.; Savrasov, S.Y.; Humphreys, C.J.; Sutton, A.P. Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. B 1998, 57, 1505–1509. [Google Scholar] [CrossRef]
  55. Krukau, A.V.; Vydrov, O.A.; Izmaylov, A.F.; Scuseria, G.E. Influence of the exchange screening parameter on the performance of screened hybrid functionals. J. Chem. Phys. 2006, 125, 224106. [Google Scholar] [CrossRef]
  56. Blöchl, P.E.; Jepsen, O.; Andersen, O.K. Improved tetrahedron method for Brillouin-zone integrations. Phys. Rev. B 1994, 49, 16223–16233. [Google Scholar] [CrossRef]
  57. Evans, D.S.; Baynor, G.V. Lattice spacings in thorium—Yttrium alloys. J. Nucl. Mater. 1960, 2, 209–215. [Google Scholar] [CrossRef]
  58. Pechan, M.J.; Stassis, C. Magnetic structure of holmium. J. Appl. Phys. 1984, 55, 1900–1902. [Google Scholar] [CrossRef]
  59. Kerker, G.P. Efficient iteration scheme for self-consistent pseudopotential calculations. Phys. Rev. B 1981, 23, 3082–3084. [Google Scholar] [CrossRef]
  60. Sun, W.; Ceder, G. Efficient creation and convergence of surface slabs. Surf. Sci. 2013, 617, 53–59. [Google Scholar] [CrossRef]
  61. Tran, R.; Xu, Z.; Radhakrishnan, B.; Winston, D.; Sun, W.; Persson, K.A.; Ong, S.P. Surface energies of elemental crystals. Sci. Data 2016, 3, 160080. [Google Scholar] [CrossRef]
  62. Rahm, J.M.; Erhart, P. WulffPack: A Python package for Wulff constructions. J. Open Source Softw. 2020, 5, 1944. [Google Scholar] [CrossRef]
  63. Momma, K.; Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  64. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  65. Rogers, E.; Dorenbos, P.; van der Kolk, E. Systematics in the optical and electronic properties of the binary lanthanide halide, chalcogenide and pnictide compounds: An overview. New J. Phys. 2011, 13, 093038. [Google Scholar] [CrossRef]
  66. Bader, R.F.W. Atoms in molecules. Acc. Chem. Res. 1985, 18, 9–15. [Google Scholar] [CrossRef]
  67. Henkelman, G.; Arnaldsson, A.; Jónsson, H. A fast and robust algorithm for Bader decomposition of charge density. Comp. Mat. Sci. 2006, 36, 354–360. [Google Scholar] [CrossRef]
  68. Sanville, E.; Kenny, S.D.; Smith, R.; Henkelman, G. Improved grid-based algorithm for Bader charge allocation. J. Comput. Chem. 2007, 28, 899–908. [Google Scholar] [CrossRef] [PubMed]
  69. Tang, W.; Sanville, E.; Henkelman, G. A grid-based Bader analysis algorithm without lattice bias. J. Phys. Condens. Matter 2009, 21, 84204. [Google Scholar] [CrossRef] [PubMed]
  70. Yu, M.; Trinkle, D.R. Accurate and efficient algorithm for Bader charge integration. J. Chem. Phys. 2011, 134, 64111. [Google Scholar] [CrossRef]
  71. Brown, P.J.; Forsyth, J.B.; Hansen, P.C.; Leask, M.J.M.M.; Ward, R.C.C.C.; Wells, M.R. Neutron diffraction determination of magnetic order in holmium trifluoride, HoF3. J. Phys. Condens. Matter 1990, 2, 4471–4484. [Google Scholar] [CrossRef]
  72. Boettger, J.C. Nonconvergence of surface energies obtained from thin-film calculations. Phys. Rev. B 1994, 49, 16798–16800. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The orthorhombic unit cell of β -YF3 with Pnma symmetry. Lattice constants are given in Table 1. The distorted tricapped trigonal prisms formed by nine fluorides (green) around each yttrium (cyan) are visualized by transparent, cyan planes.
Figure 1. The orthorhombic unit cell of β -YF3 with Pnma symmetry. Lattice constants are given in Table 1. The distorted tricapped trigonal prisms formed by nine fluorides (green) around each yttrium (cyan) are visualized by transparent, cyan planes.
Materials 15 06048 g001
Figure 2. Most stable terminations of the relaxed surface structures: the coordination number of the surface metals (CNsurf) and the surface energies in J m−2 (Esurf) are given. The first entry corresponds to YF3 and the second to HoF3. The mean of both values corresponds to the given order from top left to bottom right. Each (hkl) slab is rotated in a way to show the surface coordination best. For (111), two surfaces are given, as (111)-2 is preferred by YF3 and (111)-3 by HoF3.
Figure 2. Most stable terminations of the relaxed surface structures: the coordination number of the surface metals (CNsurf) and the surface energies in J m−2 (Esurf) are given. The first entry corresponds to YF3 and the second to HoF3. The mean of both values corresponds to the given order from top left to bottom right. Each (hkl) slab is rotated in a way to show the surface coordination best. For (111), two surfaces are given, as (111)-2 is preferred by YF3 and (111)-3 by HoF3.
Materials 15 06048 g002
Figure 3. Relaxed surface energies of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core) of all 24 terminations. The surfaces are color-coded to their nominal surface net charge ( q surf ) in e of 0 (blue), + 1 (green) or + 2 (orange). This magnified plot does not show the HoF3 (110)-3 value of 2.09 J m 2 .
Figure 3. Relaxed surface energies of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core) of all 24 terminations. The surfaces are color-coded to their nominal surface net charge ( q surf ) in e of 0 (blue), + 1 (green) or + 2 (orange). This magnified plot does not show the HoF3 (110)-3 value of 2.09 J m 2 .
Materials 15 06048 g003
Figure 4. Wulff Plots of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core) from relaxed surfaces. The percentage shows the relative abundance of each surface, which is also given in Table 2.
Figure 4. Wulff Plots of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core) from relaxed surfaces. The percentage shows the relative abundance of each surface, which is also given in Table 2.
Materials 15 06048 g004
Figure 5. Bader charges rounded to 0.1 e for all slabs of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core): the Bader charges of the bulk are highlighted by the dotted line (Y: 2.4 e, Ho: 2.3 e, F: 0.8 e). Terminations are differently colored by their formal surface net charge ( q surf ) in e for 0 (blue cycle), + 1 (green diamond) and + 2 (orange cross).
Figure 5. Bader charges rounded to 0.1 e for all slabs of YF3 (left, PBE) and HoF3 (right, PBE+U d /3 eV/4f-in-core): the Bader charges of the bulk are highlighted by the dotted line (Y: 2.4 e, Ho: 2.3 e, F: 0.8 e). Terminations are differently colored by their formal surface net charge ( q surf ) in e for 0 (blue cycle), + 1 (green diamond) and + 2 (orange cross).
Materials 15 06048 g005
Table 1. Comparison of the relaxed unit cells to experiment including standard deviation in parentheses [22,23]. Given is the absolute difference ( Δ ), as well as the deviation from experiment in percentage ( Δ % ).
Table 1. Comparison of the relaxed unit cells to experiment including standard deviation in parentheses [22,23]. Given is the absolute difference ( Δ ), as well as the deviation from experiment in percentage ( Δ % ).
YF3 (PBE)HoF3 (PBE+U d /3 eV/4f-in-Core)
a (Å) b (Å) c (Å) V 3 ) a (Å) b (Å) c (Å) V 3 )
calc.6.32156.80594.3300186.296.41646.87964.3440191.76
exp. [22]6.353(3)6.850(3)4.393(3)191.26.404(3)6.875(3)4.379(3)192.8
Δ 0.0320.0440.0634.90.0120.0050.0351.0
Δ % 0.5%0.6%1.4%2.6%0.2%0.1%0.8%0.5%
exp. [23]6.3537(7)6.8545(7)4.3953(5)191.42
Δ 0.03220.04860.06535.13
Δ % 0.5%0.7%1.5%2.7%
Table 2. The YF3 (PBE) and HoF3 (PBE+U d /3 eV/4f-in-core) surfaces with respective terminations (term.), slab thickness in layers of formula units without terminal F-deficit ( L MF 3 ), nominal surface net charge ( q surf ) in e, surface energies of relaxed ( E surf ) and unrelaxed slabs ( E surf unrel . ) in J m 2 , as well as the relaxed surface metal coordination number (CN surf ). The lowest surface energies per ( h k l ) cut are highlighted in bold. For these, also the abundance obtained by the Wulff plot (% surf ) is given.
Table 2. The YF3 (PBE) and HoF3 (PBE+U d /3 eV/4f-in-core) surfaces with respective terminations (term.), slab thickness in layers of formula units without terminal F-deficit ( L MF 3 ), nominal surface net charge ( q surf ) in e, surface energies of relaxed ( E surf ) and unrelaxed slabs ( E surf unrel . ) in J m 2 , as well as the relaxed surface metal coordination number (CN surf ). The lowest surface energies per ( h k l ) cut are highlighted in bold. For these, also the abundance obtained by the Wulff plot (% surf ) is given.
L MF 3 CN surf E surf ( E surf unrel . )% surf
( h k l )term. q surf YF3HoF3YF3HoF3YF3HoF3YF3HoF3
(100)1020245,91.61 (2.87)0.93 (1.48)
2022266,91.03 (2.02)0.58 (0.96)7%25%
3+120245,81.24 (1.61)0.62 (0.68)
4+222264,71.79 (2.14)0.87 (0.90)
(010)1010128,80.58 (0.84)0.47 (0.49)26%34%
2+210126,61.80 (2.05)1.52 (1.52)
(001)1020245,8,8,91.23 (2.45)1.37 (2.25)
2022266,7,8,90.58 (1.39)0.67 (1.16)10%6%
3+222264,5,8,91.27 (1.70)1.23 (1.29)
(110)1020246,8,81.01 (1.80)0.99 (1.59)5%0%
2022266,8,81.00 (2.41)1.00 (2.18)
3+222264,6,94,6,81.42 (1.73)2.09 (1.36)
(101)1020246,7,8,80.82 (1.48)0.89 (1.33)
2020246,6,8,80.82 (3.34)0.88 (3.17)
3+120246,7,8,80.76 (1.16)0.69 (0.89)20%14%
4+122265,6,7,95,6,8,81.07 (2.10)1.03 (1.70)
5+220244,5,8,85,6,8,80.98 (1.39)0.99 (0.99)
(011)1010127,7,9,90.78 (1.30)0.81 (1.14)
2010126,6,8,80.61 (1.32)0.68 (1.15)22%13%
3+210124,4,8,81.25 (1.68)1.35 (1.38)
(111)1020246,7,7,87,7,8,81.02 (3.46)0.87 (3.29)
2+120245,6,8,80.83 (1.30)0.82 (1.04)10%
3+122266,6,7,91.05 (1.70)0.75 (1.11) 7%
4+220245,5,7,70.93 (1.22)0.95 (1.13)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anders, J.; Limberg, N.; Paulus, B. First Principle Surface Analysis of YF3 and Isostructural HoF3. Materials 2022, 15, 6048. https://doi.org/10.3390/ma15176048

AMA Style

Anders J, Limberg N, Paulus B. First Principle Surface Analysis of YF3 and Isostructural HoF3. Materials. 2022; 15(17):6048. https://doi.org/10.3390/ma15176048

Chicago/Turabian Style

Anders, Jennifer, Niklas Limberg, and Beate Paulus. 2022. "First Principle Surface Analysis of YF3 and Isostructural HoF3" Materials 15, no. 17: 6048. https://doi.org/10.3390/ma15176048

APA Style

Anders, J., Limberg, N., & Paulus, B. (2022). First Principle Surface Analysis of YF3 and Isostructural HoF3. Materials, 15(17), 6048. https://doi.org/10.3390/ma15176048

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

Article Metrics

Back to TopTop