Structures and Properties of Osmium Hydrides under Pressure from

The pressure-induced new structures and properties of osmium hydrides were systematically explored in a wide pressure range 0–300 GPa using ab initi...
0 downloads 6 Views 4MB Size
Article pubs.acs.org/JPCC

Structures and Properties of Osmium Hydrides under Pressure from First Principle Calculation Yunxian Liu, Defang Duan, Xiaoli Huang, Fubo Tian, Da Li, Xiaojing Sha, Chao Wang, Huadi Zhang, Ting Yang, Bingbing Liu, and Tian Cui* State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun 130012, P. R. China S Supporting Information *

ABSTRACT: The pressure-induced new structures and properties of osmium hydrides were systematically explored in a wide pressure range 0−300 GPa using ab initio methods. Three stable stoichiometries, that is, OsH, OsH3, and OsH6, are predicted above 50 GPa. The above hydrides exhibit metallic character with the notable band structures exception of OsH6. It is interesting to note that the phase P21/c of hydrogen-rich OsH6 adopts intriguing structures with H2 units. The electron−phonon coupling calculations indicate that the superconducting critical temperature (Tc) values of Fm-3m-OsH is 2.1 K at 100 GPa. Comparing to pure Os, the addition of hydrogen is in favor of improving the superconducting temperature.



INTRODUCTION Hydrides have attracted a great deal of attention for several years due to the fact that they not only could be hightemperature superconductor at a relatively low pressure, but also can be used as a potential material for hydrogen storage. For example, sulfur hydrides was found having a high Tc with 190 K.1−3 And (H2)4CH4 can be used as a hydrogen storage material with 33.4 wt % molecular hydrogen.4 As we know, many transition metals can react with hydrogen to form hydrides under ambient conditions.5,6 However, platinum group metals, including Pt, Pd, Os, Ir, Ru, and Rh, do not form hydrides at ambient pressure, which owes to the fact that platinum-group metals are normally inert and resistant to react with other elements at normal conditions. High pressure can provide a different route to new material synthesis and manipulating novel physical properties. In the condensed and high pressure environments, the atoms close together and bonding patterns change, resulting in new materials that cannot be synthesized at normal conditions. With the pressure increasing, the chemical potential of hydrogen rises steeply and then hydrogen reacts with metals to form metal hydrides.7−10 The platinum-group metal hydrides, such as Rh, Pd, Ir, and Pt hydrides have been experimentally synthesized under high pressure.7,8,11−13 In addition, some of them have also been explored in theory.14−17 In recent years, the platinum-group metal hydrides have been of great scientific interest, which owes to their applications in the hydrogen economy and the potential substances of achieving metallic hydrogen at modest pressure. For example, RhH2 is a high volumetric hydrogen density compound, which theoretical volumetric hydrogen densities reach 163.7 g H2/L.11 And PtH is a superconductor with Tc = 12−16 K at 105 GPa.14 © XXXX American Chemical Society

As a platinum group metal, osmium is the densest and stiffest metal among transition metals, having a very large bulk modulus comparable to that of diamond. And it maintains the close-packed structure (P63/mmc) up to high pressures. To date, the various aspects of osmium and its compounds have been experimentally and theoretically investigated.18−20 However, no osmium hydrides are known at high pressures. Based on the applications of platinum-group metal hydrides materials, as well as the novel physical properties motivated us to perform a thorough search on osmium hydrides under high pressure. In addition, for most known D-metal hydrides, the metal atoms generally form a closed-packed host lattice in which hydrogen atoms occupy the octahedral or tetrahedral interstitial sites. A further motivation factor for the present investigations is to check whether the osmium hydrides are the above structural type or form some exotic structures. In the present work, we employed a first-principles global structural optimization method to extensively explore and search for the structures of hitherto unknown osmium hydrides at zero temperature. The phase stabilities of several stoichiometric binary osmium hydrides OsHn (n = 1−8) under various pressures were investigated. Then we obtain thermodynamically stable structures and identify a large number of stable phases which may be synthesized under proper pressure conditions. Particularly, we examined in detail their structures and calculated their dynamical stability, the corresponding electronic band structures, bonding patterns, as well as the superconductivity of Fm-3m-OsH. Received: April 20, 2015 Revised: June 11, 2015

