BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Yanming Ma

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

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.supr-conphysics.comp-phphysics.chem-phcond-mat.mes-hallcond-mat.str-elphysics.geo-phquant-ph

Catalog footprint

What is connected

40works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

40 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlocal Pseudopotential Energy Density Functional for Orbital-Free Density Functional Theory

Orbital-free density functional theory (OF-DFT) runs at low computational cost that scales linearly with the number of simulated atoms, making it suitable for large-scale material simulations. It is generally considered that OF-DFT strictly requires the use of local pseudopotentials, rather than orbital-dependent nonlocal pseudopotentials, for the calculation of electron-ion interaction energies, as no orbitals are available. This is unfortunate situation since the nonlocal pseudopotentials are known to give much better transferability and calculation accuracy than local ones. We report here the development of a theoretical scheme that allows the direct use of nonlocal pseudopotentials in OF-DFT. In this scheme, a nonlocal pseudopotential energy density functional is derived by the projection of nonlocal pseudopotential onto the non-interacting density matrix (instead of 'orbitals') that can be approximated explicitly as a functional of electron density. Our development defies the belief that nonlocal pseudopotentials are not applicable to OF-DFT, leading to the creation of an alternate theoretical framework of OF-DFT that works superior to the traditional one.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nonlocal Kinetic Energy Density Functionals for Isolated Systems via Local Density Approximation Kernels

Despite a large number of nonlocal kinetic energy density functionals (KEDFs) available for large-scale calculations, most of those nonlocal KEDFs designed for the extended systems cannot be directly applied to isolated systems. In this manuscript, we proposed a generalized scheme to construct nonlocal KEDFs via the local density approximation kernels and construct a family of KEDFs for simulations of isolated systems within orbital-free density functional theory. The performance of KEDFs has been demonstrated by several clusters encompassing Mg, Si and GaAs. The results show that our constructed KEDFs can achieve high numerical accuracy and stability for random clusters, therefore, making orbital-free density functional theory accessible for practical simulations of isolated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A linear scaling method to evaluate the ion-electron potential of crystalline solids

We propose a simple linear scaling expression in reciprocal space for evaluating the ion--electron potential of crystalline solids. The expression replaces the long-range ion--electron potential with an equivalent localized charge distribution and corresponding boundary conditions on the unit cell. Given that no quadratic scaling structure factor is required---as used in traditional methods---the expression shows inherent linear behavior, and is well suited to simulating large-scale systems within orbital-free density functional theory. The scheme is implemented in the ATLAS software package and benchmarked by using a solid Mg bcc lattice containing tens of thousands of atoms in the unit cell. The test results show that the method can efficiently model large crystals with high computational accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Design of ternary alkaline-earth metal Sn(II) oxides with potential good p-type conductivity

Oxides with good p-type conductivity have been long sought after to achieve high performance all-oxide optoelectronic devices. Divalent Sn(II) based oxides are promising candidates because of their rather dispersive upper valence bands caused by the Sn-5s/O-2p anti-bonding hybridization. There are so far few known Sn(II) oxides being p-type conductive suitable for device applications. Here, we present via first-principles global optimization structure searches a material design study for a hitherto unexplored Sn(II)-based system, ternary alkaline-earth metal Sn(II) oxides in the stoichiometry of MSn2O3 (M = Mg, Ca, Sr, Ba). We identify two stable compounds of SrSn2O3 and BaSn2O3, which can be stabilized by Sn-rich conditions in phase stability diagrams. Their structures follow the Zintl behaviour and consist of basic structural motifs of SnO3 tetrahedra. Unexpectedly they show distinct electronic properties with band gaps ranging from 1.90 (BaSn2O3) to 3.15 (SrSn2O3) eV, and hole effective masses ranging from 0.87 (BaSn2O3) to above 6.0 (SrSn2O3) m0. Further exploration of metastable phases indicates a wide tunability of electronic properties controlled by the details of the bonding between the basic structural motifs. This suggests further exploration of alkaline-earth metal Sn(II) oxides for potential applications requiring good p-type conductivity such as transparent conductors and photovoltaic absorbers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High-pressure Phase Stability and Superconductivity of Pnictogen Hydrides and Chemical Trends for Compressed Hydrides

