<h1>The surface effect on the thermodynamic stability, half-metallic and optical properties of Co2MnGa(001) films: a DFT study</h1>
			<ul class="item-list">
	        	        <li>
	            R. Ashtari Faregh	            	            	            	            <sup aria-label="Affiliated with Department of Physics, Hamedan Branch, Islamic Azad University, Hamedan, IR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            A. Boochani	            	            <abbr title="This is the corresponding author for this article">*</abbr>
	            	            	                <a href="mailto:arash_bch@yahoo.com" class="tiny-icon email-link mx-1" title="Email A. Boochani">
	                    Email
	                </a>
	            	            	            <sup aria-label="Affiliated with Department of Physics, Kermanshah Branch, Islamic Azad University, Kermanshah, IR">
	                2	            </sup>
	            	        </li>
	        	        <li>
	            S. R. Masharian	            	            	            	            <sup aria-label="Affiliated with Department of Physics, Hamedan Branch, Islamic Azad University, Hamedan, IR">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            F. H. Jafarpour	            	            	            	            <sup aria-label="Affiliated with Physics Department, Bu-Ali Sina University, Hamedan, 65174-4161, IR">
	                3	            </sup>
	            	        </li>
	        	    </ul>
	    	    <ol class="affiliations" aria-hidden="true">
	        <li>Department of Physics, Hamedan Branch, Islamic Azad University, Hamedan, IR</li><li>Department of Physics, Kermanshah Branch, Islamic Azad University, Kermanshah, IR</li><li>Physics Department, Bu-Ali Sina University, Hamedan, 65174-4161, IR</li>	    </ol>
	    
<h2>Abstract</h2>
<p>In this study, the half-metallic properties, thermodynamic stability and optical parameters of the full-Heusler Co
<sub>2</sub>
MnGa compound and its four different terminations of Co–Co, Co–Mn, Mn–Ga and Co–Ga from the surface of Co
<sub>2</sub>
MnGa (001) have been calculated based on the density functional theory (DFT). The results confirm the ferromagnetic half-metallic behavior with a magnetic moment of 4.08 
<span id="IEq1"><mml:math id="IEq1_Math"><mml:msub><mml:mi>μ</mml:mi><mml:mtext>B</mml:mtext></mml:msub></mml:math><tex-math id="IEq1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \mu_{\text{B}} $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq1.gif"/></span>
 and a gap of 0.32 eV at the Fermi level of Co
<sub>2</sub>
MnGa bulk phase having a Cu
<sub>2</sub>
MnAl-type structure. The density of states curves showed that all possible terminations from the Co
<sub>2</sub>
MnGa (001) surface eliminate the half-metallic behavior except the termination of Mn–Ga case. Moreover, the results indicate that the termination of Mn–Ga with the lowest surface energy is the most stable termination for the application in spintronics. The optical coefficients such as real and imaginary dielectric function, refraction, extinction, energy loss function, optical conductivity and reflections of the bulk and Mn–Ga termination have been calculated and compared.</p><hr/><section><h2>Introduction</h2>
<p>Many of novel aspects of the Heusler alloys in industries such as spintronics, microelectronics and electromechanics [
<a href="#CR1"><sup>1</sup></a>
, 
<a href="#CR2"><sup>2</sup></a>
–
<a href="#CR3"><sup>3</sup></a>
] have revealed the necessity of producing these compounds in the films form. The growth of Heusler films is a serious challenge in research due to the emergence of the structural disorders. Local disorders are major barriers to achieving 100% spin polarization in nanostructures made of layers of Heusler compounds. Hence, the theoretical and experimental study of the half-metallic behavior of the Heusler thin films has been the sketch of many researches [
<a href="#CR4"><sup>4</sup></a>
, 
<a href="#CR5"><sup>5</sup></a>
–
<a href="#CR6"><sup>6</sup></a>
]. Yet, controlling the structural disorders is still an important issue to achieve applicable films with unique properties, just like samples with perfect order.</p>
<p>The half-Heusler compound NiMnSb with half-metallic behavior was discovered in about 1983 [
<a href="#CR7"><sup>7</sup></a>
]. However, cobalt-based full-Heusler compounds showing full spin polarization have only been known for a few years. Full-Heusler alloys with 
<em>X</em>
<sub>2</sub>
<em>YZ</em>
 chemical formula and various physical properties, depending on the composition, might be metal, semiconductor or semi-metal, but many of them are ferromagnetic. High-spin polarized Heusler ferromagnetisms seem to be ideal for applications in giant magneto-resistance (GMR) and tunneling magneto-resistance in magnetic tunneling junctions (MTJs). [
<a href="#CR8"><sup>8</sup></a>
] Among these, Co
<sub>2</sub>
MnAl films may be introduced as ideal candidates for high-sensitive Hall sensors due to the unusual inherent Hall conductance [
<a href="#CR9"><sup>9</sup></a>
, 
<a href="#CR10"><sup>10</sup></a>
–
<a href="#CR11"><sup>11</sup></a>
]. So far, reports have been presented on the study of layered structure samples of Co
<sub>2</sub>
YZ (
<em>Y</em>
 = Mn, Cr, Fe, 
<em>Z</em>
 = Al, Si, Ga). For example, the effects of structural order on the transport and magnetic properties of Co
