Pressure-induced phase transitions of lead iodide

Lihua Yanga, Yu Zhang*a, Jun Wangb, Yiding Wanga and William W. Yu*ac
aState Key Laboratory on Integrated Optoelectronics, College of Electronic Science and Engineering, Jilin University, Changchun 130012, China. E-mail: yuzhang@jlu.edu.cn; wyu6000@gmail.com
bState Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, China
cDepartment of Chemistry and Physics, Louisiana State University, Shreveport, LA 71115, USA

Received 26th June 2016 , Accepted 31st August 2016

First published on 1st September 2016


Pressure is a fundamental thermodynamic variable that can efficiently modify crystal or electronic structure leading to new material states of interest. Here we search the crystal structures of lead iodide (PbI2) up to 200 GPa using the swarm-intelligent CALYPSO structure prediction method combined with ab initio calculations. Four new stable high-pressure phases (orthorhombic Pnma, tegragonal I4/mmm, orthorhombic Immm and monoclinic C2/c) have been uncovered. A layered to three-dimensional structure transition from the ambient-pressure 2H phase to the orthorhombic Pnma phase is proposed to be taking place at ∼2.5 GPa. The newly proposed Pnma phase is a semiconductor with an indirect band gap of 1.68 eV. It is isosymmetric to the well-known cotunnite phase but with a very different bonding pattern. The semiconductor–metal transition of PbI2 occurs at ∼27.2 GPa, where the semiconducting Pnma phase transforms to the metallic I4/mmm phase due to pressure-induced broadening of both the valance and conduction bands. The metallic nature persists in the higher-pressure phases of Immm and C2/c up to about 169.1 GPa, above which PbI2 decomposes to its constituent elements. The current results represent a significant step toward the understanding of structural and electronic properties of PbI2 under compression.


Introduction

Lead iodide (PbI2) has received considerable research effort for more than several decades because of its important scientific and technological applications.1 Being a typical layered semiconductor with a direct bandgap around 2.3–2.5 eV,2 it is a prototypical system for investigating the fundamental properties of excitions3,4 and has been traditionally used as a stable nuclear radiation detector,5–7 which converts γ- and X-ray radiation directly into current carriers. Recently, inspired by the emergence and rapid growth of two-dimensional (2D) materials (such as graphene,8,9 transition-metal dichalcogenides,10,11 and phosphorene12,13), single layered and ultrathin 2D PbI2 crystals have also been intensively explored by both experiment and theory as potential complementary compounds of existing 2D materials, which show promising applications in nanoelectronics and optoelectronics.14,15 In addition, PbI2 can serve as a precursor for low-cost synthesis of hybrid organic-inorganic perovskites, which are recent super star materials for photovoltaics.16,17 Due to the versatile applications, various forms of PbI2 including nanoclusters, nanotubes and thin film have been intensively explored.18–23

At ambient condition, PbI2 adopts a layered crystalline structure with repeating units of [I–Pb–I] stacking along the z direction. Each [I–Pb–I] layer consists of a hexagonally closed-packed layer of lead ions sandwiched between two layers of iodine ions [Fig. 1(a) and (b)], which are tightly hold together by ionic or covalent bonding.24 Adjacent layers are bound via the relatively weak van der Waals forces typical of many two dimensional materials (e.g., graphite). Different stacking pattern of the [I–Pb–I] layers lead to a series of polytypes, such as 2H, 4H, 6H, 6R, 12R, etc, among which the 2H phase with P[3 with combining macron]m1 symmetry (denote 2H-P[3 with combining macron]m1 below) is the most stable form.24


image file: c6ra16487k-f1.tif
Fig. 1 Crystal structures of various PbI2 phases, (a) and (b) 2H-P[3 with combining macron]m1, (c) Pnma, (d) I4/mmm, (e) Immm and (f) C2/c.

