Calcium (Ca) is a very abundant element in the Universe and a crucial component of many important materials, e.g., animal bones, signaling molecules in cells, and fluorite.1–3 Like other Group II alkaline-earth elements, calcium has a silvery metallic appearance in its pure form, but is not found in this condition in nature because of its high reactivity. Again like other alkaline earth metal atoms, calcium atoms have very low electronegativity and easily lose two valence electrons when bonded to other elements. This corresponds to valence II or oxidation state +2, which is perfectly consistent with the octet rule. Can one imagine a situation in which a calcium atom gives away only one of its valence electrons, i.e., is monovalent? Stable compounds with low-valence calcium have long been sought.4 For instance, low-valence calcium has been explored in clusters, such as LCa(N₂)CaL (where L is a bulky β-diketiminate ligand), which can be characterized crystallographically.5 Two-dimensional CaCl crystals with monovalent calcium ions have been synthesized on reduced graphene oxide (rGO) membranes, revealing many interesting properties.6 The atomic valence is the number of electrons that an atom gives or takes in the process of chemical bonding. It substantially affects the physical and chemical properties of compounds—in particular, their color, stability, electrical conductivity, and redox properties. Among the many examples are the oxides of chromium, including the black semiconducting chromium(II) oxide CrO with strongly reducing properties, the green insulating Cr₂O₃ with Cr(III), the half-metallic CrO₂ with Cr(IV), and the red insulating CrO₃ with Cr(VI) and strong oxidizing properties. In contrast to most transition metal atoms, calcium has long been assumed to always be in a divalent state, Ca(II), in its compounds. The discovery of compounds with low-valence calcium will not only improve our knowledge of chemical laws, but should also be valuable for the facilitation of unconventional chemical reactions and the synthesis and exploration of exotic materials and their properties (e.g., for catalysis, ferromagnetism, and hydrogen storage). For example, variable valence of calcium could lead to new catalysts, and if the remaining electron is localized, the monovalent calcium atom Ca(I) will be magnetic.6

Unexpected stoichiometries can be synthesized under extreme conditions (high pressures and temperatures). Recently, pressure–composition phase diagrams have been predicted for many binary systems, including Na–Cl, Ca–C, Na–He, Zn–N, Ti–O, Ta–N, S–O, Rb–B, Fe–He, Ce–N, Fe–Bi, Fe–F, and others,7–18 and these theoretical predictions have provided the basis for the successful experimental synthesis of previously unexpected compounds such as Na₃Cl, Ca₂C₃, and Na₂He under high pressures. Such exotic compounds may exhibit unusual electronic, magnetic, optical, and mechanical properties due to their unusual electronic distributions and peculiar bonding,19 for example, the best known high-temperature superconductors belonging to this class of exotic high-pressure compounds, including LaH₁₀, YH₆, ThH₁₀, H₃S, BaH₁₂, and Sr(Hf, Zr)H₁₈.20–25

In the present work, systematic searches for stable compounds in the calcium–chlorine system are carried out at pressures of up to 100 GPa using the ab initio evolutionary algorithm USPEX.26–28 Besides the anticipated calcium(II) compound CaCl₂, we find three new compounds, CaCl, Ca₅Cl₆, and Ca₃Cl₄, which are predicted to be stable at high pressures according to our calculations. CaCl is stable at pressures above 18 GPa with monovalent Ca cations. It is predicted to undergo a phase transition from Fm-3m-CaCl to Pm-3m-CaCl at 40 GPa. Ca₅Cl₆ is stable above 37 GPa with space group P-1. Ca₃Cl₄ is stable above 73 GPa with space group R-3. We also analyze the band structures, densities of states (DOS), electron localization functions (ELFs),29 and many other properties of these compounds, which exhibit many interesting characteristics. Following our predictions, we successfully synthesize Pm-3m-CaCl in laser-heated diamond anvil cell (DAC) experiments. Our results indicate that the emergence of an unexpected valence state at high pressure has great potential for creating entirely new materials.

Stable stoichiometries and crystal structures were predicted using USPEX26–28 in combination with the Vienna Ab initio Simulation Package (VASP) code30 within the framework of density functional theory (DFT) and the generalized gradient approximation (GGA), using the exchange–correlation functional of Perdew, Burke, and Ernzerhof (PBE).31 We employed the projector-augmented wave method32 with 3s²3p⁶4s² and 3s²3p⁵ electrons treated as the valence electrons for Ca and Cl atoms, respectively.

