SCIENTIFIC REPpRTS

Received: 02 April 2015 Accepted: 16 July 2015 Published: 12 August 2015

Hydrogen segregation and its roles in structural stability and metallization: silane under pressure

Wenwen Cui1, Jingming Shi1, Hanyu Liu2, Yansun Yao2'3, Hui Wang1, Toshiaki litaka4 & Yanming Ma1,s

We present results from first-principles calculations on silane (SiH4) under pressure. We find that a three dimensional P-3 structure becomes the most stable phase above 241 GPa. A prominent structural feature, which separates the P-3 structure from previously observed/predicted SiH4 structures, is that a fraction of hydrogen leaves the Si-H bonding environment and forms segregated H2 units. The H2 units are sparsely populated in the system and intercalated with a polymeric Si-H framework. Calculations of enthalpy of formation suggest that the P-3 structure is against the decomposition into Si-H binaries and/or the elemental crystals. Structural stability of the P-3 structure is attributed to the electron-deficient multicenter Si-H-Si interactions when neighboring silicon atoms are linked together through a common hydrogen atom. Within the multicenter bonds, electrons are delocalized and this leads to a metallic state, possibly also a superconducting state, for SiH4. An interesting outcome of the present study is that the enthalpy sum of SiH4 (P-3 structure) and Si (fcc structure) appears to be lower than the enthalpy of disilane (Si2H6) between 200 and 300 GPa (for all previously predicted crystalline forms of Si2H6), which calls for a revisit of the stability of Si2H6 under high pressure.

The pursuit of metallic hydrogen and its high-temperature superconductivity is of great interest to the condensed-matter community1'2. Producing metallic hydrogen in bulk, however, requires a significant static pressure that exceeds the limit of current experimental capacities. Alternative approaches may be possible, as suggested by Ashcroft3, that some simple hydrides, with the average valence density comparable to pure hydrogen, could become metallic at significantly reduced pressures. Similar to pure hydrogen, the metallic states of these hydrides should also have large band overlapping, high Debye temperature, and strong electron-phononcoupling (EPC) that can all lead to phonon-mediated superconductivity. According to the Goldhammer-Herzfeld criterion, Group-IV hydrides, such as methane (CH4) and silane (SiH4), are potential candidates for this behavior3.

At ambient temperature, SiH4 solidifies at around 1.7 GPa into a four-fold molecular crystal. Upon increasing the pressure, one expects higher coordinatedstructures of SiH4 to become stable, by virtue of multicenter bonding. Previously, several high-pressure structures of SiH4 have been predicted by different groups, with the metallization pressure estimated from as low as 20 GPa to above 220 GPa4-8. Some of the predicted phases have been suggested to be superconducting, as well, with the estimate of

1State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China. 2Department of Physics

and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan, S7N 5E2, Canada. 3Canadian Light Source, Saskatoon, Saskatchewan, S7N 0x4, Canada. Computational Astrophysics Laboratory, RIKEN, 2-1

Hirosawa, Wako, Saitama 351-0198, Japan. 5Beijing Computational Science Research Center, Beijing 10084, China.

Correspondence and requests for materials should be addressed to Y.Y. (email: yansun.yao@usask.ca) or H.W.

(email: huiwang@jlu.edu.cn)

Tc ~ 16 K to 260 K. Experimental determinations of high-pressure structures of SiH4 have been carried out using various techniques. Degtyarevaet al. observed from powder x-ray diffraction patterns that SiH4 remains four-fold to around 25 GPa, and likely transforms to a new phase at higher pressures9. Chen et al. observed that SiH4 sample turns opaque above 27 GPa, which, in combination with measured Raman and IR reflectivity, led to the suggestion that SiH4 metallizes near 60 GPa10. Eremets et al. later reported the metallization of SiH4 to occur between 50 and 65 GPa based on electrical resistance measurement, and also found this phase to be superconducting(maximum Tc ~ 17.5 K)11. The crystal structure of the metallic phase (P63), however, was unusual: its volume shrinks nearly 60% from its precursor the P2llc structure at the phase transition. According to Degtyareva et al., this implies that the observed metallization actually arisesfrom contaminations in the system rather than from the SiH4 sample12. This was supported by the x-ray diffractionstudy of Hanfland et al.13 in which SiH4 was found to undergo pressure-inducedamorphization at pressures above 60 GPa recrystallizing at 90 GPainto an insulating/4jlastructure original proposed from ab initio calculations by Pickard and Needs14,15. At about the same time, Strobel et al. reported a comprehensive experimental study which confirmed Hanfland et al.'s observations16. Strobel et al. observed that SiH4 darkens above ~50 GPa, and attributed it to a partial loss of crystallinity. Above 100 GPa, SiH4 recrystallizes intothe I4jla (or I-42d) structure, and remainsin this structure to at least 150 GPa.

