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

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 shear modulus 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.


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 spintronics and magnetoelectronics [1][2][3][4][5][6]. Half-metallicity is characterized by the spin polarized band structures with a band gap in one spin direction, similar to semiconductors, 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][8][9][10], superconductivity [11][12][13], topological insulativity [14,15], magnetic shape memory [16][17][18][19][20][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][23][24][25][26][27][28][29][30][31][32][33]. Compared with ternary materials, quaternary ones are more promising because a larger compositional flexibility can be realized [34][35][36][37][38][39][40][41][42][43][44]. A large and extensively studied quaternary material group is the A I 2 B II C IV D VI 4 compounds [45][46][47][48][49], in which A I is Cu or Ag; B II is from Mn, Fe, Ni, Zn and Cd; C IV is from IV main group elements Si, Ge and Sn; and D VI is from VI main group elements S, Se and Te. They are originated from tetrahedrally coordinated derivatives of B II -D VI 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, Cu 2 FeSnS 4 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 Cu 2 FeSnS 4 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 transition metal 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 Cu 2 FeSnS 4 and serve as a helpful guide for its future applications.

Computational Methodology
The electronic, magnetic and mechanical properties of the quaternary compound Cu 2 FeSnS 4 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 3d 10 4s 1 , 3d 6 4s 2 , 5s 2 5p 2 and 3s 2 3p4 2 , 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.

