*Article* **Quantum-Mechanical Study of Nanocomposites with Low and Ultra-Low Interface Energies**

#### **Martin Friák 1,\*, David Holec <sup>2</sup> and Mojmír Šob 3,1,4**


Received: 27 November 2018; Accepted: 12 December 2018; Published: 15 December 2018

**Abstract:** We applied first-principles electronic structure calculations to study structural, thermodynamic and elastic properties of nanocomposites exhibiting nearly perfect match of constituting phases. In particular, two combinations of transition-metal disilicides and one pair of magnetic phases containing the Fe and Al atoms with different atomic ordering were considered. Regarding the disilicides, nanocomposites MoSi2/WSi2 with constituents crystallizing in the tetragonal C11*<sup>b</sup>* structure and TaSi2/NbSi2 with individual phases crystallizing in the hexagonal C40 structure were simulated. Constituents within each pair of materials exhibit very similar structural and elastic properties and for their nanocomposites we obtained ultra-low (nearly zero) interface energy (within the error bar of our calculations, i.e., about 0.005 J/m2). The interface energy was found to be nearly independent on the width of individual constituents within the nanocomposites and/or crystallographic orientation of the interfaces. As far as the nanocomposites containing Fe and Al were concerned, we simulated coherent superlattices formed by an ordered Fe3Al intermetallic compound and a disordered Fe-Al phase with 18.75 at.% Al, the *α*-phase. Both phases were structurally and elastically quite similar but the disordered *α*-phase lacked a long-range periodicity. To determine the interface energy in these nanocomposites, we simulated seven different distributions of atoms in the *α*-phase interfacing the Fe3Al intermetallic compound. The resulting interface energies ranged from ultra low to low values, i.e., from 0.005 to 0.139 J/m2. The impact of atomic distribution on the elastic properties was found insignificant but local magnetic moments of the iron atoms depend sensitively on the type and distribution of surrounding atoms.

**Keywords:** MoSi2; WSi2; TaSi2; NbSi2; elasticity; ab initio; interface energies; Fe3Al; disorder

#### **1. Introduction**

Ever increasing demand for energy-conversion units exhibiting a higher efficiency leads to increasing operating temperatures in these systems and, therefore, new materials, which would sustain such conditions, are needed. Because the development of these materials is highly complex and multi-faceted, combinations of often mutually-conflicting properties are rarely found in a single-phase structures. Composites then represent a critically important class of materials. In particular, coherent nanocomposites require optimum matching of constituting phases for their stability and rather low interface energies. In our study, we addressed three nanocomposites combining materials intended for high or elevated temperature applications: two pairs of transition-metal disilicides that are predicted

to possess ultra-low interface energies and then a pair of two different phases from the Fe-Al system, which are expected to have the interface energies ranging from ultra low to low values.

Regarding the transition-metal silicides, they are currently considered as very promising bases for future high-temperature structural materials (see, e.g., Refs. [1–4]), in particular for operational temperatures above those of Ni-based superalloys. At high temperatures, transition-metal silicides are known to combine the ductility and thermal conductivity of metals with high strength and corrosion resistance of ceramics. As far as composites combining them are concerned, MoSi2/WSi2 composite powders with different phase compositions are fabricated via a self-propagating high-temperature synthesis (SHS) method [5]. This approach is widely recognized as an effective manufacturing strategy for the fabrication of materials applied in high-temperature fields, in particular for refractories such as transition-metal carbides, nitrides, silicides, and borides [6]. Preparation of MoSi2/WSi2 composites using elemental powders of Mo, W and Si by the thermal explosion mode of SHS have been theoretically calculated and investigated by experiments in Ref. [7]. Phase composition and microchemical area analyses were conducted by XRD, SEM and EDAX methods. Pure MoSi2/WSi2 composites are fabricated by the thermal explosion mode of SHS, and MoSi2/WSi2 exists as a solid solution of (Mo*x*,W1−*<sup>x</sup>*)Si2 but the chemical elements inside of individual grains are not uniformly distributed. As an alternative processing route, five kinds of WSi2/MoSi2 composites are fabricated by mechanical alloying in [8]. WSi2-reinforced MoSi2 composites are successfully prepared also by mechanical activation followed by in situ reactive spark plasma sintering of Mo, Si, and W elemental powders in [9]. The addition of W to the reactants leads to a finer microstructure than that obtained using pure MoSi2, resulting in a significant improvement of mechanical properties. The Vicker's hardness of 20 vol % WSi2/MoSi2 is as high as 16.47 GPa. Nanocomposite of (Mo,W)Si2/WSi2 was synthesized via mechanical alloying (MA) and heat treatment in Ref. [10]. Increasing the milling time to 80 h followed by the post-annealing at 1000 ◦C caused the complete formation of (Mo,W)Si2/WSi2 nanocomposite. MoSi2/WSi2 composites were successfully prepared by pressureless sintering from mechanically-assisted combustion synthesized powders in [11].