It is known that pressure is a fundamental thermodynamic variable that can efficiently alter the binding pattern of compounds and thus leads to new phase with intriguing chemical and physical properties. Studies concerning behaviors of PbI2 under compression have long been carried out since 1937 and 1948, where two phase transitions from 2H-P[3 with combining macron]m1 structure were observed at 0.4 GPa and 2.5 GPa by Bridgman as characterized by volume changes of 4% and 2%, respectively.25,26 Later measurements on the optical absorption coefficient have also observed the first phase transition at 0.4 GPa.27 Jayaraman et. al. performed Raman and optical absorbtion experiments, where the two phase transition was observed at 0.5 and 3 GPa, respectively.28 They deduced that the first phase is likely to possess a three-dimensional (3D) structure. In 2000, Saitoh et. al. investigated the pressure dependence of the Raman and luminescence spectra of 2H-P[3 with combining macron]m1 structure up to 1.5 GPa.29 Based on the systematical analysis of experimental data, they found that the first phase transition occurs at 0.56 GPa and the new phase is a 2 × 2 supercell of 2H-P[3 with combining macron]m1 with relative atom position and bondings hardly changed. Though pressure-induced phase transitions of PbI2 have long been found for several decades, the detailed structural information of these high-pressure phases are still elusive, impeding further understanding of the physical and chemical properties. Therefore, it is highly desired to determine the crystal structure of PbI2 at high pressure.

Here, we have extensively explored the structures of PbI2 at a wide pressure range of 0–200 GPa through elaborative swarm structure searching30,31 combined with first-principles calculations. Four high pressure phases, orthorhombic Pnma, tetragonal I4/mmm, orthorhombic Immm and monoclinic C2/c, were uncovered to be stable at pressure intervals of 2.5–27.2, 27.2–74.5, 74.5–92.5 and 92.5–169.1 GPa, respectively. The predicted orthorhombic Pnma structure is a semiconductor with an intriguing 3D structure, which is isosymmetry to the well-known cotunnite structure for AB2 compounds but with very different bonding patterns. The other three high pressure phases also adopt 3D structural patterns but show metallic nature. Above 169.1 GPa, PbI2 is predicted to be unstable and dissociate into the constituent elements (Pb + 2I).

Computational method

Our structure searching simulations were performed through the swarm-intelligence CALYPSO method,30,31 which enables a global minimization of free-energy (reduced to enthalpy at 0 K) surface merging first-principles total-energy calculations. The method is specially designed for global search of ground-state structures at given chemical composition and external condition (e.g. pressure) unbiased by any known structural information, and has been successfully applied to predictions of high-pressure phases of various system, such as Li,32 Bi2Te3,33 Xe–Fe,34 B3NO35 and mercury chalcogenides,36 etc. The underlying ab initio structural relaxations and electronic band structure calculations were performed in the framework of density functional theory (DFT) using plane-wave pseudopotential method as implemented in the VASP code.37 The electron–ion interaction is described by means of projector-augmented wave method.38 Configurations of 6s25d106p2 for Pb, 5s25p5 for I are considered as valence electrons, respectively. We use the Perdew–Burke–Ernzerhof (PBE) parameterization of generalized gradient approximation as the exchange-correlation functional.39 An energy cutoff of 400 eV for the plane-wave expansion and Monkhorst–Pack k-point meshes with spacing of 2π × 0.03 Å − 1 was used to sample the Brillouin zone, which is shown to yield excellent convergence for total energies. It is known that DFT-PBE approach seriously underestimates the band gaps of all the materials due to the self-interaction error. In order to given an accurate description of the electronic structures, the hybrid HSE06 functional40 (with the standard 25% non-local Fock exchange) is used for band structure and density of states calculations. Furthermore, relativistic effect of spin–orbit coupling (SOC), which is especially important to high-Z element of Pb, is also included for band structure and density of states calculations. The phonon calculations were carried out by using a supercell approach as implemented in the PHONOPY code.41

Results and discussions