In these calculations, all stoichiometries were allowed, with the only constraint being that randomly produced structures contained no more than 50 atoms in the primitive cell. The initial population of 120 structures and 20 compositions evolved toward the fittest states (with the fitness being calculated as the vertical distance of each candidate solution from the thermodynamic convex hull), and new structures/compositions were produced from the fittest individuals of the previous generation automatically using heredity, transmutation, and softmutation. Furthermore, randomly produced structures were added to each generation using symmetric28 and topological33 random structure generators.

Our searches proceeded in two steps. First, we performed full searches in the Ca–Cl system and found that CaCl₂ is a stable compound at all pressures studied here. We then performed searches on the calcium-rich side, i.e., for the Ca–CaCl₂ system. The enthalpy of formation per atom of Ca_nCl_m is defined asΔHCanClm=HCanClm−n−m2HCa−m2HCaCl2/(m+n),where all enthalpies H are given at the same pressure and zero temperature. At a given pressure, compounds located on the convex hull are thermodynamically stable rather than decomposing into any isochemical mixture of compounds, while the compounds located above the convex hull are metastable or unstable.

The most interesting structures were further relaxed using the VASP code at high pressures, and their properties were calculated to higher precision with a basis set cutoff of 800 eV and a uniform k-point grid with k-point spacing below 2π × 0.02 Å⁻¹. The energy differences, stress tensors, and structural parameters showed excellent convergence in our calculations. We computed the phonon spectra of all the predicted stable structures using the PHONOPY code to check their dynamical stability.34 We used VASP and the PBE functional to determine the band structures and DOS of these compounds. Bader analysis was performed to explore chemical bonding and analyze the charge transfer.35 The VASPKIT code was used to post-process of the data produced by the VASP code.36 Crystal orbital Hamilton population (COHP) analysis was performed using LOBSTER.37 Crystal structures were visualized using VESTA.38

High-purity CaCl₂ (Alfa Aesar, ultradry, 99.99%) powder and a calcium plate (Alfa Aesar, 99.5%) were loaded into a DAC with culets of 200 µm inside an inert Ar glovebox. The sample chamber was a 100 μm-diameter hole drilled in a pre-indented rhenium gasket using a drilling system at the Center for High Pressure Science and Technology Advanced Research (HPSTAR). The calcium plate in the sample chamber was surrounded by CaCl₂ powder positioned on the diamond anvil. The samples were compressed to 50 GPa and then heated to about 2000 K for 3–5 min over the entire sample by a laser-heating (1064 nm) system at HPSTAR. Double-sided laser heating was used to apply sufficient heating to the samples. During laser heating, the temperatures of the samples were calculated by collecting the emitted thermal radiation, correcting for the optical system response, and fitting the spectral data to the Wien function. The pressure was calibrated from the fluorescence of ruby balls placed inside the sample chamber.39

After laser heating, the X-ray diffraction (XRD) data for samples at different pressures (20–60 GPa) were collected. Synchrotron XRD data were collected on the synchrotron beamlines 4W2 at the Beijing Synchrotron Radiation Facility and BL15U1 at the Shanghai Synchrotron Radiation Facility (λ = 0.6199 Å). The XRD patterns were integrated using Dioptas,40 and full profile refinements were performed using FullProf.41

We explored the pressure–composition phase diagram of the Ca–CaCl₂ system at pressures of 0, 30, 60, and 100 GPa. By evaluating the enthalpy of formation ΔH for each composition at different pressures, the convex hulls for the Ca–CaCl₂ system at different pressures were constructed and are shown in Fig. 1(a). We then carried out more careful calculations of the enthalpy under different pressures for the stable structures on the convex hulls, which enabled us to construct the pressure–composition diagram of the Ca–CaCl₂ system [Fig. 1(b)].