Binary hydrides formed by the pnictogens of phosphorus, arsenic and antimony are studied at high pressures using first principles methods. Stable structures are predicted and their electronic, vibrational and superconducting properties are investigated. We predict that SbH$_{4}$ and AsH$_{8}$ will be high-temperature superconductors at megabar pressures, with critical temperatures in excess of 100 K. The highly symmetric hexagonal SbH$_{4}$ phase is predicted to be stabilized above about 150 GPa, which is readily achievable in diamond anvil cell experiments. We find that all phosphorus hydrides are metastable with respect to decomposition into the elements within the pressure range studied. Trends based on our results and literature data reveal a connection between the high-pressure behaviors and ambient-pressure chemical quantities which provides insight into understanding which elements may form hydrogen-rich high-temperature superconducting phases at high pressures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inverse Design of Inorganic Electrides

Electrides are ionic solids that consist of cationic frameworks and anionic electrons trapped in the voids of lattices. Organic electrides exist in a large abundance, but the thermal instability at room temperature and sensitivity to moisture are bottlenecks that limit their practical uses. Known inorganic electrides are rare but appear to have high thermal and chemical stability and exhibit promising applications as electron-emitting materials, superior catalysts and strong reducing agents. Here, we report a developed inverse-design method that can be used to search for a large variety of inorganic electrides. Our method utilizes the intrinsic property of interstitial electron localization of electrides as the global variable function being incorporated into the swarm-intelligence based structure searches. Through screening 99 binary ionic compounds, we have designed 89 new inorganic electrides that are classified into three-, two-, and zero-dimensional species according to the way that the interstitial electrons are localized and the conductive properties of the systems. Our work reveals the rich abundance of inorganic electrides by extending them into more general forms and provides new structure types for electrides that are not thought of as before.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sn(II)-containing phosphates as optoelectronic materials

We theoretically investigate Sn(II) phosphates as optoelectronic materials using first principles calculations. We focus on known prototype materials Sn$_n$P$_2$O$_{5+n}$ (n=2, 3, 4, 5) and a previously unreported compound, SnP$_2$O$_6$ (n=1), which we find using global optimization structure prediction. The electronic structure calculations indicate that these compounds all have large band gaps above 3.2 eV, meaning their transparency to visible light. Several of these compounds show relatively low hole effective masses ($\sim$2-3 m$_0$), comparable the electron masses. This suggests potential bipolar conductivity depending on doping. The dispersive valence band-edges underlying the low hole masses, originate from the anti-bonding hybridization between the Sn 5s orbitals and the phosphate groups. Analysis of structure-property relationships for the metastable structures generated during structure search shows considerable variation in combinations of band gap and carrier effective masses, implying chemical tunability of these properties. The unusual combinations of relatively high band gap, low carrier masses and high chemical stability suggests possible optoelectronic applications of these Sn(II) phosphates, including p-type transparent conductors. Related to this, calculations for doped material indicate low visible light absorption, combined with high plasma frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

ATLAS: A Real-Space Finite-Difference Implementation of Orbital-Free Density Functional Theory

Orbital-free density functional theory (OF-DFT) is a promising method for large-scale quantum mechanics simulation as it provides a good balance of accuracy and computational cost. Its applicability to large-scale simulations has been aided by progress in constructing kinetic energy functionals and local pseudopotentials. However, the widespread adoption of OF-DFT requires further improvement in its efficiency and robustly implemented software. Here we develop a real-space finite-difference method for the numerical solution of OF-DFT in periodic systems. Instead of the traditional self-consistent method, a powerful scheme for energy minimization is introduced to solve the Euler--Lagrange equation. Our approach engages both the real-space finite-difference method and a direct energy-minimization scheme for the OF-DFT calculations. The method is coded into the ATLAS software package and benchmarked using periodic systems of solid Mg, Al, and Al$_{3}$Mg. The test results show that our implementation can achieve high accuracy, efficiency, and numerical stability for large-scale simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constructing optimal local pseudopotentials from first principles

Local pseudopotential (LPP) is an important component of the orbital free density functional theory (OF-DFT), which is a promising large scale simulation method that can still maintain information of electron state in materials. Up to date, LPP is usually extracted from the solid state DFT calculations. It is unclear how to assess its transferability while applying to a much different chemical environment. Here we reveal a fundamental relation between the first principles norm-conserving PP (NCPP) and the LPP. Using the optimized effective potential method developed for exchange functional, we demonstrate that the LPP can be constructed optimally from the NCPP for a large number of elements. Our theory also reveals that the existence of an LPP is intrinsic to the elements, irrespective to the parameters used for the construction. Our method provides a unified method in constructing and assessing LPP in the framework of first principles pseudopotentials.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetism in Na-filled Fe-based skutterudites

