Unexpected Trend in Stability of Xe–F Compounds ... - ACS Publications

Oct 24, 2016 - College of Physics and Electronic Information, Luoyang Normal University, Luoyang 471022, P. R. China. § ... Department of Chemistry a...
2 downloads 6 Views 1016KB Size
Subscriber access provided by Northern Illinois University

Letter

Unexpected Trend in Stability of Xe–F Compounds Under Pressure Driven by Xe–Xe Covalent Bonds Feng Peng, Jorge Botana, Yanchao Wang, Yanming Ma, and Mao-Sheng Miao J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b01922 • Publication Date (Web): 24 Oct 2016 Downloaded from http://pubs.acs.org on October 25, 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 Letters 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 7

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 Letters

Unexpected Trend in Stability of Xe–F Compounds Under Pressure Driven by Xe–Xe Covalent Bonds Feng Peng,†‡§ Jorge Botana,§ Yanchao Wang,† Yanming Ma,†* and Maosheng Miao#§* †



State Key Lab of Superhard Materials, Jilin University, Changchun 130012, P. R. China College of Physics and Electronic Information, Luoyang Normal University, Luoyang 471022, P. R. China § Beijing Computational Science Research Center, Beijing 10084, P. R. China # Department of Chemistry and Biochemistry, California State University Northridge, CA 91220, USA * Corresponding author. Email: [email protected] and [email protected] Abstract: Xenon difluoride is the first and the most stable of hundreds of noble-gas (Ng) compounds. These compounds reveal the rich chemistry of Ngs. No stable compound that contains a Ng–Ng bond has been reported previously. Recent experiments have shown intriguing behaviors of this exemplar compound under high pressure, including increased coordination numbers and an insulator-to-metal transition. None of the behaviors can be explained by electronic-structure calculations with fixed stoichiometry. We therefore conducted a structure search of xenon–fluorine compounds with various stoichiometries, and studied their stabilities under pressure using first principles calculations. Our results revealed, unexpectedly, that pressure stabilizes xenon–fluorine compounds selectively, including xenon tetrafluoride, xenon hexafluoride, and the xenon-rich compound Xe2F. Xenon difluoride becomes unstable above 81 GPa, and yields metallic products. These compounds contain xenon–xenon covalent bonds, and may form intercalated graphitic xenon lattices, which stabilize xenon-rich compounds and promote the decomposition of xenon difluoride.

For Table of Contents Only

Predicted phase diagram and stable phases of Xe-F compounds.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

(vdW) interactions were examined using the PBE and vdWDF2 density functional and found to have only minor effects (in the order of meV/atom) on the energetics of XeFx at high pressures (> 50 GPa), and the vdW corrections are negligible under such conditions. To ensure the validity of the pseudopotentials used here, we also performed full-potential allelectron calculations for the equation of states of the XeF2 compound using the WIEN2K code22. Our all-electron calculations yield almost identical results to the Vienna Ab-initio Simulation Package (Figure S1), which validates our pseudopotentials.