Structure predictions through the CALYPSO code were performed with 1–4 formula unit (f.u.) per simulation cell at 0, 50, 80, 100, 150, 200 GPa and 1–8 f.u. at 10 GPa. The detailed structure information and calculated enthalpies of low-lying structures from the structure searches are given in Table SI–IV in the ESI. It is seen that the energy differences between these low-lying structures is usually several or several tens meV, which is substantial enough to identify the lowest-energy structures. At 0 GPa, our simulations readily reproduced the experimental 2H-P[3 with combining macron]m1 structure as the global minimum structure for PbI2, validating the method adopted here. At higher pressures, we uncovered a group of new structures. Among them, four structures were found to be most stable at certain pressure range. At 10 GPa, simulations revealed an orthorhombic Pnma structure with 4 f.u. per unit cell (Z = 4). As depicted in Fig. 1(c), this structure is 3-dimensional, which is in stark contrast to that of layered 2H-P[3 with combining macron]m1 [Fig. 1(a) and (b)]. It contains one crystallographically distinct Pb atom and two I atoms, all of which locate on the fourfold 4c site. The Pb atom is coordinated to eight I atoms, while each I atom is coordinated to four Pb atoms. It should be noteworthy that the current orthorhombic Pnma phase is different from the isosymmetric cotunnite phase (α-PbCl2 phase), which is commonly adopted for many AB2 compounds.41 As seen in Fig. 2, the structural topology of the Pb–I network in the two phases are essentially different, and the cations are 9-coordinated in the well-known cotunnite phase.42 It should be further notice that another Pnma structure has also been found for the high-pressure phases of FeF2 and CoF2 previously.43,44 However, this structure is different from the current one, since it is about 415 meV higher in energy the current one for PbI2 at 10 GPa (Table SI in the ESI). To the best of our knowledge, the current orthorhombic Pnma phase has not been observed in AB2 compounds previously. At 50 GPa, a tetragonal I4/mmm structure with Z = 2 was uncovered. As shown in Fig. 1(d), this structure is also a 3-dimensional structure composed of one crystallographically distinct Pb atom occupying 2b site and one I atom occupying 4e site. Within this structure, each Pb atom is coordinated to ten I atoms forming a double-capped square prism and each I atom is five-coordinated. The distance between central Pb atom and I atom at the vertex of the prism is 3.015 Å, and the distance between central Pb atom and capped I atoms is 3.311 Å. The double-capped square prisms stacked along the c axis and share the triangular facets each other. The structure uncovered at 80 GPa through the CALYPSO searches adopts an orthorhombic Immm symmetry with Z = 2 [Fig. 1(e)], which can be derived from symmetry lowing of the lower-pressure I4/mmm structure. As a result, on pair of facets of the prism formed by I atoms change from regular square to rectangle, but the coordination number of the Pb and I atoms are same as those in I4/mmm structure. Upon increasing pressure to 100 and 150 GPa, a monoclinic C2/c structure with Z = 4 is found as show in Fig. 1(f). This structure contains one distinct Pb atom and two I atoms, all of which occupy the 4e site. Within this structure, Pb atoms is still 10-coordinated, but the coordination number of the two distinct I atoms increase to 5 and 7, respectively.
image file: c6ra16487k-f2.tif
Fig. 2 Frameworks of (a) current orthorhombic Pnma structure and (b) the cotunnite structure (α-PbCl2 structure).