Elemental Ca is predicted to undergo several phase transitions induced by pressure: fcc-bcc-I4₁/amd-P4₃2₁2. In our calculations, under ambient conditions, Ca adopts the fcc structure, while bcc-Ca is stable from 7.5 to 33 GPa, followed by transitions to I4₁/amd and P4₃2₁2 at 33 and 78 GPa, respectively. These results agree with those of previous studies.42–45 Furthermore, the pressure–composition phase diagram predicted by our calculations contains a surprisingly large number of new stable compounds. The phase diagram shows that CaCl₂ remains stable over the whole pressure range investigated here (0–100 GPa), and many new compounds become stable at elevated pressures. We found that CaCl, Ca₅Cl₆, and Ca₃Cl₄ are thermodynamically stable at high pressures [Fig. 1(b)], and are all located on the convex hull from 18, 37, and 73 GPa, respectively, to at least 100 GPa. Calculations reveal that Ca, CaCl₂, and CaCl undergo structural phase transitions at high pressures. The crystallographic parameters of all stable phases are listed in Table S1 (supplementary material), and crystal structures are shown in Fig. 2. Lattice dynamics calculations suggest that all the phases shown in the phase diagram [Fig. 1(b)] have no phonons with imaginary frequencies (i.e., are dynamically stable) in their predicted pressure ranges of thermodynamic stability (Fig. S1, supplementary material). To investigate the electronic properties of these compounds, we calculated their ELFs, electronic band structures, Bader charges, integrated crystal orbital Hamiltonian populations (ICOHPs), and DOS. As shown in Fig. S2 (supplementary material), the electronic DOS shows that all these new compounds (CaCl, Ca₅Cl₆, and Ca₃Cl₄) are metallic.

CaCl₂ is the only known thermodynamically stable compound in the Ca–Cl system. From standard chemical considerations such as the octet rule, CaCl₂ is the only expected compound of calcium and chlorine. CaCl₂ has been studied at high pressures, and experiments have found that within a narrow pressure range of 0–10 GPa, it crystallizes in five different phases with the CaCl₂, α-PbO₂, EuI₂, SrI₂, and PbCl₂ (cotunnite) structure types, with the PbCl₂-type structure being stable up to at least 25 GPa.46 In our calculations, CaCl₂ is a stable stoichiometry in the entire region of explored pressures, and we find a new post-cotunnite P-62m structure [Fig. 2(a)] to be stable in the 18–82 GPa pressure range. Figure S3(a) (supplementary material) shows the band structure and the projected density of states (PDOS) of P-62m-CaCl₂. It has a direct DFT bandgap of 2.71 eV, in contrast to the zero-pressure structure of Pnnm-CaCl₂, which has a direct DFT bandgap of 5.34 eV. The DOS near the Fermi level is dominated by the p electronic states of chlorine. Interestingly, above 82 GPa, P-62m-CaCl₂ transforms back into the cotunnite structure, which remains stable up to at least 100 GPa.

From our calculations, a new weakly metallic compound, CaCl, becomes stable at 18 GPa. It has two stable phases, which have a rocksalt (or B1) structure (space group Fm-3m) and a CsCl-type (or B2) structure (space group Pm-3m), as shown in Figs. 2(b) and 2(c). The phase transition from the rock salt to the CsCl-type structure, accompanied by an increase in the coordination number from 6 to 8, is commonly observed at high pressures.7 For CaCl, B1–CaCl is stable at pressures from 18 to 40 GPa, while B2–CaCl is stable over a wide range of pressures from 40 GPa to at least 100 GPa. Phonon calculations show that the dynamical stability of B1–CaCl is retained as the pressure decreases to 10 GPa and it is not quenchable to ambient pressure, but B2–CaCl retains its dynamical stability when quenched to 0 GPa, which indicates that once formed at high pressure, this phase might be quenchable to ambient pressure.

To determine the valence of Ca in CaCl, we computed the Bader charges of the Ca and Cl atoms (Table S1, supplementary material). To interpret these values, we recall that in all cases known so far, Bader charges have smaller absolute values than the corresponding chemical valences. For example, Bader charges of 1.3–1.4 on Ca atoms in CaCl₂ are compatible with the classical divalent state, whereas the charge of 0.9 in Pm-3m-CaCl can be interpreted as evidence of the Ca(I) state. The intermediate charge of 1.1 in Fm-3m-CaCl is possibly a sign of the mixed-valence state. These unusual valence states are expected to carry many new properties. As previous studies have indicated, the pressure-induced change of orbital energies leads to occupation of orbitals with a higher angular momentum at high pressures.47,48 As the Ca atom is compressed, the occupation of 3d orbitals increases at the expense of 4s orbitals, and this 4s → 3d electron transfer in calcium under pressure makes calcium an “incipient transition metal.” The pressure-driven s–d transfer enables monovalent calcium to exist and stabilizes these new compounds. This is proven by the electronic band structure and PDOS, which indicate that the metallicity of CaCl phases is predominantly due to the 3d electrons of calcium [Figs. 3(a) and 3(b)]. In forming CaCl, a calcium atom gives one electron to chlorine, while the second valence electron of the Ca atom can either remain on the atomic d orbital of calcium (in this case, calcium is monovalent, giving Ca⁺Cl⁻) or delocalize over the whole crystal as an electron gas (in this case, the valence of calcium is two). Below, we argue that in CaCl, we have a mixture of these two situations, but with an essential contribution from the monovalent configuration.

