Initio Study of Optoelectronic and Magnetic Properties of Ternary Chromium Chalcogenides

Using first-principles calculations, we investigate the magnetic order in the ground state of several ternary chromium chalcogenide compounds. Electronic band structure calculations indicate that these compounds are either metallic or semiconductors with relatively low bandgap energies. &e large optical absorption coefficients, predicted by our calculations, suggest that some of these compounds may be useful as light harvesters in solar cells or as infrared detectors.


Introduction
Ternary compounds of the form ABX 3 , where A and B are metal atoms and X is a halogen or a chalcogen, are exciting candidates as semiconductors for photovoltaics and other applications.ese materials are highly tunable due to the large variety of possible elements that can be included.In addition to tunability, their low cost, facile synthesis, and long carrier lifetimes have gained the attention of the research community [1][2][3].In particular, perovskite solar cells made from lead halide materials (APbX 3 ) have enjoyed rapid development, with measured power conversion efficiencies rivaling state-of-the-art silicon devices [4][5][6][7].As a result, there is interest in commercial application of lead-halide perovskites for photovoltaics; however, concerns regarding the toxicity of lead and high instability in the presence of moisture have hindered these efforts [2].
us, it would be of great interest to find materials with better stability and nontoxic composition, while still keeping the attractive properties of these leadhalide compounds.
Recently, ternary metal chalcogenide materials have been proposed as alternatives to these lead-halide perovskites [8,9].Materials of the form ABO 3 have been extensively studied and have been shown to be resilient in the presence of moisture [8,10].Unfortunately, the energy difference between transition metal d orbitals and the 2p orbitals of oxygen results in a large bandgap; in fact, BiFeO 3 , with a bandgap of 2.7 eV, has the lowest bandgap among this class of materials [10].Solid solutions have been reported with bandgaps between 1.1 and 3.8 eV [8]; however, these materials have high intrinsic disorder.
To lower the bandgap of these ABO 3 materials, substitution of oxygen with chalcogens such as sulfur or selenium has been recently explored.Although reports of synthesis of these ternary transition metal chalcogenides have been around for a long time [11], there are still few studies focusing on their optoelectronic properties.eoretical studies have been reported on substitution of oxygen with sulfur [12,13].Additionally, experimental reactions of perovskite oxides with gaseous H 2 S or CS 2 [14], as well as solid-state reactions have been recently reported as methods of substituting oxygen with sulfur [15,16].It has been shown that transition metal chalcogenides are widely tunable, ranging from metallic to semiconducting with gaps exceeding 2 eV [9,16].Several materials are predicted to have bandgaps corresponding to theoretical photovoltaic power conversion efficiencies above 28 percent [9], and zirconium sulfide perovskites have been experimentally measured to have strong optical absorption and a bandgap suitable for photovoltaic applications [16].
One interesting subset of these ternary transition metal chalcogenides are those of the form ACrX 3 , where X � S, Se, or Te.Because of the presence of chromium, these materials tend to exhibit magnetic properties and have recently seen a resurgence of interest [17,18].Several of these materials, such as SbCrS 3 [19], SbCrSe 3 [18][19][20], and SiCrTe 3 [21,22], have been the subjects of recent experimental studies as magnetic materials with metallic or semiconducting behavior.However, many of these materials, including many of the lanthanide chromium sulfides, have not been thoroughly studied even though their syntheses have been reported a long time ago [11].
In this work, we present first-principles calculations on some chromium chalcogenides where experimentally determined crystal data are available, but electronic or magnetic properties have not been thoroughly investigated.Using total energy calculations, we determine the groundstate magnetic ordering at 0 K and then calculate optoelectronic properties to determine their usefulness in various applications.
With the exception of PCrSe 3 and SiCrTe 3 , all the structures examined have the space group 62 (Pnma) and they adopt the NH 4 CdCl 3 prototype structure.is structure contains dimers of CrX 6 (X � S, Se, or Te) octahedra and ninefold coordination between the A site and the chalcogen atoms.Unlike the perovskite structure, which contains corner-sharing octahedra [23,24], the octahedra in this structure share edges, as shown in Figure 1.
e structure of PCrSe 3 is monoclinic with space group 12 (C2/m).is structure contains layers of a single P atom surrounded by six CrSe 6 octahedra in the ab plane, as pictured in Figure 2.
e structure of SiCrTe 3 is rhombohedral with space group 148 (R3) and is shown in Figure 3.In this structure, edge-sharing CrTe 6 octahedra form layers on the ab plane, while the Si atoms have threefold coordination with the Te atoms.