Motivated by the above-mentioned studies of MoSi2/WSi2 (nano)composites containing structurally and elastically very similar pair of materials crystallizing in the C11*<sup>b</sup>* structure, we also assessed another pair of matching materials, TaSi2 and NbSi2, which are crystallizing in the C40 structure.

Finally, as a system with rather low interface energy, we studied nanocomposites formed by two phases from the Fe-Al binary system [12–15]. These materials are also considered as promising for elevated-temperature applications and as such they are intensively studied [16–40]. A sub-class of Fe-Al-based materials is represented by composites consisting of an ordered Fe3Al with the D03 structure and a disordered Fe-Al solid solution with about 18-19 at.% Al. The existence of these composites can be experimentally proved using, for example, Mössbauer spectroscopy [41] or transmission electron microscopy (TEM) techniques. The latter are sensitive to anti-phase boundaries (APBs), which have a different character in Fe3Al and the Fe-Al phase [42–45].

#### **2. Methods**

Our quantum-mechanical calculations were performed within the framework of density functional theory [46,47] using the Vienna Ab initio Simulation Package (VASP) [48,49] and projector augmented wave (PAW) pseudopotentials [50,51]. When studying the transition-metal disilicides, the exchange and correlation energy was treated in the local density approximation (LDA) [52] but, in the case of phases containing the Fe and Al atoms, the generalized gradient approximation (GGA) parameterized by Perdew and Wang [53] (PW91) with the Vosko–Wilk–Nusair correction [54] was necessary to correctly reproduce the ground-state D03 structure of Fe3Al. Regarding the MoSi2 and WSi2, we used a plane-wave energy cut-off of 400 eV and the k-point Monkhorst–Pack [55] meshes contained 20 × 20 × 2 (10 × 10 × 8) k-points in the case of 24-atom supercells modeling the superlattices with stacking along the [001] ([100] and [110]) directions. The calculations for TaSi2 and NbSi2 required the cut-off energy of 480 eV and the 12 × 12 × 4 k-point mesh in the case of 18-atom supercells. When computing Fe-Al-based nanocomposites, the cut-off energy was equal to 350 eV and the sampling of the Brillouin zone was done using Monkhorst–Pack grids with 10 × 10 × 6 k-points in the case of computational supercells containing 32 atoms. All calculations had an estimated error-bar of about 0.001 eV/atom. Finally, the second-order elastic constants were determined using the stress-strain method [56].

#### **3. Results**

The first type of studied nanocomposite containing transition-metal disilicides is visualized in Figure 1a. WSi2 and MoSi2, which crystallize in the tetragonal C11*<sup>b</sup>* structure, form a coherent nanocomposite where two conventional cells of each materials are stacked one on top of the other along the [001] direction (the interfaces are perpendicular to this direction) and alternate. It should be emphasized that, due to the periodic boundary conditions, which are applied to all nanocomposites in our calculations, the simulated nanocomposites form so-called superlattices [57–78] when the atomic planes continue from one phase into another.

