Co 2 CrAl Heuslerene: Mechanical, Thermodynamic and Electronic Properties

: This work investigates the ground state’s stability of the bulk and three Heuslerene Co 2 CrAl compounds, named as α , β , and γ phases, by density functional theory (DFT) with the generalized gradient approximation (GGA), GGA+U, and GGA+U+mBJ approximations. The results demonstrate the ground state stability of all mentioned cases since they pass the thermodynamic, elastic, and phonon stability tests. All three structures are more stable in the ferromagnetic phase than the antiferromagnetic phase. In the β phase, Young’s and Shear’s moduli were 73.97 GPa and 24.83 GPa, respectively. The thermodynamic diagram has shown existence of the accessible region, which indicates the possibility of making this structure. For all three structures, the phonon branches in the symmetry paths are positive, which represent the complete dynamic stability of these compositions in the presence of mechanical stresses and thermal vibrations. According to the electronic calculations, the bulk phase of Co 2 CrAl is a half-metal with 3 µ B magnetic moment and 100% spin polarization at the Fermi level. Furthermore, all imposed approximations approve that α and γ Heuslerenes are metal for both spin directions, while the GGA+U+mBJ approximation indicates that β phase is a ferromagnetic half-metal of 1 µ B magnetic moment. Based on the electron density diagrams, the highest (lowest) amount of electron density is created on the α ( β ) phase surface.


