Abstract
Room-temperature superconductivity has been attracting increasing attention in recent years. Recent studies have proved the potential of compressed H-rich materials for achieving room-temperature superconductivity. In this paper, we study the phononic thermal transport in the rare earth yttrium hydrides allotropes under 0, 50, and 300 GPa by using Boltzmann transport equation. We find that the lattice thermal conductivity of yttrium hydrides increases with the pressure among different allotropes, which is attributed to the increase of bond strength and the decrease of phonon-phonon scattering due to structural compression. Yttrium hydrides structure at high pressure of 300 GPa is the superconducting phase, and has high thermal conductivity around 1,360 Wm−1K−1 at room temperature. Comparison of phonon properties with existing high thermal conductivity materials further uncovers the origin for the observed high thermal conductivity. For the zero pressure allotrope, a large number of optical flat bands mix with the low-frequency acoustic phonons, which significantly increases the phonon scattering channel and effectively suppresses the phonon lifetime. As for yttrium hydrides allotropes under 50 and 300 GPa, there are two obvious band gaps in the phonon dispersion relation, and the band gap of the structure at 300 GPa is significantly wider. The occurrence of the band gap effectively inhibits the absorption and emission process of the three-phonon interactions, leading to the decrease of phonon scattering and thus the increase of the phonon lifetime and thermal conductivity at high pressure. Our work reveals the physical mechanism of the thermal transport behaviors in yttrium hydrides structures under different pressures.
Introduction
Because of the fascinating electrical transport ability, the high transition temperature (Tc) superconductivity has attracted lots of research attention (Wigner and Huntington, 1935; ). Previous studies demonstrated that the light element riched metallic solids at high pressure possess promising high-Tc phonon-mediated superconductivity due to the high Debye temperature and strong electron-phonon coupling, such as the H-riched materials (; ; Zhong et al., 2016). H-rich materials have been considered as a potential superconductor at high pressure (; Wang H. et al., 2012) because they can become metallic and superconducting at much lower pressure than hydrogen due to “chemical compression” (). Extensive explorations of the superconducting high pressure phase diagrams of rare earth (Sc, Y, La, Ce, Pr, etc.) hydrides are revealed (; ; ), concentrating on H-rich materials by performing swarm intelligence algorithm based on Crystal structure Analysis by Particle Swarm Optimization structure searches (Wang et al., 2010; Wang Y. et al., 2012). For example, reported some high pressure phases of yttrium hydrides (YH3) that are more likely to be superconducting materials.
On the other hand, the thermal transport property is crucial for the performance and reliability of various devices, which can have promising applications in both thermal management and thermoelectric devices (; ; ; Xie et al., 2018; Zhang and Chen, 2018; Zhang et al., 2018; ; Zhu et al., 2019; Xu et al., 2020; ; ). Especially the materials with high thermal conductivity and excellent mechanical properties are of great importance in solving the heat dissipation problem of highly integrated electronic devices (). Meanwhile, the thermal transport properties of superconductors have also been investigated (; ; ). For example, found the thermal conductivity of the yttrium hydrides can reach about 80 W/mK, which is larger than that of metal and depends significantly on temperatures. Moreover, previous studies also demonstrated that lattice thermal conduction can be largely tuned by external pressure or strains (; ; ). Remarkably, due to the change of the valence electronic state of Y atom under high pressure, there is more possibility of bonding between Y atom and H atom (). Among these different yttrium hydrides, YH3 allotropes provides a broad platform for the study of the structural behavior and pressure tunable thermal properties (; ; ). For example, YH3 is insulating in nature at ambient pressure (0 GPa) (), while the insulator to superconductor transition in YH3 occurs at pressure around 20 GPa (; ). More interestingly, with increasing external pressure, the system becomes normal metallic (i.e., not superconductivity) around 35 GPa, eventually, the superconducting phase re-emerges in YH3 when pressure is greater than 50 GPa (). As lots of new allotropes of YH3 at different pressures have been proposed recently, it would be interesting to study the thermal conductivity of these new superconductors and also the effect of pressure on their thermal transport properties. Here, we focus on the effect of pressure on thermal conductivity of yttrium hydrides allotropes, so we choose the stable structures for YH3 under different pressures (0, 50, and 300 GPa) as research objects.
In this work, based on the Boltzmann transport equation and ab initio calculations, we study the phononic thermal transport properties of yttrium hydrides allotropes at different pressures. We limit our discussion on the lattice thermal conductivity only, as the electron contribution and electron-phonon coupling are neglected. The simulation results indicate that the lattice thermal conductivity of yttrium hydrides is strongly dependent on the pressure. An unexpected high thermal conductivity, around 1,360 Wm−1K−1 at room temperature, is observed in the structure under high pressure (300 GPa). The underlying mechanisms are studied from the structural compressibility, group velocity phonon dispersion, and phonon lifetime. Our study uncovers thermal transport behaviors in yttrium hydrides, and also provides the perspective to achieve high lattice thermal conductivity via pressure.
Computational Method
The computational studies in this paper are based on the generalized gradient approximation by adopting the Perdew-Burke-Ernzerhof parametrization as the exchange-correlation for density functional theory. The projector augmented wave potentials are adopted for Y and H in Y-H system with the Vienna ab initio simulation package (). Due to the application of high pressure phase hydrogen, the plane wave basis set is employed with the energy cutoff of 1,000 eV in calculations at 50 and 300 GPa, whereas the energy cutoff adopted in YH3 at 0 GPa is 450 eV. Brillouin zone sampling is done with a 16 × 16 × 16 Monkhorst-Pack k-points mesh () for the YH3 structures. In order to ensure the accuracy of the results and decrease the calculation scale, the self-consistent convergence of the total energy is set to 10–7 eV/atom. Combining the Vienna ab initio simulation package and PHONOPY package (), the phonon dispersion relation and density of states are obtained for all structures.
The lattice thermal conductivity () of YH3 was calculated by using ShengBTE code (; ), which implements a fully iterative solution of Peierls-Boltzmann transport equation (PBTE). Based on the solution of PBTE, the lattice thermal conductivity can be described as (; ).where , , , and N are Boltzmann constant, temperature, the volume of the unit cell, the reduced Plank constant and the number of q-points, respectively. Here, is the distribution function, is the phonon frequency of the phonon mode λ defined by the phonon branch, and the sum goes over all the phonon modes that are composed of wave vector and phonon branch. is the phonon group velocity in the -direction, and can be expressed as (; ; )where is phonon lifetime of mode which is obtained by relaxation time approximation (RTA), and is a correction term obtained by fully iterative solution of PBTE that merges the inelastic three-phonon scattering processes ().
The anharmonic three-phonon scattering process is the dominating phonon scattering mechanism in pristine crystals, and the resulting phonon lifetime can be calculated based on the Fermi’s Golden Rule as the inversion of the intrinsic scattering rate ()where and represent the second and third phonon mode scattering corresponding to the first phonon mode . and represent the absorbing () and emitting () three-phonon scattering processes, respectively. The linewidth can be obtained by Fermi's Golden Rule and the third-order force-constant matrix. Apart from the three-phonon scatterings, we have put the phonon-isotopic scattering into consideration in this work, and the scattering rate can be expressed as (; )where is the Pearson deviation coefficient for the mass of isotopes s for atom i (found with relative frequency ) and represents the average mass. The natural concentration of isotope in literature () is adopted in our work. Meanwhile, all scattering rates follow the Matthiessen rule ().
The phase space of all possible three-phonon scattering processes that conserve both energy and quasi-lattice momentum based on phonon dispersions and Fermi's Golden Rule is determined. The allowed three-phonon scattering processes must conform towhere is the reciprocal lattice vector. represents the normal (N) processes, while denotes the umklapp (U) processes. The parameter denote the phase space corresponding to absorption and emission of the phonon mode , which are defined as the sum of frequency-dependent factors in the three-phonon transition probabilities expressions, as follows ()
Results and Discussion
Phonon Dispersion and Group Velocity
The unit cells of YH3 under different pressures are displayed in Figure 1. As reported by , these crystals have strong thermal stability. The three types of YH3 are all three-dimensional structures with different space groups and number of atoms in their unit cells (). The optimized lattice constants and atomic position coordinates of YH3 in our calculations are obtained from the work by , as listed in Table 1. These different phases of YH3 provide a valuable comparison for revealing the deep relationship among structures, electrical properties and thermal conductivities (; ; ). Due to the same stoichiometry of YH3 crystals, we denote these structures with different pressures for convenience in the following sections.
FIGURE 1
TABLE 1
| Space group | Lattice Parameters(Å,°) | Atom | X | Y | Z | |
|---|---|---|---|---|---|---|
| YH3 (0 GPa) | P3c1 | a = b = 6.358 | Y(6f) | 0.66312 | 0.00000 | 0.25000 |
| c = 6.606 | H(12g) | 0.34872 | 0.02500 | 0.09332 | ||
| α = β = 90.00° | H(2a) | −0.25000 | 0.37872 | 0.35212 | ||
| γ = 120.00° | H(4d) | 0.33333 | 0.66667 | 0.18101 | ||
| YH3 (50 GPa) | I4/mmm | a = b = 3.29560 | Y(2a) | 0.00000 | 0.00000 | 0.00000 |
| c = 4.72730 | H(4d) | 0.00000 | 0.50000 | 0.25000 | ||
| α = β = γ = 90.00° | H(2b) | 0.00000 | 0.50000 | 0.00000 | ||
| YH3 (300 GPa) | Cmcm | a = 2.55080 | Y(4c) | 0.00000 | −0.11989 | 0.75000 |
| b = 5.97010 | H(4c) | 0.00000 | −0.82205 | 0.75000 | ||
| c = 4.03880 | H(8f) | 0.00000 | −0.60333 | 0.41492 | ||
| α = β = γ = 90.00° | — | — | — | — |
Structural parameters of YH3 for 0, 50 and 300 GPa, respectively, adapted from Ref. ().
Figure 2 shows the phonon dispersions for YH3 structures calculated based on the harmonic force constant matrix under different pressures. For the zero pressure structure, the phonon band below 6 THz is notably flattened, which provides more channels for phonon-phonon scattering. In contrast, the phonon dispersion of YH3 at 50 and 300 GPa have two obvious band gaps. The first band gap is between acoustic and medium optical branches, and the second band gap is between medium optical and highest optical branches. With the increase of pressure, the first band gap becomes wider. For instance, the first band gap of YH3 at 300 GPa is 12.5 THz, which is 42% wider than that of YH3 at 50 GPa. The appearance of phonon band gap can effectively inhibit the absorption and emission three-phonons processes and further affect the thermal conductivity (), which implies the high thermal conductivity of YH3 under 300 GPa.
FIGURE 2
In addition, YH3 structure becomes even harder with the increase of pressure, which is manifested by the increases of frequency for both acoustic and optical phonon branches. By comparing their lattice constants between these three YH3 structures in Table 1, we can find that the unit cell volume is compressed by more than 70% when the pressure increases from 0 to 300 GPa. As the unit cell volume of YH3 is compressed, all bond lengths in YH3 decrease significantly, and the bond energy increases sharply, especially the H-H bonds with weak covalent interactions. The enhancement of band energy will increase the harmonic force constant, which will bring a vast improvement in phonon frequency and group velocity of YH3. Consequently, as shown in Figure 3, this dependence leads to the increase of phonon group velocities with pressure, particularly for the substantial enhancement in group velocity for low frequency phonons (below 15 THz). The phonon group velocity of YH3 at 300 GPa is the largest, which is significantly higher than the case at 0 GPa.
FIGURE 3
Thermal Conductivity
Figure 4 shows the calculated temperature dependent thermal conductivities of YH3 structures, which decreases significantly as temperature increases for all pressures. In our calculation, anharmonic phonon-phonon scattering acts as a major scattering mechanism in pristine crystals without defect or boundary scattering. At high temperature, more phonons are excited due to the enhanced phonon population with temperature. At the same time, the anharmonic effect is enhanced, which greatly increases the scattering rate and limits the phonon mean free path. Moreover, we find that thermal conductivities of YH3 under different pressures follow the ĸ∼T−1 law quite well at high temperatures, indicating the dominance of U-process in the three-phonon scattering process (). More significantly, the enhancement of thermal conductivity is observed with increasing pressure. The room temperature thermal conductivity of YH3 reaches 1,360 Wm−1K−1 at 300 GPa. Such high thermal conductivity is comparable to the superior thermal conductive materials, such as graphene (∼3,000 Wm−1K−1) (), boron arsenide (∼1,400 Wm−1K−1) (; ; ; ), boron phosphide (580 Wm−1K−1) (Zheng et al., 2018), BC2N (1,200 Wm−1K−1) (). In addition, the room temperature thermal conductivities of YH3 at 0 and 50 GPa are around 10 and 70 Wm−1K−1, respectively. Previous studies reveal that thermal conductivity in solids is positively associated with external pressure, and the relationship is almost linear when the lattice structure is not changed (). In this work, the observed dependence of thermal conductivity in YH3 on pressure is highly nonlinear, suggesting that the influences of phase transformation and structure variation are quite significant to the thermal conductivity of YH3. For the whole temperature range, the thermal conductivity of YH3 at 300 GPa is at least one order of magnitude higher than that at low pressure. However, it can be seen from Figure 3 that the phonon group velocities for three structures, although different, are still on the same order of magnitude. Obviously, the phonon group velocity is not the dominant factor responsible for the observed significantly different thermal conductivity.
FIGURE 4
To further explore the different thermal conductivity among three structures, the normalized accumulative thermal conductivities at 300 K vs. phonon frequency are presented in Figure 5. It is noteworthy that accumulative thermal conductivities in all YH3 structures increase rapidly with frequency in the low-frequency acoustic region (<15 THz), indicating the dominant contributions from acoustic phonon branches. With the increased frequency of phonon, the accumulative thermal conductivity of YH3 at 0 GPa converges after ∼4.8 THz, while the saturation frequency is ∼9.8 THz and ∼12.6 THz (vertical lines in Figure 5) for the 50 GPa structure and 300 GPa structure, respectively. This increase of saturation frequency with pressure among different structures is consistent with the lifting of phonon frequency upon pressure in the dispersion relation observed in Figure 2.
FIGURE 5
Phonon Lifetime and Phase Space
Considering that the harmonic property in terms of phonon group velocity is on the same order of magnitude among three YH3 structures, the anharmonic phonon-phonon scatterings should have a significant impact on thermal transport in YH3. In this regard, we further compare in Figure 6A the phonon lifetime of three YH3 structures at 300 K. The phonon lifetime of YH3 at 300 GPa is significantly higher than that in other two structures for the whole frequency range, while the phonon lifetime in the structure at 0 GPa is the smallest. Particularly for the low-frequency phonons below 15 THz, the phonon lifetime in the structure at 300 GPa is at least one order of magnitude higher than that in other two structures, which is consistent with the observed notable difference in thermal conductivity. This result further indicates that the large difference in thermal conductivities of structures under different pressures is mainly caused by anharmonic phonon-phonon scattering. Moreover, we found that the phonon relaxation time of YH3 is on the same order of magnitude as that of the high thermal conductivity materials such as graphene (). Therefore, the large phonon group velocity and relaxation time, especially in the low frequency region (<15 THz), eventually lead to the high thermal conductivity of YH3 at 300 GPa.
FIGURE 6
The suppression of phonon lifetime should be originated from the detailed phonon-phonon scattering channels, which can be evaluated by the phase space parameter W in Eq. 6. The phase space parameter W essentially quantifies the existing scattering channels for various phonons, with a large value indicating frequent phonon-phonon scattering process thus a small phonon lifetime ().
As shown in Figure 6B, the phase space of the structure at 0 GPa is the largest for the whole frequency range, followed by the structure at 50, and 300 GPa. This observation also agrees well with the obtained phonon lifetime in Figure 6A. As we mentioned before, there are two obvious band gaps in the phonon dispersion for the YH3 structures at 50 and 300 GPa, which is absent in the structure at 0 GPa. Moreover, the band gap for the structure at 300 GPa is notably wider than that at 50 GPa. The existence of band gap makes the phonon absorption process () and emission process () more difficult (). On the one hand, when the pressure varies from 0 to 300 GPa, YH3 undergoes a phase transition from a hexagonal to a cubic crystal system. Table 1 reveals that the unit cell volume is compressed by more than 70% when the pressure increases from 0 to 300 GPa. As the unit cell volume of YH3 is compressed, all bond lengths in YH3 decrease significantly, and the bond energy increases sharply, especially the H-H bonds with weak covalent interactions. For crystalline materials, the enhancement of bond energy will increase the force constant matrix of lattice vibration, which makes the phonon frequency and group velocity increase rapidly (; ), especially for the acoustic phonon branch. Moreover, the increase of pressure can weaken the phonon-phonon scattering and increase the phonon relaxation time, thus enhancing the material’s thermal conductivity (). As the number of three phonon-phonon scattering channels is greatly reduced at high pressure, the phonon lifetime of low-frequency acoustic phonons and further thermal conductivity is significantly enhanced. As the band gap is smaller in the structure at 50 GPa, the enhancement of phonon lifetime and thermal conductivity is suppressed. Moreover, because of the mixed acoustic and optical branches for the structure at 0 GPa, the pronounced phonon-phonon scattering processes are observed, manifested by the large phase space parameter. Thus, the suppressed phonon lifetime and thermal conductivity are observed in the structure at 0 GPa.
Conclusion
We have studied the lattice thermal conductivities of YH3 structures under different pressures by using Boltzmann transport equation and first-principle calculations. It is found that thermal conductivity of YH3 increases nonlinearly with pressure. Strikingly, the YH3 structure at 300 GPa, which is the superconducting phase, has remarkably high thermal conductivity around 1,360 W m−1K−1 at room temperature. For the YH3 structure at 0 GPa, a large number of flat bands appear at the low-frequency region in the phonon dispersion relation, and are mixed with the acoustic branches, which causes the strong phonon-phonon scattering and thus effectively suppresses the phonon lifetime and thermal conductivity. With increasing pressure, the unit cell volume of YH3 is compressed by more than 70% and YH3 structure changes from a hexagonal to a cubic crystal system. The increase of group velocity due to the enhancement of bond energy and the weakened phonon-phonon scattering are both responsible for the ultrahigh thermal conductivity of YH3 at 300 GPa. Our study uncovers the phononic thermal transport properties of yttrium hydrides allotropes, and also provides a new opportunity to achieve high thermal dissipation ability with external pressure.
Funding
This project is supported in part by the grants from the National Natural Science Foundation of China (Grant Nos. 11890703 and 11775158), Science and Technology Commission of Shanghai Municipality (Grant Nos. 19ZR1478600 and 18JC1410900), the Fundamental Research Funds for the Central Universities (Grant No. 22120200069), and Open Fund of Hunan Provincial Key Laboratory of Advanced Materials for New Energy Storage and Conversion (Grant No. 2018TP1037_201901).
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
WR performed the calculations. JC supervised the project. All authors discussed the results and wrote the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlvaG.LinY.FangG. (2018). An overview of thermal energy storage systems. Energy144, 341–378. 10.1016/j.energy.2017.12.037
2
AshcroftN. W. (1968). Metallic hydrogen: a high-temperature superconductor?. Phys. Rev. Lett.21, 1748–1749. 10.1103/physrevlett.21.1748
3
AshcroftN. W. (2004). Hydrogen dominant metallic alloys: high temperature superconductors?. Phys. Rev. Lett.92, 187002. 10.1103/physrevlett.92.187002
4
BalandinA. A. (2011). Thermal properties of graphene and nanostructured carbon materials. Nature Mater10, 569–581. 10.1038/nmat3064
5
BeldjoudiK.SahraouiF. A.BouhemadouA. (2019). Pressure effect on the structural, elastic and thermodynamic properties of the BeP2N4 compound: first-principles investigation. Comput. Condens. Matter21, e00408. 10.1016/j.cocom.2019.e00408
6
ChenJ.WaltherJ. H.KoumoutsakosP. (2014). Strain engineering of kapitza resistance in few-layer graphene. Nano Lett.14, 819–825. 10.1021/nl404182k
7
ErreaI.CalandraM.PickardC. J.NelsonJ. R.NeedsR. J.LiY.et al (2016). Quantum hydrogen-bond symmetrization in the superconducting hydrogen sulfide system. Nature532, 81–84. 10.1038/nature17175
8
FengT.LindsayL.RuanX. (2017). Four-phonon scattering significantly reduces intrinsic thermal conductivity of solids. Phys. Rev. B96, 161201. 10.1103/physrevb.96.161201
9
FengT.QiuB.RuanX. (2015). Coupling between phonon-phonon and phonon-impurity scattering: a critical revisit of the spectral Matthiessen's rule, Phys. Rev. B92, 235206. 10.1103/physrevb.92.235206
10
GhoshS.CalizoI.TeweldebrhanD.PokatilovE. P.NikaD. L.BalandinA. A.et al (2008). Extremely high thermal conductivity of graphene: prospects for thermal management applications in nanoelectronic circuits. Appl. Phys. Lett.92, 151911. 10.1063/1.2907977
11
GrishakovK. S.DegtyarenkoN. N.MazurE. A. (2019). Electron, phonon, and superconducting properties of yttrium and sulfur hydrides under high pressures. J. Exp. Theor. Phys.128, 105–114. 10.1134/s1063776119010072
12
HeilC.di CataldoS.BacheletG. B.BoeriL. (2019). Superconductivity in sodalite-like yttrium hydride clathrates, Phys. Rev. B99, 220502. 10.1103/physrevb.99.220502
13
HuS.ZhangZ.JiangP.ChenJ.VolzS.NomuraM.et al (2018). Randomness-induced phonon localization in graphene heat conduction. J. Phys. Chem. Lett.9, 3959–3968. 10.1021/acs.jpclett.8b01653
14
HummelP.LechnerA. M.HerrmannK.BiehlP.RösselC.WiedenhöftL.et al (2020). Thermal transport in ampholytic polymers: the role of hydrogen bonding and water uptake. Macromolecules53, 5528–5537. 10.1021/acs.macromol.0c00596
15
ItoM.MatsunagaJ.SetoyamaD.MutaH.KurosakiK.UnoM.et al (2005). Thermal properties of yttrium hydride. J. Nucl. Mater.344, 295–297. 10.1016/j.jnucmat.2005.04.058
16
JarosikM. W.SzczȩśniakR.WronaI. A.KostrzewaM. (2017). Non-BCS superconducting state in yttrium hydride at a record low value of the external pressure. Solid State Commun.250, 5–8. 10.1016/j.ssc.2016.11.002
17
KangJ. S.LiM.WuH.NguyenH.HuY. (2018). Experimental observation of high thermal conductivity in boron arsenide. Science361, 575–578. 10.1126/science.aat5522
18
KimD. Y.ScheicherR. H.AhujaR. (2009). Predicted high-temperature superconducting state in the hydrogen-dense transition-metal hydride YH3 at 40 K and 17.7 GPa. Phys. Rev. Lett.103, 077002. 10.1103/physrevlett.103.077002https://link.aps.org/doi/10.1103/PhysRevLett.103.077002
19
KimD. Y.ScheicherR. H.MaoH.-k.KangT. W.AhujaR. (2010). General trend for pressurized superconducting hydrogen-dense materials. Proc. Natl. Acad. Sci. U. S. A.107, 2793–2796. 10.1073/pnas.0914462107
20
KresseG.FurthmüllerJ. (1996). Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set. Phys. Rev. B54, 11169–11186. 10.1103/physrevb.54.11169
21
KunduA.MingoN.BroidoD. A.StewartD. A. (2011). Role of light and heavy embedded nanoparticles on the thermal conductivity of SiGe alloys. Phys. Rev. B84, 125426. 10.1103/physrevb.84.125426
22
LiR.ChengX.XieQ.SunY.LiD.LiY.et al (2015). Topological metal of NaBi with ultralow lattice thermal conductivity and electron-phonon superconductivity. Sci. Rep.5, 8446. 10.1038/srep08446
23
LiS.ZhengQ.LvY.LiuX.WangX.HuangP. Y.et al (2018). High thermal conductivity in cubic boron arsenide crystals. Science361, 579–581. 10.1126/science.aat8982
24
LiW.CarreteJ.KatchoN. A.MingoN. (2014). ShengBTE: a solver of the Boltzmann transport equation for phonons. Comput. Phys. Commun.185, 1747–1758. 10.1016/j.cpc.2014.02.015
25
LiW.LindsayL.BroidoD. A.StewartD. A.MingoN. (2012). Thermal conductivity of bulk and nanowire Mg2SixSn1-x alloys from first principles. Phys. Rev. B86, 174307. 10.1103/physrevb.86.174307
26
LiW.MingoN. (2014). Thermal conductivity of fully filled skutterudites: role of the filler. Phys. Rev. B89, 184304. 10.1103/physrevb.89.184304
27
LiW.MingoN. (2015). Ultralow lattice thermal conductivity of the fully filled skutterudite YbFe4Sb12 due to the flat avoided-crossing filler modes. Phys. Rev. B91, 144304. 10.1103/physrevb.91.144304
28
LiX.PengF. (2017). Superconductivity of pressure-stabilized vanadium hydrides. Inorg. Chem.56, 13759–13765. 10.1021/acs.inorgchem.7b01686
29
LindsayL.BroidoD. A. (2008). Three-phonon phase space and lattice thermal conductivity in semiconductors. J. Phys. Condens. Matter20, 165209. 10.1088/0953-8984/20/16/165209
30
LindsayL.BroidoD. A.ReineckeT. L. (2013). First-Principles determination of ultrahigh thermal conductivity of boron arsenide: a competitor for diamond?. Phys. Rev. Lett.111, 025901. 10.1103/physrevlett.111.025901
31
LiuH.NaumovI. I.HoffmannR.AshcroftN. W.HemleyR. J. (2017). Potential high-Tc superconducting lanthanum and yttrium hydrides at high pressure. Proc. Natl. Acad. Sci. U. S. A.114, 6990–6995. 10.1073/pnas.1704505114
32
LiuH.QinG.LinY.HuM. (2016). Disparate strain dependent thermal conductivity of two-dimensional penta-structures. Nano Lett.16, 3831–3842. 10.1021/acs.nanolett.6b01311
33
LiuL.-L.SunH.-J.WangC. Z.LuW.-C. (2017). High-pressure structures of yttrium hydrides. J. Phys. Condens. Matter29, 325401. 10.1088/1361-648x/aa787d
34
MannS.RaniP.KumarR.DubeyG. S.JindalV. K. (2016). Thermodynamic properties of pure and doped (B, N) graphene. RSC Adv.6, 12158–12168. 10.1039/c5ra25239c
35
MichaelB.MichaelE. W. (2011). Isotopic compositions of the elements 2009 (IUPAC technical report). Pure Appl. Chem.83, 397–410. 10.1351/PAC-REP-10-06-02
36
MonkhorstH. J.PackJ. D. (1976). Special points for Brillouin-zone integrationsPhys. Rev. B13, 5188–5192. 10.1103/physrevb.13.5188
37
OminiM.SparavignaA. (1995). An iterative approach to the phonon Boltzmann equation in the theory of thermal conductivity. Phys. B Condens. Matter212, 101–112. 10.1016/0921-4526(95)00016-3
38
OuyangT.HuM. (2015). Competing mechanism driving diverse pressure dependence of thermal conductivity of XTe (X = Hg, Cd, and Zn). Phys. Rev. B92, 235204. 10.1103/physrevb.92.235204
39
OuyangY.ZhangZ.LiD.ChenJ.ZhangG. (2019). Emerging theory, materials, and screening methods: new opportunities for promoting thermoelectric performance. Ann. Phys.531, 1800437. 10.1002/andp.201800437
40
PengF.SunY.PickardC. J.NeedsR. J.WuQ.MaY. (2017). Hydrogen clathrate structures in rare earth hydrides at high pressures: possible route to room-temperature superconductivity. Phys. Rev. Lett.119, 107001. 10.1103/physrevlett.119.107001
41
ShafiqueA.ShinY.-H. (2019). Ultrahigh and anisotropic thermal transport in the hybridized monolayer (BC2N) of boron nitride and graphene: a first-principles study. Phys. Chem. Chem. Phys.21, 17306–17313. 10.1039/c9cp02068c
42
ShiL.ChenJ.ZhangG.LiB. (2012). Thermoelectric figure of merit in Ga-doped [0001] Z00nO nanowires. Phys. Lett.376, 978–981. 10.1016/j.physleta.2011.12.040
43
SorokaO.SturmJ. M.van de KruijsR. W. E.LeeC. J.BijkerkF. (2018). Control of YH3 formation and stability via hydrogen surface adsorption and desorption, Appl. Surf. Sci.455, 70–74. 10.1016/j.apsusc.2018.05.134
44
TamuraS.-i. (1983). Isotope scattering of dispersive phonons in Ge. Phys. Rev. B27, 858–866. 10.1103/physrevb.27.858
45
TanakaK.TseJ. S.LiuH. (2017). Electron-phonon coupling mechanisms for hydrogen-rich metals at high pressure. Phys. Rev. B96, 100502. 10.1103/physrevb.96.100502
46
TianF.SongB.ChenX.RavichandranN. K.LvY.ChenK.et al (2018). Unusual high thermal conductivity in boron arsenide bulk crystals. Science361, 582–585. 10.1126/science.aat7932
47
TogoA.ObaF.TanakaI. (2008). First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-typeSiO2 at high pressures. Phys. Rev. B78, 134106. 10.1103/physrevb.78.134106
48
WangY.LvJ.ZhuL.MaY. (2010). Crystal structure prediction via particle-swarm optimization. Phys. Rev. B82, 094116. 10.1103/physrevb.82.094116
49
WangH.TseJ. S.TanakaK.IitakaT.MaY. (2012). Superconductive sodalite-like clathrate calcium hydride at high pressures. Proc. Natl. Acad. Sci. U. S. A.109, 6463–6466. 10.1073/pnas.1118168109
50
WangY.LvJ.ZhuL.MaY. (2012). CALYPSO: a method for crystal structure prediction. Comput. Phys. Commun.183, 2063–2070. 10.1016/j.cpc.2012.05.008
51
WignerE.HuntingtonH. B. (1935). On the possibility of a metallic modification of hydrogen. J. Chem. Phys.3, 764–770. 10.1063/1.1749590
52
XieG.JuZ.ZhouK.WeiX.GuoZ.CaiY.et al (2018). Ultra-low thermal conductivity of two-dimensional phononic crystals in the incoherent regime. Npj Comput. Mater.4, 21. 10.1038/s41524-018-0076-9
53
XuX.ZhouJ.ChenJ. (2020). Thermal transport in conductive polymer-based materials. Adv. Funct. Mater.30, 1904704. 10.1002/adfm.201904704
54
ZhangZ.ChenJ. (2018). Thermal conductivity of nanowires. Chin. Phys. B27, 035101. 10.1088/1674-1056/27/3/035101
55
ZhangZ.ChenJ.LiB. (2017a). Negative Gaussian curvature induces significant suppression of thermal conduction in carbon crystals. Nanoscale9, 14208–14214. 10.1039/c7nr04944g
56
ZhangZ.HuS.ChenJ.LiB. (2017b). Hexagonal boron nitride: a promising substrate for graphene with high heat dissipation. Nanotechnology28, 225704. 10.1088/1361-6528/aa6e49
57
ZhangZ.HuS.NakayamaT.ChenJ.LiB. (2018). Reducing lattice thermal conductivity in schwarzites via engineering the hybridized phonon modes. Carbon139, 289–298. 10.1016/j.carbon.2018.06.057
58
ZhangZ.HuS.XiQ.NakayamaT.VolzS.ChenJ.et al (2020a). Tunable phonon nanocapacitor built by carbon schwarzite based host-guest system. Phys. Rev. B101, 081402(R). 10.1103/physrevb.101.081402
59
ZhangZ.OuyangY.ChengY.ChenJ.LiN.ZhangG. (2020b). Size-dependent phononic thermal transport in low-dimensional nanomaterials. Phys. Rep.860, 1–26. 10.1016/j.physrep.2020.03.001
60
ZhengQ.LiS.LiC.LvY.LiuX.HuangP. Y.et al (2018). High thermal conductivity in isotopically enriched cubic boron phosphide. Adv. Funct. Mater.28, 1805116. 10.1002/adfm.201805116
61
ZhongX.WangH.ZhangJ.LiuH.ZhangS.SongH.-F.et al (2016). Tellurium hydrides at high pressures: high-temperature superconductors, Phys. Rev. Lett.116, 057002. 10.1103/physrevlett.116.057002
62
ZhuX.-L.LiuP.-F.ZhangJ.ZhangP.ZhouW.-X.XieG.et al (2019). Monolayer SnP3: an excellent p-type thermoelectric material. Nanoscale11, 19923–19932. 10.1039/c9nr04726c
Summary
Keywords
Peierls-Boltzmann transport equation, lattice thermal conductivity, superconductor, phonon lifetime, first-principle calculations
Citation
Ren W, Zhang Z, Chen C, Ouyang Y, Li N and Chen J (2020) Phononic Thermal Transport in Yttrium Hydrides Allotropes. Front. Mater. 7:569090. doi: 10.3389/fmats.2020.569090
Received
03 June 2020
Accepted
27 October 2020
Published
17 November 2020
Volume
7 - 2020
Edited by
Weifeng Li, Shandong University, China
Reviewed by
Ibrahim Sayed Hussein, King Khalid University, Saudi Arabia
Yinchang Zhao, Yantai University, China
Yandong Ma, Shandong University, China
Updates
Copyright
© 2020 Ren, Zhang, Chen, Ouyang, Li and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Chen, jie@tongji.edu.cn
This article was submitted to Thin Solid Films, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.