Methods
To determine the ground-state magnetic ordering of the chromium chalcogenides, we used density functional theory (DFT) as implemented in the Vienna Ab Initio Simulation Package (VASP) [25][26][27][28] code, which uses the projector-augmented wave (PAW) pseudopotentials and a plane-wave basis set to expand the electronic wave function.
e Perdew-Burke-Ernzerhof implementation of the generalized gradient approximation (PBE-GGA) was used to model the exchange-correlation effects [29,30].We performed spin-polarized calculations, where varying spin orderings were imposed on the chromium ions.A Hubbard on-site Coulomb term (U � 4 eV) and Hund's exchange term (J � 0.7 eV) for the chromium ions are included in the calculation.Experimentally determined crystal structures from the Inorganic Crystal Structure Database were used in our calculations, with the exception of SiCrTe 3 , where the X-ray diffraction measurements at 292 K of Casto et al. were used [21].e total energies of different magnetic orders were compared to determine the ground state.e electronic band structures, density of states, and optical properties were calculated using the full-potential, linearized, augmented plane wave (FP-LAPW) method, as implemented in the WIEN2k code [31].Here, space is divided into an interior region of nonoverlapping muffin tin spheres centered on each atom and an exterior region comprising the interstitial space between these spheres.Within the spheres, the electronic wave function is treated as atomic like, with an expansion in spherical harmonics up to ℓ max � 10.Within the interstitial space, the wave function is expanded in plane waves, with a maximum wave vector of k max .e maximum wave vector is determined by setting R mt k max � 8, where R mt is the radius of the smallest sphere.
e charge density is Fourier expanded up to a maximum wave vector of 14/a 0 , where a 0 is the Bohr radius.On-site Coulomb interaction was taken into account, using the same U and J values from the VASP total energy calculation.We performed self-consistent field calculation until energy convergence within 10 −4 Ry and charge convergence within 10 −3 e were achieved.e calculations were spin polarized, and we used the magnetic order in the ground state which was determined in the VASP calculation.e optical properties of the chromium chalcogenides were calculated using the wave functions obtained from the electronic calculation.Specifically, the imaginary part of the dielectric function ε 2 (ω) was obtained using the following momentum matrix elements [32]: where p is the momentum operator, |kn〉 is a crystal wave function with energy eigenvalue E kn , Zω is the photon energy, and f is the Fermi-Dirac distribution function.A Kramers-Kronig transformation was used to obtain the real part of the dielectric function ε 1 (ω), and the absorption coefficient α(ω) was calculated as