<sub>2</sub>
MnGa have been investigated, as well as the process of growth of these films with amorphous structures, and also their structural order improvement. [
<a href="#CR12"><sup>12</sup></a>
, 
<a href="#CR13"><sup>13</sup></a>
, 
<a href="#CR14"><sup>14</sup></a>
–
<a href="#CR15"><sup>15</sup></a>
]. In this group of Heusler compounds, Manganese is the main element, because one can say that the most known ferromagnetic Heuslers, except CrO
<sub>2</sub>
, are Mn-based compounds. Recent studies on Co
<sub>2</sub>
MnX (X = Si, Ge, Sn) and their films by Raphael et al. [
<a href="#CR5"><sup>5</sup></a>
, 
<a href="#CR6"><sup>6</sup></a>
, 
<a href="#CR16"><sup>16</sup></a>
, 
<a href="#CR17"><sup>17</sup></a>
–
<a href="#CR18"><sup>18</sup></a>
] predicted that the thin films of these compounds are not half-metallic ferromagnetic, or that their full spin polarization is limited by antisite defects [
<a href="#CR19"><sup>19</sup></a>
, 
<a href="#CR20"><sup>20</sup></a>
, 
<a href="#CR21"><sup>21</sup></a>
, 
<a href="#CR22"><sup>22</sup></a>
, 
<a href="#CR23"><sup>23</sup></a>
, 
<a href="#CR24"><sup>24</sup></a>
–
<a href="#CR25"><sup>25</sup></a>
]. Hence, in this work, the half-metallic behavior of Co
<sub>2</sub>
MnGa thin films is discussed, and the results of the study of the electrical and optical properties and thermodynamic stability of the 4 different terminations from Co
<sub>2</sub>
MnGa (001) surface are presented as well as a comparison with its bulk state. In Sect. 
<a href="#Sec2"><sup>2</sup></a>
, we express the calculation methods in details. Sections 
<a href="#Sec4"><sup>3.1</sup></a>
–
<a href="#Sec8"><sup>3.3</sup></a>
 discuss the properties of the bulk state of the Heusler compound Co
<sub>2</sub>
MnGa, its (001) surface, and the optical properties of the bulk and surface cases, respectively; and finally, Sect. 
<a href="#Sec9"><sup>4</sup></a>
 is the conclusions.</p></section>
<section><h2>Calculation methods and models</h2>
<p>The calculations were performed based on the density functional theory (DFT) implemented in Wien2K computational package, using the electronic full potential linear augmented plane waves method (FP-LAPW) [
<a href="#CR26"><sup>26</sup></a>
, 
<a href="#CR27"><sup>27</sup></a>
], and the generalized gradient approximation (GGA) was applied to obtain the exchange–correlation interaction potential [
<a href="#CR28"><sup>28</sup></a>
]. The energy and forces on each atom were converged with 1 × 10
<sup>4</sup>
 eV/atom and 0.02 eV/
<span id="IEq2"><mml:math id="IEq2_Math"><mml:mi mathvariant="normal">Å</mml:mi></mml:math><tex-math id="IEq2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ {\AA} $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq2.gif"/></span>
. Also, the optimized cutoff energy, RK
<sub>max</sub>
, lmax and KPoint were selected as − 6, 8 Ry, 10 and 3000, respectively.</p></section>
<section><h2>Results</h2>
<section><h2>The bulk Heusler Co2MnGa properties</h2>
<p>After optimizing the lattice parameters through energy minimization, the full-Heusler compound Co
<sub>2</sub>
MnGa was crystallized in a cubic structure L2
<sub>1</sub>
 with the Cu
<sub>2</sub>
MnAl-type Fm-3m space group and the equilibrium lattice constant of 5.695 Å. We also found that Co
<sub>2</sub>
MnGa is a ferromagnetic compound with a total magnetic momentum 4.08 
<span id="IEq3"><mml:math id="IEq3_Math"><mml:msub><mml:mi>μ</mml:mi><mml:mtext>B</mml:mtext></mml:msub></mml:math><tex-math id="IEq3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \mu_{\text{B}} $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq3.gif"/></span>
 according to Slider–Pauling’s relation 
<em>μ</em>
<sub>t</sub>
 = 
<em>Z</em>
<sub>t</sub>
<span id="IEq4"><mml:math id="IEq4_Math"><mml:mo>-</mml:mo></mml:math><tex-math id="IEq4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ - $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq4.gif"/></span>
 18 [
<a href="#CR29"><sup>29</sup></a>
, 
<a href="#CR30"><sup>30</sup></a>
]. These values are in agreement with the results of others’ calculations [
<a href="#CR4"><sup>4</sup></a>
, 
<a href="#CR31"><sup>31</sup></a>
].</p>
<p>The primitive L2
<sub>1</sub>
 structure consists of four FCC sub-lattices merged together with the atomic positions as: Co-I (0.25, 0.25, 0.25), Co-II (0.75, 0.75, 0.75), Ta (0.5, 0.5, 0.5) and, Z (0,0,0) depicted in Fig. 