Introduction
The unique properties of spin devices based on Heusler alloys motivated scientists to develop their fundamental knowledge about applying these materials for industrial and research applications.Ferromagnetic materials with Heusler compounds have significant applications in various modern technologies like data storage, sensors, and energy conversion [1,2].Simplicity in manufacturing, high mechanical resistance, low-cast and availability of materials as well as high Curie temperature of Heusler compounds have made them good candidates for using in spintronics, opto-electronic, giant magneto resistance (GMR), tunneling magneto resistance (TMR), magnetic thermoelectric applications, magneto-optic devices [3][4][5][6][7][8][9].
The synthesis of Co-based Heusler compounds began in 1970s [10].Kübler et al. [11] found the disappearance of the minority spin density at the Fermi level of Co 2 MnSn and Co 2 MnAl, resulted from the specific transport of the majority spin density towards meeting the Fermi level.Meanwhile, de Groot et al. [12] introduced half-metals (HM) ferromagnets (HMFMs) with 100% spin polarization at Fermi level, as ideal candidates for spintronic devices with spin injection [13].Band structure (BS) of Heusler compounds has the important role in describing their magnetic characterization and prediction of HMFMs.The first attempts to calculate the BS of some Co-based Heusler compounds (Co 2 MnSn, Co 2 TiSi, and Co 2 TiAl) did not show HM ferromagnetism since the basis of calculations were spherical potentials.Moreover, the simple form of exchange-correlation potential was approximated by local spin density approximation (LSDA) [14][15][16][17][18]. Ishida et al. revealed the first signs of HM ferromagnetism in Co-based Heusler compounds for Ru 2 MnC and Co 2 MnC (C = Si, Al, Sb, and Sn) [19,20].Mohn et al. [21] could only find the magnetic ground state of Co 2 TiC (C = Sn and Al) by full symmetry potentials, not the HM state.Also, Galanakis et al. [22] reported HM behavior in various A 2 BC compounds which was in good agreement with those found by Picozi.et al. for the Mn-based compounds.They used the GGA because, despite the pure LSDA function, the GGA technique considers the gradients of density [23].Verifying the HM state of Co 2 FeAl is not possible by spherical potentials or GGA approximation and for the complete series of Co 2 Cr 1−xx Fe x Al, the HM ferromagnetic ground state can be found when both GGA approximation and full symmetry potentials are applied [24,25].
The ferromagnetic properties with the significant magnetic moment and high Curie temperature were observed in the majority spin of Co-based full Heusler compounds with cubic L 21 structure [26].The Curie temperature and magnetic moment of Co 2 MnSi are 985 K and 4.96 µB, respectively.Furthermore, many of them were predicted to be half-metal and hence, are of particular interest in spintronic devices [27][28][29][30].Recently, the electronic BS and magnetic properties of full Heusler compounds have been investigated on a large scale, both theoretically and experimentally [26][27][28].For instance, the total magnetic moments in these compounds followed the Slater-Pauling-type behavior and were explained by their electronic structures [23].Also, the high Curie temperatures of Co-based Heulser compounds were calculated by ab initio calculations related to the electronic structures [27].
The popularity of two-dimensional (2D) materials with few atoms thickness started in 2004 by discovery of graphene.These materials are considered as promising candidates for the next generation of optoelectronics and electronics devices, because of their mechanical, thermal, and electrical properties [28,31,32].Graphene is a single layer of sp 2 -bonded carbon atoms with high carrier mobility [33] which is almost chemically inert [34].Its mechanical flexibility [35] makes it to be applied in flexible electronic and ultrathin 2D devices.Therefore, graphene has a wide variety of applications in transistors as well as chemical and strain sensors [36,37].For group-IV 2D materials like silicone, graphene, stanene, and germanene, the zero band gap can be opened by the adsorption of special species or vertical electric field.However, their application in transistors is restricted at room temperature because in general, the value of opened band gap is smaller than 0.4 eV [38,39].Recently, materials such as nitrogenene [40], black phosphorene (BP) [41], arsenene [42], and antimonene [43] in group-V 2D have attracted great attention.The direct band gap of monolayer and few-layer BP is almost 0.59-1.51eV [44] and its unique electronic properties like ON/OFF current ratio and high carriers mobility (∼1000 cm 2 /Vs) indicate its necessity in electronic and optoelectronic applications [45][46][47].The discovery of semiconductors and insulators has expanded 2D materials and graphene over the past few years [48,49].Using the 2D structures of the direct band gap semiconductors in transition metal dichalcogenides (TMDs) like monolayer WS 2 and MoS 2 have significantly improved optoelectronics devices [50,51].However, further improvements require continuous investigations of new 2D materials.Current methods for creating 2D monolayers cannot be used in bulk-layered materials since their access to crystals has not been explored yet [52].
The high compactness with mechanical stability of the FCC thin films and 2D hexagonal structures along their [111] crystallographic direction make them a suitable candidate for electronic devices.In our last works, the graphene-like monolayers of Ti Heusler compounds were made along the [111] direction and named Heuslerenes [53,54].According to the DFT-based calculations, the two mentioned Heuslerenes were stable from the static and dynamic viewpoins.In addition, electronic calculations by GGA and mBJ approximations revealed them to be half-metals.The main goal of the present work is to expand the idea of Heuslerenes based on Co 2 CrAl and dealing with their expected properties in some details.
Considering the half-filled d orbitals of two Co and Cr atoms in the Co 2 CrAl compound, we expect to see a half-metallic behavior by cutting the Co 2 CrAl FCC structure and creating 2D compounds named as Co 2 CrAl Heuslerenes.Therefore, for feasibility of making them, at first step, we calculate the formation energy and show that it is possible to form theses combinations.Then we look for the ground state point to show the existence of mechanical stability in the static mode.In the following, stability of these structures under stress and strain is checked by evaluating their elastic coefficients.At the next step, the mechanical stability will be evaluated from the vibration point of view by means of phonon dispersion diagrams.Furthermore, to complete the possibility of making these compounds in reality, thermodynamic phase diagrams are obtained, the electronic properties of these compounds in the ferromagnetic phase are discussed by drawing the density of states (DOS), BS and electron density diagrams.

Computational Methods
The electronic structure, thermodynamic phase diagram and mechanical stability of Co 2 CrAl in the bulk and Heuslerene phases were studied by density functional theory (DFT) and full potential linearized augmented plane waves (FP-LAPW) method [55] such as our previous research [56].The exchange-correlation potential in the Wien2K [57] code were calculated by generalized gradient approximation (GGA) [58], GGA+U [59,60], and GGA+U+mBJ [61,62] approximations.Our calculations considered by spin polarization mode and Hubbard potential (U) of the Co and Cr atoms were optimized as 0.235 Ry and 0.0, respectively.The optimized R mt K max parameter is 8 and the Brillouin zone KPoint for the bulk and Heuslerens are 17 × 17 × 17 and 19 × 19 × 5, respectively.The muffin tin radius of Co and Cr atoms is 2.15 a.u. and for Al is 2.00 a.u.Elastic constants of the bulk and Heuslerene crystals were proposed by Jamal [56,63,64] considering the Hubbard U parameter (GGA+U) in the Wien2K framework.In order to show the crystal dynamical stability, phonon dispersions have been extracted from the Quantum Espresso package [65,66].For increasing the accuracy of phonon calculations, and other electronic results the atomic forces were relaxed up to the best optimization of 10 −6 (a.u./dyn) and the optimized k-point and q-point meshes for phonon calculations were selected as 24 × 24 × 3 and 10 × 10 × 2, respectively.To calculate the mechanical, thermodynamic phase diagram, and electronic properties, we used from the 40 core super computer in six month.As mentioned above, the layers cut along its [111] crystallographic direction must be hexagonal, dense enough, and stable.Thus, considering the atomic termination of the slices, we have produced three ultrathin films of about 8 Bohr thick, named as α, β and γ phases (see Figure 1) with p3m1 space group.It is important to note that the stoichiometry of each product is different from the main source, as Co 2 CrAl 2 , Co 2 Cr 2 Al, and Co 3 CrAl for the α, β, and γ shapes, respectively.We will see that this extra atom in each monolayer plays a positive role in its stability.In order to maintain a system, the first step is to find equilibrium volume at minimum energy.The total energy of each system has been calculated by fitting it to the Brich-Mornaghan equation for different unit cell volumes and the results are summarized in energy-volume (E-V) plots of Figure 2. In these calculations, after finding the equilibrium volume and the ground state point by E-V curves, the atomic forces has been optimized up to 10 −6 (a.u./dyn).