In addition to pure silane, the mixtures of silane and molecular hydrogen have also been investigated. In mixtures, hydrogen molecules are perturbed by SiH4. Such intermolecular interactions can destabilize the covalent bonds in H2, showing promise of obtaining a metallic state at moderate pressures. Almost concurrently, Strobel et al.17 and Wang et al.18 reported the synthesis of novel SiH4-H2complexes under high-pressure conditions. Using power X-ray diffraction, Strobel et al. observedthe formation of a stoichiometric SiH4(H2)2 crystal near 7 GPaand identified its structure as a cubic F-43m structure. Using Raman spectroscopy, Wang et al. studied the behavior of SiH4-H2 fluids and established a binary eutec-tic point at 72 mol% H2 and 6.1 GPa. Above the eutectic point the fluid crystallizes in a solid solution. The SiH4-H2 phase diagrams revealed in both studies are very similar. Detailed electronic structures and chemical bonding of the SiH4-H2 mixtures have been subsequently examined through the use of first-principles methods19-23.

While recent studies revealed a great deal of information on the behavior of SiH4 under pressure, the jury is still out on the initial proposal. As the research progresses, the reported metallization of SiH4 near 55 GPa is now facing challenges. Theory suggests that SiH4 reaches a metallic state at a much greater pressure (220 GPa), through a phase transition to a Pbcn structure8. Motivated by the significant interests in this subject, in the present study we employ theoretical structural searches and property predictions to investigate SiH4 in the post-I4jla region. We systematically searched for stable structures of SiH4 at pressures above 200 GPa, using a heuristic algorithm based on particle swarm optimization (calypso method)24,25. We predict a new structure of SiH4 above 241 GPa. This structure (P-3) is computed to be metallic and perhaps also superconducting as inferred from the electron-phonon coupling. Our prediction helps to address the phase diagram of SiH4 at high pressure; while its significance goes beyond. Disilane (Si2H6) is the next member in silicon hydride series (SinH2n+2), which was suggested in theory to become stable under pressure, stable with respect to decompositions into other binaries andlor the elemental crystals26. Using similar argument as for SiH4, solid Si2H6 was predicted to be metallic and superconducting. With the addition of the P-3 structure in the SiH4 phase diagram, however, according to our convex hull calculations all previously predicted Si2H6 structures appear to be thermodynamically unstable with respect to the decomposition of Si and SiH4 to at least 300 GPa. This therefore suggests a revisit of previous studies on solid Si2H6 and a reconsidering of its very stability at high pressure.

Methods

Our structure prediction is based on a global minimum search of the enthalpy surfaces obtained by ab initio calculations at a constant pressure, through CALYPSO (crystal structure analysis by particle swarm optimization) methodology and its same-name code24,25. For silane, disilane, silylene and lili-cane, structure predictions with upto 6 formula units (f.u.) in the simulation cell were performed at 200 and 300 GPa, respectively, however, for the complex Si5H18, simulation cells containing only 1 f.u. were considered. The ab initio calculations were performed using density functional theory27 within the Perdew-Burke-Ernzerhof (PBE) parameterization28 of the generalized gradient approximation (GGA) as implemented in the Vienna ab initio simulation package (VASP) code29. The all-electron projector-augmented wave (PAW) method30 was adopted with 3s23p2 and 1s1 as valence electrons for Si and H atoms, respectively. Planewave energy cutoffs of 500 eV and 700 eV, and uniform Monkhorst-Pack (MP) meshes31 for Brillouin zone (BZ) sampling with resolutions of 2n x 0.06 A-1 and 2n x 0.03 A-1 were employed in the structure predictions and subsequent calculations (e.g. of thermodynamic stability), respectively. The errors in the enthalpy differences of the structures studied in the structure predictions and the subsequent calculations were found to be in the order of 1 meVlatom and less than 1 meVlatom, respectively. The contribution of the electronic entropy to the enthalpy differences aroused by the MP method is less than 0.1 meVlatom, and negligible. The phonon density of states was calculated in harmonic approximation by the finite displacement method32 as implemented in the PHONOPY program33, and the vibrational free energy was estimated therefrom34. The underlying force constants were calculated in 2 x 2 x 1 supercells with 4 x 4 x 5 MP meshes for BZ sampling for the P-3 and Pbcn

Figure 1. Calculated enthalpies per formula unit (f. u.) for candidate structures of SiH4. Enthalpies of the Pbcn structure (Ref. [8]) were taken as the reference enthalpies. (Inset) Relative energetics for the I41/a (Ref. [14]), Pbcn and P-3 structures with the inclusion of the vibrational free energy.

Lattice parameters (A, deg) Atomic coordinates (Fractional)

P-3 a = 3.89 Si1 2d 0.3333 0.6667 0.3004

c = 4.73 Si2 2d 0.6667 0.3333 0.1629

Si3 1b 0 0 0.5

H1 2c 0 0 0.0793

H2 6g 0.3740 0.2957 0.9286

H3 6g 0.0886 0.6624 0.5724

H4 6g 0.7579 0.0050 0.2396

C2/m a = 5.17 Si1 4g 0.5 0.2130 0

b = 6.93 Si2 4i -0.8023 0 0.8728

c = 2.79 H1 8j -0.7719 -0.1785 0.2159

3 = 89.5 H2 8j -0.6255 0.1210 0.5903

H3 4i -0.0994 0 0.6447

H4 4i -0.5870 0 0.2256

H5 8j -0.6007 -0.3273 0.4034

Table 1. Optimized structural parameters for the P-3 and C2/m structures at 300 GPa.

structures and a 4 x4 x 4 MP mesh for the I41/a structure. Electron-phonon coupling (EPC) calculations were carried out by using density functional perturbation theory35 as implemented in the QUANTUM ESPRESSO package36 with a kinetic energy cutoff of 80 Ry. A 16 x 16 x 12 q-mesh in the first BZ was used with Gaussians of width of 0.05 Ry in the EPC calculations.