<a href="#Fig1"><sup>1</sup></a>
. Also, the Hg
<sub>2</sub>
CuTi-type structure is another schema in this phase, which consists of four non-equivalent atomic positions with a fourfold symmetry in the space group F4-3m including Co-I (0, 0, 0), Co-II (0.25, 0.25, 0.25), Ta (0.5.0.5.0.5), and Z (0.75, 0.75, 0.75).
<figure id="Fig1"><h3>Fig. 1</h3>
<figcaption><p>The Cu
<sub>2</sub>
MnAl-type crystal structure of Co
<sub>2</sub>
MnGa</p></figcaption>
<img src="40089_2019_289_Fig1_HTML.png" /></figure></p>
<p>To understand the electronic structure of the bulk Co
<sub>2</sub>
MnGa, the total and partial density of states (DOS) as well as the energy band structure are presented in Fig. 
<a href="#Fig2"><sup>2</sup></a>
. From the total DOS curve, it is found that Co
<sub>2</sub>
MnGa can almost be considered as a ferromagnetic half-metal with the spin polarization close to 100% at the Fermi level. According to the band structure curve in both the majority and minority spins, the band-width below the Fermi level from about − 4 to 0 eV is due to the strong hybridization of the Mn-3d and Co-3
<em>d</em>
 electrons in the occupied valence levels. The band structure of the majority spin strongly shows metallic behavior; while in the minority spin, the maximum valence level is almost above the Fermi level which results in a gap of 0.32 eV and approximately semiconducting behavior of this compound.
<figure id="Fig2"><h3>Fig. 2</h3>
<figcaption><p><strong>a</strong>
 The total density of states and band structure curves of the bulk Co
<sub>2</sub>
MnGa. 
<strong>b</strong>
 The partial density of states of Ga atom, and the contribution to the d orbital of Mn and Co atoms in the bulk Co
<sub>2</sub>
MnGa</p></figcaption>
<img src="40089_2019_289_Fig2_HTML.png" /></figure></p>
<p>According to the partial DOS curves, the highest electron contribution at the Fermi level which corresponds to the Slater–Pauling rule belongs to Mn and Co, and these two atoms are the origin of half-metallic behavior at the Fermi level. However, the electron contribution of Ga atoms is connected to the semi-core region, i.e., the range of (− 10 to 7 eV). Also, the localization of the states in the DOS diagram at the conduction region shows that the electronic contribution of the orbital d to the Mn and Co atoms in this range is very high, and the interaction between the 
<em>d</em>
–
<em>d</em>
 bonds of these two atoms and the 
<em>p</em>
–
<em>d</em>
 bonds of Ga and Mn atoms leads to this spin splitting at the Fermi level.</p></section>
<section><h2>Properties of the (001) surfaces</h2>
<p>The Co
<sub>2</sub>
MnGa Heusler structure was cut along the Miller planes (001), and 20 angstroms vacuum were applied in parallel to direction of the crystal growth symmetrically at the two ends of the film. Then, four ideal structures were obtained with Co–Co, Co–Mn, Mn–Ga and Co–Ga terminations containing 9 atomic layers. The lattice constants of these films are 
<span id="IEq5"><mml:math id="IEq5_Math"><mml:mfenced close=")" open="("><mml:mrow><mml:mfrac><mml:mi>a</mml:mi><mml:mn>2</mml:mn></mml:mfrac><mml:mi>√</mml:mi><mml:mn>2</mml:mn></mml:mrow></mml:mfenced></mml:math><tex-math id="IEq5_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \left( {\frac{a}{2}\surd 2} \right) $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq5.gif"/></span>
, where 
<em>a</em>
 = 5.695 Å is the equilibrium lattice constant of the bulk Co
<sub>2</sub>
MnGa. Figure 
<a href="#Fig3"><sup>3</sup></a>
 depicts the side view of these terminations. We relaxed the atomic layers by minimizing the total energy and atomic forces.
<figure id="Fig3"><h3>Fig. 3</h3>
<figcaption><p><strong>a</strong>
–
<strong>d</strong>
 Four different terminations of the (001) surface of Co
<sub>2</sub>
MnGa</p></figcaption>
<img src="40089_2019_289_Fig3_HTML.png" /></figure></p>
<section><h2>Surfaces’ stability</h2>
<p>As we know, surfaces’ stability is an important factor in the growth of thin films. A lower surface energy means that the surface will be more stable [
<a href="#CR32"><sup>32</sup></a>
]. To study the relative stability of the different terminations, we calculated the surface’s free energy as a function of atomic chemical potentials in the framework of the atomic thermodynamics theory. The surface’s free energy is calculated as follows:
<section id="Equ1"><mml:math display="block" id="Equ1_Math"><mml:mrow><mml:mi>γ</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:mn>2</mml:mn><mml:mi>A</mml:mi></mml:mrow></mml:mfrac><mml:mfenced close="]" open="["><mml:mrow><mml:mi>G</mml:mi><mml:mo>-</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>i</mml:mi></mml:munder><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msub><mml:mi>μ</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math><tex-math id="Equ1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \gamma = \frac{1}{2A}\left[ {G - \mathop \sum \limits_{i} \left( {N_{i} \mu_{i} } \right)} \right], $$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_Equ1.gif"/></section>
where 
<em>G</em>
 and 