Results and Discussion
Panels (a)-(d) of this figure indicate symmetric outline with a minimum point which confirm the static stability of the all inspected forms.Comparing panels (a) and (b) illustrates similarity between the E-V curves of the bulk state and α phase.Also, β and γ phases are similar.In the β phase, the slope of the curve increases with decreasing volume (under stress) which indicates more resistivity of this plane, while for other two planes, the curves are relatively symmetric on both sides of the equilibrium point.More curvature and slope in the E-V curve means that bulk modulus of the composition is large.Among the Heuslerne phases, β phase has the greatest E-V curve slope and thus the highest amount of bulk modulus.Due to the presence of magnetic atoms, Co and Cr, in this compound, the E-V curves are calculated in both ferromagnetic and anti-ferromagnetic phases.It is evidently concluded from Figure 2 that all Heuslerne phases are more stable in the ferromagnetic state and hence, we continue the following calculations just based on ferromagnetism.
Resulted data of the E-V curves are listed in Table 1.Clearly, the lattice constant of Co 2 VAl is in agreement with other reports.The magnetic moment of this compound is integer 3, in accordance with Slater-Pauli rule, implying that this compound is a suitable option for spintronic purposes.Besides, the bulk modulus (hardness) of the bulk phase of Co 2 VAl is twice that of steel metal, suggesting its potential application for industrial tools and parts.Derivative of the bulk modulus coefficient of this material unveils whose atomic bonds are of ionic type, resulting in a high crystalline strength.Negative values of the formation energy (E form ) for all phases in the ferromagnetic (FM) and anti-ferromagnetic (AFM) modes show that they are all quite stable from the energy point of view.The following equation is used to calculate the E form : Minimum points approve that all structures are statically stable.

Elastic Stabilities
Elasticity refers to the retrieval of the original shape of materials after removing applied force.Hence, the amount of elastic constants can specify structural stability and can create a connection between the atomic and large scales in order to predict important parameters such as hardness, brittleness/ductility nature, melting temperature, type of  The lattice constants of the three Heuslerene phases are very close to each other, indicating their structural similarity, although they are different in the arrangement of their atoms, as shown in Figure 1.The Co 2 VAl-β Heuslerene magnetic moment is an integer, 1, implying that it could be used in the spintronic industry.The bulk modulus of these 2D structures show that, due to their low thickness (a few atomic layers), they have suitable degree hardness.Also, their bulk modulus derivative verifies the existence of strong ionic bonds between their atoms.The magnetic moment in the β and γ phases have shown correct values, the partial magnetic moments have also shown that in the γ structure, we are facing a significant increase in the magnetic property of atoms and this two-dimensional structure.