Since the first discovery of noble gases (Ngs) by Ramsey in 19041, the synthesis of compounds that contain them has presented a significant challenge to chemists. These elements are inert to chemical reactions because their outermost shells are filled. The reaction of xenon (Xe) with fluorine (F) was first predicted by Pauling, and was realized in the laboratory, to form xenon difluoride (XeF2), by Bartlett and Hoppe2 some 30 years later. Since then, a few hundred Ng compounds and species have been obtained. In these compounds, Ng elements are oxidized and share their outer-shell electrons to form covalent or ionic bonds. Among them, XeF2 is the most stable and the best studied. Under ambient conditions, it is a crystalline insulator consisting of linear XeF2 molecules. Pressure may affect the atomic and electronic structure of matter. Early research suggested that the low-pressure structure of XeF2 remained stable to 50 GPa1,2. A recent experiment revealed that XeF2 may transform into two- and three-dimensional extended solids and become metallic when a higher pressure is applied3,4. However, such transitions cannot be explained by using recent density-functional theory (DFT) calculations and structure searches of XeF2, which show that the ambient-pressure I4/mmm structure persists as the thermodynamically most stable phase of XeF2 up to 105 GPa5. Above 105 GPa, the I4/mmm structure transforms to a Pnma structure that contains bent F–Xe–F molecules. Two Ng atoms that both have a filled valence shells do not engage in chemical bonding. However, once electrons have been depleted from the outer shell, they may be capable of forming chemical bonds. Drews et al.6 isolated a salt of cationic Xe2+ containing a Xe–Xe bond with an interatomic distance of 3.087(1) Å. Stein et al.7 also identified Xe2+ experimentally. Seidel et al.8 detected and characterized by calculations a blue Xe4+ cation with a linear, symmetric structure. Fernández and Frenking9 proposed the viability of a metastable HXeXeX (X = F to I) molecule with an energy barrier of 13 kcal mol−1. They predicted that it could be synthesized and studied in a low-temperature Xe matrix.9,10 Recently, Somayazulu et al. studied the Xe–H2 binary system under high pressure.11 They found that the nearest neighbor Xe–Xe distance at 4.9 GPa is close to that of neutral dimers in the gas phase (3.84 Å)12. In this work, we studied the structures and stabilities of selected stoichiometric XeFx (x = 1/4, 1/3, 2/3, and 1–8) compounds using an accurate first-principles method and an automatic structure search based on a particle swarm optimization algorithm implemented in the CALYPSO method13,14. The crystal structure predictions were performed under various pressures, including 0, 50 100, and 200 GPa. The CALYPSO method has been benchmarked in many well-studied systems, including elements and binary and ternary compounds15,16. Total-energy calculations were performed in the DFT framework within the Perdew–Burke–Ernzerhof17 parameterization of generalized gradient approximation18 as implemented in the Vienna Ab-initio Simulation Package code19. The projector-augmented wave method20 was used to describe the ionic potential. 4d105s25p6 and 2s22p5 were treated as valence electrons in the Xe and F atoms, respectively. A plane-wave kinetic energy cutoff of 1200 eV and dense k-point samplings gave excellent convergence of total energies. Phonon calculations were performed for all relevant structures using the phonopy code21. The van der Waals

Figure 1. Phase stabilities of Xe–F compounds. (a) Enthalpies of formation of XeFx under a series of pressures. The dotted lines connect the data points, and the solid lines denote the convex hull. (b) Stable pressure range for XeFx.

Our calculations reveal that Xe–F compounds with unusual stochiometries, such as XeF and Xe2F, can be stabilized under elevated pressure and temperature. We also found a number of new structures that are governed by novel chemical bonding and include the Xe–Xe covalent bonds and the intercalated Xe honeycomb lattices. These structures are dynamically stable in their respective pressure ranges by phonon calculations (Supporting information). We investigated the stabilities of the Xe–F compounds by calculating and comparing the formation enthalpies of each compound in the structure with lowest enthalpy from 0 to 200 GPa. The most efficient way to compare the stabilities of the binary compounds with variable stoichiometry is to construct a convex hull (Fig. 1a). The compounds that are located on the hull are stable; whereas those that are located above the hull are unstable, and will decompose to form the compounds on the hull23-25.

2 ACS Paragon Plus Environment

Page 2 of 7

covalent bond in Xe2F. Bader analysis28 reveals a charge transfer of 0.37 e from Xe to F. To characterize the bonding nature including Xe-Xe and Xe-F bonds, we calculated the Crystalline Orbital Hamiltonian Population (COHP) and the integrated COHP (ICOHP) for Xe-Xe and Xe-F pairs Xe2F at 200 GPa (Figure 3a). The results clearly revealed the bonding nature of neighboring Xe-Xe pairs in Xe2F. It shows that Xe-Xe pairs in Xe2F are form its covalent bonds. The COHP plots (Figure 3a) also reveal fully occupied bonding states and partially occupied anti-bonding states, leading to net covalent bonding between the two neighboring Xe atoms. We also calculated the electron localization function (ELF) 29-32 (Fig. 4a). As a rule, a typical covalent bond has an ELF value larger than 0.7. However some covalent bonds have lower ELF value (even less than 0.5)33-37. Indeed, for a covalent bonding molecule I2, the ELF value is about 0.6 as calculated by the same method (Fig. 4c). The large ELF value of ~0.5 between neighboring Xe atoms confirms the presence of covalent bonds. This is the first example of compounds that possess covalent bonds between Ng atoms. Furthermore, the nearest Xe–F (2.35 Å) and F–F (2.28 Å) distances are both much larger than the covalent Xe–F and F–F bond lengths of 2.00 Å in the XeF2 molecule38 and 1.56 Å in F2 respectively39. This suggests that Xe and F or F and F do not bond covalently in Xe2F. The calculated COHP (Figure 3b) shows that the Xe-F bonds between the nearest neighboring Xe-F pairs have smaller COHP values than Xe-Xe (Figures 3a and 3b), indicating weaker covalent bonding. It confirms the ionic character of F atoms, as shown by the Bader charge analysis and the high coordination number 8 of F atoms. In addition, the spherical ELF distribution (Fig. 4b) also confirms the states with closed-shell electron configurations of F atoms. The projected density of states (PDOS) (Fig. 4e) reveals that Xe2F is metallic. The Xe 5p and F 2p orbitals contribute substantially to the DOS at the Fermi level. a