Results and Discussion

In Fig. 1, the calculated enthalpies (H = E + pV) of the candidate structures of SiH4 are compared between 200 to 300 GPa. In this pressure range, our structural searches successfully predicted a new structure with the P-3 space group, as well as revealed the previously observed/suggested structures, i.e., structures with P21/c, Fdd2, I41/a, C2/c and Pbcnspace groups8,9,14. The experimentally identified I41/a structure11,14,16 is calculated to be the most stable phase at the low-pressure end. The I41/a to Pbcn structural transformation is calculated to occur at about 225 GPa, in a good agreement with previous study8. Near 241 GPa, the newly predicted P-3 structure becomes the most stable phase of SiH4. Enthalpy of the P-3 structure is constantly lower than that of the C2/c structure14, by more than 0.08 eV per formula unit. As well, we predict a metastable C2/m structure that also appears to be more stable than the C2/c structure. The optimized structural parameters of the P-3 and C2/m structures are listed in Table 1. To account for the temperature effects, we estimated the vibrational free energies (at 300 K) for the I41/a, Pbcn and P-3 structures at four different pressures, using a harmonic method34 (Table 2). As shown

Pressure(GPa) I4Ja Pbcn P-3 Ai A2

150 1.322 1.251 1.126 -0.196 -0.125

200 1.404 1.339 1.233 -0.171 -0.106

250 1.472 1.408 1.326 -0.146 -0.082

300 1.533 1.472 1.406 -0.127 -0.066

Table 2. Calculated vibrational free energies (eV/f. u.) of the J41/a, Pbcn, and P-3 structures. Al (A2)

is the energy difference between the P-3 structure and the I41/a (Pbcn) structure.

in the inset of Fig. 1, the inclusion of the vibrational free energies does not alter the phase transition sequence, but shift down the two transition points in pressure (see explanation later). It is worth noting that the small proton mass may induce significant quantum nuclear effects, which can deviate the system from a harmonic description, and therefore, shift the two transition pressures further. A rigorous study of this should turn to explicit calculation of the free energy with anharmonic contributions from both the thermal and the quantum nuclear effects includedx37. This could be achieved by performing ab initio path-integral molecular dynamics simulations38, however, which is beyond the computational load we can afford nowadays. Here, we just speculate that the temperature effects will further stabilize the P-3 structure with respect to the other two structures. Moreover, the appearance of a molecular bonding environment in the P-3 structure (as will be discussed below) suggests that the van der Waals (vdW) interaction may play an important role in the phase transitions. To evaluate this effect, the enthalpies of I41/a, Pbcn and P-3 structureshave been recalculated usingoptB88-vdW functional39-42. The results indicate that the vdW interaction does not change the phase transition sequence, but increase the transition pressures of I41/a ^ Pbcn and Pbcn ^ P-3 from 225 and 241GPa to 242 and 285 GPa, respectively.

The P-3 structure is shown in Fig. 2(a,b) in two different views. Its unit cell contains 5 SiH4 formula units. The extended structure can be described as a polymeric Si-H framework intercalated by quasi-hexagonal layers of H2 units. Isolated H2 units in Group-IV hydrides are not entirely new - previous studies43-48 suggested them to exist in solid germane (GeH4), stannane (SnH4), and plumbane (PbH4). In SiH4, the H2 units are sparsely populated with negligible intermolecular interactions. This is in contrast to heavier Group-IV hydrides in which the H2 units are strongly correlated with a tendency for H2 - H2 pairing. To analyze the bonding environment, we employed the electron localization function (ELF)49 as a measure of relative electron localization with respect to a uniform electron gas. The ELF value represents a probability (0 to 1), with large values identifying places in the structure where electrons of opposite spins present, e.g., in cores, bonds, or lone pairs. As Fig. 2(c) shows, the intramolecular regions in the H2 units have a very large ELF value (close to 1) resulted from a high degree of electron paring. In fact, the ELF map of the H2 units is consistent to that of a gas-phase H2 molecule, suggesting that these H2 units are inact. This conclusion is resonating with the exceptionally short H-H distance in the H2 units, i.e., 0.75 A (calculated at 300 GPa), not so different from that of a H2 molecule.

In the Si-H framework, each silicon atom is coordinated to 12 nearest neighbor hydrogen atoms, which are then bridged to other silicon atoms. Overall, a third of the bridging hydrogen atoms are shared by 4 silicon atoms and the rest are shared by 3, yielding a Si5H18 stoichiometry. Such Si-H geometry goes far beyond the ubiquitous four-fold coordination of Group-IV elements, and on the face of it, one sees the enigma of silicon atoms forming more bonds (12) than the number that seems to be sufficient for their available valence electrons (4). This 'forbidden' increment of crystal coordination, however, is a commonplace in materials under high pressure, and can be explained by the schemes of electron-deficient multicenter bonding50. A quantum Atoms in Molecules (AIM) analysis51 (Table 3) reveals that the bridging hydrogen receive ~0.7 e- per atom from neighboring silicon (which is also seen in Fig. 2(c)). In contrast, the H2 units almost have no charge transfer from/to the Si-H framework. At this point, silicon atoms have much reduced s/p mixing in its delocalized multicenter bonding, compared with that in molecular SiH4 (which favors definite coordination geometry, as we know). A similar electron-deficient scenario is also encountered in the other two candidate structures, I4j/a and Pbcn, but the difference is that in these two structures the bridging hydrogen atoms are shared by only 2 silicon atoms8,14. Each silicon atom is coordinated to 8 hydrogen atoms; pairs of which bridge to four neighboring silicon atoms, forming two three-center-two-electron (3c-2e) bonds. This bonding arrangement is a prototypical electron-deficient geometry that is commonly known for the diborane molecule (B2H6)52.