**Figure 1.** Visualization of a WSi2/MoSi2 nanocomposite (a superlattice) with the stacking along (and the interfaces perpendicular to) the [001] direction within the C11*<sup>b</sup>* lattice (**a**) accompanied with directional dependences of the Young's modulus of bulk WSi2 (**b**), bulk MoSi2 (**c**) and the composite WSi2/MoSi2 (**d**) formed out of them. The parts (**b**–**d**) were visualized by the SC-EMA [79–81] library (scema.mpie.de) based on ab initio computed elastic constants.

As seen in Table 1, both disilicides have very similar lattice parameters *a* and *c* describing their tetragonal structure and also quite similar elastic properties. Our theoretical values are in excellent agreement with both experimental data and previous calculations. The elasticity of bulk phases is conveniently visualized in Figure 1b,c in the form of directional dependences of the Young's modulus. The values for this composite with an equal amount of both phases are listed in Table 1 at the line for the composition (WSi2)4(MoSi2)4.



As a consequence of the similarity of the elastic constants of both constituents, the overall elastic properties of their nanocomposites are quite similar, too.

**Figure 2.** Schematic visualization of supercells modeling various nanocomposites with the stacking along (and the interfaces perpendicular to) the [001] direction. The mutual ratio of the constituents varies from 1:7, i.e., (WSi2)1(MoSi2)7, in the case of sub-figure (**a**), via 2:6 (**b**) and 3:5 (**c**) further for the opposite ratios 5:3 (**d**) and 6:2 (**e**) to 7:1, i.e., (WSi2)7(MoSi2)1, in the case of sub-figure (**f**).

Next, we evaluated the interface energy of the nanocomposites according to the formula:

$$\gamma((\text{WSi}\_2)\_m/(\text{MoSi}\_2)\_n) = \frac{E((\text{WSi}\_2)\_m/(\text{MoSi}\_2)\_n) - m \cdot E(\text{WSi}\_2) - n \cdot E(\text{MoSi}\_2)}{(2 \cdot A)}\tag{1}$$

using the energy *E*((WSi2)*m*/(MoSi2)*n*) of the supercell modeling the nanocomposite (WSi2)*m*/(MoSi2)*n*, which contains *m* formula units of WSi2 and *n* formula units of MoSi2, the energies of bulk phases of each constituent *E*(WSi2) and *E*(MoSi2) and the area of the interface *A*. Importantly, very probably due to the similarity of both constituents, the interface energy was found to be ultra low, essentially zero within the error bar of our calculations, i.e., the energy differences in Equation (1) are smaller than 0.001 eV/atom and the interface energies are then smaller than about 0.005 J/m2.

To examine how the ultra-low interface energies depend on the width of the layers containing individual constituents within the nanocomposite as well as on the ratio of amount of both materials, we simulated a series of seven other superlattices with varying width of the constituents, as visualized in Figure 2. The calculated structural and elastic parameters are summarized in Table 1. The lattice parameter *a* increased quite monotonously from the value calculated for bulk MoSi2 to that obtained for bulk WSi2. The changes of the lattice parameter *c* are much smaller because the values found for bulk MoSi2 and WSi2 are only very slightly different. Rather monotonous changes appear also for the elastic constants with all of them increasing when increasing the amount of elastically stiffer WSi2. Importantly, all the studied nanocomposites have again ultra low interface energies, which are essentially zero within the error-bar of our calculations. Next, as seen in the Appendix, the ultra low interface energies were nearly independent of the crystallographic orientation of the interface.

After examining the WSi2/MoSi2 nanocomposites, which were experimentally reported to exist, we next drew our attention to another class of nanocomposites containing transition-metal disilicides, which were structurally and elastically very similar. The studied TaSi2 and NbSi2 crystallize in the hexagonal C40 structure and, therefore, we simulated a superlattice based on this structure (see Figure 3).

**Figure 3.** A top-view and a side-view of the computational supercell used in our calculations as a model of TaSi2/NbSi2 nanocomposite (a superlattice) with the stacking along (and the interfaces perpendicular to) the [0001] direction within the C40 lattice (**a**) accompanied with directional dependences of the Young's modulus of bulk TaSi2 (**b**), bulk NbSi2 (**c**) and their composite TaSi2/NbSi2 (**d**).

