
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39223124
52063
10.1038/s41467-024-52063-5
Article
Carrier-phonon decoupling in perovskite thermoelectrics via entropy engineering
http://orcid.org/0000-0002-5178-2212
Zheng Yunpeng 12
http://orcid.org/0000-0001-9086-7000
Zhang Qinghua 3
Shi Caijuan 4
http://orcid.org/0000-0002-4359-1731
Zhou Zhifang 1
Lu Yang 5
Han Jian 1
Chen Hetian 1
Ma Yunpeng 1
http://orcid.org/0000-0003-4892-1100
Zhang Yujun 4
http://orcid.org/0000-0002-6405-2390
Lin Changpeng 6
Xu Wei 47
http://orcid.org/0000-0003-0178-3994
Ma Weigang 5
http://orcid.org/0000-0001-7922-598X
Li Qian 1
Yang Yueyang 1
http://orcid.org/0000-0003-0125-1760
Wei Bin 18
Yang Bingbing 19
Zou Mingchu 1
Zhang Wenyu 1
Liu Chang 1
Dou Lvye 1
Yang Dongliang 4
http://orcid.org/0000-0001-6844-0879
Lan Jin-Le 10
http://orcid.org/0000-0002-7209-1947
Yi Di 1
http://orcid.org/0000-0002-7545-0773
Zhang Xing 5
http://orcid.org/0000-0002-7504-031X
Gu Lin 11
http://orcid.org/0000-0003-2176-6403
Nan Ce-Wen 1
http://orcid.org/0000-0002-5652-6701
Lin Yuan-Hua linyh@tsinghua.edu.cn

