Abstract
In this paper, we systematically investigate the effect of hydrostatic pressure on the phononic and electronic transport properties of orthorhombic p-type GeSe using first-principles based Boltzmann transport equation approach. It is found that the lattice thermal conductivities along the a and c directions increase with pressure, whereas it experiences a decrease along the b direction. This anomalous pressure dependent lattice thermal conductivity is attributed to the combined effect of enhanced phonon group velocity and reduced phonon lifetime. Additionally, the optical phonon branches have remarkable contributions to the total lattice thermal conductivity. The electronic transport calculations indicate that the Seebeck coefficient undergoes a sign change from p-type to n-type along the a direction under pressure, and a dramatic enhancement of the power factor is observed due to the boost of electrical conductivity. The predicted ZT values along the a, b, and c directions are 1.54, 1.09, and 1.01 at 700âK and 8âGPa, respectively, which are about 14, 7.3, and 1.9 times higher than those at zero pressure at experimental carrier concentration of ~1018âcmâ3. Our study is expected to provide a guide for further optimization of the thermal and charge transport properties through hydrostatic pressure.
Similar content being viewed by others
Introduction
Thermoelectric materials, enabling a direct and reversible conversion between thermal energy and electricity based on the Seebeck effect and Peltier effect, play a vital role in the development of sustainable energy technologies. The wide applications of thermoelectric materials can effectively promote the energy utilization efficiency. The conversion efficiency of thermoelectric materials is characterized by the dimensionless figure of merit ZTâ=âSÏ2T/κ, which relies on the Seebeck coefficient S, electrical conductivity Ï, and thermal conductivity κ consisting of electronic and lattice contributions. Obviously, a high ZT value can be achieved by increasing the power factor (PF) SÏ2 or decreasing the thermal conductivity. However, owing to the coupling among Seebeck coefficient, electrical conductivity, and electronic thermal conductivity, it is extremely challenging to improve the Seebeck coefficient and electrical conductivity while reducing the electronic thermal conductivity simultaneously1. Among these thermoelectric parameters, lattice thermal conductivity is the only independent material property, thus the strategy of reducing lattice thermal conductivity is widely used to develop advanced thermoelectric materials2,3.
In recent years, SnSe, adopting a simple layered orthorhombic crystal structure, has attracted considerable attention from the scientific community as a promising thermoelectric material. Zhao et al.4 initially reported an unprecedented ZT of 2.6â±â0.3 at 923âK in undoped single crystal SnSe. Such outstanding thermoelectric performance arises from its intrinsically ultralow lattice thermal conductivity due to strong anharmonicity in bonding. As a member of IVâVI compound, GeSe possesses the same layered structure as SnSe5. Theoretical calculations have predicted that GeSe should be a potential candidate for high-efficient thermoelectric materials. Assuming the same hole concentration with SnSe in the range of 4âÃâ1019â6.5âÃâ1019âcmâ3, Hao et al.6 predicted an extremely high thermoelectric performance in hole-doped GeSe along the crystallographic b-axis with a ZT ranging from 0.8 at 300âK to 2.5 at 800âK, even superior to the SnSe compound. However, the optimization of carrier concentration by chemical doping for GeSe has proved challenging7,8. In order to increase its carrier concentration by doping, various elements, such as Cu, Ag and Na for p-type as well as Bi, Sb, La, As and I for n-type, have been tested, and only Ag-substitution on Ge site enabled a hole concentration up to ~1018âcmâ3. At a carrier concentration of ~1018âcmâ3, the maximal available ZT is only about 0.2 at 700âK for Ag-doped polycrystalline Ge0.79Ag0.01Sn0.2Se8.
Applying hydrostatic pressure has been shown as an effective way to modulate the electronic and phononic transport properties. Pressure has a strong influence on the electronic band structures, especially in the vicinity of the energy band gap. Since the Seebeck coefficient and electrical conductivity are strongly dependent on the electronic band structures, pressure can be used as a powerful tool to optimize the thermoelectric transport properties through the modification of the electronic band structures. Previous works have reported a remarkable enhancement of thermoelectric properties under high pressure either by pressure-enhanced electrical conductivity or by improved Seebeck coefficient9,10,11,12,13,14,15,16,17,18,19. Specifically, for the SnSe crystal, pressure significantly enhances the thermoelectric transport properties along all three directions in moderate temperature ranges due to the electrical transport boost15, which gives us a hint that the thermoelectric properties of GeSe may be effectively tailored by external pressure. Moreover, uniaxial or biaxial strain offer an alternative realization of pressure and have been shown to enhance the thermoelectric performance effectively20,21. In addition to the change of the electronic band structures, the phonon dispersion curves can be modified under high pressure, leading to the manipulation of phonon behaviors. Recently, the negative correlation between lattice thermal conductivity and pressure has been reported22,23,24. The decrease of the lattice thermal conductivity with the increase of pressure provides a routine to decouple electronic and phononic transport, which enlightens the strategy of minimizing the lattice thermal conductivity by pressure without deterioration in electrical conductivity for better thermoelectric performance. As a promising thermoelectric material with potential ZT value exceeding that of SnSe, the influence of pressure on the thermoelectric properties of GeSe is still unknown. Meanwhile, most previous theoretical works on the pressure dependence of the thermoelectric properties calculated the lattice thermal conductivity based on the Debye-Callaway model6,25, which ignores the contributions of optical phonon modes, and the effect of pressure on the phonon transport cannot be obtained exactly. Therefore, it is essential to conduct a systematical investigation on the variation of electronic and phononic transport properties in GeSe under pressure.
In this study, we investigate the pressure effect on the thermoelectric properties of GeSe based on the first-principles calculations. The lattice thermal conductivity is obtained through the iterative solution of the Boltzmann transport equation, while the electronic transport properties are calculated based on the semi-classical Boltzmann theory. It is observed that the lattice thermal conductivities along the three principal crystallographic directions vary anomalously with pressure. Owing to the pressure-induced enhancement of electrical conductivity, the overall ZT values of GeSe are improved significantly. The underlying mechanism is explored in detail by the phonon spectrum and electronic band structures analysis.
Results and Discussion
At low temperature and ambient pressure, the GeSe crystal adopts a hinge-like layered orthorhombic structure with a Pnma space group, as shown in Fig. 1(a). Within the plane of the two-atom-thick slab, each atom is covalently bonded to three nearest neighbors forming zigzag chains along the b axis and armchair chains along the c axis. These two-atom-thick slabs stack together along the a axis by weak van der Waals interactions. Figure 1(b) presents the variations of the lattice constants with pressure. As can be seen from Fig. 1(b), the fully relaxed lattice constants at 0âK and 0âGPa are aâ=â10.934âà , bâ=â3.878âà , and câ=â4.382âà , which show well agreement with experimental results of low-temperature Pnma GeSe crystal5,26,27. The lattice constants along all three directions decrease with increasing pressure, while the length reduction along the a direction is more significant than that along the b and c directions. The change of lattice constants with compression is consistent with previous x-ray-diffraction measurements27, which indicates that the length along the a-axis direction is more compressible than along the b and c directions. It is known that increasing temperature causes a first-order phase transition from orthorhombic structure to NaCl-type structure at ~900âK5,28. In addition, the high pressure stability of orthorhombic phase GeSe up to 9âGPa is verified by lattice dynamics computations29. Considering the phase-transition temperature and pressure above, in this work, we focus on the thermoelectric properties of low-temperature orthorhombic phase GeSe below 8âGPa.
The phonon dispersion curves obtained from the harmonic force constants are plotted in Fig. 2(a,b) for two representative pressures. The two transverse acoustic phonon branches (TA and TAâ²) and longitudinal acoustic phonon branch (LA) are highlighted by blue, green, and red lines along the three principal crystallographic directions. To accurately predict the phonon frequencies for polar materials, the effects of dipole-dipole interactions are incorporated into the dynamic matrix based on the Born effective charges and high-frequency dielectric constants, which can result in the LO-TO splitting (splitting between longitudinal and transversal optical phonon frequencies) presented at the Ð point. From a comparison of phonon dispersion curves at different pressures, we can clearly see that most phonon modes shift to higher frequencies indicating larger group velocity. Another obvious feature is that a pronounced frequency gap exists in the phonon dispersion at 0âGPa, which separates the phonon spectrum into two discernible regions, each with twelve phonon branches. Under pressure of 8âGPa, the frequency gap disappears due to the hardened low-lying optical phonon modes. This evolution of phonon dispersion with pressure agrees well with previous theoretical computation29. From the phonon densities of states shown in Fig. 2(c,d), this change can be seen more clearly. The pressure-induced close-up of frequency gap between low-lying and high-lying optical phonons and the mixing of acoustic phonons with low-lying optical phonons might bring significant modulation to the lattice thermal conductivity, rising from the strong coupling between acoustic and low-lying optical phonons and the increasing probability of optical phonons participating in three-phonon scattering process.
Calculated phonon dispersion curves of GeSe at (a) 0âGPa and (b) 8âGPa. The TA, TAâ², and LA phonon modes are highlighted by blue, green, and red lines along the crystallographic a, b, and c axes, respectively. The total and partial phonon densities of states of GeSe at (c) 0âGPa and (d) 8âGPa.
Figure 3(a) shows the temperature-dependent lattice thermal conductivity of GeSe along the a, b, and c axes under two typical pressures. From Fig. 3(a), it can be clearly seen that the intrinsic lattice thermal conductivity of GeSe is very low and and highly anisotropic. The lattice thermal conductivity decreases with the increase of temperature following the well-known κL ~ 1/T relationship. At 0âGPa, the calculated room-temperature lattice thermal conductivity of GeSe along the a, b, and c directions are 0.85, 4.85, and 2.26âW/mK, respectively, and κL reduces to 0.38, 2.15, and 1.01âW/mK along the corresponding directions as the temperature rises to 700âK. The magnitude sequence κL,bâ>âκL,câ>âκL,a is consistent with the finding in Pnma phase SnSe4,15,30, which is attributed to the different phonon group velocities along the three directions, as evidenced by the phonon dispersion. The calculated room-temperature lattice thermal conductivity at 0âGPa is higher than the theoretical predicted value of ~0.64âW/mK along the b direction at 300âK6. This is understandable that the lattice thermal conductivity obtained by Hao et al.6 is based on the Debye-Callaway model, this simple model has low computational cost but does not take into account the contributions from optical modes. As discussed later, the optical contributions are as high as those of acoustic modes and cannot be overlooked, thus, the lattice thermal conductivity is significantly underestimated by the Debye-Callaway model25. However, our results of the average lattice thermal conductivity is comparable to the experimental values (~1.8âW/mK)8 for polycrystalline GeSe samples at room temperature. At a pressure of 8âGPa, the calculated lattice thermal conductivity values along the a and c directions are larger than those at 0âGPa over the considered temperature range, whereas a surprising reduction in the lattice thermal conductivity along the b direction is observed. Specifically, the lattice thermal conductivities along the a and c directions increase by a factor of 2.31 and 1.43, respectively, while the lattice thermal conductivity along the b direction reduces by 15%. As far as we know, Meng et al.31 also reported an increase in the cross-plane lattice thermal conductivity along with a decrease in the in-plane lattice thermal conductivity in multilayer MoS2 under hydrostatic pressure by first-principles calculations. This anisotropic evolution of lattice thermal conductivity with pressure is the result of the competition between the enhanced group velocity and reduced phonon lifetime, which will be discussed in more detail below. Moreover, an important factor that must be noted is the limitations of the method adopted to calculate lattice thermal conductivity. The inputs of second-order and third-order IFCs are determined from zero-temperature ground-state systems, which actually are temperature dependent. Romero et al.32 have pointed out the necessity of the complete anharmonic behavior at finite-temperature structures when dealing with strongly anharmonic systems. Currently, it is still very challenging to treat strong anharmonic effects when higher order scattering can not be ignored32,33. As a result, our calculations only give the qualitative trends of temperature- and pressure-dependent lattice thermal conductivity.
(a) The lattice thermal conductivity of GeSe along the a, b, and c directions at two representative pressures from 300 to 700âK. The percentage contributions to the lattice thermal conductivity of GeSe from three acoustic phonon branches and total optical phonon branches along the (b) a direction, (c) b direction, and (d) c direction under different pressures.
In Fig. 3(bâd), we further show the percentage contributions from different phonon branches to the total lattice thermal conductivity along the a, b, and c directions at different pressures. It can be seen that optical phonon modes make a non-negligible contribution to the total lattice thermal conductivity, especially along the b and c directions. At a pressure of 0âGPa, optical phonons contribute ~40% to the total lattice thermal conductivity along the a direction within studied temperature range, while their contributions to the total in-plane lattice thermal conductivities reach as much as ~65%. The optical phonon contributions for GeSe are comparable to those of SnSe but larger as compared with another group IV-VI compounds PbSe and PbTe with rocksalt crystal structures30,34,35, in which the contributions of optical phonons remain about 25% for PbSe and 22% for PbTe between 300 and 700âK. Our calculations indicate that optical phonons are not negligible and play a dominant role in heat conduction of GeSe. Such high optical phonon contributions is the reflection of more optical phonon branches, because GeSe has a large primitive cell with eight basis atoms. On the other hand, the group velocities of optical phonons are large due to more dispersive phonon branches as shown in Fig. 2. At 300âK and 8âGPa, the proportions of optical phonons in the total lattice thermal conductivity along the a and c directions rise to ~44% and ~68%, respectively, but by way of comparison, the contributions of optical phonons decrease to ~62% along the b direction. The results indicate that it is the variation of optical phonon modes that governs the abnormal pressure-dependent lattice thermal conductivity along the three directions.
Next, to reveal the underlying mechanisms governing the phonon transport of GeSe, we analyze the pressure dependence of mode-level phonon properties as a function of frequency at room temperature in terms of group velocity, phonon lifetime, three-phonon phase space, and mode Grüneisen parameter. First, the changes in group velocity along the three directions are plotted (see Supplementary Fig. S1). It is can be seen that, along the a direction, the group velocities of acoustic phonons are obviously larger than those of most optical phonons, as evidenced by phonon dispersion characterized by flatted optical branches along the Î-X direction. With regard to the pressure effect, pressure enhances the group velocities in all three directions, which facilitates phonon transport in GeSe. However, the variation of group velocity along the b direction is minor as compared to those along the a and c directions.
The frequency dependent phonon lifetimes at two typical pressures are displayed in Fig. 4. It is found that most phonon lifetimes are on the order of several picoseconds, which is short enough for low phonon thermal conductivity. Moreover, many optical phonon lifetimes are comparable with those of acoustic phonons. In corporation with large group velocities and broad phonon densities of states, the contribution of optical phonons to the total lattice thermal conductivity overwhelms that of each acoustic phonon branches. Through the comparison of phonon lifetime between 0âGPa and 8âGPa, it can be found that there is a slight decrease in phonon lifetime with pressure, especially for high frequency optical phonons. The phase space for the three-phonon scattering processes, which provides a measure of the available channels for three-phonon processes, is compared between 0âGPa and 8âGPa (see Supplementary Fig. S2). Larger phase space for three-phonon processes implies a larger number of available scattering channels. Consequently, there is an inverse relationship between phase space and the intrinsic lattice thermal conductivity. Three-phonon scattering processes can be further dissociated into absorption and emission three-phonon processes (see Supplementary Fig. S3). For low frequency phonons, absorption channels like Aâ+âAâââA and Aâ+âOâââO are dominant processes, while emission channels like OâââAâ+âA and OâââAâ+âO are anticipated to play an important role at high frequency, where âAâ and âOâ refer to acoustic and optical phonon modes, respectively. As we can see, the changing trend of three-phonon phase space with pressure is consistent with the variation of phonon lifetime, suggesting the governing role affecting the phonon lifetime. In the inset of Supplementary Fig. S2, the mode-dependent phase space normalized to the zero pressure values is depicted. In the region of 1.0~3.0âTHz and 5.0~6.0âTHz, the scattering channels are significantly enhanced at high pressure, while it is restricted for the rest frequency region. The main reason for the increased scattering channels is ascribed to more dispersive optical phonons and the close-up of frequency gap between low-lying and high-lying optical phonons. At the frequency ranging from 1.0âTHz to 3.0âTHz, acoustic phonon and low-frequency optical phonon can easily emerge into one high-frequency optical phonon due to the dispersion of low-lying optical phonon modes considering the energy and momentum conservation criterion. While at the frequency ranging from 5.0âTHz to 6.0âTHz, there are more chances for a high-frequency optical phonon scattering into two low-frequency phonons. As a result, the phonon lifetime shows a remarkable decrease at high pressure in the corresponding frequency range where the three-phonon scattering channels are enhanced.
The mode Grüneisen parameters of GeSe are calculated (see Supplementary Fig. S4), which can be used to measure the strength of anharmonic phonon-phonon scattering. We can see that the effect of pressure on the mode Grüneisen parameters is not significant with the exception of slight increase for high frequency optical phonons. The mode weighted accumulative Grüneisen parameter is calculated based on \(\bar{\gamma }={\sum }_{\lambda }{c}_{{\rm{ph}},\lambda }{\gamma }_{\lambda }/{\sum }_{\lambda }{c}_{{\rm{ph}},\lambda }\), where cph is phonon volumetric specific heat, γ is mode Grüneisen parameter, and λ denotes a phonon mode with a branch index s and a wave vector q. The calculated mode averaged Grüneisen parameters are 1.59, 0.92, and 1.16 along the a, b, and c directions at 300âK and 0âGPa, respectively. Previously, Hao et al.6 has reported higher Grüneisen parameters for GeSe \({\bar{\gamma }}_{a}=4.4\), \({\bar{\gamma }}_{b}=1.9\), and \({\bar{\gamma }}_{c}=2.6\), which are averaged over only the acoustic phonons. Our results of Grüneisen parameters are comparable with those of SnSe (\({\bar{\gamma }}_{a}=2.12\), \({\bar{\gamma }}_{b}=1.55\), \({\bar{\gamma }}_{c}=1.66\)) and SnS (\({\bar{\gamma }}_{a}=2.17\), \({\bar{\gamma }}_{b}=1.44\), \({\bar{\gamma }}_{c}=1.55\)) obtained by the same method34. At 8âGPa, the mode averaged Grüneisen parameters of GeSe are \({\bar{\gamma }}_{a}=1.78\), \({\bar{\gamma }}_{b}=1.12\), and \({\bar{\gamma }}_{c}=1.38\), respectively, which are slightly larger compared with the zero pressure values. As suggested by the semi-empirical theory, the phonon lifetime due to phonon-phonon Umklapp scattering is inversely proportional to the square of Grüneisen parameter36. Therefore, the anharmonic phonon scattering is enhanced by large Grüneisen parameter, which partially contributes to the reduced phonon lifetime at high pressure.
To evaluate the relative importance of the phonon specific heat (cph), group velocity (Ï g), and phonon lifetime (Ï) to the lattice thermal conductivity, \({\kappa }_{{\rm{L}},\alpha }={\sum }_{\lambda }{c}_{{\rm{ph}},\lambda }{\upsilon }_{\alpha ,\lambda }^{2}{\tau }_{\lambda }\), we substitute these values of uncompressed state with their respective values at a pressure of 8âGPa. For example, for the case of cph replaced thermal conductivity, the phonon specific heat at 8âGPa is adopted to calculate the total thermal conductivity while fixing group velocity and phonon lifetime at 0âGPa. The replaced lattice thermal conductivity results of GeSe along different directions are shown in Fig. 5, together with unsubstituted values at 0âGPa and 8âGPa. It can be seen clearly that the reduction of phonon lifetime leads to a lower lattice thermal conductivity along all directions, while the enhanced phonon specific heat brings out the augmentation of the lattice thermal conductivity. By comparison, the positive contribution of group velocity to the a-axis and c-axis lattice thermal conductivity is more significant than that of phonon specific heat. However, the group velocity replaced b-axis lattice thermal conductivity shows a very little decrease under high pressure. This is consistent with our previous discussions about group velocity and phonon lifetime. Therefore, the increased group velocity and phonon specific heat along the a and c directions surpass the reduction of phonon lifetime, which finally results in enhanced lattice thermal conductivity along these two directions under pressure, whereas the insensitive pressure dependent group velocity and depressed phonon lifetime cancel out the increased phonon specific heat, leading to a lower lattice thermal conductivity along the b direction at 8âGPa.
Now, we turn our attention to the effect of pressure on the electronic transport properties in GeSe. Applying the semiclassical Boltzmann transport theory, we calculate the thermoelectric transport properties of GeSe under different pressures. The experimentally measured Hall transport properties of Ag-doped GeSe below 700âK suggest that the carrier concentration is ~1âÃâ1018âcmâ3â8. Therefore, in our transport calculations we consider a hole concentration of 1âÃâ1018âcmâ3.
In the output of BoltzTraP code, the electrical conductivity is scaled by the constant relaxation time of unity (1âfs). To overcome this problem and have a comparative assessment of the effect of pressure, we calculate the electron relaxation time approximated by the the single parabolic band (SPB) model at given temperature and carrier concentration. Differently from MoS2 and Mg2Sn in which the valence band maximum (VBM) is at the Î point17,19, in our model system the VBM is not located at the Î point and there are more than one deformation potentials describing the pressure/strain behavior of the electronic band structure, thus we approximately adopt the deformation potential of the VBM along the Î-Z direction for the electron relaxation time calculations. The deformation potential is estimated based on the energy change of the VBM by varying the cell volumes instead of a stretch and a contraction along different axes. The obtained deformation potential at the VBM is comparable with that of SnSe34. Another approximation we made is that the same relaxation time is used for the conduction bands as for the valence bands in the calculation. Figure 6(a) shows the evaluated electron relaxation time under pressures of 0âGPa and 8âGPa by the SPB model. Apparently, the electron relaxation time is anisotropic and inversely proportional to temperature. At high temperatures, the electron relaxation time is in the order of several tens of femtoseconds. Hao et al.6 found that the electron relaxation time along the b-axis decreases from 27 to 4âfs as the temperature rises from 300 to 800âK for the hole concentration in the range of 4âÃâ1019â6.5âÃâ1019âcmâ3. Our results seem reasonable compared to the values in ref.6 considering the lower hole concentration we adopted. Under compression, the electron relaxation time increases due to enhanced longitudinal sound velocity and reduced effective mass. It should be noted that the used SPB model only considers acoustic phonon scattering, and previous works had demonstrated the dominant role of longitudinal optical (LO) phonon scattering, such as in n-type PbTe37,38,39. We roughly evaluate the effect of LO phonon scattering based on the first-principles results of electron-phonon scattering in SnSe40, whose electron-phonon scattering mechanisms may be analogous to GeSe due to the structural and chemical similarities (see Supplementary Information for details). Assuming 60% of the total scattering rate is contributed by the LO phonons in GeSe, the electron relaxation times will reduce to 11âfs, 26âfs, and 44âfs along the a, b, and c directions at 300âK and 0âGPa. The electron relaxation times accounting for LO phonon scattering are still in the same order as the results of the SPB model. We also plot the electrical conductivity and power factor scaled by the electron relaxation time as shown in Fig. S5. The changing trends of Ï/Ïe and S2Ï/Ïe are in accordance with the results based on the SPB model. Therefore, although the SPB model relying only on acoustic phonon scattering is certainly not rigorous, it is easy to apply and can to some extent reveal the effect of temperature and pressure on the electron relaxation time and thermoelectric transport properties.
Figure 6(b) plots the Seebeck coefficient as a function of temperature at different pressures. At 0âGPa, the Seebeck coefficients of the three directions increase with increasing the temperature and the values are all above +450 μV/K, suggesting the p-type property. Under a pressure of 8âGPa, the absolute values of Seebeck coefficients of the three directions are slightly lower than those at zero pressure and decrease with the increase of temperature. Interestingly, the in-plane direction remains positive Seebeck coefficients (p-type), while the Seebeck coefficient along the a direction exhibits a sign change, implying a carrier change from p-type to n-type under pressure. In Fig. 6(c), we show the temperature dependent electrical conductivity under pressure. Opposite to negative pressure trend in Seebeck coefficient, the electrical conductivity is significantly enhanced under pressure. The decreasing Seebeck coefficient and increasing electrical conductivity with temperature at 8âGPa are due to the bipolar transport and originate from the thermal excitation of carriers, which is the result of smaller band gap under pressure as discussed below. Combining the results of Seebeck coefficient and electrical conductivity, the corresponding power factor can be assessed. As shown in Fig. 6(d), because of the boost in electrical conductivity, the power factor has a pronounced increase under pressure. For example, at 0âGPa and 700âK, the calculated power factors along the a, b, and c directions are observed to be about 0.65, 4.48, and 7.55 μW cmâ1 Kâ2, respectively, and increase to 54.15, 61.45, and 41.57 μW cmâ1 Kâ2 at 8âGPa and 700âK along the corresponding directions.
The large enhancement in thermoelectric properties of GeSe can be understood in terms of the band structure of this material. In Fig. 7, we present the band structures of GeSe at different pressures. The Fermi level is set at 0âeV and is indicated by the black dashed lines. The effect of spin-orbit coupling (SOC) is ignored. Although many works41,42 have demonstrated the significant influence of SOC on thermoelectric properties by causing band splitting of the electronic structures, no spin-orbit splitting is expected in bulk group-IV monochalcogenides, because the inversion symmetry together with time-reversal symmetry require that the bands for the two spins are degenerate43. Therefore, the SOC barely changes the band shape of GeSe, and electronic transport properties should remain unchanged. From the band structure at 0âGPa, the valence band maximum lies in the Î-Z direction while the conduction band minimum (CBM) locates at the Î point, yielding a indirect band gap of 0.97âeV, which is slightly lower than the experimental value of 1.1âeV44,45. With the increase of pressure from 0 to 8âGPa, there is a reduction in band gap along with the downward trend of conduction band at the Î point. At higher pressures we notice that the valence band along the Î-Y direction moves toward the Fermi level and the flattening of this band is reducing, which is an indication of a light band. The high dispersion in the valence band along the Î-Y direction results in the low density of states around the Fermi level but high carrier mobility, thus leading to the decrease of the Seebeck coefficient along with the increase of the electrical conductivity. Along the Î-X direction, the valence bands are all below â0.5âeV while the conduction bands become highly dispersive at high pressures. Since the reduction of band gap under pressure, this makes thermal excitation easier and the conduction bands along the Î-X direction are occupied with thermally excited electrons, thus displaying n-type property along the a direction in the real space.
Using the above calculated Seebeck coefficient, electrical conductivity, and lattice thermal conductivity, we could estimate the figure of merit by \(ZT=\frac{{S}^{2}}{{L}_{{\rm{o}}}+({\kappa }_{{\rm{L}}}/\sigma T)}\), where Loâ=â2.45âÃâ10â8âW Ω Kâ2 is the Lorentz number. The obtained ZT values along the three directions at different pressures are shown in Fig. 8 as a function of temperature. The ZT values increase as the temperature rises from 300 to 700âK and increasing the pressure causes the enhancement in ZT values in all directions, which are attributed to the temperature-dependent and pressure-dependent electrical conductivity. At 700âK and ambient pressure, the computed ZT values along the a, b, and c directions are 0.11, 0.15, and 0.53, respectively, and these results are reasonable as compared to previous experimental study, which reported a ZT value of ~0.2 for polycrystalline GeSe at the same temperature and carrier concentration8. At 700âK and 8âGPa, the estimated ZT values along the a, b, and c directions reach 1.54, 1.09, and 1.01, respectively. By the application of hydrostatic pressure, the p-type ZT along the b and c directions increase by 7.3 and 1.9 times, while the pressure-induced n-type ZT along the a direction increases by a factor of 14, which is much more prominent than the rise in the in-plane direction. For comparison, Guo et al.46 reported the thermoelectric properties of monolayer GeSe with a maximum n-type ZT of 0.8 at 900âK based on the calculated lattice thermal conductivities from Qin et al.47. In addition, a much higher ZT of 1.99 is predicted for monolayer GeSe48. It can be seen that the ZT values of bulk GeSe at high pressure are comparable to those of monolayer counterpart.
Conclusions
In summary, we have performed detailed first-principles calculations to investigate the effect of hydrostatic pressure on the electronic and phononic transport properties of GeSe. The lattice dynamics results show that the pressure has profound influence on the phonon dispersion, which exhibits decreasing frequency gap and highly dispersive phonon branches at higher pressures. The results from iterative solutions of phonon Boltzmann transport equation suggest that the lattice thermal conductivities along the three principal crystallographic directions vary differently. The lattice thermal conductivities along the a and c directions increase with pressure, while an opposite trend is observed in the b direction. A detailed analysis on the phonon related properties reveals that there is an evident increase in both phonon specific heat and group velocity under compression with the exception of pressure insensitive group velocity along the b direction. The phonon lifetime is significantly reduced at higher pressures due to the enhancement of three-phonon scattering channels and lattice anharmonicity. Along the a and c directions, the magnitude of pressure-enhanced phonon specific heat and group velocity overwhelm the reduced phonon lifetime, thus resulting in higher lattice thermal conductivities under pressure. For the lattice thermal conductivity along the b direction, the largely reduced phonon lifetime in combination with the pressure insensitive group velocity lead to a slight decrease at 8âGPa. From the discussion of the contribution from each phonon branch, it is found that optical phonon branches have considerable contributions to the total lattice thermal conductivity and their pressure variations dominate the diverse pressure-dependent lattice thermal conductivity along the three directions. The results of electronic transport properties demonstrate that the power factors are greatly enhanced due to the boost of electrical conductivity. At 8âGPa and 700âK, the estimated ZT values for p-type doping along the b and c directions are about 7.3 and 1.9 times higher than those at ambient pressure. The predicted highest ZT value for n-type doping along the a direction is greater than 1.5 at experimental carrier concentration, which is enhanced by a factor of 14 as compared with the ZT value at 0âGPa. Our calculations show that there is plenty of room for improvement of the thermoelectric performance in layered orthorhombic GeSe through the modulation of hydrostatic pressure.
Computational methods
All the first-principles calculations are carried out based on density functional theory (DFT) as implemented in the Vienna ab-initio Simulation Package (VASP)49. The projector-augmented wave (PAW) pseudopotentials50 are used to describe the interaction among atoms and the generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE)51 form is chosen as the exchange-correlation functional. The kinetic energy cutoff of the plane-wave function is set to 500âeV for all the calculations. The convergent criterion for the total energy difference between two successive self-consistency steps is 10â8âeV and the primitive cell is fully relaxed until the maximum force acting on each atom is less than 10â8âeVâà â1. A well-converged Monkhorst-Pack52 k-point grid of 5âÃâ15âÃâ13 is used to sample the first irreducible Brillouin zone. In general, GGA usually overestimates the lattice constants of weakly bonded layered solids. To include the van der Waals interaction in layered GeSe, the non-local van der Waals density functional termed as optB86b-vdW53,54 is taken into account. Since an accurate estimation of band gap is important for the prediction of thermoelectric transport properties, the Tran-Blaha modified Becke-Johnson (TB-mBJ)55 exchange-correlation potential is used to calculate the electronic band structure of bulk GeSe.
The BoltzTraP package56 is utilized to calculate the Seebeck coefficient and electrical conductivity of bulk GeSe under different pressures. The BoltzTraP code is based on semi-classical Boltzmann theory in conjunction with rigid band and constant relaxation time approximations, and has been successfully applied to various types of thermoelectric materials57,58,59,60. The electrical conductivity tensors Ï and Seebeck coefficient tensors S as a function of temperature T and chemical potential μ are expressed as
where α and β are Cartesian indices, Ω is the volume of unit cell, ε is the band energy, f is the Fermi-Dirac distribution function of carriers, âαβ(ε) is the transport distribution function and defined as
where the summation is over all bands i and over all k grid, Ï is the relaxation time, Ï is the band velocity. It is noted that the most simple approximation that the relaxation time is independent of both i and k is adopted for the results. From the above equations, it can be seen that the calculation accounts for both valence band and conduction band. Since \(\partial f(\varepsilon )/\partial \varepsilon \) has a finite value only in an energy window about a few kBT around the chemical potential μ, the carriers in this energy range contribute most to the transport properties. To get well-converged electronic transport quantities, a finer 15âÃâ45âÃâ39 k-mesh is selected.
For the phononic transport properties, an iterative self-consistent method is employed for solving the phonon Boltzmann transport equation to calculate the lattice thermal conductivity using the ShengBTE package61 with the inputs of second-order (harmonic) and third-order (anharmonic) interatomic force constants (IFCs). Both the harmonic and anharmonic IFCs are computed by the finite difference method. Several independent atomic displacements created by the PHONOPY62 program within 2âÃâ6âÃâ5 supercells are used for the calculation of harmonic IFCs. In the calculation of the anharmonic IFCs, a 2âÃâ3âÃâ3 supercell is constructed in consideration of the computational cost and the cutoff radius is considered up to the 27th nearest neighboring shell for the sufficient involvement of the long-range interactions. Finally, the lattice thermal conductivity is integrated on a well-converged 9âÃâ27âÃâ23 phonon q grid.
To evaluate the electron relaxation time more accurately at given temperature and carrier concentration, the SPB model is employed. According to the SPB model, the electron relaxation time is calculated by
where ħ is the reduced Planck constant, Ï is the mass density, Ï l is the longitudinal sound velocity, kB is the Boltzmann constant, and T is the absolute temperature. Ed is the deformation potential obtained by \({E}_{{\rm{d}}}=\partial {E}_{{\rm{v}}}/\partial (V/{V}_{0})\), which describes the energy change at the band valley with respect to the variation of volume. A series of volumes varying from 0.985V0 to 1.015V0 with a step of 0.005V0 are used to calculate the band structure, where V0 is the equilibrium volume. The calculated deformation potentials of GeSe are â11.04âeV at 0âGPa and â13.54âeV at 8âGPa for the valence band maximum, respectively. m* is the effective mass at the band extrema, which is calculated using finite difference method63. The obtained effective masses for the valence band maximum are \({m}_{kx}^{\ast }\)â=â0.58 me, \({m}_{ky}^{\ast }\)â=â0.47 me, and \({m}_{kz}^{\ast }\)â=â0.21 me at 0âGPa, respectively, and are reduced to 0.56 me, 0.19 me, and 0.15 me at a pressure of 8âGPa along the corresponding directions, where me is the free electron mass. Fx(η) is the Fermi integral given by
where ηâ=âE/kBT is the reduced chemical potential. For the given T and η, carrier concentration n is determined by
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
References
Snyder, G. J. & Toberer, E. S. Complex thermoelectric materials. Nat. Mater. 7, 105 (2008).
Chen, Z., Zhang, X. & Pei, Y. Manipulation of phonon transport in thermoelectrics. Adv. Mater. 30, 1705617 (2018).
Toberer, E. S., Zevalkink, A. & Snyder, G. J. Phonon engineering through crystal chemistry. J. Mater. Chem. 21, 15843â15852 (2011).
Zhao, L.-D. et al. Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals. Nature 508, 373â377 (2014).
Sist, M. et al. High-temperature crystal structure and chemical bonding in thermoelectric germanium selenide (GeSe). Chem. Eur. J. 23, 6888â6895 (2017).
Hao, S., Shi, F., Dravid, V. P., Kanatzidis, M. G. & Wolverton, C. Computational prediction of high thermoelectric performance in hole doped layered GeSe. Chem. Mater. 28, 3218â3226 (2016).
Shaabani, L., Aminorroaya-Yamini, S., Byrnes, J., Akbar Nezhad, A. & Blake, G. R. Thermoelectric performance of Na-Doped GeSe. ACS Omega 2, 9192â9198 (2017).
Zhang, X. et al. Thermoelectric properties of GeSe. J. Materiomics 2, 331â337 (2016).
Ovsyannikov, S. V. & Shchennikov, V. V. Pressure-tuned colossal improvement of thermoelectric efficiency of PbTe. Appl. Phys. Lett. 90, 122103 (2007).
Meng, J. F. et al. Multifold enhancement of the thermoelectric figure of merit in p-type BaBiTe3 by pressure tuning. J. Appl. Phys. 90, 2836â2839 (2001).
Polvani, D. A., Meng, J. F., Shekar, N. V. C., Sharp, J. & Badding, J. V. Large improvement in thermoelectric properties in pressure-tuned p-type Sb1.5Bi0.5Te3. Chem. Mater. 13, 2068â2071 (2001).
Ovsyannikov, S. V. et al. Enhanced power factor and high-pressure effects in (Bi, Sb)2(Te, Se)3 thermoelectrics. Appl. Phys. Lett. 106, 143901 (2015).
Ovsyannikov, S. V. et al. Giant improvement of thermoelectric power factor of Bi2Te3 under pressure. J. Appl. Phys. 104, 053713 (2008).
Ibarra-Hernandez, W., Verstraete, M. J. & Raty, J.-Y. Effect of hydrostatic pressure on the thermoelectric properties of Bi2Te3. Phys. Rev. B: Condens. Matter Mater. Phys. 90, 245204 (2014).
Zhang, Y., Hao, S., Zhao, L., Wolverton, C. & Zeng, Z. Pressure induced thermoelectric enhancement in SnSe crystals. J. Mater. Chem. A 4, 12073â12079 (2016).
Miao, N., Zhou, J., Sa, B., Xu, B. & Sun, Z. Pressure-induced semimetal-semiconductor transition and enhancement of thermoelectric performance in α-MgAgSb. Appl. Phys. Lett. 108, 213902 (2016).
Guo, H., Yang, T., Tao, P., Wang, Y. & Zhang, Z. High pressure effect on structure, electronic structure, and thermoelectric properties of MoS2. J. Appl. Phys. 113, 013709 (2013).
Javed, Y., Rafiq, M. A. & Ahmed, N. Pressure-induced changes in the electronic structure and enhancement of the thermoelectric performance of SnS2: a first principles study. RSC Adv. 7, 38834â38843 (2017).
Guo, S.-D. & Wang, J.-L. Pressure enhanced thermoelectric properties in Mg2Sn. RSC Adv. 6, 31272â31276 (2016).
Guo, S.-D. Biaxial strain tuned thermoelectric properties in monolayer PtSe2. J. Mater. Chem. C 4, 9366â9374 (2016).
Lv, H. Y., Lu, W. J., Shao, D. F., Lu, H. Y. & Sun, Y. P. Strain-induced enhancement in the thermoelectric performance of a ZrS2 monolayer. J. Mater. Chem. C 4, 4538â4545 (2016).
Ouyang, T. & Hu, M. Competing mechanism driving diverse pressure dependence of thermal conductivity of XTe (Xâ=âHg, Cd, and Zn). Phys. Rev. B: Condens. Matter Mater. Phys. 92, 235204 (2015).
Luo, Y., Wang, J., Li, Y. & Wang, J. Giant phonon anharmonicity and anomalous pressure dependence of lattice thermal conductivity in Y2Si2O7 silicate. Sci. Rep. 6, 29801 (2016).
Lindsay, L., Broido, D. A., Carrete, J., Mingo, N. & Reinecke, T. L. Anomalous pressure dependence of thermal conductivities of large mass ratio compounds. Phys. Rev. B: Condens. Matter Mater. Phys. 91, 121202 (2015).
Zhang, Y. First-principles DebyeâCallaway approach to lattice thermal conductivity. J. Materiomics 2, 237â247 (2016).
Okazaki, A. The crystal structure of germanium selenide GeSe. J. Phys. Soc. Jpn. 13, 1151â1155 (1958).
Onodera, A., Sakamoto, I., Fujii, Y., Moâri, N. & Sugai, S. Structural and electrical properties of GeSe and GeTe at high pressure. Phys. Rev. B: Condens. Matter Mater. Phys. 56, 7935â7941 (1997).
Wiedemeier, H. & Siemers, P. A. The thermal expansion and high temperature transformation of GeSe. Z. Anorg. Allg. Chem. 411, 90â96 (1975).
Deringer, V. L., Stoffel, R. P. & Dronskowski, R. Vibrational and thermodynamic properties of GeSe in the quasiharmonic approximation. Phys. Rev. B: Condens. Matter Mater. Phys. 89, 094303 (2014).
Carrete, J., Mingo, N. & Curtarolo, S. Low thermal conductivity and triaxial phononic anisotropy of SnSe. Appl. Phys. Lett. 105, 101907 (2014).
Meng, X. et al. Thermal conductivity enhancement in MoS2 under extreme strain. Phys. Rev. Lett. 122, 155901 (2019).
Romero, A. H., Gross, E. K. U., Verstraete, M. J. & Hellman, O. Thermal conductivity in PbTe from first principles. Phys. Rev. B: Condens. Matter Mater. Phys. 91, 214310 (2015).
Lindsay, L. First principles peierls-boltzmann phonon thermal transport: A topical review. Nanoscale and Microscale Thermophysical Engineering 20, 67â84 (2016).
Guo, R., Wang, X., Kuang, Y. & Huang, B. First-principles study of anisotropic thermoelectric transport properties of IVâVI semiconductor compounds SnSe and SnS. Phys. Rev. B: Condens. Matter Mater. Phys. 92, 115202 (2015).
Tian, Z. et al. Phonon conduction in PbSe, PbTe, and PbTe1âxSex from first-principles calculations. Phys. Rev. B: Condens. Matter Mater. Phys. 85, 184303 (2012).
Callaway, J. Model for lattice thermal conductivity at low temperatures. Phys. Rev. 113, 1046â1051 (1959).
Cao, J., Querales-Flores, J. D., Murphy, A. R., Fahy, S. & SaviÄ, I. Dominant electron-phonon scattering mechanisms in n-type PbTe from first principles. Phys. Rev. B: Condens. Matter Mater. Phys. 98, 205202 (2018).
Song, Q., Liu, T.-H., Zhou, J., Ding, Z. & Chen, G. Ab initio study of electron mean free paths and thermoelectric properties of lead telluride. Materials Today Physics 2, 69â77 (2017).
Bilc, D. I., Mahanti, S. D. & Kanatzidis, M. G. Electronic transport properties of PbTe and AgPbmSbTe2+m systems. Phys. Rev. B: Condens. Matter Mater. Phys. 74, 125202 (2006).
Ma, J., Chen, Y. & Li, W. Intrinsic phonon-limited charge carrier mobilities in thermoelectric SnSe. Phys. Rev. B: Condens. Matter Mater. Phys. 97, 205207 (2018).
Guo, S.-D. & Wang, J.-L. Spinâorbital coupling effect on the power factor in semiconducting transition-metal dichalcogenide monolayers. Semicond. Sci. Technol. 31, 095011 (2016).
Guo, S.-D. Importance of spinâorbit coupling in power factor calculations for half-Heusler ANiB (Aâ=âTi, Hf, Sc, Y; Bâ=âSn, Sb, Bi). J. Alloys Compd. 663, 128â133 (2016).
Gomes, L. C. & Carvalho, A. Phosphorene analogues: Isoelectronic two-dimensional group-IV monochalcogenides with orthorhombic structure. Phys. Rev. B: Condens. Matter Mater. Phys. 92, 085406 (2015).
Vaughn, D. D., Patel, R. J., Hickner, M. A. & Schaak, R. E. Single-crystal colloidal nanosheets of GeS and GeSe. J. Am. Chem. Soc. 132, 15170â15172 (2010).
Elkorashy, A. M. Photoconductivity in germanium selenide single crystals. Phys. Status Solidi B 152, 249â259 (1989).
Guo, S.-D. & Wang, Y.-H. Thermoelectric properties of orthorhombic group IVâVI monolayers from the first-principles calculations. J. Appl. Phys. 121, 034302 (2017).
Qin, G. et al. Diverse anisotropy of phonon transport in two-dimensional group IVâVI compounds: A comparative study. Nanoscale 8, 11306â11319 (2016).
Shafique, A. & Shin, Y.-H. Thermoelectric and phonon transport properties of two-dimensional IVâVI compounds. Sci. Rep. 7, 506 (2017).
Kresse, G. & Furthmuller, 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).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B: Condens. Matter Mater. Phys. 59, 1758â1775 (1999).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865â3868 (1996).
Monkhorst, H. J. & Pack, J. D. Special points for brillouin-zone integrations. Phys. Rev. B: Condens. Matter Mater. Phys. 13, 5188â5192 (1976).
Klimeš, J., Bowler, D. R. & Michaelides, A. Van der Waals density functionals applied to solids. Phys. Rev. B: Condens. Matter Mater. Phys. 83, 195131 (2011).
JiÅÃ, K., David, R. B. & Angelos, M. Chemical accuracy for the van der Waals density functional. J. Phys.: Condens. Matter 22, 022201 (2010).
Tran, F. & Blaha, P. Accurate band gaps of semiconductors and insulators with a semilocal exchange-correlation potential. Phys. Rev. Lett. 102, 226401 (2009).
Madsen, G. K. H. & Singh, D. J. BoltzTraP. A code for calculating band-structure dependent quantities. Comput. Phys. Commun. 175, 67â71 (2006).
Lee, M.-S. & Mahanti, S. D. Validity of the rigid band approximation in the study of the thermopower of narrow band gap semiconductors. Phys. Rev. B: Condens. Matter Mater. Phys. 85, 165149 (2012).
Parker, D. & Singh, D. J. High-temperature thermoelectric performance of heavily doped PbSe. Phys. Rev. B: Condens. Matter Mater. Phys. 82, 035204 (2010).
Ong, K. P., Singh, D. J. & Wu, P. Analysis of the thermoelectric properties of n-type ZnO. Phys. Rev. B: Condens. Matter Mater. Phys. 83, 115110 (2011).
Takagiwa, Y., Pei, Y., Pomrehn, G. & Jeffrey Snyder, G. Validity of rigid band approximation of PbTe thermoelectric materials. APL. Materials 1, 011101 (2013).
Li, W., Carrete, J., A. Katcho, N. A. & Mingo, N. ShengBTE: A solver of the Boltzmann transport equation for phonons. Comput. Phys. Commun. 185, 1747â1758 (2014).
Togo, A. & Tanaka, I. First principles phonon calculations in materials science. Scripta Mater. 108, 1â5 (2015).
Fonari, A. & Sutton, C. Effective Mass Calculator, https://github.com/afonari/emc (2012).
Acknowledgements
The support from National Natural Science Foundation of China [51720105007 and 51806031] and the Fundamental Research Funds for the Central Universities [DUT16RC(3)116], and the computing resources from Supercomputer Center of Dalian University of Technology and ScGrid are greatly acknowledged.
Author information
Authors and Affiliations
Contributions
K.Y. and X.Z. conceived and designed the research. K.Y. carried out the calculations. K.Y. and X.Z. co-wrote the paper with contribution from Z.S. All the authors participated in the discussions and reviewed the manuscript. D.T. supervised the whole project.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
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 license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yuan, K., Sun, Z., Zhang, X. et al. Tailoring phononic, electronic, and thermoelectric properties of orthorhombic GeSe through hydrostatic pressure. Sci Rep 9, 9490 (2019). https://doi.org/10.1038/s41598-019-45949-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-45949-8
This article is cited by
-
Microstructure, Mechanical Property and Thermal Conductivity of AlâSiâFeâMg Alloy Solidified Under High Pressures
International Journal of Metalcasting (2023)
-
Semiclassical electron and phonon transport from first principles: application to layered thermoelectrics
Journal of Computational Electronics (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.