Source author record

Hua Y. Geng

Hua Y. Geng appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Evidence for mechanical softening-hardening dual anomaly in transition metals from shock compressed vanadium

Solid usually becomes harder and tougher under compression, and turns softer at elevated temperature. Recently, compression-induced softening and heating-induced hardening (CISHIH) dual anomaly was predicted in group VB elements such as vanadium. Here, the evidence for this counterintuitive phenomenon is reported. By using accurate high-temperature high-pressure sound velocities measured at Hugoniot states generated by shock-waves, together with first-principles calculations, we observe not only the prominent compression-induced sound velocity reduction, but also strong heating-induced sound velocity enhancement, in shocked vanadium. The former corresponds to the softening in shear modulus by compression, whereas the latter reflects the reverse hardening by heat. These experiments also unveil another anomaly in Young's modulus that wasn't reported before. Based on the experimental and theoretical data, we infer that vanadium might transition from BCC into two different rhombohedral (RH1 and RH2) phases at about 79GPa and 116GPa along the Hugoniot, respectively, which implies a dramatic difference in static and dynamic loading, as well as the significance of deviatoric stress and rate-relevant effects in high-pressure phase transition dynamics.

preprint2020arXiv

Orbital localization error of density functional theory in shear properties of vanadium and niobium

It is believed that the density functional theory (DFT) describes most elements with s, p and d orbitals very well, except some materials that having strongly localized and correlated valence electrons. In this work, we find that the widely employed exchange-correlation (xc) functionals, including LDA, GGA and meta-GGA, underestimate the shear modulus and phase stability of V and Nb greatly. The advanced hybrid functional that is usually better for correlated system, on the other hand, completely fails in these two simple metals. This striking failure is revealed due to the orbital localization error in GGA, which is further deteriorated by hybrid functionals. This observation is corroborated by a similar failure of DFT+U and van der Waals functionals when applied to V and Nb. To remedy this problem, an semi-empirical approach of DFT+J is proposed which can delocalize electrons by facilitating the on-site exchange. Furthermore, it is observed that including density derivatives slightly improves the performance of the semi-local functionals, with meta-GGA outperforms GGA, and the latter is better than LDA. This discovery indicates the possibility and necessity to include higher-order density derivatives beyond the Laplacian level for the purpose to remove the orbital localization error (mainly from d orbitals) and delocalization error (mainly from s and p orbitals) completely in V and Nb, so that to achieve a better description of their electronic structures. The same strategy can be applied to other d electron system and f electron system.

preprint2020arXiv

Prediction of the Reactivity of Argon with Xenon under High Pressure

Pressure significantly modifies the microscopic interactions in condense phase, leading to new patterns of bonding and unconventional chemistry. Both argon and xenon possess closed-shell electronic structures, which renders them chemically unreactive. Using unbiased structure searching techniques combined with first-principles calculations, we demonstrate the reaction of argon with xenon at a pressure as low as 1.1 GPa, producing a novel van der Waals compound XeAr2, which crystallizes in the MgCu2-type Laves phase structure. Due to the pressure-induced delocalization of the electrons in outermost shells, the covalent Xe-Xe and Xe-Ar bonds have been detected which lead XeAr2 to be unexpectedly stable without any phase transition or decomposition at least up to 500 GPa.

preprint2020arXiv

Thermodynamic anomalies and three distinct liquid-liquid transitions in warm dense liquid hydrogen