Figure 2. Most stable structures of various Xe-F compounds at selected pressures. (a) Xe2F in an I4/mcm structure at 200 GPa; (b) Xe3F2 in a P21/c structure at 50 GPa; (c) XeF4 in an I4/m structure at 200 GPa; (d, e, f) XeF6 in C2/c structure at 50 GPa, Cmcm structure at 100 GPa, and R-3 structure at 200 GPa, respectively. The smaller green spheres are F atoms and the larger purple ones are Xe atoms.

0.5

b

0.0

E (Ry)

As shown by the convex hull, XeF2, xenon tetrafluoride (XeF4), and xenon hexafluoride (XeF6) are the only stable compounds at ambient pressure. These results compare well with experimental observations and general knowledge on Xe chemistry: Xe is stable in the +2, +4, and +6 states at ambient pressure. Although other Xe–F compounds do not locate on the convex hull at 0 GPa, the energy difference is small and is of the order of a few meV. By taking xenon trifluoride (Xe2F6) as an example, we calculate the relative enthalpies of formation at 0 GPa, including the van der Waals interaction. The results (Figure S2) show that the inclusion of this interaction can stabilize Xe2F6. The relative enthalpies of formation per atom for XeF2, Xe2F6, and XeF4 are almost identical (less than 2 meV difference), which indicates that all three compounds may become stable at 0 GPa. Pressure can change the stability of the Xe–F compounds significantly. Figure 1 shows that Xe2F and Xe3F2 become stable at 60 GPa and 25 GPa, respectively. Surprisingly, pressure destabilizes XeF2, which is the first synthesized and the best-known Ng compound26. Above 81 GPa, XeF2 decomposes to XeF4 + Xe2F. Our calculations show that Xe3F, XeF, xenon pentafluoride, xenon heptafluoride (XeF7), and xenon octafluoride (XeF8) are unstable throughout the studied pressure range (0–200 GPa).

E (Ry)

-0.5

0.00 0.75

-0.075 0.000

ICOHP

0.5

Xe–F compounds exhibit interesting structural features and electronic properties. Figure 1b shows that Xe2F is stable in the structure with I4/mcm symmetry above 60 GPa. In this structure, Xe atoms form intercalated graphitic layers (honeycomb lattice) in two perpendicular directions. The distance between neighboring layers is 3.90 Å at 100 GPa. The intercalated layers form parallel channels, in which chains of F atoms are located. At 200 GPa, all in the graphitic layer Xe-Xe distances have the same value 2.63 Å. This value is close to the summation of the covalent radius of two Xe atoms (1.31 Å)27, and suggests the formation of a Xe–Xe

0.0

-0.5

COHP

c

0 .5

-0.2 0.0 0.2

d

0.0

COHP

-0.05

0.00

ICOHP

0.5

0.0

E(Ry)

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 Letters

E (Ry)

Page 3 of 7

-0.5

-0.5

-1

0

COHP

1

-0.1 0.0

-0.75 0.00 0.75

ICOHP

COHP

3 ACS Paragon Plus Environment

-0.1

0.0

ICOHP

The Journal of Physical Chemistry Letters

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