Elastic Stabilities
Elasticity refers to the retrieval of the original shape of materials after removing applied force.Hence, the amount of elastic constants can specify structural stability and can create a connection between the atomic and large scales in order to predict important parameters such as hardness, brittleness/ductility nature, melting temperature, type of bonds, stiffness, and etc [67].These parameters can provide valuable information about the practical applications.
The theoretical background of the elastic constants along with the calculation of parameters such as bulk modulus, shear modulus, Debye temperature, and . . .for bulk materials with different symmetry have been reported in our previous reports [63,64,[67][68][69].Moreover, as represented in our recent work [56], the theoretical background of the elastic constants for 2D materials as well as their appropriate amounts for the Heuslerene Ti 2 VGe were investigated by using IRelast2D package [56,70] which worked as an interface with Wien2k code [71].It should be noted that the experimental investigation of the elastic constants of 2D materials is a challenge for scientists at present [72] and so, such a package (IRelast2D) might be the best alternative which opens a new theoretical way for investigating the mechanical properties of 2D systems.
Mechanical parameters revealed the fact that when the percentage of B/S ratio of bulk material is less/more than 175%, it is brittle/ductile [73].The Cauchy pressure sign is defined as (C 12-C 44 ) and specifies the sort of bonds.Commonly, the dominance of covalent or ionic bonds makes the sign of the bulk material's Cauchy pressure to be negative or positive, respectively [74].Moreover, Poisson's ratio determines the bond type, so that for materials with covalent (ionic) bonds it is less than 0.25 (close to 0.25 or more) [68].
Stiffness is the resistivity of the matter in reaction to an applied force so that a stiff material has a higher Young's modulus and a flexible material has less stiffness.Hence, higher values of Young's modulus are essential when deflection is disagreeable.Hardness, as another important parameter, describes the resistivity of the matter to the shape variations.Among three parameters that indicate hardness i.e., shear modulus, bulk modulus, and Vickers hardness, better results can be achieved by characterizing the Vickers hardness [68] which can be calculated by following equations provided by Chen [75] and Tian [76]: , S : shear, B : bulkmodulus) In the case of some ionic compounds for which Chen's equation detects negative values, Tian's equation reformulates it to solve this problem [76].To evaluate thermodynamic parameters like Debye temperature (θ D ) as well as transversal (V t ), longitudinal (V l ), average (V m ) sound velocities, the Hill scheme is applied as an intermediate to the Voigt and Reuss model [67].Table 2 includes the elastic constants, mechanical and thermodynamic parameters of the bulk Co 2 CrAl in comparison with the previous theoretical reports [77,78] using the same methods and exchange-correlation (FP-LAPW-PBE) [79,80].Our ab initio results are in agreement with previous theoretical calculations.The value of Poisson's ratio affirms that the dominant bonds are ionic which is compatible with our deduction from the positive sign of Cauchy pressure.The percentage of B H S H ratio for the bulk shape is 214% from which the ductile nature is inferred.The values of Vickers Hardness, based on Chen's and Tian's equations, are positive and almost close to each other.The elastic stability criterion (known as Born stability conditions) emerges when elastic energy of the compound is positive.From the mathematical standpoint, it is held when eigenvalues of the elastic constants matrix (C) are positive [81,82].For cubic symmetries, we have: Our results show that bulk Co 2 CrAl is mechanically stable and passes the elastic stability test.
We should recall that, as far as we know, no experimental or theoretical reports have been published about the 2D Co 2 CrAl yet.The directional values of Young's moduli in Table 3 show that three Heuslerene phases of Co 2 CrAl are isotropic along the x and y directions.Also, their Young's moduli are less than bulk form, meaning that all these Heuslerenes are of less stiffness than bulk.However, the α-phase has greater amounts of Y and thus, is stiffer than β and γ phases.Shear modulus (ratio of shear stress to shear strain), which is the representative of hardness, describes the resistance of the material to the shape changes.Hence, Co 2 CrAl-Al is harder than other two Heuslerenes.The hardness decreases when we move from the bulk towards the surface.

