Crystal Structure and Superconductivity of PH3 at High Pressures

Jan 20, 2016 - Dense Hydrocarbon Structures at Megabar Pressures. Hanyu Liu , Ivan I. Naumov , and Russell J. Hemley. The Journal of Physical Chemistr...
1 downloads 17 Views 1MB Size
Subscriber access provided by UNIV OSNABRUECK

Article 3

Crystal Structure and Superconductivity of PH at High Pressures Hanyu Liu, Yinwei Li, Guoying Gao, John S Tse, and Ivan Ivanovich Naumov J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b12009 • Publication Date (Web): 20 Jan 2016 Downloaded from http://pubs.acs.org on January 22, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Crystal Structure and Superconductivity of PH3 at High Pressures Hanyu Liu1,2*, Yinwei Li3, Guoying Gao4, John S. Tse2,5*, Ivan I. Naumov1 1

Geophysical Laboratory, Carnegie Institution of Washington, Washington, DC 20015, USA 2

Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2, Canada

3

School of Physics and Electronic Engineering, Jiangsu Normal University, Xuzhou 221116, China 4

State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao 066004, China

5

State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China

Abstract We have performed a systematic structure search on solid PH3 at high pressures using the particle swarm optimization method. At 100-200 GPa, the search led to two structures which along with others have P-P bonds. These structures are structurally and chemically distinct from those predicted for the high pressure superconducting H2S phase which has a different topology (i.e., does not contain S-S bonds). Phonon and electron-phonon coupling calculations indicate that both structures are dynamically stable and superconducting. The pressure dependence and critical temperature for the monoclinic (C2/m) phase of 83 K at 200 GPa are in excellent agreement with a recent experimental report.

*Corresponding author: [email protected] or [email protected]



1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction Theoretical calculations based on the Bardeen-Cooper-Schrieffer (BCS)1 phonon mediated theory of superconductivity indicate that highly compressed solid hydrogen2 and main group hydrogen-rich materials3 can possess high superconducting critical temperatures (Tc) in predicted high-pressure metallic states owing to the high frequency vibrations of the hydrogen component. This proposal has motivated active theoretical and experimental research to search for superconductivity in these hydrogen-containing materials under pressure.4-7 Theoretical calculations based on the BCS theory on a variety of hydrides show the Tc’s indeed are very high.8-13 The predictions contrast with the common view that the BCS-based Tc cannot be much higher than 40 K14. Experimentally, superconductivity had been observed in SiH4 at 120 GPa with a Tc of 17 K6 and in BaReH9, Tc = 7 K above 100 GPa.15 A theoretical calculation predicted that H2S may exhibit superconductivity with a Tc as high as ~80 K above 100 GPa.16 Subsequent experimental electrical conductivity and magnetic measurements of a material formed from compressing H2S showed a Tc of 203 K at 200 GPa17. It was proposed that the superconducting phase is not stoichiometric H2S but likely to be H3S, in which self-consistent harmonic approximation calculations predicted a Tc of 194 K at 200 GPa.18 However, another theoretical study also suggested the formation of a Im-3m (H2S)2(H2)19 structure, with a Tc of 191-204 K at 200 GPa. So far, the exact nature of the superconducting phase of compressed H2S remains elusive. Very recently PH3, an analogue of H2S, was reported to possess high-temperature superconductivity with a Tc of ~100 K at high pressures.20 The study of PH3 may provide an additional example for the understanding of the electron-phonon coupling mechanism responsible for high Tc in compressed H2S and other hydrogen-rich

2

ACS Paragon Plus Environment

Page 2 of 13

Page 3 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

