crystals

Article DFT Investigation on the Electronic, Magnetic, Mechanical Properties and Strain Effects of the Quaternary Compound Cu2FeSnS4

Jing Bai †, Jiaying Ji †, Liyu Hao †, Tie Yang and Xingwen Tan * School of Physical Science and Technology, Southwest University, Chongqing 400715, China; [email protected] (J.B.); [email protected] (J.J.); [email protected] (L.H.); [email protected] (T.Y.) * Correspondence: [email protected] Authors contributed equally. †  Received: 27 May 2020; Accepted: 10 June 2020; Published: 15 June 2020 

Abstract: The electronic, magnetic and mechanical properties of the quaternary compound Cu2FeSnS4 have been investigated with first principle calculations. Its half-metallicity has been identified with spin polarized band structures and its magnetic origination is caused by the strong spin splitting effect in the d orbitals of Fe atoms. The total magnetic moment of 4 µB is mainly contributed by the Fe atoms and the spatial distribution of the magnetic spin density and charge density difference have also been examined. Moreover, several mechanical properties of Cu2FeSnS4 have been derived and its mechanical stability is also verified. The directional dependent Young’s modulus exhibits relatively small anisotropy yet the shows strong directional anisotropy. At last, the tetragonal strain effects have been evaluated and their impact on the electronic and magnetic properties are provided. Results show the total magnetic moment stays almost unchanged while the half-metallicity can only be maintained under relatively small variations for both strains. This study can provide comprehensive information about the various properties of Cu2FeSnS4 compound and serve as a helpful reference for its future applications.

Keywords: first principles calculation; electronic band structure; half-metallicity; strain condition; quaternary compounds

1. Introduction During recent years, the family of half-metallic ferromagnets have aroused strong research interest in the field of material science and condensed matter physics, especially for and magnetoelectronics [1–6]. Half-metallicity is characterized by the spin polarized band structures with a band gap in one spin direction, similar to , and no gap in the other spin direction, similar to metals. A great deal of effort has been dedicated to this material family, both experimental characterizations and theoretical calculations, and, except the intrinsic half-metallicity, various other interesting properties have been found, such as thermoelectricity [7–10], [11–13], topological insulativity [14,15], magnetic shape memory [16–21] and so on. Research work is ongoing and devoted to further enhancing their performance or searching for new applications. In different materials, half-metallic ferromagnets are mostly found in ternary and quaternary compounds. For example, many full Heusler and quaternary Heusler materials have been confirmed to show such a property [22–33]. Compared with ternary materials, quaternary ones are more promising because a larger compositional flexibility can be realized [34–44]. A large and extensively studied I II IV VI I II quaternary material group is the A2B C D4 compounds [45–49], in which A is Cu or Ag; B is from

Crystals 2020, 10, 509; doi:10.3390/cryst10060509 www.mdpi.com/journal/crystals Crystals 2020, 10, 509 2 of 12

Mn, Fe, Ni, Zn and Cd; CIV is from IV main group elements Si, Ge and Sn; and DVI is from VI main group elements S, Se and Te. They are originated from tetrahedrally coordinated derivatives of BII-DVI binaries and have mostly diamond-like crystal structures [50]. These materials have recently received a great deal of attention due to their potential applications in thermoelectrics [51], spintronics [52], non-linear optics [53] and photocatalysts [54], and, thus, there has been a notable increase in their syntheses during the past few years. Among this family of compounds, Cu2FeSnS4 has been experimentally synthesized [55,56] and it manifests two tetragonal structure phases of the stannite type with space group I42m and the quasi-cubic type with space group P4. The two structures are very similar but assigned to different space groups because of the modified distributions of the containing metal elements. Although their magnetic susceptibility and optoelectronic properties have been investigated, a basic understanding on their the electronic, magnetic and mechanical properties is still missing. For this reason, we took the first quasi-cubic structure as an example and performed a systematic study with first principles calculation. The presence of half-metallicity in Cu2FeSnS4 was confirmed by its spin polarized band structures and the total magnetic moment of 4 µB was calculated with the main contribution from the 3d Fe atoms. In addition, its mechanical properties were assessed and the corresponding spatial Young’s modulus and shear modulus are obtained. Furthermore, the tetragonal strain conditions were also examined and their effects on the half-metallic and magnetic properties evaluated. It was found the total magnetic moment had almost no change, yet the half-metallicity could only be maintained under relatively small variations for both strains. This detailed work can provide a valuable reference regarding to the various properties for Cu2FeSnS4 and serve as a helpful guide for its future applications.

2. Computational Methodology

The electronic, magnetic and mechanical properties of the quaternary compound Cu2FeSnS4 were investigated with first principle calculations, which were performed with the pseudopotential plane wave methods [57] based on density functional theory (DFT) [58], using CASTEP package [59]. The generalized gradient approximation (GGA) [60] in the framework of the Perdew–Burke–Ernzerhof (PBE) functional [57] was adopted for the electronic exchange correlation energy. The valence electron configurations for Cu, Fe, Sn and S were treated as 3d104s1, 3d64s2, 5s25p2 and 3s23p42, respectively, and their interactions with the atomic core are described by the ultrasoft pseudopotential. Based on an initial convergence test, a cutoff energy of 500 eV was set for the plane wave and a Γ centered 10 10 10 Monkhorst-Pack grid was selected for Brillouin zone k point sampling. The crystal structure × × was fully relaxed until the total force per atom was smaller than 1 10 3 eV/Å. The self-consistent × − filed convergence tolerance was set as the total energy difference per atom less than 1 10 6 eV. × − 3. Results and Discussions

The bulk quaternary compound Cu2FeSnS4 was experimentally synthesized and crystallized with 1 a tetragonal symmetry; space group P4 (S4, No. 81). Its crystal structure is schematically illustrated in Figure1 and the unit cell contained eight atoms in total, including one Sn located at Wycko ff site 1a (0, 0, 0), two Cu at Wyckoff site 2g (0, 0.5, 0.5), four S at Wyckoff site 4h (0.255, 0.255, 0.255) and one Fe at Wyckoff site 1c (0.5, 0.5, 0). In this structure, each S anion was coordinated with four cations and they formed slightly distorted tetrahedrons, as shown by the yellow shaded areas in Figure1. The fully relaxed lattice constants for Cu2FeSnS4 were a = 5.416 Å and c = 5.419 Å, which were consistent with the experimental values very well (a = 5.4186 Å and c = 5.4202 Å) [55]. These optimized lattice constants were used for the following discussions. Crystals 2020, 10x, ,x 509 FOR PEER REVIEW 3 of 12 Crystals 2020, x, x FOR PEER REVIEW 3 of 12