C2/m structure (Figure 5d) has the lowest XeF energy. In this structure, Xe atoms can be divided into two categories: isolated Xe ions (Xe1) and Xe–Xe dimers (Xe2). At 100 GPa, the shortest Xe-F distance is 2.45 Å, which is much larger than the covalent Xe–F bond length of 2.00 Å in the XeF2 molecule38. Bader charge analysis reveals the large charge transfer from Xe2 to F indicates that XeF consists of Xe+ and F− ions, while Xe1 atoms have almost no charge transfer to F atoms. Different from most other Xe–F compounds, the F atoms do not form covalent bonds with Xe and remain as isolated anions. Xe+(Xe2) ions have seven electrons, including two that occupy 5s and five that occupy 5p orbitals. They are isoelectronic to I atoms. The bond length of Xe–Xe dimers is 2.54 Å at 100 GPa, which is smaller than two times the covalent radius of 1.31 Å for the Xe atoms27. The calculated COHP and the ICOHP for Xe-Xe pairs (Figure 3c) reveal the covalent bonding nature of neighboring Xe-Xe pairs in XeF. The calculated COHP (Figure 3d) shows that the Xe-F bonds between the nearest neighboring Xe-F pairs have smaller COHP values than Xe-Xe (Figure S3a), indicating weaker covalent bonding, which is consistent with the observation that F atoms are strongly ionic and bind with Xe atoms mainly by electrostatic forces, as shown by the Bader charge analysis. The Xe–Xe covalent bonding character of XeF is also revealed by the ELF (Figure 5c) and the charge-density difference (Figure 5d). The PDOS (Figure S5) reveals that the C2/m structure of XeF is metallic at 100 GPa.

Figure 3. Plots of the COHP and ICOHP of Xe2F and XeF. For Xe2F at 200 GPa, (a) Xe-Xe pairs separate by 2.62 Å. (b) Xe-F pairs separate by 2.38 Å. For XeF at 100 GPa, (c) Xe-Xe pairs separate by 2.59 Å. (d) Xe-F pairs separate by 2.13 Å.

Figure 1 shows that Xe3F2 is stable above 25 GPa. It is predicted to adopt a P21/c structure with two formula units per unit cell (Fig. 2b) throughout its stable pressure range (25–55 GPa). Xe3F2 is a typical molecular crystal, is composed of a linear XeF2 molecule (Xe atom named Xe1) and two isolated Xe atoms (Xe2) in a formula unit. At 50 GPa, the shortest Xe1–F distance is 1.98 Å, which is the same as that of XeF2, whereas, the next nearest Xe2–F distance is 2.55 Å, No covalent bond exists between F atoms and the Xe2 atoms. The shortest F–F distance is 2.67 Å, which is much larger than that of the F–F bond length (1.56 Å) in F2,31 and suggests that no covalent bond forms between F atoms. The calculated ELF (Fig. 4d) reveals that the Xe–F bond in the XeF2 molecule is covalent, whereas the interaction between F atoms and isolated Xe2 atoms shows no covalent bonding. Bader analysis reveals a charge transfer of 0.6 e from Xe1 to F atoms, and no charge transfer from Xe2 to other atoms. The PDOS of Xe3F2 (Fig. 4f) reveals it as a semiconductor, with a gap of 0.72 eV. This is a natural result of the fact that Xe3F2 consists of XeF2 molecules and Xe atoms, both of which have filled highest occupied molecular or atomic orbitals.

Figure 4. Electronic properties of Xe2F and Xe3F2. (a, b) Calculated ELF plots in (-1 1 0) and (1 1 0) sections of of Xe2F at 200 GPa (c, d) Calculated ELFs of I2 at 0 GPa and Xe3F2 at 50 GPa, respectively. (e, f) Electronic PDOS per atom of Xe2F at 200 GPa and Xe3F2 at 50 GPa, respectively. The dashed line indicates the Fermi energy.

Figure 5. Predicted structures, synthesis conditions and bonding properties of XeF. (a) Predicted R-3m structure for XeF at 50 GPa. (b) Gibbs free energy of formation of XeF at different temperatures at 100 GPa. (c) ELF plots in (1 0 2.41) section of C2/m structure of XeF within one unit cell at 100 GPa. (d) Calculations of charge density difference of XeF in a C2/m structure at 100 GPa.

XeF remains unstable below 200 GPa at 0 K. It tends to decompose to other Xe–F compounds, such as Xe2F and XeF4. However, the enthalpies of XeF are quite close to the convex hull. Several recent studies of other Xe compounds indicate that elevated temperature has a significant effect and may alter the stability40,41. We therefore performed quasiharmonic free-energy calculations from the phonon spectra computed using the finite-displacement method. As shown in Figure 5b, increasing the temperature suppresses the decomposition of XeF to Xe2F and XeF4. At T = 140 K, XeF becomes stable. A structure search shows that a monoclinic