materials. In this work, we investigated possible high pressure crystal structures of PH3 using first-principles electronic structure methods. Two possible candidates were found. Furthermore, the calculated Tc’s based on BCS theory on one of the structures agree well with experiment, providing additional support for the proposal that these high critical temperatures can be realized by a phonon-mediated mechanism. Methods The search for low-energy crystalline structures of PH3 at high pressure were performed using particle swarm optimization method implemented in the CALYPSO code.21 This method has been applied successfully in a broad range of crystalline systems ranging from elemental solids to binary and ternary compounds.22-26 Structure searches were performed from 50 to 200 GPa using models consisting of 1-4 formula unit. In general, the search was terminated after the generation of distinct 1500 structures. In additional, we have also performed prediction calculations on PHx (x=1, 2, 4, 5) at 200 GPa. Structural optimizations, enthalpies, electronic structures, and phonons were calculated using the density-functional electronic structure theory (DFT) with the Perdew-Burke-Ernzerhof (PBE)27 generalized gradient approximation. Phonon dispersion and electron-phonon coupling calculations were performed with density functional perturbation theory. Ultrasoft pseudopotentials for P and H were employed with a kinetic energy cutoff of 80 Ry. A q-mesh of 6×6×3 and k-mesh of 24×24×12 for C2/m structure in the first Brillouin zone was used in the electron-phonon coupling (EPC) calculations. The superconductivity calculations were performed with the Quantum-ESPRESSO package.17



3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Results and Discussion

Structure prediction calculations were performed at 100, 200 and 250 GPa. From the calculated enthalpies a structure with P212121 symmetry was found to be the most stable phase below 210 GPa but it transforms to a C2/m structure at higher pressures (Fig. 1). Since PH3 has a high concentration of hydrogen, it is important to consider the contribution of the zero point energy (ZPE) to the relative stability of the candidate structures. Using the calculated harmonic frequencies, the inclusion of ZPE energies lowers the transition pressure of P212121 to C2/m significantly, from 210 GPa to 115 GPa. Although both predicted structures contain P-P bonds (Fig. 2), details in the crystal structures are quite different. At 100 GPa, the P-P bond length in the P212121 structure is 2.13 Å. The hydrogen atoms form single P-H bonds in the P-P units and also bridged between neighboring P-P with P-H bond-lengths from 1.38-1.58 Å at 100 GPa. In comparison, in the C2/m structure the P atoms form polymeric chains decorated with P-H bonds. At 200 GPa, the P-P bond-lengths in the chain alternate between 2.06 and 2.09 Å and P-H bond-lengths are 1.41-1.42 Å. At 100 GPa, the ZPEs of P212121 and C2/m are 0.878 eV/f.u. and 0.802 eV/f.u., respectively. In the P212121 structure, the P-H stretch frequencies that contribute most the ZPE are almost 500 cm-1 or 20% higher than C2/m.(Fig. S1) There are additional P-H-P bend modes at 1250-1500 cm-1 in P212121 but absent in C2/m. Notably, both predicted structures are fundamentally different from H3S, the purportedly superconducting phase of compressed H2S at 200 GPa as in the predicted Im-3m H3S

4

ACS Paragon Plus Environment

Page 4 of 13

Page 5 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

structure consisted only of S-H bonds and there were no S-S bonds.19 This suggests that the electron-phonon interactions responsible for high Tc superconductivity in PH3 differ from that in H3S.

Fig. 1. Calculated enthalpy of the predicted lowest energy phases of PH3 with and without zero point energy corrections.

Fig. 2. Predicted crystal structures of PH3.The large and small spheres represent the P and H atoms, respectively.



5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Fig. 3. Calculated electronic band structures and total density of states (DOS) of (a) P212121 at 100 GPa and (b) C2/m at 200 GPa in PH3. We first analyze the electronic band structures of P212121 and C2/m structured PH3 (Fig. 3). The band dispersions show the metallic behavior of P212121 PH3 at 100 GPa is the result of indirect band overlap. In comparison, the C2/m structure is a genuine metal with parabolic dispersive bands crossing the Fermi level in the A-Γ-Z symmetry axis. Furthermore the electronic bands along L-M-A are also parabolic and cross the Fermi level with steep slopes indicating mobile electrons along these symmetry directions.



6

ACS Paragon Plus Environment

Page 6 of 13

Page 7 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Fig. 4. The phonon spectrum, phonon density of states and Eliashberg phonon spectral function for phase P212121 structure at 100 GPa (a) and C2/m structures at 200 GPa C2/m (b). Red solid circles in (bottom) show the electron-phonon coupling with the radius proportional to their respective strength.