Results and Discussion
e energies of different magnetic structures for the chromium chalcogenides studied are listed in Table 1.Here, A-type antiferromagnetism (A-AFM) denotes ferromagnetism within the ac plane and antiferromagnetism between planes, C-type antiferromagnetism (C-AFM) denotes ferromagnetism between planes and antiferromagnetism within the plane, and G-type antiferromagnetism (G-AFM) denotes antiferromagnetism both within and between planes.Experimental results have shown that most of these materials, such as SbCrS 3 [19], SbCrSe 3 [18], and SiCrTe 3 [21], have magnetic order at low temperatures but are paramagnetic at room temperature.
e predicted magnetic structure of SbCrS 3 [19], SbCrSe 3 [19], LaCrS 3 [33], SbCrS 3 [33], and SiCrTe 3 [21] agree with experimental measurements.For LaCrSe 3 and CeCrSe 3 , there seems to be ambiguity as to whether the structure is antiferromagnetic or weakly ferromagnetic [34], and this is reflected in our calculations by the fact that the energy differences between ferromagnetic and antiferromagnetic states are very small.e bandgaps calculated using DFT + U for the chalcogenides are given in Table 2 and compared with experimental results.ese compounds range in character from metallic to semiconducting with small gaps, the largest of which is 1.29 eV.Due to the scarcity of experimental studies on these materials, information regarding the bandgaps of most materials is not available; however, we remark that the calculated bandgap values for SbCrSe 3 and SiCrTe 3 with our chosen values of U and J are in good agreement with experiment.
e calculated electronic band structure and density of states of the majority-spin carriers in GaCrSe 3 and SbCrS 3 are shown in Figures 4 and 5, respectively.ese structures were chosen because their bandgaps are suitable for photovoltaic applications.For GaCrSe 3 , we find a direct bandgap of 1.29 eV at Γ, while for SbCrSe 3 , we find a direct bandgap of 1.47 eV at Γ and an indirect bandgap of 0.98 eV corresponding to the Γ ⟶ T transition.In both of these structures, the density of states near both conduction band and valence band edges are primarily contributed by chromium d orbitals and chalcogenide p orbitals.In GaCrS 3 , the narrowness of the bands gives rise to high density of states, reaching a maximum of around 20 states/eV above the conduction band.SbCrSe 3 has a slightly smaller density of states, with a maximum of around 15 states/eV above the conduction band.
In addition to the density of states, we also calculated the optical properties for GaCrSe 3 and SbCrS 3 , as well as CrSiTe 3 and SbCrSe 3 , which were found to be narrow-gap semiconductors.
e absorption curves, shown in Figure 6, agree with the bandgap values calculated for these chalcogenides.From this figure, we see that they exhibit strong absorption ( > 10 5 cm −1 ) in the infrared and visible regions, suggesting that these materials may be of interest for photovoltaic applications.ese large values for the absorption coefficients are comparable to those measured in Zr-based chalcogenides [16].Additionally, absorption coefficients larger than 10 5 cm −1 are observed at significantly lower energies than that of pure silicon, which is the industry standard for photovoltaics.Silicon does not achieve such a value for the absorption coefficient until a photon energy of approximately 3 eV is reached [35]; on the contrary, we observe that SbCrS 3 and GaCrSe 3 , with Advances in Materials Science and Engineering comparable bandgaps to silicon, achieve this value at 1.5 eV and 2 eV, respectively.is would allow for better capture of light in the near-infrared and visible ranges.
e narrow bandgap SbCrSe 3 and SiCrTe 3 may be useful for infrared detection applications, while SbCrS 3 and GaCrSe 3 may be possible candidates for use as light absorbers in photovoltaics.

Conclusion
In conclusion, we have carried out density functional theory calculations on various chromium chalcogenide materials.We have demonstrated that calculations within the DFT + U scheme accurately predict the magnetic structure and energy gaps of this class of materials.We have identified two   Advances in Materials Science and Engineering

Figure 1 :Figure 2 :
Figure 1: Visualization of the unit cell of the NH 4 CdCl 3 prototype structure, which contains edge-sharing octahedra.

Figure 3 :
Figure 3: Visualization of the unit cell of SiCrTe 3 .

Figure 4 :Figure 5 :Figure 6 :
Figure 4: Calculated electronic band structure and density of states for GaCrSe 3 .

Table 1 :
Calculated energies (in meV per formula unit) of different magnetic states of the chromium chalcogenides.e energy of the lowest state for each material is set to zero for ease of comparison.For PCrSe 3 , α � c � 90 ∘ and β � 107.71 ∘ .For SiCrTe 3 , α � β � c � 50.56 ∘ .For SiCrTe 3 , each unit cell contains two Cr atoms; hence, there is only one possible antiferromagnetic state.

Table 2 :
Comparison of DFT + U bandgaps with experimental results.