Modulated structure and molecular dissociation of solid chlorine at high pressures

Data-driven prediction of complex crystal structures of dense lithium

Lithium (Li) is a prototypical simple metal at ambient conditions, but exhibits remarkable changes in structural and electronic properties under compression. There has been intense debate about the structure of dense Li, and recent experiments offered fresh evidence for yet undetermined crystalline phases near the enigmatic melting minimum region in the pressure-temperature phase diagram of Li. Here, we report on an extensive exploration of the energy landscape of Li using an advanced crystal structure search method combined with a machine-learning approach, which greatly expands the scale of structure search, leading to the prediction of four complex Li crystal structures containing up to 192 atoms in the unit cell that are energetically competitive with known Li structures. These findings provide a viable solution to the observed yet unidentified crystalline phases of Li, and showcase the predictive power of the global structure search method for discovering complex crystal structures in conjunction with accurate machine learning potentials.

Multistep Dissociation of Fluorine Molecules under Extreme Compression

All elements that form diatomic molecules, such as H_{2}, N_{2}, O_{2}, Cl_{2}, Br_{2}, and I_{2}, are destined to become atomic solids under sufficiently high pressure. However, as revealed by many experimental and theoretical studies, these elements show very different propensity and transition paths due to the balance of reduced volume, lone pair electrons, and interatomic bonds. The study of F under pressure can illuminate this intricate behavior since F, owing to its unique position on the periodic table, can be compared with H, with N and O, and also with other halogens. Nevertheless, F remains the only element whose solid structure evolution under pressure has not been thoroughly studied. Using a large-scale crystal structure search method based on first principles calculations, we find that, before reaching an atomic phase, F solid transforms first into a structure consisting of F_{2} molecules and F polymer chains and then into a structure consisting of F polymer chains and F atoms, a distinctive evolution with pressure that has not been seen in any other elements. Both intermediate structures are found to be metallic and become superconducting, a result that adds F to the elemental superconductors.

Halogen molecular modifications at high pressure: the case of iodine

Metallization and dissociation are key transformations in diatomic molecules at high densities particularly significant for modeling giant planets. Using X-ray absorption spectroscopy and atomistic modeling, we demonstrate that in halogens, the formation of a connected molecular structure takes place at pressures well below metallization. Here we show that the iodine diatomic molecule first elongates by ∼0.007 Å up to a critical pressure of Pc ∼ 7 GPa, developing bonds between molecules. Then its length continuously decreases with pressure up to 15-20 GPa. Universal trends in halogens are shown and allow us to predict for chlorine a pressure of 42 ± 8 GPa for molecular bond-length reversal. Our findings contribute to tackling the molecule invariability paradigm in diatomic molecular phases at high pressures and may be generalized to other abundant diatomic molecules in the universe, including hydrogen.

Prediction of chlorine and fluorine crystal structures at high pressure using symmetry driven structure search with geometric constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Our algorithm extends the symmetry driven structure search method [R. Domingos et al., Phys. Rev. B 98, 174107 (2018)] by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa-4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P4₂/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.

Band gap closure, incommensurability and molecular dissociation of dense chlorine

Diatomic elemental solids are highly compressible due to the weak interactions between molecules. However, as the density increases the intra- and intermolecular distances become comparable, leading to a range of phenomena, such as structural transformation, molecular dissociation, amorphization, and metallisation. Here we report, following the crystallization of chlorine at 1.15(30) GPa into an ordered orthorhombic structure (oC8), the existence of a mixed-molecular structure (mC8, 130(10)-241(10) GPa) and the concomitant observation of a continuous band gap closure, indicative of a transformation into a metallic molecular form around 200(10) GPa. The onset of dissociation of chlorine is identified by the observation of the incommensurate structure (i-oF4) above 200(10) GPa, before finally adopting a monatomic form (oI2) above 256(10) GPa.

Perspective: crystal structure prediction at high pressures

Crystal structure prediction at high pressures unbiased by any prior known structure information has recently become a topic of considerable interest. We here present a short overview of recently developed structure prediction methods and propose current challenges for crystal structure prediction. We focus on first-principles crystal structure prediction at high pressures, paying particular attention to novel high pressure structures uncovered by efficient structure prediction methods. Finally, a brief perspective on the outstanding issues that remain to be solved and some directions for future structure prediction researches at high pressure are presented and discussed.

High-pressure crystallography of periodic and aperiodic crystals

More than five decades have passed since the first single-crystal X-ray diffraction experiments at high pressure were performed. These studies were applied historically to geochemical processes occurring in the Earth and other planets, but high-pressure crystallography has spread across different fields of science including chemistry, physics, biology, materials science and pharmacy. With each passing year, high-pressure studies have become more precise and comprehensive because of the development of instrumentation and software, and the systems investigated have also become more complicated. Starting with crystals of simple minerals and inorganic compounds, the interests of researchers have shifted to complicated metal-organic frameworks, aperiodic crystals and quasicrystals, molecular crystals, and even proteins and viruses. Inspired by contributions to the microsymposium 'High-Pressure Crystallography of Periodic and Aperiodic Crystals' presented at the 23rd IUCr Congress and General Assembly, the authors have tried to summarize certain recent results of single-crystal studies of molecular and aperiodic structures under high pressure. While the selected contributions do not cover the whole spectrum of high-pressure research, they demonstrate the broad diversity of novel and fascinating results and may awaken the reader's interest in this topic.

Method of increments for the halogen molecular crystals: Cl, Br, and I

Method of increments (MI) calculations reveal the n-body correlation contributions to binding in solid chlorine, bromine, and iodine. Secondary binding contributions as well as d-correlation energies are estimated and compared between each solid halogen. We illustrate that binding is entirely determined by two-body correlation effects, which account for >80% of the total correlation energy. One-body, three-body, and exchange contributions are repulsive. Using density-fitting (DF) local coupled-cluster singles, doubles, and perturbative triples for incremental calculations, we obtain excellent agreement with the experimental cohesive energies. MI results from DF local second-order Møller-Plesset perturbation (LMP2) yield considerably over-bound cohesive energies. Comparative calculations with density functional theory and periodic LMP2 method are also shown to be less accurate for the solid halogens.

Condensed astatine: monatomic and metallic

The condensed matter properties of the nominal terminating element of the halogen group with atomic number 85, astatine, are as yet unknown. In the intervening more than 70 years since its discovery significant advances have been made in substrate cooling and the other techniques necessary for the production of the element to the point where we might now enquire about the key properties astatine might have if it attained a condensed phase. This subject is addressed here using density functional theory and structural selection methods, with an accounting for relativistic physics that is essential. Condensed astatine is predicted to be quite different in fascinating ways from iodine, being already at 1 atm a metal, and monatomic at that, and possibly a superconductor (as is dense iodine).