A

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C



COMPUTATIONAL METHOD We have employed the evolutionary algorithm USPEX, implemented in the USPEX code (Universal structure predictor: Evolutionary Xtallography),21−23 to search for the most stable stoichiometries and structures of OsHn (n = 1−8). The energetic and electronic structure calculations were carried out using density functional theory within the Perdew−Burke− Ernzerhof (PBE) parametrization of the generalized gradient approximation (GGA),24 as implemented in the Vienna ab initio simulation package (VASP).25 The adopted all-electron projector-augmented wave PAW26 pseudopotentials of osmium and hydrogen treat 5d66s2 and 1s1 electrons as valence electrons. And for hydrogen and osmium, the cutoff radii were taken 0.8 and 2.5 au, respectively. A plane-wave cutoff energy of 800 eV and appropriate Monkhorst−Pack-point meshes were used for all structures to achieve that the enthalpy calculations are converged to better than 1 meV/atom. Phonon calculations were performed by the direct supercell method using a supercell approach27 with the PHONOPY code.28 Electron−phonon couplings for superconducting properties were studied by within the framework of the linear-response theory as implemented in the Quantum ESPRESSO package.29 The convergence tests concluded that suitable calculation parameters are 100 Ry for kinetic energy cutoff, 20 × 20 × 20 k-point sampling mesh, and 5 × 5 × 5 q-mesh in Brillouin zone for Fm-3m-OsH.



Figure 1. Enthalpies of formation (ΔHf, in meV/atom) of OsHn (n = 1−8) with respect to Os and H2 at selected pressures. The symbols on the solid line represent that the hydride is stable at the corresponding pressure.

RESULTS AND DISCUSSION

Here, we perform the evolutionary simulations with considering simulation sizes ranging from one to four for each stoichiometry OsHn (n = 1−8) formula units (f.u.) at 50, 100, 150, 200, 250, and 300 GPa, and find many different relaxed structures. The computed enthalpies of formation, ΔHf, of OsHn (n = 1−8) with respect to elemental solid Os and H2 at the selected pressures are provided in the convex hull, as shown in Figure 1. Moreover, the known structures of P63/mmc for Os and P63/m, C2/c, and Cmca for H230 are considered depending on the applied pressures. Structures which are thermodynamically stable relative to decomposing into other hydrides and/or pure elements are lying on the convex hull, whereas the structures above the convex hull are metastable (Figure 1). Principally, these structures can be synthesized in experiment. At the pressure of 50 GPa, the OsH6 becomes the preferred stable phase and also the only stoichiometry which can be formed, as depicted in Figure 1. At 100 GPa, OsH and OsH6 are thermodynamically stable species, and OsH6 has the most negative enthalpy of formation. With increasing pressure to 150 GPa, the enthalpy of the OsH was predicted as being the global minimum, with OsH3 and OsH6 falling on the convex hull. At the higher pressure of 200 GPa, OsH remains the most negative ΔHf, while other thermodynamically stable phase is OsH3. When the pressure reaches 300 GPa, OsH is the only most favorable hydride. We also calculated the enthalpy for the most stable structures and constructed the pressure−composition phase diagram of Os−H system, as depicted in Figure 2. It shows the stability ranges of Os, OsH, OsH3, OsH6 and H2, respectively. We can see that the P21/c of OsH6 is the most stable structure at 38 GPa. Then when the pressure reaches 155 GPa, OsH6 (space group Fdd2) becomes unstable and decomposes into OsH3 (Cmm2) and H2 (C2/c). At 140 GPa, OsH3 (space group Cmm2) is the most stable compound, and at elevated pressure

Figure 2. Pressure−composition phase diagram for Os−H system in the pressure of 0−300 GPa.

