Source author record

Hanyu Liu

Hanyu Liu 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

29works
12topics
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

29 published item(s)

preprint2026arXiv

Beyond Autoregressive RTG: Conditioning via Injection Outside Sequential Modeling in Decision Transformer

Decision Transformer (DT) formulates offline reinforcement learning as autoregressive sequence modeling, achieving promising results by predicting actions from a sequence of Return-to-Go (RTG), state, and action tokens. However, RTG is a scalar that summarizes future rewards, containing far less information than typical state or action vectors, yet it consumes the same computational budget per token. Worse, the self-attention cost of Transformers grows quadratically with sequence length, so including RTG as a separate token adds unnecessary overhead. We propose SlimDT, which removes RTG from the autoregressive sequence. Instead, we inject RTG information into the state representations before the sequential modeling step, allowing the Transformer to process only a compact (state, action) sequence. This reduces the sequence length by one-third, directly improving inference efficiency. On the D4RL benchmark, SlimDT surpasses standard DT across various tasks and achieves performance comparable to existing state-of-the-art methods. Decoupling a sparse conditioning signal from an information-rich sequence thus yields both computational gains and higher task performance.

preprint2026arXiv

Isotropic Superconductivity in Room-temperature Superconductor LaSc$_{2}$H$_{24}$

The discovery of LaSc$_{2}$H$_{24}$ represents a milestone in the quest for room-temperature superconductivity, yet the microscopic mechanism underlying its superior performance remains unclear. Through a comprehensive revisit of theoretical calculations, we uncover a pivotal transition from the anisotropic two-gap superconductivity of LaH$_{10}$ to the isotropic single-gap superconductivity in LaSc$_{2}$H$_{24}$ upon the introduction of scandium, thereby enhancing the superconducting critical temperature ($T_\mathrm{c}$). This enhancement is rooted in a critical dual role of Sc $3d$ electrons: i) the Sc-derived Jahn-Teller effect promotes hydrogen metallization via the elongation of specific interlayer H-H bonds and enhances electron-phonon coupling (EPC) through the softening of associated phonon modes; ii) Sc $3d$ electrons reconstruct the electronic structure into an MgB$_{2}$-like configuration, generating novel Sc-H-Sc $σ$- and $π$-bonding states with EPC strengths comparable to LaH$_{10}$. Crucially, the pronounced hybridization between Sc and the hydrogen cages effectively unifies these two contributions on the Fermi surface. This Sc-induced gap unification bridges the high-EPC H-H states with widespread Sc-H states, establishing an isotropic single-gap nature with a large overall EPC strength. Our findings identify this Sc-induced gap unification as the fundamental mechanism for achieving room-temperature superconductivity in LaSc$_{2}$H$_{24}$, offering a theoretical blueprint for the future design of superior superconducting hydrides.

preprint2026arXiv

MiniCPM-o 4.5: Towards Real-Time Full-Duplex Omni-Modal Interaction

Recent progress in multimodal large language models (MLLMs) has brought AI capabilities from static offline data processing to real-time streaming interaction, yet they still remain far from human-level multimodal interaction. The key bottlenecks are no longer modality coverage or latency alone, but the interaction paradigm itself. First, perception and response are still separated into alternating phases, preventing models from incorporating new inputs for timely adjustment during generation. Second, most current models remain reactive, responding only to explicit user requests instead of acting proactively in the evolving multimodal environment. We present MiniCPM-o 4.5, our latest effort towards human-like multimodal interaction, which mitigates these gaps by real-time full-duplex omni-modal interaction. It can see, listen, and speak simultaneously in real-time, while also exhibiting proactive behaviors such as issuing reminders or comments based on its continuous understanding of the live scene. The key technique behind MiniCPM-o 4.5 is Omni-Flow, a unified streaming framework that aligns omni-modal inputs and outputs along a shared temporal axis. This formulation converts conventional turn-based interaction into a full-duplex, time-aligned process, enabling simultaneous perception and response and allowing proactive behavior to arise within the same framework. With a total of 9B parameters, MiniCPM-o 4.5 approaches Gemini 2.5 Flash in vision-language capabilities, delivering state-of-the-art open-source performance at its scale. It also surpasses Qwen3-Omni-30B-A3B in omni-modal understanding and delivers better speech generation, with significantly higher computation efficiency. Driven by its efficient architecture design and inference optimization, the model can perform real-time full-duplex omni-modal interaction on edge devices with less than 12GB RAM cost.