XeF2 adopts a I4/mmm structure (Figure S6a) throughout its stable pressure range (0–100 GPa). In this structure, linear XeF2 molecules that are oriented along the c axis occupy a body-centered tetragonal lattice. The Xe–F bond lengths of the linear XeF2 units are 1.98 Å at 0 GPa. These results are in good agreement with the experimental results42 and pervi-

4 ACS Paragon Plus Environment

Page 4 of 7

Page 5 of 7

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 Letters

ous theoretical result5. Above 110 GPa, XeF2 is stabilized in a Pnma structure (Figure S6b), which consists of bent F–Xe– F molecules with unequal Xe–F bond lengths. At 150 GPa, a new closely related structure (Cmcm) (Figure S6c) approaches the enthalpy of the Pnma structure. Our results agree well with previous theoretical calculations5. However, the geometry change and the corresponding electronicstructure evolution of XeF2 under pressure cannot explain the recent experimental observations.4 Figure S7 shows the proposed Pnnm-2 and Fmmm structures to explain that the experiments are quite high in enthalpy as calculated by DFT. The Pnnm-2 and Fmmm structures exhibit substantial phonon instabilities (i.e., imaginary phonon frequencies), which suggest dynamic instability5. The resistance measurement of XeF2 under pressure showed that XeF2 became metallic above 70 GPa4. We therefore performed DFT calculations using the Heyd–Scuseria– Ernzerhof hybrid functional43 for the electron-exchange energy to investigate the electronic structure of XeF2 in the I4/mmm structure at 100 GPa. The resultant DOS (Figure 6c) shows a gap, which remains for XeF2, throughout the studied pressure range. This obvious disagreement with experiments and the inconsistency of the atomic structures between DFT and the experiments may be resolved by the trend in stability of the Xe–F compounds under pressure. The convex hull (Figure 1a) shows that XeF2 becomes unstable above 81 GPa, and decomposes to Xe2F and XeF4, of which Xe2F is metallic, and thus might contribute to the significant decrease in resistance in reduction of the XeF2 under pressure.

cal reaction. Although ∆U is constant with pressure, the ∆PV term is reduced substantially more than the constant ∆U value because of the volume reduction that is associated with a higher packing efficiency in the Xe2F and XeF4 mixture. The stability of Xe2F6 is at a critical point under atmospheric conditions. With a slight increase in pressure, Xe2F6 becomes unstable (Figure 1) and decomposes to XeF2 and XeF4. Calculations show that Xe2F6 adopts a monoclinic P21/c structure at atmospheric pressure, which is in good agreement with experimental observations44. Xe2F6 can be viewed as a molecular complex that consists of linear XeF2 molecules and planar XeF4 molecules in a 1:1 ratio. The Xe– F bond lengths in the XeF2 and XeF4 molecules are 1.99 Å and 1.93 Å, respectively. They are both slightly larger than those of pure XeF2 (1.98 Å) and XeF4 (1.91 Å) molecular crystals at atmosphere pressure. XeF4 remains stable throughout the studied pressure range up to 200 GPa. Pressure stabilizes XeF4 considerably, and is the main reason that the other compounds become unstable. A structure search at 0 GPa showed a monoclinic structure with P21/c symmetry (Figure S8), which agrees well with the experimental results45,46. At ~30 GPa, this structure transforms to a tetragonal structure with a higher symmetry of I4/m (Figure 2c), and no further structural transformation occurred to 200 GPa. Both structures consist of XeF4 planar molecules. The Xe–F bond lengths in these two structures are almost identical, at 1.91 Å at 100 GPa. The Xe–F bond length increases with increasing pressure, which results because owing to the stronger inter-molecular interactions weakening the Xe–F covalent bonds. Similar to XeF4, XeF6 is also stable throughout the studied pressure range (from 0 GPa to 200 GPa). Because its structure at ambient pressure has been well studied, we only searched for structures from 50 to 200 GPa. Our results show that the monoclinic structure with C2/c symmetry (Figure 2d) remains stable up to 60 GPa. This structure consists of XeF6 octahedron molecules, which is similar to the ambientpressure phase. The Xe coordination number is 6. The monoclinic C2/c structure transforms to a Cmcm structure (Figure 2e) at 60 GPa, and the structure features chain-like polymerization of ⋯–F–Xe–F–Xe–F⋯, which results in the connection of neighboring XeF6 octahedrons. The Xe coordination number increases to 8. At 120 GPa, XeF6 transforms to a hexagonal structure with R-3 symmetry (Figure 2f). The F atoms bonded with two Xe to form a threedimensional network structure. The Xe coordination number increases from 6 to 8 with a structural change from C2/c to Cmcm, and then to 12 in the R-3 structure. The atomic and electronic structural features indicate that the oxidation state of Xe increases monotonically with the composition ratio of F. To examine this behavior further, we calculated the charges using Bader analysis for XeFx at 100 GPa, as shown in Figure 6b. For XeFx (x ≤ 6) compounds, the Xe charge increases almost linearly with F content. In contrast, for XeFx (x > 6), the Xe charge increases much less with an increasing number of F atoms. This occurs because the 5p electrons are easier to lose than the 5s electrons, and which is also the reason that XeF7 and XeF8 cannot form below 200 GPa. In this pressure range, Xe acts as a p-block element. While sharing its 5p electrons, Xe can bond covalently with F to form XeFx molecules or extended networks. Interestingly, many different bonding features, including inter-molecular, covalent, ionic as well as metallic, can be

