Quaternary Selenides EuLnCuSe3: Synthesis, Structures, Properties and In Silico Studies

In this work, we report on the synthesis, in-depth crystal structure studies as well as optical and magnetic properties of newly synthesized heterometallic quaternary selenides of the Eu+2Ln+3Cu+1Se3 composition. Crystal structures of the obtained compounds were refined by the derivative difference minimization (DDM) method from the powder X-ray diffraction data. The structures are found to belong to orthorhombic space groups Pnma (structure type Ba2MnS3 for EuLaCuSe3 and structure type Eu2CuS3 for EuLnCuSe3, where Ln = Sm, Gd, Tb, Dy, Ho and Y) and Cmcm (structure type KZrCuS3 for EuLnCuSe3, where Ln = Tm, Yb and Lu). Space groups Pnma and Cmcm were delimited based on the tolerance factor t’, and vibrational spectroscopy additionally confirmed the formation of three structural types. With a decrease in the ionic radius of Ln3+ in the reported structures, the distortion of the (LnCuSe3) layers decreases, and a gradual formation of the more symmetric structure occurs in the sequence Ba2MnS3 → Eu2CuS3 → KZrCuS3. According to magnetic studies, compounds EuLnCuSe3 (Ln = Tb, Dy, Ho and Tm) each exhibit ferrimagnetic properties with transition temperatures ranging from 4.7 to 6.3 K. A negative magnetization effect is observed for compound EuHoCuSe3 at temperatures below 4.8 K. The magnetic properties of the discussed selenides and isostructural sulfides were compared. The direct optical band gaps for EuLnCuSe3, subtracted from the corresponding diffuse reflectance spectra, were found to be 1.87–2.09 eV. Deviation between experimental and calculated band gaps is ascribed to lower d states of Eu2+ in the crystal field of EuLnCuSe3, while anomalous narrowing of the band gap of EuYbCuSe3 is explained by the low-lying charge-transfer state. Ab initio calculations of the crystal structures, elastic properties and phonon spectra of the reported compounds were performed.

With all this in mind, in this work we directed our attention towards a series of quaternary selenides EuLnCuSe 3 (Ln = La, Sm, Gd, Tb, Dy, Ho, Y, Tm, Lu). Particularly, we focused on the efficient synthetic approach, crystal structure determination, optical and magnetic properties as well as computational analysis of the obtained structures. All the obtained results were examined with respect to the Ln 3+ radius.

Synthesis of Selenides
Quaternary selenides reported in this work were obtained from the corresponding starting materials, which were carefully prepared before the synthesis by thermolysis of the corresponding co-crystallized metal nitrates. Both Eu 2 O 3 and copper were pretreated as described in the literature [9]. Oxides of rare earth elements were annealed 1070 K to remove sorption water, hydroxides, hydrocarbonates and oxycarbonates [62,63], while Tb 4 O 7 was annealed at 470 K to keep its chemical composition and avoid a stepwise transformation: Tb 4 O 7 → Tb 11 O 20 → Tb 7 O 12 → Tb 2 O 3 [64]. The completeness of thermolysis was proved by the elimination of absorption bands in the corresponding IR spectra due to vibrations of the nitrate and hydroxyl groups. The absence of nitrogen-containing compounds was also confirmed by the distribution spectra of elements, plotted using an energy dispersive analysis system, in which no characteristic lines of nitrogen were observed, while the presence of characteristic lines of four elements O, Ln, Eu, Cu indicated the formation of oxides. According to X-ray phase analysis, the resulting samples comprised Cu 2 Ln 2 O 5 , CuO, (Ln/Eu) 2 O 3 , (Eu/Ln) 2 CuO 4 . Notably, the selenidation of complex oxides decreases the synthesis time of EuLnCuSe 3 in comparison with the selenidation of a mixture of commercial oxides CuO, Ln 2 O 3 , Eu 2 O 3 . Finally, compounds EuLnCuSe 3 in the powdered form were prepared by the reductive selenidation of the oxide mixtures, yielding aggregated particles with linear sizes up to about 30 µm ( Figure S1 in the Supplementary Materials). The energy dispersive X-ray spectrometric analysis data are in agreement with the calculated data (see the Section 3.2 Synthesis). According to the X-ray phase analysis, the content of the main phase in the samples ranges from 95.7% to 100% ( Figure S2 in the Supplementary Materials). Notably, the recently reported synthetic approach allowed the isolation of selenides with the main phase content of only 87.4-96.0% [60]. Furthermore, the described herein synthetic procedure significantly decreases the reaction time down to about 24 h at the temperature 970-1170 K.