Two interesting observations arise from a comparison of the P-3 structure to the I4j/a and Pbcnstructures. At first, the P-3 structure has higher coordination geometry which corresponds to a smaller volume. For example, at 250 GPa, the calculated volume for the I4j/a, Pbcnand P-3 structures-gradually decreases in a sequence of 13.59 A3/f.u., 13.43 A3/f.u. and 13.17 A3/f.u., respectively. Evidently, the P-3 structure would have the smallest pV work, which provides energetic advantages. Secondly, in the P-3 structure the bridging hydrogen atoms are constrained by more silicon atoms. This yields reduced frequencies for the Si...H...Si bending modes in the mid-frequency range (20-70 THz, at 150 GPa), which is the main reason why the vibrational free energy of the P-3 structure is constantly lower than the other two structures (Table 2). In the translational region, i.e., below 20 THz, the vibrational free energy

Figure 2. (a,b) Unit cell of the P-3 structure viewed from two different angles. The polymeric Si-H framework is shown as interconnected polyhedrons; Si and H atoms are located at the polyhedron centers and vertices, respectively. Isolated H2 units are represented by the dumbbells. (c) Calculated ELF map of the P-3 structure in the yz plane. Large and small balls represent the Si atoms and H atoms, respectively.

Atom Bader Charge (e )

Sil -2.7715

Si2 -2.6169

Si3 -2.5131

H1 +0.0307

H2 +0.7225

H3 +0.7606

H4 +0.7198

Table 3. Bader charges of the P-3 structure calculated at 300 GPa. Negative (positive) sign indicates an electron loss (gain) for the particular atom.

is mainly resulted from Si lattice modes, so it is very similar for all three structures. The P-3 structure also has high-frequency H-H vibron modes above 100 THz. These modes however do not have notable contributions to the vibrational free energy due to their low density of states. At 150 GPa, for example, the H-H vibrons only contribute 0.04 eV/f.u. to the vibrational free energy of the P-3 structure.

An important issue is worthy of investigation at this point. The appearance of the H2 units in the P-3 structure, as well as the short H-H contact, seems to indicate a trend of segregation toward the Si5H18 + H2 limit. As crystalline SiH4 is known to decompose between 50 ~ 100 GPa16, it is not unreasonable to speculate it to decompose again at higher pressures. To this end, we examined the phase stability of the P-3 structure with respect to the decompositions into other possible Si-H binaries and H2. Four stoichiometries were used for the binaries, namely, Si5H18, Si2H6 (disilane), SiH2 (silylene), and SiH (silicane). To find the most stable structure of these stoichiometries, we performed additional structure

100%H 0.08 *

* -0.04-

-0.08-

100%Si *

■ . ■ 200 GPa • 300 GPa

SiH, \\ ■ 1 - \\ : V .•

\\ Si,H

SiHj 1 ' I ■ 1 1 1 1

0.0 0.2 0.4 0.6

Silicon atomic content

Figure 3. Enthalpy of formation (△Hp with respect to elemental crystals of H2 and Si) for four Si-H stoichiometries calculated at 200 and 300 GPa. The abscissa x represents the atomic content (%) of Si in the stoichiometry. Stoichiometry with a negative △ Hf means that the binary phase is more stable than the elemental phases, while the convex hull of the △Hf values (solid lines) connects the thermodynamically stable stoichiometries.