Figure 6. Decomposition at high pressure, and electronic properties of XeF2. (a) Relationship between ∆H, ∆U, and ∆PV with pressure for 7XeF2 → 2Xe2F + 3XeF4. We define ∆H = 2H Xe2F + 3HXeF4 – 7HXeF2 , ∆U = 2UXe2F + 3UXeF4 – 7UXeF2, and ∆PV = 2PVXe2F + 3PVXeF4 – 7PVXeF2, where P and V are the pressure and volume per formula unit, respectively. (b) Calculated Bader charge of Xe in XeFx at 100 GPa. (c) Calculated electronic DOS of XeF2 used by normal Perdew–Burke– Ernzerhof and hybrid functional Heyd–Scuseria–Ernzerhof (HSE06) at 100 GPa, respectively.

We discuss the energetic driven mechanism on the reaction of 7XeF2 → 2Xe2F + 3XeF4 at high pressure through calculations of enthalpies (H = U + PV), PV terms, and static energies (U) of XeF2, Xe2F, and XeF4 as shown in Figure 6a. Competition between ∆U and ∆PV dominates the chemi-

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

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

media by oxidation–reduction methods: spectroscopic properties of the ion. J. Chem. Soc., Chem. Commun. 1978, 502– 504. 8. Seidel, S., Seppelt, K., van Wüllen, C. & Sun, X. Y. The Blue Xe4+ Cation: Experimental Detection and Theoretical Characterization. Angew. Chem. Int. Ed. 2007, 46, 6717– 6720. 9. Fernández, I. & Frenking, G. Neutral noble gas compounds exhibiting a Xe–Xe bond: structure, stability and bonding situation. Phys. Chem. Chem. Phys. 2012, 14, 14869–14877. 10. Jiménez Halla, C. Ó. C., Fernández, I. & Frenking, G. Is it Possible To Synthesize a Neutral Noble Gas Compound Containing a Ng-Ng Bond? A Theoretical Study of H-NgNg-F (Ng= Ar, Kr, Xe). Angew. Chem. Int. Ed. 2009, 48, 366–369. 11. Somayazulu, M. et al. Pressure-induced bonding and compound formation in xenon-hydrogen solids. Nat. Chem. 2009, 2, 50–53. 12. Hanni, M., Lantto, P., Runeberg, N., Jokisaari, J. & Vaara, J. Calculation of binary magnetic properties and potential energy curve in xenon dimer: Second virial coefficient of 129Xe nuclear shielding. J. Chem. Phys. 2004, 121, 5908– 5919. 13. Wang, Y., Lv, J., Zhu, L. & Ma, Y. Crystal structure prediction via particle-swarm optimization. Phys. Rev. B 2010, 82, 094116. 14. Wang, Y., Lv, J., Zhu, L. & Ma, Y. CALYPSO: A method for crystal structure prediction. Comput. Phys. Commun. 2012, 183, 2063–2070. 15. Zhu, L. et al. Substitutional Alloy of Bi and Te at High Pressure. Phys. Rev. Lett. 2011, 106, 145501. 16. Lv, J., Wang, Y., Zhu, L. & Ma, Y. Predicted Novel HighPressure Phases of Lithium. Phys. Rev. Lett. 2011, 106, 015503. 17. Perdew, J. P. & Wang, Y. Accurate and simple analytic representation of the electron-gas correlation energy. Phys. Rev. B 1992, 45, 13244. 18. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. 19. Kresse, G. & Furthmüller, J. Software VASP, vienna (1999). Phys. Rev. B 1996, 54, 169. 20. Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. 21. Togo, A., Oba, F. & Tanaka, I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl_ {2}-type SiO_ {2} at high pressures. Phys. Rev. B 2008, 78, 134106. 22. Blaha, P., Schwarz, K., Sorantin, P. & Trickey, S. B. Fullpotential, linearized augmented plane wave programs for crystalline systems. Comput. Phys. Commun. 1990, 59, 399– 415. 23. Feng, J., Hennig, R. G., Ashcroft, N. W. & Hoffmann, R. Emergent reduction of electronic state dimensionality in dense ordered Li-Be alloys. Nature 2008, 451, 445–448. 24. Peng, F., Miao, M., Wang, H., Li, Q. & Ma, Y. Predicted Lithium–Boron Compounds under High Pressure. J. Am. Chem. Soc. 2012, 134, 18599–18605. 25. Miao, M.-S. Caesium in high oxidation states and as a pblock element. Nat. Chem. 2013, 5, 846–852. 26. Bartlett, N. Synthesis of chemical compound XePtF6. Proc. Chem. Soc. 1962, 112, 218–218. 27. Cordero, B. et al. Covalent radii revisited. Dalton Trans. 2008, 2832-2838. 28. Bader, R. F. W. Atoms in Molecules: A Quantum Theory. (Clarendon: Oxford, UK, 1990).