The ELF provides additional information on the electronic properties of CaCl. As shown in Figs. 3(c) and 3(d), the valence ELFs of Fm-3m-CaCl and Pm-3m-CaCl are mainly concentrated on chlorine atoms. Further insight is provided by ICOHP values, which are negative for bonding interactions, and we see significant Ca–Cl and Ca–Ca bonding in both modifications of CaCl (Table S2, supplementary material). Both the metallicity and the presence of Ca–Ca bonding indicate some activity of the second valence electron of calcium (and hence a contribution of the divalent state). There is also a significant degree of localization of that electron, conferring new properties to the compound (i.e., the essential contribution of the monovalent state). As discussed above, CaCl undergoes a phase transition from a rock salt to a CsCl-type structure at 40 GPa. Interestingly, this phase transition is accompanied by a transition from a Pauli paramagnet to a ferromagnet. Usually, increasing pressure leads to a suppression of magnetism, because atoms in the high-spin state have larger radii than in low-spin states. The opposite, however, can also occur—for example, potassium has been predicted to produce s-band ferromagnetism under pressure.49 Here, ferromagnetism comes from the second valence electron of Ca localized on the d orbital. However, the magnetic moment of Pm-3m-CaCl decreases with pressure from 0.25 µB per Ca atom at 40 GPa to 0.04 µB per Ca atom at 100 GPa and eventually collapses at 110 GPa (Fig. 4).

From our calculations, as pressure increases, the atomic volumes of Ca and Cl become similar, making a Ca → Cl substitution possible. This indeed happens and makes Ca₅Cl₆ and Ca₃Cl₄ thermodynamically stable. Ca₅Cl₆ is a low-symmetry (space group P-1) phase, thermodynamically stable from 37 GPa to at least 100 GPa. It has only one ground-state structure under the whole pressure range explored here. Its crystal structure can be obtained from that of a Pm-3m-CaCl crystal by replacing Ca atoms with Cl at some sites [Fig. 2(d)]. The Bader charges of Cl atoms occupying the Ca position are smaller (−0.7|e|) than those of other Cl atoms (−1|e|). The other stable compound, Ca₃Cl₄, is predicted to have a space group R-3 and is thermodynamically stable at pressures ranging from 73 GPa to at least 100 GPa. Its crystal structure can also be obtained from that of a Pm-3m-CaCl crystal just like Ca₅Cl₆ [Fig. 2(e)]. To further investigate the properties of Ca₅Cl₆ and Ca₃Cl₄, we analyzed their electronic structures. As expected, both Ca₃Cl₄ and Ca₅Cl₆ are metallic, and their DOS at the Fermi level come mainly from Ca 3d electrons, similar to CaCl [Figs. S3(c) and S3(e), supplementary material]. Figures S3(d) and S3(f) show the similarity of the ELF distributions of Ca₃Cl₄ and Ca₅Cl₆ to that of CaCl, whose valence electrons are mainly concentrated around the Cl atoms.

Having understood these compounds theoretically, we decided to experimentally verify their stability. We performed high-pressure experiments on the Ca–Cl system in a laser-heated DAC to validate our predictions. A Ca plate and CaCl₂ powder were loaded into a Re gasket hole of the DAC, then compressed gradually to 50 GPa, and then laser-heated to a temperature of about 2000 K. A minute change in pressure in the sample chamber before and after the laser-heating process was observed and was probably caused by changes in sample volume due to chemical reaction. At this pressure, the formation of Pm-3m-CaCl is expected on the basis of its thermodynamic stability [Fig. 1(a)].

Powder XRD patterns changed greatly after laser heating. Figure 5(a) shows powder XRD data obtained at 50 GPa after laser heating. There is a good match between the experimental data and simulated diffraction patterns. Most of the reflections do indeed correspond to the predicted P-62m-CaCl₂ and Pm-3m-CaCl, and the most obvious peak of Pm-3m-CaCl belongs to (110) [Fig. 5(a)]. Minor differences between calculated and experimental XRD patterns are likely due mostly to impurities and anisotropic broadening of diffraction peaks, probably caused by inhomogeneous elastic stresses, which always arise under high pressures during laser heating. It is worth noting that some of the new diffraction rings are spotty (Fig. S4, supplementary material), indicating that we do not have a perfect powder with a very large number of randomly oriented crystallites. This texture is due to the heating and cooling process of the samples. Importantly, we did not observe elemental calcium in the products, which means that the elemental calcium had completely reacted.