preprint2026arXiv

The Impact of Ionic Anharmonicity on Superconductivity in Metal-Stuffed B-C Clathrates

Metal-stuffed B$-$C compounds with sodalite clathrate structure have captured increasing attention due to their predicted exceptional superconductivity above liquid nitrogen temperature at ambient pressure. However, by neglecting the quantum lattice anharmonicity, the existing studies may result in an incomplete understanding of such a lightweight system. Here, using state-of-the-art ab initio methods incorporating quantum effects and machine learning potentials, we revisit the properties of a series of $XY$$\text{B}_{6}\text{C}_{6}$ clathrates where $X$ and $Y$ are metals. Our findings show that ionic quantum and anharmonic effects can harden the $E_g$ and $E_u$ vibrational modes, enabling the dynamical stability of 15 materials previously considered unstable in the harmonic approximation, including materials with previously unreported ($XY$)$^{1+}$ state, which is demonstrated here to be crucial to reach high critical temperatures. Further calculations based on the anisotropic Migdal-Eliashberg equation demonstrate that the $T_\text{c}$ values for KRb$\text{B}_{6}\text{C}_{6}$ and Rb$\text{B}_{3}\text{C}_{3}$ among these stabilized compounds are 102 and 115 K at 0 and 15 GPa, respectively, both being higher than $T_\text{c}$ of 92 K of KPb$\text{B}_{6}\text{C}_{6}$ at the anharmonic level. These record-high $T_\text{c}$ values, surpassing liquid nitrogen temperatures, emphasize the importance of anharmonic effects in stabilizing B-C clathrates with large electron-phonon coupling strength and advancing the search for high-$T_\text{c}$ superconductivity at (near) ambient pressure.

preprint2023arXiv

Pressure-Induced Superconductivity in Topological Heterostructure (PbSe)5(Bi2Se3)6

Recently, the natural heterostructure of (PbSe)5(Bi2Se3)6 has been theoretically predicted and experimentally confirmed as a topological insulator. In this work, we induce superconductivity in (PbSe)5(Bi2Se3)6 by implementing high pressure. As increasing pressure up to 10 GPa, superconductivity with Tc ~ 4.6 K suddenly appears, followed by an abrupt decrease. Remarkably, upon further compression above 30 GPa, a new superconducting state arises, where pressure raises the Tc to an unsaturated 6.0 K within the limit of our research. Combining XRD and Raman spectroscopies, we suggest that the emergence of two distinct superconducting states occurs concurrently with the pressure-induced structural transition in this topological heterostructure (PbSe)5(Bi2Se3)6.

preprint2022arXiv

A Symmetry-orientated Divide-and-Conquer Method for Crystal Structure Prediction

Crystal structure prediction has been a subject of topical interest, but remains a substantial challenge, especially for complex structures as it deals with the global minimization of the extremely rugged high-dimensional potential energy surface. In this manuscript, a symmetry-orientated divide-and-conquer scheme was proposed to construct a symmetry tree graph, where the entire search space is decomposed into a finite number of symmetry-dependent subspaces. An artificial intelligence-based symmetry selection strategy was subsequently devised to select the low-lying subspaces with high symmetries for global exploration and in-depth exploitation. Our approach can significantly simplify the problem of crystal structure prediction by avoiding exploration of the most complex P1 subspace on the entire search space and have the advantage of preserving the crystal symmetry during structure evolution, making it well suitable for predicting the complex crystal structures. The effectiveness of the method has been validated by successful prediction of the candidate structures of binary Lennard-Jones mixtures and high-pressure phase of ice, containing more than one hundred atoms in the simulation cell. The work, therefore, opens up an opportunity towards achieving the long-sought goal for crystal structure prediction of complex systems.

preprint2022arXiv