Figure 1. The crystal structures of the quaternary compound Cu2FeSnS4.

Figure 1.1. The crystal structures of the quaternary compoundcompound CuCu22FeSnS4.. Based on the optimized crystal structure for Cu2FeSnS4, we calculated its electronic band Based on the optimized crystal structure for Cu FeSnS , we calculated its electronic band structures. structures.Based Sinceon the the optimized 3d transition crystal metal structure element, for2 Fe, Cu carries24FeSnS a4, largewe calculated magnetic itsmoment, electronic the band structuresstructures.Since the 3 showdSincetransition spin the polarized3 metald transition element, features metal Fe, (see carrieselement, Figure a large Fe, 2). carriesThe magnetic fermi a large moment,energy magnetic is the set bandto moment, zero structures and the only showband the energystructuresspin polarized range show near features spin the polarized Fermi (see Figure level features 2is). shown. The (see fermi InFigure the energy spin-up 2). isThe set direction,fermi to zero energy and a small only is set indirect the to energyzero band and range gaponly near was the observedenergythe Fermi range with level near a isconduction shown.the Fermi In band level the spin-upminimum is shown. direction, (CBM) In the locatedspin-up a small at direction, indirect high symmetry band a small gap point indirect was A observed and band the gap valence with was a bandobservedconduction maximum with band a conduction(VBM) minimum at G. (CBM)band In the minimum locatedspin-down at (CBM) high direct symmetry locatedion, two at point bandshigh symmetry A show and themultiple valencepoint crossingsA bandand the maximum with valence the (VBM)Fermiband maximum energy at G. In level. the (VBM) spin-down This at G.behavior In direction,the spin-down of two bandsdirection, show properties two multiple bands in show crossings one multiple spin with direction crossings the Fermi and with energymetal the propertiesFermilevel. Thisenergy in behavior the level. other ofThis spin semiconductor behavior direction of clearly semiconductor properties indicated in one that properties spin this directionmaterial in one belongs and spin metal directionto the properties family and of inmetal half- the metalpropertiesother spincompounds. directionin the other clearly spin indicateddirection thatclearly this indicated material belongsthat this tomaterial the family belongs of half-metal to the family compounds. of half- metal compounds.

Figure 2. TheThe calculated calculated electronic band structure in both spin directions forfor CuCu22FeSnS4 compound at the equilibrium lattice constants. Figurethe equilibrium 2. The calculated lattice constants. electronic band structure in both spin directions for Cu2FeSnS4 compound at the equilibrium lattice constants. To furtherfurther examine examine the the correlation correlation states states between between different different elements elements and the and the magnetism origination originationfor CuTo2FeSnS further for4 compound, Cuexamine2FeSnS 4the wecompound, havecorrelation also computedwe states have bealso thetween densitiescomputed different of statesthe elements densities and the and distributionsof thestates magnetism and of the magnetic spin density and the charge density difference, as reported in Figure3. The positive and distributionsorigination for of theCu 2magneticFeSnS4 compound, spin density we and have the chargealso computed density difference, the densities as reported of states in Figureand the 3. Thedistributions positive andof the negative magnetic values spin in density the total and and the partial charge densities density difference,of states in asFigure reported 3a,b correspondin Figure 3. The positive and negative values in the total and partial densities of states in Figure 3a,b correspond

Crystals 2020, 10, 509 4 of 12 negative values in the total and partial densities of states in Figure3a,b correspond to the spin-up and spin-down directions, respectively. A spin polarized behavior is observed: zero density at the FermiCrystals level in2020 spin-up, x, x FOR PEER direction REVIEW and nonzero density in spin-down direction, which is in accordance4 of 12 with the findings from the electronic band structures. The total density of states near the Fermi level is mostlyto the contributed spin-up and from spin-down the transition directions, metal respective elementsly. A spin Fe and polarized Cu, especially behavior is in observed: the energy zero range density at the Fermi level in spin-up direction and nonzero density in spin-down direction, which is from 2 eV to +2 eV, and the main group elements, Sn and S, have negligible contributions. However, −in accordance with the findings from the electronic band structures. The total density of states near the twothe transition Fermi level metal is mostly elements contributed exhibit from di thefferent transi densitytion metal distributions: elements Fe and an Cu, almost especially symmetric in the one for Cu element from 2 eV to 0 eV between the two spin directions; a highly asymmetric one for Fe energy range from− −2 eV to +2 eV, and the main group elements, Sn and S, have negligible elementcontributions. with two high However, protrude the two areas transition in spin-up metal directionelements exhibit at 3 different to 2 eV; density and several distributions: high peaksan in − − spin-downalmost direction symmetric at one 0 to1 for eV.Cu Especiallyelement from for −2 the eV densityto 0 eV between at the Fermi the two level spin indirections; spin-down a highly direction, it is completelyasymmetric from one for the Fe Fe element element. with Thistwo high strong protrude spin splittingareas in spin-up effect direction in the density at −3 to of −2 stateseV; and for the Fe elementseveral could high peaks result in in spin-down its strong direction magnetic at 0 moment, to1 eV. Especially which indeed for the density is confirmed at the Fermi by the level calculated in spin-down direction, it is completely from the Fe element. This strong spin splitting effect in the magnetic spin density as shown in Figure3c. The large gray shaded areas around the Fe atoms indicate density of states for the Fe element could result in its strong magnetic moment, which indeed is a strong local magnetic moment. Some shaded areas are also found around Cu atoms but they are confirmed by the calculated magnetic spin density as shown in Figure 3c. The large gray shaded areas much smalleraround the compared Fe atoms withindicate that a aroundstrong local the magnetic Fe atoms. moment. The calculated Some shaded total areas magnetic are also moment found per unit formulaaround isCu 4 atomsµB and but it they is mainly are much contributed smaller compared by the with Fe atom. that around The charge the Fe atoms. density The di calculatedfference is also computedtotal tomagnetic show themoment atomic per bonding unit formula between is 4 μ diB andfferent it is atoms mainly and contributed it is displayed by the Fe in Figureatom. The3d with the graycharge shaded density areas difference representing is also thecomputed electron to accumulation.show the atomic Webonding can clearlybetween see different that the atoms S atoms and form four bondsit is displayed with the in four Figure neighbor 3d with atomsthe gray located shaded at areas the representing corners of the the slightlyelectron distortedaccumulation. tetrahedron We and theycan haveclearly almost see that identical the S atoms electron form four excess bonds with with only the four one neighbor a little bigger atoms betweenlocated at Sthe and corners Fe atoms. of the slightly distorted tetrahedron and they have almost identical electron excess with only one a A very strong dz2 orbital shape is observed at the Fe atoms with the large dumbbell oriented to the little bigger between S and Fe atoms. A very strong dz2 orbital shape is observed at the Fe atoms with centerthe position. large dumbbell oriented to the center position.