The nanocomposites have the interfaces perpendicular to the [0001] crystallographic direction. The lattice parameters *a* and *c* and elastic constants of both constituents in their bulk phases are summarized in Table 2. Our theoretical values are in an excellent agreement with both experimental data and previous calculations. Both disilicides have all the parameters very similar. The elasticity of bulk phases is conveniently visualized in Figure 3b,c in the form of directional dependences of the Young's modulus. As the elastic properties are quite similar, the elasticity of the studied nanocomposite is not too different from that of the constituting phases (see Figure 3d). Importantly, the interface energy is ultra low for this superlattice, again zero within the error-bar of our calculations.

Similar to the case of the WSi2/MoSi2 nanocomposites studied above, we examined how the ultra-low interface energies depend on the molar ratio of the constituting phases as well as on the width of the phases forming the superlattice (see Figure 4). We performed our calculations for a series of six other nanocomposites with different ratio of the TaSi2 and NbSi2 (see Figure 4a–f). Out of the six calculated superlattices shown in Figure 4, those shown in Figure 4c,d have the same ratio of the amount of both materials but a higher number of internal interfaces, six and four, respectively.

When evaluating the interface energies, they were again found to be ultra low, namely zero within the error-bar of our calculations. The lattice parameters *a* and *c* of the C40-based structure are concerned were rather monotonously increasing from their lower values in TaSi2 to higher values in NbSi2. In contrast to this trend, the elastic constants rather monotonously decreased from their higher values in TaSi2 to lower values in NbSi2.



**Figure 4.** Schematic visualization of simulated NbSi2/TaSi2 nanocomposites with the stacking along (and the interfaces perpendicular to) the [0001] direction. The mutual ratio of the amount of constituents varies from 1:5, i.e., (NbSi2)1(TaSi2)5 in the case of (**a**) via 2:4 (**b**), 3:3 (**c**,**d**) to 4:2 (**e**) and 5:1, i.e., (NbSi2)5(TaSi2)1 in the case of (**f**). Variants shown in (**c**,**d**) have an equal amount of both constituting phases (similar to the case of Figure 3a) but a different arrangement of atomic layers. Consequently, there is a higher number of internal interfaces (6 and 4 in the case of (**c**,**d**), respectively) and the layers with different constituents have different widths.

The last systems studied here is that containing two different phases from the binary iron-aluminium system, in particular, an ordered Fe3Al intermetallic compound crystallizing in the D03 structure and a disordered solid solution of 18.75 at.% of Al with a body-centered cubic (bcc) ferromagnetic (FM) matrix, so-called *α*-or B2 phase. The structure of the former was is derived from the bcc lattice and, therefore, both materials structurally match each other. The studied nanocomposite is schematically visualized in Figure 5a. Fe3Al was modeled by a 16-atom conventional cell depicted in the upper part of in Figure 5a. The *α*-phase represents a challenge for quantum-mechanical calculations because it lacks any long-range periodicity. We used a 16-atom supercell with atoms distributed according to so-called special quasi-random structure (SQS) concept.

**Figure 5.** Visualization of a Fe3Al/Fe-Al nanocomposite (a superlattice) with the stacking along (and the interfaces perpendicular to) the [001] direction (**a**) accompanied with directional dependences of the Young's modulus of bulk Fe3Al (**b**), bulk Fe-Al (**c**) and the nanocomposite Fe3Al/Fe-Al (**d**) formed out of them. Parts (**b**–**d**) were visualized by the SC-EMA [79–81] library (scema.mpie.de) based on ab initio computed elastic constants.

