Stability of numerous novel potassium chlorides at high pressure.

K-Cl is a simple system displaying all four main types of bonding, as it contains (i) metallic potassium, (ii) elemental chlorine made of covalently bonded Cl2 molecules held together by van der Waals forces, and (iii) an archetypal ionic compound KCl. The charge balance rule, assigning classical charges of "+1" to K and "-1" to Cl, predicts that no compounds other than KCl are possible. However, our quantum-mechanical variable-composition evolutionary simulations predict an extremely complex phase diagram, with new thermodynamically stable compounds K3Cl, K2Cl, K3Cl2, K4Cl3, K5Cl4, K3Cl5, KCl3 and KCl7. Of particular interest are 2D-metallic homologs Kn+1Cln, the presence of positively charged Cl atoms in KCl7, and the predicted stability of KCl3 already at nearly ambient pressures at zero Kelvin. We have synthesized cubic -KCl3 at 40-70 GPa and trigonal -KCl3 at 20-40 GPa in a laser-heated diamond anvil cell (DAC) at temperature exceeding 2000 K from KCl and Cl2. These phases were identified using in situ synchrotron X-ray diffraction and Raman spectroscopy. Upon unloading to 10 GPa, -KCl3 transforms to a yet unknown structure before final decomposition to KCl and Cl2 at near-ambient conditions.

Recent ab initio calculations predicted the formation of unexpected novel high-pressure compounds in several simple systems, such as Li (Na)-H12 Mg-O3, and Na-Cl4. These systems were subsequently explored experimentally: while so far the predictions have not been verified for Li-H5, for Na-Cl and Mg-O, the predicted compounds (NaCl3, Na3Cl, and MgO2) have been confirmed experimentally46, while more stable NaHx compounds than originally predicted by Zurek et al.2 have been synthesized and theoretically verified7, revealing dramatic changes of chemistry under pressure. Here we study K-Cl, a system closely related to Na-Cl and find even richer chemistry and new phenomena.

The only known potassium chloride, KCl, has been extensively studied under pressure, both experimentally8910 and using ab initio simulations111213. Two crystal structures are known for KCl: the rocksalt-type (B1) structure and cesium chloride-type (B2) structure, the latter becoming stable at ~2 GPa. The same transition occurs in NaCl, but at a much higher pressure of 30 GPa1415, reflecting the general tendency for phase transitions to occur at lower pressures for compounds of heavier elements. Yet, as we find, the K-Cl system has a much richer chemistry than Na-Cl. Here we study the K-Cl system using the quantum-mechanical variable-composition evolutionary structure prediction methodology USPEX16171819, searching for stable compounds and their corresponding crystal structures (see Methods). In each of these calculations, all possible chemical compositions were allowed with up to 16 atoms in the unit cell, and calculations were performed at pressures of 1 atm, 10 GPa, 35 GPa, 50 GPa, 100 GPa, 150 GPa, 200 GPa, 250 GPa and 300 GPa. Theoretical predictions were successfully verified by experimental synthesis of two KCl3 polymorphs in a laser-heated diamond anvil cell (DAC).

Results and Discussion

The pressure-composition phase diagram (see Fig. 1, Supplementary Fig. S1) predicted in our calculations contains a surprisingly large number of new stable compounds (see Supplementary Table S1). By thermodynamically stable we mean a compound which is more stable than any isochemical mixture of the elements or other compounds – this definition leads to the convex hull construction shown in Fig. 1b. The dynamical stability of the newly predicted phases was confirmed by phonon calculations (see Supplementary Fig. S2).

Figure 1

Stability of new potassium chlorides: (a) Pressure-composition phase diagram of the K-Cl system. (b) Convex hull diagrams for the K-Cl system at selected pressures. Solid circles represent stable compounds; open circles - metastable ones.

To verify these predictions, we performed high-pressure experiments on the K-Cl system in a laser-heated DAC up to 70 GPa in the presence of excess chlorine. We specifically targeted synthesis of KCl3, which was predicted to become stable at the lowest pressures. The reaction products were examined by visual observations (see Supplementary Fig. S3), Raman confocal spectroscopy, and synchrotron x-ray diffraction (XRD) probes at room temperature. Combining experimental and theoretical approach was critical to refine the K-Cl phase diagram as several KCl3 phases showed competing enthalpies in the 0–30 GPa pressure range.

