Source author record

Kai-Ming Ho

Kai-Ming Ho 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

25works
9topics
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

25 published item(s)

preprint2025arXiv

Ab initio superionic-liquid phase diagram of Fe1-xOx under Earth's inner core conditions

The superionic state is a phase of matter in which liquid-like ionic mobility coexists with a solid crystalline lattice. Recently identified in Earth's inner core (IC), this state has attracted considerable attention for its unique kinetic behavior and geophysical implications. However, the ab initio phase diagram describing the equilibrium between the superionic phase and the liquid solution under core conditions remains largely unexplored. Here, we present a thermodynamic approach to compute the Gibbs free energy and construct the ab initio superionic-liquid phase diagram for the Fe1-xOx system under IC conditions. We find that oxygen forms superionic states in both hcp and bcc Fe phases, with a pronounced influence on cooperative diffusion of iron in the bcc lattice. The stability fields of these superionic phases are sensitive to oxygen stoichiometry. The presence of superionic states leads to a higher oxygen concentration in the IC than previously estimated. Our work establishes a framework for investigating superionic-liquid equilibria under extreme conditions.

preprint2023arXiv

Benchmarks and results of the two-band Hubbard model from the Gutzwiller conjugate gradient minimization theory

Ground-state properties, such as energies and double occupancies, of a one-dimensional two-band Hubbard model are calculated using a first principles Gutzwiller conjugate gradient minimization theory. The favorable agreement with the results from the density matrix renormalization group theory demonstrates the accuracy of our method. A rotationally invariant approach is further incorporated into the method to greatly reduce the computational complexity with a speedup of 300 times. Moreover, we investigate the Mott transition between a metal and a Mott insulator by evaluating the charge gap. With greatly reduced computational effort, our method reproduces the phase diagram in reasonable agreement with the density matrix renormalization group theory.

preprint2022arXiv

A deep machine learning potential for atomistic simulation of Fe-Si-O systems under Earth's outer core conditions

Using artificial neural-network machine learning (ANN-ML) to generate interatomic potentials has been demonstrated to be a promising approach to address the long-standing challenge of accuracy versus efficiency in molecular dynamics (MD) simulations. Here, taking the Fe-Si-O system as a prototype, we show that accurate and transferable ANN-ML potentials can be developed for reliable MD simulations of materials at high-pressure and high-temperature conditions of the Earth's outer core. The ANN-ML potential for Fe-Si-O system is trained by fitting to the energies and forces of related binaries and ternary liquid structures at high pressures and temperatures obtained by first-principles calculations based on density functional theory (DFT). We show that the generated ANN-ML potential describes well the structure and dynamics of liquid phases of this complex system. The efficient ANN-ML potential with DFT accuracy provides a promising scheme for accurate atomistic simulations of structures and dynamics of complex Fe-Si-O system in the Earth's outer core.

preprint2022arXiv

Electron-phonon coupling strength from ab initio frozen-phonon approach

We propose a fast method for high-throughput screening of potential superconducting materials. The method is based on calculating metallic screening of zone-center phonon modes, which provides an accurate estimate for the electron-phonon coupling strength. This method is complementary to the recently proposed Rigid Muffin Tin (RMT) method, which amounts to integrating the electron-phonon coupling over the entire Brillouin zone (as opposed to the zone center), but in a relatively inferior approximation. We illustrate the use of this method by applying it to MgB$_\text{2}$, where the high-temperature superconductivity is known to be driven largely by the zone-center modes, and compare it to a sister compound AlB$_\text{2}$. We further illustrate the usage of this descriptor by screening a large number of binary hydrides, for which accurate first-principle calculations of electron-phonon coupling have been recently published. Together with the RMT descriptor, this method opens a way to perform initial high-throughput screening in search of conventional superconductors via machine learning or data mining.

preprint2022arXiv

High-throughput screening of strong electron-phonon couplings in ternary metal diborides

We perform a high-throughput screening on phonon-mediated superconductivity in ternary metal diboride structure with alkali, alkaline earth, and transition metals. We find 17 ground states and 78 low-energy metastable phases. From fast calculations of zone-center electron-phonon coupling, 43 compounds are revealed to show electron-phonon coupling strength higher than that of MgB2. An anti-correlation between energetic stability and electron-phonon coupling strength is identified. We suggest two phases, i.e., Li3ZrB8 and Ca3YB8, to be synthesized, which show reasonable energetic stability and superconducting critical temperature.