The calculated enthalpies of various phases of PbI2 relative to Pb + 2I are plotted as a function of pressure in Fig. 3(a), which confirm the energetic stability of our predicted structures. The experimental 2H-P[3 with combining macron]m1 structure is most stable up to around 2.5 GPa, beyond which orthorhombic Pnma structure becomes energetically most favorable. In the pressure range 27.2–74.5 GPa, the tetragonal I4/mmm structure has the lowest enthalpies, while orthorhombic Immm structure is stable between 74.5 GPa and 92.5 GPa, the monoclinic C2/c structure is found to be stable up to 169.1 GPa, above which PbI2 becomes unstable against decomposition to constituent elements (Pb + 2I). Our phonon calculations have verified the dynamical stabilities of all these predicted structures by evidence of the absence of any imaginary frequency in the whole Brillouin zone [Fig. 4]. Thus, the current results suggest PbI2 takes the phase transition sequence of 2H-P[3 with combining macron]m1 → PnmaI4/mmmImmmC2/c under compression.


image file: c6ra16487k-f3.tif
Fig. 3 (a) Calculated enthalpies per formula unit as functions of pressure with respect to Pb + 2I and (b) calculated equation of states of various phases of PbI2.

image file: c6ra16487k-f4.tif
Fig. 4 Calculated phonon-dispersion curves for (a) Pnma, (b) I4/mmm, (c) Immm and (d) C2/c phases of PbI2. The 2 × 2 × 2, 4 × 4 × 4, 4 × 4 × 4 and 3 × 3 × 3 supercells were used for the calculations.

It is known that experiments have observed two phase transitions of PbI2, which occur at ∼0.5 and ∼3 GPa as characterized by volume changes of 4% and 2%, respectively.25–29 The current results indicate one of the high-pressure phases adopts the 3D orthorhombic Pnma structure stable above 2.5 GPa. This phase is most likely corresponding to the second one observed in experiments as evidenced by the calculated equation of states of PbI2 [Fig. 3(b)], where the volume decrease of 5.6% from the 2H-P[3 with combining macron]m1 to Pnma structure at 2.5 GPa is in good agreement with that of 6% (2% + 4%) observed by experiments.25,26 The stability field of the first high-pressure phase that do not found in the current work covers only a very narrow pressure range of about 1–2 GPa. In fact, there are some discrepancies between experiments. While Raman and optical measurements in 1986 deduce the first high-pressure phase possesses a 3D structure,28 later Raman and luminescence spectra study in 2000 suggests this phase is a 2 × 2 supercell of 2H-P[3 with combining macron]m1 structure with relative atom position and bondings hardly changed.29 Nevertheless, the first high-pressure phase should be an intermediate phase account for the layered to 3D structural transition. Thus, the complete phase transition sequence of PbI2 under compression should be 2H-P[3 with combining macron]m1 → intermediate phase → PnmaI4/mmmImmmC2/c. Determination of the detailed crystal structure of such an intermediate phase in a narrow pressure range is challenging. Intensively theoretical structure search simulations with high level of accuracy or in situ X-ray diffraction experiments are needed.

To gain insight into the electronic properties of PbI2 under compression, the electronic band structures and partial density of states (PDOS) of various phases are calculated and shown in Fig. 5. For 2H-P[3 with combining macron]m1 structure, the calculated HSE + SOC band structure shows a direct bandgap of 2.15 eV at the A point of the Brillouin zone [Fig. 5(a)], which is in reasonable agreement with the experimental value of 2.3–2.5 eV.2 For the high-Z element Pb, the Pb-6s state moves down in energy (forming lone-pair states), separating from Pb-6p state (mass-Darwin relativistic effect).45 For +2 low oxidation state of Pb in PbI2, there exists charge transfer form Pb-6p to I-5p. Thus the valence band (VB) edge is mainly formed by anti-bonding hybridization of the Pb-6s and I-5p states, while the conduction band (CB) edge is derived from anti-bonding hybridization of the Pb-6p and I-5p states. This situation can be clear seen form the PDOS as show in Fig. 5(f). The existence of Pb-6s lone-pair states renders the VB edges broaden and dispersive, leading to a small hole effective mass at the VB maximum. This is in contrast to the oxides having a rather narrow VB due to localized O-2p states, or the transition metal compounds where the narrow VB reflects localized d states, both of which show large hole effective masses. This type of physics has been elucidated earlier for binary lead chalcogenides.46 For the orthorhombic Pnma structure, the band structure [Fig. 5(b)] calculated at 10 GPa shows an indirect bandgap of 1.68 eV formed by the VB maximum at U point and CB minimum at G point of the Brillouin zone. Through the structural topology is dramatically changed, the bonding character is essentially same as that of the 2H-P[3 with combining macron]m1 structure, e.g. VB edge formed by anti-bonding hybridization of Pb-6s and I-5p states, and CB edge formed by anti-bonding hybridization of Pb-6p and I-5p states [Fig. 5(g)]. PbI2 is still a semiconductor at the stable field of Pnma phase (2.5–27.2 GPa). Band structure of the tetragonal I4/mmm structure [Fig. 5(c)] calculated at 50 GPa shows that both the VB and CB are broadened and pass through the Fermi level, indicating electron delocalization and pressure-induced metallization occurs. In higher pressure phases, Immm and C2/c, the metallic nature of PbI2 persists [Fig. 5(d) and (e)] up to the decomposition pressure of ∼169.1 GPa.


