Article
High Thermoelectric Performance ofAg9GaSe6 Enabled by Low Cutoff Frequencyof Acoustic Phonons
Siqi Lin, Wen Li, Shasha Li, ...,
Yidong Xu, Yue Chen, Yanzhong
Pei
HIGHLIGHTS
Ag9GaSe6 shows a promising
thermoelectric figure of merit
(zT � 1.5)
It has an extremely low lattice
thermal conductivity (kL� 0.15 W/m K)
Its low cutoff frequency of acoustic
phonons contributes to the low kL
The principle guides the design of
materials with expected kL
Acoustic phonons contribute most to lattice thermal conductivity (kL), due to their
high velocities. The cutoff frequency of acoustic phonons (um) indicates not only
the population but also the velocity and therefore guides the discovery of
semiconducting Ag9GaSe6 with an extremely low kL through a low um stemming
from both its large primitive cell and weakly bonded Ag. This opens new
possibilities for high-performance thermoelectrics, illustrated by the high
thermoelectric figure of merit (zT = 1.5) in this material.
Lin et al., Joule 1, 816–830
December 20, 2017 ª 2017 Elsevier Inc.
http://dx.doi.org/10.1016/j.joule.2017.09.006
Article
High Thermoelectric Performanceof Ag9GaSe6 Enabled by Low CutoffFrequency of Acoustic PhononsSiqi Lin,1 Wen Li,1 Shasha Li,2 Xinyue Zhang,1 Zhiwei Chen,1 Yidong Xu,1 Yue Chen,2
and Yanzhong Pei1,3,*
Context & Scale
Heat flow drives the charge
carriers in a material to flow along
the same direction, enabling
direct conversion from heat to
electricity, known as
thermoelectric technology. In
order to work efficiently, the
materials should be thermal
insulators to minimize heat loss.
Lattice vibrations, with long
wavelengths and high
propagation speeds (acoustic
SUMMARY
When the number of atoms in a primitive cell (N) increases, the first Brillouin
zone shrinks and folds back the high-frequency portion of acoustic vibrations
as optical ones with reduced velocities near zone boundaries. Soft bonds lead
to flattened phonon dispersion and thus a low sound velocity. This work shows
the cutoff frequency of acoustic phonons (um) as an important parameter to
include both effects on the lattice thermal conductivity (kL). A low um of
�0.5 THz in Ag9GaSe6 due to its large N and weakly bonded Ag atoms leads
this material to be one of the least thermally conductive dense solids (kL� 0.15 W/m K). This contributes to the peak thermoelectric figure of merit,
zT, as high as 1.5. The principle used here not only leads to the exploration of
Ag9GaSe6 as a promising thermoelectric but also enables a measurable param-
eter for designing materials with expected kL.
phonons), dominate heat
conduction in most thermoelectric
materials.
The cutoff frequency of acoustic
phonons (um), characterizing both
the population and the velocity, is
therefore a key measure of lattice
thermal conductivity (kL). Here, an
extremely low um in Ag9GaSe6 is
realized due to its large primitive
cell and weak chemical bonds for
both a small population and
velocity of acoustic phonons. The
resulting extremely low kL
contributes to the superior
thermoelectric figure of merit, zT,
of 1.5. The principle here guides
the design of materials with
expected kL for many potential
applications of heat, including
thermoelectrics.
INTRODUCTION
Thermoelectrics, directly converting heat into electricity, has long been considered
as a potential technology for solving the energy crisis.1 The performance of a
thermoelectric material is determined by the figure of merit zT = sS2T/(kE+kL), where
s, S, T, kE, and kL are the electrical conductivity, the Seebeck coefficient, the absolute
temperature, and the electronic and the lattice components of the thermal conduc-
tivity, respectively.
Because S, s, and kE strongly couple with each other, it is difficult to achieve a high zTby
a simple improvement of one of these parameters. Decoupling these parameters for
enhancing zT via band convergence or nestification has been demonstrated to be
successful in various materials such as IV-VI compounds,2–8 Mg2Si,9 half-Heusler,10,11
and Te.12,13 Alternatively, reducing the only independent material property, lattice
thermal conductivity kL, has also been demonstrated to effectively enhance zT. Typical
approaches include point defects,14–20 dislocations,21–23 and nanostructures.24–27
Recently, numerous efforts have been made to develop new thermoelectric
materials with superior performance. A variety of new materials with high zT, such
as b-Zn4Sb3,28 Cu2Se,
29,30 SnSe,31 and MgAgSb32–34 have been reported to be
promising for thermoelectric applications. These high zT materials have intrinsically
low lattice thermal conductivity, due to various mechanisms, including crystal
structure complexity,28 liquid-like ions,29 and lattice anharmonicity.35 This works
indicates the importance of an intrinsic low lattice thermal conductivity for exploring
new thermoelectrics.
816 Joule 1, 816–830, December 20, 2017 ª 2017 Elsevier Inc.
1Interdisciplinary Materials Research Center,School of Materials Science and Engineering,Tongji University, 4800 Caoan Road, Shanghai201804, China
2Department of Mechanical Engineering, TheUniversity of Hong Kong, Pokfulam Road, HongKong SAR, China
3Lead Contact
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.joule.2017.09.006
Intrinsically low lattice thermal conductivity usually happens in materials with a large
primitive cell. This is fundamentally associated with the reduced population of
acoustic phonons, which have a high velocity and therefore contribute dominantly
to the lattice thermal conductivity. In more detail, when the number of atoms in a
primitive cell increases, the lattice parameter increases, resulting in a reduced
volume of the first Brillouin zone. This leads the high-frequency portion of acoustic
vibration modes to be folded back into the first Brillouin zone as optical modes,
but with significantly reduced velocities near zone boundaries. Optical phonons
usually have extremely low velocities, leading to a negligible contribution to the lat-
tice thermal conductivity.36,37 In the case of a large primitive cell, the population of
acoustic phonons that dominate the lattice thermal conductivity is significantly
reduced.
Another key parameter determining the lattice thermal conductivity is the propaga-
tion of phonons (i.e., sound velocity, v). When the medium varies from a gas to a
liquid and then to a solid, sound waves travel faster and faster, roughly indicating
that a weakly bonded medium reduces the velocity. In a crystalline solid, weakly
bonded heavy constituent elements usually lead to a flattening in phonon disper-
sion. This essentially results in a low lattice thermal conductivity.
Although the above strategies are indeed effective for exploring materials
with a low lattice thermal conductivity (kL), the individual effectiveness on kL is
by definition still not well quantified. Since both a large primitive cell and soft
bonds lead to a small frequency of acoustic phonons at the boundary of the
first Brillouin zone, this work uses such a frequency as an important parameter
determining the kL of a material. This frequency is exactly the cutoff frequency of
acoustic phonons (um), which can be determined by measuring the sound velocity
(v) and the crystal structure parameters (the number of atoms in the primitive cell
[N] and the atomic volume [V ]) via um = 2pv � ðNVÞ�1=3 according to Debye
approximation.
Guided by the above concept of a low cutoff frequency for low lattice thermal
conductivities, this work focuses on a semiconducting material, Ag9GaSe6, for its
potential in thermoelectrics. The primitive cell includes a large number of atoms,
more than half of which are weakly bonded Ag+ cations. Therefore, Ag9GaSe6 is
assumed to have a low lattice thermal conductivity due to its low um for potential
thermoelectric applications. It should be noted that the weakly bonded Ag+ cations
further induce an order-disorder phase transition from b-Ag9GaSe6 to a-Ag9GaSe6at 281 K.38–40
The thermoelectric transport properties of Ag9GaSe6 above room temperature are
investigated in this work. While most of the semiconductors, including thermoelec-
tric materials, have a cutoff frequency of 1–3 THz for acoustic phonons, Ag9GaSe6shows an extremely low cutoff frequency (only �0.5 THz) of acoustic phonons.
This material successfully shows a kL as low as�0.15W/m K in the entire temperature
range, being one of the least thermally conductive crystalline solids. This work may
shed light on the origin of the extremely low lattice thermal conductivities observed
in other high-performance thermoelectrics41–48 as well. Although the thermoelectric
power factor of Ag9GaSe6 is much lower than those of conventional thermoelectrics,
a peak thermoelectric figure of merit, zT � 1.5, can be achieved in Ag9GaSe6 alloys,
of which the transport properties can be well understood by a single parabolic band
model. This work opens new possibilities of thermoelectric advancements through a
low cutoff frequency of acoustic phonons.
Joule 1, 816–830, December 20, 2017 817
Figure 1. Phase Characterization
(A–D) Crystal structure of Ag9GaSe6 in the high temperature (>281 K) phase (A). Powder X-ray
diffraction patterns (B), the lattice parameter (C), and room temperature Hall carrier concentration
(D) for Ag9Ga(Se1�xTex)6, indicating formation of a solid solution.
RESULTS AND DISCUSSION
Ag9GaSe6 shows a phase transition between b-Ag9GaSe6 (space group P213) and
a-Ag9GaSe6 (space group F43m) at 281 K.38–40 The cubic crystal structure of the
low-temperature b-Ag9GaSe6 is formed by a [GaSe4]5� and Se2� anion framework
and ordered Ag+ fully occupying its lattice sites. Above 281 K, the anion framework
becomes more symmetrical in an F43m lattice, while Ag atoms get highly disordered
with an occupancy of only 12.5% and 25%, respectively, on the 96i and 48h sites
(Figure 1A). It is interesting to note that the order-disorder transition does not induce
a change in the lattice parameter. Moreover, the disordered phase stabilizes at room
temperature or higher. Compared with previously reported materials with complex
crystal structures involving phase transitions,29,34,35,41,47,48 these features can be
advantages for thermoelectric applications.
The powder X-ray diffraction patterns (XRD) for the samples are shown in Figure 1B.
All the peaks can be indexed to the corresponding F43m structure.With increasing x,
the lattice parameter (Figure 1C) increases linearly while the carrier concentration
decreases (Figure 1D). The lattice expansion can be understood by the larger
size of Te compared with that of Se. A broad range of Hall carrier concentration
(3 � 10 3 1018 cm�3) at room temperature is achieved (Figure 1D). Other dopants
such as Cr, Cd, Zn, Ge, and S were also used to tune the carrier concentration, yet
they all seem to be not as effective (Table S1).
Scanning electron microscopy (SEM) observations and energy dispersive
spectrometer (EDS) composition mapping of Ag, Se, and Ga for Ag9GaSe6 (Figure 2)
and Ag9Ga(Se0.9Te0.1)6 (Figure S1) were carried out to further confirm the purity and
homogeneity. The optical measurements49 enabled estimation of the band gap
818 Joule 1, 816–830, December 20, 2017
Figure 2. Microstructures of Ag9GaSe6(A–D) SEM image (A) and the corresponding EDS composition mapping of Ag (B), Se (C), Ga (D) for
Ag9GaSe6.
(�0.56 eV) for Ag9GaSe6. It is shown that the band gap of Ag9Ga(Se1�xTex)6decreased from 0.56 eV for x = 0 to 0.46 eV for x = 0.17 (Figure 3). The XRD results
and the optical and transport properties for Ag9GaSe6 with other dopants
mentioned above are included in Figures S2 and S3.
Ab initio molecular dynamics simulations were performed to understand the
atomic vibrations at elevated temperatures for the low-temperature phase (Fig-
ure 4A). It is seen that all atoms vibrate in the vicinity of their equilibrium positions
at 300 K (Figure 4B), whereas the Ag atoms exhibit a diffusive behavior at 500 K
(Figure 4C). This indicates that Ag atoms are loosely bonded to the neighboring
atoms. The simulation results should rationalize the experimentally observed Ag
partial occupancy in the high-temperature phase. It should be noted that, due
to the partial occupancy of Ag in the high-temperature phase, the simulation
can be carried out on the low-temperature phase only. However, it is reasonable
to assume that Ag atoms get more diffusive in the high-temperature phase
because of the existence of unoccupied vacancies. Further considering the un-
changed lattice parameters for both low- and high-temperature cubic phases,
the main difference relies on the occupancy of Ag atoms. Therefore, it is believed
that Ag atoms are also weakly bonded (if not weaker) in the high-temperature
phase.
Further electron density distribution and Bader charges are carried out to confirm
the soft bonding in Ag9GaSe6. Directional electron density distribution between
Ga and Se atoms is seen from Figure 5, indicating strong covalent bonds. However,
the electron density distribution surrounding the Ag atoms is nearly spherical,
suggesting weak chemical bonds to neighboring atoms. Bader charge analysis
shows that each Ag and Ga atom loses 0.26 and 0.97 electrons, respectively, while
each Se atom receives 0.56 electrons. This again indicates that Ag atoms are weakly
bonded due to ionic interactions. Such bonding characteristics provide an insight in
Joule 1, 816–830, December 20, 2017 819
Figure 3. Optical Absorption
The normalized optical absorption versus photon energy at room temperature for
Ag9Ga(Se1�xTex)6.
understanding why the lattice dynamics shows a diffusive behavior of Ag atoms
(Figure 4).
It is known that soft bonding usually leads to a flattened phonon dispersion,
which is shown in Figure 6A. By projecting the phonon density of states (DOS) of
Ag9GaSe6 onto different atomic species (Figure 6B), it is found that the majority
of the low-frequency phonon modes below 3 THz are contributed by the Ag atoms.
This is consistent with the feature that Ag atoms are weakly bonded. Therefore, it
is seen that most of the optical phonon modes are nearly non-dispersive,
ensuring a negligible contribution to the lattice thermal conductivity.37
Further, due to the soft bonding, both the transverse (vt = 1,200 m/s) and longitu-
dinal (vl = 2,900 m/s) sound velocities calculated from phonon dispersions are
extremely low for Ag9GaSe6, which is well confirmed by the measurements
(vt = 1,167 m/s, vl = 2,917 m/s) at room temperature. Besides the soft bonds
emphasized in this work, other mechanisms50,51 are found to be effective for a
low sound velocity.
Because there are many atoms (high N) in the crystal cell of Ag9GaSe6, the small
volume of the first Brillouin zone and the softly bonded Ag atoms (low v) lead the
cutoff frequency of acoustic phonons to be as low as �0.5 THz (Figure 5). This cutoff
frequency of acoustic phonons is significantly lower than those of typical semi-
conductors, including thermoelectrics. The room temperature lattice thermal con-
ductivity as a function of cutoff frequency for acoustic phonons is shown in Figure 7
for semiconductors with potential as thermoelectrics (more details are given in
Table S2). It is seen that a low cutoff frequency of acoustic phonons is indeed helpful
for a low lattice thermal conductivity.
The theoretical basis relies on the fact that the cutoff frequency of acoustic
phonons, um = 2pv � ðNVÞ�1=3, defines the upper limit of the integral of the lattice
820 Joule 1, 816–830, December 20, 2017
Figure 4. Molecular Dynamics Simulations
(A–C) Crystal structure of Ag9GaSe6 in the low-temperature phase projected onto the (100) plane
(A) and the corresponding atomic trajectories at (B) 300 K and (C) 500 K, with Ag in blue, Ga in red,
and Se in green.
thermal conductivity. Equation 1 shows the expression of kL for Umklapp scat-
tering, according to a Debye phonon dispersion and exclusion of optical
phonons:36,37
kL =
ffiffiffiffiffiffiffiffi6p23
p
4p2
R
kBMV
2=3Z um
0
v2
g2T
0B@� Zu
kBT
�2
,eðZu=kBTÞ
ðeðZu=kBTÞ � 1Þ2
1CAdu; (Equation 1)
Joule 1, 816–830, December 20, 2017 821
Figure 5. Electron Density Distribution
(A and B) Electron density isosurface of 0.4 e/A3 in Ag9GaSe6 (A) and the electron density
distribution in an atomic plane crossing the Se-Ga-Se bonds (B).
where R is the gas constant, kB is the Boltzmann constant, M is the average atomic
mass, v is the average sound speed, g is the Gruneisen parameter, and �h is the
reduced Planck constant. It should be noted that a more accurate determination
on phonon dispersion such as by density functional theory (DFT) calculations50,51
would lead to a more accurate estimation of kL.
Therefore, use of the number of atoms in the primitive cell (N) for measuring the
crystal complexity and the sound velocity (v) for characterizing the bonding stiffness
in this work, a product of 2pvðNVÞ�1=3, which is essentially the cutoff frequency of
acoustic phonons according to Debye approximation, successfully includes both
effects in the lattice thermal conductivity. The cutoff frequency of acoustic phonons
um physically means the frequency of acoustic phonons at the boundary of the
Brillouin zone. In the case of a given sound velocity, the volume of the first Brillouin
zone is proportional to um.3 Therefore, this work reveals a high lattice thermal
conductivity when the volume of the first Brillouin zone is large, at a given sound
velocity (Figure 7).
According to the above discussion, a um as low as �0.5 THz in Ag9GaSe6 (Figure 6)
would lead to an expectation of extremely low lattice thermal conductivity. The
temperature-dependent total thermal conductivity (k) and its lattice component
(kL) for Ag9Ga(Se1�xTex)6 (0 % x % 0.17) are shown in Figure 8. It is found that
k for all samples is as low as 0.5W/m K over the entire temperature range. The lattice
thermal conductivity (kL) is estimated by subtracting the electronic thermal
conductivity (ke) via the Wiedemann-Franz law (ke = LT/r) from the total thermal
conductivity, where L is the Lorenz factor determined by the single parabolic
band (SPB) model with acoustic scattering (Figure S4). The observed kL of
�0.15 W m�1 K�1, over the entire temperature range, is indeed one of the
lowest among known thermoelectrics. The much lower lattice thermal conductivity
in Ag9GaSe6 (kL � 0.15 W/m K) than in Ga2Se3 with high concentration intrinsic
vacancies (kL � 0.6 W/m K52) or in GaSe (kL � 2.1 and�16 W/m K along and perpen-
dicular to the c axis at room temperature, respectively53) suggests that the loosely
bonded Ag atoms contribute to the low lattice thermal conductivity of Ag9GaSe6.
It should be noted that the total thermal conductivity (k) is as low as�0.4 W/m K and
largely (�50%) comes from its electronic component; a combination of multiple
measurement uncertainties can mathematically affect the absolute values of kL.
822 Joule 1, 816–830, December 20, 2017
Figure 6. Phonon Dispersion
(A and B) Calculated phonon dispersions (A) and the projected phonon density of states (B) for
Ag9GaSe6.
A statistical measurement on kL yields an average of 0.15 W/m K for all samples with
an SD of 0.058W/m K (Figure 8). Nevertheless, we are able to draw a safe conclusion
that an extremely low kL is achieved in Ag9GaSe6 with a low um, which supports our
main conclusion. In addition, the measured kL very weakly depends on temperature,
which is a strong indication of Ag9GaSe6 as a phonon glass by nature. This
may explain that further introduction of atomic scale defects does not enable a
clear additional reduction in kL (Figure 8), which can easily fall within the range of
its measurement uncertainty, even if additional reduction indeed exists.
The measured room temperature longitudinal (vl), transverse (vt), and mean (v)
sound velocities for all the Ag9Ga(Se1�xTex)6 samples are listed in Table 1. The
change in v does not exceed 5%, which is within the measurement uncertainty.
Importantly, it is found that the mean sound velocity v � 1,350 m/s is one of the
lowest among known thermoelectric semiconductors.41 Table 1 also lists the bulk
(B) and shear (G) modulus, Gruneisen parameter, and Debye temperature (QD)
calculated from the sound velocities.54 The Gruneisen parameter, characterizing
the anharmonicity of the lattice vibrations, is high (�2.6), which is presumably
due to the loosely bonded Ag atoms. The Debye temperature is estimated to be
�143 K with an error of <3% for Ag9Ga(Se1�xTex)6. All these features support the
low lattice thermal conductivity observed in this work.
Providing a low kL ensured by the low um in Ag9GaSe6, a generic guideline for
realizing high zT will be a synergy with electronic performance enhancement
Joule 1, 816–830, December 20, 2017 823
Figure 7. Survey of um versus kLRoom temperature lattice thermal
conductivity versus the cutoff
frequency of acoustic phonons for
semiconductors.
through first optimizing the carrier concentration and then manipulating the band
structures at an optimized carrier concentration. This gives a band structure
calculation for the low-temperature phase (Figure 9) to reveal the electronic trans-
port properties of Ag9GaSe6-based materials. The conduction band maximal was
found at the G point, therefore, the band degeneracy (Nv) in n-type is 1. This is
significantly lower than that of conventional thermoelectrics, which usually have Nv
of 4 or higher.6,7,11,55,56 For this reason, all the samples here showing an n-type
conduction are not expected to have a high thermoelectric power factor. It should
be noted that the optical measurement suggests a direct band gap in Ag9GaSe6,
which is consistent with the DFT calculations.
The modified Becke-Johnson (mBJ) exchange potential is used in our DFT calcula-
tions. The mBJ exchange potential can significantly improve the theoretical band
gap, with a similar accuracy to the hybrid functional or the GW method. As a result,
the DFT band gap obtained is in good agreement with optimal measurements. In
addition, the calculated band structure shows a single band conduction behavior
in n-type, which is confirmed by the transport property measurements. Therefore,
the currently calculated band structure is reasonably reliable for evaluating the trans-
port properties in the absence of strong bipolar conduction as is the case in the cur-
rent study.
Substituting Se by Te in Ag9GaSe6 enables effective tuning of the carrier concen-
trations in a broad range, enabling a reliable evaluation of the thermoelectric
transport properties. Hall mobility (mH) decreases with increasing temperature via
mH � T�1.5 (Figure 10A) for all the samples in the entire temperature range, indi-
cating an overall dominant charge carrier scattering by acoustic phonons.57,58
This enables a SPB model59 to understand the electronic transport properties.
Assuming the conduction band is isotropic, the estimated density-of-state effec-
tive mass (m*) and deformation potential coefficient60 (Edef) are shown in Fig-
ure 10B. Both Edef and m* show a weak dependence on temperature and carrier
concentration. The SPB model-predicted Seebeck coefficient (Figure 10C) and
Hall mobility (Figure 10D), as a function of the Hall carrier concentration, agree
well with the measurements at different temperatures. All these indicate a rigid
band behavior for this material. Moreover, Edef of �24 eV and m* of 0.16 me
(me is the free electron mass) are within the typical range for thermoelectric
semiconductors.61–64
824 Joule 1, 816–830, December 20, 2017
Figure 8. Thermal Properties
Temperature-dependent total
thermal conductivity and lattice
thermal conductivity for
Ag9Ga(Se1�xTex)6, compared with
the lattice thermal conductivity of
Ga2Se3 with intrinsic vacancies.52
The temperature-dependent Seebeck coefficient, resistivity, and zT are shown in Fig-
ure 11 for Ag9Ga(Se1�xTex)6. A negative Seebeck coefficient indicates an n-type
conduction for all the samples, which is consistent with the Hall measurements.
The majority of the samples obtained in this work show a typical degenerate semi-
conducting behavior, meaning a continuous increase in the Seebeck coefficient
and resistivity as the temperature rises. The decrease in r and S at high temperatures
for low carrier concentration samples can be ascribed to the existence of minority
carriers.
A peak zT of �1.5 is achieved at 850 K, which mainly relies on the extremely low
lattice thermal conductivity due to the low cutoff frequency of acoustic phonons.
Using an average lattice thermal conductivity of 0.15 W m�1 K�1 (Figure 8), the
SPB model enables a prediction on zT versus Hall carrier concentration, which is
shown in Figure 11D. It is found that a peak zT as high as 1.6 can be expected
when the Hall carrier concentration is reduced to about 3 3 1018 cm�3 at 800 K.
This suggests that Ag9GaSe6 is a promising thermoelectric material.
Note here that zT in n-type Ag9GaSe6 seems to have not enough scope for further
enhancement through engineering the band, due to the simplicity of the conduction
band (Figure 9). However, an even higher zT could be expected through further band
engineering, probably in p-type conduction, because of the multiple valence band
structure (Figure 9). This can be the playground for this compound as a thermoelec-
tric material. The possibility of engineering valence bands as desired can be realized
in principle by the richness in the composition of the IB9IIIVI6 compound family,
Table 1. Elastic Properties
Composition vt (m/s) vl (m/s) V (m/s) B (GPa) G (GPa) g QD (K)
x = 0 1,167 2,917 1,322 48.8 10.0 2.7 141
x = 0.05 1,222 2,916 1,382 47.5 10.9 2.6 147
x = 0.08 1195 2,953 1,353 49.8 10.4 2.6 144
x = 0.10 1177 2,948 1,333 49.9 10.1 2.7 142
x = 0.12 1,195 2,975 1,353 50.7 10.4 2.7 144
x = 0.15 1,195 2,941 1,352 49.2 10.4 2.6 143
x = 0.17 1,197 2,955 1,355 49.8 10.5 2.6 143
Room temperature transverse (vt), longitudinal (vl), and mean (v) sound velocities, bulk (B) and shear (G)
modules, Gruneisen parameter (g) as well as Debye temperature (QD) for Ag9Ga(Se1�xTex)6.
Joule 1, 816–830, December 20, 2017 825
Figure 9. Band Structure
DFT band structure of Ag9GaSe6 in P213 phase.
where IB, III, and VI can be Ag/Cu, Al/Ga/In, and S/Se/Te, respectively. Achieving
possibly higher zT is further enabled by the fact that all these materials have similarly
complex crystal structures and partial occupancy on the IB sites, thus all are likely to
show an ultra-low intrinsic lattice thermal conductivity. Furthermore, a Cu-based
analog or heavy substitution of Ag by Cu can be an effective approach to address
the issue due to the scarcity of Ag and to reduce the cost of the material.
This work shows the cutoff frequency of acoustic phonons, as an important
parameter integrating effects on the lattice thermal conductivity due to crystal struc-
ture complexity and soft bonding. This enables not only effective guidance for
Figure 10. Electronic Transport Properties
(A–D) Temperature-dependent Hall mobility, mH (A), density-of-state effective mass, m* and
deformation potential coefficient, Edef (B), Hall carrier concentration-dependent Seebeck
coefficient (C), and Hall mobility (D) at 300, 500, and 800 K for Ag9GaSe6. The solid curves in (C) and
(D) show the SPB model predictions.
826 Joule 1, 816–830, December 20, 2017
Figure 11. Thermoelectric Transport Properties
(A–C) Temperature-dependent Seebeck coefficient (A), resistivity (B), and figure of merit, zT (C) for
Ag9Ga(Se1�xTex)6.
(D) The model-predicted zT versus carrier concentration at 300, 500, and 800 K for Ag9GaSe6.
designing materials with expected lattice thermal conductivity but also effective
exploration of new thermoelectrics such as Ag9GaSe6 of zT � 1.5. The scope for
follow-up studies on IB9IIIVI6 compounds relies on the richness in the composition,
which allows sufficient degrees of freedom for electronic structure manipulation for
performance enhancements.
EXPERIMENTAL PROCEDURES
Polycrystalline Ag9Ga(Se1�xTex)6 (x % 0.1) and Ag9Ga1�xMx(Se1�ySy)6 (M = Cr, Cd,
Zn, Ge; x% 0.06; y% 0.10) samples were synthesized by melting the stoichiometric
amount of high-purity elements (>99.99%) at 1,227 K for 6 hr, quenching in cold
water, and annealing at 900 K for 3 days. The resulting ingots were hand ground
into fine powder for X-ray diffraction (XRD) and hot press. Pellet samples were
obtained by an induction heating hot press system65 at 900 K for 1 hr under a uniaxial
pressure of�60 MPa. The dense samples obtained (>98% of the theoretical density)
were about 12 mm in diameter and �1.5 mm in thickness.
The electrical transport properties, including resistivity, Seebeck coefficient, and
Hall coefficient, were measured simultaneously on the pellet samples under helium.
The Seebeck coefficient was obtained from the slope of the thermopower versus
temperature difference within 0–5 K.66 The resistivity and Hall coefficient (RH) were
measured using the van der Pauw technique under a reversible magnetic field of
1.5 T. Thermal diffusivity (D) was measured using the laser flash technique with the
Netzsch LFA457 system. The thermal conductivity was calculated via k = dCpD,
where d is the density measured using the mass and geometric volume of the pellet
and Cp is the heat capacity determined by the Dulong-Petit limit. All the transport
Joule 1, 816–830, December 20, 2017 827
property measurements were carried out in the temperature range of 300–850 K.
Hysteresis on the transport properties was initially observed but disappeared
after a few thermal cycles or annealing. The measurement uncertainty for each
transport property (S, r, and k) is about 5%. The microstructure was characterized
by SEM equipped with EDS. The sound velocity was measured using an ultrasonic
pulse receiver (Olympus-NDT) equipped with an oscilloscope (Keysight). Optical
diffusive reflectance was measured by infrared Fourier transform spectroscopy
(Bruker Tensor II).
Phonon dispersions and projected phonon DOS of Ag9GaSe6 were calculated using
the small displacement method with a 2 3 2 3 2 supercell containing 512 atoms
with Phonopy.67 The force constants were calculated based on DFT using VASP68
with an energy tolerance of 10�8 eV. The electron-ion interactions were treated
using the projector augmented wave method.69 The exchange-correlation interac-
tions were taken into account with the Perdew-Burke-Ernzerhof functional.70 The
valence electronic states were expanded in plane-wave basis sets with an energy
cutoff of 400 eV. The Brillouin zone integrations were performed with a 1 3 1 3 1
Monkhorst-Pack grid71 for the supercell containing 512 atoms. A 23 23 2 supercell
was also used for ab initiomolecular dynamics simulations in the canonical ensemble
(NVT) at 300 K and 500 K. Ab initiomolecular dynamics simulations were performed
for 4 ps with a time step of 1 fs; atomic trajectories were collected from 0.5 ps.
The mBJ method72 was applied to obtain a more realistic electronic band gap
from DFT.
SUPPLEMENTAL INFORMATION
Supplemental Information includes four figures and two tables and can be found
with this article online at http://dx.doi.org/10.1016/j.joule.2017.09.006.
AUTHOR CONTRIBUTIONS
Conceptualization, Y.P.; DFT Calculations, S. Li and Y.C.; Experiments, S. Lin, W.L.,
and Y.X.; Discussions, S. Lin, W.L., S. Li, X.Z., Z.C., Y.X., Y.C., and Y.P.; Writing, S. Lin,
W.L., and Y.P.
ACKNOWLEDGMENTS
This work is supported by the National Natural Science Foundation of China (grant
nos. 51422208, 11474219, 51772215), the National Recruitment Program of Global
Youth Experts (1000 Plan). S. Li and Y.C. are grateful for financial support from the
Early Career Scheme of RGC under Project No. 27202516 and the research
computing facilities offered by ITS, HKU.
Received: April 25, 2017
Revised: June 30, 2017
Accepted: September 8, 2017
Published: October 4, 2017
REFERENCES
1. Bell, L.E. (2008). Cooling, heating, generatingpower, and recovering waste heat withthermoelectric systems. Science 321,1457–1461.
2. Pei, Y., Shi, X., LaLonde, A., Wang, H.,Chen, L., and Snyder, G.J. (2011).Convergence of electronic bands for high
828 Joule 1, 816–830, December 20, 2017
performance bulk thermoelectrics. Nature473, 66–69.
3. Pei, Y., Wang, H., Gibbs, Z.M., LaLonde, A.D.,and Snyder, G.J. (2012). Thermopowerenhancement in Pb1-xMnxTe alloys and itseffect on thermoelectric efficiency. NPG AsiaMater. 4, e28.
4. Pei, Y., LaLonde, A.D., Heinz, N.A., Shi, X.,Iwanaga, S., Wang, H., Chen, L., and Snyder,G.J. (2011). Stabilizing the optimal carrierconcentration for high thermoelectricefficiency. Adv. Mater. 23, 5674–5678.
5. Jian, Z., Chen, Z., Li, W., Yang, J., Zhang, W.,and Pei, Y. (2015). Significant band engineering
effect of YbTe for high performancethermoelectric PbTe. J. Mater. Chem. C 3,12410–12417.
6. Li, W., Zheng, L., Ge, B., Lin, S., Zhang, X.,Chen, Z., Chang, Y., and Pei, Y. (2017).Promoting SnTe as an eco-friendly solution forp-PbTe thermoelectric via band convergenceand interstitial defects. Adv. Mater 29, 1605887.
7. Li, J., Chen, Z., Zhang, X., Sun, Y., Yang, J., andPei, Y. (2017). Electronic origin of the highthermoelectric performance of GeTe amongthe p-type group IV monotellurides. NPG AsiaMater. 9, e353.
8. Li, W., Chen, Z., Lin, S., Chang, Y., Ge, B., Chen,Y., and Pei, Y. (2015). Band and scatteringtuning for high performance thermoelectricSn1-xMnxTe alloys. J. Materiomics 1, 307–315.
9. Liu, W., Tan, X., Yin, K., Liu, H., Tang, X., Shi, J.,Zhang, Q., and Uher, C. (2012). Convergenceof conduction bands as a means ofenhancing thermoelectric performance ofn-type Mg2Si1-xSnx solid solutions. Phys. Rev.Lett. 108, 166601.
10. Fu, C., Bai, S., Liu, Y., Tang, Y., Chen, L., Zhao,X., and Zhu, T. (2015). Realizing high figure ofmerit in heavy-band p-type half-Heuslerthermoelectric materials. Nat. Commun. 6,8144.
11. Fu, C.G., Zhu, T.J., Pei, Y.Z., Xie, H.H., Wang,H., Snyder, G.J., Liu, Y., Liu, Y.T., and Zhao, X.B.(2014). High band degeneracy contributes tohigh thermoelectric performance in p-typehalf-Heusler compounds. Adv. Energy Mater.4, 1400600.
12. Lin, S.Q., Li, W., Chen, Z.W., Shen, J.W., Ge,B.H., and Pei, Y.Z. (2016). Tellurium as ahigh-performance elemental thermoelectric.Nat. Commun. 7, 10287.
13. Lin, S.Q., Li, W., Zhang, X.Y., Li, J., Chen, Z.W.,and Pei, Y.Z. (2017). Sb induces both dopingand precipitation for improving thethermoelectric performance of elemental Te.Inorg. Chem. Front 4, 1066–1072.
14. Vining, C.B., and Rowe, D.M. (1995). SiliconGermanium (CRC Press), pp. 329–337.
15. Schemer, H., and Scherrer, S. (1995). BismuthTelluride, Antimony Telluride and Their Solid(CRC Press), pp. 27–41.
16. Shen, J.W., Zhang, X.Y., Chen, Z.W., Lin, S.Q.,Li, J., Li, W., Li, S.S., Chen, Y., and Pei, Y.Z.(2017). Substitutional defects enhancingthermoelectric CuGaTe2. J. Mater. Chem. A 5,5314–5320.
17. Shen, J.W., Zhang, X.Y., Lin, S.Q., Li, J., Chen,Z.W., Li, W., and Pei, Y.Z. (2016). Vacancyscattering for enhancing the thermoelectricperformance of CuGaTe2 solid solutions.J. Mater. Chem. A 4, 15464–15470.
18. Pei, Y.Z., Zheng, L.L., Li, W., Lin, S.Q., Chen,Z.W., Wang, Y.Y., Xu, X.F., Yu, H.L., Chen, Y.,and Ge, B.H. (2016). Interstitial point defectscattering contributing to high thermoelectricperformance in SnTe. Adv. Electron. Mater. 2,1600019.
19. May, A.F., Fleurial, J.P., and Snyder, G.J. (2008).Thermoelectric performance of lanthanumtelluride produced via mechanical alloying.Phys. Rev. B 78, 125205.
20. Li, W., Lin, S.Q., Zhang, X.Y., Chen, Z.W., Xu,X.F., and Pei, Y.Z. (2016). Thermoelectricproperties of Cu2SnSe4 with intrinsic vacancy.Chem. Mater. 28, 6227–6232.
21. Kim, S.I., Lee, K.H., Mun, H.A., Kim, H.S.,Hwang, S.W., Roh, J.W., Yang, D.J., Shin, W.H.,Li, X.S., and Lee, Y.H. (2015). Dense dislocationarrays embedded in grain boundaries for high-performance bulk thermoelectrics. Science348, 109–114.
22. Chen, Z., Ge, B., Li, W., Lin, S., Shen, J., Chang,Y., Hanus, R., Snyder, G.J., and Pei, Y. (2017).Vacancy-induced dislocations within grains forhigh-performance PbSe thermoelectrics. Nat.Commun. 8, 13828.
23. Chen, Z., Jian, Z., Li, W., Chang, Y., Ge, B.,Hanus, R., Yang, J., Chen, Y., Huang, M.,Snyder, G.J., et al. (2017). Lattice dislocationsenhancing thermoelectric PbTe in addition toband convergence. Adv. Mater. 29, 1606768.
24. Hsu, K.F., Loo, S., Guo, F., Chen, W., Dyck, J.S.,Uher, C., Hogan, T., Polychroniadis, E.K., andKanatzidis, M.G. (2004). Cubic AgPbmSbTe2+m:bulk thermoelectric materials with high figureof merit. Science 303, 818–821.
25. Biswas, K., He, J., Blum, I.D., Wu, C.I., Hogan,T.P., Seidman, D.N., Dravid, V.P., andKanatzidis, M.G. (2012). High-performancebulk thermoelectrics with all-scale hierarchicalarchitectures. Nature 489, 414–418.
26. Poudel, B., Hao, Q., Ma, Y., Lan, Y., Minnich, A.,Yu, B., Yan, X., Wang, D., Muto, A., Vashaee, D.,et al. (2008). High-thermoelectric performanceof nanostructured bismuth antimony telluridebulk alloys. Science 320, 634–638.
27. Pei, Y., Lensch-Falk, J., Toberer, E.S., Medlin,D.L., and Snyder, G.J. (2011). Highthermoelectric performance in PbTe due tolarge nanoscale Ag2Te precipitates and Ladoping. Adv. Funct. Mater. 21, 241–249.
28. Snyder, G., Christensen, M., Nishibori, E.,Caillat, T., and Iversen, B. (2004). Disorderedzinc in Zn4Sb3 with phonon-glass and electron-crystal thermoelectric properties. Nat. Mater.3, 458–463.
29. Liu, H., Shi, X., Xu, F., Zhang, L., Zhang, W.,Chen, L., Li, Q., Uher, C., Day, T., and Snyder,G.J. (2012). Copper ion liquid-likethermoelectrics. Nat. Mater. 11, 422–425.
30. He, Y., Lu, P., Shi, X., Xu, F., Zhang, T., Snyder,G.J., Uher, C., and Chen, L. (2015). Ultrahighthermoelectric performance in mosaic crystals.Adv. Mater. 27, 3639–3644.
31. Peng, K., Lu, X., Zhan, H., Hui, S., Tang, X.,Wang, G., Dai, J., Uher, C., Wang, G., andZhou, X. (2016). Broad temperature plateau forhigh ZTs in heavily doped p-type SnSe singlecrystals. Energ. Environ. Sci. 9, 454–460.
32. Ying, P., Li, X., Wang, Y., Yang, J., Fu, C., Zhang,W., Zhao, X., and Zhu, T. (2017). Hierarchicalchemical bonds contributing to the intrinsicallylow thermal conductivity in a-MgAgSbthermoelectric materials. Adv. Funct. Mater.27, 1604145.
33. Li, D., Zhao, H., Li, S., Wei, B., Shuai, J., Shi, C.,Xi, X., Sun, P., Meng, S., Gu, L., et al. (2015).Atomic disorders induced by silver andmagnesium ion migrations favor highthermoelectric performance in a-MgAgSb-
based materials. Adv. Funct. Mater. 25, 6478–6488.
34. Zhao, H., Sui, J., Tang, Z., Lan, Y., Jie, Q.,Kraemer, D., McEnaney, K., Guloy, A., Chen,G., and Ren, Z. (2014). High thermoelectricperformance of MgAgSb-based materials.Nano Energy 7, 97–103.
35. Zhao, L.D., Lo, S.H., Zhang, Y., Sun, H., Tan, G.,Uher, C., Wolverton, C., Dravid, V.P., andKanatzidis, M.G. (2014). Ultralow thermalconductivity and high thermoelectric figure ofmerit in SnSe crystals. Nature 508, 373–377.
36. Morelli, D.T., and Slack, G.A. (2006). HighLattice Thermal Conductivity Solids (Springer),pp. 37–68.
37. Slack, G.A. (1979). The thermal conductivity ofnonmetallic crystals. Solid State Phys. 34, 1–71.
38. Deloume, J., and Faure, R. (1981). A newcompound, Ag9GaSe6 structure study of thealpha-phase. J. Solid State Chem. 36, 112–117.
39. Deloume, J., Faure, R., Loiseleur, H., andRoubin, M. (1978). Crystal structure of thephase Ag9GaSe6 beta. Acta Crystallogr. B 34,3189–3193.
40. Tansho, M., Wada, H., Ishii, M., and Onoda, Y.(1996). Silver ionic conductor Ag9GaSe6studied by Ag and Ga NMR. Solid State Ionics86-8, 155–158.
41. Li, W., Lin, S., Ge, B., Yang, J., Zhang, W., andPei, Y. (2016). Low sound velocity contributingto the high thermoelectric performance ofAg8SnSe6. Adv. Sci. 3, 1600196.
42. Li, L., Liu, Y., Dai, J., Hong, A., Zeng, M., Yan, Z.,Xu, J., Zhang, D., Shan, D., Liu, S., et al. (2016).High thermoelectric performance of superionicargyrodite compound Ag8SnSe6. J. Mater.Chem. C 4, 5806–5813.
43. Zhu, T.J., Zhang, S.N., Yang, S.H., and Zhao,X.B. (2010). Improved thermoelectric figure ofmerit of self-doped Ag8�xGeTe6 compoundswith glass-like thermal conductivity. Phys.Status Solidi Rapid Res. Lett. 4, 317–319.
44. Charoenphakdee, A., Kurosaki, K., Muta, H.,Uno, M., and Yamanaka, S. (2008).Reinvestigation of the thermoelectricproperties of Ag8GeTe6. Phys. Status SolidiRapid Res. Lett. 2, 65–67.
45. Fujikane, M., Kurosaki, K., Muta, H., andYamanaka, S. (2005). Thermoelectric propertiesof Ag8GeTe6. J. Alloy. Comp. 396, 280–282.
46. Charoenphakdee, A., Kurosaki, K., Muta, H.,Uno, M., and Yamanaka, S. (2009). Ag8SiTe6: anew thermoelectric material with low thermalconductivity. Jpn. J. Appl. Phys. 48, 011603.
47. Jiang, B., Qiu, P., Eikeland, E., Chen, H., Song,Q., Ren, D., Zhang, T., Yang, J., Iversen, B.B.,Shi, X., et al. (2017). Cu8GeSe6-basedthermoelectric materials with an argyroditestructure. J. Mater. Chem. C 5, 943–952.
48. Weldert, K.S., Zeier, W.G., Day, T.W.,Panthofer, M., Snyder, G.J., and Tremel, W.(2014). Thermoelectric transport in Cu7PSe6with high copper ionic mobility. J. Am. Chem.Soc. 136, 12035–12040.
49. Gibbs, Z.M., LaLonde, A., and Snyder, G.J.(2013). Optical band gap and theBurstein–Moss effect in iodine doped PbTe
Joule 1, 816–830, December 20, 2017 829
using diffuse reflectance infrared Fouriertransform spectroscopy. New J. Phys. 15,075020.
50. Jain, A., and McGaughey, A.J.H. (2014).Thermal conductivity of compoundsemiconductors: interplay of mass density andacoustic-optical phonon frequency gap.J. Appl. Phys. 116, 073503.
51. Lindsay, L., and Kuang, Y. (2017). Effects offunctional group mass variance on vibrationalproperties and thermal transport in graphene.Phys. Rev. B 95, 121404.
52. Kim, C., Kurosaki, K., Ishimaru, M., Muta, H.,and Yamanaka, S. (2011). Effect of vacancydistribution on the thermal conductivity ofGa2Te3 and Ga2Se3. J. Electron. Mater. 40,999–1004.
53. Fernelius, N.C. (1994). Properties of galliumselenide single crystal. Prog. Cryst. GrowthChar. Mater. 28, 275–353.
54. Sanditov, D.S., and Belomestnykh, V.N. (2011).Relation between the parameters of theelasticity theory and averaged bulk modulus ofsolids. Tech. Phys. 56, 1619–1623.
55. Pei, Y., Wang, H., and Snyder, G.J. (2012). Bandengineering of thermoelectric materials. Adv.Mater. 24, 6125–6135.
56. Tang, Y., Gibbs, Z.M., Agapito, L.A., Li, G., Kim,H.S., Nardelli, M.B., Curtarolo, S., and Snyder,G.J. (2015). Convergence of multi-valley bandsas the electronic origin of high thermoelectric
830 Joule 1, 816–830, December 20, 2017
performance in CoSb3 skutterudites. Nat.Mater. 14, 1223–1228.
57. Ravich, Y.I., Efimova, B.A., and Smirnov, I.A.(1970). Semiconducting Lead Chalcogenides(Plenum Press).
58. Pei, Y., LaLonde, A.D., Wang, H., and Snyder,G.J. (2012). Low effective mass leading to highthermoelectric performance. Energ. Environ.Sci. 5, 7963–7969.
59. Shen, J.W., Chen, Z.W., Lin, S.Q., Zheng, L.L.,Li, W., and Pei, Y.Z. (2016). Single parabolicband behavior of thermoelectric p-typeCuGaTe2. J. Mater. Chem. C 4, 209–214.
60. Bardeen, J., and Shockley, W. (1950).Deformation potentials and mobilities in non-polar crystals. Phys. Rev. 80, 72–80.
61. Pei, Y., Gibbs, Z.M., Balke, B., Zeier, W.G., andSnyder, G.J. (2014). Optimum carrierconcentration in n-type PbTe thermoelectrics.Adv. Energy Mater. 4, 1400486.
62. LaLonde, A.D., Pei, Y., and Snyder, G.J. (2011).Reevaluation of PbTe1-xIx as high performancen-type thermoelectric material. Energ. Environ.Sci. 4, 2090–2096.
63. Wang, H., Pei, Y., LaLonde, A.D., and Snyder,G.J. (2012). Weak electron-phonon couplingcontributing to high thermoelectricperformance in n-type PbSe. Proc. Natl. Acad.Sci. USA 109, 9705–9709.
64. Wang, H., Schechtel, E., Pei, Y., and Snyder,G.J. (2013). High thermoelectric efficiency ofn-type PbS. Adv. Energy Mater. 3, 488–495.
65. LaLonde, A.D., Ikeda, T., and Snyder, G.J.(2011). Rapid consolidation of powderedmaterials by induction hot pressing. Rev. Sci.Instrum. 82, 025104.
66. Zhou, Z.H., and Uher, C. (2005). Apparatus forSeebeck coefficient and electrical resistivitymeasurements of bulk thermoelectric materialsat high temperature. Rev. Sci. Instrum. 76,023901.
67. Togo, A., and Tanaka, I. (2015). First principlesphonon calculations in materials science. Scr.Mater. 108, 1–5.
68. Kresse, G., and Furthmuller, J. (1996). Efficiencyof ab-initio total energy calculations for metalsand semiconductors using a plane-wave basisset. Comp. Mater. Sci. 6, 15–50.
69. Blochl, P.E. (1994). Projector augmented-wavemethod. Phys. Rev. B 50, 17953.
70. Perdew, J.P., Burke, K., and Ernzerhof, M.(1996). Generalized gradient approximationmade simple. Phys. Rev. Lett. 77, 3865.
71. Monkhorst, H.J., and Pack, J.D. (1976). Specialpoints for Brillouin-zone integrations. Phys.Rev. B 13, 5188.
72. Koller, D., Tran, F., and Blaha, P. (2011). Meritsand limits of the modified Becke-Johnsonexchange potential. Phys. Rev. B 83, 195134.
JOUL, Volume 1
Supplemental Information
High Thermoelectric Performance
of Ag9GaSe6 Enabled by Low Cutoff
Frequency of Acoustic Phonons
Siqi Lin, Wen Li, Shasha Li, Xinyue Zhang, Zhiwei Chen, Yidong Xu, YueChen, and Yanzhong Pei
Table S1. Elastic properties for Ag9Ga1-xMx(Se1-ySy)6. Room temperature Hall carrier concentration
(nH), transverse (vt), longitudinal (vl) and mean (v) sound velocities, bulk (B) and shear (G) modules,
Gruneisen parameter (γ) and Debye temperature (D) for Ag9GaSe6 with various dopants.
Compounds nH
(cm-3)
vt
(m/s)
vl
(m/s)
v
(m/s)
B
(GPa)
G
(GPa) γ
D
(K)
Ag9GaSe6 7.0e18 1168 2899 1323 48.1 10.0 2.7 142
Ag9Ga0.98Cr0.02Se6 7.6e18 1147 2920 1300 49.4 9.6 2.7 139
Ag9Ga0.96Cr0.04Se6 7.1e18 1206 2857 1364 45.4 10.6 2.5 146
Ag9Ga0.94Cr0.06Se6 6.3e18 1185 2800 1340 43.6 10.3 2.5 143
Ag9Ga0.94Cd0.06Se6 1.1e19 1140 2757 1290 42.8 9.5 2.6 138
Ag9Ga0.94Zn0.06Se6 9.1e18 1145 2851 1297 46.6 9.6 2.7 138
Ag9Ga0.94Ge0.06Se6 9.0e18 1173 2890 1328 47.6 10.0 2.6 142
Ag9Ga(Se0.95S0.05)6 9.0e18 1179 2913 1335 48.4 10.1 2.6 143
Ag9Ga(Se0.90S0.10)6 8.9e18 1230 3097 1394 55.2 11.0 2.7 150
Table S2. Relationship between cut-off frequency and lattice thermal conductivity. Cut-off
frequency for transverse (TA1,TA2) and longitudinal (LA) acoustic phonon branches, the mean
cut-off frequency (m) and the room temperature lattice thermal conductivity(L) for semiconductors.
Compounds TA1 TA2 LA m L Refs.
Argyrodite
semiconductors
Ag9GaSe6 0.51 0.51 0.64 0.54 0.2 This work
Ag8SnSe6 0.44 0.44 0.56 0.47 0.2 1
Cu8GeSe6 0.42 0.42 0.67 0.46 0.3 2
Elemental
semiconductors
Diamond 24.50 24.50 36.40 26.75 2000.0 3, 4
Si 4.41 4.41 12.26 5.01 156.0 4, 5
Ge 2.40 2.40 7.29 2.73 60.0 4, 5
SiC 11.19 11.19 19.20 12.41 490.0 4, 6
Se 1.81 2.16 2.16 2.01 4.5 7, 8
Te 1.04 1.47 1.47 1.25 2.3 8
III-V
semiconductors
BN 9.17 21.55 31.04 12.81 391.0 9, 10
BP 9.59 9.59 16.14 10.61 360.0 11, 12
BAs 5.93 5.93 9.72 6.55 210.0 11, 12
AlN 10.23 10.23 17.78 11.36 319.0 13, 14
AlP 3.98 3.98 10.58 4.52 90.0 11, 15
AlAs 2.80 2.80 6.40 3.16 91.0 5, 16
AlSb 2.19 2.19 4.64 2.47 56.0 5, 14
GaN 5.84 5.84 10.72 6.51 130.0 17, 18
GaP 3.11 3.11 7.65 3.52 100.0 19, 20
GaAs 2.46 2.46 6.83 2.80 44.0 5, 21
GaSb 1.71 1.71 5.00 1.95 33.0 5, 22
II-VI
semiconductors
ZnS 2.66 2.66 6.47 3.01 27.0 23, 24
ZnSe 2.07 2.07 5.59 2.35 19.0 24, 25
ZnTe 1.71 1.71 4.15 1.94 18.0 24, 26
CdSe 1.34 1.34 4.57 1.52 9.0 27, 28
CdTe 1.18 1.18 3.79 1.34 7.5 24, 26
I-III-VI2
semiconductors
CuInSe2 1.53 1.53 1.80 1.60 3.7 29, 30
AgGaSe2 1.02 1.02 1.57 1.12 1.0 31, 32
CdGeAs2 1.99 1.99 2.66 2.14 4.0 30, 33
I2-IV-VI3
&I3-V-VI4
semiconductors
Cu3SbS4 0.76 0.76 0.94 1.74 4.1 34
Cu3SbSe3 1.64 1.77 1.83 0.80 0.5 35
Cu3SbSe4 1.31 1.38 1.42 1.37 2.9 34
IIx-IVy
semiconductors
Mg2Ge 2.27 2.27 3.63 2.50 6.6 36, 37
Mg2Sn 2.25 2.25 3.73 2.49 13.0 30, 38
IVx-VIy
semiconductors
GeSe 0.45 0.45 0.76 1.35 2.0 28, 39
PbS 1.23 1.23 1.63 1.32 2.5 28, 39
PbSe 0.97 0.97 1.49 1.06 1.6 40, 41
PbTe 0.73 0.73 0.96 0.78 2.0 41, 42
SnTe 1.49 1.49 1.03 1.26 2.8 43, 44
GeTe 1.52 1.67 2.04 1.69 3.0 45, 46
Typical
thermoelctrics
MgAgSb 0.73 0.92 1.10 0.87 0.6 47
CaAgSb 0.99 0.99 1.44 1.08 1.7 47, 48
CoSb3 1.57 1.57 2.62 1.74 10.0 49, 50
Ba8Ga16Ge30 1.06 1.06 1.44 1.14 1.6 51, 52
Bi2Te3 0.50 0.50 0.79 0.55 1.3 53
The mean cut-off frequency (m) for acoustic phonons is obtained by:
𝜔𝑚 = 1
3 1
𝜔𝑇𝐴13 +
1
𝜔𝑇𝐴23 +
1
𝜔𝐿𝐴3
−13
Ag10 mm Ga
Se Te
10 mm
10 mm10 mm
10 mm
(a) (b) (c)
(d) (e)
Ag9Ga(Se0.9Te0.1)6
Ag: Ga: Se: Te53.79: 7.60: 35.41: 3.19
(f)
Figure S1. Microstructure of Ag9Ga(Se0.9Te0.1)6. SEM image (a), EDS composition mapping of Ag
(b), Ga (c), Se (d), Te(e) and the energy spectrum (f) for Ag9Ga(Se0.9Te0.1)6.
0.5 0.6 0.7 0.8
x=0, y=0
x=0.02, y=0 Cr
x=0.04, y=0 Cr
x=0.06, y=0 Cr
x=0.06, y=0 Cd
x=0.06, y=0 Zn
x=0.06, y=0 Ge
x=0, y=0.05
x=0, y=0.10
Ag9Ga
1-xM
x(Se
1-yS
y)6 (M=Cr,Cd,Zn,Ge)
A
bso
rpti
on
(N
orm
alize
d)
hv (eV)
Eg=0.56ev
20 30 40 50 60 70 80 90
x=0, y=0.10
Ag9Ga
1-xM
x(Se
1-yS
y)6 (M=Cr,Cd,Zn,Ge)
x=0, y=0.05
x=0.06, y=0Ge
x=0.06, y=0Zn
x=0.06, y=0Cd
x=0.06, y=0Cr
x=0.04, y=0Cr
x=0.02, y=0Cr
Inte
ns
ity
(a
.u.)
2 (deg.)
x=0, y=0
ICSD #15252
(a) (b)
Figure S2. XRD and optical measurements for Ag9Ga1-xMx(Se1-ySy)6. Powder X-ray diffraction
patterns (a) and the normalized optical absorption versus photon energy at room temperature (b) for
Ag9Ga1-xMx(Se1-ySy)6. All the peaks can be well indexed to the corresponding F43m structure,
indicating the high phase purity. The optical measurements shows that the band gap remains nearly
unchanged for various doping/alloying for Ag9GaSe6.
1017
1018
1019
0
1
2
3
4
5
6
7
8; ; ; Prediction
; ; ; Exp. Ag9Ga(Se
1-xTe
x)
; ; ; Exp. other doping&alloys
n
H (cm
-3)
PF
(m
W/c
m-K
2)
300K
500K
800K
m* ~ 0.16me
Edef ~ 24eV
300 400 500 600 700 800 900
-150
-100
-50
Ag9Ga
1-xM
x(Se
1-yS
y)
6 (M=Cr,Cd,Zn,Ge)
S (m
V/K
)
T (K)300 400 500 600 700 800 900
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
4.5Ag
9Ga
1-xM
x(Se
1-yS
y)
6 (M=Cr,Cd,Zn,Ge)
x=0, y=0
x=0.02, y=0 Cr x=0.04, y=0 Cr
x=0.06, y=0 Cr x=0.06, y=0 Cd
x=0.06, y=0 Zn x=0.06, y=0 Ge
x=0, y=0.05 x=0, y=0.10
(
m
cm
)
T (K)
300 400 500 600 700 800 9000.0
0.2
0.4
0.6
0.8
1.0
1.2
Ag9Ga
1-xM
x(Se
1-yS
y)6 (M=Cr,Cd,Zn,Ge)
x=0, y=0
x=0.02, y=0 Cr
x=0.04, y=0 Cr
x=0.06, y=0 Cr
x=0.06, y=0 Cd
x=0.06, y=0 Zn
x=0.06, y=0 Ge
x=0, y=0.05
x=0, y=0.10
zT
T (K)300 400 500 600 700 800 900
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
L~0.15W/m-K
Standard deviation ~ 0.058 W/m-K
,
L (
W/m
-K)
T (K)
(a) (b)
(c) (d)
(e) (f)
300 400 500 600 700 8000
10
20
30
m*
(me)
Ed
ef (
ev)
T (K)
0.0
0.2
0.4
0.6
0.8
Figure S3. Transport properties for Ag9Ga1-xMx(Se1-ySy)6. Temperature dependent Seebeck
coefficient (a), resistivity (b),density of state effective mass, m* and deformation potential coefficient,
Edef (c) for Ag9Ga1-xMx(Se1-ySy)6. The comparison between the measurements and prediction on power
factor versus Hall carrier concentration (d). Temperature dependent total thermal conductivity, lattice
thermal conductivity (e) and figure of merit (zT) for Ag9Ga1-xMx(Se1-ySy)6.
300 400 500 600 700 800 9001.5
1.6
1.7
1.8
1.9
2.0
2.1
2.2
2.3 x=0 1.2E19
x=0.05 8.5E18
x=0.08 6.6E18
x=0.10 3.1E18
x=0.12 4.1E18
x=0.15 4.0E18
x=0.17 2.3E18
Ag9Ga(Se
1-xTe
x)6
T (K)
L ×
10
8 (
V2/K
2)
Figure S4. Lorenz factor for Ag9Ga(Se1-xTex)6. Temperature dependent Lorenz factor for
Ag9Ga(Se1-xTex)6.
Supplemental References
1. Li, W., Lin, S., Ge, B., Yang, J., Zhang, W., Pei, Y. (2016). Low Sound Velocity Contributing to the High Thermoelectric
Performance of Ag8SnSe6. Adv. Sci. 3, 1600196.
2. Jiang, B., Qiu, P., Eikeland, E., Chen, H., Song, Q., Ren, D., Zhang, T., Yang, J., Iversen, B.B., Shi, X.et al. (2017).
Cu8GeSe6-based thermoelectric materials with an argyrodite structure. J. Mater. Chem. C. 5, 943-952.
3. Warren, J.L., Yarnell, J.L., Dolling, G., Cowley, R.A. (1967). Lattice Dynamics of Diamond. Phy. Rev. 158, 805-808.
4. Slack, G.A. (1964). Thermal Conductivity of Pure and Impure Silicon, Silicon Carbide, and Diamond. J. Appl. Phys. 35, 3460.
5. Giannozzi, P., de Gironcoli, S., Pavone, P., Baroni, S. (1991). Ab initiocalculation of phonon dispersions in semiconductors.
Phys. Rev. B 43, 7231-7242.
6. Hofmann, M., Zywietz, A., Karch, K., Bechstedt, F. (1994). Lattice dynamics of SiC polytypes within the bond-charge model.
Phys. Rev. B 50, 13401-13411.
7. Hansen, F.Y., McMurry, H.L. (1980). Normal mode calculations of trigonal selenium. J. Chem. Phys. 72, 5550-5564.
8. E. Gerlach, P. Grosse (1979). The Physics of Selenium and Tellurium (Springer Berlin).
9. Kern, G., Kresse, G., Hafner, J. (1999). Ab initio calculation of the lattice dynamics and phase diagram of boron nitride. Phys.
Rev. B 59, 8551-8559.
10. Sichel, E.K., Miller, R.E., Abrahams, M.S., Buiocchi, C.J. (1976). Heat capacity and thermal conductivity of hexagonal pyrolytic
boron nitride. Phys. Rev. B 13, 4607-4611.
11. Pletl, T., Pavone, P., Engel, U., Strauch, D. (1999). First-principles study of lattice-dynamical and elastic trends in tetrahedral
semiconductors. Phys. B 263, 392-395.
12. Slack, G.A. (1973). Nonmetallic crystals with high thermal conductivity. J. Phys. and Chem. Solids 34, 321-335.
13. Nipko, J.C., Loong, C.K. (1998). Phonon excitations and related thermal properties of aluminum nitride. Phys. Rev. B 57,
10550-10554.
14. Slack, G.A., Tanzilli, R.A., Pohl, R.O., Vandersande, J.W. (1987). The intrinsic thermal conductivity of AIN. J. Phys. Chem.
Solids 48, 641-647.
15. Steigmeier, E.F. (1963). The Debye Temperatures of Iii‐V Compounds. Appl. Phys. Lett. 3, 6-8.
16. Afromowitz, M.A. (1973). Thermal conductivity of Ga1−xAlxAs alloys. J. Appl. Phys. 44, 1292-1294.
17. Karch, K., Wagner, J.M., Bechstedt, F. (1998). Ab initiostudy of structural, dielectric, and dynamical properties of GaN. Phys.l
Rev. B 57, 7043-7049.
18. Sichel, E., Pankove, J. (1977). Thermal conductivity of GaN, 25–360 K. J. Phys. Chem. Solids 38, 330.
19. Yarnell, J., Warren, J., Wenzel, R., Dean, P. (1968) Lattice Dynamics of Gallium Phosphide. In: Neutron Inelastic Scattering Vol.
I. Proceedings of a Symposium on Neutron Inelastic Scattering (International Atomic Energy Agency ViennaPress).
20. Muzhdaba, V., NASHELSKY, A., Tamarin, P., Shalyt, S. (1969). Thermal conductivity and thermo-emf of AlSb and GaP at low
temperatures. Sov. Phys. Solid State 10, 9.
21. Blakemore, J.S. (1982). Semiconducting and other major properties of gallium arsenide. J. Appl. Phys. 53, R123-R181.
22. Maycock, P.D. (1967). Thermal conductivity of silicon, germanium, III–V compounds and III–V alloys. Solid-State Electron. 10,
161-168.
23. Talwar, D.N., Vandevyver, M., Kunc, K., Zigone, M. (1981). Lattice dynamics of zinc chalcogenides under compression:
Phonon dispersion, mode Grüneisen, and thermal expansion. Phys. Rev. B 24, 741-753.
24. Slack, G.A. (1972). Thermal Conductivity of II-VI Compounds and Phonon Scattering by Fe2+ Impurities. Phys. Rev. B 6,
3791-3800.
25. Kushwaha, M.S., Kushwaha, S.S. (1980). Lattice dynamics of II–VI and III–V compounds. J. Phys. Chem. Solids 41, 489-494.
26. Kushwaha, M.S., Kushwaha, S.S. (1980). Lattice dynamics of ZnTe, CdTe, GaP, and InP. CAN. J. Phys. 58, 351-358.
27. Dal Corso, A., Baroni, S., Resta, R., de Gironcoli, S. (1993). Ab initiocalculation of phonon dispersions in II-VI semiconductors.
Phys. Rev. B 47, 3588-3592.
28. Madelung, O. (2012). Semiconductors: data handbook (Springer Science & Business Media).
29. Fouret, R., Hennion, B., Gonzalez, J., Wasim, S.M. (1993). Elastic stiffness constants of copper indium diselenide determined
by neutron scattering. Phys. Rev. B 47, 8269-8272.
30. Spitzer, D.P. (1970). Lattice thermal conductivity of semiconductors: A chemical bond approach. J. Phys. Chem. Solids 31,
19-40.
31. Artus, L., Pujol, J., Pascual, J., Camassel, J. (1990). Lattice dynamics ofAgGaSe2. II. Theoretical aspects. Phys. Rev. B 41,
5727-5734.
32. Aggarwal, R.L., Fan, T.Y. (2005). Thermal diffusivity, specific heat, thermal conductivity, coefficient of thermal expansion, and
refractive-index change with temperature in AgGaSe2. Appl. Optics 44, 2673.
33. Tyuterev, V.G., Skachkov, S.I. (1992). Lattice dynamics, thermodynamic and elastic properties of CdGeAs2. II Nuovo Cimento
D 14, 1097-1103.
34. Suzumura, A., Watanabe, M., Nagasako, N., Asahi, R. (2014). Improvement in Thermoelectric Properties of Se-Free Cu3SbS4
Compound. J. Electron. Mater. 43, 2356-2361.
35. Qiu, W., Xi, L., Wei, P., Ke, X., Yang, J., Zhang, W. (2014). Part-crystalline part-liquid state and rattling-like thermal damping in
materials with chemical-bond hierarchy. Proc. Natl. Acad. Sci.USA 111, 15031-15035.
36. Chung, P.L., Whitten, W.B., Danielson, G.C. (1965). Lattice dynamics of Mg2Ge. J. Phys. Chem. Solids 26, 1753-1760.
37. LaBotz, R.J., Mason, D.R., O'Kane, D.F. (1963). The Thermoelectric Properties of Mixed Crystals of Mg2GexSi1−x. J.
Electrochem. Soc. 110, 127.
38. Kearney, R.J., Worlton, T.G., Schmunk, R.E. (1970). Lattice dynamics of magnesium stannide at room temperature. Journal of
Phys. Chem. Solids 31, 1085-1097.
39. Okhotin, A.S., Krestovnikov, A.N., Aivazov, A.A., Pushkarskii, A.S. (1969). Thermal Conductivity of GeS and GeSe. Phys. Status
Solid (b) 31, 485-487.
40. Tian, Z., Garg, J., Esfarjani, K., Shiga, T., Shiomi, J., Chen, G. (2012). Phonon conduction in PbSe, PbTe, and PbTe1−xSex from
first-principles calculations. Phys. Rev. B 85,184303.
41. Wang, H., Schechtel, E., Pei, Y., Snyder, G.J. (2013). High Thermoelectric Efficiency of n-type PbS. Adv. Energy Mater. 3,
488-495.
42. Cochran, W., Cowley, R.A., Dolling, G., Elcombe, M.M. (1966). The Crystal Dynamics of Lead Telluride. Proceedings of the
Royal Society A: Mathematical, Phys. Eng. Sci. 293, 433-451.
43. Cowley, E.R., Darby, J.K., Pawley, G.S. (1969). The lattice dynamics of tin telluride. J. Phys. C: Solid State Phys. 2, 1916-1925.
44. Pei, Y., Zheng, L., Li, W., Lin, S., Chen, Z., Wang, Y., Xu, X., Yu, H., Chen, Y., Ge, B. (2016). Interstitial Point Defect Scattering
Contributing to High Thermoelectric Performance in SnTe. Adv. Electron. Mater. 2, 1600019.
45. Shaltaf, R., Gonze, X., Cardona, M., Kremer, R.K., Siegle, G. (2009). Lattice dynamics and specific heat of α-GeTe: Theoretical
and experimental study. Phys. Rev. B 79,075204.
46. Li, J., Chen, Z., Zhang, X., Sun, Y., Yang, J., Pei, Y. (2017). Electronic origin of the high thermoelectric performance of GeTe
among the p-type group IV monotellurides. NPG Asia Mater. 9, e353.
47. Ying, P., Li, X., Wang, Y., Yang, J., Fu, C., Zhang, W., Zhao, X., Zhu, T. (2017). Hierarchical Chemical Bonds Contributing to the
Intrinsically Low Thermal Conductivity in α-MgAgSb Thermoelectric Materials. Adv. Funct. Mater. 27, 1604145.
48. Wang, J., Liu, X.C., Xia, S.Q., Tao, X.T. (2013). Ca1-XRExAg1-ySb (RE = La, Ce, Pr, Nd, Sm; 0 ≤x≤1; 0 ≤y≤1): Interesting
Structural Transformation and Enhanced High-Temperature Thermoelectric Performance. J. Am. Chem. Soc. 135,
11840-11848.
49. Rotter, M., Rogl, P., Grytsiv, A., Wolf, W., Krisch, M., Mirone, A. (2008). Lattice dynamics of skutterudites: Inelastic x-ray
scattering onCoSb3. Phys. Rev. B 77, 144301.
50. Caillat, T., Borshchevsky, A., Fleurial, J.P. (1996). Properties of single crystalline semiconducting CoSb3. J. Appl. Phys. 80,
4442-4449.
51. Christensen, M., Abrahamsen, A.B., Christensen, N.B., Juranyi, F., Andersen, N.H., Lefmann, K., Andreasson, J., Bahl, C.R.,
Iversen, B.B. (2008). Avoided crossing of rattler modes in thermoelectric materials. Nat. Mater. 7, 811-815.
52. May, A.F., Toberer, E.S., Saramat, A., Snyder, G.J. (2009). Characterization and analysis of thermoelectric transport in n-type
Ba8Ga16−xGe30+x. Phys. Rev. B 80,105825.
53. Qiu, B., Ruan, X. (2009). Molecular dynamics simulations of lattice thermal conductivity of bismuth telluride using two-body
interatomic potentials. Phys. Rev. B 80,165203.