The nanocomposite was formed by stacking the Fe3Al on top of the *α*-phase along the [001] direction (the interfaces are perpendicular to this direction). In contrast to the above-discussed cases of superlattices formed by pairs of ordered transition-metal disilicides, when the both interfaces in the simulation supercell were identical, the supercells modeling the Fe-Al-based nanocomposites have different interfaces due to different distribution of atoms in the disordered *α*-phase. The interface energies were then averaged values related to the two interfaces. As another difference between the pairs of structurally and elastically nearly identical disilicides discussed above, the two constituting phases have clearly distinguishable elastic properties. Again, they are conveniently visualized in the form of directional dependences of the Young's modulus in Figure 5b,c for the Fe3Al and the Fe-Al *α*-phase, respectively. The calculated values of elastic constants for Fe3Al compound are *C*<sup>11</sup> = 211 GPa, *C*<sup>12</sup> = 161 GPa and *C*<sup>44</sup> = 139 GPa. The elastic constants calculated for the disordered Fe-Al *α*-phase were projected onto a set of elastic constants possessing a cubic symmetry according to the rigorous mathematical theory by Moakher and Norris [87]. Similar concepts are often used in case of systems with any form of disorder (see, e.g., Refs. [59,88–91]). The resulting cubic-symmetry elastic constants are *C*<sup>11</sup> = 217 GPa, *C*<sup>12</sup> = 131 GPa and *C*<sup>44</sup> = 120 GPa. Both phases exhibit 001 directions as elastically soft and 111 directions as elastically hard (i.e., with the minimum and maximum values of the Young's modulus, respectively). The Fe3Al is also apparently elastically more anisotropic. The overall elasticity of their composite is then shown in Figure 5d.

**Figure 6.** Schematic visualizations of different computed Fe3Al/Fe-Al nanocomposites. The computed variants shown in sub-figures (**a**–**f**) differ by the distribution of atoms in the disordered Fe-Al phase.

Importantly, the interface energy was found to be ultra low again, only 0.005 J/m2, which represents an energy difference appearing in Equation ((1)) smaller than 0.001 eV/atom, i.e., within the error-bar of our calculations. To determine an impact of distribution of atoms in the disordered Fe-Al *α*-phase on the interface energies, we performed a series of six additional calculations for supercells which have the same stoichiometry but differ in distributions of atoms in the Fe-Al *α*-phase (see Figure 6). In fact, the structure of coordination shells of atoms remain the same. For example, considering the Al atoms, their distribution in the part corresponding to the Fe-Al *α*-phase in the structural variants in Figure 6 are the same but the part of the supercell corresponding to the Fe-Al *α*-phase is either rotated and/or the atomic planes are permuted. As far as the latter process is concerned, if the atomic planes perpendicular to the [001] direction in the *α*-phase part of Figure 5a would be numbered 1, 2, 3, 4, then by a permutation is meant, e.g., an arrangement 2, 3, 4, 1. Importantly, if being in a single-phase bulk, such permutation of atomic planes within the periodically repeated cell or rotations of the whole cell would not change the energy because the position of the origin

of coordinates (and the attached coordinate frame) can be arbitrarily shifted with respect to the crystal lattice.

**Table 3.** Calculated elastic constants of Fe3Al/Fe-Al nanocomposites with the stacking along (and the interfaces perpendicular to) the [001] direction. The nanocomposites are visualized in Figure 6.


The situation is, on the other hand, different in nanocomposites because the interface represents a reference point not existing in the single-phase bulk. The calculated energies are covering a broader range: 0.055 J/m2 (Figure 6a), 0.021 J/m2 (Figure 6b), 0.032 J/m2 (Figure 6c), 0.006 J/m2 (Figure 6d), 0.139 J/m2 (Figure 6e), and 0.014 J/m<sup>2</sup> (Figure 6f). While this sensitivity on the atomic distribution is rather significant, the elastic properties of nanocomposites shown in Figure 6 are very similar. The computed values of elastic constants are summarized in Table 3.

As both phases appearing in the studied nanocomposites are magnetic, we further examined local magnetic moments of the iron atoms in configurations visualized in Figure 6a–f. The magnitude of local magnetic moments are shown in Figure 7 by the diameter of the spheres representing the Fe atoms. The lowest and the highest value (1.8 μ<sup>B</sup> and 2.44 μB) are explicitly mentioned in Figure 7a to demonstrate a scaling, by which the magnitude of the local magnetic moment is indicated by the diameter of the spheres. Importantly, the magnetic properties of the Fe atoms turned out to be very sensitive to the distribution of atoms (they are reduced when the Al atoms are nearby).