image file: c6ra16487k-f5.tif
Fig. 5 Electronic band structures (top row) and partial density of states (bottom row) for 2H-P[3 with combining macron]m1 [(a) and (f)], Pnma [(b) and (g)], I4/mmm [(c) and (h)], Immm [(d) and (i)] and C2/c [(e) and (j)] phases of PbI2, respectively.

Conclusions

In summary, through systematic swarm-intelligent CALYSPO structure searches combined with first-principles calculations, we report the phase transition sequence of 2H-P[3 with combining macron]m1 → intermediate phase → PnmaI4/mmmImmmC2/c for PbI2 under compression, in which four new stable high-pressure phases have been uncovered. Within this phase transition sequence, layered to 3D structure transition was proposed to be taking place at 2.5 GPa from the 2H-P[3 with combining macron]m1 structure to an intriguing orthorhombic Pnma structure that have not been seen in other AB2 compounds. Pressure-induced metallization is observed at the phase transition from Pnma to I4/mmm (27.1 GPa). The metallicity of PbI2 persists in higher-pressure phases of Immm and C2/c up to about 169.1 GPa, above which PbI2 decomposes to constituent elements. The current work presents significant results of fundamental structural and electronic properties of PbI2 under compression and will inevitably stimulate future experimental and theoretical study of PbI2 at high pressure.

Acknowledgements

The authors appreciate the financial supports form the National Natural Science Foundation of China (51272084, 61225018, 61475062), the Jilin Province Key Fund (20140204079GX) and the Postdoctoral Science Foundation of China (2015M571360).