To investigate the superconductive behaviour, phonon and electron-phonon coupling calculations were performed. The phonon band structures reveal no imaginary frequency in P212121 and C2/m structures indicating both are dynamically stable. Phonon band structures and projected density of vibration states within the stable pressure range of these structures are shown in Fig. 4. At 100 GPa (Fig. 4a), the electron-phonon coupling parameter (λ) for the P212121 structure is 0.53 with an average phonon frequency ωln of 990 K defined by the Allen Dynes equation,28 an extension of the McMillan theory for strong coupling.29

Using a nominal Coulomb

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

pseudopotential parameter (µ*) of 0.13 the estimated superconducting critical temperature, Tc, is 9 K for the P212121 structure. At 200 GPa, the calculated λ for the C2/m structure of 1.45 is much larger than that of the P212121 phase. With an average phonon frequency ωln of 826 K, again using µ*=0.13, the estimated Tc is 83 K for the C2/m phase.

Fig. 5. The calculated Tc in PH3 along experimental data.20

We investigated the pressure dependence of the critical temperature. The results summarized in Fig. 5 show the Tc for PH3 increases almost linearly with pressure within the range 100 to 200 GPa; a trend in good agreement with experimental data.20 It is clearly seen that calculated Tc of C2/m is better in agreement with experiment results than P212121, suggesting the C2/m phase may be a good candidate of the experimental structure. This finding suggests the calculated Tc for C2/m phase agrees better with the measurements than for P212121 phase and is consistent with the result in Fig. 1 indicating enhanced stability of C2/m phase owing to ZPE. The larger λ and

8

ACS Paragon Plus Environment

Page 8 of 13

Page 9 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

higher Tc of the C2/m phase can be traced to the strong Kohn anomalies that considerably soften the P-P phonons along the high symmetry paths (L-M-A-Γ-Z-V).The absence of similar (associated with the S-S bonds) vibrations in H3S means that the particular electron-phonon interactions leading to a relatively high Tc in each case differ, but they still are more effective in H3S. The calculated λ and ωln for the Im-3m structure of H3S are, respectively, 2.19 and 1334 K at 200 GPa.19 These values are to be compared with C2/m PH3 where λ and ωln are 1.46 and 826 K. A smaller λ and ωln in PH3 as compared to H3S leads to a lower Tc in PH3.

Fig. 6. The formation enthalpy of PHx (X=1, 2, 3, 4, 5) at 200 GPa. We also calculated the formation enthalpy of P-H system at high pressures. It is interesting to note that no thermodynamically stable phase relative to elemental P and solid H2 is found until 200 GPa (Fig. 6). The result is similar to a recent preprint arXiv paper on pnictogen hydrides30 but differs from another report.31 In the latter work, the stability of P-H system was investigated systematically and the results suggest PH2 is the most stable phase with a negative formation enthalpy at high

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

pressure.31 Furthermore, both enthalpy and superconducting calculations suggest the PH2 may be a possible contestant for the experimentally observed phase.31 A more recent theoretical study32 reported the high pressure phase diagram of PHn (n=1, 2, 3, 4, 5, 6) from structures determined by the Minima Hopping Method.33, 34 The results suggested the lowest energy phases of PH1,2,3 are also superconducting with TC comparable to experiments. In conclusion, theoretical structural search of the crystal structures of PH3 at 100-200 GPa reveal two low energy structures containing P-P bonds. The structures found differ from those predicted for H3S in which contain no heavier element bonds. Electronic band structures and phonon calculations show both the structures are metallic and dynamically stable. Electron-phonon coupling calculations confirm that the C2/m structure is a superconductor with a Tc of 83 K at 200 GPa.

