

Thermoelectric Properties of Ag-doped In4Se2.95 Polycrystalline Compounds
English
Thermoelectric Properties of Ag-doped In4Se2.95 Polycrystalline Compounds
-
Key words:
- In4Se2.95
- / Ag doping
- / thermoelectric property
-
1. INTRODUCTION
Thermoelectric (TE) materials can convert thermal energy from a temperature gradient into electrical energy and vice versa. The TE devices composed of n- and p-type TE materials have been the alternative source of global sustainable energies because of its many excellent advantages, such as no moving parts, long operating lifetimes, reliability and simplicity[1-3]. The high TE efficiency of TE devices should require both a large working temperature gradient and a high TE figure-of-merit ZT = S2σ T /κ , where S, σ, κ, and T are Seebeck coefficient, electrical conductivity, total thermal conductivity and absolute temperature, respectively. Therefore, the TE materials having high ZT values paly a key role in the TE application. Among the huge family of TE materials, several TE materials containing heavy metals such as Pb, Bi, Sb and Te have been developed so far and they are mostly toxic or rare in the earth. Therefore, the researches on the environment friendly and earth-abundant TE materials are of great important such as the lead-free In4Se3-based TE materials.
As shown in Fig. 1a, the quasi two-dimensional layered crystal structure of In4Se3 consists of In+ ions and [In3Se3]- corrugated layers, that is, In4Se3 = [In]+[In3]5+[Se2-]3. The anionic layers of [In3Se3]- are interconnected by In+ ions (In4 atoms) via ionic bonding interactions and in each layer, the (In3)5+ multivalent clusters which comprise In1, In2 and In3 atoms are connected by Se ions (Se 1~3) with covalent bond[4, 5]. As a new kind of n-type chalcogenide TE material, In4Se3-based TE materials have become the promising medium temperature (500~900 K) n-type TE material candidates since the remarkable ZT = 1.48 at 705 K (In4Se2.35 single crystal) along the bc plane[6] was reported. The intrinsically low lattice thermal conductivity (κlat) and high Seebeck coefficient of In4Se3-δ crystals (with Se-deficiency, δ) originate from the quantum confinement effect and the charge density wave induced by Periels distortion along the quasi-1D indium chain[6]. Nevertheless, the In4Se3-δ single crystals suffer strong anisotropy[6, 7] and the cleavage habit, which are not perfect for application, so it is very necessary to research on the In4Se3-δ polycrystalline compounds.
It is a feasible approach to improve TE properties of In4Se3-based TE materials by Se-deficiency[8-10] and dual even multi-doping[11-13]. It aims at increasing the power factor (PF = S2σ) through optimizing conductivity performance because of the unsatisfied carrier concentration (ne ~ 1017 cm-3) in In4Se3 polycrystalline compound[14-16]. In the past years, the doping effect of many cationic dopants on the In4Se3-based TE materials have been reported and the results indicate that doping in different lattice sites have distinct influence on the conductivity performance. For example, cationic mono-doping at In2 (e. g. Sn)[17, 18] or In3 (e. g. Zn and Ga)[18] is reported to be non-effective for optimizing conductivity performance, while substitution at In4 (e. g. Na, Ca and Pb)[12, 18] or intercalation sites (e. g. Cu)[12] is verified to be more effective to optimizing conductivity performance because of the obviously increased ne. Herein, given that the atomic radius of Ag is smaller than that of In, and most importantly, Ag should be easier to ionize than Cu as a donor when doping at intercalation sites due to its lower first ionization energy, we investigated experimentally and theoretically the effect of Ag mono-doping on the TE properties of In4Se2.95 polycrystalline compounds. The maximum ZT of 0.92 at 723 K is obtained byIn3.98Se983 Ag0.02 compound which is 15% enhancement over that of the un-doped one.
2. EXPERIMENTAL
2.1. Synthesis and characterization
In4-xSe2.95Agx (x = 0.0, 0.01, 0.02, 0.03 and 0.05) polycrystalline compounds were synthesized by the solid state methods and hot-press technique. Stoichiometric amounts of In (granular, 99.995%), Se (granular, 99.999%) and Ag (granular, 99.999%) were weighed and mixed in quartz tubes in an Ar-filled glove box. The tubes were flame sealed under 0.1 Pa and then slowly heated to 1173 K in tube furnace and kept for 1 day followed by water quenching. The melted ingots were annealed at 743 K for 5 days and then ground into fine powders and subsequently sintered by hot-press technique under uniaxial pressure of 57 MPa at 793 K for 50 minutes under vacuum. The experimental relative densities of all hot-press samples were above 98%.
The powder X-ray diffraction (PXRD) data were obtained on a MiniFlex II powder diffractometer with a CuKα radiation in a 2θ range from 10° to 60° with a step size of 0.02°. The electrical resistivity (ρ) and Seebeck coefficient (S) were measured simultaneously from room temperature to 723 K under an argon atmosphere by using a standard four-probe method with the aid of a ULVAC-RIKO ZEM-3 instrument system. The total thermal conductivity (κtot) was calculated byκ tot =αCpd , where α, Cp and d were thermal diffusivity, specific heat and density, respectively. The thermal diffusivity was measured on a sample disk with about 10.0 mm in diameter and about 2.0 mm in thickness by a laser flash method (LFA-457, NETZSCH, Germany) under an argon atmosphere in the temperature range of 300~723 K. The specific heat was directly derived from Dulong-Petit value (0.2508 J·K-1·g-1)[6, 19]. The measured density was determined using the dimensions and mass of the sample and then reconfirmed by measurements using Archimedes' principle on a home-built device.
2.2. Details of the first-principle calculations
The calculations were based on density functional theory (DFT)[20], with the aid of Vienna ab-initio simulation package (VASP). The generalized gradient approximation (GGA)[21] was implemented by Perdew, Burke, and Ernzerhof (PBE) for the exchange correlation energy functional and a plane-wave basis with the projector augmented wave (PAW)[22] potentials was used. We took into account In4Se3 structure with unit cell composed of 28 atoms including 16 indium and 12 selenium atoms to investigate the effect of metal substitution in an indium site. There are seven different sites in the unit cell of In4Se3 (e.g. four sites for indium and three sites for selenium) and each site occupies four equivalent atomic positions. Each site of indium atoms at four different indium sites was chosen for the metal substitution, coupled with the interstitial position. In the geometrical optimization, an energy cutoff of 400 eV and a k-point spacing with 0.03 Å-1 in the Brillouin zone was used. Atomic positions were fully relaxed until the Hellman- Feynman force on each atom was less than 0.02 eV·Å-1. In the static self-consistent-field calculation, a plane-wave cutoff energy of 500 eV and threshold of 10-5 eV were set. A k-point spacing with 0.02 Å-1 in the Brillouin zone was used by the tetrahedron method[23].
3. RESULTS AND DISCUSSION
Fig. 1b shows PXRD patterns of the In4-xSe2.95Agx (x = 0~0.5) samples which are well indexed as Pnnm In4Se3 (PDF-83-0039) with no noticeable impurity peaks except for the x = 0.03 and 0.05 samples. The arrival of detectable impurity phase InSe (PDF-34-1431) indicates that the maximum solubility of Ag (xm) is reached (xm < 0.03). It is worth to mention that the reason why the element Ag can not be found in PXRD patterns is that the amount of Ag dopant is below the detective limit of XRD technique (usually 5 wt%). Besides, it is known that the atom radius of Ag (1.75 Å) is smaller than that of In (2.0 Å). Thus, to evaluate the stability of Ag atom occupying the In1, In2, In3, In4 and intercalating (Agint) sites, the formation energy calculation was carried out and shown as follows:
In16Se12-In + Ag → In15AgSe12
(substituting the In site)
In16Se12 + Ag → In16AgSe12
(entering the intercalating site)
Ef = Et (In15Ag1Se12) - Et (In16Se12) +
Et (In) - Et (Ag)
Ef = Et (In16Ag1Se12) - Et (In16Se12) - Et (Ag)
Figure 1
Figure 1. a) Crystal structure of In4Se3 viewed along the c axis with selected atom numbers and unit cell marked. b) PXRD of the In4-xSe2.95Agx (x = 0~0.05) polycrystalline compounds. Column chart in the bottom and asterisks (*) refer to In4Se3 and InSe phases, respectivelywhere Ef is the formation energy, and Et the total energy of each material. The formation energies for Ag occupying each site are shown in Table 1. Obviously, the formation energies in the intercalating and In4 sites are lower than the values in that of other In 1~3 sites. The result shows that Ag occupying the intercalating or In4 sites is more energetically stable than in other In sites.
Table 1
In1 In2 In3 In4 Agint Ag 1.504 1.351 1.58 0.587 0.556 The temperature dependence electrical resistivity ρ(T) of In4-xSe2.95Agx samples is shown in Fig. 2a. The ρ of all samples decreases and tends to convergence with the elevated temperature, while it is unsatisfied because the electrical conductivity behavior is not drastically optimized by cationic Ag-doping. Based on the density of states (DOS) calculation shown in Fig. 2d, after Ag doping, the Fermi level moves into the conduction band for Ag inserting the intercalating site while it is nearly unchanged for Ag substituting in the In4 site. Therefore, it predicts an n-type transport behavior for Ag in intercalation site as the Cu-doping[12] in the intercalation site. However, the maximum solubility of Ag (xm < 0.03) is so low that the Ag-doping effect on the electrical conductivity behavior is not significant. Besides, the negative values of S (Fig. 2b) confirm that all samples are n-type TE materials. It is noted that there exists a peak in the S (T) curve as similar with many previous studies[10-13], which is related to the electron-hole bipolar transport. In the bipolar region, the contribution to the final S value should include positive and negative charge carriers and the S can be expressed by the conductivity weighted average[24, 25]:
$S=\frac{{{\sigma }_{p}}{{S}_{p}}-{{\sigma }_{n}}\left. {{S}_{n}} \right.}{{{\sigma }_{p}}+{{\sigma }_{n}}}$ $ Figure 2
At low temperature, the population of minority carriers (holes for n-type TE materials) is so small that the contribution to the overall S value mainly stems from the electron conductivity (σn); while at higher temperature, the minority carrier conductivity (σp) will increase rapidly, thus resulting in the reduction in S. Therefore, the peak appears. The electron-hole bipolar effect mainly happens at the narrow band gap semiconductor, and the band gap can be estimated via the relation developed by Goldsmid and Sharp[25]:
${{E}_{g}}=2e{{\left. S \right.}_{\max }}{{T}_{\max }}$ $ The estimated band gaps of all samples are located at the range of 0.425~0.435 eV which indicates that the amount of Ag doping is too small to change the band gap of In4Se2.95 compound. According to the above analysis, it can be clearly seen that owing to the slight effect of Ag doping on optimizing electrical conductivity behavior, the power factors (PF = S2σ) of all samples increase with increasing temperature but not enhanced significantly as displayed in Fig. 2c. In detail, the maximum value of PF = 5.38 μ·W·K-2·cm-1 at 723 K is obtained by In3.98Se983 Ag0.02 compound.
Fig. 3a shows the temperature-dependent total thermal conductivity κtot (T) of all samples. This conventional 1/T behavior is mainly caused by an acoustic phonon contribution to thermal transport. In the temperature range of 300 to 723 K, when x < 0.03, the κtot values of all samples are lower that of un-doped one, and then increase with the increased Ag doping concentration. Interestingly, the κtot (T) curves of In4-xSe2.95Agx (x = 0.01, 0.02) samples are nearly coincident and the same phenomenon happens in the In4-xSe2.95Agx (x = 0.03, 0.05) samples. At 723 K, the κtot of In4-xSe2.95Agx (x = 0.01, 0.02) and In4-xSe2.95Agx (x = 0.03, 0.05) samples are 0.42 and 0.47 W·m-1·K-1, respectively. In general, the κtot is the sum of electrical thermal conductivity (κe) and lattice thermal conductivity (κlat). The κe can be calculated by the Wiedemann- Franz relationκ e = LT / ρ , where L is the Lorenz number. In metallic system, the L0 = π2/3(kB/e)2 = 2.45×10-8 W·Ω·K-2. However, the Lorenz number is incorrect for a degenerate semiconductor. Herein, to obtain a more reliable Lorenz number, the L is calculated based on the single parabolic band (SPB) model dominated by acoustic phonon scattering (r = -1/2) (Fig. 3b) as the following expressions[26-28]:
$S=\frac{{{k}_{B}}}{e}=\left[ \frac{2{{F}_{1}}\left( \eta \right)}{{{F}_{0}}\left( \eta \right)}-\eta \right]$ $ ${{F}_{n}}\left( \eta \right)=\int_{0}^{\infty }{\frac{{{x}^{n}}}{1+{{e}^{x-\eta }}}}dx$ $ $L={{\left( \frac{{{k}_{B}}}{e} \right)}^{2}}\frac{3{{F}_{0}}\left( \eta \right){{F}_{2}}\left( \eta \right)-4F_{1}^{2}\left( \eta \right)}{F_{0}^{2}\left( \eta \right)}$ $ Figure 3
where η, Fn (η), kB and h are the reduced Fermi energy (η = EF / kB T ), nth order Fermi integral, Boltzmann constant and Planck constant, respectively. Therefore, using the calculated L and ρ, we can obtain κe (T) and then κlat (T) by subtracting κe from κtot, as shown in Fig. 3c. All doped samples exhibit similar κe but differ in κlat which results in various κtot. Therefore, the lower κtot mainly comes from the reduced κlat, especially for x = 0.01 and 0.02 samples. Basically, the low κlat of In4Se3-based TE materials is closely related to the positive effect of Peierls distortion[6].
As manifested in Fig. 3d, the ZT values increase with increasing the temperature for all samples. Obviously, a maximum ZT value reaching 0.92 at 723 K inIn3.98Se983 Ag0.02 sample is about 15% enhancement over the pristine In4Se2.95 sample (ZT = 0.8, 723 K). The enhanced ZT for x = 0.01 and 0.02 samples mainly originates from the reduced κtot because of the reduced κlat. A download trend of ZT value when x > 0.03 may be also connected with the existence of InSe impurity phase because the negative influence of InSe phase was also found in the In4Pb0.01Sn0.03Se2.9Cl0.06 compound[27]. The results indicate that Ag mono-doping is not obviously effective to improve the TE property of In4Se2.95 polycrystalline compounds.
4. CONCLUSION
In summary, In4-xSe2.95Agx (x = 0~0.05) polycrystalline compounds are synthesized successfully and the effect of Ag mono-doping on the electrical and thermal properties are investigated theoretically and experimentally. The maximum solubility of Ag-dopant (xm) in In4Se2.95 is very low (xm < 0.03). According to the formation energy calculation, Ag-occupation at intercalating or In4 sites is more energetically stable than at other In sites. However, the electrical transport behavior is nearly unchanged by the cationic Ag mono-doping. Consequently, the maximum ZT of 0.92 (723 K) obtained by In3.98Se983 Ag0.02 polycrystalline compound is 15% enhancement over the pristine one. The enhancement on ZT mainly benefits from the slightly enhanced PF and the reduced κtot. Considering the dual even multi-doping effect and self-doping effect of Se-deficient, further researches on the In4Se3- based TE materials can focus on the synergistic effect of Se deficiency and Ag-doping or the co-doping of Ag and another element.
-
-
[1]
Elsheikh M. H, Shnawah D. A, Sabri M. F. M, Said S. B. M, Hassan M. H, Bashir M. B. A, Mohamad M. A review on thermoelectric renewable energy: principle parameters that affect their performance[J]. Renew. Sust. Energ. Rev, 2014, 30: 337-355. doi: 10.1016/j.rser.2013.10.027
-
[2]
Sales B. C. Smaller is cooler[J]. Science, 2002, 295: 1248-1249. doi: 10.1126/science.1069895
-
[3]
DiSalvo F. J. Thermoelectric cooling and power generation[J]. Science, 1999, 285: 703-706. doi: 10.1126/science.285.5428.703
-
[4]
Losovyj Y. B, Makinistian L, Albanesi E. A, Petukhov A. G, Liu J, Galiy P, Dveriy O. R, Dowben P. A. The anisotropic band structure of layered In4Se3 (001)[J]. J. Appl. Phys, 2008, 104: 083713. doi: 10.1063/1.3000453
-
[5]
Hogg J. H. C, Sutherland H. H, Williams D. J. Crystallographic evidence for the existence of the phases In4Se3 and In4Te3 which contain the homonuclear triatornic cation (In3)5+[J]. Chem. Commun, 1971, : 1568-1569.
-
[6]
Rhyee J. S, Lee K. H, Lee S. M, Cho E, Kim S. II, Lee E, Kwon Y. S, Shim J. H, Kotliar G. Peierls distortion as a route to high thermoelectric performance in In4Se3-δ crystals[J]. Nature, 2009, 459: 965-968. doi: 10.1038/nature08088
-
[7]
Zhai Y. B, Zhang Q. S, Jiang J, Zhang T, Xiao Y. K, Yang S. G, Xu G. J. Thermoelectric performance of the ordered In4Se3-In composite constructed by monotectic solidification[J]. J. Mater. Chem. A, 2013, 1: 8844-8847. doi: 10.1039/c3ta01599h
-
[8]
Zhu G. H, Lan Y. C, Wang H, Joshi G, Chen Q, Ren Z. F. Effect of selenium deficiency on the thermoelectric properties of n-type In4Se3-x compounds[J]. Phys. Rev. B, 2011, 83: 115201. doi: 10.1103/PhysRevB.83.115201
-
[9]
Rhyee J. S, Cho E, Lee K. H, Lee S. M, Kim S. II, Kim H. S, Kwon Y. S, Kim S. J. Thermoelectric properties and anisotropic electronic band structure on the In4Se3-x compounds[J]. Appl. Phys. Lett, 2009, 95: 212106. doi: 10.1063/1.3266579
-
[10]
Kim J. H, Kim M. J, Oh S, Rhyee J. S. Thermoelectric properties of Se-deficient and Pb-/Sn-codoped In4Pb0.01Sn0.03Se3-x polycrystalline compounds[J]. J. Alloys Compd, 2014, 615: 933-936. doi: 10.1016/j.jallcom.2014.06.196
-
[11]
Luo Y. B, Yang J. Y, Liu M, Xiao Y, Fu L. W, Li W. X, Zhang D, Zhang M. Y, Cheng Y. D. Multiple heteroatoms induced carrier engineering and hierarchical nanostructure for high thermoelectric performance of polycrystalline In4Se2.5[J]. J. Mater. Chem. A, 2015, 3: 1251-1257. doi: 10.1039/C4TA05508J
-
[12]
Luo Y. B, Yang J. Y, Li G, Liu M, Xiao Y, Fu L. W, Li W. X, Zhu P. W, Peng J. Y, Gao S. Z, Zhang J. Q. Enhancement of the thermoelectric performance of polycrystalline In4Se2.5 by copper intercalation and bromine substitution[J]. Adv. Energy Mater, 2014, 4: 1300599. doi: 10.1002/aenm.201300599
-
[13]
Lin Z. S, Chen L, Wang L. M, Zhao J. T, Wu L. M. A promising mid-temperature thermoelectric material candidate: Pb/Sn-codoped In4PbxSnySe3[J]. Adv. Mater, 2013, 25: 4800-4806. doi: 10.1002/adma.v25.34
-
[14]
He S. H, Lin Z. X, Muhammad A. K, Liu P. F, Chen L, Wu L. M. Thermoelectric properties of Ni-substituted polycrystalline In4Se3[J]. Chin. J. Struct. Chem, 2015, 34: 1217-1223.
-
[15]
Shi X, Cho J. Y, Salvador J. R, Yang J, Wang H. Thermoelectric properties of polycrystalline In4Se3 and In4Te3[J]. Appl. Phys. Lett, 2010, 96: 162108. doi: 10.1063/1.3389494
-
[16]
Rhyee J. S, Cho E, Ahn K, Lee K. H, Lee S. M. Thermoelectric properties of bipolar diffusion effect on In4Se3-xTex compounds[J]. Appl. Phys. Lett, 2010, 97: 152104. doi: 10.1063/1.3493269
-
[17]
Lim Y. S, Jeong M, Seo W. S, Lee J. H, Park C. H, Sznajder M, Kharkhalis L. Y, Bercha D. M, Yang J. Condenson state and its effects on thermoelectric properties in In4Se3[J]. J. Phys. D: Appl. Phys, 2013, 46: 275304. doi: 10.1088/0022-3727/46/27/275304
-
[18]
Ahn K, Cho E, Rhyee J. S, Kim S. II, Lee S. M, Lee H. K. Effect of cationic substitution on the thermoelectric properties of In4-xMxSe2.95 compounds (M = Na, Ca, Zn, Ga, Sn, Pb; x = 0.1)[J]. Appl. Phys. Lett, 2011, 99: 102110. doi: 10.1063/1.3637053
-
[19]
Rhyee J. S, Ahn K, Lee K. H, Ji H. S, Shim J. H. Enhancement of the thermoelectric figure-of-merit in a wide temperature range in In4Se3-xCl0.03 bulk crystals[J]. Adv. Mater, 2011, 23: 2191-2194. doi: 10.1002/adma.v23.19
-
[20]
Kresse G, Furthmüller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Phys. Rev. B, 1996, 54: 11169-11186. doi: 10.1103/PhysRevB.54.11169
-
[21]
Perdew J. P, Wang Y. Accurate and simple analytic representation of the electron-gas correlation energy[J]. Phys. Rev. B, 1992, 45: 13244-13249. doi: 10.1103/PhysRevB.45.13244
-
[22]
Kresse G, Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method[J]. Phys. Rev. B, 1999, 59: 1758-1775.
-
[23]
Bl?chl P. E, Jepsen O, Andersen O. K. Improved tetrahedron method for Brillouin-zone integrations[J]. Phys. Rev. B, 1994, 49: 16223-16233. doi: 10.1103/PhysRevB.49.16223
-
[24]
Gibbs Z. M, Kim H. S, Wang H, Snyder G. J. Band gap estimation from temperature dependent Seebeck measurement-deviations from the 2e.S.maxTmax relation[J]. Appl. Phys. Lett, 2015, 106: 022112. doi: 10.1063/1.4905922
-
[25]
Goldsmid H. J, Sharp J. W. Estimation of the thermal band gap of a semiconductor from Seebeck measurements[J]. J. Electron. Mater, 1999, 28: 869-872. doi: 10.1007/s11664-999-0211-y
-
[26]
Lee M. H, Kim K. R, Rhyee J. S, Park S. D, Snyder G. J. High thermoelectric figure-of-merit in Sb2Te3/Ag2Te bulk composites as Pb-free p-type thermoelectric materials[J]. J. Mater. Chem. C, 2015, 3: 10494-10499. doi: 10.1039/C5TC01623A
-
[27]
Kim J. H, Kim M. J, Oh S, Rhyee J. S, Park S. D, Ahn D. Thermoelectric properties and chlorine doping effect of In4Pb0.01Sn0.03Se2.9Clx polycrystalline compounds[J]. Dalton Trans, 2015, 44: 3185-3189. doi: 10.1039/C4DT03432E
-
[28]
May A. F, Toberer E. S, Saramat A, Snyder G. J. Characterization and analysis of thermoelectric transport in n-type Ba8Ga16-xGe30+x[J]. Phys. Rev. B, 2009, 80: 125205. doi: 10.1103/PhysRevB.80.125205
-
[1]
-
Table 1. Formation Energies in eV at Each Site for Ag-doped In4Se3
In1 In2 In3 In4 Agint Ag 1.504 1.351 1.58 0.587 0.556 -

计量
- PDF下载量: 0
- 文章访问数: 7372
- HTML全文浏览量: 312