Phonon Stability
Having proved the static stability of the bulk and three Heuslerene shapes of Co 2 CrAl, it is essential to investigate their dynamical stability as well.For this aim, the phonon BSs in Figure 3a-d, which concerns the crystal vibrations over the first Brillouin zone, is analyzed.All frequency branches in these plots are positive implying that all discussed phases of Co 2 CrAl are dynamically stable.For the bulk phase, steeper acoustical branches around Γ point infer the transformation of thermal energy in the 0-200 cm −1 region and so, it may act as a good heat conductor along this direction.However, in the range of 200-270 cm −1 310-390 cm −1 two big optical gaps are observed and their optical branches have low slope and are nearly flat, so that the IR energy with low mobility accumulates in this region.
observed and their optical branches have low slope and are nearly flat, so that the IR energy with low mobility accumulates in this region.
Moreover, Figure 3b-d shows that all three 2D predicted forms of Co2CrAl are strongly stable against the lattice vibrations.These figures are plotted along the -Γ K-M-Γ path within the first Brillouin zone.The acoustical branches of these Heuslerenes have a higher slope than optical branches.This is particularly observed in the α case where the slope of acoustical branches around the Γ point is more and thus, this structure can better convey the low-frequency thermal energies.There are three phonon band gaps for the β case and one phonon band gap for the γ case observed in the range of frequencies more than 240 cm −1 and 340-400 cm −1 , respectively.Because of qualitative similarity between the phonon BSs of the bulk and β phases (such as the ranges of optical gaps), more similarity in their optical behavior is expected.Further, among these 2D structures, the highest energy level (~400 cm −1 ) belongs to the γ case which is still a bit lower than that of the bulk shape.This means that any change in the type and arrangement of the comprising atoms may result in changes in the energy levels and thus optical behavior (including IR conductivity) of the ultimate compound.Moreover, Figure 3b-d shows that all three 2D predicted forms of Co 2 CrAl are strongly stable against the lattice vibrations.These figures are plotted along the -ΓK-M-Γ path within the first Brillouin zone.The acoustical branches of these Heuslerenes have a higher slope than optical branches.This is particularly observed in the α case where the slope of acoustical branches around the Γ point is more and thus, this structure can better convey the low-frequency thermal energies.There are three phonon band gaps for the β case and one phonon band gap for the γ case observed in the range of frequencies more than 240 cm −1 and 340-400 cm −1 , respectively.Because of qualitative similarity between the phonon BSs of the bulk and β phases (such as the ranges of optical gaps), more similarity in their optical behavior is expected.Further, among these 2D structures, the highest energy level (~400 cm −1 ) belongs to the γ case which is still a bit lower than that of the bulk shape.This means that any change in the type and arrangement of the comprising atoms may result in changes in the energy levels and thus optical behavior (including IR conductivity) of the ultimate compound.

Phase Stability Diagrams Stability
The phase diagram of any compound evaluates all possible formations by considering the Gibbs energies.Figure 4 displays the phase diagram of the bulk Co 2 CrAl based on the following formula for its Gibbs energy: where µ Co , µ Cr , and µ Al are chemical potentials of the comprising atoms.The accessible area which satisfies this equation is considered as the allowed region and its borders are determined by following equations: in these relations, ∆µ i is the chemical potential change of atom i and their ideal amounts are ∆µ i ≥ 0 which means that changes are assumed to be between their maximum and minimum values.The minimum value of µ i is found where atom i departs from the compound and its maximum value refers to crystallization of atom i in its bulk form with the Gibbs energy of g bulk

Electronic Properties
The DOS curves of the above-mentioned phases of Co2CrAl for both up and dow spin distributions are illustrated in Figure 5a-d.Panel (a) shows that in the bulk form, th Each corner in the triangle of Figure 4a corresponds to a point where one of the atomic components of Co 2 CrAl leaves the structure.In other words, each corner fits the chemical potential of a specific type of atom.Therefore, the chemical potential of Al atom has been fixed and over the whole range of changes in other chemical potentials, the free Gibbs energies for a variety of possible combinations of Co, Cr, and Al have been plotted.Some triangles may emerge from the intersection of these lines but one of them, which is marked as a hachured area, illustrates the region where the three atoms can be simultaneously present in Co 2 CrAl with space group of Fm − 3m.The phase diagrams for all three Heuslerene phases are compared is Figure 4b by calculating the surface free energies of these 2D structures through the following equation: In which, G slab is the surface Gibbs energy.µ i and N i are chemical potential and number of atoms in supercell.P i and T are also the partial pressure and temperature of the components, respectively.For our 2D systems, this equation can be rephrased as below: In which, the factor 2 in denominator of the right side is due to the two similar surfaces on both sides of the layer, and A is the surface area.The accessible region in this case can be found from Equations ( 4)- (6).All 2D shapes of Co 2 CrAl in Figure 4b lie into the allowed area of stability, with belonging the biggest and smallest area to the γ and α phase, respectively.Therefore, all these Heuslerenes of Co 2 CrAl are likely to be synthesized in the lab under suitable conditions.