On changing the pressure on the sample, we found that Pm-3m-CaCl was stable upon quenching down to at least 20 GPa. Additionally, some small peaks, marked with asterisks in Fig. 5(b), were observed during decompression, and these might be contributed by some minority impurities at different sampling positions. The peak position (2θ) of (110) of Pm-3m-CaCl moved from 16.59° at 50 GPa to 15.57° at 20 GPa, which means that the lattice parameter changed from 3.04 to 3.24 Å. After quenching of the sample to ambient pressure and exposing it to air, we could not observe any peaks related to the products, which means that these were unstable at ambient pressure and temperature. From the XRD data, we obtained the pressure–volume equation of state of Pm-3m-CaCl in the pressure range 20–60 GPa [Fig. 5(c)]. There was a good agreement between the experimental and theoretical equations of state of Pm-3m-CaCl. The experimental P–V data were fitted to a second-order Birch–Murnaghan equation of state with B₀ = 42.68 GPa, while B₀ was 44.29 GPa from theoretical predictions. To further exclude the possible presence of unreacted elemental calcium influencing the observation of Pm-3m-CaCl diffraction peaks, we compare the cell volumes of Pm-3m-CaCl and simple cubic Ca in Fig. 5(c). We can see that the volume differences between the two at the same pressures are sufficiently large to completely rule out the presence of elemental calcium in the sample. The fact that the reaction Ca + CaCl₂ = 2CaCl producing CaCl took place only after laser heating implies a considerable activation barrier for this reaction.

In summary, we have studied the Ca–Cl system using state-of-the-art theoretical predictions and experimental high-pressure synthesis. We have produced a pressure–composition phase diagram for this system at pressures of up to 100 GPa and synthesized a new anomalous compound, calcium monochloride (CaCl), containing calcium in a unique monovalent state, Ca(I), as predicted by theory. CaCl is predicted to undergo a pressure-induced transition from a nonmagnetic NaCl-type phase to a ferromagnetic CsCl-type phase at 40 GPa. Ca₅Cl₆ and Ca₃Cl₄ are stable at pressures above 37 and 73 GPa with space groups P-1 and R-3, respectively. The predicted Pm-3m-CaCl was synthesized in a laser-heated DAC at 50 GPa and shown to be quenchable down to at least 20 GPa. The application of pressure enables the creation of a wider palette of valence states of atoms, greatly increasing the range of accessible properties. For example, making calcium a variable-valence element could lead to completely unexpected calcium-based catalysts. This new side of the seemingly simple chemistry of calcium may provide important surprises in geochemistry and planetary sciences, since calcium is one of the most abundant planet-forming elements, and high pressures are typical of planetary interiors. Even under ambient conditions, the Ca(I) state may be more accessible than commonly thought (in intermediate states of chemical reactions, for example). Upon rigorous theoretical and experimental analysis, what was initially thought to be a fundamentally simple chemical system has revealed itself to be a tremendously rich system with novel physics and chemistry. This work demonstrates that by revisiting simple systems, it is possible to develop new chemical principles and identify remarkable new materials and phenomena.

See the supplementary material for the crystal structure and Bader charges of stable Ca–Cl compounds (Table S1), calculated ICOHPs for Ca–Ca, Ca–Cl, and Cl–Cl pairs in CaCl (Table S2), phonon dispersion curves of Fm-3m-CaCl, Pm-3m-CaCl, Ca₅Cl₆, and Ca₃Cl₄ (Fig. S1), electronic DOS of calcium chlorides (Fig. S2), electronic structures of Ca₅Cl₆ at 40 GPa, Ca₃Cl₄ at 80 GPa, and CaCl₂ at 60 GPa (Fig. S3), and a two-dimensional image of the XRD pattern after laser heating (Fig. S4).

Acknowledgment. This work was supported by the National Science Foundation of China (Grant Nos. 92263101, 12174200, 21627802, 51722209, and 21273206), the Science Challenge Project (Grant No. TZ2016001), the Key Research Project of Higher Education (Grant Nos. 15A140016 and 2010GGJS-110), and the National Key R&D Program of China (Grant No. YS2018YFA070119). The work of A.R.O. was supported by the Russian Science Foundation (Grant No. 24-43-00162). The calculations were performed and supported by Tianhe II in Guangzhou and the Supercomputing Center of Nankai University (NKSC).