{"id":96822,"date":"2025-10-01T12:08:59","date_gmt":"2025-10-01T12:08:59","guid":{"rendered":"https:\/\/www.europesays.com\/ie\/96822\/"},"modified":"2025-10-01T12:08:59","modified_gmt":"2025-10-01T12:08:59","slug":"the-effect-of-non-stoichiometry-in-sb2se3-films-on-their-phonon-and-electronic-properties-in-the-infrared-range","status":"publish","type":"post","link":"https:\/\/www.europesays.com\/ie\/96822\/","title":{"rendered":"The effect of non-stoichiometry in Sb2Se3 films on their phonon and electronic properties in the infrared range"},"content":{"rendered":"<p>First principle calculations of electron and phonon properties of SbxSey films<\/p>\n<p>The calculated dielectric permittivities \u03b5\u221e (electronic contribution) and \u03b5ion (ionic contribution) of Sb2Se3 crystal (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>) are in acceptable agreement with the experimental ones presented in Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Petzelt, J. &amp; Grigas, J. Far infrared dielectric dispersion in Sb2S3, Bi2S3 and Sb2Se3 single crystals. Ferroelectrics 5, 59&#x2013;68 (1973).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR21\" id=\"ref-link-section-d20448370e2403\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>. To our knowledge, Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Petzelt, J. &amp; Grigas, J. Far infrared dielectric dispersion in Sb2S3, Bi2S3 and Sb2Se3 single crystals. Ferroelectrics 5, 59&#x2013;68 (1973).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR21\" id=\"ref-link-section-d20448370e2407\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>. provides a detailed experimental study of the reflection spectra and corresponding dielectric functions of Sb2Se3 single crystals in the far-infrared range, as well as a theoretical study of the phonon modes for the crystal at the orthorhombic space group of symmetry Pnma (no. 62).<\/p>\n<p>In the present study, instead, the experimental reflection spectra of Sb2Se3 polycrystalline films deposited on soda-lime-glass substrates are measured in the far-infrared range. The stoichiometry and crystal structure of these polycrystalline films may differ from the corresponding features in the above-mentioned Sb2Se3 single crystals. Furthermore, first-principles calculations based on density functional theory of the Sb\u2082Se\u2083 single crystal were performed, taking into account the van der Waals interaction, and the corresponding dielectric function was obtained. The imaginary part of this dielectric function is used for comparison with the experimental reflection spectra of the studied Sb\u2082Se\u2083 films. The present study aims to provide more information on the electron and phonon properties of both stoichiometric and non-stoichiometric Sb2Se3 crystals, which may be utilized for the practical applications of Sb2Se3 films in solar cell technology.<\/p>\n<p>According to the analysis of the accuracy of the experimental dielectric function curves \u03b5(\u03bd) performed in Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Petzelt, J. &amp; Grigas, J. Far infrared dielectric dispersion in Sb2S3, Bi2S3 and Sb2Se3 single crystals. Ferroelectrics 5, 59&#x2013;68 (1973).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR21\" id=\"ref-link-section-d20448370e2441\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>., the accuracy of determining the electronic dielectric permittivity \u03b5\u221e was estimated to be about 20%, while the corresponding error for the ionic dielectric permittivity, \u0394\u03b5ion\/\u03b5ion, can be larger. The relative differences between the averaged theoretical (row 5 of Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>) and experimental (row 6 of Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>) dielectric permittivities are presented in row 7 of Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>. The obtained relative differences 2\u00d7|\u03b5t &#8211; \u03b5e|\/(\u03b5t\u2009+\u2009\u03b5e), appear acceptable when considering that the theoretical and experimental accuracies of determining the dielectric permittivity are both approximately 20%. The characteristic feature of Sb2Se3 crystal is the large anisotropy of the calculated dielectric permittivity tensor \u03b5: the contributions to the averaged dielectric permittivity from the y- and z-directions of the crystal are much larger than that from the x- one (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). The spontaneous electric polarization in Sb2Se3 crystal and SbxSey deposited films is not expected because the space group of symmetry Pnma is nonpolar. The absence of polar symmetry in SbxSey films may be partly caused by the relatively small difference of antimony and selenium electronegativities, \u03c7(Se) &#8211; \u03c7(Sb)\u2009=\u20090.5 (Pauling scale), which corresponds to the relatively small ionic character (IC)<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 40\" title=\"Callister, W. D. &amp; Rethwisch, D. G. Jr. Materials Science and Engineering: An Introduction 43 (Wiley, 2014).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR40\" id=\"ref-link-section-d20448370e2507\" rel=\"nofollow noopener\" target=\"_blank\">40<\/a> of Sb Se bonding, IC = %. Small ionic character of Sb Se bonds may cause small Sb Se dipole moments in SbxSey films, which can decrease the probability of polarized parts arising in these films during their deposition on a substrate.<\/p>\n<p><b id=\"Tab2\" data-test=\"table-caption\">Table 2 Calculated cartesian components of the macroscopic static dielectric permittivities, \u03b5\u221e (electronic contribution) and \u03b5ion (ionic contribution), of Sb2Se3 crystal at the space group of symmetry Pnma (no. 62) calculated for GGA-PE, GGA-PS, and GGA-MK\u2009+\u2009vdW exchange-and-correlation functionals.<\/b><\/p>\n<p>It is known<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 13\" title=\"Mavlonov, A. et al. Qiao, L. A review of Sb2Se3 photovoltaic absorber materials and thin-film solar cells. Sol Energy. 201, 227&#x2013;246 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR13\" id=\"ref-link-section-d20448370e2842\" rel=\"nofollow noopener\" target=\"_blank\">13<\/a> that the interatomic bonding in Sb2Se3 crystal is strongly anisotropic because the covalently bonded (Sb4Se6)n ribbon-like structure elements are realized along one of the crystallographic directions ([001]) and weaker forces take place in the two other orthogonal directions. Very high anisotropy of the dielectric properties in Sb2Se3 is observed for the ionic part of the dielectric permittivity, \u03b5ion_x\u2009ion_y, \u03b5ion_x\u2009ion_z (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). Similar but not so high anisotropy of the dielectric properties in Sb2Se3 takes place for the electronic part of the dielectric permittivity, \u03b5\u221ex\u2009\u221ey, \u03b5\u221ex\u2009\u221ez (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). The observed anisotropy of the dielectric permittivity tensor corresponds to the layer-like crystal structure of Sb2Se3, for which the corresponding quasi-layers are perpendicular to the [100] direction (a-direction) (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). Thus, the comparatively weak van-der-Waals interaction is realized between the quasi-layers with boundaries indicated in the b-cut view of Sb2Se3 crystal structure by the thick horizontal lines (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). The relatively weak van-der-Waals interaction in [100] direction of the crystal is a reason for the much smaller ionic and smaller electronic parts of dielectric permittivity in this direction in comparison to [010] and [001] ones (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). No layer-like crystal structure is found in the [010] direction (b-direction).<\/p>\n<p><b id=\"Fig5\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 5<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/5\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig5\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig5_HTML.png\" alt=\"figure 5\" loading=\"lazy\" width=\"685\" height=\"446\"\/><\/a><\/p>\n<p>The b-cut view of Sb2Se3 crystal structure at Pnma (no. 62) space group of symmetry. Thick horizontal lines indicate planes of boundaries between structural layers perpendicular to the a-axis of the crystal\u2019s unit cell.<\/p>\n<p>Dispersion of the phonon energy bands of Sb2Se3 crystal is observed in the frequency range of 0\u2013220\u00a0cm\u22121, where the acoustic modes are realized in the range of 0\u201350\u00a0cm\u22121 and the optical modes are observed in the range of 50\u2013220\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig6\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>a, b). The phonon dispersions of Sb2Se3 calculated for PBE and PBEsol exchange-and-correlation functional are close to one another (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig6\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>a, b). Negative frequencies were observed in the dispersion dependence \u03bd(Q), calculated without taking into account the vdW interaction and this happened in the Q-segments \u0393 X and Z \u2013 U, where the Q[100] component was changing. This result is caused by the relatively weak interatomic interaction in the [100] direction of the crystal (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). Another manifestation of the relatively weak interatomic interaction in the [100] direction of the crystal is the flat dispersion dependences \u03bd(Q) at the segments S \u2013 Y, Z \u2013 U, and R \u2013 T, where the changes of the type, [0QyQz] \u2194 [\\(\\frac{1}{2}\\)QyQz], of the phonon Q-points in the Brillouin zone take place.<\/p>\n<p><b id=\"Fig6\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 6<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/6\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig6\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig6_HTML.png\" alt=\"figure 6\" loading=\"lazy\" width=\"685\" height=\"259\"\/><\/a><\/p>\n<p>Dispersions of phonon energy bands of Sb2Se3 calculated with using van-der-Waals interaction of DFT-D3 method of Grimme<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 41\" title=\"Grimme, S., Antony, J., Ehrlich, S. &amp; Krieg, S.&#xA0;A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu.&#xA0;J. Chem. Phys. 132, 154104 (2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR41\" id=\"ref-link-section-d20448370e3035\" rel=\"nofollow noopener\" target=\"_blank\">41<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 42\" title=\"Grimme, S., Ehrlich, S. &amp; Goerigk, L.&#xA0;Effect of the damping function in dispersion corrected density functional theory.&#xA0;J. Comput. Chem. 32, 1456 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR42\" id=\"ref-link-section-d20448370e3038\" rel=\"nofollow noopener\" target=\"_blank\">42<\/a> at (a) PBE and (b) PBEsol exchange-and-correlation functionals. The correspondence of the Brillouin zone symbols and Q-points is following: \u0393[000], X[\\(\\frac{1}{2}\\)00], S[\\(\\frac{1}{2}\\)\\(\\frac{1}{2}\\)0], Y[0\\(\\frac{1}{2}\\)0], Z[00\\(\\frac{1}{2}\\)], U[\\(\\frac{1}{2}\\)0\\(\\frac{1}{2}\\)], R[\\(\\frac{1}{2}\\)\\(\\frac{1}{2}\\)\\(\\frac{1}{2}\\)], T[0\\(\\frac{1}{2}\\)\\(\\frac{1}{2}\\)].<\/p>\n<p><b id=\"Fig7\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 7<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/7\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig7\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig7_HTML.png\" alt=\"figure 7\" loading=\"lazy\" width=\"685\" height=\"259\"\/><\/a><\/p>\n<p>(a) Vibration density of states vDOS and (b) vibration partial density of states vPDOS of antimony and selenium atoms in Sb2Se3 crystal calculated at PBE exchange-and-correlation functional.<\/p>\n<p>The calculated vibration total and partial densities of states (vDOS and vPDOS) of Sb2Se3 crystal (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig7\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>) are characterized by numerous maxima, which correspond to the numerous groups of phonon energy bands (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig6\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>). Clear hybridization of the antimony and selenium vibration states is observed, which is manifested by the equality of the frequency positions of vPDOS(Sb) and vPDOS(Se) maxima (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig7\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>b). The characteristic feature of these vPDOS is a dominance of the antimony density of states in the frequency range 0\u201390\u00a0cm\u22121 and a dominance of the selenium density of states in the range 90\u2013190\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig7\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>b). The latter feature is expected because of the larger atomic mass of antimony atom (M(Sb)\u2009=\u2009121.76) in comparison to the mass of selenium one (M(Sb)\u2009=\u200978.971).<\/p>\n<p>The results of our DFT-based calculations of lattice dynamics of Sb2Se3 crystal have shown that the corresponding infrared phonon excitations take place in the frequency range of 50\u2013220\u00a0cm\u22121. According to these results, the spectral range above 220\u00a0cm\u22121 should be the range of crystal transparency, if some other mechanisms of radiation absorption are absent. However, from the SEM and EDX results of the deposited films, we also expect the metallic phase of the redundant antimony, which can be the source of free electrons. These free electrons may cause the inverse proportional dependence of the reflection coefficient R on the frequency \u03bd<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Yu, P. Y. &amp; Cardona, M. Fundamentals of Semiconductors 310 (Springer, 2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR43\" id=\"ref-link-section-d20448370e3260\" rel=\"nofollow noopener\" target=\"_blank\">43<\/a>.<\/p>\n<p>Reflection spectra in the far-infrared spectral range SbxSey films<\/p>\n<p>Analysis of SbxSey reflection spectra R(\u03bd) shows the presence of an interference-like structure in the frequency range of \u03bd\u2009&gt;\u20091200\u00a0cm\u22121 for several samples (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>a). The distances between neighboring interference extrema for samples no. 3, 4, and 7 correspond to film thicknesses of about 1\u00a0\u03bcm, which are close to the values determined by the other method (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>). The absence of a clear interference-like structure for other films in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>a can be caused by the relatively strong radiation scattering on their rough surfaces. This suggestion is in agreement with the observed quasi-monotonous decrease of the reflection coefficient R with the increase of the frequency \u03bd, caused by the surface roughness, which is expected for the mean surface roughness of about 0.5\u00a0\u03bcm (see the cross-sections of films in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>).<\/p>\n<p><b id=\"Fig8\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 8<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/8\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig8\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig8_HTML.png\" alt=\"figure 8\" loading=\"lazy\" width=\"685\" height=\"262\"\/><\/a><\/p>\n<p>(a) Reflection spectra R(\u03bd) of eight SbxSey films on soda-lime glass substrates (SLG) and SLG measured at room temperature. The legend values correspond to the index of compound relation x\/y determined from EDX measurements and sample number (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>b, d). The values of R correspond to the absolute reflection coefficients. The experimental calibration of the R scale was performed using the gold mirror and its known optical constants. (b) Reflection spectra R(\u03bd) of eight SbxSey films on soda-lime glass substrates (SLG) and SLG measured at room temperature and theoretical spectra of imaginary part of dielectric permittivity \u03b52(\u03bd) (narrow black, red and green peaks). Here, for \u03b52(\u03bd) spectra, black, red and green colors correspond to the light polarizations along x-, y-, and z-axis of Sb2Se3 crystal at the space group of symmetry Pnma (no. 62). The legend values correspond to the index of compound relation x\/y determined from EDX measurements and sample number (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>b, d). The values of R correspond to the absolute reflection coefficients. The experimental calibration of the R scale was performed using the gold mirror and its known optical constants.<\/p>\n<p>The reflection spectra R(\u03bd) in the range \u03bd\u2009\u22121 are presented in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b for SbxSey films studied. The largest reflection coefficients are observed for the samples 1, 2 and 5 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). This is caused by the increased amount of metal phase of antimony in the deposited films and, consequently, by the increased reflection coefficients due to Drude-like behavior in the far-infrared range<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Yu, P. Y. &amp; Cardona, M. Fundamentals of Semiconductors 310 (Springer, 2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR43\" id=\"ref-link-section-d20448370e3409\" rel=\"nofollow noopener\" target=\"_blank\">43<\/a>. The mean value\/level of the reflection coefficient of SbxSey films in the range of 25 230\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) is proportional to the concentration of residual antimony atoms (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>b), being in the metal-like state, compared to their bonded states in the Sb2Se3 polycrystals.<\/p>\n<p>The reflection spectra of SbxSey films in the far-IR range are similar because of the almost identical positions of R(\u03bd) peaks (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). The highest similarity is observed for the spectra of samples 4\u20138, which may confirm the formation of the same crystalline phase Sb2Se3 in these samples. These positions are in good quantitative agreement with the positions of maxima of the imaginary part of the dielectric function \u03b52(\u03bd) obtained from our first-principle calculations of Sb2Se3 crystal within the framework of DFT (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). The structure of the reflection spectra R(\u03bd) measured for SbxSey films (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) is rather similar to that obtained experimentally for Sb2Se3 single crystals<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Petzelt, J. &amp; Grigas, J. Far infrared dielectric dispersion in Sb2S3, Bi2S3 and Sb2Se3 single crystals. Ferroelectrics 5, 59&#x2013;68 (1973).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR21\" id=\"ref-link-section-d20448370e3472\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>. The latter was also observed in a similar spectral range, 25\u2013210\u00a0cm\u22121, and was characterized by a similar number of reflection maxima. One of the clear differences between the experimental reflection spectra of the Sb2Se3 films studied (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) and the reflection spectrum of an Sb2Se3 single crystal is the reflection maximum in Sb2Se3 films in the range 221\u2013225\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b), which is absent in this frequency range in Sb2Se3 single crystal<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Petzelt, J. &amp; Grigas, J. Far infrared dielectric dispersion in Sb2S3, Bi2S3 and Sb2Se3 single crystals. Ferroelectrics 5, 59&#x2013;68 (1973).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR21\" id=\"ref-link-section-d20448370e3505\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a>. This should be caused by the vibration mode associated with the special states of selenium ions in Sb2Se3 films, which are absent in the corresponding single crystal. This should be associated with the selenium ions, first of all, because the vibration states of selenium ions are characterized by a larger vPDOS compared to the vPDOS of antimony ones in the high-frequency part of the Sb2Se3 phonon spectrum (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). Such special vibration states of the selenium ions may be the states of Se \u2013 Sb bonds created on the surface of Sb2Se3 nanocrystals in SbxSey films. The corresponding frequency shifts of the Raman-active and the IR-active modes for sized and layered two-dimensional structures were observed in different nanomaterials<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Zi, J. et al. Raman shifts in Si nanocrystals. Appl. Phys. Lett. 69 (2), 200&#x2013;202 (1996).\" href=\"#ref-CR44\" id=\"ref-link-section-d20448370e3529\">44<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Yang, X. X. et al. Raman spectroscopic determination of the length, strength, compressi4bility, Debye temperature, elasticity, and force constant of the C-C bond in graphene. Nanoscale 4 (2), 502&#x2013;510 (2012).\" href=\"#ref-CR45\" id=\"ref-link-section-d20448370e3529_1\">45<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Fujii, M., Kanzawa, Y., Hayashi, S. &amp; Yamamoto, K. Raman scattering from acoustic phonons confined in Si nanocrystals. Phys Rev. B 54(12), R8373&#x2013;R8376 (1996).\" href=\"#ref-CR46\" id=\"ref-link-section-d20448370e3529_2\">46<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Cheng, W. &amp; Ren, S. F. Calculations on the size effects of Raman intensities of silicon quantum Dots. Phys. Rev. B. 65 (20), 205305 (2002).\" href=\"#ref-CR47\" id=\"ref-link-section-d20448370e3529_3\">47<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Lee, C. et al. Anomalous lattice vibrations of single and few-layer MoS2. ACS Nano. 4 (5), 2695&#x2013;2700 (2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR48\" id=\"ref-link-section-d20448370e3532\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>. The reflection maximum in Sb2Se3 films in the range 221\u2013225\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) may originate, in principle, from the Se \u2013 Se vibrations of the defect structure of SbxSey films studied. This possibility is mentioned in Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Razykov, T. M. et al. Characteristics of thin Sb2Se3 films obtained by the chemical molecular beam deposition method for thin-film solar cells. Thin Solid Films. 774, 139844 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR49\" id=\"ref-link-section-d20448370e3550\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>. (see the selected Raman frequencies in Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>). On the basis of the experimental reflection spectra obtained in the far-IR range for eight SbxSey films (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b), we tried to study the correlation between the relative value of the reflection maximum, Rmax &#8211; Rmin, in the range 221\u2013225\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) and the relation c(Se)\/c(Sb) of the concentrations of selenium and antimony atoms (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>c). We have found a weak positive correlation of these values (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig9\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>a), which supports the statement of the prevalent selenium origin of the reflection maximum in the frequency range 221\u2013225\u00a0cm\u22121.<\/p>\n<p><b id=\"Fig9\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 9<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/9\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig9\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig9_HTML.png\" alt=\"figure 9\" loading=\"lazy\" width=\"685\" height=\"260\"\/><\/a><\/p>\n<p>(a) Dependence of the relative value of the reflection maximum, Rmax &#8211; Rmin, in the range 221\u2013225\u00a0cm\u22121 (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) and the relation of the concentrations of selenium and antimony atoms c(Se)\/c(Sb), (see Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>c). The red direct line is the linear fit of this dependence. (b) Raman spectrum IRaman(\u03bd) of Sb2Se3 crystal at the space group symmetry of Pnma (No. 62) calculated using CASTEP code at PBEsol exchange-and-correlation functional.<\/p>\n<p><b id=\"Tab3\" data-test=\"table-caption\">Table 3 Peak positions of Raman modes (in cm\u22121) of SbxSey films deposited at different source temperatures in the range 830\u20131000\u00a0\u00b0C and Raman modes of different crystalline substances containing Sb and se atoms<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Razykov, T. M. et al. Characteristics of thin Sb2Se3 films obtained by the chemical molecular beam deposition method for thin-film solar cells. Thin Solid Films. 774, 139844 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR49\" id=\"ref-link-section-d20448370e3668\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>.<\/b><\/p>\n<p>We have calculated the Raman spectrum of the Sb2Se3 crystal at the space group symmetry of Pnma (No. 62) (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig9\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>b) using CASTEP code<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 50\" title=\"Clark, S. J. et al. First principles methods using CASTEP, zeitschrift fur kristallographie -. Crystalline Mater. 220, 567&#x2013;570 (2005).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR50\" id=\"ref-link-section-d20448370e4276\" rel=\"nofollow noopener\" target=\"_blank\">50<\/a> and have found that it is very similar to the experimental Raman spectra obtained for Sb2Se3 films deposited by the same method (CMBD) on the same substrates (soda-lime glass) by several authors of the present study Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Razykov, T. M. et al. Characteristics of thin Sb2Se3 films obtained by the chemical molecular beam deposition method for thin-film solar cells. Thin Solid Films. 774, 139844 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR49\" id=\"ref-link-section-d20448370e4285\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>. Besides, the calculated frequency positions of the Raman spectrum of Sb2Se3 crystal (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig9\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>b) agree satisfactorily with the corresponding experimental Raman spectrum<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 51\" title=\"Vidal-Fuentes, P. et al. Izquierdo roca. V. Multiwavelength excitation Raman scattering study of Sb2Se3 compound: fundamental vibrational properties and secondary phases detection. 2D Mater. 6, 045054 (2019).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR51\" id=\"ref-link-section-d20448370e4297\" rel=\"nofollow noopener\" target=\"_blank\">51<\/a> and with Raman spectra of Sb2Se3 films<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 52\" title=\"Shongalova, A. et al. On the identification of Sb2Se3 using Raman scattering. MRS Commun.&#xA0;&#010;                  https:\/\/doi.org\/10.1557\/mrc.2018.94&#010;                  &#010;                \" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR52\" id=\"ref-link-section-d20448370e4305\" rel=\"nofollow noopener\" target=\"_blank\">52<\/a>.<\/p>\n<p>Discussion<\/p>\n<p>The increased reflection coefficients of SbxSey films for several x\/y relations in the wide spectral range 30\u20135000\u00a0cm\u22121 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>) may be caused by the over-stoichiometric antimony content, leading to the metal-like character of the radiation reflection. To verify this supposition we have performed DFT calculations of the electron energy structure of the contents Sb25Se36 and Sb26Se36, which contain one (Sb25Se36) and two (Sb26Se36) residual antimony atoms in the 1\u2009\u00d7\u20093\u2009\u00d7\u20091 supercell of the elementary unit cell of Sb8Se12. The optimized positions of these residual antimony atoms Sb1 and Sb2 are shown in Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig10\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>.<\/p>\n<p><b id=\"Fig10\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 10<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/10\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig10\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig10_HTML.png\" alt=\"figure 10\" loading=\"lazy\" width=\"685\" height=\"385\"\/><\/a><\/p>\n<p>Positions of the residual antimony atoms Sb1 and Sb2 in the over-stoichiometric contents Sb25Se36 (a) and Sb26Se36 (b) corresponding to the supercell 1\u2009\u00d7\u20093\u2009\u00d7\u20091 of the elementary unit cell 1\u2009\u00d7\u20091\u2009\u00d7\u20091 (Sb8Se12). The distance d\u2009=\u20098.32 \u00c5 between two residual antimony atoms is indicated in the part (b) of the figure.<\/p>\n<p>For the calculations of the density of states and bandgap Eg of the crystals Sb24Se36, Sb25Se36, and Sb26Se36, we have utilized the HSEsol hybrid exchange-and-correlation functional<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 53\" title=\"Schimka, L., Harl, J. &amp; Kresse, G. Improved hybrid functional for solids: the HSEsol functional. J. Chem. Phys. 134, 024116 (2011).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR53\" id=\"ref-link-section-d20448370e4413\" rel=\"nofollow noopener\" target=\"_blank\">53<\/a>, which is often used for the calculations of the energy bandgap Eg of crystals with high reliability. The obtained bandgap of the stoichiometric crystal Sb24Se36 is found to be Eg = 1.28\u00a0eV. In the calculated density of states of the content Sb25Se36, one separate maximum of the valence states is at the energy of 4.8\u00a0eV, which is below the Fermi level of the crystal EF = 5.414\u00a0eV (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig11\" rel=\"nofollow noopener\" target=\"_blank\">11<\/a>a). The maximum of the partial density of states of Sb25Se36 at 4.8\u00a0eV (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig11\" rel=\"nofollow noopener\" target=\"_blank\">11<\/a>b) is formed mainly by the electronic states of the residual antimony ion. In the over-stoichiometric content Sb25Se36, the bandgap is smaller, Eg = 0.923\u00a0eV, and numerous partly-occupied electronic states are in the upper part of the forbidden energy gap (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig12\" rel=\"nofollow noopener\" target=\"_blank\">12<\/a>). The electron states at the Fermi energy, E\u2009\u2265\u2009EF\u200b, and the thermally stimulated transitions of electrons from the upper donor energy states of the residual antimony atoms to the conduction band of Sb\u2082\u2085Se\u2083\u2086 can form a free electron plasma. This plasma increases the optical reflection coefficient of the Sb\u2082\u2085Se\u2083\u2086 layer in the range of 25\u2013250\u00a0cm\u207b\u00b9 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b) according to the Drude relation<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Yu, P. Y. &amp; Cardona, M. Fundamentals of Semiconductors 310 (Springer, 2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR43\" id=\"ref-link-section-d20448370e4472\" rel=\"nofollow noopener\" target=\"_blank\">43<\/a>. The electron states just below the Fermi level EF\u200b (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig12\" rel=\"nofollow noopener\" target=\"_blank\">12<\/a>) may also be excited by photons to the lowest states of the conduction band (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig11\" rel=\"nofollow noopener\" target=\"_blank\">11<\/a>), thus forming the infrared absorption spectrum in the photon energy range of 1200\u20135000\u00a0cm\u207b\u00b9 (0.15\u20130.62\u00a0eV).<\/p>\n<p><b id=\"Fig11\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 11<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/11\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig11\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig11_HTML.png\" alt=\"figure 11\" loading=\"lazy\" width=\"685\" height=\"274\"\/><\/a><\/p>\n<p>Calculated density of states (a) and partial density of states (b) of antimony selenide Sb25Se36 with one additional antimony atom in 1\u2009\u00d7\u20093\u2009\u00d7\u20091 supercell (Sb24Se36) of the elementary unit cell (Sb8Se12) of Sb2Se3.<\/p>\n<p><b id=\"Fig12\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 12<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/12\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig12\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig12_HTML.png\" alt=\"figure 12\" loading=\"lazy\" width=\"685\" height=\"479\"\/><\/a><\/p>\n<p>Energy of electron states E, in the range of the bandgap Eg, versus their electron occupancy calculated for non-stoichiometric Sb25Se36 (Sb\u2009+\u2009Sb24Se36) with one additional antimony atom in 1\u2009\u00d7\u20093\u2009\u00d7\u20091 supercell (Sb24Se36) of the elementary unit cell (Sb8Se12) of Sb2Se3.<\/p>\n<p>The optical bandgap Eg of SbxSey films of thickness d was determined from the measurements of their transmittance T as a function of the photon energy h\u03bd (Fig. SD-9 of Supplementary designation). The bandgap Eg was obtained from the slope of the linear part of the dependence of the value (-h\u03bd\u22c5ln(T)\/d)2 on the photon energy h\u03bd, corresponding to the direct optical transitions. It was found that the dependence of the bandgap Eg of SbxSey films on their content index x\/y (x\/y\u2009=\u2009NSb\/NSe) is close to a linear one with a negative slope coefficient dEg\/d(x\/y) and is similar to that obtained in the previous study<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"Razykov, T. M. et al. Structural and optical properties of SbxSey thin films obtained by chemical molecular beam deposition method from Sb and se precursors. Sol Energy. 254, 67&#x2013;72 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR54\" id=\"ref-link-section-d20448370e4644\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a> (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig13\" rel=\"nofollow noopener\" target=\"_blank\">13<\/a>a). However, the absolute value of the experimental slope coefficient |dEg\/d(x\/y)| = 0.292 (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig13\" rel=\"nofollow noopener\" target=\"_blank\">13<\/a>a) is substantially smaller than the theoretical one |dEg\/d(x\/y)| = 10.07 (13b). This is likely caused by the fact that the majority of the residual antimony atoms are placed out of the polycrystals of SbxSey films in the form of a metallic phase in SbxSey films as already seen by XRD analysis. The latter circumstance also explains the experimentally observed increase of the reflection coefficient of SbxSey films with an increase of the NSb\/NSe relation (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>).<\/p>\n<p><b id=\"Fig13\" class=\"c-article-section__figure-caption\" data-test=\"figure-caption-text\">Fig. 13<\/b><a class=\"c-article-section__figure-link\" data-test=\"img-link\" data-track=\"click\" data-track-label=\"image\" data-track-action=\"view figure\" href=\"https:\/\/www.nature.com\/articles\/s41598-025-17551-8\/figures\/13\" rel=\"nofollow noopener\" target=\"_blank\"><img decoding=\"async\" aria-describedby=\"Fig13\" src=\"https:\/\/www.europesays.com\/ie\/wp-content\/uploads\/2025\/10\/41598_2025_17551_Fig13_HTML.png\" alt=\"figure 13\" loading=\"lazy\" width=\"685\" height=\"278\"\/><\/a><\/p>\n<p>Dependences (a) of the optical bandgap Eg of SbxSey films determined by the transmittance measurements T(h\u03bd) (red points relate to the corresponding experimental study in Ref<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"Razykov, T. M. et al. Structural and optical properties of SbxSey thin films obtained by chemical molecular beam deposition method from Sb and se precursors. Sol Energy. 254, 67&#x2013;72 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR54\" id=\"ref-link-section-d20448370e4721\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a>. and (b) of the theoretical optical bandgap Eg of Sb24Se36, Sb25Se36, and Sb26Se36 crystals on their atomic content NSb\/NSe.<\/p>\n<p>To confirm the supposition that the majority of the non-stoichiometric antimony atoms can be aggregated outside the polycrystals Sb2Se3 in SbxSey films we have performed comparative calculations of the total energy of (a) one and two residual antimony atoms inside Sb24Se36 supercell and (b) one residual antimony atom adsorbed on the surface of Sb48Se72 supercell. In the latter case, a twice larger supercell dimension perpendicular to the surfaces of the Sb48Se72 slab was chosen to minimize the atomic interaction between atoms of two parallel slabs\u2019 surfaces. For the residual antimony atom localized inside the Sb24Se36 supercell, the binding energy Eb(Sb) ranges from 3.10 to 3.51\u00a0eV (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>). For the non-stoichiometric antimony atom localized on the slab\u2019s surfaces [100], [010], and [001] originated from the Sb48Se72 supercell, the binding energy Eb(Sb) is found to be in the range of 3.35 4.19\u00a0eV (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). These results indicate that the probability of localization of the non-stoichiometric antimony atoms on the [100], [010], and [001] facets of the polycrystals Sb2Se3 in SbxSey films is greater than that inside these polycrystals.<\/p>\n<p><b id=\"Tab4\" data-test=\"table-caption\">Table 4 Total energy E for supercell contents Sb24Se36, Sb25Se36, and Sb26Se36 and binding energy Eb(Sb) for the residual antimony atom in Sb25Se36 and Sb26Se36 supercells.<\/b><b id=\"Tab5\" data-test=\"table-caption\">Table 5 Total energy E of Sb24Se36 slabs [100], [010], and [001] and one antimony atom on the surfaces of these slabs and binding energy of the residual antimony atom Eb(Sb) on the surfaces of these slabs.<\/b><\/p>\n<p>On the other hand, we have found that the calculated cohesive (or binding) energies of the metallic antimony (Sb) and antimony selenide (Sb2Se3) crystals are respectively equal to Ec(Sb)\u2009=\u20093.12\u00a0eV and Ec(Sb2Se3)\u2009=\u20093.34\u00a0eV. This, in turn, indicates that the probability of crystal growth of Sb2Se3 is higher than that of the metallic antimony. These results indicate that during SbxSey film deposition, the growth of Sb2Se3 crystals is preferable to the independent growth of metallic antimony crystals. However, the growth of metallic antimony crystals becomes more probable on the surfaces of Sb2Se3 crystals due to the high bonding energy of the antimony atoms Eb(Sb) on the facets of Sb2Se3 crystals (Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>). Moreover, in the late stages of SbxSey film deposition, the over-stoichiometric antimony atoms of Sb2Se3 crystals can move due to diffusion from the inner parts of the Sb2Se3 crystal to its surface, because of the larger binding energy of the antimony atoms on the crystal facets (3.35\u20134.19\u00a0eV, Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab5\" rel=\"nofollow noopener\" target=\"_blank\">5<\/a>) in comparison to its body (3.10\u20133.51\u00a0eV, Table\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"table anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Tab4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>). These factors may promote the possibility of depositing the over-stoichiometric antimony on the facets of Sb2Se3 polycrystals in the form of a thin metallic antimony layer. In turn, such a metallic antimony layer on SbxSey films studied may be a reason for the increasing reflection coefficient in the far-infrared range with increasing their antimony content due to the free electron plasma (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). The proposed explanation agrees with the XRD measurements of the films studied (see Chap.\u00a02), in which two main substances, the antimony selenide Sb2Se3 and the residual metallic Sb, were detected.<\/p>\n<p>Summarizing the results obtained for Sb\u2082Se\u2083-based films, one can state that the experimental and theoretical studies of the optical properties of the films in the wide infrared range of 30\u20138000\u00a0cm\u207b\u00b9 can be useful for checking the necessary properties important for the photoelectric applications of these and similar semiconductor films.<\/p>\n<p>In the spectral range of 500\u20138000\u00a0cm\u207b\u00b9 (0.06\u20130.99\u00a0eV), one can obtain information mainly on the possible local and distributed electron states caused by doping and structural imperfections in the films deposited by different methods. The local electron states can be detected in the form of spectral bands, which are frequently observed at relatively low temperatures. In turn, the distributed electron states are revealed as quasi-free electrons or holes, which can be fitted by the known Drude-Lorentz relation<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 43\" title=\"Yu, P. Y. &amp; Cardona, M. Fundamentals of Semiconductors 310 (Springer, 2010).\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#ref-CR43\" id=\"ref-link-section-d20448370e5289\" rel=\"nofollow noopener\" target=\"_blank\">43<\/a>. In the studied Sb\u2082Se\u2083-based films, both types of electron states occur. The local electron states in these films can cause the dependence of their bandgap Eg\u200b on the antimony content (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig13\" rel=\"nofollow noopener\" target=\"_blank\">13<\/a>). In turn, the quasi-free electrons or holes in the films are revealed in the increased level of the reflection coefficient for the frequencies, which are proportional to the concentration of these carriers (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b).<\/p>\n<p>In the spectral range of 30\u2013240\u00a0cm\u207b\u00b9 (3.72\u201330 meV), phonon excitations in Sb\u2082Se\u2083-based films are revealed in the form of well-defined spectral peaks (Fig.\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41598-025-17551-8#Fig8\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>b). The phonon spectral structure is closely related to the spatial symmetry and interatomic distances in the Sb\u2082Se\u2083 polycrystals that form the deposited films. Thus, these characteristics, when obtained for the deposited Sb\u2082Se\u2083-based films, may also be useful for the optimization of the corresponding deposition conditions.<\/p>\n","protected":false},"excerpt":{"rendered":"First principle calculations of electron and phonon properties of SbxSey films The calculated dielectric permittivities \u03b5\u221e (electronic contribution)&hellip;\n","protected":false},"author":2,"featured_media":96823,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[271],"tags":[22299,62101,914,18,62107,35155,1306,1099,19,17,62103,1100,8408,7072,62106,62102,452,1098,62105,62104,133],"class_list":{"0":"post-96822","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-physics","8":"tag-chemical-physics","9":"tag-computational-methods","10":"tag-condensed-matter-physics","11":"tag-eire","12":"tag-electron-properties","13":"tag-electronic-structure","14":"tag-electronics","15":"tag-humanities-and-social-sciences","16":"tag-ie","17":"tag-ireland","18":"tag-materials-for-energy-and-catalysis","19":"tag-multidisciplinary","20":"tag-nanoscale-materials","21":"tag-optical-physics","22":"tag-phonon-spectra","23":"tag-photonics-and-device-physics","24":"tag-physics","25":"tag-quantum-physics","26":"tag-sb2se3-crystal","27":"tag-sbxsey-films","28":"tag-science"},"share_on_mastodon":{"url":"","error":""},"_links":{"self":[{"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/posts\/96822","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/comments?post=96822"}],"version-history":[{"count":0,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/posts\/96822\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/media\/96823"}],"wp:attachment":[{"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/media?parent=96822"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/categories?post=96822"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.europesays.com\/ie\/wp-json\/wp\/v2\/tags?post=96822"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}