Acknowledgements Work at Carnegie was supported by EFree, an Energy Frontier Research Center funded by the DOE, Office of Science, Basic Energy Sciences under Award No. DE-SC-0001057 (salary support for H.L). The infrastructure and facilities used at Carnegie were supported by NNSA Grant No. DE-NA-0002006, CDAC. A portion of the calculations were performed using computing resources provided by WestGrid, Compute Canada and XSEDE. J.T. and H.L. acknowledge the National Science Foundation of China (11474126) and support from the University of Saskatchewan research computing group and the use of the HPC resources (Plato machine). Y.L acknowledges support from the National Natural Science Foundation of China under Grant Nos. 11204111 and 11404148, the Natural Science Foundation of Jiangsu

10

ACS Paragon Plus Environment

Page 10 of 13

Page 11 of 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

province under Grant No. BK20130223, and the PAPD of Jiangsu Higher Education Institutions.

References: 1.

Bardeen, J.; Cooper, L. N.; Schrieffer, J. R., Microscopic Theory of Superconductivity. Phys. Rev.

1957, 106, 162-164. 2.

Ashcroft, N. W., Metallic Hydrogen: A High-temperature Superconductor? Phys. Rev. Lett. 1968,

21, 1748-1749. 3. Ashcroft, N., Hydrogen Dominant Metallic Alloys: High Temperature Superconductors? Phys. Rev. Lett. 2004, 92, 187002-187005. 4.

Strobel, T. A.; Goncharov, A. F.; Seagle, C. T.; Liu, Z.; Somayazulu, M.; Struzhkin, V. V.; Hemley, R.

J., High-pressure Study of Silane to 150 GPa. Phys. Rev. B 2011, 83, 144102-144115. 5.

Strobel, T. A.; Ganesh, P.; Somayazulu, M.; Kent, P.; Hemley, R. J., Novel Cooperative Interactions

and Structural Ordering in H2S-H2. Phys. Rev. Lett. 2011, 107, 255503-255506. 6.

Eremets, M.; Trojan, I.; Medvedev, S.; Tse, J.; Yao, Y., Superconductivity in Hydrogen Dominant

Materials: Silane. Science 2008, 319, 1506-1509. 7.

Chen, X.-J.; Wang, J.-L.; Struzhkin, V. V.; Mao, H.-k.; Hemley, R. J.; Lin, H.-Q., Superconducting

Behavior in Compressed Solid SiH4 with a Layered Structure. Phys. Rev. Lett. 2008, 101, 077002-077005. 8.

Li, Y.; Hao, J.; Liu, H.; John, S. T.; Wang, Y.; Ma, Y., Pressure-Stabilized Superconductive Yttrium

Hydrides. Sci. Rep. 2015, 5, 9948-9952. 9.

Gao, G.; Hoffmann, R.; Ashcroft, N. W.; 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-184107. 10. 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. Nati. Acad. Sci. USA 2010, 107, 15708-15711. 11. Wang, S.; Mao, H.-k.; Chen, X.-J.; Mao, W. L., High Pressure Chemistry in The H2-SiH4 System. Proc. Nati. Acad. Sci. USA 2009, 106, 14763-14767. 12. Gao, G.; Oganov, A. R.; Bergara, A.; Martinez-Canales, M.; Cui, T.; Iitaka, T.; Ma, Y.; Zou, G., Superconducting High Pressure Phase of Germane. Phys. Rev. Lett. 2008, 101, 107002-107005. 13. Tse, J.; Yao, Y.; Tanaka, K., Novel Superconductivity in Metallic SnH4 under High Pressure. Phys. Rev. Lett. 2007, 98, 117004-117007. 14. Schilling, A.; Cantoni, M.; Guo, J.; Ott, H., Superconductivity above 130 K in The Hg-Ba-Ca-Cu-O System. Nature 1993, 363, 56-58. 15. Muramatsu, T.; Wanene, W. K.; Somayazulu, M. S.; Vinitsky, E.; Chandra, D.; Strobel, T. A.; Struzhkin, V. V.; Hemley, R. J., Metallization and Superconductivity in The Hydrogen-rich Ionic Salt BaReH9. J. Phys. Chem. C 2015, 119, 18007-18013. 16. Li, Y.; Hao, J.; Liu, H.; Li, Y.; Ma, Y., The Metallization and Superconductivity of Dense Hydrogen Sulfide. J. Chem. Phys 2014, 140, 174712-174719.

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