<em>A</em>
 are the film’s total energy and total surface’s area, respectively. 
<em>N</em>
<sub><em>i</em></sub>
 and 
<em>µ</em>
<sub><em>i</em></sub>
 are the number chemical potential of the 
<em>i</em>
th atom in the sliced layers. The sum of these chemical potentials is equal to the Gibbs free energy of the bulk Co
<sub>2</sub>
MnGa, since it is assumed that the films’ surfaces are in the thermodynamic equilibrium with the middle layers, i.e.,
<section id="Equ2"><mml:math display="block" id="Equ2_Math"><mml:mrow><mml:mn>2</mml:mn><mml:msub><mml:mi>μ</mml:mi><mml:mtext>Co</mml:mtext></mml:msub><mml:mo>+</mml:mo><mml:mi>μ</mml:mi><mml:mi>E</mml:mi><mml:mo>+</mml:mo><mml:msub><mml:mi>μ</mml:mi><mml:mtext>Ga</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:msub><mml:mtext>Co</mml:mtext><mml:mn>2</mml:mn></mml:msub><mml:mtext>MnGa</mml:mtext></mml:mrow></mml:msub></mml:mrow></mml:math><tex-math id="Equ2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ 2\mu_{\text{Co}} + \mu E + \mu_{\text{Ga}} = G_{{{\text{Co}}_{ 2} {\text{MnGa}}}} $$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_Equ2.gif"/></section></p>
<p>By removing 
<em>μ</em>
<sub>Ga</sub>
 from Eq. 
<a href="#Equ2"><sup>2</sup></a>
, we can obtain the surface’s free energy relation in terms of 
<em>μ</em>
<sub>Mn</sub>
 and 
<em>μ</em>
<sub>Co</sub>
. The maximum chemical potentials of the Mn and Co atoms are obtained from their bulk structures, and their minimum chemical potentials are determined in such a way that by decreasing 
<em>μ</em>
<sub>Mn</sub>
 and 
<em>μ</em>
<sub>Co</sub>
, Co and Mn leave the structure and Co
<sub>2</sub>
Ga and MnGa structures are produced.
<section id="Equ3"><mml:math display="block" id="Equ3_Math"><mml:mrow><mml:mfrac><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:msub><mml:mtext>Co</mml:mtext><mml:mn>2</mml:mn></mml:msub><mml:mtext>MnGa</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mtext>MnGa</mml:mtext></mml:msub></mml:mrow></mml:mfenced><mml:mo>≤</mml:mo><mml:msub><mml:mi>μ</mml:mi><mml:mtext>Co</mml:mtext></mml:msub><mml:mo>≤</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mtext>Co</mml:mtext></mml:msub></mml:mrow></mml:math><tex-math id="Equ3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \frac{1}{  2} \left( {G_{{{\text{Co}}_{2} {\text{MnGa}}}} - G_{\text{MnGa}} } \right) \le \mu_{\text{Co}} \le G_{\text{Co}} $$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_Equ3.gif"/></section>
<section id="Equ4"><mml:math display="block" id="Equ4_Math"><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:msub><mml:mtext>Co</mml:mtext><mml:mn>2</mml:mn></mml:msub><mml:mtext>MnGa</mml:mtext></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mrow><mml:msub><mml:mtext>Co</mml:mtext><mml:mn>2</mml:mn></mml:msub><mml:mtext>Ga</mml:mtext></mml:mrow></mml:msub><mml:mo>≤</mml:mo><mml:msub><mml:mi>μ</mml:mi><mml:mtext>Mn</mml:mtext></mml:msub><mml:mo>≤</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mtext>Mn</mml:mtext></mml:msub></mml:mrow></mml:math><tex-math id="Equ4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ G_{{{\text{Co}}_{2} {\text{MnGa}}}} - G_{{{\text{Co}}_{2} {\text{Ga}}}} \le \mu_{\text{Mn}} \le G_{\text{Mn}} $$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_Equ4.gif"/></section></p>
<p>The results of phase diagram show that all of the four terminations are achievable. Since the phase transition point of these terminations is within the allowed range of the phase diagram, it is concluded that the Co
<sub>2</sub>
MnGa bulk structure is also thermodynamically stable. Among the terminations of Co
<sub>2</sub>
MnGa (001) surface, the Mn–Ga termination with the lowest surface free energy is the most stable one, and the Co–Co termination is the most unstable surface with the highest free energy (see Fig. 
<a href="#Fig4"><sup>4</sup></a>
).
<figure id="Fig4"><h3>Fig. 4</h3>
<figcaption><p>The phase diagram of the surface terminations Co
<sub>2</sub>
MnGa (001)</p></figcaption>
<img src="40089_2019_289_Fig4_HTML.png" /></figure></p></section>
<section><h2>The electronic properties of surfaces</h2>
<p>The total density of states curves of Co
<sub>2</sub>
MnGa (001) surfaces with Co–Co, Co–Mn, Mn–Ga and Co–Ga terminations are shown in Fig. 
<a href="#Fig5"><sup>5</sup></a>
a. As can be seen, among the four terminations, only the Mn–Ga termination conserves an approximately half-metallic behavior with a spin polarization of about 94 percent. However, other three terminations of Co–Co, Co–Mn and Co–Ga have spin polarization rates of about 83, 51, and 75%, respectively. Because the distance of the intermediate layers in these films are similar to the bulk system, the atomic density of states of these middle layers matches the corresponding atoms in the bulk structure, and it is observed that at each of the four terminations, the electronic localization below the Fermi level corresponds to the orbits 
<em>d</em>
 of Co, Mn and orbital 