The properties of hydrogen at high pressure have wide implications in astrophysics and high-pressure physics. Its phase change in the liquid is variously described as a metallization, H2-dissociation, density discontinuity or plasma phase transition. It has been tacitly assumed that these phenomena coincide at a first-order liquid-liquid transition (LLT). In this work, the relevant pressure-temperature conditions are thoroughly explored with first-principles molecular dynamics. We show there is a large dependency on exchange-correlation functional and significant finite size effects. We use hysteresis in a number of measurable quantities to demonstrate a first-order transition up to a critical point, above which molecular and atomic liquids are indistinguishable. At higher temperature beyond the critical point, H2-dissociation becomes a smooth cross-over in the supercritical region that can be modelled by a pseudo-transition, where the H2-2H transformation is localized and does not cause a density discontinuity at metallization. Thermodynamic anomalies and counter-intuitive transport behavior of protons are also discovered even far beyond the critical point, making this dissociative transition highly relevant to the interior dynamics of Jovian planets. Below the critical point, simulation also reveals a dynamic H2-2H chemical equilibrium with rapid interconversion, showing that H2 and H are miscible. The predicted critical temperature lies well below the ionization temperature. Our calculations unequivocally demonstrate that there are three distinct regimes in the liquid-liquid transition of warm dense hydrogen.

preprint2016arXiv

Compression and phase diagram of lithium hydrides at elevated pressures and temperatures by first-principles calculations

High pressure and high temperature properties of AB (A = $^6$Li, $^7$Li; B = H, D, T) are investigated with first-principles method comprehensively. It is found that the H$^{-}$ sublattice features in the low-pressure electronic structure near the Fermi level of LiH are shifted to that dominated by the Li$^{+}$ sublattice in compression. The lattice dynamics is studied in quasi-harmonic approximation, from which the phonon contribution to the free energy and the isotopic effects are accurately modelled with the aid of a parameterized double-Debye model. The obtained equation of state (EOS) matches perfectly with available static experimental data. The calculated principal Hugoniot is also in accordance with that derived from shock wave experiments. Using the calculated principal Hugoniot and the previous theoretical melting curve, we predict a shock melting point at 56 GPa and 1923 K. In order to establish the phase diagram for LiH, the phase boundaries between the B1 and B2 solid phases are explored. The B1-B2-liquid triple point is determined at about 241 GPa and 2413 K. The remarkable shift in the phase boundaries by isotopic effect and temperature reveal the significant role played by lattice vibrations. Furthermore, the Hugoniot of the static-dynamic coupling compression is assessed. Our EOS suggests that a precompression of the sample to 50 GPa will allow the shock Hugoniot passing through the triple point and entering the B2 solid phase. This transition leads to a discontinuity with 4.6% volume collapse, about four times greater than the same B1-B2 transition at zero temperature.

preprint2016arXiv

Lattice stability and high pressure melting mechanism of dense hydrogen up to 1.5 TPa

Lattice stability and metastability, as well as melting, are important features of the physics and chemistry of dense hydrogen. Using ab initio molecular dynamics (AIMD), the classical superheating limit and melting line of metallic hydrogen are investigated up to 1.5 TPa. The computations show that the classical superheating degree is about 100 K, and the classical melting curve becomes flat at a level of 350 K when beyond 500 GPa. This information allows us to estimate the well depth and the potential barriers that must be overcome when the crystal melts. Inclusion of nuclear quantum effects (NQE) using path integral molecular dynamics (PIMD) predicts that both superheating limit and melting temperature are lowered to below room temperature, but the latter never reach absolute zero. Detailed analysis indicates that the melting is thermally activated, rather than driven by pure zero-point motion (ZPM). This argument was further supported by extensive PIMD simulations, demonstrating the stability of Fddd structure against liquefaction at low temperatures.

preprint2016arXiv

Stability of rhombohedral phases in vanadium at high-pressure and high-temperature: first-principles investigations