Experimental phase transition mapping for hydrogen above 300 K up to 300 GPa

A vast amount of Raman spectroscopic data, obtained using diamond anvil cell technique under in situ high pressure-temperature conditions, has been used for mapping the phase transitions of hydrogen in the temperature range of 300 K to 900 K and pressure up to 300 gigapascal, a pressure-temperature region in which most phase information was unknown previously. The isothermal and isobaric dependencies in frequency and peak width for intramolecular vibration Raman mode of hydrogen are established based on thousands of data files which are measured during hundreds of independent high pressure-temperature experimental runs. Discontinuities of pressure and temperature dependencies in frequency and peak-width, together with appearance or disappearance of observed Raman modes were obtained. Those transition data show self-consistency and have been used for outlining 13 phase transition boundaries and 15 possible phases. This large, diversified behavior of spectra demonstrates that surprisingly rich and complicated phase transitions may exist in this pressure-temperature region.

preprint2022arXiv

Novel Valence Transition in Elemental Metal Europium around 80 GPa

Valence transition could induce structural, insulator-metal, nonmagnetic-magnetic and superconducting transitions in rare-earth metals and compounds, while the underlying physics remains unclear due to the complex interaction of localized 4f electrons as well as their coupling with itinerant electrons. The valence transition in the elemental metal europium (Eu) still has remained as a matter of debate. Using resonant x-ray emission scattering and x-ray diffraction, we pressurize the states of 4f electrons in Eu and study its valence and structure transitions up to 160 GPa. We provide compelling evidence for a valence transition around 80 GPa, which coincides with a structural transition from a monoclinic (C2/c) to an orthorhombic phase (Pnma). We show that the valence transition occurs when the pressure-dependent energy gap between 4f and 5d electrons approaches the Coulomb interaction. Our discovery is critical for understanding the electrodynamics of Eu, including magnetism and high-pressure superconductivity.

preprint2021arXiv

En route to high Tc superconductivity via Rb substitution of guest metal atoms in SrB3C3 clathrate

Recently, a host/guest clathrate SrB3C3 with sp3-bonded boron-carbon framework was synthesized at around 50 GPa. On the basis of electron count, the structure is understood as guest Sr2+ cations intercalated in the (B3C3)3- framework. Previous calculations suggest that SrB3C3 is a hole conductor with an estimated superconducting critical temperature (Tc) of 42 K at ambient pressure. If atoms with similar radius, such as Rb, can substitute Sr2+ in the lattice, the electronic as well as superconductivity properties of this material will be modified significantly. Here, we perform extensive simulations on the stability and physical properties of Rb-Sr-B3C3 system using first-principles density functional calculation in combination with cluster expansion and CALYPSO structure prediction method. We predict a phonon-mediated superconductor Rb0.5Sr0.5B3C3 with a remarkably high Tc of 78 K at ambient pressure, which is a significant improvement from the estimated value (42 K) in SrB3C3. The current results suggest that substitution of alkali atom in synthesized clathrate SrB3C3 is a viable route toward high-Tc compounds.

preprint2021arXiv

Pressure-induced Superconductivity in dual-topological semimetal Pt2HgSe3

Recently monolayer jacutingaite (Pt2HgSe3), a naturally occurring exfoliable mineral, discovered in Brazil in 2008, has been theoretically predicted as a candidate quantum spin Hall system with a 0.5 eV band gap, while the bulk form is one of only a few known dual-topological insulators which may host different surface states protected by symmetries. In this work, we systematically investigate both structure and electronic evolution of bulk Pt2HgSe3 under high pressure up to 96 GPa. The nontrivial topology persists up to the structural phase transition observed in the high-pressure regime. Interestingly, we found that this phase transition is accompanied by the appearance of superconductivity at around 55 GPa and the critical transition temperature Tc increases with applied pressure. Our results demonstrate that Pt2HgSe3 with nontrivial topology of electronic states displays new ground states upon compression and raises potentials in application to the next-generation spintronic devices.

preprint2020arXiv

Emotion Recognition for In-the-wild Videos