Electronic Properties
The DOS curves of the above-mentioned phases of Co 2 CrAl for both up and down spin distributions are illustrated in Figure 5a-d.Panel (a) shows that in the bulk form, the spin-up electronic states have overlapped over the whole positive and negative energies and cut the Fermi level, suggesting a strong metallic nature.However, for the spin-down channel, a p-type semiconductor outline is observed whose values of the energy gap and spin-flip gap are 0.64 eV and 0.14 eV, respectively.In half-metallic compounds and in minority spins, the distance from the Fermi level to the edge of electronic states is called spin-flip gap, and is obtained based on Kohn-Sham equations as follows [90]: In which ε L↑ (ε L↓ ) is the upper (lower) energy of unoccupied orbitals and ε H↑ (ε H↓ ) is the upper (lower) energy of occupied orbitals.
Consequently, the bulk phase of Co 2 CrAl is a half-metal with the magnetic moment of 3 µB and 100% spin polarization at the Fermi level, which makes it applicable in spintronics and spin injection applications.Plotted in panels (b-f) of Figure 5 are the electronic DOS of the 2D phases for the majority and minority spins using GGA, GGA+U, and GGA+U+mBJ approximations.It is evident that DOS curves of all Heuslerene phases, in comparison to the bulk, are more serrated, owing to their nano-scale entity and the important role of free electrons on their surfaces.The α and γ phases seem to be metal for both spin directions in all imposed approximations, however, they exhibit the most anisotropy by GGA+U+mBJ.Salient is the β phase which exhibits different magnetic properties for different approximations.For the region of −3 eV to 2 eV, it is magnetically anisotropic in GGA, isotropic in GGA+U, but extremely anisotropic at the Fermi level in GGA+U+mBJ.In the last case, the β phase is a p-type semiconductor with 0.27 eV energy gap for spin-up.But for spin-down, a strong metal with a continuous spectrum of states is observed.In Figure 5, the partial DOS curve of the β phase with GGA+U+mBJ approximation is plotted for both up and down spins.It is observable that below and above the Fermi level, d orbital of the Cr atom has the highest value and is the main reason for conduction and transport of the electrons and holes.Overlap of the d orbitals of Cr and Co with p orbitals of Al in the energy range of −0.5 eV to −3.5 eV gives rise to the d-d and p-d bonds.
For more details, the electronic BS of the aforementioned phases of Co 2 CrAl have been calculated and plotted in Figure 6.Clearly, for minority spins, the bulk phase is a p-type semiconductor with the direct gap of 0.64 eV around the Γ point, which is compatible with the DOS results and other reports.However, for majority spins, high slope of the conduction bands with energies more than 1 eV, is a testament to the enormous mobility of the related electrons, which is due to the relationship of group velocity E g = 1 ∂E ∂K and to the first derivative and the second derivative of energy relative to the momentum vector.For β case, two spin-up electronic bands meet the Fermi level, as expected, and an indirect band gap appears within the spin-down bands.Here, the maximum valance band and minimum conduction band are placed at Γ and K points, respectively.The large gradient of the electronic states is remarkable in the conduction domain, from which a huge degree of electronic mobility can be inferred.The difference between the BS profiles of the up and down spin channels means that the thermodynamic behavior of this Heuslerene compound can be adjusted by an external field.
GGA+U+mBJ.In the last case, the β phase is a p-type semiconductor with 0.27 eV gap for spin-up.But for spin-down, a strong metal with a continuous spectrum is observed.In Figure 5, the partial DOS curve of the β phase with GGA approximation is plotted for both up and down spins.It is observable that bel above the Fermi level, d orbital of the Cr atom has the highest value and is th reason for conduction and transport of the electrons and holes.Overlap of the d of Cr and Co with p orbitals of Al in the energy range of −0.5 eV to −3.5 eV gives ri d-d and p-d bonds.
For more details, the electronic BS of the aforementioned phases of Co2Cr been calculated and plotted in Figure 6.Clearly, for minority spins, the bulk ph p-type semiconductor with the direct gap of 0.64 eV around the Γ point, w compatible with the DOS results and other reports.However, for majority spi slope of the conduction bands with energies more than 1 eV, is a testamen enormous mobility of the related electrons, which is due to the relationship o velocity  = ℏ and effective mass  * = ℏ to the first derivative second derivative of energy relative to the momentum vector.For β case, two electronic bands meet the Fermi level, as expected, and an indirect band gap within the spin-down bands.Here, the maximum valance band and m conduction band are placed at Γ and K points, respectively.The large gradien electronic states is remarkable in the conduction domain, from which a huge d electronic mobility can be inferred.The difference between the BS profiles of the down spin channels means that the thermodynamic behavior of this Heu compound can be adjusted by an external field.In Figure 7, electron density diagrams are plotted on the Heuslerene Co 2 CrAl plates.It is clear that Al, Cr and Co atoms are present at the surface of the α, β and γ phases, respectively.The presence of electron density and dangling bonds on the surface of films and 2D structures gives rise to new physical properties in them.According to the figures, the highest presence of electron density occurs at the α phase, and the lowest at β case.The magnitude of electron density in the α phase relative to β one is in the order of ten.Accumulation of electron density at the α and γ phases have caused them to exhibit strong metallic properties.On the other hand, presence of Cr atoms in the surface of the β phase causes more polarization of the electron density than other two phases.Besides, this electron density is due to half-full d-orbitals that have given a magnetic behavior to the electronic properties of this phase.As we saw in Table 1, the magnetic moment of this phase is the integer 1.0 µB.Therefore, the presence of more free electrons in the two-dimensional structure of α and γ phases have increased the metallic behavior in them, and in the β phase, we saw full half-metallic behavior with controllable magnetic behavior.The magnitude of electron density in the α phase relative to β one is in the order o Accumulation of electron density at the α and γ phases have caused them to ex strong metallic properties.On the other hand, presence of Cr atoms in the surface of phase causes more polarization of the electron density than other two phases.Be this electron density is due to half-full d-orbitals that have given a magnetic behav the electronic properties of this phase.As we saw in Table 1, the magnetic moment o phase is the integer 1.0 μB.Therefore, the presence of more free electrons in two-dimensional structure of α and γ phases have increased the metallic behav them, and in the β phase, we saw full half-metallic behavior with controllable mag behavior.In Figure 7, electron density diagrams are plotted on the Heuslerene Co2CrAl plates.It is clear that Al, Cr and Co atoms are present at the surface of the α, β and γ phases, respectively.The presence of electron density and dangling bonds on the surface of films and 2D structures gives rise to new physical properties in them.According to the figures, the highest presence of electron density occurs at the α phase, and the lowest at β case.The magnitude of electron density in the α phase relative to β one is in the order of ten.Accumulation of electron density at the α and γ phases have caused them to exhibit strong metallic properties.On the other hand, presence of Cr atoms in the surface of the β phase causes more polarization of the electron density than other two phases.Besides, this electron density is due to half-full d-orbitals that have given a magnetic behavior to the electronic properties of this phase.As we saw in Table 1, the magnetic moment of this phase is the integer 1.0 μB.Therefore, the presence of more free electrons in the two-dimensional structure of α and γ phases have increased the metallic behavior in them, and in the β phase, we saw full half-metallic behavior with controllable magnetic behavior.