of 246 GPa it separates into OsH (Fm-3m) and H2 (C2/c). For OsH, the structure Fm-3m structure becomes favorable in the pressure range from 94 to 300 GPa. We can see that, for the stable stoichiometries, the general trend of n (n denotes the number of H in OsHn) is declining with pressure, which is similar to the case of Ir−H system.16 That is to say, the H concentration appears to decrease with increasing pressure. As we know, the electronegativity value (2.18) of hydrogen is comparable or smaller than those of Os and Ir (2.2), which may lead to the gradual instability of stable H-rich stoichiometries with the pressure increasing. The high pressure crystal structures of OsH, OsH3 and OsH6 were plotted in Figure 3. All the structures are represented in the form of a conventional cell with the addition of symmetry. At 300 GPa, OsH adopts a face-centered cubic Fm-3m (4 f.u./cell) structure, which can be B

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 3. Crystal structures of the energetically stable Os−H system. (a) Fm-3m-OsH at 300 GPa, (b) Cmm2-OsH3 at 250 GPa, (c) P21/c-OsH6 at 100 GPa, and (d) Fdd2-OsH6 at 200 GPa. Purple and blue balls denote Os and H atoms, respectively. Yellow balls represent H2 units. The unit cell is drawn with black solid lines.

Figure 4. Phonon dispersion curves and phonon density of states (PDOS) projected for Fm-3m-OsH, Cmm2-OsH3, P21/c-OsH6, and Fdd2-OsH6 at 100, 200, 100, and 150 GPa, respectively.

viewed as NaCl-type (see Figure 3a). At the pressure of 250 GPa, OsH3 prefers the stable structure of orthorhombic Cmm2 (8 f.u./cell) and we carried out three times expanding cell in caxis direction in the Cmm2 phase to better describe its structural feature, as shown in Figure 3b. We note that the structure segregates into layers, in which one is consistent with

Os and H atoms, and the other is formed by only H atoms. For OsH6, we predict monoclinic P21/c and orthorhombic Fdd2 at 100 and 150 GPa, respectively. Remarkably, P21/c (4 f.u./cell) can be viewed as being composed of H2 units, H atoms and Os atoms, as shown in Figure 3c. The Fdd2 (8 f.u./cell) structure consists of H atoms and Os atoms, in which Os atoms can be C

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 5. (a−c) Electronic band structures and (e−h) calculated total and partial density of states (PDOS) of the Fm-3m-OsH, Cmm2-OsH3, P21/cOsH6, and Fdd2-OsH6 at different pressures. The vertical dashed line at zero is the Fermi energy level.

addition of a certain element may make OsH6 become a good hydrogen storage material at mild condition. Dynamic property of the phase is one of the basic criteria when considering the structure stability. Figure 4 presents the calculated phonon dispersion curves and projected phonon density of states (PHDOS) for Fm-3m (OsH), Cmm2 (OsH3), P21/c, and Fdd2 (OsH6) at accessible pressures. There is no sign of imaginary frequencies in the entire Brillouin zone, which indicates that the above phases are dynamical stability in their corresponding studied pressure range. A striking feature of PHDOS is the two separate regions: the high frequency region and the low frequency modes (below 10 THz) are related to H and Os atoms, which attributes to osmium much heavier than hydrogen atom. To investigate the electronic properties of Os−H system, we explore the electronic band structures and partial density of states (PDOS) for Fm-3m (OsH), Cmm2 (OsH3), P21/c, and Fdd2 (OsH6) at 100, 200, 100, and 150 GPa, respectively, as plotted in Figure 5. Note that, for Fm-3m OsH and Cmm2

seen to form a distorted simple orthorhombic structure, as seen in Figure 3d. Detail parameters of these structures at different pressures are listed in Table S1 in the Supporting Information. Figure S1 in the Supporting Information shows the enthalpy curves of OsH6 with respect to Os and H2 as a function of pressure. OsH6 become stable at 38 GPa, meaning that it can be synthesized in the recent diamond anvil cell experiment. Under compressing, Fdd2 stands out as lowest-enthalpy phase in 104−300 GPa. It is noted that RhH2 was synthesized above 8 GPa and released H at 4 GPa.11 Based on Rh and Os belong to platinum group metals, the OsH6 with hydrogen content of 3.1 wt % may be have same property in terms of hydrogen storage. In addition, Mg base hydrogen storage alloys perform excellent properties on hydrogen storage at proper conditions, for example, Mg3Pr, Mg (In, Al) and CeH2.73−MgH2−Ni.31−33 The studying of structures and deuterium desorption for Ca3Mg2Ni13 can provide a blueprint for the mechanism of hydrogen desorption.34 The above reaches inspire us that the D

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

point, the Tc was estimated by using the standard Allen−Dynes modified McMillan equation.35