FigureFigure 3. The 3. calculatedThe calculated total total (a )(a and) and partial partial ((bb)) densitiesdensities of of states states in inboth both spin spin directions directions and the and the calculated magnetic spin density (c) and charge density difference (d) with isosurface value of 0.1 for calculated magnetic spin density (c) and charge density difference (d) with isosurface value of 0.1 for Cu2FeSnS4 compound at the equilibrium lattice constants. The locations of different atoms in (c) and Cu FeSnS compound at the equilibrium lattice constants. The locations of different atoms in (c,d) can 2 (d) can4 be referred to Figure 1. be referred to Figure1.

Next, we turn to the mechanical properties of Cu2FeSnS4 compound. Since it has a P4 tetragonal Next, we turn to the mechanical properties of Cu FeSnS compound. Since it has a P4 tetragonal structure, there should be seven independent elastic 2constants4 in total: C11, C12, C13, C16, C33, C44 and structure,C66. By there applying should the be stress–strain seven independent method [61], elastic we successfully constants calculated in total: theirC11, values,C12, C13 as, summarizedC16, C33, C44 and C66. Byin applying Table 1. With the stress–strain these values, methodwe can immediately [61], we successfully examine the calculated mechanical their stability values, of asCu summarized2FeSnS4 in Tablecompound1. With theseby employing values, the we Born–Huang can immediately generalized examine elastic thecriteria mechanical [62], as follow: stability of Cu 2FeSnS4 compound by employing the Born–Huang generalized elastic criteria [62], as follow:

Crystals 2020, 10, 509 5 of 12 Crystals 2020, x, x FOR PEER REVIEW 5 of 12