Conclusions
The F.C.C structure of the Heusleres provides a compact arrangement of atoms along their [111] direction and hence, the thin films and 2D structures extracted by cutting the bulk form in this direction may be of high compactness and stability.In this work, Heuslerenes were expanded to the novel 2D structures by introducing a new member based on Co 2 CrAl and dealing with its expected properties in details.The crystal elastic constants for the bulk structure are obtained using the method proposed by Jamal in the Wien2K framework.
The structural stability was estimated by amounts of the elastic constants.Also, thermodynamic characteristics were evaluated by Hill scheme as an intermediate to the Voigt and Reuss evaluations.Our ab initio results are in good agreement with previous theoretical calculations.The value of Poisson's ratio confirmed the ionic dominant bonds which is compatible with our deduction from the positive sign of Cauchy pressure.Besides, the ductile nature was inferred from the percentage of B H S H ratio in bulk Co 2 CrAl which was 214%.The values of Vickers Hardness, based on Chen's and Tian's equations, were positive and almost close to each other.The directional values of Young's moduli revealed the isotropic nature of all three Heuslerene phases of Co 2 CrAl along the x and y directions.Moreover, while moving from bulk scale towards the sheet, hardness decreased.Comparing to other Co 2 CrAl Heuslerene phases, the elastic constants (C 11 and C 12 and C 66 ) of Co 2 CrAl-Al were larger and hence, it was mechanically stronger.The positive value of all frequency branches of phonon distribution diagrams implied that all phases of Co 2 CrAl were dynamically stable.Also, all 2D predicted forms of Co 2 CrAl were strongly stable against the lattice vibrations.
According to our phase diagram calculation for the bulk shape, a hachured area illustrated the region where all three atoms could be simultaneously placed in Co 2 CrAl with the space group of Fm − 3m.All 2D phases of Co 2 CrAl lied into the allowed area of stability and the biggest and smallest area belonged to the γ and α phase, respectively.
Moreover, the half-metallic nature of the bulk phase was demonstrated by extracting its electronic structure.Besides, the value of magnetic moment and the percentage of spin polarization at the Fermi level were 3 µ B and 100%, respectively, meaning that it could be employed for spintronics and spin injection applications.α and γ phases seemed to be metal in both spin directions and all approximations but the most anisotropy was exhibited in the GGA+U+mBJ approximation.Finally, β phase was noticeable for exhibiting different magnetic properties in different approximations, such that in the −3 eV to 2 eV region, it was magnetically anisotropic in GGA, but isotropic in GGA+U, and extremely anisotropic at the Fermi level in GGA+U+mBJ.Also, it was a p-type semiconductor with a 0.27 eV energy gap for spin-up, but a strong metal with a continuous spectrum of states for spin-down.As a prospect, mechanical exfoliation, as the most popular methods to extract 2D materials from bulk crystals, may be applied to bring this Heuslerene structure into existence.

