Shuang Zheng
Shoutao Zhang
Yu Sun1
Guochun Yang- 1Centre for Advanced Optoelectronic Functional Materials Research and Key Laboratory for UV Light-Emitting Materials and Technology of Ministry of Education, Northeast Normal University, Changchun, China
- 2Departmento de Física de la Materia Condensada, Universidad del País Vasco, UPV/EHU, Bilbao, Spain
- 3Donostia International Physics Center, Donostia, Spain
- 4Centro de Física de Materiales, Centro Mixto CSIC-UPV/EHU, Donostia, Spain
Unveiling the relation between crystal structure and superconductivity of hydrides becomes a fascinating research area in chemistry and condensed-mater physics. Although much efforts have been made to study chemical reaction between tungsten and hydrogen, the crystal structures, superconductivity, and phase diagram of tungsten hydrides under high pressure have not been fully explored and built thus far. In this work, we carried out extensive structural search on W-H binary compounds through first-principles swarm-intelligence structural search calculations. Besides reproducing the known W-H compounds, a new stoichiometry WH5 with P6mm symmetry becomes stable above 230.2 GPa. Intriguingly, P6mm WH5 shows a critical temperature (Tc) value of 60.8 K, which is much higher than 31.6 K in WH6. This finding is different from the knowledge of compounds with higher H content exhibiting higher Tc, which might be attributed to the appearance of unique H network and tetrahedron H units in WH5. Electronic property and superconductivity of the other tungsten hydrides are also investigated. The built pressure-composition phase diagram provides some useful information for experimental synthesis.
Introduction
An important objective in condensed-matter physics and chemistry is to prepare hydrogen-rich (H-rich) compounds because of their broad applications and scientific significance [1–5]. In more detail, H-rich compounds are the excellent hydrogen energy storage materials [6–8] and high-temperature superconductors [9–11]. On the other hand, H-rich compounds, especially at high pressures, contain various kinds of chemical bonds (e.g., covalent, ionic, and van der Waals bonds) [12] and diverse hydrogen structure units [e.g., unpaired H [13], H2 [14], H3 units [15], and 3D hydrogen network resembling metallic atomic hydrogen [16, 17]], providing opportunity for us to fully understand physical and chemical properties of hydrogen.
Wigner et al. predicted that metalized solid hydrogen is the potential room-temperature superconductor due to its high phonon frequencies and strong electron-phonon coupling [18]. However, hydrogen metallization is rather difficult such as extreme high pressure [19, 20]. Fortunately, Ashcroft proposed that metallization of hydrogen-rich compounds needs relatively lower pressures [10], which also show high-temperature superconductivity [1, 2, 13, 21, 22]. Among the various H-rich compounds, metal hydrides as the potential superconductors have drawn much attentions because of large density of states at the Fermi level, high phonon frequencies, and strong electron-phonon coupling [23, 24]. Great progress has been made in metal hydrides [11, 25–28]. Interestingly, some of metal hydrides under high pressure show large Tc values [2, 29, 30]. For example, the Tc for LiH6 reaches 82 K at 300 GPa [31]. The C2/c phase of KH6 is a superconductor with an estimated Tc of 58.66–69.84 K at 166 GPa [5]. Recent theoretical calculation revealed that the potential high-temperature superconductor of CaH6 has a Tc of 220–235 K at 150 GPa [2]. For Mg-H system, MgH6 shows a Tc value of 263 K at 300 GPa [29]. Moreover, a theoretical work indicated that Tc value of 34.1–40.5 K for R-3 InH3 at 200 GPa [15]. The P6/mmm SnH4 was predicted to be a superconductor with a high Tc (73 K at 120 GPa) [21].
Among the various hydrides, transition metal hydrides have drawn a growing interest due to their lower synthesis pressure with respect to the other hydrides, multiple stoichiometries, catalytic transformations, and superconductivity [32–34]. These unique properties can be attributed to the diverse d electron configurations of transition metals. More interestingly, transition metal hydrides with higher hydrogen content become stabilized with increasing pressure [35, 36], consisting of diverse hydrogen structure units and exhibiting intriguing properties [6, 11–16]. For instance, FeHx (x = 1–4) compounds have been reported [37–40]. Very recently, Pepin et al. reported synthesis of FeH5 containing non-bonded hydrogen atoms above 130 GPa, which opens up a strategy to access metallic atomic hydrogen [41]. It is worth noting that finding metallic atomic hydrogen at a low pressure is the long-pursuit targets. In addition, FeH5 shows a superconducting transition temperature of 51 K [42]. Inspired by these findings, Kvashnin et al. [43] and Zhang et al. [44] independently studied phase diagram and electronic properties of Fe-H binary compounds with diverse stoichiometries at high pressures. A new H-rich compound, FeH6, was identified. Notably, Cmmm FeH6 is metallic [43], whereas C2/c FeH6 is non-metallic [44], which becomes rare example of pressure-induced non-metallicity among the transition-metal hydrides [45]. On the other hand, parts of transition metal hydrides especially at high pressures are metallization and exhibit superconductivity [9, 11, 13, 22, 33, 35]. In more detail, Li et al. revealed that the YH4 and YH6 are potential superconductors with the estimated Tc values of 84–95 K and 251–264 K at 120 GPa [13]. It is found that TaH6 with an estimated Tc of 124.2–135.8 K under pressure is a candidate as a high-temperature superconductor [35]. As for vanadium hydrides, the high-pressure structures of VH5 and VH8 have been explored, their Tc values were predicted to be 18.5 and 71.4 K at 200 GPa due to the strong electron-phonon coupling (EPC) [36]. Consequently, these studies motivated us to search high-temperature superconductors of transition metal hydrides through compression.
Tungsten (W) is one of transition metal elements. At atmospheric pressure, molecular tungsten hydrides such as WH2, WH4, and molecular complex WH4(H2)4 were identified using spectroscopy at low temperatures and in neon matrix [46, 47]. However, they are unstable with respect to hydrogen and bulk W at ambient conditions. In other words, there is no stable solid tungsten hydride at ambient conditions. Once, it was thought that W hardly reacts with H2 under pressure [48], and was used as a sealing material in high-pressure experiment [19]. However, finding near-stoichiometric WH above ~30 GPa arouses much attention to studying tungsten hydrides at high pressures [49]. Specifically, variable composition structural prediction was performed to investigate structures and phase stability of WHn (n = 1–6, 8, 10, and 12) up to 150 GPa [46, 50]. Three W-H compounds (i.e., WH, WH2, and WH6) are found to be stable [46, 51]. Subsequently, Gregoryanz et al. studied the formation, decomposition, and mechanical properties of tungsten hydrides in a diamond anvil cell up to 115 GPa [52]. Notably, these studies mainly focus on structures and stabilities of tungsten hydrides below 150 GPa. To the best of our knowledge, their electronic properties and superconductivity have not been explored so far.
Generally speaking, most of metal hydrides exhibit superconductivity, and high pressure is in favor of stabilizing H-rich compounds [e.g., LaH8 [11], BaH6 [53], YH10 [16], and LaH10 [17]]. On the other hand, H-rich compounds show larger Tc value [54]. Therefore, we expect that tungsten hydrides with a higher hydrogen content become stable under high pressure. Moreover, understanding electronic structures and superconductivity of tungsten hydrides is important to explore their potential applications. In this work, we studied the structural stabilities of tungsten hydrides containing various WxHy (x = 1, y = 1–12; x = 2, y = 1) compositions at pressures of 100, 200, and 300 GPa via first-principle structure search calculations. Besides reproducing the reported tungsten hydrides, we find a new stable stoichiometry WH5 consisting of unique H network and tetrahedron H units, exhibiting the highest Tc value among the tungsten hydrides.
Computational Details
The structural prediction of W-H system was performed via particle swarm optimization as implemented in CALYPSO code [55, 56] at the selected pressures of 100, 200, and 300 GPa. The main merit of this code is capable of finding the ground-state structures only depending on the given chemical composition. Its efficiency and validity have been proven by applying different kinds of systems, ranging from elemental solids to binary and ternary compounds [2, 57–62]. Considering these advantages, CALYPSO becomes a leading structure prediction method in the field. Then, the structural relaxations for the considered W-H compounds were performed through the Broyden-Fletcher-Goldfarb-Shanno (BFGS) [63] method as implemented in Cambridge Serial Total Energy Package (CASTEP) program [64]. The Perdew-Burke-Ernzerh (PBE) of generalized gradient approximation (GGA) was selected for the exchange-correlation functional [65]. The electron-ion interaction was described by using the ultrasoft pseudopotential [66] plane-wave method, in which 5s25p65d46s2 and 1s1 were taken as valence electrons for W and H atoms, respectively. The kinetic cutoff energy of 400 eV and the Monkhorst-Pack scheme [67] with a k-point grid of 2π × 0.03 Å−1 were chosen to ensure that the total energy well-converged. The phonon spectra and the electron-phonon coupling were calculated within the framework of the linear-response theory as carried out in the Quantum-ESPRESSO package [68]. The kinetic energy cutoff for wave-function expansion is chosen as 80 Ry. A k-point mesh of 15 × 15 × 14 and a q-point mesh of 5 × 5 × 2 for WH5 were used in the first Brillouin zone. More detailed calculated information can be found in the Supporting Information.
Results and Discussion
To obtain higher Tc value, high hydrogen content in compounds is one of the necessary prerequisites [9, 11, 54]. As a consequence, we carried out the crystal structure search on WxHy (x = 1, y = 1–12; x = 2, y = 1) stoichiometries at the selected pressures of 100, 200, and 300 GPa through first-principles swarm structure searching calculations. First, based on the predicted stable structures, we studied the relative phase stability of each stoichiometric compound by calculating its enthalpy of formation (ΔH = [H(WxHy) – xH(W) – yH(H2)/2]/(x + y)) relative to elemental W and H2 solids at the selected pressures. Here, H(WxHy) is the enthalpy of considered structure, H(W) and H(H2) are the enthalpies of elemental W and H2 solids, respectively. Then, convex hull data for considered W-H compositions at selected pressures are constructed using enthalpy of formation per atom ΔH for each composition to determine its stability (Figure 1A). In general, a structure whose enthalpy of formation sits on the convex hull (i.e., solid line) is thermodynamically and is experimentally synthesizable. At 100 GPa, the known WH, WH2, and WH4 are well-reproduced by our structural search, validating that our adopted method is suitable to the W-H system. At 200 GPa, WH, WH2, and WH6 are stable, whereas WH4 becomes unstable with respect to WH2 and WH6. With further compression, a new stoichiometry WH5, is thermodynamically stable in addition to WH and WH2 at 300 GPa, whereas WH6 decomposes into WH5 and H2. It is noted that WH is the most stable composition in the whole considered pressure range. To provide more information for experimental synthesis, we determine the pressure stability region of W-H compounds, as displayed in Figure 1B. Intriguingly, these stable compounds do not undergo any structural transitions. This is in sharp contrast with other transition metal hydrides, such as VH3 [71] and ScH6 [9].
Figure 1. (A) Phase stabilities of various W-H compounds at 100, 200, and 300 GPa. The elemental W solid phase with symmetry Im-3m [69], and P63/m, C2/c, Cmca-12 elemental H2 solids [70] are adopted to calculate the formation of enthalpy. (B) Pressure-composition phase diagram for the W-H compounds.
WH5 becomes stable above 230.2 GPa and stabilizes into a hexagonal structure (space group P6mm, one formula unit per cell, see Figure 2). This structure contains two inequivalent H's occupying 2b (0.3333, 0.6667, 0.6985) and 3c (0.0000, 0.5000, 0.0825) sites and one equivalent W's sitting at 1a (0.0000, 0.0000, 0.4728) position. Each W atom forms a twelve-fold coordination with H atoms (Figure 2A), having two different W–H distances of 1.70 and 1.71 Å. Mulliken charge analysis [72] shows W atom transfers 0.17 e to each H atom, indicating that W-H bonding is ionic. There are two kinds of H coordination. The H atom at 2b site is coordinated by three W atoms, while H atom at 3c site is shared by two W atoms (Figure 2B). Notably, H atoms form intriguing H network with alternate triangle and hexagon rings along c axis with equal H-H distance of 1.38 Å (Figure 2C), which is called a kagome lattice [73–75]. Moreover, the H atom in the ab plane constitutes vertex-sharing tetrahedron units having H-H distance of 1.28~1.38 Å (Figure 2D). These H-H distances in WH5 is much larger than H2 (0.74 Å) [16] at ambient conditions and monatomic H (0.98 Å) at 500 GPa [76], and comparable to lanthanum and yttrium hydrides (1.02~1.54 Å) [11]. indicating that WH5 might has high Tc value as observed in CaH6 [2], LaH8 [11], and YH10 [13].
Figure 2. Crystal structure of WH5 with P6mm symmetry at 300 GPa. (A) View of W-H polyhedrons in P6mm WH5. (B) Two inequivalent H atoms with three- and two-fold W atoms. (C) H network with alternate triangle and hexagon rings along c axis is a kagome lattice. (D) Vertex-sharing tetrahedron unit consisting of four H atoms in the ab plane.
Dynamic stability of a compound strongly correlates with its practical application. Thus, we calculated phonon spectra of P6mm WH5 as implemented in Quantum-ESPRESSO package (Figure 3A). The absence of any imaginary frequency in the entire Brillouin zone obviously confirms its dynamic stability. Based on the analysis of projected phonon density of states (PHDOS, Figure 3B), there are two separate regions. The high frequency modes (above ~30 THz) mainly originate from the contribution of H atoms, whereas the low-frequency modes (below ~15 THz) associate with W atoms, which is attributed to that H atom has a lighter atomic mass than W. Moreover, P63/mmc WH, Pnma WH2, P4/nmm WH4, and C2/m WH6 are also dynamically stable (Figure S2).
Figure 3. (A) Phonon dispersion curve, (B) Eliashberg spectral function α2F(ω), electron-phonon coupling integration λ(ω), and projected phonon density of states (PHDOS) on W and H atoms, (C) electronic band structure, and (D) PDOS for P6mm WH5 at 300 GPa.
Subsequently, we probe the electronic band structure and corresponding projected density of states (PDOS) of the stable W-H compounds as depicted in Figures 3C,D and Figure S1. As we expected, all the W-H compounds exhibit metallic character due to several bands across their Fermi levels. For these W-H compounds, there is a common character of their electron properties (i.e., the main contribution to PDOS near the Fermi level comes from W 5d orbital). Moreover, there is large overlap between H 1s orbital and W 6s or 6p orbital. The appearance of “flat band-steep band” [3, 77, 78] at Fermi level and large N(Ef) indicate that W-H compounds are the potential high-temperature superconductors [79].
Subsequently, we investigate their superconductivity based on Bardeen-Cooper-Schrieffer (BCS) theory [80]. Among the stable tungsten hydrides, WH5 has the largest electron-phonon coupling parameter, λ, which is comparable to those predicted for H2S [59] and H3Se [78], and slightly smaller that of H3S [77], indicating that there is strong EPC in WH5 (Table 1). The combined projected phonon density of states (PHDOS) with Eliashberg spectral function α2F(ω) (Figure 3B) can be used to determine the origin of λ. For WH5, the low-frequency vibrations from the heavy W atoms (below ~15 THz) contribute 38.7% of the total λ, whereas the high-frequency vibrations (above ~30 THz) resulting from H atoms account for 61.3% of λ. For WH6, the low-frequency vibrations (below ~15 THz) from the W atoms contribute 50% of the total λ (Figure S3), which is in sharp contrast with WH5. The λ (0.74) in WH6 is smaller than 0.93 in WH5, indicating a weak electron-phonon coupling in WH6. On the other hand, compared PHDOS (Figure 3B and Figure S3) between WH5 and WH6, the high-frequency and low-frequency modes of WH5 are clearly separate, unlike the situation of WH6 in which the H-stretching vibrations mix together within the mid-frequency region. These observations can be attributed to the unique H network and the tetrahedron H units in WH5, whereas hydrogen atoms in WH6 directly bond W atoms (Figure S4). Based above results, WH5 with unique H network shows good superconductivity, as will be discussed below. Similar observation has been found in CaH6 [2], LaH8 [11], and YH10 [16].
Table 1. The logarithmic average phonon frequency ωlog, the electronic density of states at the Fermi level N(Ef), the electron-phonon coupling parameter λ, and the superconducting transition temperature Tc of W-H compounds at selected pressures.
The superconducting Tc values of the stable tungsten hydrides are evaluated through the Allen-Dynes modified McMillan equation [81] by using the calculated logarithmic average frequency, ωlog. The typical Coulomb pseudopotential μ* usually adopts the value in the range of 0.1 to 0.13 for many metal hydrides [10]. Thus, μ* with 0.1 and 0.13 was adopted here. The calculated Tc values are shown in Table 1. As expected, WH5 shows the highest Tc value of 51.3–60.8 K, which is attributed to its large λ and high ωlog. The Tc values of W-H compounds decrease when μ* increases from 0.1 to 0.13, similar to V-H [36] and In-H compounds [15]. Meanwhile, we calculated the pressure-dependent Tc value of WH5. As shown in Figure 4, the effect of pressure on Tc values is rather small. This is different from the decreased Tc for CrH3 [82] and ThH10 [83] and the increased Tc in ScH8 [9] with increasing pressure. In more detail, as the pressure increases, ωlog becomes larger, whereas λ becomes smaller, and N(Ef) remains nearly unchanged (Figure 4).
Figure 4. The logarithmic average phonon frequency ωlog, the electronic density of states at the Fermi level N(Ef) (states/Ry/cell), the electron-phonon coupling parameter λ, and the superconducting transition temperature Tc of P6mm WH5 at the selected pressures of 230, 250, 270, and 300 GPa.
Conclusions
To determine phase diagram and superconductivity of W-H binary compounds, we have explored their crystal structures and stabilities with various chemical compositions with the aid of first-principles swarm structural calculations. We identify a hitherto unknown WH5 compound, consisting of unprecedented H network and tetrahedron H units, which can be synthesized from compression of C2/m WH6 and Pnma WH2 above 230.2 GPa. All the stable tungsten hydrides are metallic in their stable pressure region. The P6mm WH5 shows high Tc value up to 60.8 K at 300 GPa, which is much higher than WH6. This work provides an opportunity for full understanding of crystal structures, electron properties, and superconducting mechanism of tungsten hydrides at megabar pressures.
Author Contributions
GY and AB contributed conception and design of the study. ShoZ organized the database. YS, JZ, and JL performed the statistical analysis. ShuZ wrote the first draft of the manuscript. GY and AB wrote sections of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.
Conflict of Interest Statement
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.
Acknowledgments
This work is supported by the Natural Science Foundation of China under Nos. 21573037, 11704062, and 51732003; the Natural Science Foundation of Jilin Province (No. 20150101042JC); the Postdoctoral Science Foundation of China (under Grant No. 2013M541283); and the Fundamental Research Funds for the Central Universities (2412017QD006).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2018.00101/full#supplementary-material
References
1. Drozdov AP, Eremets MI, Troyan IA, Ksenofontov V, Shylin SI. Conventional superconductivity at 203 Kelvin at high pressures in the sulfur hydride system. Nature (2015) 525:73–6. doi: 10.1038/nature14964
2. Wang H, Tse JS, Tanaka K, Iitaka T, Ma Y. Superconductive sodalite-like clathrate calcium hydride at high pressures. Proc Natl Acad Sci USA. (2012) 109:6463–6. doi: 10.1073/pnas.1118168109
3. Li Y, Gao G, Xie Y, Ma Y, Cui T, Zou G. Superconductivity at ~100 K in dense SiH4(H2)2 predicted by first principles. Proc Natl Acad Sci USA. (2010) 107:15708–11. doi: 10.1073/pnas.1007354107
4. Zhong G, Zhang C, Chen X, Li Y, Zhang R, Lin H. Structural, electronic, dynamical, and superconducting properties in dense GeH4(H2)2. J Phys Chem C (2012) 116:5225–34. doi: 10.1021/jp211051r
5. Zhou D, Jin X, Meng X, Bao G, Ma Y, Liu B, et al. Ab initio study revealing a layered structure in hydrogen-rich KH6 under high pressure. Phys Rev B Condens Matter Mater Phys. (2012) 86:14118. doi: 10.1103/PhysRevB.86.014118
6. Mao WL, Mao H-k. Hydrogen storage in molecular compounds. Proc Natl Acad Sci USA. (2004) 101:708–10. doi: 10.1073/pnas.0307449100
7. Wolverton C, Siegel DJ, Akbarzadeh AR, Ozoli V. Discovery of novel hydrogen storage materials: an atomic scale computational approach. J Phys Condens Matter (2008) 20:64228. doi: 10.1088/0953-8984/20/6/064228
8. Graetz J, Reilly JJ, Yartys VA, Maehlen JP, Bulychev BM, Antonov VE, et al. Aluminum hydride as a hydrogen and energy storage material: past, present and future. J Alloys Comp. (2011) 509S:S517–28. doi: 10.1016/j.jallcom.2010.11.115
9. Qian S, Sheng X, Yan X, Chen Y, Song B. Theoretical study of stability and superconductivity of ScHn (n = 4–8) at high pressure. Phys Rev B (2017) 96:94513. doi: 10.1103/PhysRevB.96.094513
10. Ashcroft NW. Hydrogen dominant metallic alloys: high temperature superconductors? Phys Rev Lett. (2004) 92:187002. doi: 10.1103/PhysRevLett.92.187002
11. Liu H, Naumov II, Hoffmann R, Ashcroft NW, Hemley RJ. Potential high-Tc superconducting lanthanum and yttrium hydrides at high pressure. Proc Natl Acad Sci USA. (2017) 114:6990–5. doi: 10.1073/pnas.1704505114
12. Wang H, Li X, Gao G, Li Y, Ma Y. Hydrogen-Rich superconductors at high pressures. Wiley Interdiscip Rev Comput Mol Sci. (2018) 8:1–13. doi: 10.1002/wcms.1330
13. Li Y, Hao J, Liu H, Tse JS, Wang Y, Ma Y. Pressure-stabilized superconductive yttrium hydrides. Sci Rep. (2015) 5:9948. doi: 10.1038/srep09948
14. Liu Y, Duan D, Huang X, Tian F, Li D, Sha X, et al. Structures and properties of osmium hydrides under pressure from first principle calculation. J Phys Chem C (2015) 119:15905–11. doi: 10.1021/acs.jpcc.5b03791
15. Liu Y, Duan D, Tian F, Liu H, Wang C, Huang X, et al. Pressure-Induced structures and properties in indium hydrides. Inorg Chem. (2015) 54:9924–8. doi: 10.1021/acs.inorgchem.5b01684
16. Peng F, Sun Y, Pickard CJ, Needs RJ, Wu Q, Ma Y. Hydrogen clathrate structures in rare earth hydrides at high pressures: possible route to room-temperature superconductivity. Phys Rev Lett. (2017) 119:107001. doi: 10.1103/PhysRevLett.119.107001
17. Geballe ZM, Liu H, Mishra AK, Ahart M, Somayazulu M, Meng Y, et al. Synthesis and stability of lanthanum superhydrides. Angew Chem Int Ed Engl. (2018) 57:688–92. doi: 10.1002/anie.201709970
18. Wigner E, Huntington HB. On the possibility of a metallic modification of hydrogen. J Chem Phys. (1935) 3:764–70. doi: 10.1063/1.1749590
19. Narayana C, Luo H, Orloff J, Ruoff AL. Solid hydrogen at 342 GPa: no evidence for an alkali metal. Nature (1998) 393:46–9. doi: 10.1038/29949
20. Dalladay-Simpson P, Howie RT, Gregoryanz E. Evidence for a new phase of dense hydrogen above 325 gigapascals. Nature (2016) 529:63. doi: 10.1038/nature16164
21. Tse JS, Yao Y, Tanaka K. Novel superconductivity in metallic SnH4 under high pressure. Phys Rev Lett. (2007) 98:117004. doi: 10.1103/PhysRevLett.98.117004
22. Ye X, Zarifi N, Zurek E, Hoffmann R, Ashcroft NW. High hydrides of scandium under pressure: potential superconductors. J Phys Chem C (2018) 122:6298–309. doi: 10.1021/acs.jpcc.7b12124
23. Ashcroft NW. Metallic hydrogen: a high-temperature superconductor? Phys Rev Lett. (1968) 21:1748. doi: 10.1103/PhysRevLett.21.1748
24. Bi T, Zarifi N, Terpstra T, Zurek E. The search for superconductivity in high pressure hydrides. arXiv [Preprint] arXiv (2018) 1806:163.
25. Gao G, Hoffmann R, Ashcroft NW, Liu H, Bergara A, Ma Y. Theoretical study of the ground-state structures and properties of niobium hydrides under pressure. Phys Rev B (2013) 88:184104. doi: 10.1103/PhysRevB.88.184104
26. Liu Y, Huang X, Duan D, Tian F, Liu H, Li D, et al. First-principles study on the structural and electronic properties of metallic HfH2 under pressure. Sci Rep. (2015) 5:11381. doi: 10.1038/srep11381
27. Shanavas KV, Lindsay L, Parker DS. Electronic structure and electron-phonon coupling in TiH2. Sci Rep. (2016) 6:28102. doi: 10.1038/srep28102
28. Wang L, Duan D, Yu H, Xie H, Huang X, Ma Y, et al. High-Pressure formation of cobalt polyhydrides: a first-principle study. Inorg Chem. (2018) 57:181–6. doi: 10.1021/acs.inorgchem.7b02371
29. Feng X, Zhang J, Gao G, Liu H, Wang H. Compressed sodalite-like MgH6 as a potential high-temperature superconductor. RSC Adv. (2015) 5:59292. doi: 10.1039/C5RA11459D
30. Zhang H, Jin X, Lv Y, Zhuang Q, Liu Y, Lv Q, et al. A novel stable hydrogen-rich SnH8 under high pressure. RSC Adv. (2015) 5:107637. doi: 10.1039/c5ra20428c
31. Xie Y, Li Q, Oganov AR, Wang H. Superconductivity of lithium-doped hydrogen under high pressure. Acta Crystallogr Sect C Struct Chem. (2014) 70:104–11. doi: 10.1107/S2053229613028337
32. Wiedner ES, Chambers MB, Pitman CL, Bullock RM, Miller AJM, Appel AM. Thermodynamic hydricity of transition metal hydrides. Chem Rev. (2016) 116:8655–92. doi: 10.1021/acs.chemrev.6b00168
33. Abe K. Hydrogen-Rich scandium compounds at high pressures. Phys Rev B (2017) 96:144108. doi: 10.1103/PhysRevB.96.144108
34. Scheler T, Marqués M, Konôpková Z, Guillaume CL, Howie RT, Gregoryanz E. High-Pressure synthesis and characterization of iridium trihydride. Phys Rev Lett. (2013) 111:215503. doi: 10.1103/PhysRevLett.111.215503
35. Zhuang Q, Jin X, Cui T, Ma Y, Lv Q, Li Y, et al. Pressure-Stabilized superconductive ionic tantalum hydrides. Inorg Chem. (2017) 56:3901–8. doi: 10.1021/acs.inorgchem.6b02822
36. Li X, Peng F. Superconductivity of pressure-stabilized vanadium hydrides. Inorg Chem. (2017) 56:13759–65. doi: 10.1021/acs.inorgchem.7b01686
37. Pépin CM, Dewaele A, Geneste G, Loubeyre P, Mezouar M. New iron hydrides under high pressure. Phys Rev Lett. (2014) 113:265504. doi: 10.1103/PhysRevLett.113.265504
38. Sumuraya TT, Atsuura YM, Hishidou TS, Guchi TO. First-Principles study on the structural and magnetic properties of iron hydride. J Phys Soc Jpn. (2012) 81:064707. doi: 10.1143/JPSJ.81.064707
39. Bazhanova ZG, Oganov AR, Gianola O. Fe–C and Fe–H systems at pressures of the earth's inner core. Phys Uspekhi (2012) 55:489–97. doi: 10.3367/UFNe.0182.201205c.0521
40. Li F, Wang D, Du H, Zhou D, Ma Y, Liu Y. Structural evolution of FeH4 under high pressure. RSC Adv. (2017) 7:12570–5. doi: 10.1039/C6RA25591D
41. Pépin CM, Geneste G, Dewaele A, Mezouar M, Loubeyre P. Synthesis of FeH5 : a layered structure with atomic hydrogen slabs. Science (2017) 357:382–5. doi: 10.1126/science.aan0961
42. Majumdar A, Tse JS, Wu M, Yao Y. Superconductivity in FeH5. Phys Rev B (2017) 96:201107. doi: 10.1103/PhysRevB.96.201107
43. Kvashnin AG, Kruglov IA, Semenok DV, Oganov AR. Iron superhydrides FeH5 and FeH6: stability, electronic properties, and superconductivity. J Phys Chem C (2018) 122:4731–6. doi: 10.1021/acs.jpcc.8b01270
44. Zhang S, Lin J, Wang Y, Yang G, Bergara A, Ma Y. Nonmetallic FeH6 under high pressure. J Phys Chem C (2018) 122:12022–8. doi: 10.1021/acs.jpcc.7b11643
45. Kim DY, Scheicher RH, Ahuja R. Predicted high-temperature superconducting state in the hydrogen-dense transition-metal hydride YH3 at 40K and 17.7GPa. Phys Rev Lett. (2009) 103:77002. doi: 10.1103/PhysRevLett.103.077002
46. Wang X, Andrews L. Neon matrix infrared spectra and DFT calculations of tungsten hydrides WHx(x = 1-4, 6). J Phys Chem A (2002) 106:6720–9. doi: 10.1021/jp025920d
47. Wang X, Andrews L, Infante I, Gagliardi L. Infrared spectra of the WH4(H2)4 complex in solid hydrogen. J Am Chem Soc. (2008) 130:1972–8. doi: 10.1021/ja077322o
48. Arnoult William J, McLellan RB. Thermodynamics of transition metal alloys. Acta Metall. (1973) 21:1397.
49. Strobel TA, Somayazulu M, Hemley RJ. Novel pressure-induced interactions in silane-hydrogen. Phys Rev Lett. (2009) 103:65701. doi: 10.1103/PhysRevLett.103.065701
50. Zaleski-Ejgierd P, Labet V, Strobel TA, Hoffmann R, Ashcroft NW. WHn under pressure. J Phys Condens Matter (2012) 24:155701. doi: 10.1088/0953-8984/24/15/155701
51. Wang X, Andrews L. Neon matrix infrared spectrum of WH6: a distorted trigonal prism structure. J Am Chem Soc. (2002) 124:5636–7. doi: 10.1021/ja020216w
52. Scheler T, Peng F, Guillaume CL, Howie RT, Ma Y, Gregoryanz E. Nanocrystalline tungsten hydrides at high pressures. Phys Rev B (2013) 87:184117. doi: 10.1103/PhysRevB.87.184117
53. Hooper J, Altintas B, Shamp A, Zurek E. Polyhydrides of the alkaline earth metals: a look at the extremes under pressure. J Phys Chem C (2013) 117:2982–92. doi: 10.1021/jp311571n
54. Zhong X, Wang H, Zhang J, Liu H, Zhang S, Song HF, et al. Tellurium hydrides at high pressures: high-temperature superconductors. Phys Rev Lett. (2016) 116:57002. doi: 10.1103/PhysRevLett.116.057002
55. Wang Y, Lv J, Zhu L, Ma Y. Crystal structure prediction via particle-swarm optimization. Phys Rev B Condens Matter Mater Phys. (2010) 82:94116. doi: 10.1103/PhysRevB.82.094116
56. Wang Y, Lv J, Zhu L, Ma Y. CALYPSO: a method for crystal structure prediction. Comput Phys Commun. (2012) 183:2063–70. doi: 10.1016/j.cpc.2012.05.008
57. Zhu L, Wang H, Wang Y, Lv J, Ma Y, Cui Q, et al. Substitutional alloy of Bi and Te at high pressure. Phys Rev Lett. (2011) 106:145501. doi: 10.1103/PhysRevLett.106.145501
58. Lv J, Wang Y, Zhu L, Ma Y. Predicted novel high-pressure phases of lithium. Phys Rev Lett. (2011) 106:15503. doi: 10.1103/PhysRevLett.106.015503
59. Li Y, Hao J, Liu H, Li Y, Ma Y. The metallization and superconductivity of dense hydrogen sulfide. J Chem Phys. (2014) 140:174712. doi: 10.1063/1.4874158
60. Zhu L, Liu H, Pickard CJ, Zou G, Ma Y. Reactions of xenon with iron and nickel are predicted in the earth's inner core. Nat Chem. (2014) 6:644–8. doi: 10.1038/nchem.1925
61. Yang G, Wang Y, Peng F, Bergara A, Ma Y. Gold as a 6p-element in dense lithium aurides. J Am Chem Soc. (2016) 138:4046–52. doi: 10.1021/jacs.5b11768
62. Zhang S, Zhu L, Liu H, Yang G. Structure and electronic properties of Fe2SH3 compound under high pressure. Inorg Chem. (2016) 55:11434–9. doi: 10.1021/acs.inorgchem.6b01949
63. Pfrommer BG, Côté M, Louie SG, Cohen ML. Relaxation of crystals with the quasi-newton method. J Comput Phys. (1997) 131:233–40. doi: 10.1006/jcph.1996.5612
64. Clark SJ, Segall MD, Pickard CJ, Hasnip PJ, Probert MJ, Refson K, et al. First principles methods using CASTEP. Zeitschrift fuer Krist (2005) 220:567–70. doi: 10.1524/zkri.220.5.567.65075
65. Perdew J, Chevary J, Vosko S, Jackson K, Pederson M, Singh D, et al. atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation. Phys Rev B Condens Matter Mater Phys. (1992) 46:6671. doi: 10.1103/PhysRevB.48.4978.2
66. Vanderbilt D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys Rev B (1990) 41:7892–5. doi: 10.1103/PhysRevB.41.7892
67. Pack JD, Monkhorst HJ. Special points for brillouin-zone integrations. Phys Rev B (1976) 13:5188. doi: 10.1103/PhysRevB.16.1748
68. Giannozzi P, Baroni S, Bonini N, Calandra M, Car R, Cavazzoni C, et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J Phys Condens Matter (2009) 21:395502. doi: 10.1088/0953-8984/21/39/395502
69. Deshpande VT, Pawar R. X-Ray Determination of the thermal expansion of tungsten. Curr Sci. (1962) 31:497–9.
70. Pickard CJ, Needs RJ. Structure of phase III of solid hydrogen. Nat Phys. (2007) 3:473–6. doi: 10.1038/nphys625
71. Zhuang Q, Jin X, Lv Q, Li Y, Shao Z, Liu Z, et al. Investigation of superconductivity in compressed vanadium hydrides. Phys Chem Chem Phys. (2017) 19:26280. doi: 10.1039/C7CP03435K
72. Segall M, Shah R, Pickard C, Payne M. Population analysis of plane-wave electronic structure calculations of bulk materials. Phys Rev B (1996) 54:16317–20. doi: 10.1103/PhysRevB.54.16317
73. Marck SC. Site percolation and random walks on d-Dimensional kagomé lattices. J Phys Math Gen. (1998) 31:3349–460. doi: 10.1088/0305-4470/31/15/010
74. Atwood JL. Kagomé lattice: a molecular toolkit for magnetism. Nat Mater. (2002) 1:91. doi: 10.1038/nmat740
75. Schlickum U, Decker R, Klappenberger F, Zoppellaro G, Klyatskaya S, Auwärter W, et al. Chiral kagomé lattice from simple ditopic molecular bricks. J Am Chem Soc. (2008) 130:11778–82. doi: 10.1021/ja8028119
76. McMahon JM, Ceperley DM. Ground-State structures of atomic metallic hydrogen. Phys Rev Lett. (2011) 106:165302. doi: 10.1103/PhysRevLett.106.165302
77. Duan D, Liu Y, Tian F, Li D, Huang X, Zhao Z, et al. Pressure-Induced metallization of dense (H2S)2H2 with high-Tc superconductivity. Sci Rep. (2015) 4:6968. doi: 10.1038/srep06968
78. Zhang S, Wang Y, Zhang J, Liu H, Zhong X, Song HF, et al. Phase diagram and high-temperature superconductivity of compressed selenium hydrides. Sci Rep. (2015) 5:15433. doi: 10.1038/srep15433
80. Bardeen J, Cooper LN, Schrieffer JR. Theory of superconductivity. Phys Rev. (1957) 108:1175. doi: 10.1103/PhysRev.108.1175
81. Allen PB, Dynes RC. Transition temperature of strong-couple superconductors reanalyzed. Phys Rev B (1975) 12:905–22. doi: 10.1103/PhysRevB.12.905
82. Yu S, Jia X, Frapper G, Li D, Oganov AR, Zeng Q, et al. Pressure-Driven formation and stabilization of superconductive chromium hydrides. Sci Rep. (2015) 5:17764. doi: 10.1038/srep17764
Keywords: tungsten hydride, superconductivity, high pressure, structure prediction, density functional theory
Citation: Zheng S, Zhang S, Sun Y, Zhang J, Lin J, Yang G and Bergara A (2018) Structural and Superconducting Properties of Tungsten Hydrides Under High Pressure. Front. Phys. 6:101. doi: 10.3389/fphy.2018.00101
Received: 20 June 2018; Accepted: 23 August 2018;
Published: 19 September 2018.
Edited by:
Gang Zhang, Institute of High Performance Computing (A*STAR), SingaporeReviewed by:
Chengyong Zhong, Chongqing University of Posts and Telecommunications, ChinaWeifeng Li, Shandong University, China
Copyright © 2018 Zheng, Zhang, Sun, Zhang, Lin, Yang and Bergara. 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: Guochun Yang, yanggc468@nenu.edu.cn orcid.org/0000-0003-3083-472X
Aitor Bergara, a.bergara@ehu.eus