The phase diagram shows that KCl remains stable in the whole pressure range investigated here, but many new compounds become stable at elevated pressures. Perhaps most unexpected is the prediction that KCl3 is stable already at 1 atm and 0 K. The structure belongs to P31m space group with 3 formula units (f.u.) in the unit cell (Fig. 2a), and contains exotic trichloride-ions Cl3−. The P31m phase is a semiconductor, with a DFT band gap of 2.60 eV. The first phase transition is to the Pnma structure at 1.3 GPa. Pnma structure has 4 f.u. in the unit cell and also contains trichloride-ions (Fig. 2b). Bader analysis gives the charge configuration K+0.83[Cl−0.28Cl−0.04Cl−0.51]−0.83, nearly the same as for Pnma-NaCl3 [ref. 4]. [Cl3]− ion is an isoelectronic analogue of the well-known triiodide-ion [I3]− (for example, compound KI3 is well known), Br3− and ClICl− ions, and can be also related to the known [Li3]− [ref. 20] and hypothetical [H3]− [ref. 5] ions. At 9.3 GPa, P  structure of KCl3 with 6 f.u. (Fig. 2c) in the unit cell becomes stable. The DFT band gap of P -KCl3 is 1.78 eV at 20 GPa. Interestingly, at P > 160 GPa, P -KCl3 turns metallic due to the band gap closure. Metallic -KCl3, isostructural with stable P -NaCl3, is also energetically competitive under pressure. K3Cl5 and KCl7 become stable at pressure above 140 GPa and 225 GPa, respectively.

Figure 2

Crystal structures of (a) P31m-KCl3 at 1atm. (b) Pnma-KCl3 at 5 GPa. (c) P -KCl3 at 20 GPa. (d) -KCl3 at 240 GPa. (e) Pm3-KCl7 at 240 GPa. (f) P-K3Cl5 at 240 GPa.

We synthesized KCl3 at elevated pressures and temperatures at conditions of excess of Cl2. To overcome the kinetic barriers, the reagents were laser-heated above 2000 K. The temperature was determined radiometrically. This heating procedure also promotes better mixing of reagents as chlorine melts and becomes highly diffusive. Pressures in excess of 20 GPa were necessary to initiate a chemical reaction between KCl and Cl2. A set of new Bragg peaks was observed after laser heating at 20–40 GPa with intensities and angular positions in agreement with P -KCl3 (Fig. 3a). Rietveld refinement, however, was not possible because of the apparent texturing of new reflections (see Supplementary Fig. S4). A rich Raman spectrum, with at least 15 peaks (Fig. 3b), was observed for the synthesized compound, which is consistent with group theory allowing 16 Raman active modes (Г = 5A1g + 11 Eg) for P -KCl3. Likewise, Raman shift of the experimentally observed bands agrees well with that computed for P -KCl3 at corresponding pressures (see Supplementary Fig. S5). Typical agreement of vibrational frequencies from DFT and experiment for materials with well-established structure is up to 10%, e.g. ref. 21. Therefore, experiments confirm that P -KCl3 is the most stable phase in the 20–40 GPa pressure range.

Figure 3

Experimental evidence for KCl3: (a) XRD pattern of P -KCl3 and Cl2 at 21 GPa. (b) Raman spectrum of P -KCl3 and Cl2 at 22 GPa. Blue bars show computed spectral positions and intensities; Red ticks represent tentative assignment of the theoretically predicted Raman modes with the experimental data. (c) XRD pattern of -KCl3, KCl, and Cl2 at 57 GPa. (d) Experimental and theoretical pressure-volume equations of state of P  and -KCl3. In (a,c) black lines show the intensity difference (Iobs − Icalc), Le Bail refinement residuals are Rwp = 0.139 and Rexp = 0.096 in (a) and Rwp = 0.233 and Rexp = 0.151 in (c). X-ray wavelengths are 0.3100 Ǻ in (a) and 0.3344 Ǻ in (c).

XRD of quenched samples prepared at P > 35–40 GPa shows a mixture of two space-separated phases of KCl3: P  and . Larger yields of the phase were achieved at P > 50–60 GPa (Fig. 3c) in qualitative agreement with theoretical predictions showing that the energy difference between P and -KCl3 decreases with increasing pressure. We could only use the Le Bail refinement of -KCl3 because of the spotty character of XRD images (see Supplementary Fig. S6). The agreement between the experimentally measured and computed equations of state of KCl3 (Fig. 3d) further validates theoretical predictions. It is remarkable that in a number of experiments KCl reacted completely, forming KCl3, with the only remaining material in the probed area being Cmca-chlorine, which was easily identified based on experimental22 and our theoretically calculated lattice parameters.