found in these compounds. Some bonds, for example the XeXe in Xe2F and XeF, show mixed features of covalent and metallic bonding; and Xe-F in Xe2F show mixed features of ionic and covalent bonding. In summary, using an unbiased structure-search method based on particle-swarm optimization algorithms in combination with DFT calculations, we investigated the phase stabilities and structural changes of various Xe–F systems under high pressure. We identified the formation of six stoichiometric Xe fluorides (Xe2F, Xe3F2, XeF, XeF2, XeF4, and XeF6) with unforeseen structural features. These fluorides may be synthesized experimentally over a wide range of pressures. Surprisingly, we found that XeF2 tends to decompose to XeF4 + Xe2F above 81 GPa. XeF4 remains as a molecular crystal up to 200 GPa, whereas XeF6 can polymerize with a higher coordination number under increasing pressure. On the Xe-rich side, Xe fluorides (Xe3F2, Xe2F, and XeF) with an ultralow oxidation state can form at lower pressure, and Xe–Xe dimers and graphite-like layers that contain the novel Xe–Xe covalent bond were discovered. Here, Xe behaves like a p-block element, and it can form compounds with molecular, covalent, ionic, and metallic features under pressure. Acknowledgements We thank the Natural Science Foundation of China under 11304141 and 11534003, China Postdoctoral Science Foundation under 2016M590033. This work is also sponsored by the Program for Science and Technology Innovation Talents in University of Henan Province Grant No. 17HASTIT015 and Open Project of the State Key Laboratory of Superhard Materials (Jilin University No. 201602). We also acknowledge the support of NSF-funded XSEDE resources (TG-DMR130005) including Stampede cluster maintained by Texas Advanced Computing Center (TACC). Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI:xxx. Equation of state, the enthalpies of formation, crystal structure, electronic properties, phonon dispersion curves, and ELF curves for various stoichiometry under different pressures.