preprint2022arXiv

Structure and motifs of iron oxides from 1 to 3 TPa

Iron oxides are fundamental components of planet-forming materials. Understanding the Fe-O system's behavior and properties under high pressure can help us identify many new phases and states possible in exoplanetary interiors, especially terrestrial ones. Using the adaptive genetic algorithm (AGA), we investigate the structure of iron oxides for a wide range of stoichiometries ($0.25\leq x_O \leq 0.8$) at 1, 2, and 3 TPa. Five unreported ground-state structures with Fe$_2$O, FeO, Fe$_3$O$_5$, FeO$_2$, and FeO$_4$ compositions are identified. The calculated density of states (DOS) suggests that, except for FeO$_4$, all phases are metallic, but their carrier densities decrease with increasing pressure and oxygen content. The cluster alignment analysis of stable and metastable phases shows that several motifs may co-exist in a structure of iron oxides with low O content. In contrast, most iron oxides with high O content adopt a simple BCC motif at TPa pressures. Our results provide a crystal structure database of iron oxides for modeling and understanding the interiors of exoplanets.

preprint2022arXiv

Two-step nucleation of the Earth's inner core

It has long been assumed the Earth's solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earth's core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.

preprint2020arXiv

Efficient step-merged quantum imaginary time evolution algorithm for quantum chemistry

We develop a resource efficient step-merged quantum imaginary time evolution approach (smQITE) to solve for the ground state of a Hamiltonian on quantum computers. This heuristic method features a fixed shallow quantum circuit depth along the state evolution path. We use this algorithm to determine binding energy curves of a set of molecules, including H$_2$, H$_4$, H$_6$, LiH, HF, H$_2$O and BeH$_2$, and find highly accurate results. The required quantum resources of smQITE calculations can be further reduced by adopting the circuit form of the variational quantum eigensolver (VQE) technique, such as the unitary coupled cluster ansatz. We demonstrate that smQITE achieves a similar computational accuracy as VQE at the same fixed-circuit ansatz, without requiring a generally complicated high-dimensional non-convex optimization. Finally, smQITE calculations are carried out on Rigetti quantum processing units (QPUs), demonstrating that the approach is readily applicable on current noisy intermediate-scale quantum (NISQ) devices.

preprint2020arXiv

Evidence for large Rashba splitting in PtPb4 from Angle-Resolved Photoemission Spectroscopy

We studied the electronic structure of PtPb$_{4}$ using laser angle-resolved photoemission spectroscopy(ARPES) and density functional theory(DFT) calculations. This material is closely related to PtSn$_{4}$, which exhibits exotic topological properties such as Dirac node arcs. Fermi surface(FS) of PtPb$_{4}$ consists of two electron pockets at the center of the Brillouin zone(BZ) and several hole pockets around the zone boundaries. Our ARPES data reveals significant Rashba splitting at the $Γ$ point in agreement with DFT calculations. The presence of Rashba splitting may render this material of potential interest for spintronic applications.

preprint2020arXiv

Ground state properties of one-dimensional and two-dimensional Hubbard model from Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

preprint2020arXiv

Gutzwiller Hybrid Quantum-Classical Computing Approach for Correlated Materials

Rapid progress in noisy intermediate-scale quantum (NISQ) computing technology has led to the development of novel resource-efficient hybrid quantum-classical algorithms, such as the variational quantum eigensolver (VQE), that can address open challenges in quantum chemistry, physics and material science. Proof-of-principle quantum chemistry simulations for small molecules have been demonstrated on NISQ devices. While several approaches have been theoretically proposed for correlated materials, NISQ simulations of interacting periodic models on current quantum devices have not yet been demonstrated. Here, we develop a hybrid quantum-classical simulation framework for correlated electron systems based on the Gutzwiller variational embedding approach. We implement this framework on Rigetti quantum processing units (QPUs) and apply it to the periodic Anderson model, which describes a correlated heavy electron band hybridizing with non-interacting conduction electrons. Our simulation results quantitatively reproduce the known ground state quantum phase diagram including metallic, Kondo and Mott insulating phases. This is the first fully self-consistent hybrid quantum-classical simulation of an infinite correlated lattice model executed on QPUs, demonstrating that the Gutzwiller hybrid quantum-classical embedding framework is a powerful approach to simulate correlated materials on NISQ hardware. This benchmark study also puts forth a concrete pathway towards practical quantum advantage on NISQ devices.