<em>p</em>
 of Ga atom.
<figure id="Fig5"><h3>Fig. 5</h3>
<figcaption><p><strong>a</strong>
 The total density of states curve of the terminations Mn–Ga, Co–Mn, Co–Co, and Co–Ga. 
<strong>b</strong>
 The energy band structure of the Mn–Ga termination</p></figcaption>
<img src="40089_2019_289_Fig5_HTML.png" /></figure></p>
<p>From the band structure that is shown in Fig. 
<a href="#Fig5"><sup>5</sup></a>
b, for the termination Mn–Ga, the strong metallic behavior in spin up, and the partial presence of energy levels above the Fermi level and the reception of a small indirect gap of 0.2 eV in spin down confirm the nearly half-metallic behavior in the bulk state. To more precisely study the surface electronic behavior, the atomic DOS curves of the surface and intermediate layers of this termination are shown in Fig. 
<a href="#Fig6"><sup>6</sup></a>
. The electronic states contribution to the Mn atom in the surface layer far exceed the one in the intermediate layer at the Fermi level, and basically, the surface effects on the Mn atoms in the surface layer have led to this magnetic anisotropy, so that splitting of the states is observed in the minority spin. However, it is further observed that the contribution of Ga to magnetic anisotropy at the Fermi level is very small, and the main density of the states of these atoms belongs to the core region (the energy range from − 7 to − 10 eV). An even more interesting point is that in the valence band in the range of (− 3 to − 5 eV) of the middle and sub-surface Co atoms, the electronic states’ behavior are almost homogeneous and completely non-magnetic. However, in the Fermi region, we see an increase in the electronic states, which is due to the localization of Co electronic states at the film’s surface.
<figure id="Fig6"><h3>Fig. 6</h3>
<figcaption><p>The partial density of states curve of Co, Mn, and Ga atoms in the Mn–Ga termination</p></figcaption>
<img src="40089_2019_289_Fig6_HTML.png" /></figure></p></section></section>
<section><h2>Optical properties of the bulk and surfaces</h2>
<p>In the following, we study the optical parameters for Co
<sub>2</sub>
MnGa bulk and the Mn–Ga termination of Co2MnGa (001) surface in Figs. 
<a href="#Fig7"><sup>7</sup></a>
a–g. Since in the bulk mode, due to symmetry, the crystal is completely homogeneous in different states, so we only plot the 
<em>x</em>
-axis direction.
<figure id="Fig7"><h3>Fig. 7</h3>
<figcaption><p><strong>a</strong>
, 
<strong>b</strong>
 The real and imaginary parts of dielectric function, 
<strong>c</strong>
 the electron energy loss function, 
<strong>d</strong>
 refraction index, 
<strong>e</strong>
, 
<strong>f</strong>
 reflection and extinction index, 
<strong>g</strong>
 the optical conductivity coefficient of the bulk Co
<sub>2</sub>
MnGa and Mn–Ga termination</p></figcaption>
<img src="40089_2019_289_Fig7_HTML.png" /></figure></p>
<p>In the curve of the real part of dielectric function (
<em>ε</em>
(
<em>ω</em>
)) in Fig. 
<a href="#Fig7"><sup>7</sup></a>
a, the similarity of the optical behavior between the bulk state and the Mn–Ga termination is quite clear. The static dielectric function of the bulk Co
<sub>2</sub>
MnGa tends to negative infinity, which is a proof of the metallic behavior of this compound in this interval of energy, while in the Mn–Ga termination along 
<em>x</em>
, it reached a constant value of − 18. Along 
<em>x</em>
-axis, the general form of 
<em>Reε</em>
(
<em>ɷ</em>
) of Mn–Ga termination is similar to the bulk state, except that the response rate of this mono-layer to the incident light is far greater than its bulk state, and it is even more in the areas where the sign of 
<em>ε</em>
<sub><em>x</em></sub>
(
<em>ɷ</em>
) is positive. In the 
<em>z</em>
 direction, the response of Mn–Ga termination to the incident light is quite positive and different from the 
<em>x</em>
-axis. In the region of 0–12 eV, the response of the matter to the light decreases with increasing energy of the incident photon, and there is no distinction between 
<em>x</em>
 and 
<em>z</em>
 directions in 