Crystal Structures
The crystal structures of EuLnCuSe 3 were best solved in the orthorhombic system from powder X-ray diffraction patterns ( Figure 1, Table 1, Figure S2 and Tables S2-S5 in the Supplementary Materials). The calculated unit cell parameters obtained using the PBE0 density functional theory (DFT) functional are in good agreement with the experimental values ( Figure 2, Table 1 and Table S6 in the Supplementary Materials). Notably, the PBE0 functional was recently shown to be efficient to reproduce crystal structures, phonon spectra and elastic properties of compounds with ionic and covalent bonds [65]. Furthermore, it was also shown that the use of pseudopotentials for rare-earth ions replacing the inner shells, including 4f orbitals, allows the successful description of the structure and dynamics of the lattice of compounds with a sublattice of rare-earth ions [66]. [66].
Two types of orthorhombic Pnma and one type of orthorhombic Cmcm crystal structures were revealed for EuLnCuSe3 (Figures 1 and 2). Compound EuLaCuSe3 belongs to the structural type Ba2MnS3, while compounds EuLnCuSe3 (Ln = Sm, Eu [31], Gd, Tb, Dy, Ho, Y) are isostructural to Eu2CuS3. These crystal structures have the same type of symmetry but differ in the system of bonds and coordination of Ln 3+ . Finally, selenides EuLnCuSe3 (Ln = Tm, Yb, Lu) belong to the structural type KZrCuS3. In all compounds, the Eu 2+ and Ln 3+ cations each occupy two crystallographically independent positions.  Table S4 in the Supplementary Materials. Color code: magenta polyhedron = EuSe7 for EuLnCuSe3 (Ln = La-Y) and EuS6 for EuLnCuSe3 (Ln = Tm-Lu), green polyhedron = CuS4, blue polyhedron = LaSe7, purple polyhedron = LnSe6 (Ln = Sm-Y) and grey polyhedron = LnSe6 (Ln = Tm-Lu).
The distorted tetrahedra CuSe4 are linked by vertex atoms yielding polymeric chains along the b axis for the structural types Ba2MnS3 and Eu2CuS3, and along the a axis for the structural type KZrCuS3. In the structure of EuLaCuSe3, the polyhedron EuSe7 is formed from two atoms of Se1 i , two atoms of Se1 ii , two atoms of Se2 i and one atom of Se3 (Table  S4 in the Supplementary Materials) with all the Eu-Se distances being within the  Table S4 in the Supplementary Materials. Color code: magenta polyhedron = EuSe 7 for EuLnCuSe 3 (Ln = La-Y) and EuS 6 for EuLnCuSe 3 (Ln = Tm-Lu), green polyhedron = CuS 4 , blue polyhedron = LaSe 7 , purple polyhedron = LnSe 6 (Ln = Sm-Y) and grey polyhedron = LnSe 6 (Ln = Tm-Lu).
Two types of orthorhombic Pnma and one type of orthorhombic Cmcm crystal structures were revealed for EuLnCuSe 3 (Figures 1 and 2). Compound EuLaCuSe 3 belongs to the structural type Ba 2 MnS 3 , while compounds EuLnCuSe 3 (Ln = Sm, Eu [31], Gd, Tb, Dy, Ho, Y) are isostructural to Eu 2 CuS 3 . These crystal structures have the same type of symmetry but differ in the system of bonds and coordination of Ln 3+ . Finally, selenides EuLnCuSe 3 (Ln = Tm, Yb, Lu) belong to the structural type KZrCuS 3 . In all compounds, the Eu 2+ and Ln 3+ cations each occupy two crystallographically independent positions. eighth selenide ion is located at a distance that significantly exceeds the sum of ionic radii (>4 Å). The layered structures of EuLnCuSe3 (Ln = Eu-Lu) are formed by two-dimensional layers [LnCuSe3] in the ab plane for EuLnCuSe3 (Ln = Eu-Y) and in the ac plane for EuLnCuSe3 (Ln = Tm-Lu), between which the Eu 2+ ions are located.   The distorted tetrahedra CuSe 4 are linked by vertex atoms yielding polymeric chains along the b axis for the structural types Ba 2 MnS 3 and Eu 2 CuS 3 , and along the a axis for the structural type KZrCuS 3 . In the structure of EuLaCuSe 3 , the polyhedron EuSe 7 is formed from two atoms of Se1 i , two atoms of Se1 ii , two atoms of Se2 i and one atom of Se3 (Table S4 in the Supplementary Materials) with all the Eu-Se distances being within the calculated sum of ionic radii 3.15-3.18 Å (rSe 2− = 1.98 Å (coordination number = 6); rEu 2+ = 1.17 Å (coordination number = 6); rEu 2+ = 1.20 Å (coordination number = 7) [61]). The eighth selenide ion is located at a distance that significantly exceeds the sum of ionic radii (>4 Å). The layered structures of EuLnCuSe 3 (Ln = Eu-Lu) are formed by two-dimensional layers [LnCuSe 3 ] in the ab plane for EuLnCuSe 3 (Ln = Eu-Y) and in the ac plane for EuLnCuSe 3 (Ln = Tm-Lu), between which the Eu 2+ ions are located.
In the EuLnCuSe 3 , as the radius rLn 3+ decreases, the average Ln-Se distance gradually decreases from 3.09(7) Å down to 2.82(1) Å, while the Cu-Se distance remains intact and of 2.46(2)-2.48(3) Å (Table S4 in the Supplementary Materials). Notably, since the distorted octahedra LnSe 6 are linked into layers with the CuSe 4 tetrahedra, as rLn 3+ decreases, the crystal-chemical compression of these layers is observed. A decrease in the coordination saturation of the lanthanide cation leads to changes in the coordination polyhedron, structural type and space group ( Figure 1). Particularly, changing the coordination polyhedron from the monocapped trigonal prism LnSe 7 in EuLaCuSe 3 to the LnSe 6 octahedron (Ln = Eu-Y) leads to the change in the structural type Ba 2 MnS 3 to Eu 2 CuS 3 , and the monocapped trigonal prism EuSe 7 in EuLnCuSe 3 (Ln = Eu-Y) to the trigonal prism EuSe 6 in EuLnCuSe 3 (Ln = Tm-Lu) leads to the change in the structural type Eu 2 CuS 3 to KZrCuS 3 . A similar change in the structural types was observed for SrLnCuSe 3 (Ln = La-Lu) [26][27][28]36,42] and EuLnCuS 3 [10]. Thus, the change in the structural type in EuLnCuSe 3 is dictated by the size of the Ln 3+ cation.
In the structures of EuLnCuSe 3 (Ln = Sm, Gd, Tb, Dy, Ho, Y; structural type Eu 2 CuS 3 ), six Eu-Se distances, namely two Eu-Se1 i , two Eu-Se2 i and two Eu-Se3 i , are shorter than 3.2 Å, while the seventh and the longest distance Eu-Se3 ii gradually increases in the following sequence: 3.331 Å (EuSmCuSe 3 Table S4 in the Supplementary Materials). A similar trend was observed for isostructural compounds EuLnCuS 3 [10] and SrLnCuS 3 [67].
In the structures of EuLnCuSe 3 (Ln = Tm, Yb, Lu; structural type KZrCuS 3 ), six Eu-Se distances, namely four Eu-Se1 ii and two Eu-Se2 ii , are also shorter than 3.2 Å, while two more separations Eu-Se1 are longer than 3.55 Å and, thus, are not included in the coordination sphere of the metal ion ( Figure 1, Table S4 in the Supplementary Materials).
In compounds EuLnCuSe 3 (Ln = La-Lu), the unit cell parameters b (space group Pnma) and a (space group Cmcm) decrease. The parameters a (space group Pnma) and c (space group Cmcm) also decrease within each structural type of EuLnCuSe 3 (Ln = Sm-Ho) and EuLnCuSe 3 (Ln = Tm-Lu). The unit cell parameters a and c (space group Pnma) upon the changing of the structural type from Ba 2 MnS 3 to Eu 2 CuS 3 abruptly increase and decrease, respectively. The unit cell volume changes from EuHoCuSe 3 to EuTmCuSe 3 smoothly, despite the change in symmetry. A similar process was observed for isomorphic compounds ALnMX 3 (A = Sr, Eu, Ba; M = Cu, Ag; X = S, Se, Te) [25,29,68]. All these changes clearly testify to morphotropy in a series of the reported EuLnCuSe 3 compounds.
Despite that EuYCuSe 3 belongs to the same structural type as EuLnCuSe 3 (Ln = Sm-Ho), and the yttrium ion radius is less than that of holmium, the unit cell parameters of the yttrium-based selenide are larger than those for EuDyCuSe 3 and EuHoCuSe 3 , which is in agreement with the data obtained for isostructural quaternary sulfides ALnCuS 3 (A = Sr, Eu) [10,42,68,69]. This can be explained by the electronic structure of yttrium, which is a d-element.
In compounds EuLnCuSe3 (Ln = La-Lu), the unit cell parameters b (space group Pnma) and a (space group Cmcm) decrease. The parameters a (space group Pnma) and c (space group Cmcm) also decrease within each structural type of EuLnCuSe3 (Ln = Sm-Ho) and EuLnCuSe3 (Ln = Tm-Lu). The unit cell parameters a and c (space group Pnma) upon the changing of the structural type from Ba2MnS3 to Eu2CuS3 abruptly increase and In general, the EuLnCuSe 3 compounds are isostructural to EuLnCuS 3 [10]. When selenium is included instead of sulfur, the anionic sublattice expands, since selenium is a larger anion. A comparison of sulfide and selenide analogues EuLnCuX 3 (X = S, Se, Ln = La-Lu, Y) allows the reveal of the formation of different types of crystal structures with the same rare-earth metals ( Figure 4). Particularly, SrYCuS 3 and SrHoCuS 3 crystallize in space group Pnma, while their selenide derivatives crystallize in space group Cmcm ( Figure 4) [28,42,53,69].
less pronounced ionic character of atoms, with the mm2 symmetry (less symmetric coordination). As the ionic radius of the lanthanide rLn 3+ decreases in ALnCuSe3, the formation of a highly symmetric coordination of A 2+ occurs earlier in SrLnCuSe3 than in EuLnCuSe3. Particularly, SrSe7 (Pnma) transforms into SrSe6 (Cmcm) from Dy to Ho [42], while the change in the coordination polyhedron of europium occurs later, viz. from Y to Er. A similar process was an observed trend for ALnCuS3 (A = Eu, Sr) [10,42]. In SrLnCuS3 and SrLnCuS3, for example, the change in the strontium and europium coordination polyhedra is observed from Y to Er and from Er to Tm, respectively [10,42]. This is due to the more pronounced ionic characteristics of the alkaline earth elements in comparison with the rare earth metals. Thus, with an increase in the degree of ionicity of the A 2+ -X bond, as well as the A 2+ radius, the number of compounds, crystallized in the Cmcm space group, increases, and the A 2+ ions become six-coordinated earlier.
To separate orthorhombic space groups Cmcm and Pnma in sulfides ALnCuS3 (A = Sr, Eu, Ba) clearly, the tolerance factor t′ = [rA 2+ ] × [rCu + ] + a × [rLn 3+ ] 2 (a = −0.6, critical value of t′ = 0.28) was applied [10,70]. The experimental data on selenides ALnCuSe3 obtained here allow for checking the formula for these compounds. Considering the data on both sulfides and selenides, the adjusted value of the tolerance factor is 0.273. A slight change in the value of t′ confirms the previous assumption that the change in the space group of compounds ALnCuX3 is almost independent on the ionic radius of the chalcogen, but depends on the ionic radii of A 2+ , Ln 3+ and Cu + [70]. We assume that studying new  In the structures of SrYCuSe 3 and EuYCuSe 3 , the Sr 2+ and Eu 2+ ions exhibit different coordination polyhedra. Compound SrYCuSe 3 contains the SrSe 6 polyhedron with the most pronounced ionic character of the included atoms, having a local symmetry of 43m (highly symmetric coordination), while EuYCuSe 3 contains the EuSe 7 polyhedron with a less pronounced ionic character of atoms, with the mm2 symmetry (less symmetric coordination). As the ionic radius of the lanthanide rLn 3+ decreases in ALnCuSe 3 , the formation of a highly symmetric coordination of A 2+ occurs earlier in SrLnCuSe 3 than in EuLnCuSe 3 . Particularly, SrSe 7 (Pnma) transforms into SrSe 6 (Cmcm) from Dy to Ho [42], while the change in the coordination polyhedron of europium occurs later, viz. from Y to Er. A similar process was an observed trend for ALnCuS 3 (A = Eu, Sr) [10,42]. In SrLnCuS 3 and SrLnCuS 3 , for example, the change in the strontium and europium coordination polyhedra is observed from Y to Er and from Er to Tm, respectively [10,42]. This is due to the more pronounced ionic characteristics of the alkaline earth elements in comparison with the rare earth metals.
Thus, with an increase in the degree of ionicity of the A 2+ -X bond, as well as the A 2+ radius, the number of compounds, crystallized in the Cmcm space group, increases, and the A 2+ ions become six-coordinated earlier.
To separate orthorhombic space groups Cmcm and Pnma in sulfides ALnCuS 3 (A = Sr, Eu, Ba) clearly, the tolerance factor t = [rA 2+ ] × [rCu + ] + a × [rLn 3+ ] 2 (a = −0.6, critical value of t = 0.28) was applied [10,70]. The experimental data on selenides ALnCuSe 3 obtained here allow for checking the formula for these compounds. Considering the data on both sulfides and selenides, the adjusted value of the tolerance factor is 0.273. A slight change in the value of t confirms the previous assumption that the change in the space group of compounds ALnCuX 3 is almost independent on the ionic radius of the chalcogen, but depends on the ionic radii of A 2+ , Ln 3+ and Cu + [70]. We assume that studying new compounds ALnCuX 3 (X = S, Se; A is an element with an ionic radius in the range of 1.2 Å (rEu 2+ ) < (rA 2+ ) < 1.35 Å (rBa 2+ )), the value of t can vary in the range from 0.273 to 0.285.
Thus, our results on the crystal structures of the reported quaternary selenides are in full agreement with the previously established structural types for compounds of the type A 2+ Ln 3+ Cu +1 X 3 (A = Eu, Sr, Pb, Ba; X = S, Se) [8]. Furthermore, the formation of the most symmetric structure of the type KZrCuS 3 with the decreasing of rLn 3+ and decreasing of the distortion of two-dimensional [LnCuSe 3 ] layers in a series of the discussed compounds is also confirmed.
In compounds EuLnCuS 3 , the europium ion is Eu 2+ [11,68,71]. An insignificant structural difference also anticipates a similarity of charges in EuLnCuX 3 (X = S, Se). The bond valence sum calculation is frequently used as an efficient tool to estimate the oxidation states of atoms.
The oxidation states of the Eu, Ln and Cu ions in EuLnCuSe 3 are close to 2, 3 and 1, respectively ( Table S7 in the Supplementary Materials). This is in agreement with the data for Eu 2 CuSe 3 on the presence of both Eu 2+ and Eu 3+ [31]. Moreover, the Eu 3+ ion occupies the same position in Eu 2+ Eu 3+ CuSe 3 as Ln 3+ in EuLnCuSe 3 .