Results and Discussions
The bulk quaternary compound Cu 2 FeSnS 4 was experimentally synthesized and crystallized with a tetragonal symmetry; space group P4 (S 1 4 , No. 81). Its crystal structure is schematically illustrated in Figure 1 and the unit cell contained eight atoms in total, including one Sn located at Wyckoff 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 Figure 1. The fully relaxed lattice constants for Cu 2 FeSnS 4 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. Based on the optimized crystal structure for Cu2FeSnS4, we calculated its electronic band structures. Since the 3d transition metal element, Fe, carries a large magnetic moment, the band structures show spin polarized features (see Figure 2). The fermi energy is set to zero and only the energy range near the Fermi level is shown. In the spin-up direction, a small indirect band gap was observed with a conduction band minimum (CBM) located at high symmetry point A and the valence band maximum (VBM) at G. In the spin-down direction, two bands show multiple crossings with the Fermi energy level. This behavior of semiconductor properties in one spin direction and metal properties in the other spin direction clearly indicated that this material belongs to the family of halfmetal compounds. To further examine the correlation states between different elements and the magnetism origination for Cu2FeSnS4 compound, we have also computed the densities of states and the distributions of the magnetic spin density and the charge density difference, as reported in Figure 3. The positive and negative values in the total and partial densities of states in Figure 3a,b correspond Based on the optimized crystal structure for Cu 2 FeSnS 4 , we calculated its electronic band structures. Since the 3d transition metal element, Fe, carries a large magnetic moment, the band structures show spin polarized features (see Figure 2). The fermi energy is set to zero and only the energy range near the Fermi level is shown. In the spin-up direction, a small indirect band gap was observed with a conduction band minimum (CBM) located at high symmetry point A and the valence band maximum (VBM) at G. In the spin-down direction, two bands show multiple crossings with the Fermi energy level. This behavior of semiconductor properties in one spin direction and metal properties in the other spin direction clearly indicated that this material belongs to the family of half-metal compounds. Based on the optimized crystal structure for Cu2FeSnS4, we calculated its electronic band structures. Since the 3d transition metal element, Fe, carries a large magnetic moment, the band structures show spin polarized features (see Figure 2). The fermi energy is set to zero and only the energy range near the Fermi level is shown. In the spin-up direction, a small indirect band gap was observed with a conduction band minimum (CBM) located at high symmetry point A and the valence band maximum (VBM) at G. In the spin-down direction, two bands show multiple crossings with the Fermi energy level. This behavior of semiconductor properties in one spin direction and metal properties in the other spin direction clearly indicated that this material belongs to the family of halfmetal compounds. To further examine the correlation states between different elements and the magnetism origination for Cu2FeSnS4 compound, we have also computed the densities of states and the distributions of the magnetic spin density and the charge density difference, as reported in Figure 3. The positive and negative values in the total and partial densities of states in Figure 3a,b correspond To further examine the correlation states between different elements and the magnetism origination for Cu 2 FeSnS 4 compound, we have also computed the densities of states and the distributions of the magnetic spin density and the charge density difference, as reported in Figure 3. The positive and negative values in the total and partial densities of states in Figure 3a,b correspond to the spin-up and spin-down directions, respectively. A spin polarized behavior is observed: zero density at the Fermi level in spin-up direction and nonzero density in spin-down direction, which is in accordance with the findings from the electronic band structures. The total density of states near the Fermi level is mostly contributed from the transition metal elements Fe and Cu, especially in the energy range from −2 eV to +2 eV, and the main group elements, Sn and S, have negligible contributions. However, the two transition metal elements exhibit different density distributions: an almost symmetric one for Cu element from −2 eV to 0 eV between the two spin directions; a highly asymmetric one for Fe element with two high protrude areas in spin-up direction at −3 to −2 eV; and several high peaks in spin-down direction at 0 to1 eV. Especially for the density at the Fermi level in spin-down direction, it is completely from the Fe element. This strong spin splitting effect in the density of states for the Fe element could result in its strong magnetic moment, which indeed is confirmed by the calculated magnetic spin density as shown in Figure 3c. The large gray shaded areas around the Fe atoms indicate a strong local magnetic moment. Some shaded areas are also found around Cu atoms but they are much smaller compared with that around the Fe atoms. The calculated total magnetic moment per unit formula is 4 µ B and it is mainly contributed by the Fe atom. The charge density difference is also computed to show the atomic bonding between different atoms and it is displayed in Figure 3d with the gray shaded areas representing the electron accumulation. We can clearly see that the S atoms form four bonds with the four neighbor atoms located at the corners of the slightly distorted tetrahedron and they have almost identical electron excess with only one a little bigger between S and Fe atoms. A very strong d z 2 orbital shape is observed at the Fe atoms with the large dumbbell oriented to the center position. to the spin-up and spin-down directions, respectively. A spin polarized behavior is observed: zero density at the Fermi level in spin-up direction and nonzero density in spin-down direction, which is in accordance with the findings from the electronic band structures. The total density of states near the Fermi level is mostly contributed from the transition metal elements Fe and Cu, especially in the energy range from −2 eV to +2 eV, and the main group elements, Sn and S, have negligible contributions. However, the two transition metal elements exhibit different density distributions: an almost symmetric one for Cu element from −2 eV to 0 eV between the two spin directions; a highly asymmetric one for Fe element with two high protrude areas in spin-up direction at −3 to −2 eV; and several high peaks in spin-down direction at 0 to1 eV. Especially for the density at the Fermi level in spin-down direction, it is completely from the Fe element. This strong spin splitting effect in the density of states for the Fe element could result in its strong magnetic moment, which indeed is confirmed by the calculated magnetic spin density as shown in Figure 3c. The large gray shaded areas around the Fe atoms indicate a strong local magnetic moment. Some shaded areas are also found around Cu atoms but they are much smaller compared with that around the Fe atoms. The calculated total magnetic moment per unit formula is 4 μB and it is mainly contributed by the Fe atom. The charge density difference is also computed to show the atomic bonding between different atoms and it is displayed in Figure 3d with the gray shaded areas representing the electron accumulation. We can clearly see that the S atoms form four bonds with the four neighbor atoms located at the corners of the slightly distorted tetrahedron and they have almost identical electron excess with only one a little bigger between S and Fe atoms. A very strong d z 2 orbital shape is observed at the Fe atoms with the large dumbbell oriented to the center position. Next, we turn to the mechanical properties of Cu2FeSnS4 compound. Since it has a P4 tetragonal structure, there should be seven independent elastic constants in total: C11, C12, C13, C16, C33, C44 and C66. By applying the stress-strain method [61], we successfully calculated their values, as summarized in Table 1. With these values, we can immediately examine the mechanical stability of Cu2FeSnS4 compound by employing the Born-Huang generalized elastic criteria [62], as follow: Next, we turn to the mechanical properties of Cu 2 FeSnS 4 compound. Since it has a P4 tetragonal structure, there should be seven independent elastic constants in total: C 11 , C 12 , C 13 , C 16 , C 33 , C 44 and C 66 . By applying the stress-strain method [61], we successfully calculated their values, as summarized in Table 1. With these values, we can immediately examine the mechanical stability of Cu 2 FeSnS 4 compound by employing the Born-Huang generalized elastic criteria [62], as follow: (1) It is found that all conditions of Equation (1) are satisfied, verifying the mechanical stability of Cu 2 FeSnS 4 compound. In addition, several other mechanical parameters can also be derived with the Voigt-Reuss-Hill approximation [63,64], such as bulk modulus B, Young's modulus E, shear modulus G, compressibility β, Poisson's ratio υ and Pugh's ratio (B/G), and their values are also included in Table 1. In general, the ductile/brittle behavior of material can also be referred to its Poisson's ratio υ.
For ductile materials, υ should be larger than 1/3, otherwise the material behaves in a brittle manner. Here, the calculated υ of Cu 2 FeSnS 4 is 0.314, which is very close to the threshold value. Furthermore, Pugh's ratio (B/G) can also be used to describe the brittle or ductile property of materials. When it is larger than 1.75, the material is ductile and vice versa. For Cu 2 FeSnS 4 compound, its derived Pugh's ratio is 2.36, which is much larger than the threshold value. In combination, Cu 2 FeSnS 4 is categorized as a ductile material.  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, and it is also accessed for Cu 2 FeSnS 4 compound by calculating the spatial elastic moduli [65,66]. The directional dependent Young's modulus and shear modulus are shown in Figures 4 and 5, respectively. The 3D contour plot is accompanied with the 2D projections in different planes. It is found that the Young's modulus exhibits relatively small anisotropy because its spatial shapes are quite closed to the ideal isotropic sphere, yet the shear modulus shows strong directional anisotropy. This mechanical anisotropy is related with the crystal symmetry. For example, the Young's modulus exhibits fourfold central rotoinversion symmetry simply because the crystal system has the same symmetry operation. It should be noted here that although Cu 2 FeSnS 4 compound has been synthesized in the experiment, there are still no direct mechanical parameters reported yet. Thus, the mechanical properties obtained in current work can provide a very helpful reference for its further study. (1) It is found that all conditions of Equation (1) are satisfied, verifying the mechanical stability of Cu2FeSnS4 compound. In addition, several other mechanical parameters can also be derived with the Voigt-Reuss-Hill approximation [63,64], such as bulk modulus B, Young's modulus E, shear modulus G, compressibility β, Poisson's ratio υ and Pugh's ratio (B/G), and their values are also included in Table 1. In general, the ductile/brittle behavior of material can also be referred to its Poisson's ratio υ. For ductile materials, υ should be larger than 1/3, otherwise the material behaves in a brittle manner. Here, the calculated υ of Cu2FeSnS4 is 0.314, which is very close to the threshold value. Furthermore, Pugh's ratio (B/G) can also be used to describe the brittle or ductile property of materials. When it is larger than 1.75, the material is ductile and vice versa. For Cu2FeSnS4 compound, its derived Pugh's ratio is 2.36, which is much larger than the threshold value. In combination, Cu2FeSnS4 is categorized as a ductile material.  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, and it is also accessed for Cu2FeSnS4 compound by calculating the spatial elastic moduli [65,66]. The directional dependent Young's modulus and shear modulus are shown in Figures 4 and 5, respectively. The 3D contour plot is accompanied with the 2D projections in different planes. It is found that the Young's modulus exhibits relatively small anisotropy because its spatial shapes are quite closed to the ideal isotropic sphere, yet the shear modulus shows strong directional anisotropy. This mechanical anisotropy is related with the crystal symmetry. For example, the Young's modulus exhibits fourfold central rotoinversion symmetry simply because the crystal system has the same symmetry operation. It should be noted here that although Cu2FeSnS4 compound has been synthesized in the experiment, there are still no direct mechanical parameters reported yet. Thus, the mechanical properties obtained in current work can provide a very helpful reference for its further study.   As is commonly known, a crystal structure cannot maintain an ideal equilibrium lattice in practical application and strain conditions often happen, especially during material processing and different structure connections in device construction. For some cases, lattice strain is also used on purpose to engineer materials' physical properties [67][68][69][70]. Consequently, we have also investigated the strain effects of Cu2FeSnS4. As this material has a tetragonal structure, there are two independent lattice constants, namely a and c, and they are oriented perpendicular to each other. Thus, we introduced lattice strain in a tetragonal manner, that is to say, one lattice is varied when the other is kept constant. The considered lattice variation range is ±5% from the equilibrium value for both two lattice constants.