**Figure 7.** Calculated local magnetic moments of iron atoms (the magnitude is indicated by the diameter of the spheres representing the Fe atoms). The local magnetic moments shown in sub-figures (**a**–**f**) correspond to atomic distributions of nanocomposites visualized in sub-figures (**a**–**f**) of Figure 6, respectively.

#### **4. Discussion**

The above-discussed ultra low interface energies in the MoSi2/WSi2 nanocomposites may also indicate that both constituents are prone to mixing even at the atomic level. Indeed, it seems that longer annealing times lead to formation of solid-solution phases rather than (nano)composites. For example, MoSi2/WSi2 powders are synthesized by means of self-propagating high temperature combustion in [92] but solid solutions of MoSi2/WSi2 and Mo5Si3/W5Si3 are found. In [93], it is also reported that it is hard to distinguish MoSi2 and WSi2 phases and (W,Mo)Si2 mainly in solid solution is found in Ref. [94]. The difficulties to distinguish MoSi2 and WSi2 Bragg peaks can be attributed to tetragonal MoSi2 and WSi2 phase having a long-range structure with very similar lattice constants (*a* is equal to 0.3202 nm and 0.3211, *c* amounts to 0.7855nm and 0.7835nm, respectively) [83]. It is also confirmed that MoSi2/WSi2 solid solution powder with nanometric (Mo,W)Si2 structure forms via combustion synthesis method from the mechanical activated powder mixture [95].

To test the scenario of a random solid solution of Mo and W within a C11*<sup>b</sup>* lattice, we performed a series of calculations for supercells modeling these states (see Figure 8a–d). The corresponding enthalpies of mixing (evaluated with respect to the energy of MoSi2 and WSi2 as reference end-members) are shown in Figure 8e and all of them are between zero and −0.001 eV/atom, i.e., within the error-bar of our calculations and comparable to the energy differences obtained when simulating the MoSi2/WSi2 nanocomposites.

**Figure 8.** Schematic visualization of the supercells modeling random solid solutions of Mo and W within a C11*<sup>b</sup>* lattice in the case of Mo:W ratio equal to 1:7 (**a**), 2:6 (**b**), 3:5 (**c**), and 4:4 (**d**) together with the correspondingenthalpies of mixing (**e**). The supercells for the Mo:W ratios equal to 5:3, 6:2 and 7:1 were obtained by swapping Mo and W atoms in the supercells shown in sub-figures (**a**–**c**). The atoms in the 32-atom supercells (2 × 2 × 1 multiple of 6-atom conventional cell of the C11*<sup>b</sup>* structure) were distributed according to the special quasi-random structure (SQS) concept [96].

The above-discussed competition between formation of two-phase nanocomposites on the one hand and single-phase solid solutions on the other hand probably explains why a suitable preparation route is still being searched for in the case of TaSi2/NbSi2 nanocomposites when Nb solubility in TaSi2 extremely large [97]. Our results related to the TaSi2/NbSi2 nanocomposites are intended as a motivation for future studies of this interesting system.

#### **5. Conclusions**

We performed a first-principles study of structural, thermodynamic and elastic properties of nanocomposites exhibiting ultra low or low interface energies. As examples of systems with predominantly covalent interatomic bonds, we studied two combinations of transition-metal disilicides: (i) MoSi2/WSi2 nanocomposites with individual constituents crystallizing in the tetragonal C11*<sup>b</sup>* structure; and (ii) TaSi2/NbSi2 with the two components crystallizing in the hexagonal C40 structure. The constituents within each pair of materials exhibit very similar structural and elastic properties and we obtained ultra low (nearly zero) interface energy for their nanocomposites (within the error bar of our calculations, i.e., about 0.005 J/m2). The interface energy was found to be nearly independent on the width of individual constituents within the nanocomposites and/or crystallographic orientation of the interfaces.