On decompression to below 10 GPa, Raman bands of P -KCl3 disappeared completely, while new strong pressure-dependent bands appeared near 450 cm−1 (see Supplementary Fig. S7). We tentatively assigned these bands to stretching vibrations of the linear Cl3− ions23. Changes in XRD also suggest a phase transition, although the quality of the diffraction pattern was not sufficient to index new peaks and pinpoint the structure. At room temperature, this new phase becomes unstable below 2 GPa: Raman spectroscopy, X-ray diffraction and visual observations showed only the presence of KCl and Cl2 in the decompressed sample cavity.

In the K-Cl system, in contrast with Na-Cl, there is yet another chlorine-rich phase, P m2-K3Cl5, which has a pseudocubic cell with 1 formula unit. The K atom in the center of the unit cell is surrounded by 4 K atoms and 10 Cl atoms, together forming a bicapped hexagonal antiprism (Fig. 2f). The electronic structure (Fig. 4a,c) shows that it is a poor metal with a deep pseudogap of width ~4.6 eV at 240 GPa. In Fig. 4, we compare the total and atom-projected electronic densities of states of P m2-K3Cl5, -KCl3 and Pm3-KCl7. All these phases are poor metals with pronounced pseudogaps at the Fermi level, implying electronic stabilization. The main contribution at the Fermi level comes from chlorine atoms, and one can observe that different chlorine sites play very different roles – for example, in Pm2-K3Cl5 only p-orbitals of Cl (4j) contribute at the Fermi level, and are thus responsible for its metallicity. Due to excess of chlorine atoms, which act as electron acceptors, -KCl3 has DOS similar to p-type semiconductors. The central, positively charged Cl (1b) donating electrons to the system in Pm3-KCl7, makes the DOS at the Fermi level in KCl7 much higher than that in Pm2-K3Cl5 and -KCl3. Distributions of valence electron localization function (ELF, e.g., Fig. 4b) also show that crystallographically inequivalent Cl atoms have very different ELF distributions – from spherical (around atoms with the most negative Bader charge, indicating a closed-shell configuration, and also around the positively charged Cl atom in KCl7) to toroidal (around atoms with small negative charges).

Figure 4

Electronic structure: (a) band structure and electronic density of states of P-K3Cl5 at 240 GPa. (b) electron localization function of P-K3Cl5 at 240 GPa with isosurface ELF = 0.77. (c) total and atom-projected densities of states of P-K3Cl5. (d) total and atom-projected densities of states of -KCl3. (e) total and atom-projected densities of states of Pm3-KCl7 at 240 GPa.

Comparing Bader charges of K-Cl phases (see Supplementary Table S2) with those of Na-Cl phases, we see higher charges on Na atoms in Pnma-NaCl3, -NaCl3 and Pm3-NaCl7 (about +0.8) than in their K-counterparts (about +0.65). This is counterintuitive, but consistent with our finding24 that under pressure K has higher electronegativity and lower reactivity than Na, due to the well-known s → d electronic transition in K atoms under pressure. Related to this is the observation that the depth of the convex hull (i.e. the enthalpy of formation of KCl or NaCl) in the K-Cl system (Fig. 1) changes from −2.9 eV/atom at 40 GPa to −1.5 eV/atom at 300 GPa, whereas for the Na-Cl system3, it changes from −2.5 eV/atom at 40 GPa to −2.9 eV/atom at 300 GPa.

In the studied pressure range, besides the known B1 and B2 phases, we find a new phase of KCl: I41/amd–KCl, stable above 201 GPa, shown in Fig. 5a. This structure is a derivative of the fcc structure. Figure 5b shows K3Cl, the other fcc-derived superstructure compound stable in the K-Cl system (above 149 GPa) – square planar layers with stoichiometry KCl alternate with similar layers of stoichiometry K2 along the c-axis, leading to the total stoichiometry K3Cl. These two compounds can be described as fcc-based homologs.

Figure 5

Crystal structures of (a) high-pressure I41/amd-KCl, (b) fcc-derived I4/mmm-K3Cl, and bcc-derived Kn+1Cln homologs: (c) I4/mmm-K2Cl, (d) I4/mmm-K3Cl2, (e) I4/mmm-K4Cl3, (f) I4/mmm-K5Cl4. (g) Spatial distribution of electrons (shown by isosurfaces and density contours) at the Fermi level in I4/mmm-K5Cl4, showing clear 2D-metallic character.