preprint2020arXiv

Shallow-circuit variational quantum eigensolver based on symmetry-inspired Hilbert space partitioning for quantum chemical calculations

Development of resource-friendly quantum algorithms remains highly desirable for noisy intermediate-scale quantum computing. Based on the variational quantum eigensolver (VQE) with unitary coupled cluster ansatz, we demonstrate that partitioning of the Hilbert space made possible by the point group symmetry of the molecular systems greatly reduces the number of variational operators by confining the variational search within a subspace. In addition, we found that instead of including all subterms for each excitation operator, a single-term representation suffices to reach required accuracy for various molecules tested, resulting in an additional shortening of the quantum circuit. With these strategies, VQE calculations on a noiseless quantum simulator achieve energies within a few meVs of those obtained with the full UCCSD ansatz for $\mathrm{H}_4$ square, $\mathrm{H}_4$ chain and $\mathrm{H}_6$ hexagon molecules; while the number of controlled-NOT (CNOT) gates, a measure of the quantum-circuit depth, is reduced by a factor of as large as 35. Furthermore, we introduced an efficient "score" parameter to rank the excitation operators, so that the operators causing larger energy reduction can be applied first. Using $\mathrm{H}_4$ square and $\mathrm{H}_4$ chain as examples, We demonstrated on noisy quantum simulators that the first few variational operators can bring the energy within the chemical accuracy, while additional operators do not improve the energy since the accumulative noise outweighs the gain from the expansion of the variational ansatz.

preprint2020arXiv

Theoretical prediction of a highly responsive material: Spin fluctuations and superconductivity in FeNiB2 system

By analyzing Fe-Ni-B compositional diagram we predict an energetically and dynamically stable FeNiB2 compound. This system belongs to the class of highly responsive state of material, as it is very sensitive to the external perturbations. This state is also characterized by a high level of spin fluctuations which strongly influence possible magnetic long- and short-range orders. Furthermore, we demonstrate that these antiferromagnetically dominating fluctuations could lead to the appearance of spin mediated superconductivity. The obtained results suggest a promising avenue for the search of strong spin fluctuation systems and related superconductors.

preprint2018arXiv

Spatially-correlated Site Occupancy in the Nonstoichiometric Meta-stable ε-Al60Sm11 Phase during Devitrification of Al-10.2 at.% Sm Glasses

A metastable ε-Al60Sm11 phase appears during the initial devitrification of as-quenched Al-10.2 at.% Sm glasses. The ε phase is nonstoichiometric in nature since Al occupation is observed on the 16f Sm lattice sites. Scanning transmission electron microscopic images reveal profound spatial correlation of Sm content on these sites, which cannot be explained by the "average crystal" description from Rietveld analysis of diffraction data. Thermodynamically favorable configurations, established by Monte Carlo (MC) simulations based on a cluster-expansion model, also give qualitatively different correlation functions from experimental observations. On the other hand, molecular dynamics simulations of the growth of ε-Al60Sm11 in undercooled liquid show that when the diffusion range of Sm is limited to ~ 4 Å, the correlation function of the as-grown crystal structure agrees well with that of the STEM images. Our results show that kinetic effects, especially the limited diffusivity of Sm atoms plays the fundamental role in determining the nonstoichiometric site occupancies of the ε-Al60Sm11 phase during the crystallization process.

preprint2017arXiv

A self-contained algorithm for determination of solid-liquid equilibria in an alloy system

We describe a self-contained procedure to evaluate the free energy of liquid and solid phases of an alloy system. The free energy of a single-element solid phase is calculated with thermodynamic integration using the Einstein crystal as the reference system. Then, free energy difference between the solid and liquid phases is calculated by Gibbs-Duhem integration. The central part of our method is the construction of a reversible alchemical path connecting a pure liquid and a liquid alloy to calculate the mixing enthalpy and entropy. We have applied the method to calculate the free energy of solid and liquid phases in the Al-Sm system. The driving force for fcc-Al nucleation in Al-Sm liquid and the melting curve for fcc-Al and Al3Sm are also calculated.