Notes and references

  1. M. Tubbs, Phys. Status Solidi, 1972, 11, 11–50 CrossRef .
  2. R. Ahuja, H. Arwin, A. F. Da Silva, C. Persson, J. M. Osorio-Guillén, J. Souza De Almeida, C. M. Araujo, E. Veje, N. Veissid, C. Y. An, I. Pepe and B. Johansson, J. Appl. Phys., 2002, 92, 7219–7224 CrossRef CAS .
  3. M. Watanabe and T. Hayashi, J. Phys. Soc. Jpn., 1994, 63, 785–794 CrossRef CAS .
  4. B. Dorner, R. E. Ghosh and G. Harbeke, Phys. Status Solidi, 1976, 73, 655–659 CrossRef CAS .
  5. T. Yanagida, Y. Fujimoto, A. Yoshikawa, Y. Yokota, K. Kamada, J. Pejchal, V. Chani, N. Kawaguchi, K. Fukuda, K. Uchiyama, K. Mori, K. Kitano and M. Nikl, Appl. Phys. Express, 2010, 3, 056202 CrossRef .
  6. R. A. Street, S. E. Ready, K. Van Schuylenbergh, J. Ho, J. B. Boyce, P. Nylen, K. Shah, L. Melekhov and H. Hermon, J. Appl. Phys., 2002, 91, 3345–3355 CrossRef CAS .
  7. M. Nikl, Meas. Sci. Technol., 2006, 17, R37–R54 CrossRef CAS .
  8. K. S. Novoselov, Science, 2004, 306, 666–669 CrossRef CAS PubMed .
  9. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V Grigorieva, S. V. Dubonos and A. A. Firsov, Nature, 2005, 438, 197–200 CrossRef CAS PubMed .
  10. J. N. Coleman, M. Lotya, A. O'Neill, S. D. Bergin, P. J. King, U. Khan, K. Young, A. Gaucher, S. De, R. J. Smith, I. V. Shvets, S. K. Arora, G. Stanton, H.-Y. Kim, K. Lee, G. T. Kim, G. S. Duesberg, T. Hallam, J. J. Boland, J. J. Wang, J. F. Donegan, J. C. Grunlan, G. Moriarty, A. Shmeliov, R. J. Nicholls, J. M. Perkins, E. M. Grieveson, K. Theuwissen, D. W. McComb, P. D. Nellist and V. Nicolosi, Science, 2011, 331, 568–571 CrossRef CAS PubMed .
  11. Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman and M. S. Strano, Nat. Nanotechnol., 2012, 7, 699–712 CrossRef CAS PubMed .
  12. S. Zhang, M. Xie, F. Li, Z. Yan, Y. Li, E. Kan, W. Liu, Z. Chen and H. Zeng, Angew. Chem., Int. Ed., 2015, 55, 1666–1669 CrossRef PubMed .
  13. L. Li, Y. Yu, G. J. Ye, Q. Ge, X. Ou, H. Wu, D. Feng, X. H. Chen and Y. Zhang, Nat. Nanotechnol., 2014, 9, 372 CrossRef CAS PubMed .
  14. M. Zhou, W. Duan, Y. Chen and A. Du, Nanoscale, 2015, 7, 15168–15174 RSC .
  15. P. Wangyang, H. Sun, X. Zhu, D. Yang and X. Gao, Mater. Lett., 2016, 168, 68–71 CrossRef CAS .
  16. A. Kojima, K. Teshima, Y. Shirai and T. Miyasaka, J. Am. Chem. Soc., 2009, 131, 6050–6051 CrossRef CAS PubMed .
  17. T. C. Sum and N. Mathews, Energy Environ. Sci., 2014, 7, 2518–2534 CAS .
  18. G. K. Kasi, N. R. Dollahon and T. S. Ahmadi, J. Phys. D: Appl. Phys., 2007, 40, 1778 CrossRef CAS .
  19. R. Kreizman, S. Y. Hong, J. Sloan, R. Popovitz-Biro, A. Albu-Yaron, G. Tobias, B. Ballesteros, B. G. Davis, M. L. H. Green and R. Tenne, Angew. Chem., Int. Ed., 2009, 48, 1230–1233 CrossRef CAS PubMed .
  20. R. Kreizman, A. N. Enyashin, F. L. Deepak, A. Albu-Yaron, R. Popovitz-Biro, G. Seifert and R. Tenne, Adv. Funct. Mater., 2010, 20, 2459–2468 CrossRef CAS .
  21. Y. Liu, C. Hu, A. Comotti and M. D. Ward, Science, 2011, 333, 436–440 CrossRef CAS PubMed .
  22. X. Liu, S. T. Ha, Q. Zhang, M. de la Mata, C. Magen, J. Arbiol, T. C. Sum and Q. Xiong, ACS Nano, 2015, 9, 687–695 CrossRef CAS PubMed .
  23. Y. Wang, Y.-Y. Sun, S. Zhang, T.-M. Lu and J. Shi, Appl. Phys. Lett., 2016, 108, 013105 CrossRef .
  24. P. A. Beckmann, Cryst. Res. Technol., 2010, 45, 455–460 CrossRef CAS .
  25. P. W. Bridgman, Proc. Am. Acad. Arts Sci., 1937, 72, 45 CrossRef CAS .
  26. P. W. Bridgman, Proc. Am. Acad. Arts Sci., 1948, 76, 71 CrossRef CAS .
  27. C. Carillon and G. Martinez, Nuovo Cimento B, 1977, 38, 496–502 CrossRef .
  28. A. Jayaraman, R. G. Maines and T. Chattopadhyay, Pramana, 1986, 27, 449–457 CrossRef CAS .
  29. A. Saitoh, T. Komatsu and T. Karasawa, Phys. Rev. B: Condens. Matter Mater. Phys., 2000, 62, 11398–11404 CrossRef CAS .
  30. Y. Wang, J. Lv, L. Zhu and Y. Ma, Phys. Rev. B: Condens. Matter Mater. Phys., 2010, 82, 94116 CrossRef .
  31. Y. Wang, J. Lv, L. Zhu and Y. Ma, Comput. Phys. Commun., 2012, 183, 2063–2070 CrossRef CAS .
  32. J. Lv, Y. Wang, L. Zhu and Y. Ma, Phys. Rev. Lett., 2011, 106, 15503 CrossRef PubMed .
  33. L. Zhu, H. Wang, Y. Wang, J. Lv, Y. Ma, Q. Cui, Y. Ma and G. Zou, Phys. Rev. Lett., 2011, 106, 145501 CrossRef PubMed .
  34. L. Zhu, H. Liu, C. J. Pickard, G. Zou and Y. Ma, Nat. Chem., 2014, 6, 644 CAS .
  35. Q. Li, J. Wang, M. Zhang, Q. Li and Y. Ma, RSC Adv., 2015, 5, 35882–35887 RSC .
  36. Y. Yan, S. Zhang, Y. Wang, G. Yang and Y. Ma, RSC Adv., 2015, 5, 104426–104432 RSC .
  37. G. Kresse and J. Furthmüller, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 11169–11186 CrossRef CAS .
  38. P. E. Blöchl, Phys. Rev. B: Condens. Matter Mater. Phys., 1994, 50, 17953–17979 CrossRef .
  39. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868 CrossRef CAS PubMed .
  40. A. V. Krukau, O. A. Vydrov, A. F. Izmaylov and G. E. Scuseria, J. Chem. Phys., 2006, 125, 224106 CrossRef PubMed .
  41. A. Togo, F. Oba and I. Tanaka, Phys. Rev. B: Condens. Matter Mater. Phys., 2008, 78, 134106 CrossRef .
  42. J. Haines, J. M. Léger and O. Schulte, Phys. Rev. B: Condens. Matter Mater. Phys., 1998, 57, 7551–7555 CrossRef CAS .
  43. J. A. Barreda-Argüeso, S. López-Moreno, M. N. Sanz-Ortiz, F. Aguado, R. Valiente, J. González, F. Rodríguez, A. H. Romero, A. Muñoz, L. Nataf and F. Baudelet, Phys. Rev. B: Condens. Matter Mater. Phys., 2013, 88, 214108 CrossRef .
  44. S. López-Moreno, A. H. Romero, J. Mejía-López, A. Muñoz and I. V. Roshchin, Phys. Rev. B: Condens. Matter Mater. Phys., 2012, 85, 134110 CrossRef .
  45. B. Thaller, The Dirac Equation, Springer, Berlin, Heidelberg, 1992 Search PubMed .
  46. S.-H. Wei and A. Zunger, Phys. Rev. B: Condens. Matter Mater. Phys., 1997, 55, 13605–13610 CrossRef CAS .

Footnote

Electronic supplementary information (ESI) available: The detailed structure information and calculated enthalpies of ground-state and low-lying structures form the structure searches. See DOI: 10.1039/c6ra16487k

This journal is © The Royal Society of Chemistry 2016