The interplay of superconductivity and magnetism is a subject of ongoing interest, stimulated most recently by the discovery of Fe-based superconductivity and the recognition that spin-fluctuations near a magnetic quantum critical point may provide an explanation for the superconductivity and the order parameter. Here we investigate magnetism in the Na filled Fe-based skutterudites using first principles calculations. NaFe4Sb12 is a known ferromagnet near a quantum critical point. We find a ferromagnetic metallic state for this compound driven by a Stoner type instability, consistent with prior work. In accord with prior work, the magnetization is overestimated, as expected for a material near an itinerant ferromagnetic quantum critical point. NaFe4P12 also shows a ferromagnetic instability at the density functional level, but this instability is much weaker than that of NaFe4Sb12, possibly placing it on the paramagnetic side of the quantum critical point. NaFe4As12 shows intermediate behavior. We also present results for skutterudite FeSb3, which is a metastable phase that has been reported in thin film form.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Metallic Icosahedron Phase of Sodium at Terapascal Pressures

Alkali metals exhibit unexpected structures and electronic behavior at high pressures. Compression of metallic sodium (Na) to 200 GPa leads to the stability of a wide-band-gap insulator with the double hexagonal hP4 structure. Post-hP4 structures remain unexplored, but they are important for addressing the question of the pressure at which Na reverts to a metal. Here we report the reentrant metallicity of Na at the very high pressure of 15.5 terapascal (TPa), predicted using first-principles structure searching simulations. Na is therefore insulating over the large pressure range of 0.2-15.5 TPa. Unusually, Na adopts an oP8 structure at pressures of 117-125 GPa, and the same oP8 structure at 1.75-15.5 TPa. Metallization of Na occurs on formation of a stable and striking body-centered cubic cI24 electride structure consisting of Na12 icosahedra, each housing at its center about one electron which is not associated with any Na ions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Prediction of Silicon-Based Layered Structures for Optoelectronic Applications

A method based on the particle swarm optimization (PSO) algorithm is presented to design quasi-two-dimensional (Q2D) materials. With this development, various single-layer and bi-layer materials in C, Si, Ge, Sn, and Pb were predicted. A new Si bi-layer structure is found to have a much-favored energy than the previously widely accepted configuration. Both single-layer and bi-layer Si materials have small band gaps, limiting their usages in optoelectronic applications. Hydrogenation has therefore been used to tune the electronic and optical properties of Si layers. We discover two hydrogenated materials of layered Si8H2 and Si6H2 possessing quasi-direct band gaps of 0.75 eV and 1.59 eV, respectively. Their potential applications for light emitting diode and photovoltaics are proposed and discussed. Our study opened up the possibility of hydrogenated Si layered materials as next-generation optoelectronic devices.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stable Xenon Nitride at High Pressures

Nitrides in many ways are fascinating since they often appear as superconductors, high energy density and hard materials. Though there exist a large variety of nitrides, noble gas nitrides are long missing in nature. Pursuit of noble gas nitrides has therefore become the subject of topical interests, but remains as a great challenge since molecular nitrogen (N2, a major form of nitrogen) and noble gases are both inert systems and do not interact at normal conditions. We show through a swarm structure searching simulation that high pressure can lift the reactivity of both N2 and xenon (Xe), making chemical reaction of them possible. The resultant nitride has a peculiar stoichiometry of XeN6, possessing a high-energy-density of approximately 2.4 kJg-1, rivaling that of the modern explosives. Chemically, XeN6 is more intriguing with the appearance of chaired N6 hexagons and an emergent 12-fold Xe by acceptance of unprecedentedly 12 Xe-N weak covalent bonds. Our work opens up the possibility of achieving Xe nitrides whose formation is long sought as impossible.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tuning Optical Properties of Transparent Conducting Barium Stannate by Dimensional Reduction

We report calculations of the electronic structure and optical properties of doped $n$-type perovskite BaSnO3 and layered perovskites. While doped BaSnO$_3$ retains its transparency for energies below the valence to conduction band onset, the doped layered compounds exhibit below band edge optical conductivity due to transitions from the lowest conduction band. This gives absorption in the visible for Ba2SnO4. Thus it is important to minimize this phase in transparent conducting oxide (TCO) films. Ba3Sn2O7 and Ba4Sn3O10 have strong transitions only in the red and infrared, respectively. Thus there may be opportunities for using these as wavelength filtering TCO.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ultrahigh Pressure Superconductivity in Molybdenum Disulfide