<em>Reε</em>
(
<em>ɷ</em>
) after 12 eV.</p>
<p>Figure 
<a href="#Fig7"><sup>7</sup></a>
b shows the 
<em>Imε</em>
(
<em>ɷ</em>
) curve for the bulk and Mn–Ga termination of Co
<sub>2</sub>
MnGa. In the bulk state, the static values of this function tend to infinity, which is a solid proof of its metallic behavior and that it has intra-band optical transitions. As seen in the figure, due to the metallic nature of this compound, most of the transitions occur in the infrared and visible regions, and by increasing the energy of the incident photon at the edge of UV, the intensity of transitions decreases, and practically after 12 eV, no transition is observed for this material and it passes light entirely. But in Mn–Ga termination, 
<em>Imε</em>
(
<em>ɷ</em>
) is completely anisotropic in two directions 
<em>x</em>
 and 
<em>z</em>
 at low energies, so that in the 
<em>x</em>
 direction, its behavior is similar to the bulk state with lesser intensities of transitions in the energies of 1.3, 2.2 eV; while in the 
<em>z</em>
 direction, some peaks are emerged with increasing the energy of incident photon at 0.9 eV, 1.5 eV, and 2.5 eV. In both directions, 
<em>Imε</em>
(
<em>ɷ</em>
) tends to zero at high energies. In Fig. 
<a href="#Fig7"><sup>7</sup></a>
c, we see the energy loss function (
<em>L</em>
(
<span id="IEq6"><mml:math id="IEq6_Math"><mml:mi>ω</mml:mi></mml:math><tex-math id="IEq6_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \omega $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq6.gif"/></span>
)) curve, which represents the amount of energy loss of the incident photon. For the bulk Co
<sub>2</sub>
MnGa, all optical regions practically do not show a significant drop in the incident energy other than the energy of about 22 eV, which has a very sharp Dirac peak. In the curve of the loss function for the Mn–Ga termination, it is observed that the Dirac peak of the bulk state is destroyed and in both directions the curve has a broadening. In both 
<em>x</em>
 and 
<em>z</em>
 directions, more than one peak is observed with lower intensities than the bulk Co
<sub>2</sub>
MnGa, and these peaks represent the roots of the 
<em>Reε</em>
(
<em>ɷ</em>
) in the 
<em>x</em>
 and 
<em>z</em>
 directions.</p>
<p>We consider the refractive index in Fig. 
<a href="#Fig7"><sup>7</sup></a>
d. For the bulk Co
<sub>2</sub>
MnGa, according to our expectation of metallic behavior of this compound, the static refractive index tends to large quantities, because in these energies and in the infrared region, electrons of this metal can easily be excited by incident light. By increasing the energy of the incident photon and entering the visible region, the refractive index decreases. Moreover, in the refractive index curve of the Mn–Ga termination, we see that in low energies, the refractive index is anisotropic in the two directions of 
<em>x</em>
, 
<em>z</em>
. It is worth considering the above issues of the metallic behavior of the Mn–Ga termination along 
<em>x</em>
, which is also clear here, so that, in low energies, the high refractive index confirms this metallic behavior, but with increasing intensity of the incident photon, the refraction of light decreases in both 
<em>x</em>
 and 
<em>z</em>
 directions. In the 
<em>z</em>
 direction, the static refractive index is about 6, which suggests the semiconductivity behavior of the compound in this direction. After 10 eV energy, the refraction indexes are lower than one that indicated the super luminance phenomena.</p>
<p>The reflection index curve is observed in Fig. 
<a href="#Fig7"><sup>7</sup></a>
e, the general form of which for Mn–Ga termination with less intensity is similar to the bulk state. The static value of the reflection index for the bulk Co
<sub>2</sub>
MnGa and the Mn–Ga termination along 
<em>x</em>
 direction is about 95%, i.e., a very strong metal in this direction. After about 10 eV, the reflection coefficient is shifted toward long peaks; so in these energies, most of the incident light is reflected.</p>
<p>The extinction coefficient is a measure of light absorption in matter. In Fig. 
<a href="#Fig7"><sup>7</sup></a>
f, the extinction coefficient curves of bulk Co
<sub>2</sub>
MnGa and its Mn–Ga termination similarly show that the maximum extinction coefficient occurs in the infrared region and at zero energy limit, which is due to the metallic nature of the matter. However, with the increase in the energy of the incident photon, the amount of extinction coefficient intensely decreases to the range of 1 eV, then continue decreasing but with a lower rate. As shown in the reflection coefficient curve, the reflectance in this region is about 50%. At the energy level of 3.2 eV, An characteristic peak for both the bulk and Mn–Ga termination is observed along 
<em>x</em>
, which is equivalent to the root of 
<em>Reε</em>
(
<em>ɷ</em>
) in this energy. Moreover, in the 
<em>z</em>
 direction, for the Mn–Ga termination, we see three peaks in the visible area and the UV edge (at 0.8 eV, 2.1 eV and 3.1 eV, respectively). After 9 eV, the curves of the extinction coefficient are symmetric in both directions, overlap and tend to zero. In Fig. 
<a href="#Fig7"><sup>7</sup></a>
g, it is observed that the optical conductivity coefficient has a very high value, due to the occurrence of metallic behavior in the bulk and the Mn–Ga termination along 
<em>x</em>
 at zero energy, which is in agreement with the graph of 