17. Drozdov, A.; Eremets, M.; Troyan, I.; Ksenofontov, V.; Shylin, S., Conventional Superconductivity at 203 K at High Pressures. Nature 2015, 525, 73-76. 18. Errea, I.; Calandra, M.; Pickard, C. J.; Nelson, J.; Needs, R. J.; Li, Y.; Liu, H.; Zhang, Y.; Ma, Y.; Mauri, F., High-Pressure Hydrogen Sulfide from First Principles: A Strongly Anharmonic Phonon-mediated Superconductor. Phys. Rev. Lett. 2015, 114, 157004-157007. 19. 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-6971. 20. Drozdov, A. P.; Eremets, M. I.; Troyan, I. A., Superconductivity above 100 K in PH3 at High Pressures. arXiv 2015. 1508, 06224-06227 21. Wang, Y.; Lv, J.; Zhu, L.; Ma, Y., Crystal Structure Prediction Via Particle-swarm Optimization. Phys. Rev. B 2010, 82, 094116-094119. 22. Zhong, X.; Wang, Y.; Peng, F.; Liu, H.; Wang, H.; Ma, Y., Pressure Stabilization of Long-missing Bare C6 Hexagonal Rings in Binary Sesquicarbides. Chem. Sci. 2014, 5, 3936-3940. 23. Liu, H.; Wang, H.; Ma, Y., Quasi-molecular and Atomic Phases of Dense Solid Hydrogen. J. Phys. Chem. C 2012, 116, 9221-9226. 24. Zhu, L.; wang, H.; Wang, Y.; Lv, J.; Ma, Y.; Cui, Q.; Ma, Y.; Zou, G., Substitutional Alloy of Bi and Te at High Pressure. Phys. Rev. Lett 2011, 106, 145501-145504. 25. Lv, J.; Wang, Y.; Zhu, L.; Ma, Y., Predicted Novel High-pressure Phases of Lithium. Phys. Rev. Lett 2011, 106, 015503-015507. 26. Zhang, X.; Qin, J.; Liu, H.; Zhang, S.; Ma, M.; Luo, W.; Liu, R.; Ahuja, R., Pressure-induced Zigzag Phosphorus Chain and Superconductivity in Boron Monophosphide. Scientific Reports 2015, 5, 8761-8764. 27. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 28. Allen, P. B.; Dynes, R. C., Transition Temperature of Strong-coupled Superconductors Reanalyzed. Phy. Rev. B 1975, 12, 905-922. 29. Dynes, R., McMillan's Equation and The Tc of Superconductors. Solid State Commun. 1972, 10, 615-618. 30. Fu, Y.; Du, X.; Zhang, L.; Peng, F.; Zhang, M.; Pickard, C. J.; Needs, R. J.; Singh, D. J.; Zheng, W.; Ma, Y., High-Pressure Phase Stability and Superconductivity of Pnictogen Hydrides and Chemical Trends for Compressed Hydrides. arXiv preprint arXiv 2015. 1510, 04415-04418 31. Shamp, A.; Terpstra, T.; Bi, T.; Falls, Z.; Avery, P.; Zurek, E., Decomposition Products of Phosphine under Pressure: PH2 Stable and Superconducting? arXiv 2015. 1509, 05455-05468 32. Flores-Livas, J. A.; Amsler, M.; Heil, C.; Sanna, A.; Boeri, L.; Profeta, G.; Wolverton, C.; Goedecker, S.; Gross, E., Superconductivity in Metastable Phases of Phosphorus-hydride Compounds under High Pressure. arXiv 2015, 1512, 02132-02139. 33. Goedecker, S., Minima hopping: An Efficient Search Method for the Global Minimum of The Potential Energy Surface of Complex Molecular Systems. J. Chem. Phys. 2004, 120, 9911-9917. 34. Amsler, M.; Goedecker, S., Crystal Structure Prediction using The Minima Hopping Method. J. Chem. Phys. 2010, 133, 224104-224116.



12

ACS Paragon Plus Environment

Page 12 of 13

Page 13 of 13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC



13

ACS Paragon Plus Environment