preprint2016arXiv

Ferromagnetic Quantum Critical Point Avoided by the Appearance of Another Magnetic Phase in LaCrGe$_3$ under Pressure

The temperature-pressure phase diagram of the ferromagnet LaCrGe$_3$ is determined for the first time from a combination of magnetization, muon-spin-rotation and electrical resistivity measurements. The ferromagnetic phase is suppressed near $2.1$~GPa, but quantum criticality is avoided by the appearance of a magnetic phase, likely modulated, AFM$_Q$. Our density functional theory total energy calculations suggest a near degeneracy of antiferromagnetic states with small magnetic wave vectors $Q$ allowing for the potential of an ordering wave vector evolving from $Q=0$ to finite $Q$, as expected from the most recent theories on ferromagnetic quantum criticality. Our findings show that LaCrGe$_3$ is a very simple example to study this scenario of avoided ferromagnetic quantum criticality and will inspire further study on this material and other itinerant ferromagnets.

preprint2016arXiv

Large magnetic anisotropy predicted for rare-earth free Fe16-xCoxN2 alloys

Structures and magnetic properties of Fe16-xCoxN2 are studied using adaptive genetic algorithm and first-principles calculations. We show that substituting Fe by Co in Fe16N2 with Co/Fe ratio smaller than 1 can greatly improve the magnetic anisotropy of the material. The magnetocrystalline anisotropy energy from first-principles calculations reaches 3.18 MJ/m3 (245.6 μeV per metal atom) for Fe12Co4N2, much larger than that of Fe16N2 and is one of the largest among the reported rare-earth free magnets. From our systematic crystal structure searches, we show that there is a structure transition from tetragonal Fe16N2 to cubic Co16N2 in Fe16-xCoxN2 as the Co concentration increases, which can be well explained by electron counting analysis. Different magnetic properties between the Fe-rich (x < 8) and Co-rich (x > 8) Fe16-xCoxN2 is closely related to the structural transition.

preprint2015arXiv

A coefficient average approximation towards Gutzwiller wavefunction formalism

Gutzwiller wavefunction is a physically well motivated trial wavefunction for describing correlated electron systems. In this work, a new approximation is introduced to facilitate evaluation of the expectation value of any operator within the Gutzwiller wavefunction formalism. The basic idea is to make use of a specially designed average over Gutzwiller wavefunction coefficients expanded in the many-body Fock space to approximate the ratio of expectation values between a Gutzwiller wavefunction and its underlying noninteracting wavefunction. To check with the standard Gutzwiller approximation (GA), we test its performance on single band systems and find quite interesting properties. On finite systems, we noticed that it gives superior performance than GA, while on infinite systems it asymptotically approaches GA. Analytic analysis together with numerical tests are provided to support this claimed asymptotic behavior. At the end, possible improvements on the approximation and its generalization towards multiband systems are illustrated and discussed.

preprint2015arXiv

Electronic structure and energetics of praseodymium and $α$-plutonium under pressure

We develop a new implementation of the Gutzwiller approximation (GA) and interface it with the local density approximation (LDA). This formulation enables us to study complex $4f$ and $5f$ systems. We perform calculations of praseodymium and $α$-plutonium under pressure, which compare very well with the experiments. Our study of praseodymium indicates that both structure change and $f$-delocalization are important to obtain the correct phase diagram and, in particular, the pressure-induced volume-collapse transition. Our calculations of $α$-plutonium indicate that, even though the $f$ electrons are delocalized in this phase, the electron-correlations affect substantially its electronic structure and thermodynamical properties.

preprint2015arXiv

Exploration of tetrahedral structures in silicate cathodes using a motif-network scheme

Using a motif-network search scheme, we studied the tetrahedral structures of the dilithium/disodium transition metal orthosilicates A2MSiO4 with A = Li or Na and M = Mn, Fe or Co. In addition to finding all previously reported structures, we discovered many other different tetrahedral-network-based crystal structures which are highly degenerate in energy. These structures can be classified into structures with 1D, 2D and 3D M-Si-O frameworks. A clear trend of the structural preference in different systems was revealed and possible indicators that affect the structure stabilities were introduced. For the case of Na systems which have been much less investigated in the literature relative to the Li systems, we predicted their ground state structures and found evidence for the existence of new structural motifs.