OsH3, the conduction and the valence bands are overlapped, suggesting that the two phases are metallic (as shown in Figure 5a and b), while for P21/c and Fdd2 of OsH6, they are semiconducting phases due to their energy band gaps (Figure 5c and d). A visual inspection of the PDOS of Fm-3m OsH and Cmm2 OsH3 (Figure 5e and f) shows that Os-d orbits mainly dominate the density of states at the Fermi level. Whereas for P21/c and Fdd2 OsH6, the PDOSs show significant hybridization between the Os-d and H-s orbits in the energy range (Figure 5g and h). To further visualize the interactions between the atoms in the polymeric hydrogenic sublattice for P21/cOsH6, we calculated electron localization function (ELF), as depicted in Figure 6c. The ELF value located between H2 unit

Tc =

⎡ ⎤ 1.04(1 + λ) exp⎢ ⎥ ⎣ λ − μ*(1 + 0.62λ) ⎦ 1.2

ω log

This equation has been found to be highly accurate for materials with λ < 1.5. The electron−phonon coupling strength (EPC) λ, the logarithmic average phonon frequency ωlog, and the electronic DOS at the Fermi level N(Ef) were calculated to explore the possible superconductivity. According to our calculations, with the Coulomb parameter set μ* = 0.1, the EPC parameter λ becomes 0.42, ωlog is 347.84 K, N(Ef) is 4.22 (states/spin/Ry/unit cell), and the Tc is 2.1 K (at 100 GPa). Moreover, to our best knowledge, the Tc of Os is only 0.65 K. Therefore, we considered that the addition of hydrogen is in favor of improving the superconducting critical temperature. Besides the above stable structures, we also observed some interesting metal-stable phases, and their structural motifs are presented in Figure 7. At 100 and 150 GPa, we predicted two monoclinic phases with Cm (4 f.u./cell) and P21 (4 f.u./cell) for OsH5, which are composed of Os, H atoms and H3 units. However, the morphology of H3 is different: the H3 is linear in Cm phase, but it becomes bent in P21 structure. OsH7 adopts a monoclinic C2/m (4 f.u./cell) and a orthorhombic Ima2 (2 f.u./cell) at 150 and 300 GPa, respectively. Both phases consist of H, Os atoms and H2 units, in which the H−H distances are 0.87 and 0.79 Å. For OsH8, a structure with P21/c symmetry containing four formula units (Z = 4) was predicted at 200 GPa. The P21/c contains H, Os atoms and H3 units, and the H3 units form triangles with the distance of 0.87 Å, 0.90 Å, and 0.96 Å. The map of the electron localization function (ELF) for P21/c-OsH8 (Figure 6d) shows localized electron density within the H3 ring (ELF over 0.9), implying that localized spin paired electron density exist.



CONCLUSION In summary, we have employed the evolutionary algorithm USPEX in combination with first-principles calculations to investigate the high-pressure crystal structures of Os−H system (H-rich). Three stoichiometries (OsH, OsH3, and OsH6) are found stable, and OsH6 become stable relative to pure elements above 38 GPa. The hydrogenic sublattices of osmium hydrides are composed of H atoms and H2 and H3 units. The band structures of stable phases reveal that Fm-3m (OsH) and Cmm2 (OsH3) are metallic phases. Moreover, the estimated super-

Figure 6. Calculated ELF of (a, b) P21/c-OsH6 with isosurface value of 0.8 and (c, d) P21/c-OsH8 with isosurface value of 0.8.

is close to 1.0, which is significantly larger than the ELF value between Os and the nearest H atom (ELF < 0.5), implying that the former form strong covalent bonding in H2 unit, while no bonds were present between the Os and H atoms. The density of states for the Fm-3m OsH phase at the Fermi level has a significant value, suggesting that it may be a superconductor at accessible conditions. Further examining this

Figure 7. Crystal structures of the metastable Os−H system at different pressures. Purple and blue balls denote Os and H atoms, respectively. Yellow balls represent H2, H3 units. E

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