C12 (GPa)
Firstly, the total energy difference under different tetragonal strains was calculated and is plotted in Figure 6. Three green curves for the strain a condition and three blue curves for the strain c condition are shown, and they are perpendicularly orientated. In addition, two biaxial strain curves with black color are considered and they correspond to the positive biaxial strain condition with the same variation between a and c and the negative biaxial strain condition with the opposite variation between a and c. It can be clearly observed from the figure that all strains lead to a total energy increase, implying that strain distortion makes Cu2FeSnS4 energetically unstable. Different variations are found: a much smaller total energy change with strain c at zero strain a compared with that with strain a at zero strain c, see the central two perpendicular curves in Figure 6. This can be simply understood as follows: at the same strain ratio, the volume of the Cu2FeSnS4 crystal structure is varied linearly with strain c but quadratically with strain a. It was also found that the reverse biaxial strain gave rise to a much smaller energy increase than the same biaxial strain, meaning the total energy variation of Cu2FeSnS4 compound was strongly correlated with its volume change.  As is commonly known, a crystal structure cannot maintain an ideal equilibrium lattice in practical application and strain conditions often happen, especially during material processing and different structure connections in device construction. For some cases, lattice strain is also used on purpose to engineer materials' physical properties [67][68][69][70]. Consequently, we have also investigated the strain effects of Cu 2 FeSnS 4 . As this material has a tetragonal structure, there are two independent lattice constants, namely a and c, and they are oriented perpendicular to each other. Thus, we introduced lattice strain in a tetragonal manner, that is to say, one lattice is varied when the other is kept constant. The considered lattice variation range is ±5% from the equilibrium value for both two lattice constants.
Firstly, the total energy difference under different tetragonal strains was calculated and is plotted in Figure 6. Three green curves for the strain a condition and three blue curves for the strain c condition are shown, and they are perpendicularly orientated. In addition, two biaxial strain curves with black color are considered and they correspond to the positive biaxial strain condition with the same variation between a and c and the negative biaxial strain condition with the opposite variation between a and c. It can be clearly observed from the figure that all strains lead to a total energy increase, implying that strain distortion makes Cu 2 FeSnS 4 energetically unstable. Different variations are found: a much smaller total energy change with strain c at zero strain a compared with that with strain a at zero strain c, see the central two perpendicular curves in Figure 6. This can be simply understood as follows: at the same strain ratio, the volume of the Cu 2 FeSnS 4 crystal structure is varied linearly with strain c but quadratically with strain a. It was also found that the reverse biaxial strain gave rise to a much smaller energy increase than the same biaxial strain, meaning the total energy variation of Cu 2 FeSnS 4 compound was strongly correlated with its volume change.
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.
As is commonly known, a crystal structure cannot maintain an ideal equilibrium lattice in practical application and strain conditions often happen, especially during material processing and different structure connections in device construction. For some cases, lattice strain is also used on purpose to engineer materials' physical properties [67][68][69][70]. Consequently, we have also investigated the strain effects of Cu2FeSnS4. As this material has a tetragonal structure, there are two independent lattice constants, namely a and c, and they are oriented perpendicular to each other. Thus, we introduced lattice strain in a tetragonal manner, that is to say, one lattice is varied when the other is kept constant. The considered lattice variation range is ±5% from the equilibrium value for both two lattice constants.
Firstly, the total energy difference under different tetragonal strains was calculated and is plotted in Figure 6. Three green curves for the strain a condition and three blue curves for the strain c condition are shown, and they are perpendicularly orientated. In addition, two biaxial strain curves with black color are considered and they correspond to the positive biaxial strain condition with the same variation between a and c and the negative biaxial strain condition with the opposite variation between a and c. It can be clearly observed from the figure that all strains lead to a total energy increase, implying that strain distortion makes Cu2FeSnS4 energetically unstable. Different variations are found: a much smaller total energy change with strain c at zero strain a compared with that with strain a at zero strain c, see the central two perpendicular curves in Figure 6. This can be simply understood as follows: at the same strain ratio, the volume of the Cu2FeSnS4 crystal structure is varied linearly with strain c but quadratically with strain a. It was also found that the reverse biaxial strain gave rise to a much smaller energy increase than the same biaxial strain, meaning the total energy variation of Cu2FeSnS4 compound was strongly correlated with its volume change.  Then, we moved on to examine the strain effects on the electronic properties for Cu 2 FeSnS 4 compound. Since it shows half-metallic band structure with a small gap in the spin-up direction, we focused on the variations of CBM and VBM in the same spin under different tetragonal strains. The obtained results are reported in Figure 7, where the half-transparent gray area indicated the presence of band gap. We can see that CBM and VBM always descended with the strain increase from -5% to +5% for both a and c. The band gap can be preserved under a much larger variation of strain c at zero strain a than strain a at zero strain c, see the two plots in the central column. At 105% strain a, the band gap completely disappeared throughout the entire variation of strain c, which led to a total transformation from half-metal to metal.
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 compound. Since it shows half-metallic band structure with a small gap in the spin-up direction, we focused on the variations of CBM and VBM in the same spin under different tetragonal strains. The obtained results are reported in Figure 7, where the half-transparent gray area indicated the presence of band gap. We can see that CBM and VBM always descended with the strain increase from -5% to +5% for both a and c. The band gap can be preserved under a much larger variation of strain c at zero strain a than strain a at zero strain c, see the two plots in the central column. At 105% strain a, the band gap completely disappeared throughout the entire variation of strain c, which led to a total transformation from half-metal to metal. Finally, the total and atom-resolved partial magnetic moments were investigated under different tetragonal strains and their variations are shown in Figure 8. Because the total moment is mainly contributed to by the Fe atoms, different scales were applied in the vertical axis to have a better visualization. It is shown that the total magnetic moment was almost maintained constant over all strain variations, except for a very small drop at the compressive sides for both a and c of 95% strains, see the two plots in the left column. This unchanged total moment was due to the reverse variations between the Fe atoms and the other three atoms. The partial magnetic moment of the Fe atoms shows a relatively larger variation with strain a than strain c, and its variation with strain a (or strain c) 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. Finally, the total and atom-resolved partial magnetic moments were investigated under different tetragonal strains and their variations are shown in Figure 8. Because the total moment is mainly contributed to by the Fe atoms, different scales were applied in the vertical axis to have a better visualization. It is shown that the total magnetic moment was almost maintained constant over all strain variations, except for a very small drop at the compressive sides for both a and c of 95% strains, see the two plots in the left column. This unchanged total moment was due to the reverse variations between the Fe atoms and the other three atoms. The partial magnetic moment of the Fe atoms shows a relatively larger variation with strain a than strain c, and its variation with strain a (or strain c) 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 Cu 2 FeSnS 4 .