This paper is a brief introduction to our submission to the seven basic expression classification track of Affective Behavior Analysis in-the-wild Competition held in conjunction with the IEEE International Conference on Automatic Face and Gesture Recognition (FG) 2020. Our method combines Deep Residual Network (ResNet) and Bidirectional Long Short-Term Memory Network (BLSTM), achieving 64.3% accuracy and 43.4% final metric on the validation set.

preprint2020arXiv

Experimental Syntheses of Sodalite-like Clathrate EuH$_6$ and EuH$_9$ at Extreme Pressures

The recent discovery of a class of sodalite-like clathrate superhydrides (e.g., YH6, YH9, ThH9, ThH10, and LaH10) at extreme pressures, which exhibit commonly a high-temperature superconductivity with the highest Tc approaching 260 K for LaH10, opened up a new era in search of high-temperature superconductors in metal superhydrides. There is a high interest towards the finding of alternative clathrate superhydrides that might witness the long-dreamed room-temperature superconductivity. Here, we target on the experimental synthesis of strongly-correlated europium (Eu) superhydrides where theory can fail for the prediction of superconductivity. We pressurized and laser-heated the mixture of metal Eu and ammonia borane (NH3BH3) in a diamond anvil cell and successfully synthesized the sodalite-like clathrate EuH6 and EuH9 at conditions of 152 GPa and 1,700 K, and 170 GPa and 2,800 K, respectively. Two non-clathrate structured phases of EuH5 and EuH6 were also synthesized that are not reported in lanthanide superhydrides. Calculated large H-derived electronic density of states at the Fermi level in clathrate EuH6 implies the potential of high temperature superconductivity. Our work created a model superhydride platform for subsequent investigation on how strongly-correlated effect in electronic structure can affect the superconductivity of superhydrides, a phenomenon that is not known thus far.

preprint2015arXiv

Dissociation products and structures of solid H2S at strong compression

Hydrogen sulfides have recently received a great deal of interest due to the record high superconducting temperatures of up to 203 K observed on strong compression of dihydrogen sulfide (H2S). A joint theoretical and experimental study is presented in which decomposition products and structures of compressed H2S are characterized, and their superconducting properties are calculated. In addition to the experimentally known H2S and H3S phases, our first-principles structure searches have identified several energetically competitive stoichiometries that have not been reported previously; H2S3, H3S2, and H4S3. In particular, H4S3 is predicted to be thermodynamically stable within a large pressure range of 25-113 GPa. High-pressure room-temperature X-ray diffraction measurements confirm the presence of H3S and H4S3 through decomposition of H2S that emerge at 27 GPa and coexist with residual H2S, at least up to the highest pressure studied in our experiments of 140 GPa. Electron-phonon coupling calculations show that H4S3 has a small Tc of below 2 K, and that H2S is mainly responsible for the observed superconductivity of samples prepared at low temperature (<100K).

preprint2015arXiv

High-energy-density and superhard nitrogen-rich B-N compounds

The pressure-induced transformation of diatomic nitrogen into non-molecular polymeric phases may produce potentially useful high-energy-density materials. We combine first-principles calculations with structure searching to predict a new class of nitrogen-rich boron nitrides with a stoichiometry of B3N5 that are stable or metastable relative to solid N2 and h-BN at ambient pressure. The most stable phase at ambient pressure has a layered structure (h-B3N5) containing hexagonal B3N3 layers sandwiched with intercalated freely rotating N2 molecules. At 15 GPa, a three-dimensional C2221 structure with single N-N bonds becomes the most stable. This pressure is much lower than that required for triple-to-single bond transformation in pure solid nitrogen (110 GPa). More importantly, C2221-B3N5 is metastable, and can be recovered under ambient conditions. Its energy density of 3.44 kJ/g makes it a potential high-energy-density material. In addition, stress-strain calculations estimate a Vickers hardness of 44 GPa. Structure searching reveals a new clathrate sodalite-like BN structure that is metastable under ambient conditions.

preprint2015arXiv

Hydrogen sulphide at high pressure: a strongly-anharmonic phonon-mediated superconductor