(10) Ye, X.; Hoffmann, R.; Ashcroft, N. W. Theoretical Study of Phase Separation of Scandium Hydrides under High Pressure. J. Phys. Chem. C 2015, 119, 5614−5625. (11) Li, B.; Ding, Y.; Kim, D. Y.; Ahuja, R.; Zou, G.; Mao, H. K. Rhodium dihydride (RhH2) with high volumetric hydrogen density. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 18618−18621. (12) Scheler, T.; Marqués, M.; Konôpková, Z.; Guillaume, C. L.; Howie, R. T.; Gregoryanz, E. High-Pressure Synthesis and Characterization of Iridium Trihydride. Phys. Rev. Lett. 2013, 111, 215503− 215507. (13) Scheler, T.; Degtyareva, O.; Marqués, M.; Guillaume, C. L.; Proctor, J. E.; Evans, S.; Gregoryanz, E. Synthesis and properties of platinum hydride. Phys. Rev. B 2011, 83, 214106−214110. (14) Zhang, C.; Chen, X. J.; Lin, H. Q. Phase transitions and electron-phonon coupling in platinum hydride. J. Phys.: Condens. Matter 2012, 24, 035701−035707. (15) Zhou, X.-F.; Oganov, A. R.; Dong, X.; Zhang, L.; Tian, Y.; Wang, H.-T. Superconducting high-pressure phase of platinum hydride from first principles. Phys. Rev. B 2011, 84, 054543−054547. (16) Zaleski-Ejgierd, P. High-pressure formation and stabilization of binary iridium hydrides. Phys. Chem. Chem. Phys. 2014, 16, 3220− 3229. (17) Gao, G.; Wang, H.; Zhu, L.; Ma, Y. Pressure-Induced Formation of Noble Metal Hydrides. J. Phys. Chem. C 2012, 116, 1995−2000. (18) Armentrout, M. M.; Kavner, A. Incompressibility of osmium metal at ultrahigh pressures and temperatures. J. Appl. Phys. 2010, 107, 093528. (19) Liu, K.; He, D.-W.; Zhou, X.-L.; Chen, H.-H. First-principles study of structural and thermodynamic properties of osmium. Phys. B 2011, 406, 3065−3069. (20) Minisini, B.; Roetting, J.; Tsobnang, F. Elastic and thermodynamic properties of OsSi, OsSi2 and Os2Si3. Comput. Mater. Sci. 2008, 43, 812−817. (21) Oganov, A. R.; Glass, C. W. Crystal structure prediction using ab initio evolutionary techniques: principles and applications. J. Chem. Phys. 2006, 124, 244704−244718. (22) Oganov, A. R.; Lyakhov, A. O.; Valle, M. How Evolutionary Crystal Structure Prediction Works and Why. Acc. Chem. Res. 2011, 44, 227−237. (23) Lyakhov, A. O.; Oganov, A. R.; Stokes, H. T.; Zhu, Q. New developments in evolutionary structure prediction algorithm USPEX. Comput. Phys. Commun. 2013, 184, 1172−1182. (24) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (25) Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15−50. (26) Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758−1775. (27) Parlinski, K.; Li, Z.; Kawazoe, Y. First-principles determination of the soft mode in cubic ZrO2. Phys. Rev. Lett. 1997, 78, 4063−4067. (28) Togo, A.; Oba, F.; Tanaka, I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures. Phys. Rev. B 2008, 78, 134106−134114. (29) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; Dal Corso, A.; de Gironcoli, S.; Fabris, S.; Fratesi, G.; Gebauer, R.; Gerstmann, U.; Gougoussis, C.; Kokalj, A.; Lazzeri, M.; Martin-Samos, L.; Marzari, N.; Mauri, F.; Mazzarello, R.; Paolini, S.; Pasquarello, A.; Paulatto, L.; Sbraccia, C.; Scandolo, S.; Sclauzero, G.; Seitsonen, A. P.; Smogunov, A.; Umari, P.; Wentzcovitch, R. M. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys.: Condens. Matter 2009, 21, 395502. (30) Pickard, C. J.; Needs, R. J. Structure of phase III of solid hydrogen. Nat. Phys. 2007, 3, 473−476. (31) Ouyang, L. Z.; Yang, X. S.; Dong, H. W.; Zhu, M. Structure and hydrogen storage properties of Mg3Pr and Mg3PrNi0.1 alloys. Scr. Mater. 2009, 61, 339−342.