The pressure-induced transition of vanadium from BCC to rhombohedral structures is unique and intriguing among transition metals. In this work, the stability of these phases is revisited by using density functional theory. At finite temperatures, a novel transition of rhombohedral phases back to BCC phase induced by thermal electrons is discovered. This reentrant transition is found not driven by phonons, instead it is the electronic entropy that stabilizes the latter phase, which is totally out of expectation. Parallel to this transition, we find a peculiar and strong increase of the shear modulus C44 with increasing temperature. It is counter-intuitive in the sense that it suggests an unusual harding mechanism of vanadium by temperature. With these stability analyses, the high-pressure and finite-temperature phase diagram of vanadium is proposed. Furthermore, the dependence of the stability of RH phases on the Fermi energy and chemical environment is investigated. The results demonstrate that the position of the Fermi level has a significant impact on the phase stability, and follows the band-filling argument. Besides the Fermi surface nesting, we find that the localization/delocalization of the d orbitals also contributes to the instability of rhombohedral distortions in vanadium.

preprint2014arXiv

Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model, the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4 fold for a two-level implementation, and can be increased to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of $r_{s}=0.912$.

preprint2013arXiv

Theoretical assessment on the possibility of constraining point defect energetics by pseudo-phase transition pressures

Making use of the energetics and equations of state of defective uranium dioxide that calculated with first-principles method, we demonstrate a possibility of constraining the formation energy of point defects by measuring the transition pressures of the corresponding pseudo-phase of defects. The mechanically stable range of fluorite structure of UO2, which dictates the maximum possible pressure of relevant pseudo-phase transitions, gives rise to defect formation energies that span a wide band and overlap with the existing experimental estimates. We reveal that the knowledge about pseudo-phase boundaries can not only provide important information of energetics that is helpful for reducing the scattering in current estimates, but also be valuable for guiding theoretical assessments, even to validate or disprove a theory. In order to take defect interactions into account and to extrapolate the physical quantities at finite stoichiometry deviations to that near the stoichiometry, we develop a general formalism to describe the thermodynamics of a defective system. We also show that it is possible to include interactions among defects in a simple expression of point defect model (PDM) by introducing an auxiliary constant mean-field. This generalization of the simple PDM leads to great versatility that allows one to study nonlinear effects of stoichiometry deviation on materials' behavior. It is a powerful tool to extract the defect energetics from finite defect concentrations to the dilute limit. Besides these, the full content of the theoretical formalism and some relevant and interesting issues, including reentrant pseudo-transition, multi-defect coexistence, charged defects, and possible consequence of instantaneous defective response in a quantum crystal, are explored and discussed.

preprint2012arXiv

Anomalies in non-stoichiometric uranium dioxide induced by pseudo-phase transition of point defects

A uniform distribution of point defects in an otherwise perfect crystallographic structure usually describes a unique pseudo phase of that state of a non-stoichiometric material. With off-stoichiometric uranium dioxide as a prototype, we show that analogous to a conventional phase transition, these pseudo phases also will transform from one state into another via changing the predominant defect species when external conditions of pressure, temperature, or chemical composition are varied. This exotic transition is numerically observed along shock Hugoniots and isothermal compression curves in UO2 with first-principles calculations. At low temperatures, it leads to anomalies (or quasi-discontinuities) in thermodynamic properties and electronic structures. In particular, the anomaly is pronounced in both shock temperature and the specific heat at constant pressure. With increasing of the temperature, however, it transforms gradually to a smooth cross-over, and becomes less discernible. The underlying physical mechanism and characteristics of this type of transition are encoded in the Gibbs free energy, and are elucidated clearly by analyzing the correlation with the variation of defect populations as a function of pressure and temperature. The opportunities and challenges for a possible experimental observation of this phase change are also discussed.

preprint2012arXiv

High-pressure behavior of dense hydrogen up to 3.5 TPa from density functional theory calculations

