Rasmus S.
Nielsen
*a,
Ángel Labordet
Álvarez
abc,
Axel G.
Medaille
de,
Ivan
Caño
de,
Alejandro
Navarro-Güell
de,
Cibrán L.
Álvarez
ef,
Claudio
Cazorla
ef,
David R.
Ferrer
de,
Zacharie J.
Li-Kao
de,
Edgardo
Saucedo
de and
Mirjana
Dimitrievska
*a
aNanomaterials Spectroscopy and Imaging, Transport at Nanoscale Interfaces Laboratory, Swiss Federal Laboratories for Material Science and Technology (EMPA), Ueberlandstrasse 129, 8600 Dübendorf, Switzerland. E-mail: rasmus.nielsen@empa.ch; mirjana.dimitrievska@empa.ch
bDepartment of Physics, University of Basel, 4056 Basel, Switzerland
cSwiss Nanoscience Institute, University of Basel, 4056 Basel, Switzerland
dUniversitat Politècnica de Catalunya (UPC), Photovoltaic Lab – Micro and Nano Technologies Group (MNT), Electronic Engineering Department, EEBE, Av. Eduard Maristany 10-14, Barcelona 08019, Spain
eUniversitat Politècnica de Catalunya (UPC), Barcelona Centre for Multiscale Science & Engineering, Av. Eduard Maristany 10-14, Barcelona 08019, Spain
fGroup of Characterization of Materials, Departament de Física, Universitat Politècnica de Catalunya (UPC), Campus Diagonal-Besòs, Av. Eduard Maristany 10-14, Barcelona 08019, Spain
First published on 27th August 2025
Chalcohalides are an emerging family of semiconductors with tunable material properties, shaped by the intricate interplay between their unique structural chemistry and vibrational dynamics. Despite their promise for next-generation optoelectronics and solar energy conversion devices, their intrinsic optoelectronic properties remain largely unexplored. Here, we focus on the (Sb,Bi)(S,Se)(Br,I) system, a subset of compounds that share the same quasi-1D crystal structure. Using a two-step physical vapor deposition (PVD) process, we synthesize the eight ternary chalcohalide compounds, demonstrating bandgaps ranging from 1.38 to 2.08 eV with sharp, single-component photoluminescence (PL) peaks. In a parallel exploration of carrier dynamics and intrinsic electron–phonon interactions – comprehensively studied using power-, temperature-dependent, and time-resolved PL measurements – we map their direct impact on optoelectronic performance. Supported by first-principles density functional theory (DFT) defect calculations, we establish clear structure–property relations, identifying solid-solutions engineering as an effective means to fine-tune the native phonon structures and further suppress non-radiative recombination. This study provides a blueprint for optimizing chalcohalides as high-efficiency materials across a wide range of optoelectronic applications.
Chalcohalides are an emerging family of inorganic semiconductors that combine both chalcogen and halogen anions.16–21 The instability of halide perovskites is often attributed to the high electronegativity and ionic bonding character of halogen anions. In chalcohalides, however, the partial substitution of halogens with chalcogen anions – sometimes referred to as a “split-anion approach”22 – introduces a more covalent chemical bonding environment, which is predicted to enhance both chemical and thermal stability.23–25 The coexistence of divalent chalcogen and monovalent halogen anions also enables a high degree of tunability in structural and optoelectronic properties. Moreover, the incorporation of trivalent metal cations offers an electronic configuration similar to Pb2+, where the unique chemistry of the ns2 lone-pair state is known to play a crucial role in the exceptional defect tolerance of lead-halide perovskites.26–28 This suggests that chalcohalides, with their analogous electronic structure, may exhibit similar defect-tolerant properties. These distinctive electronic structures also feature highly dispersive band edges, leading to high carrier mobilities.29,30 Finally, by mitigating issues related to halogen volatility, chalcohalides are more thermodynamically stable than halide perovskites, further supporting their potential for integration into a wide range of optoelectronic devices.
Despite combining many of the most favorable material properties for optoelectronic applications, chalcohalides are still in the early stages of research. Unlocking their full potential and identifying the most promising candidates requires a deeper understanding of their fundamental properties. This includes not only advancing synthesis strategies to realize and optimize thin films, but also addressing the current lack of experimental studies on their optoelectronic properties. Moreover, given the vast compositional space accessible through different combinations of chalcogen and halogen anions with various metal cations, it would be more efficient to explore smaller subsets of chalcohalides – such as those sharing similar crystal structures – in parallel. Such targeted studies can accelerate the mapping of synthesis–structure–property relations, enabling precise functionality engineering.31 These mapped relations may ultimately serve as a blueprint for the inverse design of next-generation optoelectronic materials with tailored properties.
In this work, we investigate the optoelectronic properties of a combinatorial subset of heavy pnictogen-based ternary chalcohalides with the general formula (Sb,Bi)(S,Se)(Br,I). These eight compounds all crystallize in the same quasi-one-dimensional orthorhombic structure with space group Pnma, providing an ideal platform for comparative studies. We present a two-step synthesis route, in which chalcogenide precursors are first formed via thermal co-evaporation, followed by high-pressure reactive annealing to incorporate the halogen and form the final ternary phase. We first characterize their structural and morphological properties, alongside optical measurements to evaluate light absorption and emission behavior. Based on these results, we identify the most promising compounds for more in-depth photoluminescence (PL) studies, including power- and temperature-dependent measurements, as well as excitation intensity-dependent time-resolved PL spectroscopy (TRPL). These measurements are used to assess the nature of the emissive states, charge-carrier recombination dynamics, and electron–phonon interactions, offering insight into the fundamental processes that govern the intrinsic optoelectronic performance in this class of materials.
Following deposition, the chalcogenide precursor was placed in a Petri dish along with the corresponding halide source – SbI3 (Thermo Scientific, 99.999%), SbBr3 (Thermo Scientific, 99.995%), BiI3 (Puratronic, 99.999%), or BiBr3 (Thermo Scientific, 99%). For example, Sb2Se3 was combined with SbI3 to form SbSeI. The Petri dish was then sealed within a steel tubular furnace under an inert Ar atmosphere and subjected to high-temperature, high-pressure annealing. The annealing conditions were optimized for each chalcohalide composition, as detailed in Table 1. All purities specified are metals basis. The synthesis process has been schematically illustrated in Fig. S1. See ref. 32 for further synthesis details.
Compound | T (°C) | p (bar) | t (min) |
---|---|---|---|
BiSBr | 400 | 2.5 | 15 |
BiSI | 425 | 2.5 | 15 |
BiSeBr | 450 | 4.0 | 15 |
BiSeI | 500 | 4.0 | 15 |
SbSBr | 325 | 2.5 | 15 |
SbSI | 300 | 2.0 | 15 |
SbSeBr | 450 | 4.0 | 15 |
SbSeI | 450 | 3.0 | 15 |
Photoluminescence (PL) spectra were measured using a Horiba LabRam confocal microscope in a backscattering configuration with 538 nm and 488 nm excitation lasers. The excitation beam was focused onto the sample surface using a long-range 50× microscope objective with a numerical aperture of 0.55, yielding a beam diameter of approximately 1.1 μm for the 488 nm laser and 1.2 μm for the 532 nm laser. Temperature-dependent PL spectra were recorded from 70 K to 300 K using a MicroStat HiRes optical cryostat from Oxford Instruments with liquid nitrogen cooling, while excitation power-dependent PL studies were conducted at 70 K. The backscattered light was analyzed using a spectrometer equipped with a 150 g mm−1 holographic grating and a thermoelectrically cooled CCD.
Time-resolved photoluminescence (TRPL) measurements were conducted using a MicroTime 100 time-resolved confocal microscope coupled with a PicoQuant detection unit. Samples were excited using a pulsed 488 nm laser with a pulse duration of <100 ps and a repetition rate of 64 MHz. The excitation beam was focused onto the sample surface through a long-range 20× objective with a numerical aperture of 0.45. Emitted photons were collected through the same objective and guided via a 50 μm diameter optical fiber to the detection unit, which included a FluoTime 300 photospectrometer, a monochromator, and a photomultiplier detector. A cutoff filter was employed to suppress reflected laser photons. TRPL measurements were performed at room temperature.
Scanning electron microscopy (SEM) images were acquired using a Zeiss Gemini 460 equipped with an in-lens secondary electron detector, operated at an acceleration voltage of 2 kV. The SEM was also equipped with an Ultim Max 170 EDS system from Oxford Instruments, which was used to perform energy-dispersive X-ray spectroscopy (EDS) at an acceleration voltage of 10 kV. X-ray diffraction (XRD) patterns were measured using a Bruker D8 Advance in Bragg–Brentano geometry with Cu Kα radiation. UV-vis measurements were conducted using a PerkinElmer Lambda 950 UV/vis/NIR spectrophotometer equipped with an integrating sphere.
Ab initio calculations based on density functional theory (DFT)33 were carried out as implemented in the VASP code.34 The projector augmented-wave method35 was used to represent the ionic cores36 and for each element the maximum possible number of valence electronic states was considered. Wave functions were represented in a plane-wave basis typically truncated at 600 eV. By using these parameters and a dense k-point grid for Brillouin zone integration of 8 × 4 × 3 centered at Γ, the resulting zero-temperature energies were converged to within 1 meV per formula unit. In the geometry relaxations, a tolerance of 0.005 eV Å−1 was imposed in the atomic forces.
We employed the Perdew–Burke–Ernzerhof exchange–correlation functional revised for solids.37–39 Long-range dispersion interactions were taken into account through the van der Waals D3 correction scheme.40 Phonon calculations were performed with the small-displacement method and the PHONOPY software.41 Large supercells (324 atoms) constructed as 3 × 3 × 3 unit cell replications were employed along with a k-point grid of 4 × 2 × 2 for sampling of the Brillouin zone.
The selenides, BiSeBr and BiSeI, exhibit compositions closer to the expected stoichiometry, and their surface morphologies appear more homogeneous compared to their sulfur-based counterparts. As with the sulfides, all XRD reflections match the reference patterns from the ICSD database, suggesting phase-pure samples. Unlike prior reports on BiSeI, which describe crystallographic morphologies as platelet-like structures or rod-like microcrystals obtained via vapor transport and solvothermal synthesis,46,47 our synthesis approach resulted in compact thin films. This distinction is particularly relevant for optoelectronic device applications, where planar architectures are typically used. Given the limited number of studies on BiSeI and BiSeBr, and the fact that we obtained compact thin film morphologies for both of these compounds, investigating their optoelectronic properties is particularly interesting, as these materials can be readily integrated into various optoelectronic applications.
Among the Bi-based compounds, BiSI is the most extensively studied and is the only material to have been explored in various energy-related applications, including solar cells,48–51 photocatalysis,52 and photodetectors.53 However, despite its widespread study and use in energy devices, relatively little has been reported on its intrinsic optoelectronic properties. The large needle-like single crystals grown here offer a unique opportunity to study these optoelectronic properties in greater detail.
SEM images show that only SbSBr forms a compact thin film; however, it is also the most halogen-deficient and exhibits the strongest XRD reflections corresponding to the binary precursor phase. Although reflections corresponding to the desired ternary phase are observed, the lack of phase purity complicates the assessment of intrinsic material properties and their viability for device applications. This underscores the need to further optimize the synthesis process of the Sb-based chalcohalides to achieve complete halogenation of the binary precursors and obtain compact, phase-pure thin films.
It is important to note that no prior reports of SbSeBr exist in the ICSD database. Consequently, the reference powder pattern for this compound was generated using the computationally predicted structure retrieved from the Materials Project.54 However, the presence of secondary binary Sb2Se3 phases in the XRD patterns and significant uncertainty in the peak identification of the desired ternary phase indicate that results obtained from this sample should be interpreted with caution. Further studies are needed to confirm the structure of SbSeBr and optimize the synthesis process.
Among the Sb-based chalcohalides, SbSI is the most extensively studied. It was the first compound in this family identified as both ferroelectric and photoconductive, with its photoelectric properties investigated by Nitsche and Merz,55 and its ferroelectric behavior first reported by Fatuzzo et al.56 The discovery of ferroelectricity in SbSI sparked significant interest, particularly for applications in photodetectors,16,57–60 solar cells,61,62 photocatalysts,63–66 and nanogenerators.67 In contrast, SbSeI, SbSBr, and SbSeBr remain largely unexplored for energy devices, despite their structural similarity to SbSI. Investigating their optoelectronic properties and potential for functional devices could uncover new opportunities for this emerging materials class.
Across all eight chalcohalide compounds, we observe direct bandgaps ranging from 1.38 eV to 2.08 eV, a range that is highly desirable for various semiconductor and optoelectronic applications. While most chalcohalides are known to exhibit indirect fundamental bandgaps,16 the energy difference between the indirect and lowest direct transitions is typically so small that the indirect edge has minimal, if any, impact on most optical spectroscopy measurements. The absorption coefficients of the Bi-based compounds reach values on the order of 105 cm−1, highlighting their strong potential as efficient light-absorbing materials.
By comparing the absorption coefficients and PL emission peaks, we find that the absorption onset and PL peak position match closely for both the Bi-based and Sb-based chalcohalides. However, a small anti-Stokes shift is observed in SbSeI and SbSBr, where the emission peak is slightly shifted toward higher photon energies. This shift may be attributed to carrier localization effects, where shallow defect states or compositional inhomogeneities create a preferential recombination pathway at slightly higher energies. Alternatively, excitonic effects could contribute, as variations in exciton binding energies may introduce subtle differences between absorption and emission. The higher bandgap materials, BiSBr and SbSI, exhibit broader and more asymmetric emission peaks, which may arise from a combination of electron–phonon coupling and the influence of defect-related or band tail states. The presence of secondary phases in the Sb-based chalcohalides is not immediately apparent from the absorption and emission characteristics, as both Sb2Se3 and Sb2S3 have lower bandgaps of approximately 1.2 eV68 and 1.7 eV,69 respectively. However, their presence could contribute to the observed peak broadening.
With the exception of SbSeBr, the substitution of sulfur with selenium generally lowers the bandgap, a trend attributed to variations in the lattice parameters. However, as previously mentioned, the structural determination of SbSeBr was inconclusive, and therefore, the observed deviation from the sulfur–selenium trend should be treated with caution. These measurements may not accurately reflect the intrinsic properties of the desired ternary compound, and as such, the results have been labeled with an asterisk. Since all eight chalcohalide compounds share the same crystal structure, it is also worth exploring solid solutions of both the chalcogenide and halogen anions. Similar to other optoelectronic materials, such as Cu(In,Ga)Se2 (CIGS) and Cu2ZnSn(S,Se)4 (CZTS), where solid solutions have been widely used, chalcohalide solid solutions could provide a means to fine-tune material properties, including the bandgap, further enhancing their versatility for device applications.
First, we investigated the nature of the emission peaks of the two thin films through power-dependent PL measurements at 70 K using both 488 nm and 532 nm excitation lasers. This analysis helps determine whether the significantly weaker PL intensity observed in the thin films, relative to the needle-like BiSI crystals, arises from excitonic effects or defect-related transitions. Fig. 4(a) shows the power-dependent PL spectra of BiSeBr, while Fig. 4(b) and (c) present the integrated PL intensity, determined by fitting the emission peak with a single Gaussian, as a function of excitation power density for BiSeBr and BiSeI, respectively. The measurements are fitted using the power-law relation I ∝ Pk, where I is the integrated PL intensity, P is the excitation power density, and k is the power-law exponent.
![]() | ||
Fig. 4 Photoluminescence characterization of BiSI, BiSeI, and BiSeBr. (a) Power-dependent PL spectrum of BiSeBr at T = 70 K. (b and c) Power dependence of the integrated PL intensity for the two selenide-based thin films, measured using 532 nm and 488 nm excitation lasers, respectively, and fitted to a power-law relation. (d) Temperature-dependent PL spectra of BiSI and BiSeBr (inset), highlighting the significant difference in relative PL yield between the needle-like crystal and the thin films. BiSeI exhibits similar PL yields to BiSeBr (see Fig. S3 in the SI). (e) Temperature dependence of the fitted PL emission peak positions, modeled using eqn (1). (f) Temperature dependence of the peak width, modeled using eqn (2). (g) Time-resolved PL (TRPL) at room temperature, with transient decays fitted using a bi-exponential function. (h and i) Excitation intensity dependence of the fast and slow decay components, respectively. BiSeBr does not exhibit a slow decay component. |
For both thin films, we consistently obtain k values close to 1 under both excitation wavelengths, indicating that the PL emission peaks originate from band-to-band-like transitions. This behavior rules out defect-mediated recombination, which typically results in sub-linear (k < 1) power dependencies, as well as excitonic transitions, which are expected to yield super-linear (k > 1) responses.70,71
Having established that the PL emission originates from intrinsic, band-to-band-like transitions, we next investigate the temperature-dependent PL. As shown in Fig. 4(d), the relative PL yield is significantly higher for the needle-like BiSI crystals compared to the BiSeBr and BiSeI thin films (the latter shown in Fig. S3 in the SI) across all measured temperatures.
Se-based compounds typically exhibit lower PL yields than their S-based counterparts due to a combination of intrinsic electronic and structural factors. The larger atomic radius and weaker metal–selenium bonds in Se compounds tend to increase the density of structural defects and non-radiative recombination centers, both of which can suppress radiative emission. Recent work by López et al.72 demonstrated that Se-based pnictogen chalcohalides (e.g., BiSeI and BiSeBr) exhibit higher equilibrium concentrations of native defects – particularly selenium vacancies – than their sulfur-based analogues (e.g., BiSI and BiSBr). This trend is attributed to the lower formation energies and high volatility of selenium under typical synthesis conditions. Furthermore, the reduced bandgaps of Se-based compounds favor thermal population of non-radiative states, which can enhance phonon-assisted decays.73,74 The stronger spin–orbit coupling in selenium, due to its higher atomic mass, can also facilitate exciton relaxation into dark states, further reducing radiative efficiency. These effects have been widely observed in other systems. For instance, monolayer MoSe2 exhibits significantly lower PL intensity than MoS2 under identical conditions, primarily due to increased non-radiative decay pathways and dark exciton formation.75 Similarly, CdSe quantum dots without core–shell passivation generally show lower PL quantum yields than CdS counterparts,76 highlighting a broader trend of diminished PL efficiency in Se-based semiconductors.
This difference in PL yield may also be partly attributed to a deeper excitation penetration depth in BiSI, coupled with a much higher density of grain boundaries in the thin films, which are expected to act as non-radiative recombination centers. This suggests that the absolute PL intensity and thus the overall optoelectronic quality of the selenide-based thin films could be improved by increasing the crystal grain size or through grain boundary passivation.
To quantify the temperature-dependence of the PL emission, we fit the peak at each temperature with a single Gaussian (R2 > 0.95) and track the extracted peak position and spectral width. These trends are shown in Fig. 4(e) and (f), respectively. As both the thermal shift of the bandgap and the broadening of the emission peak are typically attributed to electron–phonon interactions, we model them using Bose–Einstein statistics. The temperature dependence of the bandgap is described by:
![]() | (1) |
![]() | (2) |
The temperature-dependent PL behavior of BiSeBr and BiSI is well described by the Bose–Einstein models, as reflected by their modest coupling strengths (Γph = 50 meV and 36 meV, respectively) and characteristic phonon energies (ℏωph = 31 meV and 26 meV). These relatively low values suggest that PL broadening in these materials is primarily governed by moderate electron–phonon interactions, consistent with weakly perturbed radiative recombination. The small values of λ and ℏ further indicate that the thermal redshift of the PL peak is gradual, and that the underlying phonon modes responsible for broadening are within a physically reasonable energy range. Additionally, the higher bandgap energies (1.54 eV for BiSeBr and 1.63 eV for BiSI) help suppress the thermal activation of non-radiative states, consistent with the defect tolerance of these materials recently reported by López et al.72
In contrast, the fitted parameters for BiSeI show a significantly higher phonon coupling strength (Γph = 161 meV) and an unusually large characteristic phonon energy (ℏωph = 71 meV), well beyond the typical energy range of lattice vibrations. This suggests that while the Bose–Einstein model generally describes the temperature dependence of PL in BiSeI, it likely oversimplifies the true broadening mechanisms, potentially involving anharmonic lattice dynamics or multiple phonon modes. These anomalies align with the recent insights into the native defect chemistry of this materials system from López et al.,72 who showed that the equilibrium concentration of selenium vacancies (VSe) in BiSeI is relatively high due to a low formation enthalpy under Se-poor conditions. These vacancies introduce deep charge transition states that drive non-radiative recombination, significantly reducing photocarrier lifetimes and radiative efficiency. Although the capture coefficients of these defects are calculated to be modest, their high concentrations can substantially limit device performance. The detrimental impact of these vacancies may be mitigated through cation-poor synthesis conditions and strategic anion substitutions, offering potential pathways to improve the optoelectronic performance of BiSeI. In contrast, BiSeBr benefits from higher defect formation enthalpies – particularly for antisite defects – due to larger ionic radius mismatches, which effectively suppress defect formation. Collectively, these results illustrate that the stronger apparent phonon coupling and enhanced PL broadening observed for BiSeI not only reflect more intense electron–phonon interactions but also highlight the severe impact of high defect densities on non-radiative carrier decay.
To assess whether the carrier dynamics in these compounds reflect similar behavior as suggested by the temperature-dependent PL results, we performed time-resolved photoluminescence (TRPL) measurements, shown in Fig. 4(g). The PL decay in BiSeBr can be modeled using a mono-exponential function, indicating that the carrier dynamics in this sample are relatively straightforward to describe. In contrast, the transient decays for BiSeI and the needle-like BiSI crystal require at least bi-exponential fits, suggesting multiple recombination pathways occurring in parallel and/or carrier localization effects.
To better understand these processes, we performed excitation intensity-dependent TRPL measurements to study how the different decay components evolve with increasing injection levels. The results are shown in Fig. 4(h) and (i) for the fast and slow decay component, respectively. As non-radiative recombination centers become saturated at higher carrier densities, the effective carrier lifetime is expected to increase asymptotically.77 This trend is indeed observed for the fast decay component in both BiSI and BiSeBr, suggesting that this component reflects the effective carrier lifetime in these materials. Recent work by Yuan et al. showed that shallow defects typically limit the recombination lifetime in triple-cation perovskites, leading to a decrease in the measured lifetime with increasing carrier concentration.78 Since none of the decay components in any of our samples exhibit this trend, shallow traps are unlikely to be the dominant non-radiative recombination mechanism. This implies that the effective carrier lifetime may be significantly improved by identifying and mitigating deep-level defects.
In contrast, the fast decay component in BiSeI increases approximately linearly with excitation intensity, while the slower component shows little dependence on injection level. This behavior aligns with the temperature-dependent PL results, indicating that the recombination dynamics and electron–phonon interactions in this sample are not well-captured by simple models and point to relatively poor optoelectronic quality. Nevertheless, the fast decay component remains below 1 ns for all three samples, highlighting the need for defect mitigation strategies to enhance carrier lifetimes.
In contrast, BiSeI stands out for its anomalously strong phonon coupling (Γph = 161 meV) and a fitted characteristic phonon energy (ℏωph = 71 meV) that is well beyond the physical range of known vibrational modes. While the Bose–Einstein model captures the general temperature trend, these values suggest that the model oversimplifies the underlying recombination mechanisms. The exaggerated coupling suggests a combination of anharmonic lattice effects, multi-mode phonon interactions, and substantial carrier trapping by deep defects.
Recent defect calculations by López et al.72 shed light on the native defect chemistry, identifying chalcogen vacancies as the killer point defects in these systems. BiSeI is particularly susceptible to form Se vacancies (VSe) due to a relatively low formation enthalpy under Se-poor conditions. These vacancies act as deep non-radiative centers, consistent with the poor PL performance and sub-nanosecond carrier lifetimes observed in our measurements. In contrast, BiSeBr exhibits higher antisite formation energies, partially attributed to larger ionic radius mismatches, which enhance its defect tolerance and radiative efficiency.
TRPL measurements further confirm these findings. BiSeBr exhibits clean mono-exponential decay, indicative of a single dominant recombination pathway consistent with efficient band-to-band emission. In contrast, BiSeI and BiSI show bi-exponential decay, suggesting the presence of multiple recombination channels and possible carrier localization. In BiSeI, the fast decay component increases linearly with excitation density and remains below 1 ns across all measured conditions, consistent with a high density of unsaturable non-radiative traps. BiSeBr and BiSI, however, exhibit an asymptotic increase in their fast decay components, consistent with saturation of defects. These trends support the notion that deep-level defects govern recombination in BiSeI, while BiSeBr may be further optimized through surface or grain boundary passivation.
In addition to defect chemistry, vibrational effects also play a central role in the optoelectronic performance of these chalcohalides. Fig. 5 shows the total phonon density of states (DOS) for BiSI, calculated using density functional theory (DFT), revealing a pronounced phonon gap just below the sulfur-related phonon band (25–30 meV), which closely aligns with the experimentally extracted characteristic phonon energy ℏωph. This gap originates from the large atomic mass differences between the heavy pnictogen and the much lighter sulfur atoms and is consistently observed across all sulfide-based compounds (Fig. S4 in the SI). While this gap can reduce phonon scattering pathways – potentially improving radiative efficiency – it may also enhance energy dissipation through specific high-energy modes, promoting non-radiative recombination. In contrast, the corresponding selenium-based compounds lack this phonon gap, as selenium-related phonon modes extend into the previously forbidden vibrational range. This filled gap intensifies carrier-phonon scattering, providing a compelling vibrational explanation for the anomalous electron–phonon coupling strength observed in BiSeI.
Forming solid solutions, such as Bi(S,Se)I, may help mitigate these issues. Partial substitution of S with Se or Br with I offers a means to engineer the phonon landscape, reintroducing beneficial gaps or shifting scattering bands to suppress non-radiative recombination pathways. Notably, López et al.72 reported that such substitutions also improve intrinsic point defect tolerance, providing a dual strategy for optimizing optoelectronic performance through combined phonon and defect engineering.
Altogether, these results highlight the importance of both intrinsic material design and extrinsic processing conditions in tuning the optoelectronic properties of chalcohalides. While BiSeI is currently limited by strong electron–phonon coupling and a high equilibrium concentration of native defects, BiSeBr and BiSI present promising platforms for further development. Strategies such as targeted alloying, grain boundary passivation, and optimization of growth conditions will be key to advancing this materials class toward high-efficiency applications in photovoltaics and photodetection.
Our results establish a clear structure–property relationship within this material family, where defect formation energies and phonon dynamics play a critical role in determining optoelectronic performance. We suggest that BiSeBr and BiSI hold promise for further development, while BiSeI requires defect management to improve its radiative efficiency. Strategies such as solid-solutions engineering (e.g., Bi(S,Se)I), optimizing growth conditions, and implementing surface or grain boundary passivation can mitigate defect-induced losses and tune the native phonon landscape. By demonstrating how intrinsic defect chemistry and lattice dynamics govern light emission and carrier lifetimes, this study provides a roadmap for optimizing chalcohalides for high-performance applications in photovoltaics, photodetectors, and beyond.
The SI provides experimental details, synthesis procedures, structural and morphological characterization, photoluminescence, phonon density of states, and bandgap data for the eight (Sb,Bi)(S,Se)(Br,I) chalcohalides studied. See DOI: https://doi.org/10.1039/d5ta05011a.
This journal is © The Royal Society of Chemistry 2025 |