conducting critical temperature Tc of OsH reaches 2.1 K at 100 GPa by application of the Allen−Dynes modified McMillan equation. Our findings help to shed light on understanding the rich and complex crystal structures of osmium hydrides, which have broad implications for further exploring other platinumgroup metal hydrides.



ASSOCIATED CONTENT

S Supporting Information *

Calculated structural parameters of our predicted stable structures for Os−H system at their corresponding pressures; enthalpy curves (relative to Os and H2) for OsH6 in the pressure range from 5 to 300 GPa. The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.5b03791.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Tel/fax: +86-431-85168825. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Basic Research Program of China (No. 2011CB808200), Program for Changjiang Scholars and Innovative Research Team in University (No. IRT1132), National Natural Science Foundation of China (Nos. 11204100, 51032001, 11074090, 10979001, 51025206, 11104102, and 11404134), National Found for Fostering Talents of basic Science (No. J1103202), China Postdoctoral Science Foundation (2012M511326, 2013T60314, and 2014M561279), Graduate Innovation Fund of Jilin University (2015102). Parts of calculations were performed in the High Performance Computing Center (HPCC) of Jilin University.



REFERENCES

(1) Duan, D.; Liu, Y.; Tian, F.; Li, D.; Huang, X.; Zhao, Z.; Yu, H.; Liu, B.; Tian, W.; Cui, T. Pressure-induced metallization of dense (H2S)2H2 with high-Tc superconductivity. Sci. rep. 2014, 4, 6968− 6973. (2) Drozdov, A.; Eremets, M.; Troyan, I. Conventional superconductivity at 190 K at high pressures. 2014, arXiv:1412.0460. arXiv.org e-Print archive. http://arxiv.org/abs/1412.0460. (3) Duan, D.; Huang, X.; Tian, F.; Li, D.; Liu, Y.; Ma, Y.; Liu, B.; Cui, T. Pressure-induced decomposition of solid hydrogen sulfide. Phys. Rev. B 2015, 91, 180502. (4) Mao, W. L.; Mao, H. K. Hydrogen storage in molecular compounds. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 708−710. (5) Smithson, H.; Marianetti, C.; Morgan, D.; Van der Ven, A.; Predith, A.; Ceder, G. First-principles study of the stability and electronic structure of metal hydrides. Phys. Rev. B 2002, 66, 144107− 144116. (6) Miwa, K.; Fukumoto, A. First-principles study on 3d transitionmetal dihydrides. Phys. Rev. B 2002, 65, 155114−155120. (7) Antonov, V. Phase transformations, crystal and magnetic structures of high-pressure hydrides of d-metals. J. Alloys Compd. 2002, 330, 110−116. (8) Degtyareva, O.; Proctor, J. E.; Guillaume, C. L.; Gregoryanz, E.; Hanfland, M. Formation of transition metal hydrides at high pressures. Solid State Commun. 2009, 149, 1583−1586. (9) 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−2992. F

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (32) Wang, H.; Zhong, H.; Ouyang, L.; Liu, J.; Sun, D.; Zhang, Q.; Zhu, M. Fully Reversible De/hydriding of Mg Base Solid Solutions with Reduced Reaction Enthalpy and Enhanced Kinetics. J. Phys. Chem. C 2014, 118, 12087−12096. (33) Ouyang, L. Z.; Yang, X. S.; Zhu, M.; Liu, J. W.; Dong, H. W.; Sun, D. L.; Zou, J.; Yao, X. D. Enhanced Hydrogen Storage Kinetics and Stability by Synergistic Effects of in Situ Formed CeH2.73and Ni in CeH2.73-MgH2-Ni Nanocomposites. J. Phys. Chem. C 2014, 118, 7808−7820. (34) Zhang, Q.; Sun, D.; Zhang, J.; Latroche, M.; Ouyang, L.; Zhu, M. Structure and Deuterium Desorption from Ca3Mg2Ni13 Deuteride: A Neutron Diffraction Study. J. Phys. Chem. C 2014, 118, 4626−4633. (35) Allen, P. B. Transition temperature of strong-coupled superconductors reanalyzed. Phys. Rev. B 1975, 12, 905−922.

G

DOI: 10.1021/acs.jpcc.5b03791 J. Phys. Chem. C XXXX, XXX, XXX−XXX