Computational prediction of new magnetic materials.

The discovery of new magnetic materials is a big challenge in the field of modern materials science. We report the development of a new extension of the evolutionary algorithm USPEX, enabling the search for half-metals (materials that are metallic only in one spin channel) and hard magnetic materials. First, we enabled the simultaneous optimization of stoichiometries, crystal structures, and magnetic structures of stable phases. Second, we developed a new fitness function for half-metallic materials that can be used for predicting half-metals through an evolutionary algorithm. We used this extended technique to predict new, potentially hard magnets and rediscover known half-metals. In total, we report five promising hard magnets with high energy product (|BH|MAX), anisotropy field (Ha), and magnetic hardness (κ) and a few half-metal phases in the Cr-O system. A comparison of our predictions with experimental results, including the synthesis of a newly predicted antiferromagnetic material (WMnB2), shows the robustness of our technique.

A stable compound of helium and sodium at high pressure.

Helium is generally understood to be chemically inert and this is due to its extremely stable closed-shell electronic configuration, zero electron affinity and an unsurpassed ionization potential. It is not known to form thermodynamically stable compounds, except a few inclusion compounds. Here, using the ab initio evolutionary algorithm USPEX and subsequent high-pressure synthesis in a diamond anvil cell, we report the discovery of a thermodynamically stable compound of helium and sodium, Na2He, which has a fluorite-type structure and is stable at pressures >113 GPa. We show that the presence of He atoms causes strong electron localization and makes this material insulating. This phase is an electride, with electron pairs localized in interstices, forming eight-centre two-electron bonds within empty Na8 cubes. We also predict the existence of Na2HeO with a similar structure at pressures above 15 GPa.

Diverse Chemistry of Stable Hydronitrogens, and Implications for Planetary and Materials Sciences.

Nitrogen hydrides, e.g., ammonia (NH3), hydrazine (N2H4) and hydrazoic acid (HN3), are compounds of great fundamental and applied importance. Their high-pressure behavior is important because of their abundance in giant planets and because of the hopes of discovering high-energy-density materials. Here, we have performed a systematic investigation on the structural stability of N-H system in a pressure range up to 800 GPa through evolutionary structure prediction. Surprisingly, we found that high pressure stabilizes a series of previously unreported compounds with peculiar structural and electronic properties, such as the N4H, N3H, N2H and NH phases composed of nitrogen backbones, the N9H4 phase containing two-dimensional metallic nitrogen planes and novel N8H, NH2, N3H7, NH4 and NH5 molecular phases. Another surprise is that NH3 becomes thermodynamically unstable above ~460 GPa. We found that high-pressure chemistry of hydronitrogens is much more diverse than hydrocarbon chemistry at normal conditions, leading to expectations that N-H-O and N-H-O-S systems under pressure are likely to possess richer chemistry than the known organic chemistry. This, in turn, opens a possibility of nitrogen-based life at high pressure. The predicted phase diagram of the N-H system also provides a reference for synthesis of high-energy-density materials.

Combined Theoretical and in Situ Scattering Strategies for Optimized Discovery and Recovery of High-Pressure Phases: A Case Study of the GaN-Nb2O5 System.

The application of pressure in solid-state synthesis provides a route for the creation of new and exciting materials. However, the onerous nature of high-pressure techniques limits their utility in materials discovery. The systematic search for novel oxynitrides-semiconductors for photocatalytic overall water splitting-is a representative case where quench high-pressure synthesis is useful and necessary in order to obtain target compounds. We utilize state of the art crystal structure prediction theory (USPEX) and in situ synchrotron-based X-ray scattering to speed up the discovery and optimization of novel compounds using high-pressure synthesis. Using this approach, two novel oxynitride phases were discovered in the GaN-Nb2O5 system. The (Nb2O5)0.84:(NbO2)0.32:(GaN)0.82 rutile structured phase was formed at 1 GPa and 900 °C and gradually transformed to a α-PbO2-related structure above 2.8 GPa and 1000 °C. The low-pressure rutile type phase was found to have a direct optical band gap of 0.84 eV and an indirect gap of 0.51 eV.

A novel phase of beryllium fluoride at high pressure.