3 . 1 .
Stability of the Bulk and 2D Sheets of Co 2 CrAl 3.1.1.Structural Stability Co 2 CrAl is a Co-based full-Heusler compound crystallized in F.C.C. structure with space group of Fm − 3m.

Figure 1 .
Figure 1.Crystal structure of the Co2CrAl compound from top and side view for: (a) α-, (b) β-and, (c) γ-phases (d) shows the thickness of the atomic layers from the side.

Figure 1 .
Figure 1.Crystal structure of the Co 2 CrAl compound from top and side view for: (a) α-, (b) βand, (c) γphases (d) shows the thickness of the atomic layers from the side.

i.Figure 4 .
Figure 4.The stability phase diagrams of Co2CrAl in the (a) bulk, and (b) Heuslerene phase Chemical potential changes are expressed in Rydberg unit.

Figure 4 .
Figure 4.The stability phase diagrams of Co 2 CrAl in the (a) bulk, and (b) Heuslerene phases.Chemical potential changes are expressed in Rydberg unit.

Figure 6 .
Figure 6.Electronic BSs for two spin distributions of just two phases of Co 2 CrAl: (a) Bulk phase by GGA (b) β phase by GGA, (c) β phase by GGA+U, and (d) β phase by GGA+U+mbJ.

Figure 6 .
Figure 6.Electronic BSs for two spin distributions of just two phases of Co2CrAl: (a) Bulk phase by GGA (b) β phase by GGA, (c) β phase by GGA+U, and (d) β phase by GGA+U+mbJ.

Figure 7 .
Figure 7. Three Heuslerene monolayers of Co2CrAl, named α, β and γ.The difference arises from the termination of the surface atoms.

Figure 7 .
Figure 7. Three Heuslerene monolayers of Co 2 CrAl, named α, β and γ.The difference arises from the termination of the surface atoms.

Table 1 .
Structural and mechanical properties of Co 2 CrAl compound in the bulk and three Heuslerene phases, including: crystal lattice constants (a, c), formation energy (E for ), total magnetic moment (M tot ), partial magnetic moment (M p ), Bulk modulus (B), and derivative bulk modulus (B ' ).Al are the total energy of the bulk or Heuslerene unit cell, and energies of the Co, Cr, and Al crystals, respectively.
total , E Co , E Cr , and E

Table 2 .
The values of elastic constants (C ij in (GPa)), Bulk moduli (B), shear moduli (S), and Young's moduli (Y) in GPa in the Voigt-Reuss-Hill (VRH) approximation, Poisson's ratios (ν), as well as transverse (V t ), longitudinal (V L ) and average wave velocities (V m in m/s), Debye temperature (Θ D in (K)) and Vickers hardness base on the Chen and Tian equations (H V in (GPa)) for the bulk form of Co 2 CrAl.

Table 3 .
TheAs another result, all Poisson's ratios of the Heuslerenes are positive, with largest value for the γ-phase, and independent of the x and y directions.With reference to Table3, the largest and smallest amounts of the elastic parameters, including C 11 , C 22 , C 66 , Y x , Y y , and S xy belong to the α and γ phase, respectively, while for the Poisson's ratio an inverse behavior is observed.Obviously, Co 2 CrAl-Al is mechanically stronger than other two Heuslerene phases, since its elastic constants (C 11 , C 12 and C 66 ) are larger.From the viewpoint of elasticity, all phases pass the elastic stability test and hence, they are mechanically stable.
values of elastic constants (C ij in (N/m)), shear moduli (S), and Young's moduli (Y) in N/m and Poisson's ratios (ν) for three Heuslerene phases of Co 2 CrAl.