References 1. Räsänen, M. Argon out of thin air. Nat. Chem. 2014, 6, 82– 82. 2. Schwarz, U. & Syassen, K. The behaviour of solid XeF2 under pressure. High Pressure Res. 1992, 9, 47–50. 3. Hoppe, R., Dähne, W., Mattauch, H. & Rödder, K. Fluorination of xenon. Angew. Chem. Int. Ed. 1962, 1, 599–599. 4. Kim, M., Debessai, M. & Yoo, C.-S. Two- and threedimensional extended solids and metallization of compressed XeF2. Nat. Chem. 2010, 2, 784–788. 5. Kurzydłowski, D., Zaleski-Ejgierd, P., Grochala, W. & Hoffmann, R. Freezing in Resonance Structures for Better Packing: XeF 2Becomes (XeF +)(F −) at Large Compression. Inorg. Chem. 2011, 50, 3832–3840. 6. Drews, T. & Seppelt, K. The Xe 2+ ion—preparation and structure. Angew. Chem. Int. Ed. 1997, 36, 273–274. 7. Stein, L., Norris, J. R., Downs, A. J. & Minihan, A. R. Formation of the dixenon cation, Xe 2+, in fluoroantimonate (V)

6 ACS Paragon Plus Environment

Page 6 of 7

Page 7 of 7

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 Letters

29. Becke, A. D. & Edgecombe, K. E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 1990, 92, 5397. 30. Savin, A., Nesper, R., Wengert, S. & Fassler, T. E. The Electron Localization Function. Angeii Chem Int Ed Engl. 1997, 36, 1808–1832. 31. Burdett, J. K. & McCormick, T. A. Electron Localization in Molecules and Solids: The Meaning of ELF. J. Phys. Chem. A 1998, 102, 6366–6372. 32. Steinmann, S. N., Mo, Y. & Corminboeuf, C. How do electron localization functions describe π-electron delocalization? Phys. Chem. Chem. Phys. 2011, 13, 20584. 33. Li, X. & Cai, J. Electron density properties and metallophilic interactions of gold halides AuX 2and Au 2X (X = FI): Ab Initiocalculations. Int. J. Quantum Chem. 2016, 116, 1350–1357. 34. Kisowska, K., Berski, S. & Latajka, Z. The structure and chemical bonding in the N2-CuX and N 2···XCu (X = F, Cl, Br) systems studied by means of the molecular orbital and Quantum Chemical Topology methods. J. Comput. Chem. 2008, 29, 2677–2692. 35. Chen, W., Chen, G.-H., Wu, D. & Wang, Q. BNg 3F 3: the first three noble gas atoms inserted into mono-centric neutral compounds – a theoretical study. Phys. Chem. Chem. Phys. 2016, 18, 17534–17545. 36. Makarewicz, E., Gordon, A. J. & Berski, S. Nature of the Bonding in the AuNgX (Ng = Ar, Kr, Xe; X = F, Cl, Br, I) Molecules. Topological Study on Electron Density and the Electron Localization Function (ELF). J. Phys. Chem. A 2015, 119, 2401–2412. 37. Liu, H.-T. et al. Probing the nature of gold–carbon bonding in gold–alkynyl complexes. Nat. Commun. 2013, 4, 1–7. 38. Levy, H. A. & Agron, P. A. The crystal and molecular structure of xenon difluoride by neutron diffraction. J. Am. Chem. Soc. 1963, 85, 241–242. 39. Meyer, L., Barrett, C. S. & Greer, S. C. Crystal Structure of α‐Fluorine. J. Chem. Phys. 1968, 49, 1902. 40. Zhu, L., Liu, H., Pickard, C. J., Zou, G. & Ma, Y. Reactions of xenon with iron and nickel are predicted in the Earth's inner core. Nat. Chem. 2014, 6, 644–648. 41. Miao, M.-S. et al. Anionic Chemistry of Noble Gases: Formation of Mg–NG (NG= Xe, Kr, Ar) Compounds under Pressure. J. Am. Chem. Soc. 2015, 137, 14122–14128. 42. Siegel, S. & Gebert, E. Crystallographic Studies of XeF2 and XeF4. J. Am. Chem. Soc. 1963, 85, 240–240. 43. Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207–8215. 44. Burns, J. H., Ellison, R. D. & Levy, H. A. The crystal structure of the molecular addition compound xenon difluoridexenon tetrafluoride. Acta Cryst. 1965, 18, 11–16. 45. Ibers, J. A. & Hamilton, W. C. Xenon tetrafluoride: crystal structure. Science 1963, 139, 106–107. 46. Templeton, D. H., Zalkin, A., Forrester, J. D. & Williamson, S. M. Crystal and Molecular Structure of Xenon Tetrafluoride. J. Am. Chem. Soc. 1963, 85, 242–242.

7 ACS Paragon Plus Environment