1. Introduction
Heusler alloys belong to a wide family of materials, which can be metals, semiconductors, superconductors, topological insulators, or other classes of compounds [
1]. Discovery of novel Heusler alloys is supported by theoretical modeling [
2,
3]. Heusler alloys have gained significant attention in recent years due to their remarkable properties, including high-temperature stability, half-metallicity, and ferromagnetic shape-memory behavior, and they can be used as functional materials in spin-dependent devices, etc. [
4,
5]. Among the various Heusler alloys, the Mn-based Heusler alloys have been extensively studied due to their potential applications in spintronics and magnetic storage devices [
6]. In a number of theoretical papers, it was shown that the local spin density approximation (LSDA) or GGA schemes are not sufficient for describing the electronic structure correctly [
7]. For transition-metal-based half-metallic compounds, the calculated Hubbard U parameters lie between 1.7 and 3.8 eV, being the smallest for MnAs (Mn-3d orbitals) [
8] and largest for Cr
CoGa (Co-3d orbitals) [
9]. However, Mn-based full Heusler compounds such as Co
MnSi lose their half-metallic behavior if the correlation becomes too strong (above 2 eV) [
10]. A recent study showed that values of Hubbard U and Hund exchange J parameters for Mn ions vary between 0.55 (0.44) eV and 2.60 (0.83) eV [
11]. In particular, in Mn2-based alloys, two ferromagnetic alloys, Mn
NiGa and Mn
NiSn, were found with shape memory properties and Curie temperature values equal to 588 K [
12] and 530 K [
13], respectively. Mn
NbAl and Mn
NiGe were theoretically calculated as half-metallic alloys [
14,
15]. However, these alloys based on Mn2 in many studies were found in a complicated
-Mn structure [
8,
16] or even in the hexagonal structure in the case of Mn
PtAl [
6]. Mn
NiAl is a prototype Mn-based Heusler alloy that exhibits interesting electronic and magnetic properties, making it a promising material for various applications. The electronic structure of Mn
NiAl was studied by first-principles calculations [
17,
18], and Mn
NiAl was found to be a ferrimagnet with the total magnetic moments of 1.48 and 1.14
/f.u. for the austenite and the martensite, respectively. The total moment was estimated to be formed by the antiparallel aligned Mn moments and a small contribution of Ni moment [
17,
18], which was confirmed by recent experimental studies for Mn
NiAl and Mn
Ni
Al
ribbons [
19] and Monte Carlo simulations [
15]. An experimental study [
20] revealed that Mn
NiAl has the inverse Xa-type structure and anomalies in the temperature dependence of the magnetization, which correlate with the electrical resistivity anomalies [
20]. The experimental magnetization Mn
NiAl was found to be as low as 0.2
[
20]; however, theoretical studies substantially overestimate this small experimental value of magnetization of the studied alloy.
In order to understand why MnNiAl has such a small magnetization, in this work we evaluate the electron correlations and determine their effect on the electronic and magnetic structure of the compound under study. Our work sheds light on the role of electron correlations in determining the electronic and magnetic properties of MnNiAl. The results of our study will be useful for designing Mn-based Heusler alloys with desired properties for various applications in spintronics and magnetic storage devices.
2. Materials and Methods
The Mn
NiAl full Heusler alloy has a cubic F
3m structure of symmetry group number 216 in the list of crystallographic groups, with unit cell parameters: a = b = c = 5.795 Å,
= 90°. As shown in
Figure 1, the alloy forms the Xa-type inverse Heusler structure. Atomic positions are Mn1 (0, 0, 0), Mn2 (0.25, 0.25, 0.25), Ni (0.5, 0.5, 0.5), and Al (−0.25, −0.25, −0.25). The Mn atoms are indicated in purple, Ni atoms—gray, and Al atoms—blue. The unit cell of The Mn
NiAl contains two manganese atoms, one nickel atom, and one aluminum atom. The manganese ions Mn1 and Mn2 are in the antiparallel alignment. The crystal structure of Mn
NiAl is plotted in Vesta [
21].
In this work, the electronic structure of the Mn
NiAl alloy was computed using the Quantum ESPRESSO software package [
22]. This software package contains the most common basic and advanced exchange–correlation approximations and methods, as well as an impressive set of post-processing tools. The exchange–correlation potential was employed in a generalized gradient approximation (GGA) of Perdew–Burke–Ernzerhof (PBE) type [
23]. To take into account electron correlations, the DFT+U method [
24] was used with U
= U − J, with the values of Coulomb parameter in the range between 0.9 and 6 eV, and the calculated values for Mn
ions. The Coulomb interaction is not considered for Ni ions. The DFT+U method is an extension of the density functional theory (DFT) that is used to improve the description of strongly correlated electronic systems. In standard DFT, the exchange–correlation functional is approximated as a local density approximation (LDA) or GGA, which can lead to inaccuracies in the description of localized electronic states and strong electron–electron interactions. In turn, the DFT+U method introduces a Hubbard-like term, to the Hamiltonian that describes the onsite Coulomb repulsion U, between electrons in localized orbitals. To calculate the Coulomb interaction parameter for Mn atoms (3d shell), we use constrained DFT calculations in the form of a linear response to a certain potential applied to the system. In order to evaluate the effect of electronic correlations, we make a series of calculations varying the value of the Coulomb interaction parameter in the range from 0.9 to 6 eV. An alternative approach to take into account electron correlations was also used, namely, the hybrid functional HSE. The hybrid functional HSE (Heyd–Scuseria–Ernzerhof) is a density functional theory method that combines the traditional exchange–correlation functionals with a fraction of Hartree–Fock (HF) exchange. In the HSE method, the exchange–correlation functional is split into two parts: a short-range part that is described by LDA or GGA, and a long-range part that is described by a fraction of HF exchange. The weight of the HF exchange is determined empirically or from ab initio calculations. In our calculations, for the correct estimation of the Coulomb interaction parameter, we use the weight of the HF exchange equal to the standard value of 0.3. K-grid and q-grid dimensions consist of 4 × 4 × 4 points. For GGA and DFT+U calculation, the k-grid dimension consists of 8 × 8 × 8 points. The calculations used the standard ultrasoft potentials from the pseudopotential library of Quantum ESPRESSO [
25].
3. Results and Discussions
It is worth noting that the studies were carried out as in the case of the initial ferromagnetic or antiferromagnetic ordering of manganese atoms. In the course of self-consistent calculations, it was found that all starting orderings converged to the ferrimagnetic one. Thus, the antiparallel alignment of the Mn ions was found to be energetically favorable. Therefore, below, we consider it.
In our investigation, the GGA calculation for experimental lattice parameter value 5.795 Å (
Figure 2) was performed first. As can be seen, Mn
NiAl exhibits semimetallic properties in the “minority” spin projection and metallic properties in the “majority” spin projection. Let us consider this in detail (
Figure 3). In the “minority” spin projection, indirect bandgaps are observed in the vicinity of the Fermi energy. The Mn band intersects the Fermi energy in the conduction band, but no bands from the valence band cross it. In the energy range between −4 and −1 eV, the strong hybridization between the Mn and Ni electronic states with the high peaks in the density of states is observed. In the case of the “minority” spin projection, the hybridization between the Mn and Ni states with the high peaks of the density of states is observed too. Mostly, this hybridization and peaks are observed in the valence band, but in the conduction band the narrow peaks of the Mn 3d states are observed at 1.3 eV in the “majority” spin projection and near 0.8 eV. It is important to note that the densities of electronic states for different spin projections differ from each other and are not mirrored. Based on this, it can be concluded that Mn
NiAl is a ferrimagnetic material.
Using constrained DFT calculations, we obtain the Coulomb interaction parameter of Mn atoms (3d shell) to be 0.9 eV. In
Figure 4, the density of states within the framework of the DFT+U method with U = 0.9 eV is presented. With the U parameter included, significant changes are observed. There is a sliding of 3d states of manganese. The Mn and Ni 3d states are strongly localized, forming three peaks in “majority” and “minority” spin projection. It can be observed that hybridization between the Ni and Mn 3d states both in “minority” spin projection on −1.6 eV and −3.6 eV energy and “majority” spin projection on −2.1 eV and −3.2 eV energy decreases.
In
Figure 5, the density of states of 3d manganese orbitals are shown more particularly for the Mn1 and Mn2 ions. One can see that both Mn1 and Mn2 have double and triple degenerate e
and t
orbitals, respectively. A strong hybridization between xy, zx, zy orbitals and between z
and x
-y
is observed. Let us focus on
Figure 5a, where the density of states of 3d Mn1 orbitals is shown. A major part of states is localized between −5 and −1.2 eV in the valance band with the “majority” spin projection and from 0 to 2 eV in the conduction band with the “minority” spin projection. The t
orbitals have several peaks at −4, −3.2, and 1.8 eV energy values, the intensity of which reaches 0.8 eV. In case of the e
orbitals, there are also several peaks localizing in vicinity −3.8, 1.2, and 2.1 eV. The intensity of these peaks is slightly less and more smeared than the t
peaks intensity, and reach the value at about 0.7 eV. Also, there are some hybridized states. Let us consider
Figure 5b, where the density of states of 3d Mn2 orbitals is shown. A major part of states is localized between the same ranges as in the case of Mn1, but only in the opposite spin projections. The e
orbitals have two main peaks localizing at −3.8 eV energy in the valence band with the intensity about 1.6 eV, and on 1.8 eV energy in the conduction band with the intensity about 0.6 eV. The t
orbitals have three pronounced peaks of the density of states. These peaks localize at −4.1, −1.8 eV energies in the valence band with intensities 0.6 and 0.5 eV, respectively, and at 1.8 eV energy in the conduction band with the intensity about 1.2 eV. There are also some hybridized states. It is worth noting that the value of peaks intensity of Mn2 is significantly more than the value of peaks intensity of Mn1.
With the increase of U to 3 eV (
Figure 6), the 3d-band become less hybridized Ni and Mn. The Mn 3d states shifted to the low energy levels in the valence band and high energy levels in the conduction band. The Ni 3d states are localized in the vicinity of the −2 eV energy level. In the case of “minority” spin projection, manganese 3d states localize in the valence band in the vicinity of the −5 eV energy level, and in the conduction band in the vicinity of the 1.4 eV energy level. For “majority” spin projection, the same situation is observed except for the fact that the Mn 3d states peak in the conduction band are shifted to the higher energy relative to the Mn 3d states of “majority” spin projection.
After increasing the Coulomb interaction parameter to 6 eV (
Figure 7), a significant displacement of the manganese 3d states to the lower energy is observed in the valence band in each spin projection, but shifting of the manganese 3d states in the conduction band is not obvious and constitutes only 1 eV in the direction of the higher energies. There is also a strict localization of the manganese 3d states. For the Mn1 states, the peaks are localized in the ranges of energies between −8 and −6.75 eV in the valence band and from 1 to 3.5 eV in the conduction band. For the Mn2 states, the peaks are localized in ranges of energies from −8 to −6.5 eV in the valence band and from 2 to 4 eV in the conduction band. It is also worth noting the splitting of the peaks of the 3d states of manganese in the valence band. The Ni 3d states do not shift and localize in the same range of energies between −3.5 and −1 eV in the valence band. Despite these changes, the compound under study retains the metallic properties.
Figure 8 demonstrates the density of the electronic states of the compound Mn
NiAl in the framework of the HSE method. It can be seen that there are numerous localized maximum peaks throughout the figure, but if we compare the obtained HSE graph with previous DFT+U results, we find that the obtained HSE curve is matched with the DFT+U one for U = 3 eV.
Let us now focus on
Table 1 and
Figure 9, where the partial and total magnetic moments are listed. In the case of GGA approximation calculation, the Mn magnetic moments are −2.83 and 3.12
, which correspond to antiparallel ferrimagnetic. The total magnetization is 1.10
and does not match with experimental value. Considering the Coulomb parameter and using the DFT+U method, Mn magnetic moments become equal to −3.26 and 3.60
when we include calculated U = 0.9 eV, which is much more than the case of calculation with the GGA approximation. The total magnetization is reduced by half and equal to 0.54
, which is consistent with the experimental data [
20] and in contrast with the previous DFT calculations [
17,
18]. This decreasing occurs by reduction of the difference between the Mn magnetic moments and reduction of the value of the Ni magnetic moment. When increasing the Coulomb parameter to 3 eV, the difference between Mn magnetic moments continues to decrease. The Mn magnetic moments consist of −3.75
and 3.95
and Ni magnetic moments decrease up to 0.07
, respectively. As a result, the total magnetization reduces to 0.35
. When the Coulomb parameter is 6 eV, the Mn magnetic moments increase up to −4.15 and 4.25
, Ni magnetic moments reduce to 0.00
, and the total magnetization consists of 0.18
. In the case of HSE calculation, the Mn magnetic moments consist of −3.74 and 3.87
, Ni magnetic moments consist of 0.09
, and the total magnetization consists of 0.25
. It also agrees well with the experimental value of total magnetization.
In our study, we reveal that the experimental lattice parameter does not correspond to the equilibrium minimum total energy of the ground state. Varying the lattice parameter within the range of 5.6 and 6.6 Å, we found that the more energetically favorable lattice parameter is 6.1 Å (
Figure 10). Below, we present the results of calculations for Mn
NiAl with lattice parameter corresponding to 6.1 Å and compare the results obtained.
As in the case of calculation with GGA for Mn
NiAl with the experimental lattice parameter, let us consider the electronic and band structure of Mn
NiAl with the energetically favorable lattice parameter in the case of GGA (
Figure 11). It can be seen that peaks of density of electronic states on the energy level of −1.5 eV in all spin projections become discernible. This happened as a result of the removal of a strong d-band hybridization between Mn and Ni (
Figure 12). In the case of “majority” spin projection, shifting of the 3d-Mn states to the Fermi level is observed. The main peaks of density of Ni 3d states are localized at −1.3 eV and −1.5 eV energy levels in the valence band. The main peaks of the density of Mn1 3d states are localized at the −2.5 eV energy level in the valence band and on 1 eV energy level in the conduction band. The main peaks of density of Mn1 3d states are localized at the −3 eV energy level in the valence band and on 1.75 eV energy level in the conduction band. In general, the type of electronic structure is comparable with the case of the experimental lattice parameter with U = 0.9 eV. It is important to note changes in band structure. Now in the case of “minority” spin projection, the Mn band intersects the Fermi level both in the conduction band and the valence band. However, the intersection between the conduction band and the valence band is still not observed. In the case of “majority” spin projection, shifting of the 3d-Mn states to the Fermi level is observed. In spite of this, Mn
NiAl still exhibits metallic properties.
Similar to the case of the experimental value of the lattice, we calculated the Coulomb interaction parameter using the constrained DFT calculations, and its value is equal to 1 eV. The obtained results are shown in
Figure 13. The behavior of Mn
NiAl with the lattice parameter of 6.1 Å is similar to the behavior of Mn
NiAl with the lattice parameter of 5.795 Å. There is a decrease in hybridization between 3d states of Ni and Mn. The main peaks of the density of Mn1 3d states localize in the vicinity of the −3.5 eV energy level in the valence band and in the vicinity of the 1.5 eV energy level in the conduction band. In turn, the main peaks of the density of Mn2 3d states localize in the similar range of energy levels in the valence band, but the peak of the density of Mn2 3d states in the conduction band is shifted relative to another peak of the density of Mn1 3d states towards higher energies, and is localized at the 2.3 eV energy level. However, unlike Mn
NiAl with the lattice parameter of 5.795 eV, the peak of hybridized 3d states of Ni and Mn is observed on the Fermi level. It is also observed that the weakly intense peaks of 3d-Mn states localize in the vicinity of the Fermi level.
With the subsequent increase of the Coulomb parameter to 3 eV (
Figure 14), a shifting of the electronic states of manganese in proportion to the increase in the Coulomb parameter occurs. However, the main peaks of density of Ni 3d states are localized at the same energy and their intensity does not change. The main peaks of the density of Mn 3d states shift to the −5 eV energy level in the valence band. In the conduction band, there is a shift of 0.25 eV towards high energies, and also the decrease of the intensity of the density of Mn1 3d states occurs. The system under study reacts to a change in the Coulomb parameter in the same way as a system with the experimental value of the lattice constant. The compound under the study retains the metallic properties.
After increasing the Coulomb interaction parameter up to 6 eV (
Figure 15), displacement of the manganese 3d states to the lower energy occurs in the valence band in each spin projection. Their intensity increases to 6.2 eV due to localization. There is also a splitting of peaks. The peak of the density of Mn1 3d states in the conduction band splits into two peaks with intensity of about 2.3 eV. The peak of the density of Mn2 3d states shifts to the 3.5 energy level, and its intensity is about 6 eV. It is worth noting that the overall intensity of the peaks decreases due to the removal of hybridization between the 3d states of manganese and nickel. Despite these changes, the compound Mn
NiAl retains the metallic properties.
Similar to the HSE calculation presented above for Mn
NiAl with a lattice constant of 5.795 Å, a figure of the density of the electronic states of the Mn
NiAl alloy with a lattice constant of 6.1 Å was obtained (
Figure 16). As can be seen from the figure below, the density of electronic states corresponds to the density of electronic states in the case of DFT+U calculation with the Coulomb parameter of 3 eV.
Unlike the electronic structure, there is a significant change in the magnetic properties (
Table 2). The values of the magnetic moments of manganese are noticeably greater (
Figure 17). Even when calculated with the GGA approximation, the values of the magnetic moments of manganese are −3.24 and 3.60
, but their change with an increase in the Coulomb parameter remains the same as for Mn
NiAl with the lattice parameter of 5.795 Å. With the growth of the Coulomb parameter, the magnetic moments of the manganese increase, reaching −4.43 and 4.46
, respectively, and the difference between them decreases. The behavior of the nickel magnetic moment differs from the previous one. With the growth of the Coulomb parameter, its value decreases, and when the Coulomb parameter reaches 6 eV, the magnetic moment of nickel changes its orientation and takes the value −0.04
. The total magnetization of the unit cell at the lattice parameter of 6.1 Å is slightly lower. Without taking into account the Coulomb parameter, its value is 0.59
. With an increase in the Coulomb parameter, the total magnetization decreases, reaching a value of 0.02
. In the case of the HSE calculation, the obtained values of the total magnetization agree well with the experimental data [
20].