A previously unknown thermodynamically stable high-pressure phase of BeF2 has been predicted using the evolutionary algorithm USPEX. This phase occurs in the pressure range 18-27 GPa. Its structure has C2/c space group symmetry and contains 18 atoms in the primitive unit cell. Given the analogy between BeF2 and SiO2, silica phases have been investigated as well, but the new phase has not been observed to be thermodynamically stable for this system. However, it is found to be metastable and to have comparable energy to the known metastable phases of SiO2, suggesting a possibility of its synthesis.

The phase diagram and hardness of carbon nitrides.

Novel superhard materials, especially those with superior thermal and chemical stability, are needed to replace diamond. Carbon nitrides (C-N), which are likely to possess these characteristics and have even been expected to be harder than diamond, are excellent candidates. Here we report three new superhard and thermodynamically stable carbon nitride phases. Based on a systematic evolutionary structure searches, we report a complete phase diagram of the C-N system at 0-300 GPa and analyze the hardest metastable structures. Surprisingly, we find that at zero pressure, the earlier proposed graphitic-C3N4 structure () is dynamically unstable, and we find the lowest-energy structure based on s-triazine unit and s-heptazine unit.

Novel hydrogen hydrate structures under pressure.

Gas hydrates are systems of prime importance. In particular, hydrogen hydrates are potential materials of icy satellites and comets, and may be used for hydrogen storage. We explore the H2O-H2 system at pressures in the range 0-100 GPa with ab initio variable-composition evolutionary simulations. According to our calculation and previous experiments, the H2O-H2 system undergoes a series of transformations with pressure, and adopts the known open-network clathrate structures (sII, C0), dense "filled ice" structures (C1, C2) and two novel hydrate phases. One of these is based on the hexagonal ice framework and has the same H2O:H2 ratio (2:1) as the C0 phase at low pressures and similar enthalpy (we name this phase Ih-C0). The other newly predicted hydrate phase has a 1:2 H2O:H2 ratio and structure based on cubic ice. This phase (which we name C3) is predicted to be thermodynamically stable above 38 GPa when including van der Waals interactions and zero-point vibrational energy, and explains previously mysterious experimental X-ray diffraction and Raman measurements. This is the hydrogen-richest hydrate and this phase has a remarkable gravimetric density (18 wt.%) of easily extractable hydrogen.

Pressure-induced stabilization and insulator-superconductor transition of BH.

Diborane (B(2)H(6)), a high energy density material, was believed to be stable in a wide P, T interval. A systematic investigation of the B-H system using the ab initio variable-composition evolutionary simulations shows that boron monohydride (BH) is thermodynamically stable and can coexist with solid B, H(2), and B(2)H(6) in a wide pressure range above 50 GPa. B(2)H(6) becomes unstable and decomposes into the Ibam phase of BH and H(2) (C2/c) at 153 GPa. The semiconducting layered Ibam structure of BH at 168 GPa transforms into a metallic phase with space group P6/mmm and a 3D topology with strong B-B and B-H covalent bonds. The Ibam-P6/mmm transformation pathway suggests the possibility of obtaining the metastable Pbcm phase on cold decompression of the P6/mmm phase. The electron-phonon coupling calculations indicate that P6/mmm-BH is a phonon-mediated superconductor with a critical temperature of superconductivity (T(c)) of 14.1-21.4 K at 175 GPa.

An ab initio study on the transition paths from graphite to diamond under pressure.

We calculate and compare the transition paths from graphite to two types of diamond using the variable cell nudged elastic band method. For the phase transition from graphite to cubic diamond, we analyze in detail how the π bonds transit to the σ bonds in an electronic structure. Meanwhile, a new transition path with a lower energy barrier for the transformation from graphite to hexagonal diamond is discovered. The path has its own peculiar sp(2)-sp(3) bonding configurations, serving as a transition state. Further calculation suggests that the sp(2)-sp(3) transition state represents an expected general phenomenon for cold-compressed graphite.

First-principles determination of the structure of magnesium borohydride.

The energy landscape of Mg(BH(4))(2) under pressure is explored by ab initio evolutionary calculations. Two new tetragonal structures, with space groups P4 and I4(1)/acd, are predicted to be lower in enthalpy by 15.4 and 21.2 kJ/mol, respectively, than the earlier proposed P4(2)nm phase. We have simulated x-ray diffraction spectra, lattice dynamics, and equations of state of these phases. The density, volume contraction, bulk modulus, and simulated x-ray diffraction patterns of I4(1)/acd and P4 structures are in excellent agreement with the experimental results.