Open Access Article
Oleksandr Dobrozhan†
ab,
Sajida Kousar†a,
Girish C. Tewari
a and
Maarit Karppinen
*a
aDepartment of Chemistry and Materials Science, Aalto University, FI-00076 Espoo, Finland. E-mail: maarit.karppinen@aalto.fi
bDepartment of Electronics and Computer Technology, Sumy State University, 40007 Sumy, Ukraine
First published on 2nd April 2026
Metal chalcogenides represent state-of-the-art thermoelectric (TE) materials, both in practical TE energy conversion devices and in new TE material research. To further push the performance boundaries of these materials, the major challenge is to suppress the thermal conductivity as much as possible while maintaining high electrical transport properties. The strategies employed to address this challenge include band engineering, chemical composition tuning, and nano- and microstructural modification. Regarding chemical composition tuning, an effective approach specific to metal chalcogenide thermoelectrics is anion mixing through solid solution formation at a chalcogen site, where sulfur, selenium, and tellurium atoms can occupy the same crystallographic site. This perspectival review aims to provide first a wider overview of such possibilities among the metal chalcogenides in general and then highlight the recent research on three currently strongly emerging TE metal chalcogenide families: (i) different binary, ternary and quaternary copper chalcogenides, (ii) ternary cobalt antimony chalcogenides, and (iii) binary silver chalcogenides.
In recent years, various new metal chalcogenide families have been highlighted as prominent TE material candidates. Compared to metal oxides, the higher bond covalency in metal chalcogenides is beneficial for enhancing the power factor (PF), while the heaviness of the chalcogen atoms results in lower thermal conductivity, such that overall, the TE figure-of-merit (ZT) is enhanced. Also importantly, metal chalcogenides can be doped into both n-type and p-type conductors and engineered into different kinds of structures; these are essential requirements for the fabrication of actual TE devices. Moreover, by controlling the anion composition in these materials, their structural, electrical, and thermal properties can be optimized.
Indeed, an effective approach specific to metal chalcogenide thermoelectrics is anion mixing via solid solution formation at the chalcogen site occupied by sulfur, selenium or tellurium, as opposed to mixed-anion frameworks where anions occupy distinct crystallographic positions. In the resultant metal chalcogenides, the cations are bonded to more than one kind of anion, inducing disorder in the lattice and suppressing thermal conductivity, while also providing a powerful tool for fine-tuning the electronic structure for improved electronic transport, thereby enhancing the ZT value. In Fig. 1, we summarize the expected benefits achieved by the anion mixing approach in the metal chalcogenides.
This perspectival review aims to provide, after a brief introduction to the basics of thermoelectrics in Section 2, first a wider overview of commonly used engineering approaches for tuning charge and heat transport in state-of-the-art metal chalcogenides (Section 3) and then highlight recent research on three currently strongly emerging TE metal chalcogenide families: different binary, ternary and quaternary copper chalcogenides (Section 4), ternary cobalt antimony chalcogenides (Section 5), and binary silver chalcogenides (Section 6). Finally, we make a brief outlook in Section 7.
The interrelationship between electrical and thermal transport properties poses a major challenge in thermoelectric materials. Electrical conductivity depends on the charge carrier concentration (n), as expressed by: σ = neμ, where e is the elementary charge and μ is the carrier mobility. Meanwhile, the electronic contribution to thermal conductivity follows the Wiedemann–Franz law: κe = L0 σT, where L0 is the Lorenz number. The total thermal conductivity can then be written as: κ = κl + L0neμT. While increasing the carrier concentration enhances electrical conductivity, it also raises thermal conductivity, thus reducing the ZT value. This trade-off highlights the difficulty of designing high-performance TE materials. To improve the ZT value, researchers typically aim to enhance the power factor (PF = S2σ), while minimizing the lattice thermal conductivity through band structure engineering and/or nano- and microstructural modifications.2–4
One strategy to enhance the thermoelectric performance is to reduce the lattice thermal conductivity without negatively impacting the electrical transport properties. This has been achieved by incorporating heavy elements such as lead, thallium, bismuth, antimony, and tellurium into the TE materials, as well as designing complex crystal structures that increase phonon scattering. However, complex structures may reduce carrier mobility, thereby negatively affecting electrical conductivity.1,5
Thermoelectric materials have found applications in various fields, spanning from small-scale power generation to advanced aerospace and defense technologies. One of their most promising uses is in waste heat recovery, where they enhance energy efficiency by capturing excess heat from industrial processes, automotive engines, and power plants, converting it into usable electrical energy.6–9 Similarly, wearable thermoelectric generators (TEGs) have been developed for self-powered electronics, providing sustainable energy solutions for biomedical sensors, smart textiles, and personal gadgets.10 Thermoelectric technology is vital in aerospace and defense, providing reliable power for deep-space missions and nuclear detectors operating in extreme environments.11,12
An important advantage of TE modules is their solid-state operation, which makes them compact and durable.13 Moreover, TE devices offer environmental benefits, as they operate without producing greenhouse gas emissions or relying on harmful chemicals, making them a sustainable and eco-friendly option for energy conversion. Their efficiency, reliability, and versatility continue to drive research efforts, particularly in the field of chalcogenide-based TE materials, which are being optimized for even greater performance with high mechanical and chemical stability in next-generation energy applications.
Extensive research has been devoted to optimizing the TE performance of PbTe through band convergence,15,16 nanostructuring,23 dislocation control,24 doping,25,26 and alloying.27,28 The NaCl-type structure allows doping/alloying with compatible elements/compounds. In particular, cation doping with Mg and Se improved the ZT to 2.2 at 820 K,28 with Na, Eu, and Sn to 2.5 at 900 K,29 and with Zn and I to 1.35 at 650 K.30
On the other side, PbSe shows promising TE properties, outperforming PbTe at 900 K.31 Even though intrinsic PbSe suffers from low n and high κ, these parameters can be improved by doping polycrystalline PbSe with Na (ZT = 1.2 at 850 K),32 In (ZT = 1.6 at 850 K),33 Cu and Te (ZT = 1.7 at 773 K),34 or combining Br-doping and Cu2Se-nanomixing (ZT = 1.8 at 723 K).35 Recently, lillianite-type minerals have shown great promise as TE materials, reaching a peak ZT of 1.35 at 800 K for polycrystalline n-(Pb0.95Ga0.05)7Bi4Se13 samples, attributed to unusually high valley degeneracy due to alignment of nested conduction bands, strong phonon anharmonicity, and nanoscale lattice imperfection.36
Similar strategies have been applied to enhance the TE performance of polycrystalline PbS, which is also characterized by a low n and high κ. For example, adding CdS and ZnS nano-inclusions enhanced ZT up to 1.3 at 923 K,37 Cu2S, SrS and CaS in Na-doping – up to 1.2 at (823–923) K,38,39 GeS and Sb-doping – up to 1.32 at 923 K,40 (Sb, Cu) co-doping – up to 1.23 at 923 K,41 (Sb, Cl) co-doping – up to 1.0 at 823 K,42 and Cl-doping – up to 0.7 at 850 K.43 Compositing is also an effective strategy for balancing electronic and thermal transport in PbS-based materials. For example, incorporating Cu or Ag can strongly increase the charge carrier concentration through spillover from metallic inclusions, while the presence of multiple phases reduces the κl, yielding ZT values of 1.1 at 750 K and 1.7 at 900 K.44,45
Anion mixing has been shown to be effective in enhancing the TE performance of polycrystalline Pb-based materials. For instance, anion mixing in p-type Pb(Te,Se) led to a ZT of 1.8 at 850 K due to high valley degeneracy, enhanced alloy scattering, and increased solubility limit for cation doping.46 Similarly, Se substitution at the S site has been shown to be beneficial as it was found to increase bond covalency, modify effective mass, and weaken electron–phonon coupling, thereby improving charge carrier mobility. At the same time, it introduces point defects that scatter high frequency phonons. As a result, a ZT value of 1.7 at 900 K was achieved for (Pb0.93Sb0.05)(S0.5Se0.5).47 Notably, entropy-driven engineering in Na-doped (PbTe)0.84(PbSe)0.07(PbS)0.07 alloys outperforms some endotaxial nanostructured counterparts due to the efficient reduction of κl through alloy scattering and point defects, together with simultaneous valence band modification, yielding a ZT of 2.0 at 800 K.48 Moreover, the incorporation of PbSe into the ternary PbTe–PbS system broadens the solubility limit for anion doping, enabling a wider carrier concentration range. It also provides an effective means of tuning the band gap, while the strong atomic mass contrast enhances phonon scattering, resulting in an excellent ZT of 1.1 at 800 K for (PbTe)0.75(PbS)0.15(PbSe)0.1,49 1.52 at 700 K for (PbTe0.9295Se0.07Cl0.0005)0.93(PbS)0.07,50 and 1.6 at 750 K for (Pb0.98Na0.02)(Se0.1S0.05Te0.85).51 Another notable example is (Pb0.98Na0.02)(Te0.75Se0.2S0.05), in which high-temperature band convergence and nanoprecipitations resulted in an exceptionally high ZT of 2.3 at 800 K.52
Various doping strategies have been explored to enhance the TE efficiency of Bi2Te3, particularly effective being the cationic Sb-for-Bi and anionic Se-for-Te substitutions:58 (Bi,Sb)2Te3 compositions typically exhibit p-type behavior, while Bi2(Te,Se)3 compositions favor n-type conduction.59 Advancements in nano- and microstructural engineering have significantly improved the TE properties, with ZT reaching 1.86 at 373 K in p-type (Bi,Sb)2Te3 owing to enhanced phonon scattering and increased carrier mobility.60,61 Additionally, n-type nanoplates Bi2(Te0.9Se0.1)3 synthesized using a microwave-assisted method have demonstrated ZT = 1.23 at 480 K, with a remarkably reduced thermal conductivity of 0.98 W m−1 K−1 due to grain boundary scattering, dislocations, and point defects.62 Notably, Bi2Te3-based thin films have shown significantly enhanced ZT values compared to bulk Bi2Te3 alloys, with the highest values reported at 300 K of 2.75 for n-type nanoparticle-based Bi2(Se0.4Te0.6)3 thin films and 2.4 for p-type Bi2Te3/Sb2Te3 superlattice structures.63,64
The Se-counterpart, Bi2Se3, crystallizes in a layered trigonal structure with a narrow bandgap (0.3 eV) but different electronic transport behavior in comparison with Bi2Te3 due to differences in defect chemistry and carrier mobility. The TE properties of Bi2Se3 are significantly influenced by intrinsic defects, primarily Se vacancies, which result in n-type conductivity.65 Various cation dopants (e.g. Cu, In, Pb) have been explored to modulate n and optimize the TE performance in single crystals of Bi2(Se,Te)3.66–68 However, Bi2Se3 has still remained inferior to Bi2Te3, due to its higher lattice thermal conductivity.
Although SnSe has a high ZT, its optimization possibilities are limited by the low-temperature Pnma phase, which resists effective doping due to its layered structure and strongly covalent, locally distorted bonding. However, alkali metals can act as effective acceptor dopants in SnSe; e.g. Na-doping in single-crystal SnSe increased hole concentration by two orders of magnitude and ZT up to 2.0 at 773 K.73,74 Practically high performance in p-type SnSe stems from its complex electronic band structure, while n-type SnSe achieves ultrahigh ZT through phonon–electron interactions and 3D-charge/2D-phonon transport.75,76
Due to the high cost and slow growth of single-crystal SnSe, research has shifted toward polycrystalline SnSe.77 However, its high resistivity from low intrinsic carrier concentration makes doping essential. Alkali metal (Na, K) doping enhanced conductivity and the Seebeck coefficient via band convergence, achieving ZT = 0.8 at 800 K with just 1% doping, while Zn doping further increased ZT to 0.96 at 873 K.78 Cu interstitials effectively reduced κl, yielding ZT = 1.2 at 873 K.79 For highly textured polycrystals, ZT was increased to 1.05 at 873 K.80 Textured SnSe ingots with large single-crystal domains showed suppressed κl, achieving a ZT of 1.1 at 873 K.81
For the sulfide, SnS, a ZT value of 0.16 at 823 K was reported for polycrystals.82 In SnS single crystals, doping increased carrier concentration, boosting the ZT to 1.1 at 870 K with Na-doping and up to 1.6 at 873 K with Se alloying, mainly due to tailoring the electronic band structure, particularly regulating effective mass and mobility, and activating Fermi pockets.83–85
The telluride counterpart, SnTe, is a narrow-bandgap (0.13 eV) semiconductor; it is considered a promising, eco-friendly alternative to PbTe due to its similar band structure. However, pristine SnTe suffers from high κl and a large energy gap (0.35 eV at 300 K) between its light-hole (L) and heavy-hole (Σ) bands. Thus, achieving effective band convergence requires high dopant concentrations.86 To address this, various doping strategies have been explored. While In, Cd and Hg showed limited solubility,87–89 Mn proved to be more effective. Mn-doped polycrystalline SnTe exhibited enhanced S (80 µV K−1 at 300 K; 230 µV K−1 at 900 K), increased power factors, and a ZT of 1.3 at 900 K, due to band convergence and point defect scattering.90 Alloying AgBiTe2 into polycrystalline SnTe compensates the hole concentration via Bi electron donors and reduces thermal conductivity through point defect scattering and Ag-rich nanostructuring, yielding a ZT value of 1.1 at 775 K.91
Among them, anion mixing represents one of the most effective strategies for defect engineering and band manipulation. For instance, isoelectronic doping with S atoms at the Se site was demonstrated to suppress the annihilation of donor defects (VIn and Ini), yielding a peak ZT of 0.67 in polycrystalline α-In2S0.05Se2.95 at 923 K.93 Furthermore, Se substitution in polycrystalline In2.67S4 thiospinels enhanced the charge carrier concentration – a well-known limiting factor in In–X materials – and promoted stabilization of the thermoelectrically favorable α-phase. This was attributed to the simultaneous reduction of both the band gap and the effective electron mass. Additionally, structural disorder in the ternary thiospinel contributed to a reduction in κ. Collectively, these effects resulted in an approximately 2.5-fold improvement in ZT, reaching a maximum of 0.25 at 760 K in In2.67S3.5Se0.5.94
Besides that, anion substitution of Se by Te in polycrystalline In0.96Bi0.04Se0.97Te0.03 introduced effective phonon scattering centers, thereby reducing the κl. The synergistic effect of cation–anion co-doping in this system yielded a maximum ZT of 0.13 at 630 K.95 Moreover, Te substitution can also improve the PF through an increase in the effective mass; as a result, ZT was enhanced more than twofold, reaching 0.096 at 700 K in In0.9Si0.1Se0.93Te0.07.96 Anion substitution with 7% S or Te in polycrystalline InSe alloys effectively increased σ with a more pronounced effect observed for Te-doped samples owing to their smaller band gap. Both dopants also reduced κ via point defect phonon scattering. Consequently, ZT reaches maxima of 0.13 and 0.28 at 735 K for InSe0.93S0.07 and InSe0.93Te0.07.97
Binary Cu2X and ternary CuCrX2 chalcogenides containing monovalent copper have been widely investigated because of their high ionic conductivity and ultralow κl, arising from the liquid-like behavior of Cu+ ions. Structurally, binary Cu2X compounds comprise a rigid X-atom framework responsible for electronic transport and a disordered, mobile Cu+ sublattice that suppresses transverse and shear vibrations, thereby reducing κl.104–107 Indeed, binary Cu2X (X = S, Se, Te) chalcogenides exhibit phase transitions and ultralow lattice thermal conductivity (<0.8 W m−1 K−1) despite the low atomic mass of Cu. Cu2S, for example, transforms from monoclinic γ to hexagonal β at 370 K and then to superionic cubic α at 700 K.108 The α-phase contains an fcc S2− framework and highly disordered Cu+ ions, resulting in superionic conductivity. Cu2S is an intrinsic semiconductor, while Cu deficiency readily induces p-type behavior. Many Cu2S-based compositions display excellent TE properties with ZT > 1.0, and Cu2S:3%Pb reaches ZT = 2.0 at 880 K.109 Stability-oriented strategies, such as structural stabilization and suppression of sulfur loss and Cu segregation, also promote high TE performance and PF.110 Alloying with Te and introducing nanoscale mosaicity in polycrystalline Cu2(S0.52Te0.48) produced a high ZT of 2.1 at 1000 K, mainly due to its exceptionally low κl of 0.3–0.4 W m−1 K−1. Mosaicity offers a distinct route for tuning TE properties by enabling efficient electron transport along quasi-single-crystalline frameworks while strongly scattering phonons at strained lattice regions and mosaic nanograin interfaces.111
In contrast, Cu2Se is an intrinsically p-type conductor with a bandgap of 1.2–2.3 eV.112 Near 400 K, α-Cu2Se transforms to the fast ion-conducting β-phase, where an fcc Se framework and disordered Cu+ ions give rise to liquid-like thermal conductivity (0.34–0.60 W m−1 K−1) while maintaining high electrical transport.113–119 Its TE performance has been markedly enhanced since 2012, reaching a ZT of 1.8 by self-propagating synthesis and melt-quenching and a ZT of 2.6 at 850 K for Cu2Se:1 mol% In.120–122 Yet, instability under strong electric fields and large temperature gradients still limits applications, although operation below a critical voltage and Cu excess can improve stability.123,124 By employing a melt-solidification technique, Zhao et al.125 improved ZT up to 1.9 at 973 K in polycrystalline samples, though stability issues persisted due to high ion mobility. To address these challenges, Cu2(Se,S) solid solutions were developed, where anion mixing enhances structural stability. Besides the intrinsically low thermal conductivity originating from the liquid-like Cu+ sublattice, mixed-anion solid solution formation further reduces lattice thermal conductivity via defect-induced phonon scattering and simultaneously optimizes carrier concentration. A ZT of 2.0 at 1000 K was achieved for polycrystalline Cu2(Se0.92S0.08) with improved stability, owing to stronger bonding that fixed Cu atoms in the lattice for suppressed Cu vacancy formation and reduced the transverse sound velocity, thus leading to optimal carrier concentration and lowered thermal conductivity.126 Similar behavior was observed by Zhao et al. in polycrystalline Cu2(Se0.8S0.2), where point-defect phonon scattering in the solid solution resulted in extremely low lattice thermal conductivity (0.2–0.4 W m−1 K−1) and a peak ZT of 1.65 at 950 K.127 A peak ZT of 2.3 at 1000 K was achieved in polycrystalline Cu1.94(Se0.5S0.5).128
The crystal structure of CuxTe remains debated due to the presence of closely related Cu2Te, CuTe, and CuTe2 phases at different temperatures.129,130 However, its cubic ε-phase above 850 K exhibits promising TE performance, reaching ZT = 0.29 in undoped Cu2Te and 1.0 in Ag-doped Cu2Te at 900 K.129 An ultralow κl (0.17 W m−1 K−1) has been observed in Cu5TeS3I3 polycrystals via coordination design for partial ionization of the Cu–I bonds, leading to reduced transverse speed of sound (νT = 839 m s−1) and high anharmonicity (γ = 2.76).131
Similarly, ternary CuCrX2 (X = S, Se) compounds exhibit κ owing to Cu+ ion mobility.132 They adopt a hexagonal layered structure with Cu+ ions located at tetrahedral sites between CrX2 layers.133–140 In CuCrSe2, Cu+ redistribution above 365 K induces kinetic disorder, markedly lowering κl while maintaining high σ; a similar effect occurs in CuCrS2 above 675 K.134,135,138,141–143 CuCrSe2 further shows promising transport properties and reaches ZT = 1.0 at 773 K.132,138,144–146 For CuCrS2, a high ZT value of 2.0 at RT was initially reported.147 The subsequent attempts by Kaltzoglou et al.148 to reproduce these results under analogous synthesis conditions yielded a ZT of only 0.023, suggesting that the originally reported value may be sensitive to differences in the microstructure or measurement methodology. Despite their promising TE properties, phase-pure CuCrSe2 is difficult to synthesize by conventional methods because of secondary phases such as CuCr2Se4.146,149 Accordingly, high-temperature treatment with water quenching, mechanical alloying followed by annealing, and spark plasma sintering (SPS) have been used to suppress impurity formation.139,150,151 CuCrX2 compounds are also highly sensitive to synthesis temperature because of native defect formation; sintering at elevated temperatures improves TE performance via partial Cr redistribution.133,135–137 We further found that Se/S anion mixing in CuCr(S,Se)2 drives a semiconductor-to-metal transition, narrows the bandgap, increases n, and lowers κ, all of which favor high TE performance.152 Moreover, Se-for-S substitution changes the magnetic state from long-range antiferromagnetism to a spin-glass-like regime with gapless spin-liquid-like excitations.153
Beyond the Cu2X and CuCrX2 families, several other multicomponent Cu-based chalcogenides have emerged as promising TE materials. For example, BiCuSeO is a notable quaternary p-type thermoelectrics.154 Its low κ enables high TE performance, with ZT reaching 1.1 through co-doping with (Ba, Ca, Pb)155–157 and anion mixing with Te.158 Moreover, argyrodite-like materials with multiple chalcogenide anions, e.g. bulk Cu8GeSe4Te2,159 have emerged as promising PLEC (phonon-liquid electron-crystal) candidates, characterized by dynamically disordered crystal structures that give rise to ultralow κl. Another example is polycrystalline low-symmetry Cu2Hf0.8Ti0.2Te3, in which electronic band nesting and hierarchical bonding formation collectively enhance the PF and suppress κl, yielding a ZT peak of 1.11 at 873 K.160 Also, anion mixing in argyrodite-type copper chalcogenides is an effective strategy for stabilizing the high-symmetry cubic phase at lower temperatures through configurational entropy engineering.161 This structural stabilization promotes a multivalley electronic structure of the valence band, modulates carrier concentration, and reduces κl up to 0.2–0.3 W m−1 K−1, collectively enhancing the TE performance. Representative examples include Cu8−xGeS4.36Te1.64 (ZT = 0.20 at 670 K),162 Cu7P(S,Se)6 (ZT = 0.75 at 673 K),163 and Cu8SiS3Se3 (ZT = 1.24 at 773 K)164 among others.
Antimony-containing compounds such as Cu3SbSe4, CuSbSe2, and Cu3SbSe3 demonstrate the important role of Sb3+ lone-pair electrons in lowering κ.165 Cu3SbSe4, for example, shows strongly reduced κ at elevated temperatures, and its TE performance can be enhanced to ZT = 0.7 by Sn doping.166 Other promising Cu-based chalcogenides include CuGaTe2, CuInTe2, Cu2Zn(Sn0.9In0.1)Se, and Cu2(Sn0.9In0.1)Se3, with reported ZT values approaching or exceeding unity.167–171 Among copper sulfides, tetrahedrites, colusites, and bornite also show attractive TE performance owing to intrinsically low thermal conductivity, with the best reported ZT values of 1.0–1.2.172–175
Particular attention in multicomponent chalcogenides has been paid to the role of crystal structure and phase transitions in TE performance. For example, Cu2SnSe3 adopts a monoclinic structure with a three-dimensional hole-conducting network, which contributes to its favorable TE properties.176 In contrast, Cu2ZnSn(S,Se)4 and Cu2CdSnSe4 possess distorted diamond-like structures, where the multielement nature suppresses phonon transport and lowers κ.177,178 CuInTe2 and CuGaTe2 undergo structural phase transitions associated with cation disorder, which also affects the transport behavior.179–181 Similarly, Cu2ZnSnS4 exhibits phase transitions at 530 K and 1149 K.182,183 Overall, further enhancement of TE performance in multicomponent copper chalcogenides relies on band engineering, doping/alloying, compositing, and entropy engineering, enabling ZT > 1.184,185
In addition, copper-based systems with alkali, alkaline-earth, and rare-earth metals are also of considerable interest for TE applications.186–191 Hodges et al. synthesized a series of bulk copper chalcogenides, ACu4.2TeS2 (A = K, Rb, Cs), in which the mixed-anion motif expands the Cu sublattice relative to ACu4S3, accommodating additional interstitial Cu atoms, resulting in high n (2 × 1021 cm−3).192 Jafarzadeh et al. synthesized polycrystalline BaCu6−xSe1−yTe6+y samples in which the Se/Te ratio adjusts Cu deficiency and thus n, with BaCu5.74(Se0.46Te6.54) achieving an electrical conductivity of 685 S cm−1. While increased Se/Te mass fluctuation increased material stability and reduced κl, the simultaneously decreased Seebeck coefficient prevented a net improvement in ZT.193 Early reports demonstrated that substitution of S by Te reduces the band gap in LaCuSTe and SmCuSTe.194 This finding is corroborated by recent studies on Cu3RETe3 (RE = Er, Ho, Tb), where the fully Te-based compositions yield promising TE performance, with peak ZT values ranging from 0.7 to 0.9 at 873 K.195
![]() | ||
| Fig. 2 Crystal structures of (a) paracostibite and (b) costibite CoSbS: Co, Sb, and S atoms are represented in blue, gold, and yellow, respectively. Paracostibite crystallizes in the space group Pbca, while costibite adopts the space group Pmn21. Reproduced with permission from AIP Publishing ©2020 (License CC BY).196 | ||
For the TE application, the paracostibite CoSbS structure is more advantageous; in this structure the Sb and S atoms form a layered arrangement within which the Co atoms occupy interstitial sites. This results in an n-type semiconductor with an indirect bandgap of 0.5 eV.197–200 In contrast, costibite CoSbS is a semimetal with the conduction band minimum overlapping with the valence band maximum.198,201 The key structural difference between the two lies in the degree of atomic distortion and bonding interactions; paracostibite exhibits a more distorted coordination environment that contributes to its semiconducting behavior, whereas the relatively symmetric atomic arrangement facilitates the metallic conduction of costibite. The more distorted paracostibite structure comprises 24 atoms and 8 formula units per unit cell (a = 5.843 Å, b = 5.956 Å, and c = 11.667 Å), whereas costibite has a unit cell consisting of 6 atoms and 2 formula units (a = 4.873 Å, b = 5.852 Å, and c = 3.608 Å).196,198,200,201 Paracostibite is the high-temperature polymorph of CoSbS, and it is stable up to 1073 K.202
Paracostibite CoSbS was studied as a TE material for the first time in 2011 by Carlini et al.203 Since then, it has gained continuous research interest owing to its flexibility for chemical substitutions, no intrinsic doping limit, high band degeneracy near the Fermi level, a high Seebeck coefficient (200 µV K−1), and a large power factor (2700 µW m−1 K−2).196 Electronic band structure calculations revealed that CoSbS has several low-lying conduction band minima and multiple electronic pockets, which are favorable for good TE performance. Interestingly, a gigantic low-temperature (40 K) S of 2.5 mV K−1 was reported for a single crystal, being notably higher than the values measured for polycrystalline samples.203–205
Paracostibite CoSbS is readily fabricated using various synthesis techniques: ball milling,206–208 a vapor phase technique,203 solid-state synthesis followed by SPS for pellet densification,201,209–211 and a heat-up method followed by two-step hot-pressing for nanoparticle synthesis.196 Chemically synthesized nanoparticles offer advantages such as enhanced phonon scattering to reduce κ.196,212 The thermal conductivity of CoSbS is unfortunately high (8 W m−1 K−1) compared to state-of-the-art TE materials due to stiff bonds and light anions.
Density functional theory (DFT) calculations predicted that a high ZT value up to 1.0 at 600 K could be achieved for optimized CoSbS if the high κ can be significantly reduced.203,213 To reduce κ and to optimize the TE performance, various cation substitutions with Ni, Te, Se, Zn, Mn, Fe, Cu, Pd, and Cr have been investigated. Partial Ni-for-Co substitution decreased the κ to 3.5 W m−1 K−1 at 600 K,206,214 while the Cu-for-Co and Cr-for-Co substitutions reduce κ to 2–3 W m−1 K−1 (at 700–900 K), respectively.207,215 Selenium is an efficient anionic dopant as it decreases the formation energy of point defects, suppressing κ to 3 W m−1 K−1 at 600–900 K through mass and strain fluctuations; this leads to enhanced Umklapp scattering.211,216–218
Similarly, Te-for-Sb substitution could reduce κ to 2.23 W m−1 K−1 due to point defect scattering, Umklapp scattering, mass fluctuation (mass difference in Sb and Te), and strain fluctuation (size difference in Sb and Te).209,210,219 On the other hand, simultaneous introduction of Se and Ni could not decrease κ below 2.73 W m−1 K−1.201,208
Recently, our group succeeded in massively suppressing κ down to 1.32 W m−1 K−1 for the Co(Sb0.9Ge0.1)(S0.95Se0.05) composition due to softening of the low energy acoustic phonon modes and strong scattering of phonons from Ge and Se born defects.220 DFT calculations on Ge-doped CoSbS involved replacing 12.5% of the Sb atoms with Ge. This substitution leads to a slight (2%) reduction in unit cell volume and a modest (1–4%) decrease in bond lengths. Importantly, Ge acts as an acceptor dopant by introducing metallic behavior: it shifts the Fermi level into the valence band, as shown in Fig. 3, thereby enhancing electrical conductivity.
![]() | ||
| Fig. 3 Electronic band structures of CoSbS (black) and Ge-substituted Co(Sb0.875Ge0.125)S (blue). Reproduced with permission from The Royal Society of Chemistry ©2024 (License CC BY 3.0).220 | ||
Phonon dispersion analysis confirmed that both non-substituted and Ge-for-Sb substituted CoSbS are dynamically stable. However, Ge substitution alters the vibrational properties by lifting the degeneracy of the transverse acoustic (TA) modes and softening the TA1 branch along the Γ–X–S–U–R path. A nearly dispersionless phonon mode also appears at around 33 meV. Projected phonon density of states (PDOS) reveals that Ge predominantly affects low-energy phonons, reducing their group velocity and enhancing phonon scattering, as shown in Fig. 4.
![]() | ||
| Fig. 4 DFT-calculated phonon dispersion curves for (a) CoSbS and (b) Ge-doped CoSbS (Co(Sb0.875Ge0.125)S). Black lines represent the acoustic phonon modes, while red lines indicate the optical phonon modes. Reproduced with permission from The Royal Society of Chemistry ©2024 (License CC BY 3.0).220 | ||
These computational insights explain the experimental observation that upon the Ge-for-Sb substitution, lattice thermal conductivity drops dramatically from 5.17 W m−1 K−1 to 1.32 W m−1 K−1, approaching the theoretical minimum. Fig. 5a shows the thermal conductivity data for Ge-for-Sb and Se-for-S co-substituted CoSbS. The reduction seen for κ is attributed to increased phonon–defect interactions and phonon dispersion analysis. The phonon mean free path (MFP) shortens considerably with Ge incorporation, confirming enhanced phonon scattering. Fig. 5b shows the thermal conductivity plot for different compositions of CoSbS-based materials.
![]() | ||
| Fig. 5 (a) Temperature-dependent thermal conductivity of Co(Sb,Ge)(S,Se) samples. (b) Comparison of thermal conductivity at 400 K for CoSbS-based materials studied in this work and in previous reports. Reproduced with permission from The Royal Society of Chemistry ©2024 (License CC BY 3.0).220 | ||
CoSbS is widely reported as an n-type semiconductor in the 200–900 K temperature range.201,205,206,209,210,214 However, it has also been reported as a p-type material with a gigantic PF below RT.204 Upon substitutions with Ni-for-Co, Te-for-Sb, and Se-for-S, n-type conduction is observed.201,206,209,214,221,222 In contrast, a transition in conduction from p-type to n-type is observed at 475 K when substituted with Cu-for-Co and Se-for-S.207,216,219,223 These findings suggest the bipolar nature of CoSbS, which needs systematic study for deeper understanding.
Recently we confirmed the bipolarity of CoSbS by slightly changing the vapor pressure of sulfur during synthesis; this was shown to change the type of conductivity from p- to n-type.224 Fig. 6 illustrates how sulfur stoichiometry influences the conduction type in CoSbS. When the material is sulfur-deficient, it exhibits p-type behavior, indicated by positive Seebeck and Hall coefficients. In contrast, sulfur-rich compositions result in n-type conduction, where both the Seebeck and Hall coefficients are negative.
![]() | ||
| Fig. 6 The p-type to n-type conductivity transition as a function of sulfur content in CoSbS. The Seebeck coefficient (blue, left axis) and the Hall coefficient (purple, right axis) of CoSbS plotted against nominal sulfur stoichiometry. Reproduced with permission from AIP Publishing ©2022 (License CC BY 4.0).224 | ||
Predictions based on Wiedemann–Franz law considerations and DFT calculations have suggested that ZT values higher than unity at high temperatures should be possible for properly optimized CoSbS.203,205 Early experimental attempts for partially Ni-for-Co substituted samples reached ZT values up to 0.50 at 873 K when using ball milling to synthesize nanostructured samples.206,214 By means of the anionic Se substitutions only, the ZT value could not be increased beyond 0.35.211,216,225
However, by combining the Se-for-S and Ni-for-Co substitutions, a record-high ZT value of 0.58 was realized at 900 K.215,223 The reason behind this is that the Ni-for-Co substitution increases n and m*, while the larger atomic radius of Se as compared to S softens the bonding, hence enhancing the phonon scattering. In line with this, Te-for-S substitution has been found to remarkably reduce κ owing to point defect scattering and also to increase n, leading to a jump in the PF to the highest reported value of 2.7 W m−1 K−1 and an enhanced ZT value to 0.47 at 725 K.208,210 Co-substitutions with Te and Ni or Te and Se have been found effective as well, increasing the m* and reaching the highest reported ZT value of 0.65 at 873 K208 and 0.62 at 730 K.201
The selenium- and tellurium counterparts of CoSbS are somewhat different from the sulfur-based compound. Early reports suggest that CoSbSe exhibits an orthorhombic structure with the space group Pnn2. However, Chmielowski et al.201 later reported that the exact space group is Pnm21, with a band gap of 0.4 eV. The phase shows low resistivity (<10−3 Ω cm),199,201,226 high carrier density (1021 cm−3), and metallic behavior with a small negative Seebeck coefficient of −14 µV K−1 at RT. Despite the high carrier density, the κ of CoSbSe is mainly governed by lattice contributions.227 The corresponding characteristics of CoSbTe are a pseudo-orthorhombic structure (Pnn2) and metallic transport behavior with a carrier density of the order of 1022 cm−3. The RT Seebeck coefficient is negative (−7.5 µV K−1) and the thermal conductivity behavior resembles that of metallic alloys, with a significant electronic contribution.227,228 In Table 1, we summarize the early progress in optimizing the TE characteristics of various CoSb(S,Se,Te) compositions.
| Year | Composition | Synthesis | κ | ZT | T | Ref. |
|---|---|---|---|---|---|---|
| The units for κ and T are W m−1 K−1 and K, respectively. Abbreviations for synthesis and processing: SSS = solid-state synthesis; BM = ball milling; SPS = spark plasma sintering; HP = hot pressing. | ||||||
| 2024 | Co(Sb0.8S0.8Se0.4) | SSS + SPS | 2.6 | 0.3 | 823 | 217 |
| 2024 | Co(Sb0.9Ge0.1)(S0.95Se0.05) | SSS | 1.32 | 0.1 | 400 | 220 |
| 2023 | Co(Sb0.9Se0.1)S | SSS + SPS | 3.02 | 0.22 | 823 | 218 |
| 2022 | (Ni0.10Co0.90)Sb(S0.85Se0.15) | SSS | 3.0 | 0.52 | 876 | 215 |
| 2022 | (Co0.94Ni0.06)Sb(S0.96Te0.04) | BM | 2.0 | 0.65 | 873 | 208 |
| 2020 | CoSbS nanoparticles | HP | 2.0 | 0.05 | 661 | 196 |
| 2020 | Co(Sb0.96Te0.04)(S[Te]0.28) | SSS + SPS | 2.73 | 0.45 | 773 | 210 |
| 2019 | (Co0.06Cu0.04)SbS | BM | 2.0 | 0.12 | 723 | 207 |
| 2018 | (Co0.93Ni0.07)Sb(S0.93Se0.07) | SSS + SPS | 2.4 | 0.58 | 900 | 223 |
| 2018 | Co(Sb0.93Te0.07)S | SSS + SPS | 2.23 | 0.43 | 900 | 219 |
| 2017 | (Co0.95Zn0.05)Sb(S0.85Se0.15) | SSS + SPS | 3.8 | 0.34 | 875 | 225 |
| 2017 | Co(Sb0.96Te0.04) (S0.75Se0.25) | SSS + SPS | 3.48 | 0.62 | 730 | 201 |
| 2017 | CoSb(S0.85Se0.15) | SSS + SPS | 4.72 | 0.35 | 923 | 211 |
| 2017 | CoSb(S0.99Se0.01) | SSS + SPS | 3.0 | 0.26 | 900 | 216 |
| 2016 | Co(Sb0.96Te0.04)S | SSS + SPS | 4.0 | 0.47 | 725 | 209 |
| 2015 | (Co0.92Ni0.08)SbS | BM | 3.6 | 0.5 | 873 | 206 |
| 2013 | (Co0.94Ni0.06)SbS | SSS | 3.5 | 0.35 | 773 | 214 |
Binary silver chalcogenides (Ag2S, Ag2Se, Ag2Te) with n-type conductivity crystallize in multiple polymorphic forms and undergo several phase transitions upon heating.233 At room temperature, α-Ag2S and β-Ag2Te adopt a monoclinic structure (space group P21/c), while β-Ag2Se exhibits an orthorhombic structure (space group P212121). Upon heating, all three compounds transition into distinct cubic crystal structures, while Ag2S and Ag2Te undergo at least two phase transitions. The Ag2S semiconductor is distinguished by its remarkable ductility, attributed to the rearrangement of atoms within the unit cell without bond breaking under pressure, as confirmed by DFT simulations.237 Ag2S materials can withstand up to 50% compressive strain, 20% bending strain, and 4.2% tensile strain.
Meanwhile, the monoclinic Ag2Te and orthorhombic Ag2Se phases appear brittle at RT, in contrast to their more ductile sulfur-containing analogs. The analysis of TE transport in Ag2S reveals that this compound, with a relatively wide band gap (Eg ≈ 1.0 eV),238 exhibits a near-zero ZT at RT, primarily due to its low electrical conductivity (σ ≈ 10 S m−1). This low σ results from the low carrier concentration (n ≈ 1014–1015 cm−3).233 Moreover, the high volatility of sulfur can alter the Ag-to-S stoichiometry, leading to a non-optimal concentration of interstitial Ag defects, which, in turn, significantly impacts the carrier concentration even at low concentrations.239–242 In contrast, β-Ag2Se and β-Ag2Te, with Eg values of 0.05–0.20 eV,243,244 exhibit significantly higher TE performance, reaching ZT values of 0.60–0.96 at RT, primarily due to the high carrier mobilities (μ).245,246 Notably, Ag2(S,Se,Te) compounds exhibit low lattice thermal conductivity (κl < 0.6 W m−1 K−1), attributed to structural disorder and liquid-like motion of Ag ions.234 Thus, the development of silver-based chalcogenide materials that combine the high ductility of Ag2S with the excellent TE properties of Ag2Se and Ag2Te represents a promising approach for enhancing the performance and durability of flexible TEGs. This can be achieved through anion substitution in Ag2X via doping, alloying, and composite formation, all of which have been shown to enhance ZT by optimizing the n and μ through precise band structure engineering and reducing the formation energy of Ag defects. Also, anion substitution can reduce κl by forming multiscale phonon scattering centers. Table 2 summarizes the TE properties of the Ag-based materials with mixed chalcogens reported since 2020.
| Year | Composition | ZT | T | σ | S | PF | κ | κl | Synthesis and processing | Ref. |
|---|---|---|---|---|---|---|---|---|---|---|
| The units for T, σ, S, PF, and κ and κl are K, S cm−1, μV K−1, μW cm−1 K−2, and W m−1 K−1, respectively. Abbreviations for synthesis and processing: MeltAnneal = melting–annealing; Quench = quenching; MechAl = mechanical alloying (ball milling); ZonMelt = zone melting; WetMechAl = wet-mechanical alloying; HydS = hydrazine solution-based synthesis; SolvoS = solvothermal synthesis; Rol = rolling; SPS = spark plasma sintering; Cut = cutting; HP = hot pressing. | ||||||||||
| 2025 | (Ag1.99Au0.01)(S0.5Te0.5) | 0.95 | 550 | 221 | −170 | 6.4 | 0.37 | 0.19 | MeltAnneal, Cut | 255 |
| 2025 | Ag2Te–Ag2S | 0.42 | 373 | 700 | −118 | 9.7 | 0.86 | 0.40 | SolvoS, SPS | 272 |
| 2024 | (Ag1.9Sn0.1)(S0.7Se0.3) | 0.42 | 343 | 500 | −100 | 5.0 | 0.41 | 0.18 | MeltAnneal, Quench, Cut | 253 |
| 2024 | Ag1.9992(S0.7Te0.3) | 0.85 | 550 | 329 | −150 | 7.4 | 0.48 | — | MeltAnneal, Cut | 240 |
| 2024 | Ag2(S0.7Te0.3) | 0.59 | 723 | 449 | −118 | 6.3 | 0.77 | 0.24 | MeltAnneal, Quench, SPS | 258 |
| 2024 | Ag2Se–1.0%Bi2S3 | 0.96 | 370 | 1040 | −160 | 26.6 | 1.03 | 0.45 | HydS, HP | 262 |
| 2024 | Ag2(Te0.5S0.5) | 1.34 | 523 | 592 | −114 | 7.7 | 0.30 | — | MeltAnneal, Quench, Cut | 265 |
| 2024 | Ag2(Te0.9S0.1) | 0.60 | 600 | 402 | −130 | 6.8 | 0.68 | 0.26 | MeltAnneal, Cut | 267 |
| 2024 | Ag1.98(S0.34Se0.33Te0.33) | 0.85 | 500 | 236 | −180 | 7.6 | 0.45 | 0.25 | MeltAnneal, Cut | 239 |
| 2024 | Ag1.98(S1/3Se1/3Te1/3) | 0.80 | 460 | 168 | −190 | 6.1 | 0.35 | 0.20 | MeltAnneal, Cut, Rol | 276 |
| 2023 | (Ag1.97V0.03)(S0.55Se0.45) | 0.71 | 350 | 200 | −239 | 11.4 | 0.56 | 0.43 | MeltAnneal, SPS | 250 |
| 2023 | Ag2(S0.55Te0.457) | 0.39 | 300 | 431 | −120 | 6.2 | 0.48 | 0.25 | MeltAnneal, Cut | 256 |
| 2023 | Ag2(Te0.6S0.4) | 0.80 | 623 | 205 | −186 | 7.1 | 0.55 | — | MeltAnneal, Cut | 269 |
| 2023 | Ag2(S0.5Se0.5)–0.5%Ag2Te | 0.43 | 323 | 325 | −135 | 5.9 | 0.45 | 0.28 | MeltAnneal, MechAl, SPS | 278 |
| 2022 | Ag2(S0.7Te0.3) | 0.80 | 600 | 185 | −190 | 6.7 | 0.50 | 0.27 | MeltAnneal, Cut | 257 |
| 2022 | Ag2(Te0.6S0.4) | 0.80 | 573 | 357 | −140 | 7.0 | 0.50 | 0.18 | MeltAnneal, Quench, Cut | 268 |
| 2022 | Ag2((Se0.6Te0.4)0.6S0.4)–0.03%Se | 0.68 | 470 | 248 | −166 | 6.8 | 0.47 | 0.30 | MeltAnneal, Cut | 277 |
| 2021 | Ag20(S7Te3) | 0.80 | 600 | 242 | −169 | 7.0 | 0.51 | 0.28 | MeltAnneal, Cut | 287 |
| 2021 | Ag2(S0.5Se0.5) | 0.27 | 330 | 294 | −118 | 4.1 | 0.50 | 0.30 | MeltAnneal, MechAl, SPS | 252 |
| 2021 | Ag2(Se0.9Te0.1) | 0.80 | 390 | 1100 | −132 | 19.2 | 0.94 | 0.18 | WetMechAl, SPS | 263 |
| 2021 | Ag1.9(S1/3Se1/3Te1/3) | 0.55 | 423 | 188 | −160 | 4.8 | 0.37 | 0.24 | MechAl, SPS | 242 |
| 2020 | Ag2SeS0.01 | 0.80 | 300 | 1224 | −140 | 24.0 | 0.90 | 0.60 | MeltAnneal, Cut | 241 |
| 2020 | Ag2(Se0.9S0.1) | 0.80 | 358 | 1190 | −137 | 22.3 | 1.00 | 0.53 | MeltAnneal, SPS | 260 |
| 2020 | Ag2(Te0.6S0.4) | 0.70 | 573 | 567 | −118 | 7.9 | 0.65 | — | MeltAnneal, Cut, HP | 270 |
The first-principles calculations by Nam et al.247 revealed that the single substitution of Ag with transition metals has little impact on the TE performance of binary Ag2S. In addition, Wuliji et al.248 examined 17 aliovalent dopants for α-Ag2S, which shows their extremely low doping limit (<0.002%), highlighting the crucial role of lattice manipulation at the anion site. In contrast, isovalent Se and Te dopants effectively modulate both the electronic band structure and the formation energies of Ag interstitials, enabling carrier concentration adjustment by 2–3 orders of magnitude.
This was further supported by Nam et al.249 and Sato et al.,250 whose first-principles calculations of the electronic transport properties in the crystal structure of Ag2(S,Se) (Fig. 7a) revealed that Se alloying reduces the band gap (Fig. 7b–d), lowers the electron effective mass (m*), and decreases the formation energy of Ag interstitials.
![]() | ||
| Fig. 7 (a) Crystal structure of Ag2(S,Se) solid solutions with random substitution of S and Se (yellow/green atoms) and high-symmetry points in the first Brillouin zone. (b)–(d) Projected electronic band structures and partial density-of-states of Ag2(S,Se) for different Se concentrations. (e)–(g) Phonon band structures, phonon DOS, and Bose–Einstein distribution functions (orange curves) of the non-alloyed and alloyed Se systems. Reproduced with permission from AIP Publishing ©2022.249 | ||
Se alloying maintains a reasonably high Seebeck coefficient (S), increases n, and thereby enhances the power factor (PF). The expected increase in electronic thermal conductivity (κe) is compensated by a decrease in κl due to the new scattering centers for phonons (Fig. 7e–g). Indeed, Hirata et al.251 showed that phonons in polycrystalline Se-alloyed Ag2S, synthesized via a melting and annealing method, can reach the strongest scattering limit of the Cahill model, resulting in a remarkably low κl of 0.4 W m−1 K−1 at RT. Notably, an additional electron density distribution appears between Ag2 sites in Ag2(S0.5Se0.5), which suggests the presence of interstitial Ag sites. This observation implies that the number of interstitial Ag defects increases with higher Se content. Consequently, such an increase results in a higher n, an increase in κe, and a reduction in κl.
Li et al.252 demonstrated that Se alloying in polycrystalline Ag2S optimizes the n and enhances TE performance. Furthermore, the co-alloying of Cu and Se in Ag2S not only alters the phase composition, shifting the phase transition to 370 K, but also introduces a secondary Cu-based phase and enhances the liquid-like behavior, leading to a reduction in κl. As a result, Cu–Se co-alloying effectively modulates phase transition in Ag2S, resulting in a ZT of 0.20 at 450 K. A peak ZT of 0.27 at 330 K was observed for the Ag2(S0.5Se0.5).
Moreover, Se-alloying in Ag2S can also benefit the extension of the cation substitution limit. For example, as shown by Wu et al.,253 dual Se and Sn alloying in polycrystalline (Ag1.9Sn0.1)(S0.7Se0.3) composition with a monoclinic Ag2S crystal structure reduced Eg and shifted the Fermi level (EF) into the conduction band, facilitating n-type charge transport. The computed crystal orbital Hamilton population (–COHP) for Ag–S bonding in Ag2S and Ag2(S0.7Se0.3) indicates that Se alloying slightly decreases this value from 1.566 eV in Ag2S to 1.482 eV in the Se-alloyed compound. The (Ag1.9Sn0.1)(S0.7Se0.3) composition with a high n ≈ 7.29 × 1018 cm−3 and μ ≈ 321 cm2 V−1 s−1 exhibited a PF of 5 µW cm−1 K−2 and a ZT of 0.42 at 343 K, while maintaining exceptional ductility exceeding 90% due to the formation of a biphasic structure that induces multiscale phonon scattering centers, resulting in a reduced κl of 0.18 W m−1 K−1. Adjustment of the S/Se ratio and V-doping in polycrystalline (Ag1.97V0.03)(S0.55Se0.45) composition with an α-Ag2S crystal structure resulted in a ZT of 0.71 at 350 K.250
High-entropy materials possess a distorted lattice and diversified chemical bonding, which may significantly extend the element doping/alloying limit.254 Indeed, the Ag2(S0.5Te0.5) solid solution with a highly distorted crystal structure due to the large atomic size and electronegativity mismatches between S and Te extends the cation substitution limit, i.e., effectively accommodates the foreign Au atoms at cation sites up to 0.07.255 The (Ag,Au)2(S0.5Te0.5) were composed of a mixture of cubic and amorphous phases. Substituting Ag with Au in (Ag,Au)2(S0.5Te0.5) reduced n from 8.7 × 1018 cm−3 at x = 0 to 7.7 × 1017 cm−3 at x = 0.07. Therefore, the reduced σ and κ, combined with increased S and PF, resulted in a ZT of 0.95 at 550 K in (Ag1.99Au0.01)(S0.5Te0.5). Notably, Au alloying had a negligible effect on the mechanical properties of (Ag1.99Au0.01)(S0.5Te0.5), preserving a bending strain of 15%, a compressive strain of 20%, and a Vickers hardness of 20 kgf mm−2.
Li et al.256 synthesized an Ag2(S0.55Te0.45) solid solution with a slight excess of Te to reduce the n. This modification enhanced the electrical transport properties and reduced the κ, resulting in a high PF of 6.2 μW cm−1 K−2 and a ZT of 0.39 at 300 K. Notably, Ag2(S0.55Te0.457) also exhibited a high average ZT among ductile silver chalcogenides. Moreover, Te-alloying in Ag2(S0.7Te0.3) solid solution enhanced μ, leading to a ZT of 0.80 at 600 K (Fig. 8a and b).257 An adjustment of Ag enabled further optimization of n in Ag1.9992(S0.7Te0.3), yielding a ZT of 0.85 at 550 K (Fig. 8c–f).240 The Ag2(S,Te) solid solutions exhibited a body-centered cubic structure at 300–600 K and demonstrated superior ductility compared to Ag2S.
![]() | ||
| Fig. 8 (a) Te content dependence of the Hall carrier concentration (nH) and mobility (μH) of Ag2(S,Te) at 300 K. ZT values of (b) Ag2(S,Te). Reproduced with permission from Elsevier ©2022.257 (c) Ag2+x(S0.7Te0.3) as a function of temperature, (d) Ag1.9992(S0.7Te0.3) and reported Ag2(S,Te), and (e) Ag1.9992(S0.7Te0.3) at 300–550 K, and (f) Ag1.9992(S0.7Te0.3) as a function of carrier concentration at 300 K and 550 K. Reproduced with permission from Elsevier ©2024.240 | ||
Efforts to stabilize the cubic phase of Ag2S at low temperatures have been undertaken by Zhong et al.258 The study demonstrated that Ag2(S,Te) solid solutions can maintain a stable cubic phase at RT. These materials also exhibited excellent mechanical properties, and the Ag2(S0.8Te0.2) composition showed a bending yield strength of 46.52 MPa at 673 K, surpassing that of Ag2S, which declined from 80.15 MPa at 300 K to 12.66 MPa at 673 K. Notably, the Ag2(S0.7Te0.3) solid solution achieved a high ZT of 0.59 at 723 K. Jin et al.259 synthesized an Ag2(S0.7Te0.3) solid solution via a zone melting technique, confirming that S and Te atoms randomly occupied the anionic sites. This composition exhibited a high μ ≈ 410 cm2 V−1 s−1 and a low κl ≈ 0.2 W m−1 K−1, resulting in a peak ZT value of 0.30 at 600 K.
Orthorhombic Ag2Se is brittle, whereas monoclinic Ag2S exhibits notable ductility. Liang et al.260 demonstrated that the compositional transition boundary between monoclinic and orthorhombic crystal structures in polycrystalline Ag2(Se1−xSx) occurs at x = 0.3. Substituting Se with S modifies the valence band maximum rather than the conduction band minimum, leading to an unaffected effective mass (m*). The relatively small m* (0.25me) in Ag2(S, Se) contributes to their high μ, thereby yielding a high ZT of 0.8 at 358 K for Ag2(Se0.9S0.1) (Fig. 9a). Additionally, orthorhombic Ag2(Se0.9S0.1) exhibits limited bending deformation (1.5%), whereas monoclinic Ag2(Se0.6S0.4) demonstrates significantly higher bending deformation (∼10%) without cracking. Thus, S-alloying in Ag2Se enables the controlled formation of either an orthorhombic or monoclinic structure, offering an effective strategy for tuning its mechanical and TE properties.233,235 A minute anion excess in polycrystalline Ag2SeS0.01 can also stabilize the orthorhombic structure, preventing the formation of a metastable phase in the temperature range of 300–400 K, just before the transition to the cubic phase at ∼407 K.241,261 As a result, the high n (3.5 × 1018 cm−3) and μ (2030 cm2 V−1 s−1) contribute to a high PF of 2.4 mW m−1 K−2 at 300 K. The relatively low κ leads to a ZT of 0.8 at 300 K in the Ag2SeS0.01 composition (Fig. 9b). However, during the phase transition, the orthorhombic structure becomes unstable, leading to the re-emergence of a metastable phase. Although additional cation doping with Cu further enhances the PF (∼2.6 mW m−1 K−2), the resulting increase in κe prevents improvement in ZT.
![]() | ||
| Fig. 9 Temperature-dependent ZT of (a) Ag2(Se,S). Reproduced with permission from the American Association for the Advancement of Science ©2020.260 (b) (Ag,Cu)2(SeS0.01). Reproduced with permission from American Chemical Society ©2020.241 (c) Ag2Se–x wt% Bi2S3. Reproduced with permission from American Chemical Society ©2024.262 (d) Ag2(Se,Te). Reproduced with permission from American Chemical Society ©2021.263 | ||
The Ag2Se/Bi2S3 composites, synthesized in an aqueous solution with hydrazine as a reducing agent, exhibited an improved S and a slightly altered σ due to a decrease in n and an increase in μ, resulting from the partial substitution of Se with S.262 The PF of the polycrystalline Ag2Se–1.0 wt% Bi2S3 composite was enhanced to 2.66 mW m−1 K−2 at 370 K. Furthermore, with a decreased κ due to the effective phonon scattering both at point defects and interphase boundaries, ZT values of the Ag2Se–1.0 wt% Bi2S3 composites were improved up to 0.96 at 370 K (Fig. 9c).
DFT calculations and single parabolic band (SPB) modeling of Te-doped Ag2Se indicate that the reduced contribution of Se 4p orbitals to the total density of states decreases the carrier effective mass (m*) with increasing Te content. This reduction, along with dense phonon dispersion, is expected to enhance the theoretical ZT.263 The decrease in m* reduces the S but increases the σ. However, the deformation potential (Edef) also rises with increasing Te content, suppressing μ. Despite this, the reduced m* still contributes to an enhancement of PF. Thus, the theoretical peak ZT exceeds 1.0 in the Ag2(Se0.9Te0.1) composition at n = 7.0 × 1017 cm−3.
Moreover, SPB modeling predicts that theoretical ZT values above 1.1 at 300 K are achievable for Ag2(Se0.7Te0.3) composition at n = 4.0 × 1017 cm−3. Ternary Ag2(Se,Te) solid solutions, synthesized via wet-mechanical alloying and SPS, also exhibited high Te solubility (∼50%).263 The as-sintered ternary Ag2(Se,Te) solid solutions demonstrated complex microstructures, including dislocations, nanograins, grain boundaries, TeSe substitution, lattice distortions, and localized strain, which contribute to strong phonon scattering and result in ultralow κl values of 0.21–0.31 W m−1 K−1 at 300 K. The ZT of 0.65–0.80 was observed in the Ag2(Se0.9Te0.1) solid solution at 300–390 K (Fig. 9d).
The modulation of the anion ratio in Ag-based materials enhances their thermoelectric and mechanical properties. However, it also introduces phase complexity and multiple structural transitions, necessitating precise control to maintain stable TE parameters upon heating. Recently, Li et al.264 demonstrated that the Ag2(Te,S) solid solutions consist of both amorphous and crystalline phases at 20–420 K. The Ag2Te-based phase retains its amorphous state within the 20–300 K range, whereas the amorphous Ag2S-based phase begins to crystallize below 183 K. Notably, in Ag2(S0.4Te0.6), the presence of an amorphous Ag2S-based phase at 183–420 K ensures excellent ductility and low κl, yielding a ZT of 0.17 at 300 K. During the heat treatment of Ag2(Te,S) solid solutions,265 three crystal phases can emerge: a low-temperature crystalline phase, an amorphous phase, and a high-temperature cubic phase. The high-temperature cubic phase transforms into the low-temperature crystalline phase with the formation of a metastable amorphous phase, enhancing σ below the phase transition temperature. The Ag2(Te0.5S0.5) compound, annealed at 473 K, exhibited a ZT of 1.02 at 623 K, which increased to 1.34 at 523 K during the second heating cycle.
A variety of strengthening methods, including alloying, grain refinement, and secondary phase precipitation, have been employed to tune the mechanical properties of TE materials.266 In this context, understanding the relationship between phase composition and mechanical properties in Ag-based TE materials is crucial for elucidating the mechanisms underlying their exceptional ductility. The analysis of mechanical properties in Te-rich Ag2(Te,S) solid solutions by Wang et al.267 revealed a brittle–ductile transition. High ductility was observed in Ag2Te crystals with the monoclinic and cubic phases, previously considered brittle. The ductile Ag2(Te0.9S0.1) exhibited a ZT of 0.6 at 600 K. The removal of the monoclinic Ag2Te phase in the Ag2(Te,S) compositions increased compressive strain, exceeding 70%.
A high tensile elongation of 107% was achieved in Ag2(Te0.7S0.3) solid solutions.268 Moreover, the high mobility (μ ≈ 1000 cm2 V−1 s−1 at 300 K) results in a ZT value of 0.8 at 573 K for the annealed Ag2(Te0.6S0.4). The reversible brittle–ductile transition in Ag2(Te0.6S0.4) was demonstrated by tuning the phase composition through heat treatments.269 The Ag2(Te0.6S0.4), annealed at 973 K for 7 days, exhibits an amorphous/cubic structure and shows a compressive strain exceeding 80%, with an engineering stress of 600 MPa. In addition to its superior plasticity, Ag2(Te0.6S0.4) demonstrates excellent TE properties, attaining a high ZT of 0.8 at 673 K.
The flexibility and high TE performance of Ag2(Te0.6S0.4) were achieved by He et al.270 through amorphization. Amorphization in the S-doped Ag2(Te0.6S0.4) composition contributed to excellent mechanical properties, including a bending strain above 14% under a stress of 110 MPa, an elastic strain limit of 2% under 70 MPa, a compression strain of 25%, a tensile strain of 12.5%, and a Vickers hardness of 19.5 kgf mm−2. The coexistence of amorphous and crystalline phases contributed to exceptional flexibility and a low κ of 0.65 W m−1 K−1 at 573 K. A high μ (750 cm2 V−1 s−1) and n (8.6 × 1018 cm−3) led to a high PF of 7.9 μW cm−1 K−2 and a peak ZT of 0.7 at 573 K surpassing that of organic TE materials.271
The Ag2Te/Ag2S composites, synthesized via a low-cost solvothermal method and SPS, exhibited a tunable phase composition, optimizing both n and μ. In particular, Ag2(S0.8Te0.2), with a phase distribution of Ag2Te (95.06%) and Ag2S (4.94%), demonstrated ZT values of 0.42 at 373 K and 0.38 at 298 K.272 The enhanced ZT in these composites results from reduced κl, facilitated by phonon scattering at the branched vein-like microstructure, amorphous transition phase boundaries, and dislocations. The addition of Ag2S improved the plasticity of the Ag2S/Ag2Te composites, achieving a band strain of 2.5%. Also, n increased from 1014 cm−3 (in Ag2S) to 1018 cm−3 (in Ag2(S0.8Te0.2)), while μ reached 5436 cm2 V−1 s−1, 1.86 times higher than that of Ag2Te.
Doping at the anion site has been widely investigated for the synthesis of high-entropy chalcogenide compounds featuring metastable phases.273,274 The first-principles calculations by Biswas et al.275 predict an enhanced power factor in the pseudo-ternary Ag2(Se0.5Te0.25S0.25) phase, driven by spin-dependent band splitting known as the Rashba effect following S and Te substitution. Consequently, increased configurational entropy theoretically predicts a low κl (0.34 W m−1 K−1) and a high ZT of 2.1 at 400 K.
The efforts to stabilize the high-temperature cubic β-Ag2S phase through entropy engineering239 have successfully resulted in the formation of single-phase Ag2(S,Se,Te) compounds at RT, exhibiting superior mechanical properties, including ductility. Notably, the mixed-anion Ag1.98(S0.34Se0.33Te0.33) with slight Ag deficiency demonstrated an optimized σ, reduced κ, and a significantly enhanced ZT value of 0.85 at 500 K. These improvements in TE parameters are attributed to optimized configurational entropy, strong anharmonicity, and weakened ionic bonding, which collectively yield a high n ≈ 1018–1019 cm−3 and μ ≈ 826 cm2 V−1 s−1 at x = 0.33.
Chen et al.276 identified the compositional region of the cubic phase (Fig. 10) and proposed an anion-site alloyed Ag1.98(S1/3Se1/3Te1/3) with the highest configurational entropy (ΔS). This anion-site alloying not only enhanced the mechanical properties of Ag1.98(S1/3Se1/3Te1/3) but also improved its TE performance, achieving an ultralow κl of 0.2 W m−1 K−1 and a high ZT of 0.8 at 460 K. Notably, the material exhibits a tensile strain of 97%.
![]() | ||
| Fig. 10 (a) Schematic phase diagram of Ag2X (X = S, Se, Te). (b) Calculated configurational entropy ΔS considering the mixture of S/Se/Te at the anion site. (c) Phase diagram of maximum bending strain at room temperature of the Ag2S–Ag2Se–Ag2Te system. (d) Phase diagram of the room-temperature ZT value of the Ag2S–Ag2Se–Ag2Te system. Reproduced with permission from Wiley-VCH GmbH ©2024.276 | ||
Yu et al.242 synthesized an Ag1.9(Se1/3S1/3Te1/3) solid solution using high-energy ball milling followed by SPS. The anion-site alloying of Se, S, and Te significantly enhances phonon scattering, resulting in an ultralow κl of 0.2 W m−1 K−1. Additionally, sulfur volatilization during the SPS process induces porosity in the material, which further contributes to phonon scattering. Consequently, a peak ZT of 0.55 was observed at 423 K.
To extend the operating temperature range of Ag(S,Se,Te), Xie et al.277 alloyed Ag2(Se0.6Te0.4) with Ag2S to stabilize the cubic phase within the temperature range of 300–470 K. This alloying effect lowers the orthorhombic-to-cubic phase transition temperature, enabling Ag2(S0.4(Se0.6Te0.4)0.6) to maintain a stable body-centered cubic (bcc) structure. Also, Ag2(S0.4(Se0.6Te0.4)0.6) demonstrates remarkable mechanical properties, sustaining a bending strain exceeding 11% before fracture, in contrast to Ag2(Se0.6Te0.4), which exhibits a strain below 1%. Despite its stable cubic structure and excellent ductility, the Ag2(S0.4(Se0.6Te0.4)0.6) solid solution has an extremely high carrier concentration (n ≈ 3.1 × 1019 cm−3 at RT), leading to a low ZT of 0.06. To improve its TE performance, a small amount of excess Se was introduced into the Ag2(S0.4(Se0.6Te0.4)0.6), effectively reducing the carrier concentration by an order of magnitude. This modification decreased both σ and κe, while significantly enhancing S, resulting in a peak ZT of 0.68 at 470 K.
Wu et al.278 synthesized a polycrystalline Ag2(S0.5Se0.5) matrix with 0.5 mol% Ag2Te inclusions, achieving a ZT of 0.43 at 323 K and a high bending strain of 32.5%. The incorporation of Ag2Te optimized the n and μ, resulting in a high PF of 6 μW cm−1 K−2. Furthermore, Ag2Te induced the formation of Ag-poor amorphous phase boundaries, significantly enhancing the overall plasticity. These amorphous phase boundaries, combined with multiscale phonon scattering mechanisms, effectively lowered the κl to 0.28 W m−1 K−1 at 323 K. DFT calculations revealed that the intrinsic n-type conductivity of Ag2S and Ag2(S0.5Se0.5) is attributed to S and Se vacancies, as well as Ag interstitials. Additionally, Ag2Te results in amorphous phase boundaries that, along with multiscale phonon scattering centers, contribute to the reduction in κl.
The conventional method to synthesize Ag-based chalcogenides predominantly involves solid-state reactions, including melting, annealing, and subsequent quenching. The obtained ingots are typically cut into pieces or ground into powder, followed by SPS to achieve densified TE elements. However, this approach has limitations for large-scale production, whereas techniques such as hot-rolling could facilitate the fabrication of TE elements with a larger surface area. Recently, orthorhombic Ag2(Se,S) films with outstanding mechanical properties (Fig. 11) were successfully fabricated by hot-rolling bulk materials at 393 K.279 This approach promotes the formation of dense dislocations and grain refinement and increases the elastic strain while ensuring the recovery of TE properties even after 105 bending cycles. A maximal power output of 157 μW, corresponding to a specific power density of 5 μW cm−1 K−2, was achieved in the Ag2(Se0.9S0.1) solid solution at a temperature difference of 76 K. Another example is the fabrication of Ag2(Te,S) films via hot-rolling at 393 K, exhibiting an increased elastic strain of 1.8% in Ag2(Te0.8S0.2), representing a 200% enhancement compared to unalloyed Ag2Te.280
![]() | ||
| Fig. 11 (a) Softening temperature of Ag2(Se,S) alloys. (b) Dimensional evolution for Ag2(Se0.9S0.1) by multipass hot rolling. (c) Elasticity of Ag2(Se,S) alloys. (d) Stress versus strain for Ag2(Se0.9S0.1) bulks with and without hot-rolling. Reproduced with permission from American Chemical Society ©2024.279 | ||
The solid solutions between Ag2Te and Ag2Se also demonstrated remarkable mechanical and TE properties,281 a maximum compressive strain rate of 36%, a compressive yield strength of 180 MPa, and a Vickers hardness of 0.63 GPa in the Ag2(Te0.8Se0.2) bulks with a ZT of 0.3 at 300 K and 0.4 at 380 K. Wang et al.282 reported a unique room-temperature plastic deformation mechanism in inorganic semiconductors Ag2(Te,S), driven by sublattice amorphization coupled with Ag-ion diffusion.
The demonstrated high-performance Ag-based materials for flexible and wearable TEGs are promising to replace the costly and toxic Bi2Te3-based materials for thermal energy harvesting applications. A prototype of the TEG was successfully fabricated using ductile n-type Ag2(S0.55Se0.55) and over-stoichiometric p-type Cu2.075Se materials with a ZT value of 0.5 at 300 K. The monolithic TEG with negligible Ag- or Cu-ion migration exhibited a maximum power density of 0.68 mW cm−2 at ΔT = 30 K.283 The ultra-flexible TEG with 6 legs assembled with a Se-substituted Ag2S film on a nylon membrane demonstrated a power density of 0.66 mW cm−2 at a temperature difference of 28.8 K.284
Owing to its superb flexibility and high TE performance, Ag2(Te0.6S0.4) was utilized to develop a prototype of 3D wearable fabric, achieving a normalized power density of 0.4 μW m−1 K−2 at ΔT = 20 K.285 The effectiveness of orthorhombic Ag2Se and monoclinic Ag2S nanocomposites, synthesized via a solvothermal method on carbon fabric, was demonstrated in a 2-pair TE module, which showed an output voltage from 0.1 to 7.4 mV at ΔT = 3–8 K.286
The prototype TEG, consisting of 10 Ag20(S7Te3) strips, exhibited an open-circuit voltage of 69.2 mV and a maximum power output of 17.1 µW at a temperature difference of 70 K (Fig. 12a).287 The entropy-engineered Ag2(S1/3Se1/3Te1/3) can be readily rolled into flexible thin films. A six-leg in-plane device was fabricated using this material, achieving an output voltage of 13.6 mV, a maximum power of 12.8 μW, and a power density of 14.3 W m−2 at a temperature difference of 30 K (Fig. 12b).276 A flexible in-plane device with a single leg based on the amorphized Ag2(Te0.6S0.4) produced a power of 3.5 μW at a temperature difference of 50 K, demonstrating its great potential as a self-powered TE generator (Fig. 12c).270 The large-scale all-inorganic silver chalcogenide foils were produced using a rolling process by Liang et al.288 A high ZT of 0.47 at 310 K was achieved for Ag2(S0.45Se0.45Te0.1) composition. A proof-of-concept flexible TE generator based on these silver chalcogenide foils demonstrated an open-circuit voltage of 1.19 mV and an output power density of 1.8 mW m−2 with a temperature difference of 2.7 K (Fig. 12d).
![]() | ||
| Fig. 12 (a) Hetero-shaped Ag20(S7Te3)-based TEG and its output V–I, P–I characteristics. Reproduced with permission from Wiley-VCH GmbH ©2021.287 (b) Flexible in-plane Ag1.98(S1/3Se1/3Te1/3)-Pt-based TEG and power density vs. ΔT compared to typical flexible TEGs. Reproduced with permission from Wiley-VCH GmbH ©2024.276 (c) Ag2(Te0.6S0.4)-based TEG and power density vs. ΔT compared to typical flexible TEGs. Reproduced with permission from the American Association for the Advancement of Science ©2020.270 (d) Flexible Ag2(S0.6Se0.3Te0.1) foil-based TEG and its output V–I, P–I characteristics at a heat source temperature of 37 °C. Reproduced with permission from American Chemical Society ©2022.288 | ||
Beyond binary Ag-based compounds, multicomponent silver-bearing minerals deserve particular attention, including the argyrodites Ag8SnSe6,289 Ag8GeSe6,290 and Ag8SiSe6,291 as well as the pavonites AgBi3S5,292 Ag3Pb4Bi11Se22,293 and Ag1.75InSb5.75Se11.294 In particular, Ag-based superionic argyrodites have attracted considerable attention owing to their liquid-like ultralow κl supporting the PLEC concept. For instance, n-type cubic argyrodites Ag4M0.5S2Te (M = Sn, Ge) exhibit an ultralow κ of 0.30 W m−1 K−1, attributed to large lattice anharmonicity arising from triangular, rod-like, and cage-like cluster units alongside weakly bonded Ag+ ions. Optimization via generating Te deficiency enhances the PF to 2.61 μW cm−1 K−2, yielding a peak ZT of 0.74 at 823 K.295 Substituting Te for S increases the atomic mass and weakens the chemical bonds, leading to a substantial reduction in sound velocities. Through this lattice-softening approach, the thermal conductivity of Ag8Sn(S,Te)6 reaches 0.15–0.26 W m−1 K−1 at RT.296 Substituting Se for S in polycrystalline Ag8Sn(S,Se)6 shifts the orthorhombic-cubic phase transition, expanding the temperature range of the thermoelectrically favored cubic phase. A peak ZT of 0.80 at 773 K is achieved for Ag8Sn(S0.97Se0.03)6.297 Wu et al. synthesized AgBi3(Se,S)5 solid solutions. Se/S anion mixing in AgBi3(Se,S)5.08 simultaneously increases configurational entropy, enhancing phonon point defect scattering and reducing κl to 0.45 W m−1 K−1, and modifies the effective mass near the Fermi level, boosting electrical transport. As a result, polycrystalline AgBi3(Se0.9S0.1)5.08 achieves a peak ZT of 0.42 at 723 K.298
![]() | ||
| Fig. 13 Outlook of the effects of anion mixing on thermoelectric metal chalcogenides. Lead chalcogenides: (1) electronic band convergence. Adapted with permission from Springer Nature ©2011;46 (2) higher covalency and hardness, π-type TEGs. Adapted with permission from The Royal Society of Chemistry ©2020 (License CC BY 3.0);47 (3) multiscale defects (e.g., nano precipitates, strains). Adapted with permission from Elsevier ©2017.52 Bismuth chalcogenides: (1) alloying for higher anharmonicity. Adapted with permission from Wiley-VCH GmbH ©2022;58 (2) nanostructuring and interface engineering. Adapted with permission from American Chemical Society ©2016.62 Tin and indium chalcogenides: (1) tuning phase transitions and band structures; (2) softening phonon band structures. Adapted with permission from the American Association for the Advancement of Science ©2019;85 (3) band structure engineering. Adapted with permission from The Royal Society of Chemistry ©2021 (License CC BY 3.0).94 Copper chalcogenides: (1) nanoscale mosaicity. Adapted with permission from Wiley-VCH GmbH ©2015;111 (2) selectivity and control of phase composition. Adapted with permission from The Royal Society of Chemistry ©2017 (License CC BY 3.0);127 (3) semiconductor-to-metal transition. Adapted with permission from Wiley-VCH GmbH ©2023.152 Cobalt chalcogenides: (1) softening of the acoustic phonons. Adapted with permission from The Royal Society of Chemistry ©2024 (License CC BY 3.0);220 (2) tuning the p-type to n-type conductivity transition. Adapted with permission from AIP Publishing ©2022 (License CC BY).224 Silver chalcogenides: (1) tuning Ag interstitial defects. Adapted with permission from AIP Publishing ©2023 (License CC BY);251 (2) adjusting the crystal structure, phase transitions, and ductility. Adapted with permission from Elsevier ©2024;267 (3) prototyping TEGs. Adapted with permission from American Chemical Society ©2022.288 | ||
Anion mixing can also increase both anion and cation solubility limits, as was, for example, observed in entropy-driven Pb(S,Se,Te) solid solutions. This expands the accessible carrier concentration range, enables effective band gap tuning, and enhances phonon scattering through strong atomic mass contrast. The extent of S/Se/Te miscibility naturally varies across material families due to differences in the host crystal structure, site symmetry, and the tolerance of the cation sublattice to anion size mismatch, all of which govern the thermodynamic stability of the resulting solid solution.
Besides the anion mixing approach, another common chemical tuning approach applied to enhance the thermoelectric properties of metal chalcogenides is cation doping; these efforts have been briefly summarized as well in the present review for each metal chalcogenide material family. A distinct difference between the anion-mixing and cation doping/substitutions is that anion mixing at the S/Se/Te site is of the isovalent substitution type, while the most effective cation substitutions are often aliovalent. Nevertheless, isovalent anion substitution may also modify the cation-site defect equilibria, as was, for example, demonstrated in the case of In-based chalcogenides, where isovalent S-for-Se substitution was found to suppress the annihilation of donor defects and thereby enhance the electrical conductivity. In copper chalcogenides, anion mixing was found to stabilize Cu+ within its sublattice, thereby improving operational stability and fine-tuning the charge carrier transport; moreover, the resulting point defects provided additional phonon scattering centers, further reducing thermal conductivity and enhancing the thermoelectric performance.
In the present review, the anion-mixing effects were discussed in detail in separate sections both for the state-of-the-art Pb-, Bi-, Sn- and In-based systems and a number of currently strongly emerging Cu-, Co-, and Ag-based thermoelectrics. For lead chalcogenides, anion mixing was found to improve charge transport by promoting band convergence, enabling a better band gap, carrier concentration and effective mass control, and suppressing bipolar effects. Simultaneously, the resulting crystal disorder, point defects, lattice strain, and nanostructures enhance phonon scattering and reduce lattice thermal conductivity. Moreover, through tailored cation–anion bridges, decoupling of electron and phonon transport becomes possible. In bismuth chalcogenides, anion mixing enables control over the conductivity type, carrier concentration, and band gap, thereby suppressing bipolar effects. In tin and indium chalcogenides, anion mixing tunes carrier concentration and the valence band structure, increasing the density of states and improving electrical transport. At the same time, the induced defects and hierarchical microstructure strengthen phonon scattering, lower lattice thermal conductivity, and may stabilize thermoelectrically favorable phases.
The recent renewed interest in copper chalcogenides is largely due to their low lattice thermal conductivity which mainly arises from the liquid-like Cu+ sublattice; the trade-off is that this ion migration also causes phase instability. Here, impressively, the anion-mixing approach provided an effective way to not only stabilize the Cu sublattice by suppressing the Cu-vacancy formation but also to lower thermal conductivity and enhance the thermoelectric performance through strain-induced phonon scattering and an optimized electronic structure and carrier concentration. Cobalt chalcogenides, particularly the CoSbX compounds, offer a stable framework where anion mixing has been shown to induce phonon softening and strong scattering and accordingly reduce thermal conductivity. Interestingly, the electrical conductivity type of the CoSbX materials can be tuned by composition, enabling p–n module design, though challenges such as bipolar conduction still remain. In silver chalcogenides, anion mixing is an effective strategy for tuning both the thermoelectric and mechanical properties. In Ag2X compounds, low lattice thermal conductivity mainly arises from structural disorder and liquid-like Ag-ion motion, while anion substitutions further enhance phonon scattering through defects, strain, phase boundaries, and multiscale microstructures. Simultaneously, anion mixing tunes the band gap, carrier effective mass, and defect chemistry, thereby optimizing the carrier concentration, mobility, and power factor. Because cation substitution is often limited in these systems, anion mixing is valuable for expanding the compositional flexibility, stabilizing favorable phases, shifting structural transitions, and controlling crystallinity or partial amorphization. These effects make mixed-anion Ag-based chalcogenides promising for flexible devices, although further work is needed to improve the stability of practical prototypes.
Overall, metal chalcogenides provide an excellent platform in which anion mixing may function as a targeted strategy to simultaneously manipulate electron and phonon transport. Varying the S/Se/Te composition allows electronic structure tuning (control over the band gap, band alignment, band convergence, and carrier effective mass) needed to increase the power factor. In parallel, the created anion disorder may introduce mass fluctuation, local strain, and hierarchical microstructural heterogeneity, which is highly beneficial to enhance phonon scattering and suppress lattice thermal conductivity. Beyond the electronic and thermal transport tuning, anion mixing also regulates defect chemistry by modifying point defect formation energies, charge compensation, solubility limits, and secondary phase formation, while in many systems, it additionally drives phase stabilization or polymorphic transitions that can unlock the thermoelectrically favorable states. However, these benefits are not necessarily universal: excessive disorder may reduce carrier mobility, incomplete alloying may lead to phase separation, and structurally stabilized phases do not always yield superior thermoelectric performance. A deeper mechanistic understanding of how anion mixing couples electronic structure, defect chemistry, and phase stability should further accelerate the development of high-performance and especially tellurium-lean metal chalcogenide thermoelectrics so that mixed-anion metal chalcogenides could become key thermoelectric materials for sustainable energy technologies.
Footnote |
| † These authors contributed equally. |
| This journal is © The Royal Society of Chemistry 2026 |