Superconductivity commonly appears under pressure in charge density wave (CDW)-bearing transition metal dichalcogenides (TMDs), but has emerged so far only via either intercalation with electron donors or electrostatic doping in CDW-free TMDs. Theoretical calculations have predicted that the latter should be metallized through bandgap closure under pressure, but superconductivity remained elusive in pristine 2H-MoS2 upon substantial compression, where a pressure of up to 60 GPa only evidenced the metallic state. Here we report the emergence of superconductivity in pristine 2H-MoS2 at 90 GPa. The maximum onset transition temperature Tc(onset) of 11.5 K, the highest value among TMDs and nearly constant from 120 up to 200 GPa, is well above that obtained by chemical doping but comparable to that obtained by electrostatic doping. Tc(onset) is more than an order of magnitude larger than present theoretical expectations, raising questions on either the current calculation methodologies or the mechanism of the pressure-induced pairing state. Our findings strongly suggest further experimental and theoretical efforts directed toward the study of the pressure-induced superconductivity in all CDW-free TMDs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The metallization and superconductivity of dense hydrogen sulfide

Hydrogen sulfide (H2S) is a prototype molecular system and a sister molecule of water. The phase diagram of solid H2S at high pressures remains largely unexplored arising from the challenges in dealing with the looser S-H bond and larger atomic core difference between H and S. Metallization is yet achieved for water ice, but it was established for H2S above 96 GPa. However, the metallic structure of H2S remains elusive, greatly impeding the understanding of its metallicity and the potential superconductivity. We have performed an extensive structural study on solid H2S under high pressures through unbiased first-principles structure predictions based on swarm intelligence. Besides the findings of best-known candidate structures for nonmetallic phases IV and V, we are able to establish stable metallic structures violating an earlier proposal of elemental decomposition into sulfur and hydrogen [PRL 85, 1254 (2000)]. Our study unraveled a superconductive potential of metallic H2S with an estimated maximal transition temperature of ~ 80 K at 160 GPa, higher than those predicted for most archetypal hydrogen-containing compounds (e.g., SiH4 and GeH4, etc).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CALYPSO: a method for crystal structure prediction

We have developed a software package CALYPSO (Crystal structure AnaLYsis by Particle Swarm Optimization) to predict the energetically stable/metastable crystal structures of materials at given chemical compositions and external conditions (e.g., pressure). The CALYPSO method is based on several major techniques (e.g. particle-swarm optimization algorithm, symmetry constraints on structural generation, bond characterization matrix on elimination of similar structures, partial random structures per generation on enhancing structural diversity, and penalty function, etc) for global structural minimization from scratch. All of these techniques have been demonstrated to be critical to the prediction of global stable structure. We have implemented these techniques into the CALYPSO code. Testing of the code on many known and unknown systems shows high efficiency and high successful rate of this CALYPSO method [Wang et al., Phys. Rev. B 82 (2010) 094116][1]. In this paper, we focus on descriptions of the implementation of CALYPSO code and why it works.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diamondoid Structure of Polymeric Nitrogen at High Pressures

High-pressure polymeric structures of nitrogen have attracted great attention owing to their potential application as high-energy-density materials. We report the density functional structural prediction of the unexpected stabilization of a diamondoid (or N10-cage) structure of polymeric nitrogen at high pressures. The structure adopts a highly symmetric body-centered cubic form with lattice sites occupied by N10 tetracyclic cages, each of which consists of 10 atoms and is covalently bonded with its six next-nearest N10 cages. The prediction of this diamondoid structure rules out the earlier proposed helical tunnel phase and demonstrates the high-order nature of polymeric nitrogen at extreme high pressures. Diamondoid nitrogen is a wide-gap insulator and energetically more favorable than the experimental cubic gauche and previously predicted layered Pba2 phases above 263 GPa, a pressure which is accessible to high pressure experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High pressures make Hg a transition metal in a thermodynamically stable solid

The appropriateness of including Hg among the transition metals has been debated for a long time. Although the synthesis of HgF$_{4}$ molecules in gas phase was reported before, the molecules show strong instabilities and dissociate. Therefore, the transition metal propensity of Hg remains an open question. Here, we propose that high pressure provides a controllable method for preparing unusual oxidation states of matter. Using an advanced structure search method based on first-principles electronic structure calculations, we predict that under high pressures, Hg can transfer the electrons in its outmost $d$ shell to F atoms, thereby acting as a transition metal. Oxidation of Hg to the +4 state yielded thermodynamically stable molecular crystals consisting of HgF$_{4}$ planar molecules, a typical geometry for $d^{8}$ metal centers.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spiral Chain O4 Form of Dense Oxygen