Conclusions
In conclusion, we have performed a systematic investigation of the electronic, magnetic and mechanical properties of the quaternary compound Cu2FeSnS4 using first principles calculations. The fully optimized lattice structure shows very good consistency with the corresponding experimental one. The spin polarized band structures show half-metallic behavior with a small band gap in the spin-up direction and multiple bands crossing the Fermi level in the spin-down direction. The calculated total magnetic moment is 4 μB and it is mostly contributed to by the 3d transition metal element Fe, which originated from its strong spin splitting densities of states and further confirmed by the large distribution of the magnetic spin density around the Fe atoms. Moreover, the mechanical properties were further studied and several mechanical parameters were obtained, with which the mechanical stability was verified. The spatial distribution of its Young's modulus and shear modulus were computed and it showed a very minor mechanical anisotropy. Finally, the tetragonal strain conditions have been further accessed and its effects on the electronic and magnetic properties for Cu2FeSnS4 compound have been examined. Results show the total magnetic moment can be properly maintained throughout the whole variation range for both strain a and c, but the half-metallicity can only be preserved under relatively small variations. The current work can provide comprehensive information for quaternary compound Cu2FeSnS4 and a valuable reference for its future applications.

Conclusions
In conclusion, we have performed a systematic investigation of the electronic, magnetic and mechanical properties of the quaternary compound Cu 2 FeSnS 4 using first principles calculations.
The fully optimized lattice structure shows very good consistency with the corresponding experimental one. The spin polarized band structures show half-metallic behavior with a small band gap in the spin-up direction and multiple bands crossing the Fermi level in the spin-down direction. The calculated total magnetic moment is 4 µ B and it is mostly contributed to by the 3d transition metal element Fe, which originated from its strong spin splitting densities of states and further confirmed by the large distribution of the magnetic spin density around the Fe atoms. Moreover, the mechanical properties were further studied and several mechanical parameters were obtained, with which the mechanical stability was verified. The spatial distribution of its Young's modulus and shear modulus were computed and it showed a very minor mechanical anisotropy. Finally, the tetragonal strain conditions have been further accessed and its effects on the electronic and magnetic properties for Cu 2 FeSnS 4 compound have been examined. Results show the total magnetic moment can be properly maintained throughout the whole variation range for both strain a and c, but the half-metallicity can only be preserved under relatively small variations. The current work can provide comprehensive information for quaternary compound Cu 2 FeSnS 4 and a valuable reference for its future applications.

Conflicts of Interest:
The authors declare no conflicts of interest.