Haiyang Niu

Modeling hardness of polycrystalline materials and bulk metallic glasses

Though extensively studied, hardness, defined as the resist ance of a material to deformation, still remains a challengi ng issue for a formal theoretical description due to its inherent mechani c l complexity. The widely applied Teter’s empirical corre lation between hardness and shear modulus has been considered to be not alwa ys valid for a large variety of materials. Here, inspired by t he classical work on Pugh’s modulus ratio, we develop a theoret ical model which establishes a robust correlation between h ardness and elasticity for a wide class of materials, including bulk metallic glasses, with results in very good agreement with e xperiment. The simplified form of our model also provides an unambiguous theoretical evidence for Teter’s empirical correlation.

Families of superhard crystalline carbon allotropes constructed via cold compression of graphite and nanotubes.

We report a general scheme to systematically construct two classes of structural families of superhard sp(3) carbon allotropes of cold-compressed graphite through the topological analysis of odd 5+7 or even 4+8 membered carbon rings stemmed from the stacking of zigzag and armchair chains. Our results show that the previously proposed M, bct-C(4), W and Z allotropes belong to our currently proposed families and that depending on the topological arrangement of the native carbon rings numerous other members are found that can help us understand the structural phase transformation of cold-compressed graphite and carbon nanotubes (CNTs). In particular, we predict the existence of two simple allotropes, R and P carbon, which match well the experimental x-ray diffraction patterns of cold-compressed graphite and CNTs, respectively, display a transparent wide-gap insulator ground state and possess a large Vickers hardness comparable to diamond.

Hardness of T-carbon: Density functional theory calculations

We reconsider and interpret the mechanical properties of the recently proposed allotrope of carbon, T-carbon [Sheng et al., Phys. Rev. Lett. 106, 155703 (2011)], using density functional theory in combination with different empirical hardness models. In contrast with the early estimation based on Gao et al.s model, which attributes to T-carbon a high Vickers hardness of 61 GPa comparable to that of superhard cubic boron nitride (c-BN), we find that T-carbon is not a superhard material, since its Vickers hardness does not exceed 10 GPa. Besides providing clear evidence for the absence of superhardness in T-carbon, we discuss the physical reasons behind the failure of Gao et al.s and Simunek and Vackars (SV) models in predicting the hardness of T-carbon, residing in their improper treatment of the highly anisotropic distribution of quasi-sp(3)-like C-C hybrids. A possible remedy for the Gao et al. and SV models based on the concept of the superatom is suggested, which indeed yields a Vickers hardness of about 8 GPa.

Extra-electron induced covalent strengthening and generalization of intrinsic ductile-to-brittle criterion

Traditional strengthening ways, such as strain, precipitation, and solid-solution, come into effect by pinning the motion of dislocation. Here, through first-principles calculations we report on an extra-electron induced covalent strengthening mechanism, which alters chemical bonding upon the introduction of extra-valence electrons in the matrix of parent materials. It is responsible for the brittle and high-strength properties of Al12W-type compounds featured by the typical fivefold icosahedral cages, which are common for quasicrystals and bulk metallic glasses (BMGs). In combination with this mechanism, we generalize ductile-to-brittle criterion in a universal hyperbolic form by integrating the classical Pettifors Cauchy pressure with Pughs modulus ratio for a wide variety of materials with cubic lattices. This study provides compelling evidence to correlate Pughs modulus ratio with hardness of materials and may have implication for understanding the intrinsic brittleness of quasicrystals and BMGs.

Interstitial-boron solution strengthened WB3+x

By means of variable-composition evolutionary algorithm coupled with density functional theory and in combination with aberration-corrected high-resolution transmission electron microscopy experiments, we have studied and characterized the composition, structure, and hardness properties of WB3+x (x < 0.5). We provide robust evidence for the occurrence of stoichiometric WB3 and non-stoichiometric WB3+x, both crystallizing in the metastable hP16 (P63/mmc) structure. No signs for the formation of the highly debated WB4 (both hP20 and hP10) phases were found. Our results rationalize the seemingly contradictory high-pressure experimental findings and suggest that the interstitial boron atom is located in the tungsten layer and vertically interconnect with four boron atoms, thus forming a typical three-center boron net with the upper and lower boron layers in a three-dimensional covalent network, which thereby strengthen the hardness.

Electronic, optical, and mechanical properties of superhard cold-compressed phases of carbon

By means of standard and hybrid density functional theory, we analyzed the electronic, optical, and mechanical properties of the two discovered superhard orthorhombic (W) and monoclinic (M) phases of carbon, synthesized by cold compression. We demonstrated that both phases exhibit a transparent insulating behaviour with indirect band gaps of about 5.4 eV (W) and 4.5 eV (M), and highly isotropic optical spectra, substantially different to those of the related body centered tetragonal C4 phase. The analysis of the elastic constants and Vickers hardness confirmed that these phases are as hard as the second hardest material c-BC2N.

Prediction of novel stable compounds in the Mg-Si-O system under exoplanet pressures

The Mg-Si-O system is the major Earth and rocky planet-forming system. Here, through quantum variable-composition evolutionary structure explorations, we have discovered several unexpected stable binary and ternary compounds in the Mg-Si-O system. Besides the well-known SiO2 phases, we have found two extraordinary silicon oxides, SiO3 and SiO, which become stable at pressures above 0.51 TPa and 1.89 TPa, respectively. In the Mg-O system, we have found one new compound, MgO3, which becomes stable at 0.89 TPa. We find that not only the (MgO)x·(SiO2)y compounds, but also two (MgO3)x·(SiO3)y compounds, MgSi3O12 and MgSiO6, have stability fields above 2.41 TPa and 2.95 TPa, respectively. The highly oxidized MgSi3O12 can form in deep mantles of mega-Earths with masses above 20 M⊕ (M⊕:Earth’s mass). Furthermore, the dissociation pathways of pPv-MgSiO3 are also clarified, and found to be different at low and high temperatures. The low-temperature pathway is MgSiO3 ⇒ Mg2SiO4 + MgSi2O5 ⇒ SiO2 + Mg2SiO4 ⇒ MgO + SiO2, while the high-temperature pathway is MgSiO3 ⇒ Mg2SiO4 + MgSi2O5 ⇒ MgO + MgSi2O5 ⇒ MgO + SiO2. Present results are relevant for models of the internal structure of giant exoplanets, and for understanding the high-pressure behavior of materials.

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.

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.

Superconductivity and unexpected chemistry of germanium hydrides under pressure

Following the idea that hydrogen-rich compounds might be high-T