Guang-Rui Qian

Ab initio study of the formation of transparent carbon under pressure

A body-centered tetragonal carbon (bct-Carbon) allotrope has been predicted to be a transparent carbon polymorph obtained under pressure. The structural transition pathways from graphite to diamond, M-Carbon, and bct-Carbon are simulated and the lowest activation barrier is found for the graphite-bct transition. Furthermore, bct-Carbon has higher shear strength than diamond due to its perpendicular graphene-like structure. Our results provide a possible explanation for the formation of a transparent carbon allotrope via the cold compression of graphite. We also verify that this allotrope is hard enough to crack diamond.

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.

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.

Variable cell nudged elastic band method for studying solid–solid structural phase transitions

a b s t r a c t The nudged elastic band (NEB) method, as a popular technique for studying reaction paths due to its efficiency, has not been extensively used in solid state physics because of the need to deal with the variation of the unit cell during solid-solid transformations. Here we present an extended technique— a variable-cell NEB (VC-NEB) technique combined with the ab initio method, implemented in the USPEX code. Our technique is applied to reconstructive solid-solid phase transitions of GaN: from wurtzite to rocksalt (B4 → B1) and from zincblende to rocksalt (B3 → B1). As a more challenging application, we study the mechanism of the recently predicted insulator-metal phase transition of BH. The results reveal that the VC-NEB technique is an efficient and general method and should have wide applications for studying the paths and mechanisms of reconstructive phase transitions. Results of the VC-NEB method can be considered as a starting point for a more sophisticated treatment using the transition path sampling method, the main prerequisite of which is to have an initial transformation trajectory.

Chalcopyrite polymorph for superhard BC2N

A chalcopyrite structure is predicted for the potential superhard boron-carbonitride (BC2N), based on ab initio pseudopotential density functional method. Lattice constant, formation energy, bulk modulus, band structure, and electron density of states are calculated. The results show that this compound is a wide gap semiconductor with a direct band gap of about 3.3eV. The calculated Vickers hardness is about 72GPa, which is very close to the measured hardness of c-BC2N (76±4GPa).

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.

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.

Antiferromagnetic Stabilization in the Ti8O12 Cluster

Using the evolutionary algorithm USPEX and DFT+U calculations, we predicted a high-symmetry geometric structure of the bare Ti8 O12 cluster composed of 8 Ti atoms forming a cube, in which O atoms are at midpoints of all of its edges, in excellent agreement with experimental results. Using natural bond orbital analysis, adaptive natural density partitioning algorithm, electron localization function, and partial charge plots, we find the origin of the particular stability of bare Ti8 O12 cluster: unique chemical bonding where eight electrons of Ti atoms interacting with each other in antiferromagnetic fashion to lower the total energy of the system. The bare Ti8 O12 is thus an unusual molecule stabilized by d-orbital antiferromagnetic coupling.

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.

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.