We use first principles calculations to study structural, vibrational and superconducting properties of H$_2$S at pressures $P\ge 200$ GPa. The inclusion of zero point energy leads to two different possible dissociations of H$_2$S, namely 3H$_2$S $\to$ 2H$_3$S + S and 5H$_2$S $\to$ 3H$_3$S + HS$_2$, where both H$_3$S and HS$_2$ are metallic. For H$_3$S, we perform non-perturbative calculations of anharmonic effects within the self-consistent harmonic approximation and show that the harmonic approximation strongly overestimates the electron-phonon interaction ($λ\approx 2.64$ at 200 GPa) and T$_c$. Anharmonicity hardens HS bond-stretching modes and softens H--S bond-bending modes. As a result, the electron-phonon coupling is suppressed by $30\%$ ($λ\approx 1.84$ at 200 GPa). Moreover, while at the harmonic level T$_c$ decreases with increasing pressure, the inclusion of anharmonicity leads to a T$_c$ that is almost independent of pressure. High pressure hydrogen sulfide is a strongly anharmonic superconductor.

preprint2015arXiv

N$_2$H: A Novel Polymeric Hydronitrogen as a High Energy Density Material

The polymeric phase of nitrogen connected by the lower (than three) order N-N bonds has been long sought after for the potential application as high energy density materials. Here we report a hitherto unknown polymeric N$_2$H phase discovered in the high-pressure hydronitrogen system by first-principle structure search method based on particle swarm optimization algorithm. This polymeric hydronitrogen consists of quasi-one-dimensional infinite armchair-like polymeric N chains, where H atoms bond with two adjacent N located at one side of armchair edge. It is energetically stable against decomposition above ~33 GPa, and shows novel metallic feature as the result of pressure-enhanced charge transfer and delocalization of π electrons within the infinite nitrogen chains. The high energy density (~4.40 KJ/g), high nitrogen content (96.6%), as well as relatively low stabilization pressure, make it a possible candidate for high energy density applications. It also has lattice dynamical stability down to the ambient pressure, allowing for the possibility of kinetic stability with respect to variations of external conditions. Experimental synthesis of this novel phase is called for.

preprint2015arXiv

Phase Diagram and High-Temperature Superconductivity of Compressed Selenium Hydrides

Recent discovery of high-temperature superconductivity (Tc = 190 K) in sulfur hydrides at megabar pressures breaks the traditional belief on the Tc limit of 40 K for conventional superconductors, and open up the doors in searching new high-temperature superconductors in compounds made up of light elements. Selenium is a sister and isoelectronic element of sulfur, with a larger atomic core and a weaker electronegativity. Whether selenium hydrides share similar high-temperature superconductivity remains elusive, but it is a subject of considerable interest. First-principles swarm structure predictions are performed in an effort to seek for energetically stable and metallic selenium hydrides at high pressures. We find the phase diagram of selenium hydrides is rather different from its sulfur analogy, which is indicates by the emergence of new phases and the change of relative stabilities. Three stable and metallic species with stoichiometries of HSe2, HSe and H3Se are identified above ~120 GPa and they all exhibit superconductive behaviors, of which the hydrogen-rich HSe and H3Se phases show high Tc in the range of 40-110 K. Our simulations established the high-temperature superconductive nature of selenium hydrides and provided useful route for experimental verification.

preprint2015arXiv

Quantum Hydrogen-Bond Symmetrization and High-Temperature Superconductivity in Hydrogen Sulfide

Hydrogen compounds are peculiar as the quantum nature of the proton can crucially affect their structural and physical properties. A remarkable example are the high-pressure phases of H$_2$O, where quantum proton fluctuations favor the symmetrization of the H bond and lower by 30 GPa the boundary between the asymmetric structure and the symmetric one. Here we show that an analogous quantum symmetrization occurs in the recently discovered sulfur hydride superconductor with the record superconducting critical temperature $T_c=203$ K at 155 GPa. In this system, according to classical theory, superconductivity occurs via formation of a structure of stoichiometry H$_3$S with S atoms arranged on a body-centered-cubic (bcc) lattice. For $P \gtrsim 175$ GPa, the H atoms are predicted to sit midway between two S atoms, in a structure with $Im\bar3m$ symmetry. At lower pressures the H atoms move to an off-center position forming a short H$-$S covalent bond and a longer H$\cdots$S hydrogen bond, in a structure with $R3m$ symmetry. X-ray diffraction experiments confirmed the H$_3$S stoichiometry and the S lattice sites, but were unable to discriminate between the two phases. Our present ab initio density-functional theory (DFT) calculations show that the quantum nuclear motion lowers the symmetrization pressure by 72 GPa. Consequently, we predict that the $Im\bar3m$ phase is stable over the whole pressure range within which a high $T_c$ was measured. The observed pressure-dependence of $T_c$ is closely reproduced in our calculations for the $Im\bar3m$ phase, but not for the $R3m$ phase. Thus, the quantum nature of the proton completely rules the superconducting phase diagram of H$_3$S.