As an example of a magnetic system, a pair of metallic phases containing from the Fe-Al system with different atomic ordering was considered. In particular, we simulated coherent superlattices formed by an ordered Fe3Al intermetallic compound and a disordered Fe-Al phase with 18.75 at.% Al, the *α*-phase. Both constituents are structurally and elastically rather similar (but less than the two pairs of studied disilicides). To estimate the interface energy in the nanocomposite containing the disordered *α*-phase, which lacks a long-range periodicity, we simulated seven different distributions of atoms in the *α*-phase interfacing the Fe3Al intermetallic compound. The resulting interface energies were again either ultra low or low, from 0.005 to 0.139 J/m2. While the impact of atomic distribution on the elastic properties was found insignificant, the local magnetic moments of the iron atoms sensitively depended on the type and the distribution of surrounding atoms.

**Author Contributions:** Conceptualization, M.F.; Methodology, D.H.; Resources, M.Š. and M.F.; Writing—Original Draft Preparation, M.F.; Writing—Review and Editing, D.H. and M.Š.; Visualization, M.F.; Supervision, M.Š.; Project Administration, M.Š. and M.F.; and Funding Acquisition, M.Š. and M.F.

**Funding:** The authors acknowledge the Czech Science Foundation for the financial support received under the Project Nos. 16-24711S (M.Š.) and 17-22139S (M.F.). Additional resources were provided by the Academy of Sciences of the Czech Republic through the Fellowship of J. E. Purkynˇe (M.F.) and by the Ministry of Education, Youth and Sports of the Czech Republic under the Project CEITEC 2020, LQ1601 (M.Š). D.H. acknowledges financial support from the Austrian Science Fund (FWF), Project Number P30341-N36. The computational results presented have been achieved in part using the Vienna Scientific Cluster (VSC).

**Acknowledgments:** M.F. and M.Š. acknowledge the support from the Academy of Sciences of the Czech Republic (Institutional Project No. RVO:68081723) and from the Ministry of Education, Youth and Sports of the Czech Republic via the research infrastructure IPMINFRA, LM2015069. Computational resources were made available by the Ministry of Education, Youth and Sports of the Czech Republic under the Projects CESNET (Project No. LM2015042), CERIT-Scientific Cloud (Project No. LM2015085) and IT4Innovations National Supercomputer Center (Project No. LM2015070) within the program Projects of Large Research, Development and Innovations Infrastructures. Figures 1a, 2, 3a, 4, 5a, 6, 8a–d, and A1a,b were visualized using the VESTA package [98].

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**

To analyze the dependence of the ultra-low interface energies in the MoSi2/WSi2 nanocomposites on the crystallographic orientation of the interface, we computed properties of two other superlattices with the interfaces perpendicular to the [010] and [110] directions, respectively. The corresponding supercells modeling these nanocomposites are shown in Figure A1a,b, respectively. The elastic properties of these superlattices are visualized in the form of directional dependences of the Young's modulus in Figure A1c,d. It is obvious that they are again very similar to each other and also very similar to that of the superlattice with the interfaces perpendicular to the [001] direction.

**Figure A1.** Schematics of computational supercells of MoSi2/WSi2 nanocomposites (superlattices) with the stacking along (and the interfaces perpendicular to) the [010] (**a**) and [110] (**b**) directions within the C11*<sup>b</sup>* lattice, respectively. The interface planes are marked by orange color. The calculated elastic constants for these superlattices are shown as directional dependences of the Young's modulus for the composite with the [010] stacking direction (**c**) as well as that with the [110] stacking direction (**d**). Parts (**c**,**d**) were visualized by the SC-EMA [79–81] library (scema.mpie.de) based on ab initio computed elastic constants.

#### **References**


© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

MDPI St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Nanomaterials* Editorial Office E-mail: nanomaterials@mdpi.com www.mdpi.com/journal/nanomaterials

MDPI St. Alban-Anlage 66 4052 Basel Switzerland

Tel: +41 61 683 77 34 Fax: +41 61 302 89 18