CC11 > |C12| 11 | 12|  2 2 ( + 2C213C13> > 0 2 44  2C16 2> C66(C 11 C12 2C16 < C66(C11-−C12) ItIt is is found found that that all all conditions conditions ofof EquationEquation (1)(1) areare satisfied,satisfied, verifyingverifying the mechanical stability stability of of CuCu2FeSnS2FeSnS44compound. compound. InIn addition,addition, several other mechanical parameters can can also also be be derived derived with with the the Voigt–Reuss–HillVoigt–Reuss–Hill approximation approximation [63 [63,64],,64], such such as bulk as bulk modulus modulusB, Young’s B, Young’s modulus modulusE, shear E modulus, shear Gmodulus, compressibility G, compressibilityβ, Poisson’s β, ratioPoisson’sυ and ratio Pugh’s υ and ratio Pugh’s (B/G), ratio and their(B/G), values and their are alsovalues included are also in Tableincluded1. In general,in Table the 1. In ductile general,/brittle the behaviorductile/brittle of material behavior can of also material be referred can also to its be Poisson’s referred ratioto itsυ . ForPoisson’s ductile ratio materials, υ. For υductileshould materials, be larger υ than should 1/3, be otherwise larger than the 1/3, material otherwise behaves the material in a brittle behaves manner. in Here,a brittle the calculatedmanner. Here,υ of the Cu2 calculatedFeSnS4 is 0.314,υ of Cu which2FeSnS is4 very is 0.314, close which to the thresholdis very close value. to the Furthermore, threshold Pugh’svalue. ratioFurthermore, (B/G) can Pugh’s also be ratio used (B/G) to describe can also the be brittleused to or describe ductile propertythe brittle of or materials. ductile property When it of is largermaterials. than When 1.75, the it is material larger than is ductile 1.75, the and material vice versa. is ductile For Cuand2FeSnS vice versa.4 compound, For Cu2FeSnS its derived4 compound, Pugh’s ratioits derived is 2.36, whichPugh’s is ratio much is larger 2.36, thanwhich the is threshold much larger value. than In combination,the threshold Cu value.2FeSnS In4 iscombination, categorized asCu a2 ductileFeSnS4 material.is categorized as a ductile material.

TableTable 1. 1.The The calculated calculated variousvarious elasticelastic constantsconstants ((CCijij), bulk modulus BB,, Young’s Young’s modulus modulus EE, ,shear shear

modulusmodulusG G, compressibility, compressibilityβ β,, Poisson’s Poisson’s ratio ratioυ υand and Pugh’sPugh’s ratioratio ((BB//GG)) forfor CuCu22FeSnS44..

C 11 C 12 C 13 C 16 C 33 C 44 C 66 B E G β 11 C 12 C 13C 16C 33C 44C 66C B E G β υ B/G (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (TPa 1) υ B/G (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (GPa) (TPa− -1) 69.1269.12 38.84 38.84 33.55 33.55 0.16 0.16 72.69 72.69 21.9521.95 21.8621.86 46.9846.98 52.4152.41 19.9419.94 7.02 7.02 0.314 0.314 2.36 2.36

Furthermore, the mechanical anisotropy is another very important property for the practical Furthermore, the mechanical anisotropy is another very important property for the practical applications of materials, such as the epitaxial thin film growth and the heterojunction combination, applications of materials, such as the epitaxial thin film growth and the heterojunction combination, and it is also accessed for Cu2FeSnS4 compound by calculating the spatial elastic moduli [65,66]. The and it is also accessed for Cu FeSnS compound by calculating the spatial elastic moduli [65,66]. directional dependent Young’s2 modulus4 and shear modulus are shown in Figures 4 and 5, The directional dependent Young’s modulus and shear modulus are shown in Figures4 and5, respectively. The 3D contour plot is accompanied with the 2D projections in different planes. It is respectively.found that the The Young’s 3D contour modulus plot is exhibits accompanied relatively with small the 2D anisotropy projections because in different its spatial planes. shapesIt is found are thatquite the closed Young’s to modulusthe ideal isotropic exhibits relatively sphere, yet small the anisotropy shear modulus because shows its spatial strong shapes directional are quite anisotropy. closed to theThis ideal mechanical isotropic anisotropy sphere, yet is the related shear with modulus the crystal shows symmetry. strong directional For example, anisotropy. the Young’sThis mechanical modulus anisotropyexhibits fourfold is related central with therotoinversion crystal symmetry. symmetry For simply example, because the Young’s the crystal modulus system exhibits has the fourfold same central rotoinversion symmetry simply because the crystal system has the same symmetry operation. symmetry operation. It should be noted here that although Cu2FeSnS4 compound has been Itsynthesized should be noted in the here experiment, that although there are Cu 2stillFeSnS no 4direcompoundct mechanical has been parameters synthesized reported in the yet. experiment, Thus, the theremechanical are still properties no direct mechanical obtained in parameters current work reported can provide yet. Thus, a very the helpful mechanical reference properties for its obtained further instudy. current work can provide a very helpful reference for its further study.

Figure 4. The calculated spatial Young’s modulus for Cu2FeSnS4 compound at the equilibrium lattice Figure 4. The calculated spatial Young’s modulus for Cu2FeSnS4 compound at the equilibrium latticeconstants. constants.

Crystals 2020, x, x FOR PEER REVIEW 6 of 12 Crystals 2020, 10, 509 6 of 12 Crystals 2020, x, x FOR PEER REVIEW 6 of 12

Figure 5. The calculated spatial shear modulus for Cu2FeSnS4 compound at the equilibrium lattice constants. FigureFigure 5.5. TheThe calculated calculated spatial spatial shear shear modulus modulus for Cu for2FeSnS Cu2FeSnS4 compound4 compound at the atequilibrium the equilibrium lattice latticeconstants. constants. As is commonly known, a crystal structure cannot maintain an ideal equilibrium lattice in practicalAsAs is isapplication commonly commonly known, and known, strain a crystal a conditions crystal structure structure often cannot happen, cannot maintain especiallymaintain an ideal anduring equilibriumideal material equilibrium lattice processing inlattice practical andin applicationdifferentpractical applicationstructure and strain connections and conditions strain conditionsin often device happen, construction often especiallyhappen,. For especially during some cases, material during lattice material processing strain processing is and also di usedff erentand on structurepurposedifferent connectionstostructure engineer connections materials’ in device construction.physicalin device properties construction For some [67–70].. cases,For some Consequently, lattice cases, strain lattice is we also strainhave used alsois onalso investigated purpose used on to engineerthepurpose strain materials’to effects engineer of physicalCu materials’2FeSnS properties4. physicalAs this material [properties67–70]. has Consequently, [67–70]. a tetragonal Consequently, westructure, have also we there have investigated are also two investigated independent the strain elatticetheffects strain ofconstants, Cu effects2FeSnS of namely 4Cu. As2FeSnS this a and4 material. As cthis, and material has they a tetragonal arehas aoriented tetragonal structure, perpendicular structure, there are there to two areeach independent two other. independent Thus, lattice we constants,introducedlattice constants, namely lattice namelyastrainand cin, a and a and tetragonal they c, and are orientedmanner,they are that perpendicularoriented is to say,perpendicular one to eachlattice other. isto varied each Thus, other.when we introducedtheThus, other we is latticekeptintroduced constant. strain lattice in The a tetragonal strainconsidered in manner,a tetragonal lattice that variation manner, is to say, range that one is lattice to± say,5% is from variedone latticethe when equilibrium is varied the other when value is kept the for constant.otherboth twois Thelatticekept considered constant. constants. The lattice considered variation lattice range variation is 5% from range the is equilibrium ±5% from the value equilibrium for both two value lattice for both constants. two ± latticeFirstly,Firstly, constants. the the total total energy energy diff differenceerence under under different differe tetragonalnt tetragonal strains strains was calculated was calculated and is plotted and is inplotted FigureFirstly, in6. Figure Three the green total6. Three energy curves green fordifference curves the strain for under athecondition strain differe a andntcondition tetragonal three blue and strains curvesthree blue forwas the curvescalculated strain forc conditionthe and strain is arecplotted condition shown, in Figure and are shown, they 6. Three are and perpendicularly green they curvesare perpendicularly for orientated. the strain Inorientated.a condition addition, Inand two addition, three biaxial blue straintwo curves biaxial curves for strain withthe strain curves black colorwithc condition areblack considered colorare shown, are and considered and they they correspond areand perpendicularly they to correspond the positive orientated. to biaxial the positive strain In addition, condition biaxial two strain with biaxial thecondition same strain variation withcurves the betweensamewith blackvariationa and colorc betweenand are theconsiderednegative a and c andand biaxial theythe strainnegative correspond condition biaxial to the withstrain positive the condition opposite biaxial with variationstrain the condition opposite between with variationa and thec . Itbetweensame can bevariation clearlya and between observedc. It can abe fromand clearly c the and figure observedthe negative that allfrom strainsbiaxial the leadfigurestrain to conditionthat a total all energy strains with increase,the lead opposite to implyinga total variation energy that strainincrease,between distortion implyinga and c makes. It that can strainCu be2 FeSnSclearly distortion4 observedenergetically makes from Cu unstable.2 FeSnSthe figure4 energetically Di thatfferent all strains variations unstable. lead areDifferent to found:a total variations energy a much smallerareincrease, found: total implying a energymuch thatsmaller change strain withtotal distortion strainenergyc makeschangeat zero Cu strainwith2FeSnS straina compared4 energetically c at zero with strain thatunstable. a with compared strainDifferenta withat variations zero that strain with cstrain,are see found: the a at central azero much strain two smaller perpendicular c, see total the energy central curves change two in perpendicular with Figure strain6. This c at curves canzero be strain in simply Figure a compared understood 6. This with can as thatbe follows: simply with strain a at zero strain c, see the central two perpendicular curves in Figure 6. This can be simply atunderstood the same strain as follows: ratio, at the the volume same strain of the ratio, Cu2FeSnS the volume4 crystal of structurethe Cu2FeSnS is varied4 crystal linearly structure with is strain variedc butlinearlyunderstood quadratically with as strain follows: with c but strainat thequadratically samea. It wasstrain also with ratio, found strain the volume that a. It the was of reverse thealso Cu found biaxial2FeSnS that4 straincrystal the reverse gave structure rise biaxial to is avaried much strain linearly with strain c but quadratically with strain a. It was also found that the reverse biaxial strain smallergave rise energy to a much increase smaller than theenergy same increase biaxial strain,than the meaning same biaxial the total strain, energy meaning variation the of total Cu2FeSnS energy4 compoundvariationgave rise ofto was Cua much2 stronglyFeSnS smaller4 compound correlated energy was with increase strongly its volume than correlated the change. same withbiaxial its strain,volume meaning change. the total energy variation of Cu2FeSnS4 compound was strongly correlated with its volume change.

FigureFigure 6. 6.The The calculated calculated total total energy energy difference difference with respect with torespect the equilibrium to the equilibrium condition forcondition Cu FeSnS for Figure 6. The calculated total energy difference with respect to the equilibrium condition2 for 4 compoundCu2FeSnS4 undercompound different under tetragonal different strains. tetragonal strains. Cu2FeSnS4 compound under different tetragonal strains.

Crystals 2020, 10, 509 7 of 12

Crystals 2020, x, x FOR PEER REVIEW 7 of 12

Then, we moved on to examine the strain effects on the electronic properties for Cu2FeSnS4 Then, we moved on to examine the strain effects on the electronic properties for Cu2FeSnS4 compound. Since it shows half-metallic band structure with a small gap in the spin-up direction, compound. Since it shows half-metallic band structure with a small gap in the spin-up direction, we we focused on the variations of CBM and VBM in the same spin under different tetragonal strains. focused on the variations of CBM and VBM in the same spin under different tetragonal strains. The The obtained results are reported in Figure7, where the half-transparent gray area indicated the obtained results are reported in Figure 7, where the half-transparent gray area indicated the presence presence of band gap. We can see that CBM and VBM always descended with the strain increase from of band gap. We can see that CBM and VBM always descended with the strain increase from –5% to –5%+5% to for+5% both for a both and ac.and The cband. The gap band can gap be canpreserved be preserved under a under much a larger much variation larger variation of strain of c at strain zeroc at zerostrain strain a thana than strain strain a at azeroat zero strain strain c, seec ,the see two the plots two plotsin the in central the central column. column. At 105% At strain 105% a strain, the a, theband band gap gap completely completely disappeared disappeared throughout throughout the theentire entire variation variation of strain of strain c, whichc, which led ledto a to total a total transformationtransformation from from half-metal half-metal toto metal.metal.

FigureFigure 7. 7.The The calculated calculated conductionconduction band minimum minimum (CBM) (CBM) and and valence valence band band maximum maximum (VBM) (VBM) in the in the

spin-upspin-up direction direction for for Cu Cu22FeSnSFeSnS44 compound under under different different uniform uniform strains. strains. Grey Grey shade shade indicates indicates the the areaarea with with band band gap. gap.

Finally,Finally, the the total total and and atom-resolved atom-resolved partial magn magneticetic moments moments were were investigated investigated under under different different tetragonaltetragonal strains strains and and theirtheir variationsvariations are shown in in Figure Figure 8.8. Because Because the the total total moment moment is ismainly mainly contributedcontributed to to by by thethe FeFe atoms,atoms, didifferentfferent scales we werere applied applied in in the the vertical vertical axis axis to tohave have a better a better visualization.visualization. It It is is shownshown thatthat thethe totaltotal magnetic moment moment was was almost almost maintained maintained constant constant over over all all strainstrain variations, variations, except except for for aa veryvery smallsmall drop at at the the compressive compressive sides sides for for both both a anda and c ofc 95%of 95% strains, strains, seesee the the two two plots plots in in the the left left column.column. This unchange unchangedd total total moment moment was was due due to tothe the reverse reverse variations variations betweenbetween the the Fe Fe atoms atoms and and thethe otherother three atoms. Th Thee partial partial magnetic magnetic moment moment of ofthe the Fe Featoms atoms shows shows a relatively larger variation with strain a than strain c, and its variation with strain a (or strain c) only a relatively larger variation with strain a than strain c, and its variation with strain a (or strain c) manifests a small offset between different strain c (or strain a) ratios. Overall, we found that the only manifests a small offset between different strain c (or strain a) ratios. Overall, we found that the tetragonal strain has very small impact on the magnetic properties of Cu2FeSnS4. tetragonal strain has very small impact on the magnetic properties of Cu2FeSnS4.

Crystals 2020, 10, 509 8 of 12 Crystals 2020, x, x FOR PEER REVIEW 8 of 12

Figure 8. The calculated total and partial magnetic moments for Cu2FeSnS4 compound under different Figure 8. The calculated total and partial magnetic moments for Cu2FeSnS4 compound under different tetragonaltetragonal strains. strains. Note Note thethe didifferentfferent scalesscales are us useded on on the the two two sides sides of of the the break break position position in inthe the vertical axis. vertical axis.

4. Conclusions4. Conclusions InIn conclusion, conclusion, we we have have performed performed aa systematicsystematic investigation of of the the electronic, electronic, magnetic magnetic and and mechanical properties of the quaternary compound Cu2FeSnS4 using first principles calculations. The mechanical properties of the quaternary compound Cu2FeSnS4 using first principles calculations. Thefully fully optimized optimized lattice lattice structure structure shows shows very very good good consistency consistency with with the the corresponding corresponding experimental experimental one.one. The The spin spin polarized polarized band band structuresstructures showshow half-metallic behavior behavior with with a asmall small band band gap gap in inthe the spin-up direction and multiple bands crossing the Fermi level in the spin-down direction. The spin-up direction and multiple bands crossing the Fermi level in the spin-down direction. The calculated calculated total magnetic moment is 4 μB and it is mostly contributed to by the 3d transition metal total magnetic moment is 4 µ and it is mostly contributed to by the 3d transition metal element Fe, element Fe, which originatedB from its strong spin splitting densities of states and further confirmed which originated from its strong spin splitting densities of states and further confirmed by the large by the large distribution of the magnetic spin density around the Fe atoms. Moreover, the mechanical distribution of the magnetic spin density around the Fe atoms. Moreover, the mechanical properties properties were further studied and several mechanical parameters were obtained, with which the weremechanical further studied stability andwas severalverified. mechanical The spatial dist parametersribution of were its Young’s obtained, modulus with which and shear the mechanicalmodulus stabilitywere computed was verified. and it The showed spatial a very distribution minor mechanical of its Young’s anisotropy. modulus Finally, and the shear tetragonal modulus strain were computedconditions and have it showed been further a very accessed minor mechanical and its effect anisotropy.s on the electronic Finally, the and tetragonal magnetic strain properties conditions for haveCu2 beenFeSnS further4 compound accessed have and been its examined. effects on Results the electronic show the andtotal magnetic magnetic propertiesmoment can for be Cuproperly2FeSnS 4 compoundmaintained have throughout been examined. the whole Results variation show range the totalfor both magnetic strain momenta and c, but can the be half-metallicity properly maintained can throughoutonly be preserved the whole under variation relatively range small for variations. both strain Thea andcurrentc, but work the can half-metallicity provide comprehensive can only be preservedinformation under for relatively quaternary small compound variations. Cu2FeSnS The current4 and a workvaluable can reference provide comprehensivefor its future applications. information forAuthor quaternary Contributions: compound Formal Cu 2analysis,FeSnS4 Jingand Bai, a valuable Jiaying Ji referenceand Liyu Hao; for itsinvestigation, future applications. Jing Bai, Jiaying Ji and Liyu Hao; methodology, Xingwen Tan; project administration, Liyu Hao; supervision, Xingwen Tan; validation, Author Contributions: Formal analysis, J.B., J.J. and L.H.; investigation, J.B., J.J. and L.H.; methodology, X.T.; projectTie Yang; administration, visualization, L.H.; Tie supervision, Yang; writing—original X.T.; validation, draft, T.Y.; Tie visualization, Yang and Xingwen T.Y.; writing—original Tan; writing—review draft, T.Y. and and X.T.;editing, writing—review Xingwen Tan. and editing, X.T. All authors have read and agreed to the published version of the manuscript. Funding:Funding:This This research research received received no no external external funding.funding. Conflicts of Interest: The authors declare no conflicts of interest.

Crystals 2020, 10, 509 9 of 12

References

1. De Groot, R.A.; Mueller, F.M.; Engen, P.G.; Buschow, K.H.J. New Class of Materials: Half-Metallic Ferromagnets. Phys. Rev. Lett. 1983, 50, 2024–2027. [CrossRef] 2. Zutic, I.; Fabian, J.; Das Sarma, S. Spintronics: Fundamentals and applications. Rev. Mod. Phys. 2004, 76, 323. [CrossRef] 3. Gao, G.Y.; Yao, K.L.; Sasioglu, E.; Sandratskii, L.M.; Liu, Z.L.; Jiang, J.L. Half-metallic in zinc-blende CaC, SrC, and BaC from first principles. Phys. Rev. B 2007, 75, 174442. [CrossRef] 4. Ishizuka, H.; Motome, Y. Dirac Half-Metal in a Triangular Ferrimagnet. Phys. Rev. Lett. 2012, 109, 237207. [CrossRef][PubMed] 5. Li, X.; Yang, J. First-principles design of spintronics materials. Natl. Sci. Rev. 2016, 3, 365–381. [CrossRef] 6. Chadov, S.; Wu, S.-C.; Felser, C.; Galanakis, I. Stability of Weyl points in magnetic half-metallic Heusler compounds. Phys. Rev. B 2017, 96, 024435. [CrossRef] 7. Chen, X.; Huang, Y.; Yuan, H.; Liu, J.; Chen, H. Theoretical investigation on thermoelectric properties of spin gapless semiconductor Cr2ZnSi. Appl. Phys. A 2018, 124, 841. [CrossRef] 8. Patel, P.D.; Shinde, S.M.; Gupta, S.D.; Jha, P.K. A promising thermoelectric response of fully compensated ferrimagnetic spin gapless semiconducting Heusler alloy Zr2MnAl at high temperature: DFT study. Mater. Res. Express 2019, 6, 076307. [CrossRef] 9. Hasan, R.; Ur, S. Thermoelectric and Transport Properties of FeV1 xTixSb Half-Heusler System Synthesized − by Controlled Mechanical Alloying Process. Electron. Mater. Lett. 2018, 14, 725–732. [CrossRef] 10. Kara, H.; Kahaly, M.U.; Ozdogan, K. Thermoelectric response of quaternary Heusler compound CrVNbZn. J. Alloys. Compd. 2018, 735, 950–958. [CrossRef] 11. Shigeta, I.; Kubota, T.; Sakuraba, Y.; Kimura, S.; Awaji, S.; Takanashi, K.; Hiroi, M. Transport properties of epitaxial films for superconductor NbN and half-metallic Heusler alloy Co2MnSi under high magnetic fields. Phys. B. Condens. Matter 2018, 536, 310–313. [CrossRef] 12. Nakajima, Y.; Hu, R.; Kirshenbaum, K.; Hughes, A.; Syers, P.; Wang, X.; Wang, K.; Wang, R.; Saha, S.R.; Pratt, D.; et al. Topological RPdBi half-Heusler semimetals: A new family of noncentrosymmetric magnetic superconductors. Sci. Adv. 2015, 1, e1500242. [CrossRef][PubMed] 13. Sprungmann, D.; Westerholt, K.; Zabel, H.; Weides, M.; Kohlstedt, H. Evidence for triplet superconductivity in Josephson junctions with barriers of the ferromagnetic Heusler alloy Cu2MnAl. Phys. Rev. B 2010, 82, 060505. [CrossRef] 14. He, M.; Sun, H.; He, Q.L. : Spintronics and quantum computations. Front. Phys. 2019, 14, 43401. [CrossRef] 15. Chadov, S.; Qi, X.; Kubler, J.; Fecher, G.H.; Felser, C.; Zhang, S.C. Tunable multifunctional topological insulators in ternary Heusler compounds. Nat. Mater. 2010, 9, 541–545. [CrossRef][PubMed] 16. Kainuma, R.; Imano, Y.; Ito, W.; Sutou, Y.; Morito, H.; Okamoto, S.; Kitakami, O.; Oikawa, K.; Fujita, A.; Kanomata, T.; et al. Magnetic-field-induced shape recovery by reverse phase transformation. Nature 2006, 439, 957–960. [CrossRef] 17. Zhu, W.; Liu, E.K.; Feng, L.; Tang, X.D.; Chen, J.L.; Wu, G.H.; Liu, H.Y.; Meng, F.B.; Luo, H.Z. Magnetic-field- induced transformation in FeMnGa alloys. Appl. Phys. Lett. 2009, 95, 222512. [CrossRef] 18. Planes, A.; Mañosa, L.; Acet, M. Magnetocaloric effect and its relation to shape-memory properties in ferromagnetic Heusler alloys. J. Phys. Condens. Matter 2009, 21, 233201. [CrossRef] 19. Wang, J.; Sehitoglu, H. Twinning stress in shape memory alloys: Theory and experiments. Acta Mater. 2013, 61, 6790–6801. [CrossRef] 20. Aznar, A.; Gràcia-Condal, A.; Planes, A.; Lloveras, P.; Barrio, M.; Tamarit, J.-L.; Xiong, W.; Cong, D.; Popescu, C.; Mañosa, L. Giant barocaloric effect in all-d-metal Heusler shape memory alloys. Phys. Rev. Mater. 2019, 3, 044406. [CrossRef] 21. Zeng, Q.; Shen, J.; Zhang, H.; Chen, J.; Ding, B.; Xi, X.; Liu, E.; Wang, W.; Wu, G. Electronic behaviors during martensitic transformations in all-d-metal Heusler alloys. J. Phys. Condens. Matter 2019, 31, 425401. [CrossRef][PubMed] 22. Wang, X.; Khachai, H.; Khenata, R.; Yuan, H.; Wang, L.; Wang, W.; Bouhemadou, A.; Hao, L.; Dai, X.; Guo, R.; et al. Structural, electronic, magnetic, half-metallic, mechanical, and thermodynamic properties of the quaternary Heusler compound FeCrRuSi: A first-principles study. Sci. Rep. 2017, 7, 16183. [CrossRef] [PubMed] Crystals 2020, 10, 509 10 of 12

23. Cherid, S.; Benstaali, W.; Abbad, A.; Bentata, S.; Lantri, T.; Abbar, B. Theoretical prediction of half metallic ferromagnetic full-Heusler alloys Cs2CrGe. Solid State Commun. 2017, 260, 14–18. [CrossRef] 24. Bahramian, S.; Ahmadian, F. Half-metallicity and magnetism of quaternary Heusler compounds CoRuTiZ (Z=Si, Ge, and Sn). J. Magn. Magn. Mater. 2017, 424, 122–129. [CrossRef] 25. Han, Y.; Wu, Y.; Li, T.; Khenata, R.; Yang, T.; Wang, X. Electronic, Magnetic, Half-Metallic, and Mechanical Properties of a New Equiatomic Quaternary Heusler Compound YRhTiGe: A First-Principles Study. Materials (Basel) 2018, 11, 797. [CrossRef] 26. Zheng, A.M.; Huang, H.; Gao, G.Y.; Yao, K.L. Surface half-metallicity in the Heusler alloy Cr2CoGa with low magnetic moment. J. Mater. Sci. 2018, 53, 8364–8371. [CrossRef] 27. Hyun, J.M.; Kim, M. The Half-Metallicity of Co2FeGe Full Heusler Alloy in (001) Thin Film: First Principles Study. J. Korean Phys. Soc. 2018, 72, 276–282. [CrossRef] 28. Wang, J.X.; Chen, Z.B.; Gao, Y.C. Phase stability, magnetic, electronic, half-metallic and mechanical properties of a new equiatomic quaternary Heusler compound ZrRhTiIn: A first-principles investigation. J. Phys. Chem. Solids 2018, 116, 72–78. [CrossRef] 29. Wang, X.; Zhao, W.; Cheng, Z.; Dai, X.; Khenata, R. Electronic, magnetic, half-metallic and mechanical properties of a new quaternary Heusler compound ZrRhTiTl: Insights from first-principles studies. Solid State Commun. 2018, 269, 125–130. [CrossRef] 30. Liu, W.; Zhang, X.; Jia, H.; Khenata, R.; Dai, X.; Liu, G. Theoretical Investigations on the Mechanical, Magneto-Electronic Properties and Half-Metallic Characteristics of ZrRhTiZ (Z = Al, Ga) Quaternary Heusler Compounds. Appl. Sci. 2019, 9, 883. [CrossRef] 31. Chen, Y.; Chen, S.; Wang, B.; Wu, B.; Huang, H.; Qin, X.; Li, D.; Yan, W. Half-Metallicity and Magnetism of the Quaternary Heusler Compound TiZrCoIn1 xGex from the First-Principles Calculations. Appl. Sci. 2019, − 9, 620. [CrossRef] 32. Ganai, Z.S.; Yousuf, S.; Batoo, K.M.; Khan, M.; Gupta, D.C. Half-metallicity and onsite Hubbard interaction on d-electronic states: A case study of Fe2NiZ (Z = Al, Ga, Si, Ge) Heusler systems. Philos. Mag. 2019, 99, 1551–1562. [CrossRef] 33. Rani, D.; Bainsla, L.; Suresh, K.G.; Alam, A. Experimental and theoretical investigation on the possible half-metallic behaviour of equiatomic quaternary Heusler alloys: CoRuMnGe and CoRuVZ (Z = Al, Ga). J. Magn. Magn. Mater. 2019, 492, 165662. [CrossRef] 34. Galanakis, I.; Mavropoulos, P. Spin-polarization and electronic properties of half-metallic Heusler alloys calculated from first principles. J. Phys. Condens. Matter 2007, 19, 315213. [CrossRef] 35. Kudrnovsky, J.; Drchal, V.; Turek, I.; Weinberger, P. Electronic, magnetic, and transport properties and magnetic phase transition in quaternary (Cu,Ni)MnSb Heusler alloys. Phys. Rev. B 2008, 78, 054441. [CrossRef] 36. Dai, X.F.; Liu, G.D.; Fecher, G.H.; Felser, C.; Li, Y.X.; Liu, H.Y. New quarternary half metallic material CoFeMnSi. J. Appl. Phys. 2009, 105, 07E901. [CrossRef] 37. Dannenberg, A.; Siewert, M.; Gruner, M.E.; Wuttig, M.; Entel, P. Competing structural ordering tendencies in Heusler-type alloys with high Curie temperatures: Fe2CoGa1-xZnx studied by first-principles calculations. Phys. Rev. B 2010, 82, 214421. [CrossRef] 38. Gao, G.Y.; Hu, L.; Yao, K.L.; Luo, B.; Liu, N. Large half-metallic gaps in the quaternary Heusler alloys CoFeCrZ (Z=Al, Si, Ga, Ge): A first-principles study. J. Alloys. Compd. 2013, 551, 539–543. [CrossRef] 39. Özdo˘gan,K.; ¸Sa¸sıo˘glu,E.; Galanakis, I. Slater-Pauling behavior in LiMgPdSn-type multifunctional quaternary Heusler materials: Half-metallicity, spin-gapless and magnetic semiconductors. J. Appl. Phys. 2013, 113, 193903. [CrossRef] 40. Delgado, G.E.; Quintero, E.; Tovar, R.; Grima-Gallardo, P.; Quintero, M. Synthesis and crystal structure of the quaternary compound AgFe2GaTe4. J. Alloys. Compd. 2014, 613, 143–145. [CrossRef] 41. Berri, S.; Ibrir, M.; Maouche, D.; Attallah, M. First principles study of structural, electronic and magnetic properties of ZrFeTiAl, ZrFeTiSi, ZrFeTiGe and ZrNiTiAl. J. Magn. Magn. Mater. 2014, 371, 106–111. [CrossRef] 42. Berri, S.; Ibrir, M.; Maouche, D.; Attallah, M. Robust half-metallic ferromagnet of quaternary Heusler compounds ZrCoTiZ(Z = Si, Ge, Ga and Al). Comput. Condens. Matter 2014, 1, 26–31. [CrossRef] 43. Delgado, G.; Delgado, N.; Gallardo, P. Crystal structure of the quaternary semiconductor compound CuFeCrSe3. Period. Tche Quim. 2019, 16, 1. Crystals 2020, 10, 509 11 of 12

44. Hao, L.; Tan, M.; Khenata, R.; Wang, X.; Yang, T. Computational study of the electronic, magnetic, mechanical and thermodynamic properties of new equiatomic quaternary heusler compounds TiZrRuZ (Z = Al, Ga, In). Chin. J. Phys. 2019, 62, 54–64. [CrossRef]

45. Schäfer, W.; Nitsche, R. Tetrahedral quaternary chalcogenides of the type Cu2-II-IV-S4(Se4). Mater. Res. Bull. 1974, 9, 645–654. [CrossRef] 46. Wintenberger, M. Etude de la structure cristallographique et magnetique de Cu2FeGeS4 et remarque sur la structure magnetique de Cu2MnSnS4. Mater. Res. Bull. 1979, 14, 1195–1202. [CrossRef] 47. Delgado, G.E.; Sierralta, N.; Quintero, M.; Quintero, E.; Moreno, E.; Flores-Cruz, J.A.; Rincon, C. Synthesis,

structural characterization and differential thermal analysis of the quaternary compound Ag2MnSnS4. Rev. Mex. Fis. 2018, 64, 216–221. [CrossRef] 48. Delgado, G.; Sagredo, A. Synthesis and crystal structure of the quaternary semiconductor Cu2NiGeS4, a new stannite-type compound. Rev. Mex. Física 2019, 65, 355. [CrossRef] 49. Guen, L.; Glaunsinger, W.S. Electrical, magnetic, and EPR studies of the quaternary chalcogenides Cu2AIIBIVX4 prepared by iodine transport. J. Solid State Chem. 1980, 35, 10–21. [CrossRef] 50. Giriat, W.; Furdyna, J. Chapter 1 Crystal Structure, Composition, and Materials Preparation of Diluted Magnetic Semiconductors; Elsevier BV: Amsterdam, The Netherlands, 1988; Volume 25, pp. 1–34. 51. Grima-Gallardo, P.; Cárdenas, K.; Molina, L.; Quintero, M.; Ruiz, J.; Delgado, G.; Briceño, J.M. A Comparative Study of (Cu–III–Se2)x–(FeSe)1—x Alloys (III: Al, Ga, In) (0 x 1) by X-Ray Diffraction, Differential ≤ ≤ Thermal Analysis and Scanning Electron Microscopy. Phys. Status Solidi (a) 2001, 187, 395–406. [CrossRef] 52. Chambers, S.A.; Yoo, Y.K. New Materials for Spintronics. MRS Bull. 2003, 28, 706–710. [CrossRef] 53. Catella, G.C.; Burlage, D. Crystal Growth and Optical Properties of AgGaS2 and AgGaSe2. MRS Bull. 1998, 23, 28–36. [CrossRef] 54. Tsuji, I.; Shimodaira, Y.; Kato, H.; Kobayashi, H.; Kudo, A. Novel Stannite-type Complex Sulfide Photocatalysts AI2-Zn-AIV-S4 (AI = Cu and Ag; AIV = Sn and Ge) for Hydrogen Evolution under Visible-Light Irradiation. Chem. Mater. 2010, 22, 1402–1409. [CrossRef] 55. Rincón, C.; Quintero, E.; Quintero, M.; Moreno, E.; Power, C.; Morocoima, M.; Delgado, G.E. Temperature Dependence of Raman Spectra in Cu2FeSnS4 Magnetic Semiconductor Compound. Phys. Status Solidi (b) 2019, 10, 1900076. 56. Rincón, C.; Quintero, M.; Moreno, E.; Power, C.; Quintero, E.; Henao, J.A.; Macías, M.A.; Delgado, G.E.; Tovar, R.; Morocoima, M. X-ray diffraction, Raman spectrum and magnetic susceptibility of the magnetic semiconductor Cu2FeSnS4. Solid State Commun. 2011, 151, 947–951. [CrossRef] 57. Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B. Condens. Matter. 1990, 41, 7892–7895. [CrossRef] 58. Payne, M.C.; Teter, M.P.; Allan, D.C.; Arias, T.A.; Joannopoulos, J.D. Iterative minimization techniques forab initiototal-energy calculations: Molecular dynamics and conjugate gradients. Rev. Mod. Phys. 1992, 64, 1045–1097. [CrossRef] 59. Segall, M.D.; Lindan, P.J.D.; Probert, M.J.; Pickard, C.J.; Hasnip, P.J.; Clark, S.J.; Payne, M.C. First-principles simulation: Ideas, illustrations and the CASTEP code. J. Phys. Condens. Matter 2002, 14, 2717–2744. [CrossRef] 60. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [CrossRef] 61. Wang, J.; Li, J.; Yip, S.; Phillpot, S.; Wolf, D. Mechanical instabilities of homogeneous crystals. Phys. Rev. B. Condens. Matter 1995, 52, 12627–12635. [CrossRef] 62. Mouhat, F.; Coudert, F.X. Necessary and sufficient elastic stability conditions in various crystal systems. Phys. Rev. B 2014, 90, 224104. [CrossRef] 63. Hill, R. The Elastic Behaviour of a Crystalline Aggregate. Proc. Phys. Soc. Sect. A 1952, 65, 349. [CrossRef] 64. Yip, S.; Li, J.; Tang, M.J.; Wang, J.G. Mechanistic aspects and atomic-level consequences of elastic instabilities in homogeneous crystals. Mater. Sci. Eng. A Struct. Mater. Prop. Microstruct. Process. 2001, 317, 236–240. [CrossRef] 65. Gaillac, R.; Pullumbi, P.;Coudert, F.X. ELATE: An open-source online application for analysis and visualization of elastic tensors. J. Phys. Condens. Matter 2016, 28, 275201. [CrossRef][PubMed] 66. Marmier, A.; Lethbridge, Z.A.D.; Walton, R.I.; Smith, C.W.; Parker, S.C.; Evans, K.E. ElAM: A computer program for the analysis and representation of anisotropic elastic properties. Comput. Phys. Commun. 2010, 181, 2102–2115. [CrossRef] Crystals 2020, 10, 509 12 of 12

67. Yang, T.; Hao, L.Y.; Khenata, R.; Wang, X.T. Strain Conditions for the Inverse Heusler Type Fully Compensated Spin-Gapless Semiconductor Ti2MnAl: A First-Principles Study. Materials 2018, 11, 2091. [CrossRef] 68. Yang, T.; Cao, J.; Khenata, R.; Cheng, Z.; Kuang, M.; Wang, X. Strain effect for the newly discovered spin-gapless diamond-like quaternary-type semiconductor CuMn2InSe4. J. Alloys. Compd. 2019, 793, 302–313. [CrossRef] 69. Yang, T.; Hao, L.; Khenata, R.; Wang, X. Investigation of the structural competing and atomic ordering in Heusler compounds Fe2NiSi and Ni2FeSi under strain condition. R. Soc. Open Sci. 2019, 6, 191007. [CrossRef] 70. You, J.; Cao, J.; Khenata, R.; Wang, X.; Shen, X.; Yang, T. Robust Spin-Gapless Behavior in the Newly Discovered Half Heusler Compound MnPK. Materials (Basel) 2019, 12, 3117. [CrossRef]

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