preprint2015arXiv

Tellurium Hydrides at High Pressures: High-temperature Superconductors

Observation of high-temperature superconductivity in sulfur hydrides at megabar pressures has generated an irresistible wave on searching for new superconductors in other compressed hydrogen-rich compounds. An immediate effort is towards exploration of the relevant candidate of tellurium hydrides, where tellurium is isoelectronic to sulfur but it has a heavier atomic mass and much weaker electronegativity. The hitherto unknown phase diagram of tellurium hydrides at high pressures was investigated by a first-principles swarm structure search. Four energetically stable and metallic stoichiometries of H4Te, H5Te2, HTe and HTe3 were uncovered above 140 GPa, showing a distinct potential energy map of tellurium hydrides from those in sulfur and selenium hydrides. The two hydrogen-rich H4Te and H5Te2 species adopt ionic structures containing exotic quasi-molecular H2 and linear H3 units, respectively. Strong electron-phonon couplings associated with the intermediate-frequency H-derived wagging and bending modes make them good superconductors with high Tc in the range of 46-104 K.

preprint2014arXiv

Hardness of FeB4: Density functional theory investigation

A recent experimental study reported the successful synthesis of an orthorhombic FeB4 with a high hardness of 62 GPa, which has reignited extensive interests on whether transition metal borides (TRBs) compounds will become superhard materials. However, it is contradicted with some theoretical studies suggesting transition metal boron compounds are unlikely to become superhard materials. Here, we examined structural and electronic properties of FeB4 using density functional theory. The electronic calculations show the good metallicity and covalent FeB bonding. Meanwhile, we extensively investigated stress strain relations of FeB4 under various tensile and shear loading directions. The calculated weakest tensile and shear stresses are 40 GPa and 25 GPa, respectively. Further simulations (e.g. electron localized function and bond length along the weakest loading direction) on FeB4 show the weak Fe-B bonding is responsible for this low hardness. Moreover, these results are consistent with the value of Vickers hardness (11.7 to 32.3 GPa) by employing different empirical hardness models and below the superhardness threshold of 40 GPa. Our current results suggest FeB4 is a hard material and unlikely to become superhard.

preprint2013arXiv

Anomalous melting behavior of solid hydrogen at high pressures

Hydrogen is the most abundant element in the universe, and its properties under conditions of high temperature and pressure are crucial to understand the interior of of large gaseous planets and other astrophysical bodies. At ultra high pressures solid hydrogen has been predicted to transform into a quantum fluid, because of its high zero point motion. Here we report first principles two phase coexistence and Z method determinations of the melting line of solid hydrogen in a pressure range spanning from 30 to 600 GPa. Our results suggest that the melting line of solid hydrogen, as derived from classical molecular dynamics simulations, reaches a minimum of 367 K at about 430 GPa, at higher pressures the melting line of the atomics Cs IV phase regain a positive slope. In view of the possible importance of quantum effects in hydrogen at such low temperatures, we also determined the melting temperature of the atomic CsIV phase at pressures of 400, 500, 600 GPa, employing Feynman path integral simulations. These result in a downward shift of the classical melting line by about 100 K, and hint at a possible secondary maximum in the melting line in the region between 500 and 600 GPa, testifying to the importance of quantum effects in this system. Combined, our results imply that the stability field of the zero temperature quantum liquid phase, if it exists at all, would only occur at higher pressures than previously thought.