Oxygen is in many ways a unique element: the only known diatomic molecular magnet and the capability of stabilization of the hitherto unexpected O8 cluster structure in its solid form at high pressure. Molecular dissociations upon compression as one of the fundamental problems were reported for other diatomic solids (e.g., H2, I2, Br2, and N2), but it remains elusive for solid oxygen, making oxygen an intractable system. We here report the theoretical prediction on the dissociation of molecular oxygen into a polymeric spiral chain O4 structure (θ-O4) by using first-principles calypso method on crystal structure prediction. The θ-O4 stabilizes above 2 TPa and has been observed as the third high pressure phase of sulfur (S-III). We find that the molecular O8 phase remains extremely stable in a large pressure range of 0.008 - 2 TPa, whose breakdown is driven by the pressure-induced instability of a transverse acoustic phonon mode at zone boundary, leading to the ultimate formation of θ-O4. Remarkably, stabilization of θ-O4 turns oxygen from a superconductor into an insulator with a wide band gap (approximately 5.9 eV) originating from the sp3-like hybridized orbitals of oxygen and the localization of valence electrons. (This is a pre-print version of the following article: Li Zhu et al, Spiral chain O4 form of dense oxygen, Proc. Natl. Acad. Sci. U.S.A. (2011), doi: 10.1073/pnas.1119375109, which has been published online at http://www.pnas.org/content/early/2011/12/27/1119375109 .)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Superconductive "sodalite"-like clathrate calcium hydride at high pressures

Hydrogen-rich compounds hold promise as high-temperature superconductors under high pressures. Recent theoretical hydride structures on achieving high-pressure superconductivity are composed mainly of H2 fragments. Through a systematic investigation of Ca hydrides with different hydrogen contents using particle-swam optimization structural search, we show that in the stoichiometry CaH6 a body-centred cubic structure with hydrogen that forms unusual "sodalite" cages containing enclathrated Ca stabilizes above pressure 150 GPa. The stability of this structure is derived from the acceptance by two H2 of electrons donated by Ca forming a "H4" unit as the building block in the construction of the 3-dimensional sodalite cage. This unique structure has a partial occupation of the degenerated orbitals at the zone centre. The resultant dynamic Jahn-Teller effect helps to enhance electron-phonon coupling and leads to superconductivity of CaH6. A superconducting critical temperature (Tc) of 220-235 K at 150 GPa obtained from the solution of the Eliashberg equations is the highest among all hydrides studied thus far.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crystal Structure Prediction via Particle Swarm Optimization

We have developed a powerful method for crystal structure prediction from "scratch" through particle swarm optimization (PSO) algorithm within the evolutionary scheme. PSO technique is dramatically different with the genetic algorithm and has apparently avoided the use of evolution operators (e.g., crossover and mutation). The approach is based on a highly efficient global minimization of free energy surfaces merging total-energy calculations via PSO technique and requires only chemical compositions for a given compound to predict stable or metastable structures at given external conditions (e.g., pressure). A particularly devised geometrical structure factor method which allows the elimination of similar structures during structure evolution was implemented to enhance the structure search efficiency. The application of designed variable unit cell size technique has greatly reduced the computational cost. Moreover, the symmetry constraint imposed in the structure generation enables the realization of diverse structures, leads to significantly reduced search space and optimization variables, and thus fastens the global structural convergence. The PSO algorithm has been successfully applied to the prediction of many known systems (e.g., elemental, binary and ternary compounds) with various chemical bonding environments (e.g., metallic, ionic, and covalent bonding). The remarkable success rate demonstrates the reliability of this methodology and illustrates the great promise of PSO as a major technique on crystal structure determination.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Density functional study of elastic and vibrational properties of the Heusler-type alloys Fe$_2$VAl and Fe$_2$VGa

The structural and elastic properties as well as phonon-dispersion relations of the Heusler-type alloys Fe$_2$VAl and Fe$_2$VGa are computed using density-functional and density-functional perturbation theory within the generalized-gradient approximation. The calculated equilibrium lattice constants agree well with the experimental values. The elastic constants of Fe$_2$VAl and Fe$_2$VGa are predicted for the first time. From the elastic constants the shear modulus, Young's modulus, Poisson's ratio, sound velocities and Debye temperatures are obtained. By analyzing the ratio between the bulk and shear modulii, we conclude that both Fe$_2$VAl and Fe$_2$VGa are brittle in nature. The computed phonon-dispersion relation shows that both compounds are dynamically stable in the L1$_2$ structure without any imaginary phonon frequencies. The isomer shifts of Fe in the two compounds are discussed in terms of the Fe s partial density of states, which reveal larger ionicity/less hybridization in Fe$_2$VGa than in Fe$_2$VAl. For the same reason the Cauchy pressure is negative in Fe$_2$VAl but positive in Fe$_2$VGa

0 citations0 reviewsNo code linkedOpen access