Structural behavior and equation of state of atomic and molecular crystal phases of dense hydrogen at pressures up to 3.5 TPa are systematically investigated with density functional theory. The results indicate that the Vinet EOS model that fitted to low-pressure experimental data overestimates the compressibility of dense hydrogen drastically when beyond 500 GPa. Metastable multi-atomic molecular phases with weak covalent bonds are observed. When compressed beyond about 2.8 TPa, these exotic low-coordinated phases become competitive with the groundstate and other high-symmetry atomic phases. Using nudged elastic band method, the transition path and the associated energy barrier between these high-pressure phases are evaluated. In particular for the case of dissociation of diatomic molecular phase into the atomic metallic Cs-IV phase, the existent barrier might raise the transition pressure about 200 GPa at low temperatures. Plenty of flat and broad basins on the energy surface of dense hydrogen have been discovered, which should take a major responsibility for the highly anharmonic zero point vibrations of the lattice, as well as the quantum structure fluctuations in some extreme cases. At zero pressure, our analysis demonstrates that all of these atomic phases of dense hydrogen known so far are unquenchable. NOTE: In the previous version of this paper (1010.3392v1) we made a mistake when evaluating the enthalpy of Cs-IV phase, which misled us to a conclusion that taking the multi-atomic molecular phases as the ground-state. After corrected this error, however, the atomic phase of Cs-IV becomes the static structure with the lowest enthalpy. Current version not only includes a substantial improvement of the previous one, but also contains many NEW interesting topics that were not touched before.

preprint2011arXiv

First-principles study on oxidation effects in uranium oxides and high-pressure high-temperature behavior of point defects in uranium dioxide

Formation Gibbs free energy of point defects and oxygen clusters in uranium dioxide at high-pressure high-temperature conditions are calculated from first principles, using the LSDA+U approach for the electronic structure and the Debye model for the lattice vibrations. The phonon contribution on Frenkel pairs is found to be notable, whereas it is negligible for the Schottky defect. Hydrostatic compression changes the formation energies drastically, making defect concentrations depend more sensitively on pressure. Calculations show that, if no oxygen clusters are considered, uranium vacancy becomes predominant in overstoichiometric UO2 with the aid of the contribution from lattice vibrations, while compression favors oxygen defects and suppresses uranium vacancy greatly. At ambient pressure, however, the experimental observation of predominant oxygen defects in this regime can be reproduced only in a form of cuboctahedral clusters, underlining the importance of defect clustering in UO2+x. Making use of the point defect model, an equation of state for non-stoichiometric oxides is established, which is then applied to describe the shock Hugoniot of UO2+x. Furthermore, the oxidization and compression behavior of uranium monoxide, triuranium octoxide, uranium trioxide, and a series of defective UO2 at zero Kelvin are investigated. The evolution of mechanical properties and electronic structures with an increase of the oxidation degree are analyzed, revealing the transition of the groundstate of uranium oxides from metallic to Mott insulator and then to charge-transfer insulator due to the interplay of strongly correlated effects of 5f orbitals and the shift of electrons from uranium to oxygen atoms.

preprint2010arXiv

Interplay of defect cluster and the stability of xenon in uranium dioxide by density functional calculations

Self-defect clusters in bulk matrix might affect the thermodynamic behavior of fission gases in nuclear fuel such as uranium dioxide. With first-principles LSDA+U calculations and taking xenon as a prototype, we find that the influence of oxygen defect clusters on the thermodynamics of gas atoms is prominent, which increases the solution energy of xenon by a magnitude of 0.5 eV, about 43% of the energy difference between the two lowest lying states at 700 K. Calculation also reveals a thermodynamic competition between the uranium vacancy and tri-vacancy sites to incorporate xenon in hyper-stoichiometric regime at high temperatures. The results show that in hypo-stoichiometric regime neutral tri-vacancy sites are the most favored position for diluted xenon gas, whereas in hyper-stoichiometric condition they prefer to uranium vacancies even after taking oxygen self-defect clusters into account at low temperatures, which not only confirms previous studies but also extends the conclusion to more realistic fuel operating conditions. The observation that gas atoms are ionized to a charge state of Xe+ when at a uranium vacancy site due to strong Madelung potential implies that one can control temperature to tune the preferred site of gas atoms and then the bubble growth rate. A solution to the notorious meta-stable states difficulty that frequently encountered in DFT+U applications, namely, the quasi-annealing procedure, is also discussed.