preprint2013arXiv

Robust Diffusive Proton Motions in Phase IV of Solid Hydrogen

Systematic first-principles molecular dynamics (MD) simulations with long simulation times (7-13 ps) for phase IV of solid hydrogen using different supercell sizes of 96, 288, 576, and 768 atoms established that the diffusive proton motions process in the graphene-like layer is an intrinsic property and independent of the simulation cell sizes. The present study highlights an often overlook issue in first-principles calculations that long time MD is essential to achieve ergodicity, which is mandatory for a proper description of dynamics of a system. The present results contradict a recent work [Phys. Rev. B 87, 174110 (2013)] in which the analysis was relied on short time slices (1-3 ps).

preprint2013arXiv

Xenon Reacts with Iron at the Conditions of the Earth's Core

Studies of the Earth's atmosphere have shown that more than 90% of xenon (Xe) is depleted compared with its abundance in chondritic meteorites. This long-standing missing Xe paradox has become the subject of considerable interest and several models for a Xe reservoir have been proposed. Whether the missing Xe is hiding in the Earth's core has remained a long unanswered question. The key to address this issue lies in the reactivity of Xe with iron (Fe, the main constituent of the Earth's core), which has been denied by earlier studies. Here we report on the first evidence of the chemical reaction of Xe and Fe at the conditions of the Earth's core, predicted through first-principles calculations and unbiased structure searching techniques. We find that Xe and Fe form a stable, inter-metallic compound of XeFe3, adopting a Cu3Au-type face-centered cubic structure above 183 GPa and at 4470 K. As the result of a Xe -> Fe charge transfer, Xe loses its chemical inertness by opening up the filled 5p electron shell and functioning as a 5p-like element, whilst Fe is unusually negatively charged, acting as an oxidant rather than a reductant as usual. Our work establishes that the Earth's core is a natural reservoir for Xe storage, and possibly provides the key to unlocking the missing Xe paradox.

preprint2012arXiv

Hybrid Functional Study Rationalizes the Simple Cubic Phase of Calcium at High Pressures

Simple cubic (SC) phase has been long experimentally determined as the high-pressure phase III of elemental calcium (Ca) since 1984. However, recent density functional calculations within semi-local approximation showed that this SC phase is structurally unstable by exhibiting severely imaginary phonons, and is energetically unstable with respect to a theoretical body-centered tetragonal I41/amd structure over the pressure range of phase III. These calculations generated extensive debates on the validity of SC phase. Here we have re-examined the SC structure by performing more precise density functional calculations within hybrid functionals of Heyd-Scuseria-Erhzerhof (HSE) and PBE0. Our calculations were able to rationalize fundamentally the phase stability of SC structure over all other known phases by evidence of its actual energetic stability above 33 GPa and its intrinsically dynamical stability without showing any imaginary phonons in the entire pressure range studied. We further established that the long-thought theoretical I41/amd structure remains stable in a narrow pressure range before entering SC phase and is actually the structure of experimental Ca-III' synthesized recently at low temperature 14 K as supported by the excellent agreement between our simulated X-ray diffraction patterns and the experimental data. Our results shed strong light on the crucial role played by the precise electron exchange energy in a proper description of the potential energy of Ca.

preprint2012arXiv

Proton Transfer in Phase IV of Solid Hydrogen and Deuterium

The recent discovery of phase IV of solid hydrogen and deuterium consisting of two alternate layers of graphenelike three-molecule rings and unbound H2 molecules have generated great interests. However, vibrational nature of phase IV remains poorly understood. Here, we report a peculiar proton transfer and a simultaneous rotation of three molecule rings in graphenelike layers predicted by ab initio variable cell molecular dynamics simulations for phase IV of solid hydrogen and deuterium at pressure ranges of from 250 to 350 GPa and temperature range of from 300 to 500 K. This proton transfer is intimately related to the particular elongation of molecules in graphenelike layers, and it becomes more pronounced with increasing pressure at the course of larger elongation of molecules. As the consequence of proton transfer, hydrogen molecules in graphenelike layers are short lived and hydrogen vibration is strongly anharmonic. Our findings provide direct explanations on the observed abrupt increase of Raman width at the formation of phase IV and its large increase with pressure.

preprint2012arXiv

Quasi-molecular and atomic phases of dense solid hydrogen

The high-pressure phases of solid hydrogen are of fundamental interest and relevant to the interior of giant planets; however, knowledge of these phases is far from complete. Particle swarm optimization (PSO) techniques were applied to a structural search, yielding hitherto unexpected high-pressure phases of solid hydrogen at pressures up to 5 TPa. An exotic quasi-molecular mC24 structure (space group C2/c, stable at 0.47-0.59 TPa) with two types of intramolecular bonds was predicted, providing a deeper understanding of molecular dissociation in solid hydrogen, which has been a mystery for decades. We further predicted the existence of two atomic phases: (i) the oC12 structure (space group Cmcm, stable at > 2.1 TPa), consisting of planar H3 clusters, and (ii) the cI16 structure, previously observed in lithium and sodium, stable above 3.5 TPa upon consideration of the zero-point energy. This work clearly revised the known zero-temperature and high-pressure (>0.47 TPa) phase diagram for solid hydrogen and has implications for the constituent structures of giant planets.

preprint2012arXiv

Room-Temperature Structures of Solid Hydrogen at High Pressures

By employing first-principles metadynamics simulations, we explore the 300 K structures of solid hydrogen over the pressure range 150-300 GPa. At 200 GPa, we find the ambient-pressure disordered hexagonal close-packed (hcp) phase transited into an insulating partially ordered hcp phase (po-hcp), a mixture of ordered graphene-like H2 layers and the other layers of weakly coupled, disordered H2 molecules. Within this phase, hydrogen remains in paired states with creation of shorter intra-molecular bonds, which are responsible for the very high experimental Raman peak above 4000 cm-1. At 275 GPa, our simulations predicted a transformation from po-hcp into the ordered molecular metallic Cmca phase (4 molecules/cell) that was previously proposed to be stable only above 400 GPa. Gibbs free energy calculations at 300 K confirmed the energetic stabilities of the po-hcp and metallic Cmca phases over all known structures at 220-242 GPa and >242 GPa, respectively. Our simulations highlighted the major role played by temperature in tuning the phase stabilities and provided theoretical support for claimed metallization of solid hydrogen below 300 GPa at 300 K.

preprint2011arXiv

Predicting Two-Dimensional Boron-Carbon Compounds by the Global Optimization Method

We adopt a global optimization method to predict two-dimensional (2D) nanostructures through the particle-swarm optimization (PSO) algorithm. By performing PSO simulations, we predict new stable structures of 2D boron-carbon (B-C) compounds for a wide range of boron concentrations. Our calculations show that: (1) All 2D B-C compounds are metallic except for BC3 which is a magic case where the isolation of carbon six-membered ring by boron atoms results in a semiconducting behavior. (2) For C-rich B-C compounds, the most stable 2D structures can be viewed as boron doped graphene structures, where boron atoms typically form 1D zigzag chains except for BC3 in which boron atoms are uniformly distributed. (3) The most stable 2D structure of BC has alternative carbon and boron ribbons with strong in-between B-C bonds, which possesses a high thermal stability above 2000K. (4) For B-rich 2D B-C compounds, there is a novel planar-tetracoordinate carbon motif with an approximate C2v symmetry.

preprint2010arXiv

Anomalous optical and electronic properties of dense sodium

Based on ab initio density-functional-theory using generalized gradient approximation, we systematically study the optical and electronic properties of the insulating dense sodium phase (Na-hp4) reported recently [Ma \textit{et al.}, Nature \textbf{458}, 182 (2009)]. The structure is found optically anisotropic and transparent to visible light, which can be well interpreted using its electronic band structure and angular moment decomposed density of states. Through the bader analysis of Na-hp4 at different pressures, we conclude that ionicity exists in the structure and becomes stronger with increasing pressure. In addition, the absorption spectra in the energy range from 1.4 to 2.4 eV are compared with recent experimental results and found good agreement. It is found that the deep-lying valence electrons participate in the interband transition.