Pressure (GPa) Si2H6 SiH4 Si H2 A Hf! (eV/Si2H6) AHf2 (eV^HJ

200 Cmcm I4Ja Fm-3m C2/c 0.249 0.042

250 Cmcm P-3 Fm-3m C2/c 0.256 0.016

300 C2/c P-3 Fm-3m Cmcm 12 0.383 0.072

300 Pm-3m P-3 Fm-3m Cmcm 12 0.388 0.077

Table 4. Enthalpies of formation △Hf1 of Si2H6 with respect to crystals of SiH4 and Si (3/2 SiH4 + 1/2 Si) and △Hf2 with respect to SiH4, Si and H2 (SiH4 + Si + H2 ) at 200, 250 and 300 GPa. The C2/c and Pm-3m structures of Si2H6 were previously predicted in Ref. 55, while the Cmcm structurepredicted in Ref. 26.

searches at 200 and 300 GPa. The lowest-enthalpy structure that arises from the search was employed as the candidate for each stoichiometry. In Fig. 3, the calculated enthalpy of formation (△ Hf, with respect to elemental crystal of H2 and Si) is shown for all four stoichiometries. Here a negative △ Hf means that the binary phase is more stable than the elemental crystals, while the convex hull of the △ Hf values (solid lines) connects the stable phases; in this case they are H2, Si, and SiH4. Clearly, SiH4 (P-3 structure) is the most stable Si-H stoichiometry at 200 and 300 GPa. The △Hf of the P-3 structure relative to Si5H18 and H2 (Cmca1253) structure is negative which confirms its stabilities.

It is perhaps not surprising that Si5H18, SiH2, and SiH crystals are predicted to be less stable than SiH4, since these species are either reactive intermediates or not known to exist experimentally. A significant finding is that the Si2H6 crystal appears to be unstable as well, with respect to the decomposition of Si + SiH4 crystal. To our best knowledge, Si2H6 is the only known higher silanes that can readily be prepared in laboratory, usually by the reaction of silicon chloride (Si2Cl6) with lithium aluminum hydride (LiAlH4). Due to the weaker Si-Si bond (226 kJ mol-1), Si2H6 decomposes slowly even at room temperature54. On the other hand, previous theoretical studiessuggest that the stability of Si2H6 can be enhanced by applying pressure26,55. A series of crystalline Si2H6 polymorphs, with the space group Cmcm, C2/c, and Pm-3m, have been predicted to be thermodynamically stable and superconducting at pressures above 190 GPa. On the contrary, new calculations using the P-3 structure as the most stable structure of SiH4 show that these Si2H6 polymorphs would all have positive △ Hfwith respect to the crystal of SiH4 + Si between 200 and 300 GPa (Table 4). This outcome therefore suggests a revisit on the stability of crystal Si2H6 at high pressures. It should be noted here that the choices of the Si-H binaries as the decomposition products are based upon the known Si-H stoichiometries observed at ambient pressure. It is possible that other Si-H binaries, in additional to what have been considered here, can be stabilized at high pressure.

It is worth noting that another recent theoretical study56 reports an interesting result on the dis-proportionation of GeH4 at high pressure. In this study, zero-point energies were estimated using the

Figure 4. (a,b) Calculated band structure, total DOS, and site-projected DOS near the Fermi level for the P-3 structure at 300 GPa. The Fermi level is at the energy origin. (c,d) Calculated phonon dispersion relations and site-projected phonon DOS for the P-3 structure at 300 GPa. In both (b,d), H' and H2 represent the bridging hydrogen and paired hydrogen, respectively.

harmonic approximation, which is very similar to the approach adopted in the present study. Under this consideration, a new stoichiometry for the Ge-H binaries, namely GeH3, is predicted to become energetically stable with respect to the decomposition of elemental crystals of Ge and H2 near 175 GPa. The GeH4, on the hand, was suggested to be unstable with respect to the decomposition of GeH3 and H2. These results are distinctly different for that of the present study in which the SiH4 is suggested to be the most stable stoichiometry for the Si-H binaries at high pressures. Experimentally, SiH4 is the most stable binary in its hydride series (SinH2n+2) at ambient pressure54. At high pressure, SiH4 decomposes near 50 GPa but recrystallizes above 100 GPa, and remains stable to at least 192 GPa, the highest pressure attempted in the experiments11,16. The P-3 structure is predicted to become stable near 177 GPa (with harmonic zero-point corrections), which is well within the experimental stability range for SiH4. The study of GeH3 on the other hand indicates that there are maybe new stoichiometries for this group, perhaps unprecedented at ambient pressures, may become stable at high pressures and this awaits to be discovered in future experiments.

Figure 4(a,b) show the calculated electronic band structure and projected density of states (DOS) for the P-3 structure at 300 GPa. It is of considerable interest that this structure is predicted to be metallic. The DOS reveals a pseudogap developed below the Fermi energy (EF). Similar pseudogaps were also found in calculations on high-pressure GeH4 and SnH4 but not on PbH444-48. This suggests that the electronic structure of the P-3 structure is not free-electron-like. Since in this case the Fermi level lies in the pseudogap, the DOS at the EF in the P-3 structure is substantially lower than the value in a free electron gas (3n/2Ef, where n is the density of electrons), which limits the possibility of achieving high superconductivity. It perhaps merits mention that the DOS characteristics of the P-3 structure is similar to that of the AlH3 at high pressure57. Since Al has one less electron than Si, the Fermi level of AlH3 lies in the lower reach of the pseudogap which also results in a low DOS value and ultimately non-superconducting behaviors58. In Fig. 4b, hydrogen atoms (mostly bridging ones) contribute to the DOS throughout the entire energy range. This indicates that the hydrogen and silicon electronic states are strongly mixed in the Si-H framework. If the states near the EF can be effectively coupled with the phonon modes, especially the high-frequency ones from hydrogen vibrons, a credible superconductivity can be expected3.

a F(a>) X{(o)

0 20 40 60 80 100 120

Frequencies (THz)

Figure 5. Eliashberg phonon spectral function a2F(w) and integrated EPC parameter (A (u) = Ju du'2ca F (u')/u') as a function of frequency for the P-3 structure at 300 GPa.

However, the H2 units in the P-3 structure are populated sparsely in the lattice so their DOS and phonon DOS are both very low (Figs 4b,d), compared with their counterparts in SnH4, GeH4 and PbH4. We therefore do not anticipate exceptional electron-phonon coupling in the P-3 structure at this point.

Figure 4(c,d) show the calculated phonon dispersion relations and phonon DOS for the P-3 structure at 300 GPa. The absence of imaginary frequency modes suggests that this structure is mechanically stable. An outstanding feature of the 'H2' intercalation in the structure is that the vibrational modes are divided into subsets in frequency. The H2 pairs occupy primarily two subsets; one induced by roton modes around 25 THz, and the other induced vibron modes around 118 THz. In these two subsets, the contribution from the Si-H framework is minor, consistent with the fact that the H2 pairs are inactive. The Si lattice modes dominate the low-frequency subset below 20 THz, while the bridging hydrogen atoms are in the intermediate subset between 30 and 70 THz.

A brief discussion on the possibility of phonon-mediated superconductivity in the P-3 structure is of interest. The methodology was based on an extension of the BCS model in which the attractive, electron-phonon interaction for each phonon mode is treated explicitly. The strength of this interaction is characterized by the electron-phonon coupling parameterX incorporating the contributions from all participating phonon modes. In Fig. 5, we present the Eliashberg phonon spectral function a2F(w) and the integrated A (w) (A (u) = f u du'2a2F (u')/u') as a function of frequency w at 300 GPa. The overall integrated A is 0.63, which is a moderate value for hydrogen-rich materials and comparable to the zero-pressure value for MgB259-61. The Si lattice modes contribute about 0.23 to the total A. While the low frequencies of the lattice modes were considered as a disadvantage of attaining superconductivity, their reasonably strong EPC compensates it. The intermediate-frequency H modes contribute 0.36 while the high-frequency vibron modes only contribute 0.04. The logarithmic average of the phonon frequency wlog is obtained as 1320 K. The electron-screened repulsive interaction is represented by the Coulomb pseudopotential ¡Í. For typical phonon-mediated superconductors, the values of fi between 0.1-0.13 are generally considered as reasonable3. The superconducting critical temperature Tc is estimated from the McMillan formula62 along with the Allen-Dynes correction (A (u) = fu du'2a2F (u')/uO63, using an empirical value of f* = 0.12. The estimated Tc of P-3 structure turns out to be 32 K at 300 GPa. A slightly improved estimate of the Tc of 35.1 K is made by directly solving the Eliashberg equations2 with the calculated a2F(w), using the methodology previously implemented by us64. We note that these two estimates of Tc are both based on isotropic gap equations and an empirical value of f *, which may not be adequate depending on the mechanism of the electron-phonon coupling in SiH4. A more accurate estimate of the Tc may be made using a fully anisotropic treatment of the gap equations, which encourages future study in this direction. Moreover, it should be pointed out that the quantum nuclear effects may deviate the system from a harmonic description and alter the estimated Tc. As manifested by a study of AlH3, anharmonicity of atomic motions can cause renormalization of the vibrational modes and suppress the superconductivity58, On the other hand, however, anharmonic vibrations were found to enhance the electron-phonon matrix elements, as in the case of disordered materials65.

Conclusions

We present a theoretical study of high-pressure phase transitions and the metallization of crystalline silane (SiH4). A new polymorph of SiH4 (space group: P-3) was predicted from structure searches using a heuristic algorithm based on particle swarm optimization (calypso methodology). This new phase

becomes thermodynamically stable at pressures above 241 GPa, replacing the previous suggested C2/c structure in the phase diagram. The P-3 structure is calculated to be metallic and also superconducting with an estimated Tc of 32 K at 300 GPa. A unique structure feature of the P-3 structure is the presence of H2 units that are intercalated with the Si-H framework. Electron-deficient multicenter bonding along the Si...H...Si connections results in the delocalization of valence electrons and accounts for the increment of crystal coordination in the Si-H framework. Calculations of enthalpy of formation suggest that the P-3 structure is stable with respective to the decomposition to other Si-H binaries and/or the elemental crystals.

References

1. Mao, H. & Hemley, R. Ultrahigh-pressure transitions in solid hydrogen. Rev. Mod. Phys. 66, 671-692 (1994).

2. McMahon, J. M., Morales, M. A., Pierleoni, C. & Ceperley, D. M. The properties of hydrogen and helium under extreme conditions. Rev. Mod. Phys. 84, 1607-1653 (2012).

3. Ashcroft, N. Hydrogen Dominant Metallic Alloys: High Temperature Superconductors? Phys. Rev. Lett. 92, 187002 (2004).

4. Feng, J. et al. Structures and Potential Superconductivity in SiH4 at High Pressure: En Route to 'Metallic Hydrogen. Phys. Rev. Lett. 96, 017006 (2006).

5. Yao, Y., Tse, J. S., Ma, Y. & Tanaka, K. Superconductivity in high-pressure SiH4. Europhys. Lett. 78, 37003 (2007).

6. Chen, X. et al. Superconducting Behavior in Compressed Solid SiH4 with a Layered Structure. Phys. Rev. Lett. 101, 077002 (2008).

7. Kim, D. Y. et al. Crystal structure of the pressure-induced metallic phase of SiH4 from ab initio theory. Proc. Natl. Acad. Sci. USA 105, 16454-16459 (2008).

8. Martinez-Canales, M. et al. Novel Structures and Superconductivity of Silane under Pressure. Phys. Rev. Lett. 102, 087005 (2009).

9. Degtyareva, O. et al. Crystal structure of SiH4 at high pressure. Phys. Rev. B 76, 064123 (2007).

10. Chen, X. et al. Pressure-induced metallization of silane. Proc. Natl. Acad. Sci. USA 105, 20-23 (2008).

11. Eremets, M. I., Trojan, I. A., Medvedev, S. A., Tse, J. S. & Yao, Y. Superconductivity in hydrogen dominant materials: silane. Science 319, 1506-1509 (2008).

12. Degtyareva, O., Proctor, J. E., Guillaume, C. L., Gregoryanz, E. & Hanfland, M. Formation of transition metal hydrides at high pressures. Solid State Commun. 149, 1583-1586 (2009).

13. Hanfland, M., Proctor, J. E., Guillaume, C. L., Degtyareva, O. & Gregoryanz, E. High-Pressure Synthesis, Amorphization, and Decomposition of Silane. Phys. Rev. Lett. 106, 095503 (2011).

14. Pickard, C. & Needs, R. High-Pressure Phases of Silane. Phys. Rev. Lett. 97, 045504 (2006).

15. Pickard, C. J. & Needs, R. J. Structures at high pressure from random searching. Phys. status solidi 246, 536-540 (2009).

16. Strobel, T. A. et al. High-pressure study of silane to 150 GPa. Phys. Rev. B 83, 144102 (2011).

17. Strobel, T. a., Somayazulu, M. & Hemley, R. J. Novel pressure-induced interactions in silane-hydrogen. Phys. Rev. Lett. 103, 2-5 (2009).

18. Wang, S., Mao, H., Chen, X.-J. & Mao, W. L. High pressure chemistry in the H2-SiH4 system. Proc. Natl. Acad. Sci. USA 106, 14763-14767 (2009).

19. Chen, X.-Q., Wang, S., Mao, W. L. & Fu, C. L. Pressure-induced behavior of the hydrogen-dominant compound SiH4 (H2 )2 from first-principles calculations. Phys. Rev. B 82, 104115 (2010).

20. Li, Y. et al. Superconductivity at approximately 100 K in dense SiH4(H2)2 predicted by first principles. Proc. Natl. Acad. Sci. USA 107, 15708-15711 (2010).

21. Michel, K., Liu, Y. & Ozolins, V. Atomic structure and SiH4 -H2 interactions of SiH4 (H2)2 from first principles. Phys. Rev. B 82, 174103 (2010).

22. Yao, Y. & Klug, D. D. Silane plus molecular hydrogen as a possible pathway to metallic hydrogen. Proc. Natl. Acad. Sci. USA 107, 20893-20898 (2010).

23. Yim, W. L., Tse, J. S. & Iitaka, T. Pressure-induced intermolecular interactions in crystalline silane-hydrogen. Phys. Rev. Lett. 105, 1-4 (2010).

24. Wang, Y., Lv, J., Zhu, L. & Ma, Y. Crystal structure prediction via particle-swarm optimization. Phys. Rev. B 82, 094116 (2010).

25. Wang, Y., Lv, J., Zhu, L. & Ma, Y. CALYPSO: A method for crystal structure prediction. Comput. Phys. Commun. 183, 2063-2070 (2012).

26. Flores-Livas, J. A. et al. High-Pressure Structures of Disilane and Their Superconducting Properties. Phys. Rev. Lett. 108, 117004 (2012).

27. Kohn, W. & Sham, L. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133 (1965).

28. Perdew, J., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

29. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558-561 (1993).

30. Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).

31. J., M. H. & D., P. J. Special points for Brillonin-zone integrations. Phys. Rev. B 13, 5188 (1976).

32. Parlinski, K., Li, Z. & Kawazoe, Y. First-Principles Determination of the Soft Mode in Cubic ZrO2. Phys. Rev. Lett. 78, 4063-4066 (1997).

33. 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 78, 134106 (2008).

34. Pavone, P., Baroni, S. & de Gironcoli, S. phase transition in tin: A theoretical study based on density-functional perturbation theory. Phys. Rev. B 57, 10421-10423 (1998).

35. Baroni, S., de Gironcoli, S. & Dal Corso, A. Phonons and related crystal properties from density-functional perturbation theory. Rev. Mod. Phys. 73, 515-562 (2001).

36. Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).

37. Chen, J., Ren, X., Li, X.-Z., Alfe, D. & Wang, E. On the room-temperature phase diagram of high pressure hydrogen: An ab initio molecular dynamics perspective and a diffusion Monte Carlo study. J. Chem. Phys. 141, 24501 (2014).

38. Li, X.-Z. et al. Classical and quantum ordering of protons in cold solid hydrogen under megabar pressures. J. Phys. Condens. Matter 25, 85402 (2013).

39. Dion, M., Rydberg, H., Schroder, E., Langreth, D. C. & Lundqvist, B. I. Van der Waals Density Functional for General Geometries. Phys. Rev. Lett. 92, 246401 (2004).

40. Klimes, J., Bowler, D. R. & Michaelides, A. Van der Waals density functionals applied to solids. Phys. Rev. B 83, 195131 (2011).

41. Román-Pérez, G. & Soler, J. M. Efficient Implementation of a van der Waals Density Functional: Application to Double-Wall Carbon Nanotubes. Phys. Rev. Lett. 103, 096102 (2009).

42. Klimes, J., Bowler, D. R. & Michaelides, A. Chemical accuracy for the van der Waals density functional. J. Phys. Condens. Matter 22, 022201 (2010).

43. Martinez-Canales, M., Bergara, A., Feng, J. & Grochala, W. Pressure induced metallization of Germane. J. Phys. Chem. Solids 67, 2095-2099 (2006).

44. Gao, G. et al. Superconducting High Pressure Phase of Germane. Phys. Rev. Lett. 101, 107002 (2008).

45. Tse, J., Yao, Y. & Tanaka, K. Novel Superconductivity in Metallic SnH4 under High Pressure. Phys. Rev. Lett. 98, 117004 (2007).

46. Gao, G. et al. High-pressure crystal structures and superconductivity of Stannane (SnH4). Proc. Natl. Acad. Sci. USA 107, 1317-1320 (2010).

47. Gonzalez-Morelos, P., Hoffmann, R. & Ashcroft, N. W. Segregation into layers: a general problem for structural instability under pressure, exemplified by SnH4. Chemphyschem 11, 3105-3112 (2010).

48. Zaleski-Ejgierd, P., Hoffmann, R. & Ashcroft, N. W. High Pressure Stabilization and Emergent Forms of PbH4. Phys. Rev. Lett. 107, 037002 (2011).

49. Becke, A. D. & Edgecombe, K. E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 92, 5397 (1990).

50. Grochala, W., Hoffmann, R., Feng, J. & Ashcroft, N. W. The chemical imagination at work in very tight places. Angew. Chem. Int. Ed. Engl. 46, 3620-3642 (2007).

51. Bader, R. F. W. Atoms in Molecules: A Quantum Thory Vol. 22 (Clarendon Press, Oxford, UK, 1990).

52. Yao, Y. & Hoffmann, R. BH3 under pressure: leaving the molecular diborane motif. J. Am. Chem. Soc. 133, 21002-21009 (2011).

53. Pickard, C. J., Martinez-Canales, M. & Needs, R. J. Density functional theory study of phase IV of solid hydrogen. Phys. Rev. B 85, 214114 (2012).

54. Chambers, C. & Holliday, A. K. Modern inorganic chemistry: Group IV 175-176 (Butterworth, London, UK, 1975).

55. Jin, X. et al. Superconducting high-pressure phases of disilane. Proc. Natl. Acad. Sci. USA 107, 9969-9973 (2010).

56. Abe, K. & Ashcroft, N. W. Quantum disproportionation: The high hydrides at elevated pressures. Phys. Rev. B 88, 174110 (2013).

57. Goncharenko, I. et al. Pressure-Induced Hydrogen-Dominant Metallic State in Aluminum Hydride. Phys. Rev. Lett. 100, 045504 (2008).

58. Rousseau, B. & Bergara, A. Giant anharmonicity suppresses superconductivity in AlH3 under pressure. Phys. Rev. B 82, 104504 (2010).

59. Bohnen, K.-P., Heid, R. & Renker, B. Phonon Dispersion and Electron-Phonon Coupling in MgB2 and AlB2. Phys. Rev. Lett. 86, 5771-5774 (2001).

60. Choi, H., Roundy, D., Sun, H., Cohen, M. & Louie, S. First-principles calculation of the superconducting transition in MgB2 within the anisotropic Eliashberg formalism. Phys. Rev. B 66, 020513 (2002).

61. Singh, P. From E2g to Other Modes: Effects of Pressure on Electron-Phonon Interaction in MgB2. Phys. Rev. Lett. 97, 247002 (2006).

62. McMillan, W. Transition temperature of strong-coupled superconductors. Phys. Rev. 167, 331 (1968).

63. Allen, P. B. & Dynes, R. C. Superconductivity at very strong coupling. J. Phys. C Solid State Phys. 8, L158-L163 (1975).

64. Yao, Y., Tse, J., Tanaka, K., Marsiglio, F. & Ma, Y. Superconductivity in lithium under high pressure investigated with density functional and Eliashberg theory. Phys. Rev. B 79, 054524 (2009).

65. Garland, J., Bennemann, K. & Mueller, F. Effect of Lattice Disorder on the Superconducting Transition Temperature. Phys. Rev. Lett. 21, 1315-1319 (1968).

Acknowledgments

The authors appreciatively acknowledge the financial support by Natural Science Foundation of China under Grants No. 11474128, 11274136, 11025418, and 91022029, China 973 Program under Grant No. 2011CB808200, and the 2013 Program for New Century Excellent Talents in University. T.I. was supported by MEXT of Japan (Grant Nos. 20103001-20103005) and RIKEN iTHES project. A part of calculations were performed in the computing facilities at RIKEN Integrated Cluster of Clusters system (Japan) and the High Performance Computing Center of Jilin University. The work at the University of Saskatchewan was supported by Natural Sciences and Engineering ResearchCouncil of Canada (NSERC).

Author Contributions

W.C., J.S. and H.L. performed the calculation. W.C., H.L., Y.Y., H.W., T.I. and Y.M. analyzed the data. Y.Y., W.C. and H.W. wrote the manuscript. H.W. and Y.M. supervised the work.

Additional Information

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Cui, W. et al. Hydrogen segregation and its roles in structural stability and metallization: silane under pressure. Sci. Rep. 5, 13039; doi: 10.1038/srep13039 (2015).

i/S 0 I This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/