preprint2015arXiv

Nature of the Insulating Ground State of the 5d PostPerovskite CaIrO3

The insulating ground state of the 5d transition metal oxide CaIrO3 has been classified as a Mott-type insulator. Based on a systematic density functional theory (DFT) study with local, semilocal, and hybrid exchange-correlation functionals, we reveal that the Ir t2g states exhibit large splittings and one-dimensional electronic states along the c axis due to a tetragonal crystal field. Our hybrid DFT calculation adequately describes the antiferromagnetic (AFM) order along the c direction via a superexchange interaction between Ir4+ spins. Furthermore, the spin-orbit coupling (SOC) hybridizes the t2g states to open an insulating gap. These results indicate that CaIrO3 can be represented as a spin-orbit Slater insulator, driven by the interplay between a long-range AFM order and the SOC. Such a Slater mechanism for the gap formation is also demonstrated by the DFT + dynamical mean field theory calculation, where the metal-insulator transition and the paramagnetic to AFM phase transition are concomitant with each other.

preprint2015arXiv

Structures and magnetic properties of Co-Zr-B magnets studied by first-principles calculations

The structures and magnetic properties of the Co-Zr-B alloys near the Co5Zr composition were studied using adaptive genetic algorithm and first-principles calculations to guide further experimental effort on optimizing their magnetic performances. Through extensive structural searches, we constructed the contour maps of the energetics and magnetic moments of the Co-Zr-B magnet alloys as a function of composition. We found that the Co-Zr-B system exhibits the same structural motif as the "Co11Zr2" polymorphs, which plays a key role in achieving high coercivity. Boron atoms can either substitute selective cobalt atoms or occupy the interstitial sites. First-principles calculation shows that the magnetocrystalline anisotropy energies can be significantly improved through proper boron doping.

preprint2014arXiv

Energetics and Electronic Structure of Plutonium

Plutonium is the most exotic and mysterious element in the periodic table. It has 6 metallic phases and peculiar physical properties not yet understood. One of the most intriguing properties of Pu is that relatively small changes of temperature can induce transitions between different structures, that are accompanied by very large changes of equilibrium volumes. This fact has stimulated extensive theoretical and experimental studies. In spite of this, a convincing explanation of the metallurgic properties of Pu based on fundamental principles is still lacking, and none of the previous theories has been able to describe simultaneously the energetics and the $f$ electronic structure of all of the phases of Pu on the same footing. Here we provide a bird's eye view of Pu by studying the zero-temperature pressure-volume phase diagram and the $f$ electronic structure of all of its crystalline phases from first principles. In particular, we clarify the way in which the $f$-electron correlations determine its unusual energetics. Our theoretical energetics and ground-state $f$ electronic structure are both in good quantitative agreement with the experiments.

preprint2014arXiv

Interplay of spin-orbit and entropic effects in Cerium

We perform first-principle calculations of elemental cerium and compute its pressure-temperature phase diagram, finding good quantitative agreement with the experiments. Our calculations indicate that, while a signature of the volume-collapse transition appears in the free energy already at low temperatures, at larger temperatures this signature is enhanced because of the entropic effects, and originates an actual thermodynamical instability. Furthermore, we find that the catalyst determining this feature is --- in all temperature regimes --- a pressure-induced effective reduction of the $f$-level degeneracy due to the spin-orbit coupling. Our analysis suggests also that the lattice vibrations might be crucial in order to capture the behavior of the pressure-temperature transition line at large temperatures.

preprint2011arXiv

New ultrahigh pressure phases of H2O ice predicted using an adaptive genetic algorithm

We propose three new phases of H2O under ultrahigh pressure. Our structural search was performed using an adaptive genetic algorithm which allows an extensive exploration of crystal structure. The new sequence of pressure-induced transitions beyond ice X at 0 K should be ice X - Pbcm - Pbca - Pmc21 - P21 - P21/c phases. Across the Pmc21 - P21 transition, the coordination number of oxygen increases from 4 to 5 with a significant increase of density. All stable crystalline phases have nonmetallic band structures up to 7 TPa.