There is another interesting and surprisingly rich class of phases, Kn+1Cln homologs (n = 2, 3, 4 were found in our calculations, but we cannot exclude the possibility of even higher homologs) based on the B2 structure and shown in Fig. 5c–g. These have (2n + 1) layers along the c-axis, with extra K-layer serving as an antiphase boundary between B2-structured domains. All these phases have the same space group I4/mmm, and similar interatomic distances, all of them are poor metals, due to the excess of electron-donating K atoms, analogous to n-type semiconductors (see Supplementary Figs S8 and S9), and display a two-dimensional metallic character. It is surprising that phases with different n, while being structurally so similar, have quite different stability fields: e.g., K2Cl is stable at pressures above 56 GPa, whereas K5Cl4 is stable above 100 GPa. Interestingly, mobile electrons are observed only at the antiphase boundaries, whereas regions between them are insulating (Fig. 5g). These antiphase boundaries may be created as metastable growth defects also at lower pressures, with the promise of new electronic materials.

In summary, for a seemingly simple K-Cl system our calculations predict an extremely unusual behavior. Already at near-ambient conditions we have predicted and experimentally confirmed the existence of a new insulating compound KCl3, which has not been observed before. As pressure increases, a surprisingly large number of thermodynamically stable phases become stable: (1) Cl-rich metallic phases (KCl7, K3Cl5, and a metallic form of KCl3) with high coordination numbers (12–14), (2) fcc-superstructures (insulating I41/amd-KCl and metallic K3Cl), (3) layered B2-superstructures with compositions KCl (n = 3, 4, 5) and two-dimensional electronic conductivity. KCl3 decomposes into KCl and Cl2 at room temperature at pressures below 2 GPa and probably is stable at zero pressure and low temperatures, as suggested by theory. What was considered as an ultimately simple chemical system, upon careful theoretical and experimental study turned out to be a very rich system with novel physics and chemistry. Revisiting other simple systems may result in the formulation of new chemical principles that could be used for the discovery of novel materials and phenomena.

Methods

Theory

Structure/composition predictions were done using the USPEX code161718 in the variable-composition mode19. The first generation of structures was produced randomly and the subsequent generations were obtained by applying heredity, transmutation, softmutation, and lattice mutation operations, with probabilities of 60%, 10%, 20% and 10%, respectively. 60% fittest non-identical structures of each generation were used to produce the next generation. 20% new random symmetric structures were also added in each generation. All structures were relaxed using density functional theory (DFT) calculations at the generalized gradient approximation level of theory, with the Perdew-Burke-Ernzerhof (PBE)25 exchange-correlation functional, as implemented in the VASP code26. We used the all-electron projector augmented wave (PAW)27 with K [3s23p64s1], Cl [2s22p4] cores (core radii 2.20 a.u. and 1.50 a.u., respectively) and plane-wave basis sets with the 500 eV cutoff, and dense Monkhorst-Pack meshes with resolution better than 2π × 0.05 Ǻ−1. We used the normalized enthalpy of formation as fitness and visualized crystal structures and electron density distributions using the STM4 package28. Having identified the most stable compositions and structures, we relaxed them at pressures between 1 atm and 300 GPa using very accurate Brillouin zone sampling (Monkhorst-Pack meshes with resolution of better than 2π × 0.03 Ǻ−1).

Experiments

We used symmetrical diamond anvil cells (DAC) to generate static pressures up to 70 GPa. The diamond culet size was 300 μm. Sample chambers were created by laser drilling of 70–80 μm holes in Re gaskets (40 μm thick). KCl platelets (8–15 μm thick), were made by squeezing KCl single crystals (which were dried out at 130 °C) and cleaved to match the dimensions of the gasket hole. Immediately after the cleaving, KCl platelets were stacked at an angle preserving empty space (5–15 μm) between the platelets for chlorine condensation in the gap. Subsequently, DACs were cooled to 77 K with liquid nitrogen in a nitrogen-purged glove box. Chlorine gas (Linde Group, >99.8%) was injected into the gasket cavity from a capillary. Finally, the DAC was closed and brought to high pressure while still at cryogenic temperatures.

Upon compression, chlorine becomes optically opaque and can be heated directly by a 1064 nm fiber laser. Double-sided laser heating experiments were performed at GeoSoilEnviroCARS (APS, Chicago) and Extreme Conditions Beamline P02.2 at DESY (Germany). Temperature was measured spectroradiometrically. XRD data were collected both at high temperature and from quenched samples. We used a 5 μm x-ray beam to detect phase transformations (chlorine melting, KCl3 synthesis) at high temperature, while 2–3 μm beam was used to map quenched samples.

Raman radiation was excited using either 488 or 532 nm lines of a solid state laser. The laser spot on the sample was focused to 4 μm. Raman spectra were analyzed with an imaging spectrograph equipped with a charge-coupled device (CCD). The spectral resolution was 4 cm−1.

Additional Information

How to cite this article: Zhang, W. et al. Stability of numerous novel potassium chlorides at high pressure. Sci. Rep. 6, 26265; doi: 10.1038/srep26265 (2016).