Elastic Properties
Elastic constants of crystals EuLnCuSe 3 were calculated using the PBE0 functional (Table S8 in Table S8 in the Supplementary Materials), for example, for the former two structural types C 12 < C 13 , while C 12 > C 13 for the latter structural type ( compounds ALnCuX3 (X = S, Se; A is an element with an ionic radius in the range of 1.2 Å (rEu 2+ ) < (rA 2+ ) < 1.35 Å (rBa 2+ )), the value of t′ can vary in the range from 0.273 to 0.285. Thus, our results on the crystal structures of the reported quaternary selenides are in full agreement with the previously established structural types for compounds of the type A 2+ Ln 3+ Cu +1 Х3 (A = Eu, Sr, Pb, Ba; X = S, Se) [8]. Furthermore, the formation of the most symmetric structure of the type KZrCuS3 with the decreasing of rLn 3+ and decreasing of the distortion of two-dimensional [LnCuSe3] layers in a series of the discussed compounds is also confirmed.
In compounds EuLnCuS3, the europium ion is Eu 2+ [11,68,71]. An insignificant structural difference also anticipates a similarity of charges in EuLnCuX3 (X = S, Se). The bond valence sum calculation is frequently used as an efficient tool to estimate the oxidation states of atoms.
The oxidation states of the Eu, Ln and Cu ions in EuLnCuSe3 are close to 2, 3 and 1, respectively (Table S7 in the Supplementary Materials). This is in agreement with the data for Eu2CuSe3 on the presence of both Eu 2+ and Eu 3+ [31]. Moreover, the Eu 3+ ion occupies the same position in Eu 2+ Eu 3+ CuSe3 as Ln 3+ in EuLnCuSe3.

Elastic Properties
Elastic constants of crystals EuLnCuSe3 were calculated using the PBE0 functional (Table S8 in the Supplementary Materials), while the bulk modulus, Young's modulus, and shear modulus were calculated in the Voigt, Reuss and Hill approximations (Table S9 in the Supplementary Materials). The dependence of the Young's modulus on the direction demonstrates anisotropy of the elastic properties ( Figure 5). Calculations predict that the elastic properties of crystals of EuLnCuSe3 (Ln = La-Lu) strongly depend on the structural types Ba2MnS3, Eu2CuS3 and KZrCuS3 ( Figure 5, Table S8 in the Supplementary Materials), for example, for the former two structural types C12 < C13, while C12 > C13 for the latter structural type (Table S8 in   The calculated values of the shear modulus and bulk modulus allows the estimation of the Vickers hardness for SrTmCuS3 according to empirical formula HV = 0.92(G/B) 1.137 × G 0.708 , where G is the shear modulus, and B is the bulk modulus, obtained from the Hill approximation [72]. It was found that for the structures with space group Pnma, the calculated hardness decreases from La to Ho, while for the structures with space group Cmcm, the calculated hardness increases from Er to Lu (Table 2). Furthermore, using formula A U = 5(GV/GR) + BV/BR − 6 (GV and GR, BV and RR are the shear and bulk modulus, calculated in the Voigt and Reuss approximations, respectively), we also calculated the so-called "universal" anisotropy index A U to estimate anisotropy of the elastic properties of EuLnCuSe3 [73]. Deviation of the index A U from zero determines the degree of anisotropic properties of the crystal. It was found that for the crystals of EuLnCuSe3 with space group Pnma, the value of A U increases, while for the crystals with space group Cmcm, the same value decreases (Table 2). The calculated values of the shear modulus and bulk modulus allows the estimation of the Vickers hardness for SrTmCuS 3 according to empirical formula H V = 0.92(G/B) 1.137 × G 0.708 , where G is the shear modulus, and B is the bulk modulus, obtained from the Hill approximation [72]. It was found that for the structures with space group Pnma, the calculated hardness decreases from La to Ho, while for the structures with space group Cmcm, the calculated hardness increases from Er to Lu (Table 2). Furthermore, using formula A U = 5(G V /G R ) + B V /B R − 6 (G V and G R , B V and R R are the shear and bulk modulus, calculated in the Voigt and Reuss approximations, respectively), we also calculated the so-called "universal" anisotropy index A U to estimate anisotropy of the elastic properties of EuLnCuSe 3 [73]. Deviation of the index A U from zero determines the degree of anisotropic properties of the crystal. It was found that for the crystals of EuLnCuSe 3 with space group Pnma, the value of A U increases, while for the crystals with space group Cmcm, the same value decreases (Table 2).

Phonon, Raman and IR Spectra
We also applied ab initio calculations with the PBE0 functional to examine wavenumbers and types of the phonon modes. As a result, IR active modes, Raman modes and "silent" modes were revealed (Tables S10-S12 in the Supplementary Materials). Participation of each ion in a particular mode was estimated from the analysis of displacement vectors obtained from these ab initio calculations ( Figure 6 and Tables S13-S15 in the Supplementary Materials).

Phonon, Raman and IR Spectra
We also applied ab initio calculations with the PBE0 functional to examine wavenumbers and types of the phonon modes. As a result, IR active modes, Raman modes and "silent" modes were revealed (Tables S10-S12 in the Supplementary Materials). Participation of each ion in a particular mode was estimated from the analysis of displacement vectors obtained from these ab initio calculations ( Figure 6 and Tables S13-S15 in the Supplementary Materials).

EuLaCuSe3 structural type Ba2MnS3
EuTbCuSe3 structural type Eu2CuS3 EuTmCuSe3 structural type KZrCuS3 According to the obtained results, the phonon wavenumbers at the Γ-point do not exceed 230 cm −1 . In this range, not only light ions of selenium and copper but also heavy rare-earth ions are actively involved. In crystals of EuLnCuSe 3 with a structural type of Ba 2 MnS 3 , the europium ions actively participate in modes with wavenumbers up to about 150 cm −1 , while for the structural types of Eu 2 CuS 3 and KZrCuS 3 , the europium ions participate in modes with wavenumbers up to about 110 cm −1 , respectively ( Figure 6).
Both the Ln 2+ and selenium ions are involved in almost all modes for all the reported structural types ( Figure 6).
Unfortunately, contrarily to the corresponding sulfide analogues [9,11], we faced problems collecting Raman spectra for some of the reported selenides, and our numerous attempts failed. This can be explained by the narrower bandgap and sequent high absorption at the wavelengths of excitation and of the scattered Raman signal. As a result, we were able to record the Raman spectra for EuTbCuSe 3 and EuTmCuSe 3 (Figure 7). Both the experimental and calculated Raman spectra, obtained using the PBE0 functional, for three typical representatives of each structural types of the reported selenides, namely EuLaCuSe 3 (structural type Ba 2 MnS 3 ), EuTbCuSe 3 (structural type Eu 2 CuS 3 ) and EuTmCuSe 3 (structural type KZrCuS 3 ), exhibit bands exclusively in the region up to about 250 cm −1 (Figure 7). Notably, the experimental and calculated spectra of the latter two compounds are pairwise very similar with the most intensive bands at about 60 and 180 cm −1 , each corresponding to active A g modes. However, the spectrum of EuTbCuSe 3 contains an additional clearly visible band at about 15 cm −1 also due to the A g mode (Figure 7). Both spectra further exhibit bands from B 2g and B 3g modes. Notably, the experimental spectrum of EuTmCuSe 3 exhibits an additional clearly defined band at about 45 cm −1 (Figure 7). We can tentatively assign the origin of this band to the Raman active mode of the Tm 4 Se 3 O 4 impurity. Unfortunately, our numerous attempts to calculate the Raman spectrum of Tm 4 Se 3 O 4 failed and was found to be highly time-consuming, which can be explained by the nonstoichiometric occupancy of the selenium ions.
150 cm −1 , while for the structural types of Eu2CuS3 and KZrCuS3, the europium ions participate in modes with wavenumbers up to about 110 cm −1 , respectively ( Figure 6). Both the Ln 2+ and selenium ions are involved in almost all modes for all the reported structural types ( Figure 6).
Unfortunately, contrarily to the corresponding sulfide analogues [9,11], we faced problems collecting Raman spectra for some of the reported selenides, and our numerous attempts failed. This can be explained by the narrower bandgap and sequent high absorption at the wavelengths of excitation and of the scattered Raman signal. As a result, we were able to record the Raman spectra for EuTbCuSe3 and EuTmCuSe3 (Figure 7). Both the experimental and calculated Raman spectra, obtained using the PBE0 functional, for three typical representatives of each structural types of the reported selenides, namely EuLaCuSe3 (structural type Ba2MnS3), EuTbCuSe3 (structural type Eu2CuS3) and EuTmCuSe3 (structural type KZrCuS3), exhibit bands exclusively in the region up to about 250 cm −1 (Figure 7). Notably, the experimental and calculated spectra of the latter two compounds are pairwise very similar with the most intensive bands at about 60 and 180 cm −1 , each corresponding to active Ag modes. However, the spectrum of EuTbCuSe3 contains an additional clearly visible band at about 15 cm −1 also due to the Ag mode ( Figure  7). Both spectra further exhibit bands from B2g and B3g modes. Notably, the experimental spectrum of EuTmCuSe3 exhibits an additional clearly defined band at about 45 cm −1 (Figure 7). We can tentatively assign the origin of this band to the Raman active mode of the Tm4Se3O4 impurity. Unfortunately, our numerous attempts to calculate the Raman spectrum of Tm4Se3O4 failed and was found to be highly time-consuming, which can be explained by the nonstoichiometric occupancy of the selenium ions.
Interestingly, the calculated spectrum of EuLaCuSe3 is characterized by a noticeably larger number of bands ( Figure 7). An overwhelming majority of these bands are due to Ag modes, with the most intensive one at about 150 cm −1 , accompanied with two bands at about 60 and 90 cm −1 , arising from B3g and B1g modes, respectively, and one band at about 115 cm −1 , corresponding to the A1g mode (Figure 7). The calculate (black) and experimental (red) Raman spectra for the selected selenides of three different structural types. Calculations were performed at λexc = 514.5 nm and T = 300 K for EuLaCuSe3, and for λexc = 532 nm and T = 296 K for EuTbCuSe3 (Pnma) and EuTmCuSe3 (Cmcm), respectively. The wavenumber scale of the calculated spectra was multiplied by 0.9 to compensate the overestimated results obtained by the applied PBE0 functional.
Structures of EuLnCuSe3 were also studied by IR spectroscopy. Notably, all the reported selenides were found to be transparent for IR radiation in the range of 250-4000 cm −1 . However, a close inspection of the range of 85-250 cm −1 revealed a set of absorption bands in the IR spectra of compounds ( Figure 8).
Mechanical representation in the center of the Brillouin zone for EuLnCuSe3 with the structural types Ba2MnS3 and Eu2CuS3 is described as Γvibr = 12Ag + 6Au + 6B1g + 12B1u + 12B2g + 6B2u + 6B3g + 12B3u. Mechanical representation for the structural type KZrCuS3 is the sum   Figure 7. The calculate (black) and experimental (red) Raman spectra for the selected selenides of three different structural types. Calculations were performed at λ exc = 514.5 nm and T = 300 K for EuLaCuSe 3 , and for λ exc = 532 nm and T = 296 K for EuTbCuSe 3 (Pnma) and EuTmCuSe 3 (Cmcm), respectively. The wavenumber scale of the calculated spectra was multiplied by 0.9 to compensate the overestimated results obtained by the applied PBE0 functional.
Interestingly, the calculated spectrum of EuLaCuSe 3 is characterized by a noticeably larger number of bands (Figure 7). An overwhelming majority of these bands are due to A g modes, with the most intensive one at about 150 cm −1 , accompanied with two bands at about 60 and 90 cm −1 , arising from B 3g and B 1g modes, respectively, and one band at about 115 cm −1 , corresponding to the A 1g mode (Figure 7).
The structures of EuLnCuSe 3 were also studied by IR spectroscopy. Notably, all the reported selenides were found to be transparent for IR radiation in the range of 250-4000 cm −1 . However, a close inspection of the range of 85-250 cm −1 revealed a set of absorption bands in the IR spectra of compounds ( Figure 8). experimental IR spectrum is, in general, the result of overlapping of closely positioned bands without any large gaps between the IR active vibrations (Figure 8). It should also be noted that all the reported selenides, except for the La-based one, are contaminated by the Ln4O4Se3 [74] or Ln2SeO2 [75] impurities (see discussion above). Thus, additional bands from the Ln-Se vibrations of the oxyselenides are also expected in the same range as for the quaternary selenides. However, their impact is negligible as evidenced from very weak bands of the oxyselenides Ln-O vibrations found at about 300-600 cm −1 [76].  Mechanical representation in the center of the Brillouin zone for EuLnCuSe 3 with the structural types Ba 2 MnS 3 and Eu 2 CuS 3 is described as Γ vibr = 12A g + 6A u + 6B 1g + 12B 1u + 12B 2g + 6B 2u + 6B 3g + 12B 3u . Mechanical representation for the structural type KZrCuS 3 is the sum of 36 irreducible representations: 5A g + 2A u + 4B 1g + 7B 1u + B 2g + 7B 2u + 5B 3g + 5B 3u . Notably, the g and u labeled modes are Raman and IR active, respectively, except for the "silent" A u modes. Thus, the number of IR active modes is 16 and 27 for the structural types KZrCuS 3 , and Eu 2 CuS 3 and Ba 2 MnS 3 , respectively.
The experimental IR spectra of the discussed selenides are of three different profiles and in good agreement with the calculated ones ( Figure 8). Particularly, the spectra of the structural types KZrCuS 3 and Eu 2 CuS 3 (except for EuYCuSe 3 ) are very similar, and a set of bands is observed at about 80-210 cm −1 . This is also supported by the positions of the calculated IR active bands for these selenides, although the total number of these bands differs significantly (Figure 8). The experimental IR spectrum of the yttrium-based selenide exhibits bands in a broader range at about 80-240 cm −1 . Although the structure of this compound also belongs to the structural type Eu 2 CuS 3 as for the Sm-, Gd-, Tb-, Dy-and Hobased derivatives, again the different nature of yttrium (d-element) is reflected. EuLaCuSe 3 is the only described herein selenide of the structural type Ba 2 MnS 3 . Its experimental IR spectrum is, in general, the result of overlapping of closely positioned bands without any large gaps between the IR active vibrations (Figure 8). It should also be noted that all the reported selenides, except for the La-based one, are contaminated by the Ln 4 O 4 Se 3 [74] or Ln 2 SeO 2 [75] impurities (see discussion above). Thus, additional bands from the Ln-Se vibrations of the oxyselenides are also expected in the same range as for the quaternary selenides. However, their impact is negligible as evidenced from very weak bands of the oxyselenides Ln-O vibrations found at about 300-600 cm −1 [76].
It is known that, for the same compound, hybrid functionals yield a higher value of the gap width than non-hybrid ones. For example, it was shown that the hybrid functional PBE0 overestimates, while the non-hybrid functional PBE underestimates the gap in some oxides [65]. The same is true for the discussed selenides EuLnCuSe 3 (Table 2). Furthermore, the non-hybrid PBE functional also yielded an underestimated value of the gap width for selenides as reported earlier (Table 2 and Figure S3 in the Supplementary Materials) [32].
The PBE-based calculation performed in this revealed the gap width close to the results of the PBE calculations reported earlier (Table 2), where the f -shell of the rare-earth ions was considered explicitly and was not replaced by a pseudopotential [32]. We also performed a test calculation with the PBE functional for EuLaCuSe 3 , where the f -electrons of europium were considered explicitly. The obtained gap of about 1 eV is in good agreement with the reported one [32] and is close to the results, which were obtained by replacing the f -shell of europium with a pseudopotential.
The experimental band gaps of EuLnCuSe 3 were determined from the Kubelka-Munk functions ( Figure S4 in the Supplementary Materials). A recent study suggested that use of the simple Kubelka-Munk function may lead to errors in band gap values of the order of several tenths of eV [77]. Landi et al. argue that the Tauc-type modification of the Kubelka-Munk functions, though containing approximations, leads to more correct band gap values. Therefore, we tested both uncorrected ( Figure S4 in the Supplementary Materials) and modified Kubelka-Munk functions ( Table 2 and Figure S3 in the Supplementary Materials). For some samples, the Kubelka-Munk functions modified for an indirect band gap failed, and the corresponding cells in Table 2 were left empty. Most of the unmodified Kubelka-Munk spectra exhibit a shoulder at longer wavelengths preceding the onset of an evident fundamental absorption. This shoulder is associated with oxyselenide impurities ( Table 1) that are known to be semiconductors with the band gap of 1.02-1.46 eV [74]. Thus, to estimate the band gap value of the discussed selenides, an additional correction was applied with respect to impurities [78]. For the Dy-and Sm-based selenides, such a correction was less than 0.01 eV, while no correction was applied for the Y-based derivative since additional data on impurities are required. For the Gd-, Ho-, Tb-, Lu-and Yb-derived compounds, a correction of several hundredths of eV was revealed for the samples containing about 1% of impurities (Table 2).
It is known that, for the same compound, hybrid functionals yield a higher value of the gap width than non-hybrid ones. For example, it was shown that the hybrid functional PBE0 overestimates, while the non-hybrid functional PBE underestimates the gap in some oxides [65]. The same is true for the discussed selenides EuLnCuSe3 (Table 2). Furthermore, the non-hybrid PBE functional also yielded an underestimated value of the gap width for selenides as reported earlier (Table 2 and Figure S3 in the Supplementary Materials) [32]. The PBE-based calculation performed in this revealed the gap width close to the results of the PBE calculations reported earlier (Table 2), where the f-shell of the rare-earth ions was considered explicitly and was not replaced by a pseudopotential [32]. We also performed a test calculation with the PBE functional for EuLaCuSe3, where the f-electrons of europium were considered explicitly. The obtained gap of about 1 eV is in good agreement with the reported one [32] and is close to the results, which were obtained by replacing the f-shell of europium with a pseudopotential.
The experimental band gaps of EuLnCuSe3 were determined from the Kubelka-Munk functions ( Figure S4 in the Supplementary Materials). A recent study suggested that use of the simple Kubelka-Munk function may lead to errors in band gap values of The obtained values of the band gap are close to those obtained for isostructural sulfides EuLnCuS 3 ( Figure S4 in the Supplementary Materials). The band gap of BaLnCuS 3 is much higher than that of EuLnCuX 3 (X = S, Se) [70,79].
Narrowing of the band gap is well-known for the majority of ab initio calculations using the PBE functional. We applied the B3LYP and PBE0 functionals to calculate the electronic band structure, which yielded overestimated band gap values in comparison to the experimental ones, while calculations preformed earlier using the PBE functional and considering f -electrons, in general, underestimated the band gap values [32].
The most probable ion that could be responsible for this narrowing is Eu 2+ , which is present in all of the discussed compounds. We suggest that interband transitions between sub-bands originating from partially occupied orbitals of f -electrons of Eu 2+ to, e.g., unoccupied bands originating from d-orbitals of the same ion, are not correctly described by a pseudopotential, and these sub-bands are the factor that determines experimental behavior of band gaps ( Figure S3 in the Supplementary Materials). More specifically, experimental values of the direct band gap in all selenides containing divalent europium and trivalent lanthanide with f -electrons vary from 1.87 eV to 1.97 eV for the Pnma space group, and from 2.02 eV to 2.1 eV for the Cmcm space group (Tm, Lu).
The B3LYP-based calculation with the europium pseudopotential predicts direct band gaps in the range 2.38-2.63 eV and 2.52-2.59 eV for Pnma and Cmcm space groups, respectively. An indirect band gap for the Cmcm space group is predicted to be narrower than the direct one in contrast to the Pnma space group. However, for consistency, we compare direct band gaps produced by stronger direct transitions for all crystals. The discrepancy between the B3LYP-based calculations and experimental values of the order of 0.6 eV is ascribed to a lower position of excited d-bands of europium in comparison with the B3LYP calculation. This explanation is valid for all the discussed compounds except for the La-, Yand Yb-derived ones, which exhibit much narrower experimental band gaps. For the Yb containing selenide, the narrower band gap can be explained as similar to BaYbCuS 3 [79]. The Yb 3+ ion with the f 13 -electronic configuration is characterized by a low-lying charge transfer transition to the f 14 -configuration with enhanced stability. In selenides, this charge transfer transition is expected to lie at even lower energies than in sulfides, resulting in additional narrowing of the EuYbCuSe 3 band gap with respect to calculations. To explain the narrowing of the band gap of compounds where Ln 3+ completely lack f -electrons, i.e., La and Y are more challenging, one may expect that the lack of f -electrons results in certain structural changes that lead to enhancement of the crystal field acting upon europium ions and sequent narrowing of the band gap. Let us compare, e.g., a local environment of the europium ion in EuYCuSe 3 and in EuHoCuSe 3 with the ionic radius of Ln 3+ closest to that of Y 3+ . The europium ion is coordinated by seven selenium ions at distances of the order of 3 Å. Specifically, one of the selenium ions is positioned at a much longer distance of about 3.4-3.5 Å for the Ho-and Y-based selenides, respectively. Thus, asymmetry of the crystal field for EuYCuSe 3 is the largest one. Hypothetically, it can be the reason for the narrowing of the EuYCuSe 3 band gap. However, since the increase of asymmetry is rather small, this explanation must be considered as a hypothesis.

Magnetic Properties
The external field dependent magnetic moments of the samples at room temperature exhibit a linear plot, corresponding to the Curie law for paramagnets and differ only in the slope of the line (Figure 10 and Figure S5 in the Supplementary Materials).
The calculated values of µ and C for EuYCuSe 3 , EuLaCuSe 3 and EuLuCuSe 3 are the same since the cations Y 3+ , La 3+ and Lu 3+ are diamagnetic, and the final magnetic properties of these selenides are exclusively associated with the Eu 2+ cation. The experimental values of the effective magnetic moments of EuYCuSe 3 and EuLaCuSe 3 are equal within the error and differ from the calculated ones by about 1%, while the experimental magnetic moment of EuLuCuSe 3 is about 5.5% larger in comparison to the calculated one (Table 3). For the other selenides, the relative deviations of the calculated values of the magnetic moments from the calculated ones do not exceed 3%. As for the corresponding Curie constants, this value of deviation is as twice as large due to the quadratic dependence on the values of the effective magnetic moments. Figure 10. Field-dependent magnetic moments at 296 K (left), and temperature-dependent specific magnetization (middle) and reciprocal magnetic susceptibility (right) of EuTmCuSe3 (top row), EuLaCuSe3 (middle row) and EuHoCuSe3 (bottom row), respectively, at 10 Oe. The measurements of low-temperature magnetization were performed in the zero-field cooled (ZFC) and nonzero-field cooled (FC) modes. Fit line shows approximation by the Néel model of a ferrimagnet [80].
The calculated values of μ and C for EuYCuSe3, EuLaCuSe3 and EuLuCuSe3 are the same since the cations Y 3+ , La 3+ and Lu 3+ are diamagnetic, and the final magnetic properties of these selenides are exclusively associated with the Eu 2+ cation. The experimental values of the effective magnetic moments of EuYCuSe3 and EuLaCuSe3 are equal within the error and differ from the calculated ones by about 1%, while the experimental magnetic moment of EuLuCuSe3 is about 5.5% larger in comparison to the calculated one (Table 3). For the other selenides, the relative deviations of the calculated values of the magnetic moments from the calculated ones do not exceed 3%. As for the corresponding Curie constants, this value of deviation is as twice as large due to the quadratic dependence on the values of the effective magnetic moments.
Plots of both specific magnetization and reciprocal magnetic susceptibility of EuLaCuSe3 and EuLuCuSe3 at 4.2-20 K indicate that, upon cooling down to 4.2 K, the compounds remain paramagnetic with the corresponding (asymptotic) Curie temperature θp = 0.2 and 3.1 K, respectively ( Figure 10 and Figure S5 in the Supplementary  Figure 10. Field-dependent magnetic moments at 296 K (left), and temperature-dependent specific magnetization (middle) and reciprocal magnetic susceptibility (right) of EuTmCuSe 3 (top row), EuLaCuSe 3 (middle row) and EuHoCuSe 3 (bottom row), respectively, at 10 Oe. The measurements of low-temperature magnetization were performed in the zero-field cooled (ZFC) and nonzero-field cooled (FC) modes. Fit line shows approximation by the Néel model of a ferrimagnet [80]. Plots of both specific magnetization and reciprocal magnetic susceptibility of EuLaCuSe 3 and EuLuCuSe 3 at 4.2-20 K indicate that, upon cooling down to 4.2 K, the compounds remain paramagnetic with the corresponding (asymptotic) Curie temperature θ p = 0.2 and 3.1 K, respectively ( Figure 10 and Figure S5 in the Supplementary Materials). At these temperatures, the samples are most likely to be ferromagnetic. This assumption is supported by the presence of one type of magnetic ions in the compounds (Eu 2+ ), the coincidence of the data for the FC and ZFC measurements and the linearity of the temperature dependent reciprocal susceptibility down to the lowest temperatures.
The magnetic behavior of EuYCuSe 3 is similar to that of EuLaCuSe 3 and EuLuCuSe 3 . However, below 12 K, the FC and ZFC curves slightly diverge, although there are no noticeable points in these dependences ( Figure S5 in the Supplementary Materials). The paramagnetic Curie temperatures were found to be 4.2 and 3.3 K for FC and ZFC, respectively. The second value is more reliable since the dependence in this case is more linear. Probably, a transition to a ferromagnetic state occurs at this point. The thermomagnetic dependences for EuYbCuSe 3 are similar with the FC and ZFC curves being less diverged, and the transition to a ferromagnetic state occurs, most likely, very close to 4 K.
Curves for EuTbCuSe 3 , EuDyCuSe 3 (ZFC), EuHoCuSe 3 and, to a lesser extent, EuTmCuSe 3 with decreasing temperatures exhibit an abrupt drop in the linear dependence of the reciprocal susceptibility, characteristic of ferrimagnets ( Figure 10 and Figure S5 in the Supplementary Materials). Paramagnetic Curie temperatures, calculated from linear sections of the curves, are negative except for EuTmCuSe 3 . This testifies to the antiparallel coordination of the magnetic sublattices. The approximation of these curves at temperatures above the points of local minima χ −1 , according to the Néel theory for ferrimagnets, reveals the data, including the Néel temperatures, gathered in Table 4.  1 Formula applied for the model: 1 The decreasing of the χ −1 values of EuDyCuSe 3 and EuTbCuSe 3 in the FC mode in comparison to the ZFC mode at low temperatures corresponds to an increase in the magnetic moment of the sample cooled in the presence of an external field. Apparently, when heated above the Néel point, the effect of this initial magnetization persists up to some higher temperatures. A similar explanation can be proposed for the divergence of the FC and ZFC curves of EuYCuSe 3 and EuYbCuSe 3 , discussed above, although the temperature did not decrease below the points of magnetic transitions. The positive value of the paramagnetic Curie point calculated for EuTmCuSe 3 contradicts the conclusion about ferrimagnetic ordering, but this value is rather close to zero.
The magnetization of EuHoCuSe 3 below the Néel point looks significantly different, and its value is negative at temperatures from 4.2 to 4.8 K. This recently attracted increased interest for possible practical applications [81,82]. Several mechanisms of this phenomenon are considered: negative exchange coupling among ferromagnetic sublattices, negative exchange coupling among canted antiferromagnetic sublattices, negative exchange coupling among ferromagnetic/canted antiferromagnetic and paramagnetic sublattices, imbalance of spin and orbital moments, interfacial exchange coupling between ferromagnetic and antiferromagnetic phases [81]. Similar temperature dependences of the reciprocal suscep-tibility were obtained for (Tm 0.8 Mn 0.2 )MnO 3 [82]. It was established that the negative magnetization is determined by the first mechanism. Most likely, EuHoCuSe 3 exhibits similar properties and belongs to the N-type ferrimagnets according to Néel [80].

Synthesis
Powdered samples of EuLnCuSe 3 were prepared by reductive selenidation of the oxide mixtures produced by thermolysis of the co-crystallized metal nitrates. Metallic copper fragments were mechanically cleaned, treated with alcohol and dried at room temperature. Oxides of rare earth elements were annealed in a quartz glass in a muffle furnace at 1070 K for 5 h, while Tb 4 O 7 was annealed at 470 K to keep its chemical composition and avoid stepwise transformation.
The parent reagents were calculated based on the stoichiometric ratio of metals in EuLnCuSe 3 (1:1:1). Starting materials were dissolved in 65% aqueous HNO 3 (100 mL) under vigorous stirring upon heating to 330-350 K. The resulting solution was slowly evaporated, followed by the co-crystallization of nitrates. The dry residue was transferred into a quartz glass and subjected to thermolysis at 1070 K for 5 h. The baked mass was ground to a size of about 30-50 µm.
Reductive selenidation was carried out in a flow of H 2 (pressure 0.004 MPa) and H 2 Se, obtained by the interaction of hydrogen with elemental selenium at 770 K. At the outlet of the reaction zone, an excess of H 2 Se was passed through 200-300 g of activated charcoal. The resulting gas was bubbled through a 0.1 M CuSO 4 solution. Activated charcoal was replaced by a fresh portion once flakes of Cu 2 Se (K SP = 60.6 [83]) appeared. Selenidation was carried out with periodic grinding of the resulting product at 970-1040 K for 2.5-12 h and at 1170 K for 6-12 h.
The resulting products were examined by SEM-EDX, and the obtained data are collected in Table 6.

Methods
The powder X-ray diffraction data were collected at room temperature with a ДPOH 7 (Burevestnik, Saint Petersburg, Russia) powder diffractometer (Cu-K α radiation, graphite monochromator). The step size of 2θ was 0.02 • , and the counting time was 13-50 s per step. The crystal structures of EuLnCuSe 3 (Ln = La, Sm, Gd, Tb, Dy, Ho, Y, Tm, Yb, Lu) were refined by the derivative difference minimization (DDM) method [84] in the anisotropic approximation for all atoms. The data for isostructural sulfides Ba 2 MnS 3 [85], Eu 2 CuS 3 [71] and KZrCuS 3 [23] were used as initial structural models. The effects of preferred orientation, anisotropic broadening of peak and sample surface roughness and displacement were taken into account during refinement. Crystal structures were visualized in the program package Diamond 3 [86]. The bond valence method was applied to estimate the oxidation states of atoms [87]. CCDC 2125819-2125828 contain the supplementary crystallographic data. These data can be obtained free of charge via https://www.ccdc.cam.ac.uk/structures (accessed on 23 December 2021) or from the Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; fax: (+44)-1223-336-033; or e-mail: deposit@ccdc.cam.ac.uk.
Scanning electron microscopy (SEM) was performed on a JEOLJSM-6510 LV (JEOL Ltd., Tokyo, Japan) equipped with an energy dispersive spectrometer.
The Fourier-transform infrared (FTIR) absorption spectra in the range of 85-675 cm −1 were recorded on a VERTEX 80v FT-IR spectrometer (Bruker OJSC, Germany). The samples were ground in an agate mortar and then mixed with ultrahigh molecular weight polyethylene (Mitsui Petrochemical Ltd., Tokyo, Japan) and pressed into pellets of about 0.26 mm thickness and of 13 mm diameter. A Globar was used for the light source, and it was equipped with a Bruker multilayer beamsplitter and RT-DTGS FIR as a detector. The attenuated total reflectance infrared (ATR-IR) absorption spectra in the range of 400-4000 cm −1 were recorded on a Cary 630 FTIR spectrometer (Agilent Technologies Inc., Santa Clara, CA, USA) equipped with an ATR attachment and a DTGS detector.
The diffuse reflectance spectra were recorded on a UV-2600 spectrophotometer (Shimadzu OJSC, Tokyo, Japan) equipped with an ISR-2600Plus attachment with the photomultiplier PMT of the R-928 type and InGaAs detectors. BaSO 4 (99.8%) was used as a standard.
Gaseous hydrogen was obtained on a SPEKTR 16 4D hydrogen generator (Spektr Ltd., Petrozavodsk, Russia). The temperature in the electric heating furnaces used for synthesis was controlled using a Termodat-16K6 (InSAT Ltd., Moscow, Russia) temperature controller with a chromel-alumel thermocouple.
The experimental Raman spectra of EuLnCuSe 3 compounds were collected in backscattering geometry, using a triple monochromator Horiba JobinYvon T64000 Raman spectrometer (Horiba Ltd., Tokyo, Japan) operating in subtractive mode. The spectral resolution for the recorded Stokes-side Raman spectra was better than 2.5 cm −1 (this resolution was achieved by using gratings with 1800 grooves mm −1 and 100 µm slits). The single-mode radiation at 532 nm from the Spectra-Physics Excelsior laser was used as an excitation light source, the power on the sample being 1 mW.

DFT Calculations
The density functional theory (DFT) calculations were carried out using the PBE0 and B3LYP hybrid functionals, which consider both local and nonlocal Hartree-Fock exchanges. The calculations were performed using the CRYSTAL17 program designed to simulate periodic structures [91]. For Eu 3+ , the ECP53MWB quasi-relativistic pseudopotential was used with the attached valence basis set ECP53MWB [92]. For Ln 3+ , the ECPnMWB quasirelativistic pseudopotential was used with the attached valence basis set ECPnMWB-I (n = Z − 1; Z is an atomic number) [92]. Thus, the inner shells of the rare-earth ion, including 4f, were replaced by a pseudopotential. To describe the outer shells 5s 2 5p 6 , involved in chemical bonding, valence basis sets were used. Such an approach makes it possible to reconstitute both the lattice structure and dynamics in the compounds that have a lanthanide ion sublattice successfully [10]. For Y 3+ , the ECP28MHB pseudopotential was used with the attached valence basis set ECP28MWB [92]. For copper and selen, we used full-electron basis sets, known in the CRYSTAL program site as Cu_86-4111(41D)G_doll_2000 and Se_976-311d51G_towler_1995, respectively [91]. Gaussian primitives with orbital exponent values less than 0.1 were removed from basis sets since these calculations are periodic. The exponent in the outer p-orbital of the selen basis set was set to 0.1742. The exponent in the outer diffuse s-, pand dorbitals of the yttrium valence basis set was set to 0.1312. The accuracy of calculation of the self-consistent field and the two-electron integrals was set at 10 −9 a.u. and 10 −8 , respectively. Integration over the Brillouin zone was carried out according to the Monkhorst-Pack scheme with a grid of k-points equal to 8 × 8 × 8.
The sequence of the DFT calculations was as follows: the crystal structure was first optimized followed by the calculation of the phonon spectrum at the Γ point or the elastic constants.

Conclusions
In summary, we successfully fabricated novel quaternary selenides of the Eu +2 Ln +3 Cu +1 Se 3 (Ln = La, Sm, Gd, Tb, Dy, Ho, Y, Tm, Yb and Lu) composition through reductive selenidation of the oxide mixtures, produced by thermolysis of the co-crystallized metal nitrates. The applied synthetic procedure allowed a significant decrease in the reaction time and temperature, and an increase in the yield of the final product in comparison to the recently reported procedure [60]. Crystal structures of the resulting compounds were elucidated by powder X-ray diffraction. As a result, the structures were found to belong to orthorhombic space groups Pnma (structural type Ba 2 MnS for EuLaCuSe 3 and structural type Eu 2 CuS 3 for EuLnCuSe 3 , where Ln = Sm, Gd, Tb, Dy, Ho and Y) and Cmcm (structural type KZrCuS 3 for EuLnCuSe 3 , where Ln = Tm, Yb and Lu). These space groups were delimited based on the tolerance factor t , while the formation of three structural types was additionally probed by vibrational spectroscopy and further supported by ab initio calculations with hybrid functionals. It was also established that the revealed structural types for a series of the obtained selenides are the same as for the previously reported sulfide analogues.
According to calculations, the elastic properties change abruptly from one structural type to another. The wavenumbers and types of fundamental vibrations were determined from ab initio calculations, and the experimental IR spectra were interpreted.
Optical properties of the reported compounds EuLnCuSe 3 were studied by diffuse reflectance spectroscopy. The corresponding band gaps are determined by the interconfigurational transitions between sub-bands originating from the occupied f -orbitals and vacant d-orbitals of the Eu 2+ ion.
Field-dependent magnetic properties of the discussed selenides were also revealed by SQUID measurements. The experimental Curie constants and effective magnetic moments at room temperature are very close to calculated ones. Phase transitions to the ferrimagnetic state were established for EuTbCuSe 3 , EuDyCuSe 3 , EuHoCuSe 3 and EuTmCuSe 3 with the transition temperatures varying from 4.5 to 6.2 K. For EuHoCuSe 3 , the effect of negative magnetization, which is characteristic of an N-type ferrimagnet according to the Néel's classification, was found at temperatures below 4.8 K. A comparison of the magnetic properties of quaternary selenides EuLnCuSe 3 with isostructural sulfides showed a high degree of similarity.
Finally, we believe that a series of the reported selenides together with the revealed properties will be of value for the fabrication of novel materials with properties of interest. Data Availability Statement: All the data supporting the conclusions is included within the manuscript and is available on request from the corresponding authors.