1
1 grid.12527.33 0000 0001 0662 3178 State Key Laboratory of New Ceramics and Fine Processing, School of Materials Science and Engineering, Tsinghua University, Beijing, PR China
2 https://ror.org/011xvna82 grid.411604.6 0000 0001 0130 6528 Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou, PR China
3 grid.9227.e 0000000119573309 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, China
4 grid.9227.e 0000000119573309 Beijing Synchrotron Radiation Facility, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China
5 https://ror.org/03cve4549 grid.12527.33 0000 0001 0662 3178 Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Department of Engineering Mechanics, Tsinghua University, Beijing, China
6 https://ror.org/02s376052 grid.5333.6 0000 0001 2183 9049 Theory and Simulation of Materials (THEOS), École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
7 grid.499323.6 RICMASS, Rome International Center for Materials Science Superstripes, Via dei Sabelli 119A, Roma, Italy
8 https://ror.org/05vr1c885 grid.412097.9 0000 0000 8645 6375 Henan Key Laboratory of Materials on Deep-Earth Engineering, School of Materials Science and Engineering, Henan Polytechnic University, Jiaozuo, PR China
9 https://ror.org/034t30j35 grid.9227.e 0000 0001 1957 3309 Key Laboratory of Materials Physics Institute of Solid State Physics, Chinese Academy of Sciences, Hefei, China
10 grid.48166.3d 0000 0000 9931 8406 State Key Laboratory of Organic-inorganic Composite, College of Materials Science and Engineering, Beijing University of Chemical Technology, Beijing, PR China
11 https://ror.org/03cve4549 grid.12527.33 0000 0001 0662 3178 National Center of Electron Microscopy in Beijing, School of Materials Science and Engineering, Tsinghua University, Beijing, China
3 9 2024
3 9 2024
2024
15 765011 6 2024
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Thermoelectrics converting heat and electricity directly attract broad attentions. To enhance the thermoelectric figure of merit, zT, one of the key points is to decouple the carrier-phonon transport. Here, we propose an entropy engineering strategy to realize the carrier-phonon decoupling in the typical SrTiO3-based perovskite thermoelectrics. By high-entropy design, the lattice thermal conductivity could be reduced nearly to the amorphous limit, 1.25 W m−1 K−1. Simultaneously, entropy engineering can tune the Ti displacement, improving the weighted mobility to 65 cm2 V−1 s−1. Such carrier-phonon decoupling behaviors enable the greatly enhanced μW/κL of ~5.2 × 103 cm3 K J−1 V−1. The measured maximum zT of 0.24 at 488 K and the estimated zT of ~0.8 at 1173 K in (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3 film are among the best of n-type thermoelectric oxides. These results reveal that the entropy engineering may be a promising strategy to decouple the carrier-phonon transport and achieve higher zT in thermoelectrics.

Authors propose an entropy engineering strategy to realize the carrier-phonon decoupling in SrTiO3-based perovskite thermoelectrics, reducing the lattice thermal conductivity nearly to the amorphous limit and improving the weighted mobility.

Subject terms

Thermoelectrics
Thermoelectric devices and materials
https://doi.org/10.13039/501100001809 National Natural Science Foundation of China (National Science Foundation of China) 52388201 Lin Yuan-Hua issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Thermoelectrics converting electricity and heat directly are promising in waste heat harvesting and active cooling, providing solutions to the fossil fuel crisis and 5G/6G micro-device cooling1. To quantify the comprehensive thermoelectric properties, the dimensionless figure of merit, zT = S2σT/(κe + κL), was employed, where S, σ, T, κe, κL represent the Seebeck coefficient, electrical conductivity, temperature in Kelvin, carrier thermal conductivity and lattice thermal conductivity, respectively2. Accordingly, the simultaneous realization of a high Seebeck coefficient, high electrical conductivity, and low thermal conductivity, makes good thermoelectrics, requiring synergistic modulation of carrier-phonon transport3.

However, the parameters in the zT definition are strongly inter-correlated. The quality factor B=kBe28πe(2mekBT)323h3μWκLT, positively correlated to zT, indicates that the extent of carrier-phonon decoupling reflected by μWκL is the key point to the zT enhancement4. Recently, many efforts have been made to decouple carrier-phonon transport. By modulation of intrinsic crystal symmetry, e.g., Pnma to Cmcm symmetry in SnSe by Pb and Cl co-doping5, rhombohedral to near cubic phase transition in Pb and Sb co-doped GeTe6, the resulting high symmetry could benefit the improvement of the Hall mobility, while the thermal conductivity was suppressed by the dopants. Additionally, through compositing effects by adding metal particles7, band aligning precipitates4,8 in the matrix phase, energy barriers for carrier transport were lowered, and the phonons were simultaneously scattered. Furthermore, interface engineering, e.g., constructing PbTe-SrTe coherent interfaces in PbTe system9, forming Cu2Se-BiCuSeO-graphene interfaces in Cu2Se materials10, was also proved effective to concurrently achieve phonon blocking and charge transmitting9. Despite such efforts in alloys to decouple carriers and phonons effectively, carrier-phonon decoupling is still a big challenge in thermoelectric oxides.

Oxides of low cost, low pollution, high abundance, and excellent thermal stability are competitive in thermoelectrics11. SrTiO3, as a typical perovskite thermoelectric oxide, displays relatively good electrical properties12 (PF = S2σ over 1000 μW m−1 K−2) originating from special TiO6 octahedrons and high symmetry cubic phase13, but suffers from high thermal conductivity over 11 W m−1 K−1 14. Strategies such as nanostructuring15, element doping16, and vacancy modulation17, were applied in reducing lattice thermal conductivity, but the carrier mobility was sacrificed, limiting the further enhancement of zT in SrTiO3-based thermoelectrics18. Therefore, carrier-phonon decoupling is important for zT enhancement in SrTiO3-based oxides. Modifying the grain boundaries of SrTiO3-based ceramics with carbon-based materials could synergistically tune the electrical-thermal transport from the extrinsic compositing aspect19–24. However, it is not quite common to decouple carrier-phonon transport in SrTiO3 by intrinsically tuning the composition and structure.

Recently, the high-entropy strategy has been found to remarkably suppress the lattice thermal conductivity in thermal barrier coatings and thermoelectrics25–27. Normally, ion displacement is strongly related to carrier transport28, which could be regulated in perovskites-based dielectric capacitors by engineering entropy29,30. Herein, we propose that rational entropy engineering could tune the phonon and carrier transport behaviors, which could be expected to decouple the carrier-phonon transport, reaching an overall optimization of the thermoelectric performance in SrTiO3-based oxides.

Results

Design and preparation of entropy-engineered thermoelectrics

As for the crystal sites for entropy engineering, since the Ti 3d orbitals form the conduction band minimum (CBM) deciding the electron transport in n-type SrTiO3-based semiconductors31, and A-site vibrations correspond to acoustic phonon branches and low-frequency optical branches dominating lattice thermal conductivity32,33, the A-site entropy engineering could help maintain the TiO6 octahedrons undisturbed to prevent mobility deterioration and meanwhile strongly scatter phonons34–36, realizing carrier-phonon decoupling. Therefore, SrTiO3, BaTiO3, CaTiO3, and PbTiO3 were selected to form solid solutions. Starting from La-doped SrTiO3 and keeping the doping level (20%) of La fixed, the entropy was engineered by introducing Ba, Ca, and Pb, equimolarly occupying the rest 80% of A sites. (Sr0.8La0.2)TiO3 (SLTO), (Sr0.4Ba0.4La0.2)TiO3 (SBLTO), (Sr0.267Ba0.267Ca0.267La0.2)TiO3 (SBCLTO), and (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3 (SBCPLTO) films were grown on single crystal (LaAlO3)0.3(Sr2TaAlO6)0.7 (LSAT) (001) substrates by pulsed laser deposition (PLD) followed by annealing in reducing atmosphere to create O vacancies and provide electrons to make the films conductors (see below in “Methods” section). The reason for fabricating the thin films, rather than bulks, is that the high-quality thin films could reduce the influence of grain boundaries in perovskites on electrical and thermal transport37,38, making it easier to construct entropy-transport relationship and enhance the thermoelectric properties. According to the definition of the configuration entropy Sconfig.=−R((∑i=1Nxilnxi)cation−site+(∑j=1Mxjlnxj)anion−site), where R, N, M, x are ideal gas constant, the number of cation-site elements, the number of anion-site elements, and the molar ratio of elements39, the Sconfig. increased from low entropy (0.50R in SLTO) to medium entropy (1.05R, 1.38R in SBLTO, SBCLTO, respectively), and to high entropy (1.61R in SBCPLTO without taking O and Pb vacancies into account) (Supplementary Table 1) in this work.

The XRD patterns (Supplementary Figs. 1–3) showed that perovskite thin films were epitaxially grown of high quality by PLD on LSAT (001) substrates, and there were no impure phases before and after the annealing process. The unchanged position and shape of LSAT peaks (Supplementary Fig. 3) indicated chemical endurance to reduce, avoiding electrical contribution from substrates40. The shifts after sequent introduction of Ba, Ca, Pb to (Sr0.8La0.2)TiO3 were consistent with the relative size changes of introduced A-site atoms (Supplementary Fig. 1, Supplementary Table 2), illustrating the successful solid solution at A sites. The epitaxy nature could also be revealed by RSM of (103) peaks (Supplementary Fig. 4), and the peaks of the films deviated from the Qx of LSAT (103), which was a sign for in-plane strain release by thickness of around 200 nm, excluding the influence which epitaxial strain could make on transport. The interface of the thin films and substrates consisted of co-vertex connected Al/TaO6 and TiO6 octahedrons through ABF and HAADF images (Supplementary Figs. 5 and 9h), and the films showed good morphology in TEM (Supplementary Fig. 6), verifying coherently epitaxial thin films of high quality. The elements distributed uniformly at A sites in low-entropy, medium-entropy, and high-entropy samples, and no obvious ordered structure and distribution were viewed, which was reflected by EDS mapping in atomic resolution and micro-scale (Supplementary Figs. 7–11), confirming the random and homogeneous occupation at A sites to fulfill the definition of high entropy. Unavoidable Pb volatilization during target sintering, film deposition, and annealing contributed to the considerable amount of Pb vacancies (Supplementary Table 3, the molar ratio Pb/Ti = 0.073 by EPMA, nominal VPb/Ti = 0.127), adding to the entropy at A sites and tuning the A-site average radius.

Suppressing lattice thermal conductivity to amorphous limit

Since A-O skeleton dominates the transport of acoustic and low-frequency optical phonons17, high level of La doping (20%) in SLTO primarily reduced lattice thermal conductivity from 11 W m−1 K−1 of STO single crystals12 to 2.52 W m−1 K−1 due to size, mass, charge differences and the vacancies introduced (Fig. 1). The entropy engineering at A sites could further strongly suppress the lattice thermal conductivity κL, and thermal conductivity of ~1.60 W m−1 K−1 at room temperature was measured by TDTR in SBCPLTO (Fig. 1a and Supplementary Fig. 26). After increasing entropy by introducing elements, the phonon mean free path lp decreased monotonically (Supplementary Table 7), and the room temperature lattice thermal conductivity of thin films was suppressed significantly from 2.52 W m−1 K−1 for SLTO to 1.25 W m−1 K−1 for SBCPLTO, approaching the amorphous limit41 (Fig. 1b, Supplementary Fig. 12). To verify the thermal measurement of thin films, the lattice thermal conductivity of corresponding bulks was measured by LFA, and the results between TDTR and LFA, bulks and films, were close to some extent in Supplementary Fig. 24.Fig. 1 Thermal transport behaviors of entropy-engineered thin films, and the extrinsic and intrinsic origins of strong phonon scattering after entropy increase.

a Temperature-dependent total thermal conductivity (κ) of entropy-engineered perovskite oxide thin films. b Calculated lattice thermal conductivity (κL) of entropy-engineered thin films and amorphous limit thermal conductivity of SBCPLTO. c The relation between phonon mean free path (lp) and the size disorder parameter (δ) (Supplementary Table 4)43,99 at A sites. The phonon mean free path was measured on corresponding bulks at room temperature for reference (Supplementary Table 7). The definition of size disorder parameter δ can be found in supplementary materials (Supplementary Equation 1). The purple arrow is a guide for the eyes. d, e The DFT calculation of phonon dispersion of SLTO (d) and SBCPLTO (e). f Raman spectra of corresponding entropy-engineered bulks for reference. g–l The TEM, FFT, and GPA results of SLTO (g, h, i) and SBCPLTO (j, k, l). The dashed lines mark the interfaces between the films and substrates. The dislocation symbols are also shown in red.

Intrinsic structural factors and extrinsic defect factors could explain the entropy-scattering relations. The acoustic and low-frequency optical branches of phonons are mainly influenced by A sites in perovskites32, and thus for intrinsic factors, the doping at A sites with different mass, size, charge and vacancies could effectively scatter phonons25,42. The size disorder, described by a standard deviation, size disorder parameter δ (Supplementary Equation 1), could lead to lattice distortion, which is thought to be the main factor suppressing lattice thermal conductivity43. Due to the comparably large radius of Ba2+ (Supplementary Table 2), the significantly increased size disorder had already strongly scattered phonons and lowered the lp (Fig. 1c). In addition to the size disorder coming from different radius of elements, the charge difference from La3+ also could be seen as one kind of size fluctuation. Causing lattice distortion by the compensation for the local charge imbalance25, the charge disorder is usually responsible for the strong reduction in thermal conductivity16,33. The mass disorder could contribute to phonon scattering as well. As suggested by the phonon dispersion calculation in Fig. 1d, e and Supplementary Fig. 13, by introducing mass disorder to SrTiO3, the acoustic phonon branches, and low-frequency optical phonon branches became broader in dispersion, meaning the enhanced scattering rate of phonons33. Furthermore, since the A sites were nominally fully occupied, and the La3+, Sr2+, Ba2+, Ca2+, and Pb2+ have no other smaller valence states when performing pulsed laser deposition, the A-site vacancies tended to form for charge balance in low-pressure, possibly contributing to the reduction in thermal conductivity by large size, charge and mass disorder16,19,33. Consistent with the size disorder parameter δ and DFT calculation (considering mass disorder), in corresponding unannealed bulks, the Raman spectra (Fig. 1f) which are the projection of optical phonon branches at zero point (G) of the Brillouin zone, displayed broadened peaks around 100 cm−1 from SLTO to SBCPLTO, corresponding to increased linewidth of optical branches decided by A-site vibrations44. For extrinsic defects, as illustrated by FFT and GPA process from HRTEM images (Fig. 1g−l), larger amounts of dislocations (as reported in high-entropy thermal barriers27) and denser strain distribution, were found in high-entropy SBCPLTO than those in SLTO with lower entropy, functioning as phonon scatters of different scales. In addition, it is worth noting that though the entropy of SBCLTO is higher than SBLTO, the thermal capacity was not necessarily monotonically correlated with entropy, and the Ca with smaller mass could lead to larger thermal capacity (Dulong–Petit Law and Neumann–Kopp Rule45), thus larger thermal conductivity in SBCLTO than that in SBLTO (Fig. 1b and Supplementary Fig. 24).

Improving the carrier mobility by tuning electron-conducting TiO6 octahedrons

The carrier transport including electrical conductivity σ, Hall mobility μH, and weighted mobility μW46 was measured and calculated (Fig. 2a–c and Supplementary Fig. 27). As the carrier concentration was in the order of ~1021 cm−3 for the films (Supplementary Table 8), far larger than the critical concentration of degeneration47, all the films were degenerate semiconductors. For the absence of grain boundary scattering in epitaxial thin films, the σ-T, μH-T, and μW-T followed the acoustic phonon scattering dominated T−1.5 trend48. SLTO, as a typical thermoelectric material, displayed σ of ~1066 S cm−1 and μH over 2.1 cm2 V−1 s−1 at room temperature, while SBLTO, after the introduction of Ba with a much larger radius of 1.61 Å, showed deteriorated σ of ~175 S cm−1 and μH of ~0.4 cm2 V−1 s−1 at room temperature. Though the lattice thermal conductivity was low enough in SBLTO, but the coupled poor electrical properties harmed the thermoelectric performances. Substituting Ba by Ca, and Pb and simultaneously introducing Pb vacancies with lower radius, the weighted mobility recovered from SBLTO to SBCLTO to SBCPLTO, and SBCPLTO was comparable with SLTO in weighted mobility. The carrier mean free paths lc48 and the deformation potentials Ξdef5 were calculated (Supplementary Table 8), and obviously, with decreased deformation potentials, SBCPLTO obtained increased lc. The structural factors behind the recovery of carrier mobility needed exploring.Fig. 2 Electrical transport behaviors of entropy-engineered thin films and the possible structural origins of decoupled carrier transport explained by spectral methods.

a–c Temperature-dependent electrical conductivity (σ) (a), Hall mobility (μH) (b), and weighted mobility (μW) (c) of entropy-engineered thin films. The acoustic phonon scattering dominated T−1.5 trend curves were plotted (dashed gray line). d The XPS Ti 2p spectra of entropy-engineered thin films. The dashed gray line marks the peak position and serves as a guide for the eyes. e The PDF results of corresponding bulks to get referencing bond length and the fitting curves by reverse Monte Carlo methods. f The Raman spectra of LSAT substrates and entropy-engineered thin films on LSAT substrates. The typical vibration modes of perovskite thin films were highlighted by dashed lines, and the inset is the schematic of LO4 vibration mode. g SHG mapping of SBLTO, SBCLTO, and SBCPLTO in 10 μm × 10 μm.

The relation between the mobility deterioration and the TiO6 distortion was noticed. In Ti 2p XPS (Fig. 2d), there were no obvious Ti3+ shoulder peaks near Ti4+ in all films49, suggesting that Ti mainly existed as Ti4+. In this case, the shifts of Ti and O XPS binding energy peaks (Fig. 2d and Supplementary Fig. 16) were negatively correlated to the bond length. In terms of lattice parameters, SBLTO > SBCPLTO > SBCLTO > SLTO (Supplementary Fig. 14), while according to XPS, for Ti-O bond length, SBLTO > SBCLTO > SBCPLTO ≈ SLTO (Fig. 2d and Supplementary Fig. 16). SBCLTO had longer Ti-O bonds though smaller lattice than SBCPLTO, indicating that the TiO6 octahedrons were distorted in SBCLTO, and SBCPLTO had similar Ti-O chemical surroundings with centrosymmetric SLTO despite different lattice parameters. PDF (pair distribution function) by synchrotron radiation is a good method to get information on coordination and bonds. The reverse Monte Carlo simulation could give the average bond length of atom pairs by fitting the pair distribution function (Fig. 2e), and the observed tolerance factors tobs=length(A−O)2length(B−O) were calculated (Supplementary Table 6)13. The deviation of tobs from 1 in SBLTO and SBCLTO could represent lowered symmetry and distorted TiO6 octahedrons.

Furthermore, the types of TiO6 distortion could also be inferred by spectral methods. In Raman spectra of films (Fig. 2f), compared with substrates, extra peaks from 240 to 400 cm−1 in films were considered to be multiple second-order peaks typical in perovskites, and extra peaks around 795 cm−1 in SBLTO and SBCLTO corresponded to the LO4 mode, indicating off-center vibration of Ti in TiO6 octahedrons50. The absence of LO4 mode in SLTO and SBCPLTO could be thought of as a sign of Ti centrosymmetric vibrations, and the red-shifted and stronger LO4 signal in SBLTO than that in SBCLTO might be from larger displacement and longer bond length. Additionally, from the second harmonic generation (SHG) mapping (Fig. 2g and Supplementary Fig. 17), the breaking of centrosymmetry by displacement could also be inferred. The polar domains were visualized in SBLTO, and the signal weakened in SBCLTO, hence hardly detectable in SBCPLTO with recovered mobility. The Ti displacement seemed to be the type of octahedron distortion accounting for the mobility deterioration.

To validate the TiO6 distortion type, the Cs-corrected STEM was performed. The anti-direction tilting could double the unit cells, and lead to half-index diffraction in the selected area electron diffraction (SAED)51. In SLTO and SBCLTO (Supplementary Fig. 18), the presence of half-index diffraction points in SAED patterns along [110] direction ensured the tilting of TiO6, while no half-index diffraction points were obtained in SBLTO and SBCPLTO (Fig. 3b, c). The octahedron tilting could be seen in atomic resolution annular bright field (ABF) images by the displacement of oxygen ion (Supplementary Fig. 19), and the tilting in SBCLTO was stronger than that in SLTO, explaining the smaller binding energy in Ti 2p XPS of SBCLTO. In terms of B-site displacement, the displacement of Ti ions decreased from SBLTO to SBCLTO, and Ti ions hardly displaced in SBCPLTO (Fig. 3d–f). In short, the distortion types were slight TiO6 tilting in SLTO, pure large Ti displacement in SBLTO, mixed TiO6 tilting and moderate Ti displacement in SBCLTO, and minimized distortion in SBCPLTO, respectively (Supplementary Table 10). As can be seen, the TiO6 octahedrons in SLTO are tilted, and in SBCPLTO, the tilting is ignorable, but the mobility was close between SLTO and SBCPLTO. Furthermore, the TiO6 tilting is ignorable in SBLTO, but the mobility deteriorates significantly, which seems that the tilting is not the main factor affecting carrier mobility. The Ti displacement dominated the electron transport in our entropy-engineered perovskite titanates, and with an increasing displacement of Ti, the weighted mobility decreased (Fig. 3d–g). To explain, the displaced ions could create shorter and longer bonds, and the electrons could be localized in shorter bonds with overlapping electron clouds28 (Fig. 3a). From the aspect of orbitals, the breaking of the symmetry by Ti displacement could result in irregular bond lengths and O-Ti-O angles departing from 180°, and thus the split of t2g orbitals would lead to carrier localization13,52,53. Furthermore, the inhomogeneous displacement among adjacent octahedrons in disordered medium-entropy samples could contribute to inhomogeneous electronic potentials, thus scattering electrons. After improving the symmetry, the similar bond feature could enable homogeneous electron cloud overlapping among Ti-O bonds, delocalizing the electron to transport. In summary, the Ti displacement was the distortion that hindered the electron transport in this case.Fig. 3 Atomic scale electron microscopy characterization to explain the recovered mobility of entropy-engineered thin films.

a Schematic of the possible mechanism of electron scattering and localization in Ti-displaced TiO6 octahedrons. C point marked is the center of the TiO6 octahedron, and the length of Ti displacement is the length of C-Ti. b, c SAED patterns along [110] to prove the absence of octahedron tilting in SBLTO and SBCPLTO. d–f ABF images in the atomic resolution after filtering to show the Ti displacement in SBLTO (d), SBCLTO (e), and SBCPLTO (f). The orange scale bar denotes 5 Å. The orange, red, and blue spheres represent A-site atoms, O, and Ti, and the overlapped spheres were omitted. g The relation between room temperature weighted mobility (μW) and average normalized Ti displacement (d¯Ti) of SBLTO, SBCLTO, and SBCPLTO with increased entropy. The displacement was normalized by the ratio dTi = length(C-Ti)/length(C-Ob) in (a). The error bar is the standard deviation δd. The blue and red arrows are the guides to show the trends of weighted mobility and displacement. The details of the analysis and results can be found in supplementary materials.

Decoupled carrier-phonon transport and enhanced thermoelectric properties

The discussion on the carrier-phonon decoupling is needed. For phonon transport, the low-entropy SLTO displayed relatively high thermal conductivity (Fig. 1a, b), and the medium-entropy SBLTO, SBCLTO, the high-entropy SBCPLTO showed decreased phonon mean free path lp and minimized thermal conductivity (Fig. 4a and Supplementary Table 7). The increase of configuration entropy (Sconfig.) would lead to the increase in size disorder, mass disorder, and dislocation density27, and it could be summarized that the thermal diffusivity D decreased with increasing entropy, as shown in Fig. 4b. In terms of carrier transport, after the introduction of large Ba at A sites, the weighted mobility decreased. Interestingly, the weighted mobility was improved with suppressed phonon mean free path lp from SBLTO to SBCLTO to SBCPLTO by a decreased amount of large Ba cations and substitution by smaller Ca, Pb, and Pb vacancies (Fig. 4a). As discussed above, the displacement of Ti in TiO6 octahedrons was the main factor affecting the carrier transport. The TiO6 octahedrons could be tuned by changing the A-O average bond length, and the tolerance factor t=length(A−O)2length(B−O) was the reason behind13. When the tolerance factor is 1, the perovskite tends to form a cubic phase with undistorted centrosymmetric BO6 octahedrons, and SrTiO3 is a typical material with t ≈ 153. The occupation of large Ba at A sites would deviate the t from 1, distorting the TiO6, and the introduction of smaller ions or vacancies to tailor t back to 1 could realize TiO6 recovery. As concluded in Fig. 4c, with the tobs approaching 1, the TiO6 became less distorted, performing better in weighted mobility. Since entropy is a concept of disorder, and the tolerance factor t is a mean value, it is reasonable to decouple the transport behaviors by these two different math concepts. Additionally, the entropy design was at A sites, dominating phonon transport, while the carrier transport was mainly affected by TiO6 octahedrons, which was also the reason for the carrier-phonon decoupling. In summary, A-site rational entropy engineering resulted in A-site disorder, reducing lattice thermal conductivity, and meanwhile modulated the A-site radius and tolerance factor, thus reducing the Ti displacement and improving carrier mobility. Consequently, the carrier-phonon transport was decoupled in perovskite thermoelectrics via entropy engineering. Due to the decoupled carrier-phonon transport, the μW/κL increased from SBLTO to SBCPLTO, and achieved enhanced μW/κL of ~5.2 × 103 cm3 K J−1 V−1 at room temperature in SBCPLTO, outperforming the original SLTO (Fig. 4d).Fig. 4 Decoupled carrier-phonon transport and enhanced thermoelectric properties of entropy-engineered thin films.

a Schematic diagram of tuning A sites and TiO6 octahedrons to decouple the carrier-phonon transport with increasing entropy. The symbols of elements were shown in the insets. b The correlation between the thermal diffusivity (D) and the nominal configuration entropy (Sconfig.). The thermal diffusivity was measured on corresponding unannealed bulks at 923 K for reference. The purple arrow is a guide for the eyes. c The relation between weighted mobility and observed tolerance factor tobs. to explain the structural origin of mobility recovery. The tolerance factor was calculated from the bond length measured by PDF on corresponding bulks for reference. The purple arrow is a guide for the eyes. d The lattice thermal conductivity (κL), weighted mobility (μW), and μW/κL of SLTO, SBLTO, SBCLTO, and SBCPLTO to show the extent of decoupling and effects of different elements when engineering entropy at room temperature. e, f. Temperature-dependent BE (in log scale) (e), μW/κL (f) of entropy-engineered thin films. g Temperature-dependent zT of SBCPLTO, in comparison with other competitive n-type thermoelectric pure oxides12,58–62,100–102. The zT measurement was limited by the highest allowed measuring temperature of TDTR. Since the PF and κ are plateau-like with increasing T at high temperatures, the zT was extrapolated to high temperature (dashed red line).

As degenerate semiconductors, the Seebeck coefficients (S) of all samples increased linearly with temperature (Supplementary Fig. 23a). High S of −228 μV K−1 was reached in SBCPLTO at 693 K, and through Pisarenko plot at room temperature using SPB model (Supplementary Fig. 23b), it could be seen that the md* of SLTO, SBLTO and SBCLTO were close (~6.5 m0), while enhanced md* of ~8.2 m0 was calculated in SBCPLTO. The A-site doping could help tune the energy band structure and the electrical transport in SrTiO3-based thermoelectrics54,55. To explain the improved md*, the band dispersion was calculated by DFT (Supplementary Fig. 21), and SLTO, SBLTO, and SBCLTO shared similar band structures, while extra orbitals lay at the CBM of SBCPLTO. According to the partial DOS and COHP (Supplementary Figs. 22, 23c, and 29), since Sr, Ba, and Ca belong to the same alkaline-earth main group, with empty s and d orbitals of high energy, the CBM and VBM only consist of Ti-O bonding and anti-bonding orbitals. However, Pb2+ with lone pair 6s2 and empty 6p could affect band structure significantly54, and coupled with Ti 3d and A-site d orbitals, the Pb 6p orbitals with high degeneracy could effectively contribute to the DOS at CBM to enhance the md* of SBCPLTO. Additionally, the increased symmetry could also increase the degeneracy, thus the md* in SBCPLTO56. Benefiting from recovered mobility and enhanced effective mass, the PF and electronic quality factor BE57 of SBCPLTO outperformed the base material SLTO, and PF reached 1.1 × 103 μW m−1 K−2 in SBCPLTO (Fig. 4e and Supplementary Fig. 23d). In SBCPLTO, the decoupled μW/κL (Fig. 4f) and the increased effective mass combined helped significantly enhance the zT compared to low- and medium-entropy samples (Supplementary Fig. 23e), and the maximum zT of 0.24 was achieved at 488 K in SBCPLTO (limited by the TDTR measuring temperature). Considering the PF-T, κL-T plateau at the middle to high-temperature range in SBCPLTO, the zT was roughly extrapolated linearly to 1173 K, reaching a high estimated zT over 0.8, which outperforms other n-type oxide thermoelectrics (Fig. 4g)12,58–62. Due to the absence of grain boundary scattering in epitaxial thin films to enhance room temperature performances and the synergistic entropy engineering to decouple carrier-phonon transport, the average zT (RT to 473 K) of 0.14 was obtained (Supplementary Fig. 23f), which is also competitive in n-type thermoelectric oxides.

Discussion

Entropy increase usually accompanies with deteriorated carrier mobility despite lattice thermal conductivity reduction. By applying entropy engineering at phonon transport skeletons and intentionally tuning the relative sizes of ions under the guidance of factors like tolerance factor to delocalize carriers, the carrier-phonon decoupling could be realized. The strategy could be used in thermoelectrics with separate carrier and phonon transport units. Since the high entropy could suppress lattice thermal conductivity to the amorphous limit, the decoupled carrier transport could further enable competitive thermoelectric properties. This work provides solutions to applying entropy engineering in thermoelectrics and materials facing trade-offs among strongly coupled physical variables.

Methods

Sample preparation

The corresponding ceramic targets for film fabrication and the bulks for referencing characterization were fabricated by conventional solid-state reaction and pressureless sintering in air. The raw powers including La2O3 powder (99.99%, Aladdin, China), TiO2 powder (99.8%, anatase, Aladdin, China), SrCO3 powder (99.95%, Macklin, China), BaCO3 powder (99.95%, Macklin, China), CaCO3 powder (99.99%, Macklin, China), PbTiO3 powder (99.5%, Alfa Aesar, China), were stoichiometrically mixed (PbTiO3 5% excess) by ball milling with ethanol as agent. The dried mix powders were calcined for 3 h in muffle furnaces in air (1523 K for (Sr0.8La0.2)TiO3, (Sr0.4Ba0.4La0.2)TiO3, (Sr0.267Ba0.267Ca0.267La0.2)TiO3, and 1323 K for (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3). The calcined powers were ball milled and were pressed into pellets (30 mm in diameter for targets). After cold isotropic pressing at 220 MPa, the pellets were sintered for 10 h in muffle furnaces in the air (1723 K for (Sr0.8La0.2)TiO3, (Sr0.4Ba0.4La0.2)TiO3, (Sr0.27Ba0.27Ca0.27La0.2)TiO3, and 1523 K for (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3), and the dense ceramic targets were in single phases (Supplementary Fig. 25). The corresponding bulk ceramic samples for reference were prepared at the same time with targets. The corresponding bulks were not annealed to remain electrically insulated to exclude the carrier contribution to thermal transport for reference.

The thin films were grown with nominal contents (Sr0.8La0.2)TiO3, (Sr0.4Ba0.4La0.2)TiO3, (Sr0.27Ba0.27Ca0.27La0.2)TiO3, (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3 via pulsed laser deposition (PLD). The single crystal LSAT (001) substrates were used (12 × 4 × 0.5 mm3, Hefei Kejing Materials Technology Co., Ltd., China), and the high-temperature silver paste (05001-AB, 0.5 oz, Spi Supplies, the U.S.) was applied to ensure uniform temperature when depositing films. The PLD was performed with a KrF laser (wavelength = 248 nm), and the energy density, repetition rate, number of pulses, growth temperature, working oxygen partial pressure were 0.8 J cm−2, 5 Hz, 6000 pulses, 1023 K, 20 Pa, respectively. After deposition, the films were annealed in situ in oxygen (~600 Pa) for 20 min to improve the crystal quality and were cooled down to room temperature at a rate of 10 K min−1. The as-grown thin films underwent annealing in 5% H2/95% Ar mixed gas at 1173 K for 2 h to be reduced. After the annealing in reducing gas, the thin films turned into conductors, while the LSAT substrates which are not sensitive to oxygen vacancies, were tested to remain insulated, excluding the electrical contribution from the substrates40.

Measurement of transport behaviors

Electrical transport

The temperature-dependent electrical conductivity and Seebeck measurement of films were performed on the commercial apparatus MRS-3 (JouleYacht, China). The temperature-dependent carrier concentration (nH) and the Hall mobility (μH) were calculated from σH (μH = σH|RH|) and RH (nH = 1/eRH) by van der Pauw method in Hall measurement (8408, Lake Shore Cryotronics, Inc.). The effective mass m* could be derived by combining Hall data with the Seebeck coefficient. In this work, the mobility of samples varies from 0.1 to 2.2 cm2 V−1 s−1, right at the test limit of DC field measurement (1–1 × 106 cm2 V−1 s−1). Therefore, the AC field measurement developed by the Lake Shore Cryotronics63,64, using AC magnetic field rather than the traditional DC mode, was applied to remove the effect of misalignment, to enlarge the measuring range of Hall mobility (1 × 10−3 to 1 × 106 cm2 V−1 s−1) and to make the measurement more precise. According to the manual given by the Lake Shore official website, all field values in Hall measurement were nominal and could vary ±1%”, and the error of μH could be calculated by the sum of the error (1%) of five field-related variants (VH, I, B, I’, U’) to be 5%. The μH–T correlation with an error of 5% was plotted in Supplementary Fig. 27a. However, multiple factors could enlarge the error of Hall measurement, like measuring geometry, electrical contact, shape and size of metal pads, etc., and the error could even reach 20%. μH–T correlation with exaggerated errors of 20%, and 30% were plotted in Supplementary Figs. 27b and 27c. As can be seen in Supplementary Fig. 27, the μH differences between samples are still more significant than the error bar.

Thermal transport

The thermal properties of annealed LSAT(001) and annealed films were measured by the time-domain thermoreflectance (TDTR) technique65,66. TDTR utilizes a pump-probe configuration, where the pump beam heats the sample, and the probe beam monitors the transient sample surface temperature changes via thermoreflectance. The thermal properties are extracted by fitting the thermoreflectance signal with a thermal diffusion model. In this work, the pump beam was modulated with a frequency of 10.7 MHz, and focused on the sample surface together with the probe beam through a 10× objective lens, with a spot radius of 5.2 μm. Prior to the TDTR measurements, a 151 nm aluminum (Al) film was deposited on the film. The Al transducer layer thickness was determined by X-ray reflectance (XRR), and the thickness of the perovskite oxide film was determined by cross-sectional SEM images. The surface of Al film could be oxidized with a native oxide layer (Al2O3) of 2–3 nm. The native oxide layer is commonly counted as an extra Al thickness of ~3 nm in other groups67,68. The 3 nm was added to the measured Al thickness (154 nm in total in this work), and the total thickness was used in the fitting process of the TDTR method. During the temperature-dependent thermal conductivity TDTR measurement, the chamber was vacuumed, and it could be assumed that the native oxide layer would not be further oxidized during high-temperature measurement. The approximation of adding 3 nm to the Al layer still worked in the measurement. The heat capacities of the Al transducer and the annealed LSAT (001) substrate were taken from literature value and laser flash analysis (LFA), respectively, and the heat capacities of perovskite oxide films were estimated by DSC and LFA. The thermal conductivity of the annealed LSAT substrate was measured separately using LFA (κ=CpρD). It is important to clarify that the electrical transport in this work was measured in plane, while the thermal conductivity measured by the TDTR method was out of plane due to difficulty in technique. Additionally, the temperature range of the TDTR measurement was up to 473 K limited by the apparatus. Fortunately, for films whose thicknesses were ~ 200 nm, the mean free paths of phonons in entropy-engineered oxides were smaller by orders (in 100 Å), leading to little influence from transport confinement by thickness25,69. And for cubic or pseudo-cubic phases, the 2nd-order tensor-like thermal conductivity could be isotropic in the strain-released epitaxial oxides with few grain boundaries70, the zT could be roughly calculated despite measuring in different directions. Since the PF and the thermal conductivity displayed plateau-like behaviors with increasing temperature, the zT value could be estimated by fitting the data linearly to high temperature. The laser flash apparatus (LFA, LFA457, NETZSCH, Germany) measured the thermal diffusivity and the thermal capacity of substrates and corresponding ceramics for reference, and standard samples were measured simultaneously to ensure validity. The transverse and longitudinal phonon group velocity of unannealed corresponding oxide ceramics were measured by ultrasonic pulse-echo method (Olympus 5072PR, Japan).

Characterizations

HAADF-STEM and EDS and diffraction

The cross-sectional scanning transmission electron microscope (STEM) samples were acquired by focus ion beam (FIB) milling. The samples were thinned to 50 nm by FIB with a current of 240–50 pA under an accelerating voltage of 30 kV, and with 20 pA and 5 kV during polishing process. High angle annular dark field (HAADF), annular bright field (ABF), and energy dispersion spectroscopy (EDS) images along <110> zone axis with a high resolution of 0.059 nm were collected by a 200 kV JEOL ARM 200CF equipped with double aberration correctors. The convergence angle and the collection angle were set at 25 mrad, 48–2000 mrad, respectively. The exposure time when collecting EDS signals was 2 μs, and the total time was 90 min. The noise in ABF images was filtered by HREM-Filters released by HREM Research Inc.

X-ray total scattering experiment

X-ray diffraction measurements were made at beamline 3W1 of the Beijing Synchrotron Radiation (Beijing, China) using an incident X-ray beam of wavelength 0.206468 Å (60.05 keV). A Mercu 1717HS detector (2048 × 2048 pixels of a 140 × 140-µm CsI scintillator) was placed ~150 mm downstream of the sample, giving a Qmax of 25 Å−1. The setup was calibrated using the diffraction pattern from polycrystalline CeO2 powder. The unannealed ceramic bulks (Sr0.8La0.2)TiO3, (Sr0.4Ba0.4La0.2)TiO3, (Sr0.27Ba0.27Ca0.27La0.2)TiO3, (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3 were glued with Compton tapes to an aluminum alloy frame. The measurement procedure was controlled by iDetector software and 20 s exposure time was set for (Sr0.4Ba0.4La0.2)TiO3, (Sr0.27Ba0.27Ca0.27La0.2)TiO3, and (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3, 10 s for (Sr0.8La0.2)TiO3 due to its stronger signal. Background patterns were collected with the same setup and exposure time. The raw diffraction data were reduced from two-dimensional images and corrected for the effects of polarization and geometry using the program Fit2D71 and absorption, geometry, detector effects, and the normalization procedure carried out using PDFgetX272.

Reverse Monte Carlo (RMC) simulation

Using RMC simulation, the structural models were derived from the diffraction data73. To start, the configuration of ~5000 atoms was used. Constraints were set, including the connectivity of Ti-O bonds (i. e. Ti was coordinated to a reasonable number of O up to 2.5 Å) and the atom-atom approaches. The aim of setting the constraints was to avoid the results of unrealistic structures in physics. By counting the atomic configurations generated by RMC simulations, the structural information was derived.

SHG

Second harmonic generation (SHG) spectroscopy was performed under a home-designed SHG microscope74. The excitation laser was provided by a Ti: sapphire mode-locking femtosecond laser (MaiTai SP, Spectra-Physics) which generated 35 fs pulses with a repetition rate of 80 MHz and a wavelength centered at 800 nm. The fundamental laser beam was directed onto the sample at normal incidence and focused to a focal spot diameter of ≈1 μm using a 50× (numerical aperture NA = 0.55) objective. The back-scattered SHG signal was collected with the same objective and detected using a Hamamatsu photomultiplier tube. The SHG microscopy images of the square areas selected were acquired by scanning the probing position using a pair of galvanometers. When scanning the interested region, the angle between the polarization direction of the incident light and the SHG light collected was fixed the same at 90°. The rate of scanning was ~50 ms pixel−1 and ~2 μm s−1.

Other regular characterizations

The out-of-plane X-ray diffraction (XRD) in two theta-omega mode and the reciprocal space mapping (RSM) of perovskite films were conducted by X-ray diffractor (PANalytical X’Pert MRD, Netherlands) to characterize the epitaxy nature and crystal quality. The cross-sectional scanning electron microscopy (SEM, MERLIN VP Compact, ZEISS, Germany) was used to obtain the thickness of films. Cross-sectional morphology and microstructure were observed by transmission electron microscopy (TEM, JEM-2100F, JEOL, Japan) to analyze the epitaxy nature, lattice fringes, dislocation analysis by FFT, and the energy dispersive spectra (EDS) of films. The focused ion beam (FIB, TESCAN S9000X, Brno, Czech) was applied to thin the samples. Geometric phase analysis75–77 (GPA, HREM Research Inc.) plug-in helped the strain analysis and mapping from HRTEM images of films. Raman spectroscopy (HORIBA, Japan) was utilized to get information on bond vibrations and crystal symmetry, and the laser wavelength was 325 nm for thin films to minimize the signals of substrates and 532 nm for bulks. The information of bond length and charge valence were recorded by X-ray photoelectron spectroscopy (XPS, Thermo Fisher ESCALAB 250Xi, the U.S.), and the bandgaps were measured by ultraviolet-visible light spectrophotometer (UV-vis, L950, PerkinElmer, the U.S.) absorption mode on films. The rough concentration of elements in films, especially Pb, was derived from the Electron probe X-ray micro analyzer (EPMA, JXA8230, JEOL, Japan) by the Pb: Ti ratio.

Calculation

Band structure DFT

The electron density of states (DOS), crystal orbital Hamilton population (COHP)78, and band structures were calculated through density functional theory (DFT) based first-principle scheme using the Vienna ab-initio Simulation Package (VASP)79–81. The electron-ion interaction was modeled by projector-augmented wave (PAW) method82, and the Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA)83 was used to simulate the electron interactions. The doping elements of La, Ca, Ba, and Pb are treated by using virtual crystal approximation (VCA). The cut-off energy was 400 eV for all the computations. The O and Pb vacancies were not considered in the calculation.

Phonon dispersion calculation

First, we trained a deep learning potential model to describe the potential energy surface between atoms. The dataset was generated by DP-GEN84,85, using 3 × 3 × 3 and 2 × 2 × 2 supercells of the cubic primitive cells, and the 2 × 1 × 1 supercell of the tetragonal conventional cell of SrTiO3 as approximation86,87 (Supplementary Fig. 28 and CIF files in Supplementary Data 1–4), and labeled through first-principles calculations, consisting of a total of 2028 data points. The method and supercell size have been applied in phonon dispersion calculation of SrTiO3 in other works86,87. Furthermore, in the self-consistent phonon calculations by Tadano et al. in 2015, a 2 × 2 × 2 supercell is sufficient to obtain phonon dispersion close to the experimental results88. The first-principle calculations were performed using the VASP package79,80,89,90 with the projector-augmented plane wave (PAW)82 method and the PBEsol functional83,91. The energy cutoff was set to 700 eV, and the K-point grid spacing was set to 0.14 Å−1. Next, utilizing the labeled dataset, we employed DeePMD-kit92 to train a neural network. Compared to the results obtained from DFT calculations, the energy accuracy achieved 1.21 × 10−3 eV per atom. Using the deep learning potential model, we applied the SCAILD93,94 method to calculate the interatomic force constants of SrTiO3 at 300 K within a 5 × 5 × 5 supercell. To introduce mass disorder, structures were randomly generated for SrTiO3, (Sr0.8La0.2)TiO3, (Sr0.4Ba0.4La0.2)TiO3, (Sr0.27Ba0.27Ca0.27La0.2)TiO3, (Sr0.2Ba0.2Ca0.2Pb0.2La0.2)TiO3 subsequently basing on the supercell (CIF files in Supplementary Data 5–9). Finally, the UPHO95,96 package was then employed to compute the unfolded phonon spectrum by combining the IFCs and the doped structure with the lowest energy for each composition. In SCAILD and UPHO calculation, the api of phonopy package was used97,98. The O and Pb vacancies were not considered in the calculation.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Supplementary Data 8

Supplementary Data 9

Lasing Reporting Summary

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52063-5.

Acknowledgements

This work was financially supported by the Basic Science Center Project of the National Natural Science Foundation of China under grant No. 52388201 (Y.-H. Lin), and the National Science Foundation of China under grant No. 52172211 (M.Z.). This work (C.S.) was also supported by the National Key Basic Research Program of China (grant No. 2020YFA0406101). Q.Z. and L.G. were supported by the NSFC (grant No. 52025025 and No. 52072400). Z.Z. acknowledges financial support from the Shuimu Tsinghua Scholar Program and the China Postdoctoral Science Foundation under grant No. 2024T170479. W.X. was supported by the National Natural Science Foundation of China (grant No. 12075273). J.-L.L. acknowledges the National Natural Science Foundation of China under grant No. 52372075. B.W. acknowledges the support of the National Science Foundation of China (grant No. 12304040), China Postdoctoral Science Foundation (grant No. 2023M742004), and the State Key Laboratory of New Ceramic and Fine Processing Tsinghua University (grant No. KF202304). PDF measurements were conducted at beamline 3W1 of the Beijing Synchrotron Radiation (Beijing, China). We also acknowledge the contribution of 4B9B beamline of Beijing Synchrotron Radiation (Beijing, China) for the exploration on XAS measurement. We thank the Beijing PARATERA Tech CO., Ltd. (https://cloud.paratera.com). We thank Kaiqi Nie, Zunqiu Xiao, Xiaomin Jia, Jiasheng Guo, Haojie Han, and Shun Lan for fruitful discussions.

Author contributions

Y.-H. Lin and Y. Zheng conceived the study. Y. Zheng, Z. Zhou, and Y.-H. Lin wrote the manuscript. Y. Zheng conducted the experiments and analyzed the data in this study. Q. Zhang and L. Gu performed the Cs-corrected STEM. C. Shi, W. Xu, Y. Zhang, and D. Yang carried out the PDF measurement, and C. Shi analyzed the results by reverse Monte Carlo simulation. Y. Lu, W. Ma, and X. Zhang measured thermal conductivity by TDTR. J. Han and C. Lin calculated the phonon dispersion. H. Chen and D. Yi helped with the calculation of energy bands. Y. Ma and Q. Li supported by SHG mapping. Y. Yang, B. Wei, B. Yang, M. Zou, W. Zhang, C. Liu, L. Dou, J.-L. Lan, and C.-W. Nan assisted in the writing of the manuscript. All authors discussed and revised the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Data availability

The authors declare that the data supporting the findings of this study are available within the paper and its Supplementary Information files. The source data used in this study are available in the Figshare database under accession code 10.6084/m9.figshare.26797660. Any other relevant data are also available upon reasonable request from Y.-H. L.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Yunpeng Zheng, Qinghua Zhang, Caijuan Shi, Zhifang Zhou.
==== Refs
References

1. Shi XL Zou J Chen ZG Advanced thermoelectric design: from materials and structures to devices Chem. Rev. 2020 120 7399 7515 10.1021/acs.chemrev.0c00026 32614171
Shi, X. L., Zou, J. & Chen, Z. G. Advanced thermoelectric design: from materials and structures to devices. Chem. Rev. 120, 7399–7515 (2020).32614171 10.1021/acs.chemrev.0c00026
2. Ioffe. Semiconductor Thermoelements and Thermo-Electric Cooling. Solar Energy, Vol. 4 (INFOSEARCH LIMITED, London, 1960).
3. Xiao Y Zhao L-D Seeking new, highly effective thermoelectrics Science 2020 367 1196 1197 10.1126/science.aaz9426 32165572
Xiao, Y. & Zhao, L.-D. Seeking new, highly effective thermoelectrics. Science 367, 1196–1197 (2020).32165572 10.1126/science.aaz9426
4. Xiao Y Band sharpening and band alignment enable high quality factor to enhance thermoelectric performance in n-type PbS J. Am. Chem. Soc. 2020 142 4051 4060 10.1021/jacs.0c00306 32017566
Xiao, Y. et al. Band sharpening and band alignment enable high quality factor to enhance thermoelectric performance in n-type PbS. J. Am. Chem. Soc. 142, 4051–4060 (2020).32017566 10.1021/jacs.0c00306
5. Su L High thermoelectric performance realized through manipulating layered phonon-electron decoupling Science 2022 375 1385 1389 10.1126/science.abn8997 35324303
Su, L. et al. High thermoelectric performance realized through manipulating layered phonon-electron decoupling. Science 375, 1385–1389 (2022).35324303 10.1126/science.abn8997
6. Li J High-performance GeTe thermoelectrics in both rhombohedral and cubic phases J. Am. Chem. Soc. 2018 140 16190 16197 10.1021/jacs.8b09147 30360620
Li, J. et al. High-performance GeTe thermoelectrics in both rhombohedral and cubic phases. J. Am. Chem. Soc. 140, 16190–16197 (2018).30360620 10.1021/jacs.8b09147
7. Li JW Wide-temperature-range thermoelectric n-type Mg3(Sb,Bi)2 with high average and peak zT values Nat. Commun. 2023 14 7428 10.1038/s41467-023-43228-9 37973881
Li, J. W. et al. Wide-temperature-range thermoelectric n-type Mg3(Sb,Bi)2 with high average and peak zT values. Nat. Commun. 14, 7428 (2023).37973881 10.1038/s41467-023-43228-9
8. Luo ZZ Valence disproportionation of GeS in the PbS matrix forms Pb5Ge5S12 inclusions with conduction band alignment leading to high n-type thermoelectric performance J. Am. Chem. Soc. 2022 144 7402 7413 10.1021/jacs.2c01706 35420804
Luo, Z. Z. et al. Valence disproportionation of GeS in the PbS matrix forms Pb5Ge5S12 inclusions with conduction band alignment leading to high n-type thermoelectric performance. J. Am. Chem. Soc. 144, 7402–7413 (2022).35420804 10.1021/jacs.2c01706
9. Biswas K Strained endotaxial nanostructures with high thermoelectric figure of merit Nat. Chem. 2011 3 160 166 10.1038/nchem.955 21258390
Biswas, K. et al. Strained endotaxial nanostructures with high thermoelectric figure of merit. Nat. Chem. 3, 160–166 (2011).21258390 10.1038/nchem.955
10. Zhou Z Compositing effects for high thermoelectric performance of Cu2Se-based materials Nat. Commun. 2023 14 2410 10.1038/s41467-023-38054-y 37105970
Zhou, Z. et al. Compositing effects for high thermoelectric performance of Cu2Se-based materials. Nat. Commun. 14, 2410 (2023).37105970 10.1038/s41467-023-38054-y
11. Lin, Y. H., Lan, J. & Nan, C. Oxide Thermoelectric Materials: From Basic Principles to Applications. Oxide Thermoelectric Materials: From Basic Principles to Applications (John Wiley & Sons, 2019).
12. Wang J Record high thermoelectric performance in bulk SrTiO3 via nano-scale modulation doping Nano Energy 2017 35 387 395 10.1016/j.nanoen.2017.04.003
Wang, J. et al. Record high thermoelectric performance in bulk SrTiO3 via nano-scale modulation doping. Nano Energy 35, 387–395 (2017).10.1016/j.nanoen.2017.04.003
13. Tilley, R. J. D. Perovskites: Structure-Property Relationships. Perovskites: Structure-Property Relationships (John Wiley & Sons, 2016).
14. Muta H Kurosaki K Yamanaka S Thermoelectric properties of reduced and La-doped single-crystalline SrTiO3 J. Alloy. Compd. 2005 392 306 309 10.1016/j.jallcom.2004.09.005
Muta, H., Kurosaki, K. & Yamanaka, S. Thermoelectric properties of reduced and La-doped single-crystalline SrTiO3. J. Alloy. Compd. 392, 306–309 (2005).10.1016/j.jallcom.2004.09.005
15. Wang Y Interfacial thermal resistance and thermal conductivity in nanograined SrTiO3 Appl. Phys. Express 2010 3 1 4 10.1143/APEX.3.031101
Wang, Y. et al. Interfacial thermal resistance and thermal conductivity in nanograined SrTiO3. Appl. Phys. Express 3, 1–4 (2010).10.1143/APEX.3.031101
16. Daniels LM A and B site doping of a phonon-glass perovskite oxide thermoelectric J. Mater. Chem. A 2018 6 15640 15652 10.1039/C8TA03739F
Daniels, L. M. et al. A and B site doping of a phonon-glass perovskite oxide thermoelectric. J. Mater. Chem. A 6, 15640–15652 (2018).10.1039/C8TA03739F
17. Popuri SR Glass-like thermal conductivity in SrTiO3 thermoelectrics induced by A-site vacancies RSC Adv. 2014 4 33720 33723 10.1039/C4RA06871H
Popuri, S. R. et al. Glass-like thermal conductivity in SrTiO3 thermoelectrics induced by A-site vacancies. RSC Adv. 4, 33720–33723 (2014).10.1039/C4RA06871H
18. Lou Z A novel high-entropy perovskite ceramics Sr0.9La0.1(Zr0.25Sn0.25Ti0.25Hf0.25)O3 with low thermal conductivity and high Seebeck coefficient J. Eur. Ceram. Soc. 2022 42 3480 3488 10.1016/j.jeurceramsoc.2022.02.053
Lou, Z. et al. A novel high-entropy perovskite ceramics Sr0.9La0.1(Zr0.25Sn0.25Ti0.25Hf0.25)O3 with low thermal conductivity and high Seebeck coefficient. J. Eur. Ceram. Soc. 42, 3480–3488 (2022).10.1016/j.jeurceramsoc.2022.02.053
19. Lin Y Thermoelectric power generation from lanthanum strontium titanium oxide at room temperature through the addition of graphene ACS Appl. Mater. Interfaces 2015 7 15898 15908 10.1021/acsami.5b03522 26095083
Lin, Y. et al. Thermoelectric power generation from lanthanum strontium titanium oxide at room temperature through the addition of graphene. ACS Appl. Mater. Interfaces 7, 15898–15908 (2015).26095083 10.1021/acsami.5b03522
20. Feng X Graphene promoted oxygen vacancies in perovskite for enhanced thermoelectric properties Carbon N. Y. 2017 112 169 176 10.1016/j.carbon.2016.11.012
Feng, X. et al. Graphene promoted oxygen vacancies in perovskite for enhanced thermoelectric properties. Carbon N. Y. 112, 169–176 (2017).10.1016/j.carbon.2016.11.012
21. Dylla MT Kuo JJ Witting I Snyder GJ Grain boundary engineering nanostructured SrTiO3 for thermoelectric applications Adv. Mater. Interfaces 2019 6 1 7 10.1002/admi.201900222
Dylla, M. T., Kuo, J. J., Witting, I. & Snyder, G. J. Grain boundary engineering nanostructured SrTiO3 for thermoelectric applications. Adv. Mater. Interfaces 6, 1–7 (2019).10.1002/admi.201900222
22. Okhay O Thermoelectric performance of Nb-doped SrTiO3 enhanced by reduced graphene oxide and Sr deficiency cooperation Carbon N. Y. 2019 143 215 222 10.1016/j.carbon.2018.11.023
Okhay, O. et al. Thermoelectric performance of Nb-doped SrTiO3 enhanced by reduced graphene oxide and Sr deficiency cooperation. Carbon N. Y. 143, 215–222 (2019).10.1016/j.carbon.2018.11.023
23. Rahman JU Grain boundary interfaces controlled by reduced graphene oxide in nonstoichiometric SrTiO3-δ thermoelectrics Sci. Rep. 2019 9 1 12 10.1038/s41598-019-45162-7 30626917
Rahman, J. U. et al. Grain boundary interfaces controlled by reduced graphene oxide in nonstoichiometric SrTiO3 -δ thermoelectrics. Sci. Rep. 9, 1–12 (2019).30626917 10.1038/s41598-019-45162-7
24. Cao J Modulation of charge transport at grain boundaries in SrTiO3: toward a high thermoelectric power factor at room temperature ACS Appl. Mater. Interfaces 2021 13 11879 11890 10.1021/acsami.0c21699 33660994
Cao, J. et al. Modulation of charge transport at grain boundaries in SrTiO3: toward a high thermoelectric power factor at room temperature. ACS Appl. Mater. Interfaces 13, 11879–11890 (2021).33660994 10.1021/acsami.0c21699
25. Braun JL Charge-induced disorder controls the thermal conductivity of entropy-stabilized oxides Adv. Mater. 2018 30 1 8 10.1002/adma.201805004
Braun, J. L. et al. Charge-induced disorder controls the thermal conductivity of entropy-stabilized oxides. Adv. Mater. 30, 1–8 (2018).10.1002/adma.201805004
26. Jiang B High-entropy-stabilized chalcogenides with high thermoelectric performance Science 2021 371 830 834 10.1126/science.abe1292 33602853
Jiang, B. et al. High-entropy-stabilized chalcogenides with high thermoelectric performance. Science 371, 830–834 (2021).33602853 10.1126/science.abe1292
27. Han Y Ultra-dense dislocations stabilized in high entropy oxide ceramics Nat. Commun. 2022 13 2871 10.1038/s41467-022-30260-4 35610224
Han, Y. et al. Ultra-dense dislocations stabilized in high entropy oxide ceramics. Nat. Commun. 13, 2871 (2022).35610224 10.1038/s41467-022-30260-4
28. Jiang B High figure-of-merit and power generation in high-entropy GeTe-based thermoelectrics Science 2022 377 208 213 10.1126/science.abq5815 35857539
Jiang, B. et al. High figure-of-merit and power generation in high-entropy GeTe-based thermoelectrics. Science 377, 208–213 (2022).35857539 10.1126/science.abq5815
29. Yang B Engineering relaxors by entropy for high energy storage performance Nat. Energy 2023 8 956 964 10.1038/s41560-023-01300-0
Yang, B. et al. Engineering relaxors by entropy for high energy storage performance. Nat. Energy 8, 956–964 (2023).10.1038/s41560-023-01300-0
30. Zhang M Ultrahigh energy storage in high-entropy ceramic capacitors with polymorphic relaxor phase Science 2024 384 185 189 10.1126/science.adl2931 38603510
Zhang, M. et al. Ultrahigh energy storage in high-entropy ceramic capacitors with polymorphic relaxor phase. Science 384, 185–189 (2024).38603510 10.1126/science.adl2931
31. Dylla MT Kang SD Snyder GJ Effect of two-dimensional crystal orbitals on fermi surfaces and electron transport in three-dimensional perovskite oxides Angew. Chem. Int. Ed. 2019 58 5503 5512 10.1002/anie.201812230
Dylla, M. T., Kang, S. D. & Snyder, G. J. Effect of two-dimensional crystal orbitals on fermi surfaces and electron transport in three-dimensional perovskite oxides. Angew. Chem. Int. Ed. 58, 5503–5512 (2019).10.1002/anie.201812230
32. Popuri SR Phonon-glass and heterogeneous electrical transport in a-site-deficient SrTiO3 J. Phys. Chem. C. 2019 123 5198 5208 10.1021/acs.jpcc.8b10520
Popuri, S. R. et al. Phonon-glass and heterogeneous electrical transport in a-site-deficient SrTiO3. J. Phys. Chem. C. 123, 5198–5208 (2019).10.1021/acs.jpcc.8b10520
33. Daniels LM Phonon-glass electron-crystal behaviour by A site disorder in n-type thermoelectric oxides Energy Environ. Sci. 2017 10 1917 1922 10.1039/C7EE01510K
Daniels, L. M. et al. Phonon-glass electron-crystal behaviour by A site disorder in n-type thermoelectric oxides. Energy Environ. Sci. 10, 1917–1922 (2017).10.1039/C7EE01510K
34. Zheng Y Electrical and thermal transport behaviours of high-entropy perovskite thermoelectric oxides J. Adv. Ceram. 2021 10 377 384 10.1007/s40145-021-0462-5
Zheng, Y. et al. Electrical and thermal transport behaviours of high-entropy perovskite thermoelectric oxides. J. Adv. Ceram. 10, 377–384 (2021).10.1007/s40145-021-0462-5
35. Zhang P High-entropy MTiO3 perovskite oxides with glass-like thermal conductivity for thermoelectric applications J. Alloy. Compd. 2023 937 168366 10.1016/j.jallcom.2022.168366
Zhang, P. et al. High-entropy MTiO3 perovskite oxides with glass-like thermal conductivity for thermoelectric applications. J. Alloy. Compd. 937, 168366 (2023).10.1016/j.jallcom.2022.168366
36. Zhang P Reduced lattice thermal conductivity of perovskite-type high-entropy (Ca0.25Sr0.25Ba0.25RE0.25)TiO3 ceramics by phonon engineering for thermoelectric applications J. Alloy. Compd. 2022 898 162858 10.1016/j.jallcom.2021.162858
Zhang, P. et al. Reduced lattice thermal conductivity of perovskite-type high-entropy (Ca0.25Sr0.25Ba0.25RE0.25)TiO3 ceramics by phonon engineering for thermoelectric applications. J. Alloy. Compd. 898, 162858 (2022).10.1016/j.jallcom.2021.162858
37. Huang J Simultaneously breaking the double Schottky barrier and phonon transport in SrTiO3-based thermoelectric ceramics via two-step reduction ACS Appl. Mater. Interfaces 2020 12 52721 52730 10.1021/acsami.0c16084 33182996
Huang, J. et al. Simultaneously breaking the double Schottky barrier and phonon transport in SrTiO3-based thermoelectric ceramics via two-step reduction. ACS Appl. Mater. Interfaces 12, 52721–52730 (2020).33182996 10.1021/acsami.0c16084
38. Koumoto K Wang Y Zhang R Kosuga A Funahashi R Oxide thermoelectric materials: a nanostructuring approach Annu. Rev. Mater. Res. 2010 40 363 394 10.1146/annurev-matsci-070909-104521
Koumoto, K., Wang, Y., Zhang, R., Kosuga, A. & Funahashi, R. Oxide thermoelectric materials: a nanostructuring approach. Annu. Rev. Mater. Res. 40, 363–394 (2010).10.1146/annurev-matsci-070909-104521
39. Sarkar, A. et al. High-entropy oxides: fundamental aspects and electrochemical properties. Adv. Mater. 31, 1806236 (2019).
40. Scullin ML Anomalously large measured thermoelectric power factor in Sr1-x Lax TiO3 thin films due to SrTiO3 substrate reduction Appl. Phys. Lett. 2008 92 1 4 10.1063/1.2916690
Scullin, M. L. et al. Anomalously large measured thermoelectric power factor in Sr1-x Lax TiO3 thin films due to SrTiO3 substrate reduction. Appl. Phys. Lett. 92, 1–4 (2008).10.1063/1.2916690
41. Cahill DG Watson SK Pohl RO Lower limit to the thermal conductivity of disordered crystals Phys. Rev. B 1992 46 6131 6140 10.1103/PhysRevB.46.6131
Cahill, D. G., Watson, S. K. & Pohl, R. O. Lower limit to the thermal conductivity of disordered crystals. Phys. Rev. B 46, 6131–6140 (1992).10.1103/PhysRevB.46.6131
42. Lim, M. et al. Influence of mass and charge disorder on the phonon thermal conductivity of entropy stabilized oxides determined by molecular dynamics simulations. J. Appl. Phys. 125, 055105 (2019).
43. Ma Z High thermoelectric performance and low lattice thermal conductivity in lattice-distorted high-entropy semiconductors AgMnSn1-xPbxSbTe4 Chem. Mater. 2022 34 8959 8967 10.1021/acs.chemmater.2c02344
Ma, Z. et al. High thermoelectric performance and low lattice thermal conductivity in lattice-distorted high-entropy semiconductors AgMnSn1 -xPbxSbTe4. Chem. Mater. 34, 8959–8967 (2022).10.1021/acs.chemmater.2c02344
44. Rychetský I Dielectric, infrared, and Raman response of undoped (formula presented) ceramics: evidence of polar grain boundaries Phys. Rev. B Condens. Matter Mater. Phys. 2001 64 1 10
Rychetský, I. et al. Dielectric, infrared, and Raman response of undoped (formula presented) ceramics: evidence of polar grain boundaries. Phys. Rev. B Condens. Matter Mater. Phys. 64, 1–10 (2001).
45. Leitner J Voňka P Sedmidubský D Svoboda P Application of Neumann-Kopp rule for the estimation of heat capacity of mixed oxides Thermochim. Acta 2010 497 7 13 10.1016/j.tca.2009.08.002
Leitner, J., Voňka, P., Sedmidubský, D. & Svoboda, P. Application of Neumann-Kopp rule for the estimation of heat capacity of mixed oxides. Thermochim. Acta 497, 7–13 (2010).10.1016/j.tca.2009.08.002
46. Snyder GJ Weighted mobility Adv. Mater. 2020 32 1 5 10.1002/adma.202001537
Snyder, G. J. et al. Weighted mobility. Adv. Mater. 32, 1–5 (2020).10.1002/adma.202001537
47. Edwards PP Sienko MJ Universality aspects of the metal-nonmetal transition in condensed media Phys. Rev. B 1978 17 2575 2581 10.1103/PhysRevB.17.2575
Edwards, P. P. & Sienko, M. J. Universality aspects of the metal-nonmetal transition in condensed media. Phys. Rev. B 17, 2575–2581 (1978).10.1103/PhysRevB.17.2575
48. Bardeen J Shockley W Deformation potentials and mobilities in non-polar crystals Phys. Rev. 1950 80 72 80 10.1103/PhysRev.80.72
Bardeen, J. & Shockley, W. Deformation potentials and mobilities in non-polar crystals. Phys. Rev. 80, 72–80 (1950).10.1103/PhysRev.80.72
49. Watts, F. & Wolstenholme, J. Surface Analysis by XPS and AES, Vol. 1 (John Wiley & Sons, 2015).
50. Li T Strong room-temperature ferroelectricity in strained SrTiO3 homoepitaxial film Adv. Mater. 2021 33 2008316 10.1002/adma.202008316
Li, T. et al. Strong room-temperature ferroelectricity in strained SrTiO3 homoepitaxial film. Adv. Mater. 33, 2008316 (2021).10.1002/adma.202008316
51. Lu W Yang P Song WD Chow GM Chen JS Control of oxygen octahedral rotations and physical properties in SrRuO3 films Phys. Rev. B Condens. Matter Mater. Phys. 2013 88 1 11 10.1103/PhysRevB.88.214115
Lu, W., Yang, P., Song, W. D., Chow, G. M. & Chen, J. S. Control of oxygen octahedral rotations and physical properties in SrRuO3 films. Phys. Rev. B Condens. Matter Mater. Phys. 88, 1–11 (2013).10.1103/PhysRevB.88.214115
52. Bjaalie L Janotti A Himmetoglu B Van De Walle CG Turning SrTiO3 into a Mott insulator Phys. Rev. B Condens. Matter Mater. Phys. 2014 90 195117 10.1103/PhysRevB.90.195117
Bjaalie, L., Janotti, A., Himmetoglu, B. & Van De Walle, C. G. Turning SrTiO3 into a Mott insulator. Phys. Rev. B Condens. Matter Mater. Phys. 90, 195117 (2014).10.1103/PhysRevB.90.195117
53. Khomskii, D. I. Transition Metal Compounds. Transition Metal Compounds10.1017/CBO9781139096782 (2014).
54. Chen Q Effect of La3+, Ag+ and Bi3+ doping on thermoelectric properties of SrTiO3: first-principles investigation Ceram. Int. 2022 48 13803 13816 10.1016/j.ceramint.2022.01.262
Chen, Q. et al. Effect of La3+, Ag+ and Bi3+ doping on thermoelectric properties of SrTiO3: first-principles investigation. Ceram. Int. 48, 13803–13816 (2022).10.1016/j.ceramint.2022.01.262
55. Azevedo, S. A., Laranjeira, J. A. S., Martins, N. F. & Sambrano, J. R. Ag doping effect on electronic and thermoelectric properties of SrTiO3 (0 0 1) surface. Comput. Mater. Sci. 227, 112274 (2023).
56. Wang Y Enhancement of Seebeck coefficient for SrO (SrTiO3)2 by Sm substitution: crystal symmetry restoration of distorted TiO6 octahedra Appl. Phys. Lett. 2007 91 242102 10.1063/1.2820447
Wang, Y. et al. Enhancement of Seebeck coefficient for SrO (SrTiO3)2 by Sm substitution: crystal symmetry restoration of distorted TiO6 octahedra. Appl. Phys. Lett. 91, 242102 (2007).10.1063/1.2820447
57. Zhang X Electronic quality factor for thermoelectrics Sci. Adv. 2020 6 6 11 10.1126/sciadv.abc0726
Zhang, X. et al. Electronic quality factor for thermoelectrics. Sci. Adv. 6, 6–11 (2020).10.1126/sciadv.abc0726
58. Ohtaki M Araki K Yamamoto K High thermoelectric performance of dually doped ZnO ceramics J. Electron. Mater. 2009 38 1234 1238 10.1007/s11664-009-0816-1
Ohtaki, M., Araki, K. & Yamamoto, K. High thermoelectric performance of dually doped ZnO ceramics. J. Electron. Mater. 38, 1234–1238 (2009).10.1007/s11664-009-0816-1
59. Ahmad A Thermoelectric performance enhancement of vanadium doped n-type In2O3 ceramics via carrier engineering and phonon suppression ACS Appl. Energy Mater. 2020 3 1552 1558 10.1021/acsaem.9b02011
Ahmad, A. et al. Thermoelectric performance enhancement of vanadium doped n-type In2O3 ceramics via carrier engineering and phonon suppression. ACS Appl. Energy Mater. 3, 1552–1558 (2020).10.1021/acsaem.9b02011
60. Bocher L CaMn1-x NbxO3 (x ≤ 0.08) perovskite-type phases as promising new high-temperature n-type thermoelectric materials Inorg. Chem. 2008 47 8077 8085 10.1021/ic800463s 18698764
Bocher, L. et al. CaMn1-x NbxO3 (x ≤ 0.08) perovskite-type phases as promising new high-temperature n-type thermoelectric materials. Inorg. Chem. 47, 8077–8085 (2008).18698764 10.1021/ic800463s
61. Liu H Enhanced thermoelectric properties of nonstoichiometric TiO1.76 with excellent mechanical properties induced by optimizing processing parameters Ceram. Int. 2018 44 19859 19865 10.1016/j.ceramint.2018.07.246
Liu, H. et al. Enhanced thermoelectric properties of nonstoichiometric TiO1.76 with excellent mechanical properties induced by optimizing processing parameters. Ceram. Int. 44, 19859–19865 (2018).10.1016/j.ceramint.2018.07.246
62. Tan X Synergistical enhancement of thermoelectric properties in n-type Bi2O2Se by carrier engineering and hierarchical microstructure Adv. Energy Mater. 2019 9 1 7 10.1002/aenm.201900354
Tan, X. et al. Synergistical enhancement of thermoelectric properties in n-type Bi2O2Se by carrier engineering and hierarchical microstructure. Adv. Energy Mater. 9, 1–7 (2019).10.1002/aenm.201900354
63. Lindermuth DJR An introduction to AC field Hall effect measurements J. Cryst. Growth 1976 36 29 35
Lindermuth, D. J. R. An introduction to AC field Hall effect measurements. J. Cryst. Growth 36, 29–35 (1976).
64. Lindemuth, J., Dodrill, B., Meyer, J. & Vurgaftman, I. Extraction of low mobility, low conductivity carriers from field. In CS-MAX 48–50 (IOP Publishing, 2002).
65. Cahill DG Analysis of heat flow in layered structures for time-domain thermoreflectance Rev. Sci. Instrum. 2004 75 5119 5122 10.1063/1.1819431
Cahill, D. G. Analysis of heat flow in layered structures for time-domain thermoreflectance. Rev. Sci. Instrum. 75, 5119–5122 (2004).10.1063/1.1819431
66. Jiang, P., Qian, X. & Yang, R. Tutorial: time-domain thermoreflectance (TDTR) for thermal property characterization of bulk and thin film materials. J. Appl. Phys. 124, 161103 (2018).
67. Park, K. & Bayram, C. Impact of dislocations on the thermal conductivity of gallium nitride studied by time-domain thermoreflectance. J. Appl. Phys. 126, 185103 (2019).
68. Zheng Q Thermal conductivity of GaN, 71GaN, and SiC from 150 K to 850 K Phys. Rev. Mater. 2019 3 014601 10.1103/PhysRevMaterials.3.014601
Zheng, Q. et al. Thermal conductivity of GaN, 71GaN, and SiC from 150 K to 850 K. Phys. Rev. Mater. 3, 014601 (2019).10.1103/PhysRevMaterials.3.014601
69. Braun JL Size effects on the thermal conductivity of amorphous silicon thin films Phys. Rev. B 2016 93 140201 10.1103/PhysRevB.93.140201
Braun, J. L. et al. Size effects on the thermal conductivity of amorphous silicon thin films. Phys. Rev. B 93, 140201 (2016).10.1103/PhysRevB.93.140201
70. Li C Remarkably weak anisotropy in thermal conductivity of two-dimensional hybrid perovskite butylammonium lead iodide crystals Nano Lett. 2021 21 3708 3714 10.1021/acs.nanolett.0c04550 33938755
Li, C. et al. Remarkably weak anisotropy in thermal conductivity of two-dimensional hybrid perovskite butylammonium lead iodide crystals. Nano Lett. 21, 3708–3714 (2021).33938755 10.1021/acs.nanolett.0c04550
71. Hammersley AP Svensson SO Hanfland M Fitch AN Häusermann D Two-dimensional detector software: from real detector to idealised image or two-theta scan High. Press. Res. 1996 14 235 248 10.1080/08957959608201408
Hammersley, A. P., Svensson, S. O., Hanfland, M., Fitch, A. N. & Häusermann, D. Two-dimensional detector software: from real detector to idealised image or two-theta scan. High. Press. Res. 14, 235–248 (1996).10.1080/08957959608201408
72. Qiu X Thompson JW Billinge SJL PDFgetX2: a GUI-driven program to obtain the pair distribution function from X-ray powder diffraction data J. Appl. Crystallogr. 2004 37 678 10.1107/S0021889804011744
Qiu, X., Thompson, J. W. & Billinge, S. J. L. PDFgetX2: a GUI-driven program to obtain the pair distribution function from X-ray powder diffraction data. J. Appl. Crystallogr. 37, 678 (2004).10.1107/S0021889804011744
73. McGreevy RL Reverse Monte Carlo modelling J. Phys. Condens. Matter 2001 13 R877 R913 10.1088/0953-8984/13/46/201
McGreevy, R. L. Reverse Monte Carlo modelling. J. Phys. Condens. Matter 13, R877–R913 (2001).10.1088/0953-8984/13/46/201
74. Li W Delineating complex ferroelectric domain structures via second harmonic generation spectral imaging J. Mater. 2023 9 395 402
Li, W. et al. Delineating complex ferroelectric domain structures via second harmonic generation spectral imaging. J. Mater. 9, 395–402 (2023).
75. Hÿtch MJ Snoeck E Kilaas R Quantitative measurement of displacement and strain fields from HREM micrographs Ultramicroscopy 1998 74 131 146 10.1016/S0304-3991(98)00035-7
Hÿtch, M. J., Snoeck, E. & Kilaas, R. Quantitative measurement of displacement and strain fields from HREM micrographs. Ultramicroscopy 74, 131–146 (1998).10.1016/S0304-3991(98)00035-7
76. Hüe F Calibration of projector lens distortions J. Electron Microsc. 2005 54 181 190
Hüe, F. et al. Calibration of projector lens distortions. J. Electron Microsc. 54, 181–190 (2005).
77. Hÿtch MJ Plamann T Imaging conditions for reliable measurement of displacement and strain in high-resolution electron microscopy Ultramicroscopy 2001 87 199 212 10.1016/S0304-3991(00)00099-1 11334167
Hÿtch, M. J. & Plamann, T. Imaging conditions for reliable measurement of displacement and strain in high-resolution electron microscopy. Ultramicroscopy 87, 199–212 (2001).11334167 10.1016/S0304-3991(00)00099-1
78. Deringer VL Tchougréeff AL Dronskowski R Crystal orbital Hamilton population (COHP) analysis as projected from plane-wave basis sets J. Phys. Chem. A 2011 115 5461 5466 10.1021/jp202489s 21548594
Deringer, V. L., Tchougréeff, A. L. & Dronskowski, R. Crystal orbital Hamilton population (COHP) analysis as projected from plane-wave basis sets. J. Phys. Chem. A 115, 5461–5466 (2011).21548594 10.1021/jp202489s
79. Kresse G Ab initio molecular dynamics for liquid metals J. Non Cryst. Solids 1995 192–193 222 229 10.1016/0022-3093(95)00355-X
Kresse, G. Ab initio molecular dynamics for liquid metals. J. Non Cryst. Solids 192–193, 222–229 (1995).10.1016/0022-3093(95)00355-X
80. Kresse G Furthmüller J Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set Phys. Rev. B Condens. Matter Mater. Phys. 1996 54 11169 11186 10.1103/PhysRevB.54.11169
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B Condens. Matter Mater. Phys. 54, 11169–11186 (1996).10.1103/PhysRevB.54.11169
81. Kresse G Hafner J Ab initio molecular dynamics for liquid metals Phys. Rev. B 1993 47 558 561 10.1103/PhysRevB.47.558
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561 (1993).10.1103/PhysRevB.47.558
82. Blöchl PE Projector augmented-wave method Phys. Rev. B 1994 50 17953 17979 10.1103/PhysRevB.50.17953
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).10.1103/PhysRevB.50.17953
83. Perdew JP Burke K Ernzerhof M Generalized gradient approximation made simple Phys. Rev. Lett. 1996 77 3865 3868 10.1103/PhysRevLett.77.3865 10062328
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).10062328 10.1103/PhysRevLett.77.3865
84. Zhang L Lin DY Wang H Car R Weinan E Active learning of uniformly accurate interatomic potentials for materials simulation Phys. Rev. Mater. 2019 3 023804 10.1103/PhysRevMaterials.3.023804
Zhang, L., Lin, D. Y., Wang, H., Car, R. & Weinan, E. Active learning of uniformly accurate interatomic potentials for materials simulation. Phys. Rev. Mater. 3, 023804 (2019).10.1103/PhysRevMaterials.3.023804
85. Zhang Y DP-GEN: a concurrent learning platform for the generation of reliable deep learning based potential energy models Comput. Phys. Commun. 2020 253 107206 10.1016/j.cpc.2020.107206
Zhang, Y. et al. DP-GEN: a concurrent learning platform for the generation of reliable deep learning based potential energy models. Comput. Phys. Commun. 253, 107206 (2020).10.1016/j.cpc.2020.107206
86. He R Ferroelastic twin-wall-mediated ferroelectriclike behavior and bulk photovoltaic effect in SrTiO3 Phys. Rev. Lett. 2024 132 1 7 10.1103/PhysRevLett.132.176801
He, R. et al. Ferroelastic twin-wall-mediated ferroelectriclike behavior and bulk photovoltaic effect in SrTiO3. Phys. Rev. Lett. 132, 1–7 (2024).10.1103/PhysRevLett.132.176801
87. He R Structural phase transitions in SrTiO3 from deep potential molecular dynamics Phys. Rev. B 2022 105 1 10 10.1103/PhysRevB.105.064104
He, R. et al. Structural phase transitions in SrTiO3 from deep potential molecular dynamics. Phys. Rev. B 105, 1–10 (2022).10.1103/PhysRevB.105.064104
88. Tadano T Tsuneyuki S Self-consistent phonon calculations of lattice dynamical properties in cubic SrTiO3 with first-principles anharmonic force constants Phys. Rev. B 2015 92 054301 10.1103/PhysRevB.92.054301
Tadano, T. & Tsuneyuki, S. Self-consistent phonon calculations of lattice dynamical properties in cubic SrTiO3 with first-principles anharmonic force constants. Phys. Rev. B 92, 054301 (2015).10.1103/PhysRevB.92.054301
89. Kresse G Furthmüller J Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set Comput. Mater. Sci. 1996 6 15 50 10.1016/0927-0256(96)00008-0
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).10.1016/0927-0256(96)00008-0
90. Kresse G Hafner J Ab initio molecular-dynamics simulation of the liquid-metalamorphous- semiconductor transition in germanium Phys. Rev. B 1994 49 14251 14269 10.1103/PhysRevB.49.14251
Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metalamorphous- semiconductor transition in germanium. Phys. Rev. B 49, 14251–14269 (1994).10.1103/PhysRevB.49.14251
91. Perdew JP Restoring the density-gradient expansion for exchange in solids and surfaces Phys. Rev. Lett. 2008 100 136406 10.1103/PhysRevLett.100.136406 18517979
Perdew, J. P. et al. Restoring the density-gradient expansion for exchange in solids and surfaces. Phys. Rev. Lett. 100, 136406 (2008).18517979 10.1103/PhysRevLett.100.136406
92. Han J Zhang L Car R Weinan E Deep potential: a general representation of a many-body potential energy surface Commun. Comput. Phys. 2018 23 629 639 10.4208/cicp.OA-2017-0213
Han, J., Zhang, L., Car, R. & Weinan, E. Deep potential: a general representation of a many-body potential energy surface. Commun. Comput. Phys. 23, 629–639 (2018).10.4208/cicp.OA-2017-0213
93. Souvatzis P Eriksson O Katsnelson MI Rudin SP Entropy driven stabilization of energetically unstable crystal structures explained from first principles theory Phys. Rev. Lett. 2008 100 095901 10.1103/PhysRevLett.100.095901 18352725
Souvatzis, P., Eriksson, O., Katsnelson, M. I. & Rudin, S. P. Entropy driven stabilization of energetically unstable crystal structures explained from first principles theory. Phys. Rev. Lett. 100, 095901 (2008).18352725 10.1103/PhysRevLett.100.095901
94. Van Roekeghem A Carrete J Mingo N Anomalous thermal conductivity and suppression of negative thermal expansion in ScF3 Phys. Rev. B 2016 94 020303 10.1103/PhysRevB.94.020303
Van Roekeghem, A., Carrete, J. & Mingo, N. Anomalous thermal conductivity and suppression of negative thermal expansion in ScF3. Phys. Rev. B 94, 020303 (2016).10.1103/PhysRevB.94.020303
95. Körmann, F., Ikeda, Y., Grabowski, B. & Sluiter, M. H. F. Phonon broadening in high entropy alloys. npj Comput. Mater.10.1038/s41524-017-0037-8 (2017).
96. Ikeda Y Carreras A Seko A Togo A Tanaka I Mode decomposition based on crystallographic symmetry in the band-unfolding method Phys. Rev. B 2017 95 024305 10.1103/PhysRevB.95.024305
Ikeda, Y., Carreras, A., Seko, A., Togo, A. & Tanaka, I. Mode decomposition based on crystallographic symmetry in the band-unfolding method. Phys. Rev. B 95, 024305 (2017).10.1103/PhysRevB.95.024305
97. Togo, A. First-principles phonon calculations with phonopy and phono3py. J. Phys. Soc. Jpn. 92, 012001 (2023).
98. Togo A Chaput L Tadano T Tanaka I Implementation strategies in phonopy and phono3py J. Phys. Condens. Matter 2023 35 353001 10.1088/1361-648X/acd831
Togo, A., Chaput, L., Tadano, T. & Tanaka, I. Implementation strategies in phonopy and phono3py. J. Phys. Condens. Matter 35, 353001 (2023).10.1088/1361-648X/acd831
99. Wright AJ Size disorder as a descriptor for predicting reduced thermal conductivity in medium- and high-entropy pyrochlore oxides Scr. Mater. 2020 181 76 81 10.1016/j.scriptamat.2020.02.011
Wright, A. J. et al. Size disorder as a descriptor for predicting reduced thermal conductivity in medium- and high-entropy pyrochlore oxides. Scr. Mater. 181, 76–81 (2020).10.1016/j.scriptamat.2020.02.011
100. He, X. et al. Hydride anion substitution boosts thermoelectric performance of polycrystalline SrTiO3 via simultaneous realization of reduced thermal conductivity and high electronic conductivity. Adv. Funct. Mater. 33, 2213144 (2023).
101. Li JB Broadening the temperature range for high thermoelectric performance of bulk polycrystalline strontium titanate by controlling the electronic transport properties J. Mater. Chem. C. 2018 6 7594 7603 10.1039/C8TC02130A
Li, J. B. et al. Broadening the temperature range for high thermoelectric performance of bulk polycrystalline strontium titanate by controlling the electronic transport properties. J. Mater. Chem. C. 6, 7594–7603 (2018).10.1039/C8TC02130A
102. Bakhshi H Sarraf-Mamoory R Yourdkhani A Abdelnabi AA Mozharivskyj Y Highly dense Sr0.95Sm0.0125Dy0.0125□0.025Ti0.90Nb0.10O3±δ/ZrO2 composite preparation directly through spark plasma sintering and its thermoelectric properties Dalton Trans. 2019 49 17 22 10.1039/C9DT03934A 31799569
Bakhshi, H., Sarraf-Mamoory, R., Yourdkhani, A., Abdelnabi, A. A. & Mozharivskyj, Y. Highly dense Sr0.95Sm0.0125Dy0.0125□0.025Ti0.90Nb0.10O3 ±δ/ZrO2 composite preparation directly through spark plasma sintering and its thermoelectric properties. Dalton Trans. 49, 17–22 (2019).31799569 10.1039/C9DT03934A