<em>Reε</em>
(
<em>ɷ</em>
). In the energy range of (1.5–10 eV), simultaneously with the absorption of incident light, we observe the optical transitions increase in both the bulk and Mn–Ga termination, which results in the highest absorption and optical conductivity in the visible and ultraviolet spectral range. Based on the zero amount of the energy loss function in the lower energies and high 
<em>Imε</em>
(
<em>ɷ</em>
), it is shown that these compositions are good sensitive to the light with high absorption. Also, the zero amount of the 
<em>Imε</em>
(
<em>ɷ</em>
) and extinction coefficients are referred to the transparent behavior for these compounds.</p></section></section>
<section><h2>Conclusion</h2>
<p>The thermodynamic stability, the structural, electronic and optical properties of the four terminations of Co–Co, Co–Mn, Mn–Ga, and Co–Ga from Co
<sub>2</sub>
MnGa (001) surface were studied using the first principles calculations. Calculations are based on DFT framework and GGA approximation. The spin polarized calculations showed that the bulk Co
<sub>2</sub>
MnGa is nearly a half-metal with a spin polarization of less than 100 with a magnetic moment of 4.08 
<span id="IEq7"><mml:math id="IEq7_Math"><mml:msub><mml:mi>μ</mml:mi><mml:mtext>B</mml:mtext></mml:msub></mml:math><tex-math id="IEq7_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ \mu_{\text{B}} $$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2019_289_Article_IEq7.gif"/></span>
; among the all terminations, only the Mn–Ga termination retains this property. The Mn–Ga termination is the most stable thin film due to the lowest surface free energy and the most suitable for maintaining the spin injection property in spintronic devices. Examination of the bulk Co
<sub>2</sub>
MnGa and its Mn–Ga termination reflects the similarity of optical behavior in these two structures.</p></section><hr/><h2>Publisher's Note</h2>
<p>Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
<hr/><h2>References</h2>
<ol><li>Mao et al. (2016) <em>Effect of As and Nb doping on the magnetic properties for quaternary Heusler alloy FeCoZrGe</em> (pp. 1-6) <a href="https://doi.org/10.1016/j.jmmm.2015.08.123" target="_blank">10.1016/j.jmmm.2015.08.123</a></li><li>Xie et al. (2015) <em>First-principles study of four quaternary Heusler alloys ZrMnVZ and ZrCoFeZ (Z = Si, Ge)</em> (pp. 52-55) <a href="https://doi.org/10.1016/j.commatsci.2015.03.010" target="_blank">10.1016/j.commatsci.2015.03.010</a></li><li>Vasil’ev et al. (2003) <em>Shape memory ferromagnets</em> <a href="https://doi.org/10.1070/PU2003v046n06ABEH001339" target="_blank">10.1070/PU2003v046n06ABEH001339</a></li><li>Dubowik et al. (2007) <em>Films of Heusler alloys</em> (pp. 583-598)</li><li>Singh and Barber (2004) <em>Structural, magnetic, and transport properties of thin films of the Heusler alloy Co2MnSi</em> <a href="https://doi.org/10.1063/1.1690868" target="_blank">10.1063/1.1690868</a></li><li>Picozzi and Continenza (2002) <em>Co2MnX.(X = Si, Ge, Sn). Heusler compounds: an ab initio study of their structural, electronic, and magnetic properties at zero and elevated pressure</em> <a href="https://doi.org/10.1103/PhysRevB.66.094421" target="_blank">10.1103/PhysRevB.66.094421</a></li><li>De Groot et al. (1983) <em>New Class of materials: half-metallic ferromagnets</em> <a href="https://doi.org/10.1103/PhysRevLett.50.2024" target="_blank">10.1103/PhysRevLett.50.2024</a></li><li>Oogane (2006) <em>Large tunnel magnetoresistance in magnetic tunnel junctions Co2 MnX (X = Al, Si) Heusler alloys</em> (pp. 834-841) <a href="https://doi.org/10.1088/0022-3727/39/5/S09" target="_blank">10.1088/0022-3727/39/5/S09</a></li><li>Kubler and Felser (2012) <em>Berry curvature and the anomalous Hall effect in Heusler compounds</em> <a href="https://doi.org/10.1103/PhysRevB.85.012405" target="_blank">10.1103/PhysRevB.85.012405</a></li><li>Tung and Guo (2013) <em>High spin polarization of the anomalous Hall current in Co-based Heusler compounds</em> <a href="https://doi.org/10.1088/1367-2630/15/3/033014" target="_blank">10.1088/1367-2630/15/3/033014</a></li><li>Zhu and Zhao (2017) <em>Anomalous resistivity upturn in epitaxial L21-Co2 MnAl films</em> <a href="https://doi.org/10.1038/srep42931" target="_blank">10.1038/srep42931</a></li><li>Fujii et al. (1994) <em>Electronic and magnetic properties of X2Mn1-xVxSi (X = Fe and Co)</em> (pp. 1881-1888) <a href="https://doi.org/10.1143/JPSJ.63.1881" target="_blank">10.1143/JPSJ.63.1881</a></li><li>Hütten et al. (2005) <em>Half-metallic alloys</em> (pp. 241-264) Springer-Verlag</li><li>Sakuraba et al. (2006) <em>Direct observation of half-metallic energy gap in Co2 MnSi by tunneling conductance spectroscopy</em> <a href="https://doi.org/10.1063/1.2335583" target="_blank">10.1063/1.2335583</a></li><li>Kudryavtsev et al. (2007) <em>Ferromagnetic resonance in Co2MnGa films with various structural ordering</em> 310(2) (pp. 2271-2273) <a href="https://doi.org/10.1016/j.jmmm.2006.10.726" target="_blank">10.1016/j.jmmm.2006.10.726</a></li><li>Fecher et al. (2005) <em>Design of magnetic materials: the electronic structure of the ordered, doped Heusler compound Co2Cr1−xFexAl</em> <a href="https://doi.org/10.1088/0953-8984/17/46/008" target="_blank">10.1088/0953-8984/17/46/008</a></li><li>Kübler et al. (1983) <em>Formation and coupling of magnetic moments in Heusler alloys</em> <a href="https://doi.org/10.1103/PhysRevB.28.1745" target="_blank">10.1103/PhysRevB.28.1745</a></li><li>Ökoğlu and Gülseren (2010) <em>Electronic structure of half-metallic ferromagnet Co2MnSi at high-pressure</em> (pp. 321-326) <a href="https://doi.org/10.1140/epjb/e2010-00188-0" target="_blank">10.1140/epjb/e2010-00188-0</a></li><li>Boochani et al. (2015) <em>Calculation of half-metal, debye and curie temperatures of Co2VAl compound: first principles study</em> 63(5) <a href="https://doi.org/10.1088/0253-6102/63/5/641" target="_blank">10.1088/0253-6102/63/5/641</a></li><li>Majidi et al. (2017) <em>Effect of point defects on the electronic density states of SnC nanosheets: first-principles calculations</em> (pp. 3209-3215) <a href="https://doi.org/10.1016/j.rinp.2017.08.049" target="_blank">10.1016/j.rinp.2017.08.049</a></li><li>Achour et al. (2017) <em>Plasma surface functionalization of boron nitride nano-sheets</em> (pp. 110-115) <a href="https://doi.org/10.1016/j.diamond.2017.06.012" target="_blank">10.1016/j.diamond.2017.06.012</a></li><li>Solaymani et al. (2018) <em>The effect of different laser irradiation on rugometric and microtopographic features in zirconia ceramics: study of surface statistical metrics</em> (pp. 180-185) <a href="https://doi.org/10.1016/j.jallcom.2018.06.213" target="_blank">10.1016/j.jallcom.2018.06.213</a></li><li>Ţălu et al. (2016) <em>Fractal features of carbon–nickel composite thin films</em> 79(12) (pp. 1208-1213) <a href="https://doi.org/10.1002/jemt.22779" target="_blank">10.1002/jemt.22779</a></li><li>Ţălu et al. (2016) <em>Microstructure and micromorphology of Cu/Co nanoparticles: surface texture analysis</em> 12(5) (pp. 580-588) <a href="https://doi.org/10.1007/s13391-016-6036-y" target="_blank">10.1007/s13391-016-6036-y</a></li><li>Ţălu et al. (2016) <em>Microstructure and tribological properties of FeNPs@a-C: H films by micromorphology analysis and fractal geometry</em> (pp. 164-169) <a href="https://doi.org/10.1016/j.jiec.2016.08.003" target="_blank">10.1016/j.jiec.2016.08.003</a></li><li>Unknown () <em></em></li><li>Kresse and Furthmüller (1996) <em>Efficiency of ab initio total energy calculations for metals and semiconductors using a plane-wave basis set</em> (pp. 15-50) <a href="https://doi.org/10.1016/0927-0256(96)00008-0" target="_blank">10.1016/0927-0256(96)00008-0</a></li><li>Perdew et al. (1996) <em>Generalized gradient approximation made simple</em> <a href="https://doi.org/10.1103/PhysRevLett.77.3865" target="_blank">10.1103/PhysRevLett.77.3865</a></li><li>Skaftouros et al. (2013) <em>Generalized Slater-Pauling rule for the inverse Heusler compounds</em> <a href="https://doi.org/10.1103/PhysRevB.87.024420" target="_blank">10.1103/PhysRevB.87.024420</a></li><li>Jung et al. (2000) <em>Study of the 18-electron band gap and ferromagnetismin semi-Heusler compounds by non-spin-polarized electronic band structure calculations</em> (pp. 113-119) <a href="https://doi.org/10.1016/S0166-1280(00)00483-8" target="_blank">10.1016/S0166-1280(00)00483-8</a></li><li>Sargolzaei et al. (2006) <em>Spin and orbital magnetism in full Heusler alloys: a density functional theory study of Co2YZ (Y = Mn, Fe; Z = Al, Si, Ga, Ge)</em> <a href="https://doi.org/10.1103/PhysRevB.74.224410" target="_blank">10.1103/PhysRevB.74.224410</a></li><li>Han et al. (2012) <em>Preserving stable 100% spin polarization at (111) heterostructures of half-metallic Heusler alloy Co2VGa with semiconductor PbS</em> <a href="https://doi.org/10.1063/1.4759159" target="_blank">10.1063/1.4759159</a></li></ol>