Search

Electronic structure analysis is performed to study the properties of several Pb-containing 3d-intermetallics. Our study reveals that binary metastable Co3Pb and Fe3Pb intermetallic compounds exhibit very attractive intrinsic magnetic properties. We primarily focus on the magnetic anisotropic properties arising from the high spin-orbit coupling of the Pb atom. Decomposing the total anisotropy into intra- and interatomic contributions reveals a significant deviation from single ion anisotropy model with strong symmetric anisotropic pair interactions present. Furthermore, we consider magnetic properties of ternary Pb-based 3d-intermetallics which recently have been reported as stable or metastable. Giant magnetic anisotropy is found on Pb atoms in these systems. The origin of such strong anisotropy in La18Co28Pb3 appears from two sources: spin-orbit and interelectronic Breit couplings. The significance of Breit interaction for magnetic anisotropy in bulk systems has not been reported previously. It is expected that Breit coupling induced anisotropy is dominating in magnetic Pb-based magnets with lower dimensionality including thin films., Comment: 21 pages, 9 figures

The adaptive variational quantum dynamics simulation (AVQDS) method performs real-time evolution of quantum states using automatically generated parameterized quantum circuits that often contain substantially fewer gates than Trotter circuits. Here we report an improved version of the method, which we call AVQDS(T), by porting the Tiling Efficient Trial Circuits with Rotations Implemented Simultaneously (TETRIS) technique. The algorithm adaptively adds layers of disjoint unitary gates to the ansatz circuit so as to keep the McLachlan distance, a measure of the accuracy of the variational dynamics, below a fixed threshold. We perform benchmark noiseless AVQDS(T) simulations of quench dynamics in local spin models demonstrating that the TETRIS technique significantly reduces the circuit depth and two-qubit gate count. We also show a method based on eigenvalue truncation to solve the linear equations of motion for the variational parameters with enhanced noise resilience. Finally, we propose a way to substantially alleviate the measurement overhead of AVQDS(T) while maintaining high accuracy by synergistically integrating quantum circuit calculations on quantum processing units with classical calculations using, e.g., tensor networks to evaluate the quantum geometric tensor. We showcase that this approach enables AVQDS(T) to deliver more accurate results than simulations using a fixed ansatz of comparable final depth for a significant time duration with fewer quantum resources.

The discovery of novel quantum materials within ternary phase spaces containing antagonistic pair such as Fe with Bi, Pb, In, and Ag, presents significant challenges yet holds great potential. In this work, we investigate the stabilization of these immiscible pairs through the integration of Cerium (Ce), an abundant rare-earth and cost-effective element. By employing a machine learning (ML)-guided framework, particularly crystal graph convolutional neural networks (CGCNN), combined with first-principles calculations, we efficiently explore the composition/structure space and predict 9 stable and 37 metastable Ce-Fe-X (X=Bi, Pb, In and Ag) ternary compounds. Our findings include the identification of multiple new stable and metastable phases, which are evaluated for their structural and energetic properties. These discoveries not only contribute to the advancement of quantum materials but also offer viable alternatives to critical rare earth elements, underscoring the importance of Ce-based intermetallic compounds in technological applications.

Constraining the melting temperature of iron under Earth's inner core conditions is crucial for understanding core dynamics and planetary evolution. Here, we develop a deep potential (DP) model for iron that explicitly incorporates electronic entropy contributions governing thermodynamics under Earth's core conditions. Extensive benchmarking demonstrates the DP's high fidelity across relevant iron phases and extreme pressure and temperature conditions. Through thermodynamic integration and direct solid-liquid coexistence simulations, the DP predicts melting temperatures for iron at the inner core boundary, consistent with previous \textit{ab initio} results. This resolves the previous discrepancy of iron's melting temperature at ICB between the DP model and \textit{ab initio} calculation and suggests the crucial contribution of electronic entropy. Our work provides insights into machine learning melting behavior of iron under core conditions and provides the basis for future development of binary or ternary DP models for iron and other elements in the core.

We outline how pairs of strongly immiscible elements, referred to here as antagonistic pairs, can be used to synthesize ternary compounds with quasi-reduced dimensional motifs. By identifying third elements that are compatible with a given antagonistic pair, ternary compounds can be formed in which the third element segregates the immiscible atoms into spatially separated substructures. Quasi-low dimensional structural units are a natural consequence of the immiscible atoms seeking to avoid contact in the solid-state. As proof of principle, we present the discovery and physical properties of La$_4$Co$_4$X (X = Pb, Bi, Sb), a new family of intermetallics based on the antagonistic pairs Co-Pb and Co-Bi. La$_4$Co$_4$X adopts a new orthorhombic crystal structure (space group Pbam) containing quasi-2D Co slabs and La-X layers that stack along the a-axis. Consistent with our proposal, the La atoms separate the Co and X substructures, ensuring there are no direct contacts between immiscible atoms. Within the Co slabs, the atoms occupy the vertices of corner sharing tetrahedra and triangles, and this motif produces flat electronic bands near the Fermi level that favor magnetism. The Co is moment bearing in La$_4$Co$_4$X, and we show that whereas La$_4$Co$_4$Pb behaves as a three dimensional antiferromagnet with T$_N$ = 220 K, La$_4$Co$_4$Bi and La$_4$Co$_4$Sb have behavior consistent with low dimensional magnetic coupling and ordering, with T$_N$ = 153 K and 143 K respectively. In addition to the Pb, Bi, and Sb based La$_4$Co$_4$X compounds, we were likely able to produce an analogous La$_4$Co$_4$Sn in polycrystalline form, although we were unable to isolate single crystals. We anticipate that using mutually compatible third elements with an antagonistic pair represents a generalizable design principle for discovering new materials and structure types containing low-dimensional substructures.

The electronic structure and magnetic properties of a newly predicted stable ternary compound La$_{18}$Co$_{28}$Pb$_{3}$ are studied using electronic structure analysis. The ground state of this compound is ferromagnetic, with three positions of nonequivalent magnetic Co atoms. A strong dependence of magnetic properties on volume shows that this system is situated near the point of magnetic instability. A coexistence of high- and low-spin ferromagnetic states as a function of volume near equilibrium was discovered. A corresponding spin tunneling splitting was estimated. The stability of the theoretically predicted magnetic ground state was tested by varying the Hubbard parameter. The thermal spin fluctuations were added to estimate the paramagnetic moment and a Curie temperature. The necessity of experimental verification of the obtained results is emphasized.

Exploring high-temperature superconducting (high-$T_c$) material at ambient pressure holds immense significance for physics, chemistry, and materials science. In this study, we perform a high-throughput screening of strong electron-phonon interactions in X$_2$MH$_6$ compounds (X = Li, Na, Mg, Al, K, Ca, Ga, Rb, Sr, and In; M are $3d$, $4d$, and $5d$ transition metals). These compounds have a cubic structure featuring an MH$_6$ octahedron motif. Our screening calculations suggest that 26 compounds exhibit dynamic stability and strong electron-phonon coupling. Among these 26 compounds, Mg$_2$RhH$_6$, Mg$_2$IrH$_6$, Al$_2$MnH$_6$, and Li$_2$CuH$_6$ show promising energetic stability and $T_c$ of more than 50 K at ambient pressure. This study underscores promising high-$T_c$ compounds at ambient pressure with distinctive MH$_6$ motifs.

We applied the ab initio spin-polarized Correlation Matrix Renormalization Theory (CMRT) to the ferromagnetic state of the bulk BCC iron. We showed that it was capable of reproducing the equilibrium physical properties and the pressure-volume curve in good comparison with experiments. We then focused on the analysis of its local electronic correlations. By exploiting different local fluctuation-related physical quantities as measures of electronic correlation within target orbits, we elucidated the different roles of $t_{2g}$ and $e_g$ states in both spin channels and presented compelling evidence to showcase this distinction in their electronic correlation., Comment: 11 pages, 6 figures, regular article

We apply the adaptive variational quantum imaginary time evolution (AVQITE) method to prepare ground states of one-dimensional spin $S=1$ models. We compare different spin-to-qubit encodings (standard binary, Gray, unary, and multiplet) with regard to the performance and quantum resource cost of the algorithm. Using statevector simulations we study two well-known spin-1 models: the Blume-Capel model of transverse-field Ising spins with single-ion anisotropy, and the XXZ model with single-ion anisotropy. We consider system sizes of up to $20$ qubits, which corresponds to spin-$1$ chains up to length $10$. We determine the dependence of the number of CNOT gates in the AVQITE state preparation circuit on the encoding, the initial state, and the choice of operator pool in the adaptive method. Independent on the choice of encoding, we find that the CNOT gate count scales cubically with the number of spins for the Blume-Capel model and quartically for the anistropic XXZ model. However, the multiplet and Gray encodings present smaller prefactors in the scaling relations. These results provide useful insights for the implementation of AVQITE on quantum hardware., Comment: 11 pages, 6 figures

Boron carbon compounds have been shown to have feasible superconductivity. In our earlier paper [Zheng et al., Phys. Rev. B 107, 014508 (2023)], we identified a new conventional superconductor of LiB3C at 100 GPa. Here, we aim to extend the investigation of possible superconductivity in this structural framework by replacing Li atoms with 27 different cations under pressures ranging from 0 to 100 GPa. Using the high-throughput screening method of zone-center electron-phonon interaction, we find that ternary compounds like CaB3C, SrB3C, TiB3C, and VB3C are promising candidates for superconductivity. The consecutive calculations using the full Brillouin zone confirm that they have Tc < 31 K at moderate pressures. Our study demonstrates that fast screening of superconductivity by calculating zone-center electron-phonon coupling strength is an effective strategy for high-throughput identification of new superconductors.

Ni is the second most abundant element in the Earth's core. Yet, its effects on the inner core's structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe's melting point at inner core pressures. The melting temperature of Ni is shown to be 700-800 K higher than that of Fe at 323-360 GPa. hcp, bcc, and liquid phase relation differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe's crystallization at core pressures. These results suggest Ni may substantially impact the structure and formation process of the solid inner core.

Terrestrial exoplanets are of great interest for being simultaneously similar to and different from Earth. Their compositions are likely comparable to those of solar-terrestrial objects, but their internal pressures and temperatures can vary significantly with their masses/sizes. The most abundant non-volatile elements are O, Mg, Si, Fe, Al, and Ca, and there has been much recent progress in understanding the nature of magnesium silicates up to and beyond ~3 TPa. However, a critical element, Fe, has yet to be systematically included in materials discovery studies of potential terrestrial planet-forming phases at ultra-high pressures. Here, using the adaptive genetic algorithm (AGA) crystal structure prediction method, we predict several unreported stable crystalline phases in the binary Fe-Mg and ternary Fe-Mg-O systems up to pressures of 3 TPa. The analysis of the local packing motifs of the low-enthalpy Fe-Mg-O phases reveals that the Fe-Mg-O system favors a BCC motif under ultra-high pressures regardless of chemical composition. Besides, oxygen enrichment is conducive to lowering the enthalpies of the Fe-Mg-O phases. Our results extend the current knowledge of structural information of the Fe-Mg-O system to exoplanet pressures.

Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe$_5$. We show that a coherent $A_\text{1g}$ Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau-Zener-St\"uckelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe$_5$ at low temperature., Comment: 15 pages, 6 figures

Pursuit of superconductivity in light-element systems at ambient pressure is of great experimental and theoretical interest. In this work, we combine a machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary Ca-B-C compounds. Three new layered borocarbides (stable CaBC5 and metastable Ca2BC11 and CaB3C3) are predicted to be phonon-mediated superconductors at ambient pressure. The hexagonal CaB3C3 possesses the highest Tc of 26.05 K among the three compounds. The {\sigma}-bonging bands around the Fermi level account for the large electron-phonon coupling ({\lambda} = 0.980) of hexagonal CaB3C3. The ML-guided approach opens up a way for greatly accelerating the discovery of new high-Tc superconductors.

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of Fe90Si3O7 liquid close to the Earth's liquid core conditions. The simulation results reveal that the short-range structural order (SRO) in the Fe90Si3O7 liquid is very strong. About 80% of the atoms are arranged in crystalline-like SRO motifs. In particular, ~70% of Fe-centered clusters can be classified as either hexagonal-close-pack (HCP/HCP-like) or icosahedral (ICO/ICO-like) SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O diffusion coefficient is about 2-3 times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

Van der Waals (vdW) heterostructures, which are produced by the precise assemblies of varieties of two-dimensional (2D) materials, have demonstrated many novel properties and functionalities. Here we report a nano-plasmonic study of vdW heterostructures that were produced by depositing ordered molecular layers of pentacene on top of graphene. We find through nano-infrared (IR) imaging that surface plasmons formed due to the collective oscillations of Dirac fermions in graphene are highly sensitive to the adjacent pentacene layers. In particular, the plasmon wavelength declines systematically but nonlinearly with increasing pentacene thickness. Further analysis and density functional theory (DFT) calculations indicate that the observed peculiar thickness dependence is mainly due to the tunneling-type electron transfer from pentacene to graphene. Our work unveils a new method for tailoring graphene plasmons and deepens our understanding of the intriguing nano-optical phenomena due to interlayer couplings in novel vdW heterostructures., Comment: 21 pages

We report a nano-infrared (IR) imaging study of trilayer graphene (TLG) with both ABA (Bernal) and ABC (rhombohedral) stacking orders using the scattering-type scanning near-field optical microscope (s-SNOM). With s-SNOM operating in the mid-IR region, we mapped in real space the surface plasmon polaritons (SPPs) of ABA-TLG and ABC-TLG, which are tunable with electrical gating. Through quantitative modeling of the plasmonic imaging data, we found that the plasmon wavelength of ABA-TLG is significantly larger than that of ABC-TLG, resulting in a sizable impedance mismatch and hence a strong plasmon reflection at the ABA/ABC lateral junction. Further analysis indicates that the different plasmonic responses of the two types of TLG are directly linked to their electronic structures and carrier properties. Our work uncovers the physics behind the stacking-dependent plasmonic responses of TLG and sheds light on future applications of TLG and the ABA/ABC junctions in IR plasmonics and planar nano-optics., Comment: 11 pages

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.

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., Comment: 8 pages, 5 figures

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.

There has been a long debate on the stable phase of iron under the Earth's inner core conditions. Because of the solid-liquid coexistence at the inner core boundary, the thermodynamic stability of solid phases directly relates to their melting temperatures, which remains considerable uncertainty. In the present study, we utilized a semi-empirical potential fitted to high-temperature ab initio data to perform a thermodynamic integration from classical systems described by this potential to ab initio systems. This method provides a smooth path for thermodynamic integration and significantly reduces the uncertainty caused by the finite-size effect. Our results suggest the hcp phase is the stable phase of pure iron under the inner core conditions, while the free energy difference between the hcp and bcc phases is tiny, on the order of 10s meV/atom near the melting temperature., Comment: 4 figures

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.

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., Comment: 7 pages, 9 figures

Oxygen solubility in solid iron is extremely low, even at high pressures and temperatures. Thus far, no Fe-O compounds between Fe and FeO endmembers have been reported experimentally. We observed chemical reactions of Fe with FeO or Fe$_2$O$_3$ $in\ situ$ x-ray diffraction experiments at 220-260 GPa and 3,000-3,500 K. The refined diffraction patterns are consistent with a series of Fe$_n$O (n $>$ 1) compounds (e.g., Fe$_{25}$O$_{13}$ and Fe$_{28}$O$_{14}$) identified using the adaptive genetic algorithm. Like $\epsilon$-Fe in the hexagonal close-packed (hcp) structure, the structures of Fe$_n$O compounds consist of oxygen-only close-packed monolayers distributed between iron-only layers. $Ab\ initio$ calculations show systematic electronic properties of these compounds that have ramifications for the physical properties of Earth's inner core., Comment: 6 pages, 4 figures

Being a lithophile element at ambient pressure, magnesium is long believed to be immiscible with iron. A recent study by Gao et al. [1] showed that pressure turns magnesium into a siderophile element and can produce unconventional Fe-Mg compounds. Here, we extend the investigation to exoplanetary pressure conditions using an adaptive genetic algorithm-based variable-composition structural prediction approach. We identify several Fe-Mg phases up to 3 TPa. Our cluster alignment analysis reveals that most of the predicted Fe-Mg compounds prefer a BCC packing motif at terapascal pressures. This study provides a more comprehensive structure database to support future investigations of the high-pressure structural behavior of Fe-Mg and ternary, quaternary, etc. compounds involving these elements., Comment: 6 pages, 7 figures

Fe, Mg, and O are among the most abundant elements in terrestrial planets. While the behavior of the Fe-O, Mg-O, and Fe-Mg binary systems under pressure have been investigated, there are still very few studies of the Fe-Mg-O ternary system at relevant Earth's core and super-Earth's mantle pressures. Here, we use the adaptive genetic algorithm (AGA) to study ternary Fe$_x$Mg$_y$O$_z$ phases in a wide range of stoichiometries at 200 GPa and 350 GPa. We discovered three dynamically stable phases with stoichiometries FeMg$_2$O$_4$, Fe$_2$MgO$_4$, and FeMg$_3$O$_4$ with lower enthalpy than any known combination of Fe-Mg-O high-pressure compounds at 350 GPa. With the discovery of these phases, we construct the Fe-Mg-O ternary convex hull. We further clarify the composition- and pressure-dependence of structural motifs with the analysis of the AGA-found stable and metastable structures. Analysis of binary and ternary stable phases suggest that O, Mg, or both could stabilize a BCC iron alloy at inner core pressures., Comment: 9 pages, 8 figures

An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near-term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H$_4$, H$_2$O and BeH$_2$ molecules, where it yields compact variational ans\"atze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is demonstrated through a set of AVQITE calculations of quantum spin models. Finally, it is shown that quantum Lanczos calculations can also be naturally performed alongside AVQITE without additional quantum resource costs., Comment: 16 pages, 4 figures

We propose a general-purpose, self-adaptive approach to construct variational wavefunction ans\"atze for highly accurate quantum dynamics simulations based on McLachlan's variational principle. The key idea is to dynamically expand the variational ansatz along the time-evolution path such that the ``McLachlan distance'', which is a measure of the simulation accuracy, remains below a set threshold. We apply this adaptive variational quantum dynamics simulation (AVQDS) approach to the integrable Lieb-Schultz-Mattis spin chain and the nonintegrable mixed-field Ising model, where it captures both finite-rate and sudden post-quench dynamics with high fidelity. The AVQDS quantum circuits that prepare the time-evolved state are much shallower than those obtained from first-order Trotterization and contain up to two orders of magnitude fewer CNOT gate operations. We envision that a wide range of dynamical simulations of quantum many-body systems on near-term quantum computing devices will be made possible through the AVQDS framework., Comment: 12 pages, 7 figures

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 $\Gamma$ point in agreement with DFT calculations. The presence of Rashba splitting may render this material of potential interest for spintronic applications.

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., Comment: 11 pages, 5 figures

We propose a scheme of decomposition of the total relativistic energy in solids to intra- and interatomic contributions. The method is based on a variation of the speed of light from its value in relativistic theory to infinity (a non-relativistic limit). As an illustration of the method, we tested such decomposition in the case of a spin-orbit interaction variation for decomposition of the magnetic anisotropy energy (MAE) in CoPt. We further studied the {\alpha}''-Fe16N2 magnet doped by Bi, Sb, Co and Pt atoms. It has been found that the addition of Pt atoms can enhance the MAE by as large as five times while Bi and Sb substitutions double the total MAE. Using the proposed technique we demonstrate the spatial distribution of these enhancements. Our studies also suggest that Sb, Pt and Co substitutions could be synthesized by experiments.

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.

A first principle prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e. motifs, which provide a generalization of the four conventional Frank-Kasper polyhedral units. Our results show that meta-stable phases share the very same motifs as equilibrium phases. We discuss the connection with packing models, phase diagrams with repulsive potentials and the prediction of likely experimental superlattices.

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.

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., Comment: 14 pages, 5 figures

Relaxation processes significantly influence the properties of glass materials. However, understanding their specific origins is difficult, even more challenging is to forecast them theoretically. In this study, using microseconds molecular dynamics simulations together with an accurate many-body interaction potential, we predict that an Al$_{\text{90}}$Sm$_{\text{10}}$ metallic glass would have complex relaxation behaviors: In addition to the main ($\alpha$) relaxation, the glass (i) shows a pronounced secondary ($\beta$) relaxation at cryogenic temperatures and (ii) exhibits an anomalous relaxation process ($\alpha_2$) accompanying $\alpha$ relaxation. Both of the predictions are verified by experiments. Computational simulations reveal the microscopic origins of relaxation processes: while the pronounced $\beta$ relaxation is attributed to the abundance of string-like cooperative atomic rearrangements, the anomalous $\alpha_2$ process is found to correlate with the decoupling of the faster motions of Al with slower Sm atoms. The combination of simulations and experiments represents a first glimpse of what may become a predictive routine and integral step for glass physics.

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.

We explore an opportunity to induce and control tetragonal distortion in materials. The idea involves formation of a binary alloy from parent compounds having body-centered and face-centered symmetries. The concept is illustrated in the case of FeNi$_{1-x}$Co$_x$ magnetic alloy formed by substitutional doping of the L1$_0$ FeNi magnet with Co. Using electronic structure calculations we demonstrate that the tetragonal strain in this system can be controlled by concentration and it reaches maximum for $x=0.5$. This finding is then applied to create a large magnetocrystalline anisotropy (MAE) in FeNi$_{1-x}$Co$_x$ system by considering an interplay of the tetragonal distortion with electronic concentration and chemical anisotropy. In particular, we identify a new ordered FeNi$_{0.5}$Co$_{0.5}$ system with MAE larger by a factor 4.5 from the L1$_0$ FeNi magnet., Comment: 5 pages, 5 figures

Functional manipulation of graphene is an important topic in view of both fundamental researches and practical applications. In this study, we show that intercalation of 5$d$ transition metals in epitaxial graphene on SiC is a promising approach to realize topologically nontrivial phases with a finite band gap in graphene. Using first-principles calculations based on density functional theory, we show that the Re- and Ta-intercalated graphene become two-dimensional topological insulators which exhibit linear Dirac cones and quadratic bands with topological band gaps, respectively. The appearance of the topological states is attributed to the strong spin-orbit coupling strength of the intercalants. We find that topological edge states exist within the finite bulk band gap in accordance with the bulk-boundary correspondence. We also discuss the spin splitting of the band structure due to the inversion symmetry breaking and the spin-orbit coupling. Our results demonstrate that intercalation of graphene is an effective and viable method to manipulate the band gap and the topological character of graphene. Such intercalated graphene systems are potentially useful for spintronics and quantum computing applications., Comment: 11 pages, 4 figures

The formation possibility of a new (Zr0.25Nb0.25Ti0.25V0.25)C high-entropy ceramic (ZHC-1) was first analyzed by the first-principles calculations and thermodynamical analysis and then it was successfully fabricated by hot pressing sintering technique. The first-principles calculation results showed that the mixing enthalpy of ZHC-1 was 5.526 kJ/mol and the mixing entropy of ZHC-1 was in the range of 0.693R-1.040R. The thermodynamical analysis results showed that ZHC-1 was thermodynamically stable above 959 K owing to its negative mixing Gibbs free energy. The experimental results showed that the as-prepared ZHC-1 (95.1% relative density) possessed a single rock-salt crystal structure, some interesting nanoplate-like structures and high compositional uniformity from nanoscale to microscale. By taking advantage of these unique features, compared with the initial metal carbides (ZrC, NbC, TiC and VC), it showed a relatively low thermal conductivity of 15.3 + - 0.3 W/(m.K) at room temperature, which was due to the presence of solid solution effects, nanoplates and porosity. Meanwhile, it exhibited the relatively high nanohardness of 30.3 + - 0.7 GPa and elastic modulus of 460.4 + - 19.2 GPa and the higher fracture toughness of 4.7 + - 0.5 MPa.m1/2, which were attributed to the solid solution strengthening mechanism and nanoplate pullout and microcrack deflection toughening mechanism., Comment: 49 pages,6 figures, 4 tables

We investigated the homogenous nucleation of the stoichiometric B2 and B33 phases in the Ni50Zr50 alloy using the persistent embryo method and the classical nucleation theory. The two phases become very close competitors at large supercoolings, which is consistent with the experimental observations. In the case of the B2 phase, the linear temperature dependence of the solid-liquid interface (SLI) free energy extrapolated to the melting temperature leads to the same value as the one obtained from the capillarity fluctuation method (CFM). In the case of the B33 phases, the SLI free energy is also a linear function of temperature at large supercoolings but the extrapolation to the melting temperature leads to a value which is considerably different from the CFM value. This is consistent with the large anisotropy of the SLI properties of the B33 phase nearby the melting temperature observed in the simulation of the nominally flat interface migration.

The effect of dopants on the metallic glass forming ability is usually considered based on analysis of changes in the liquid structure or thermodynamics. What is missing in such considerations is an analysis of how a dopant changes the properties of the crystal phases which can form instead of the glass. In order to illuminate this aspect we performed molecular dynamics simulations to study the effects of Mg and Sm dopants on the crystal nucleation in Al. The simulation data were found to be consistent with the experimental observations that addition of Mg to Al does not lead to vitrification but addition of only 8% Sm does. The significant effect of Sm doping was related to the intolerance of Al to this dopant. This leads to increase in the solid-liquid interfacial free energy, and therefore, to increase in the nucleation barrier and to dramatic decrease in the nucleation rate. The intolerance mechanism also significantly affects the growth kinetics.

The temperature dependence of the solid-liquid interfacial free energy, {\gamma}, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed {\gamma} shows a linear dependence on the temperature. The values of {\gamma} extrapolated to the melting temperature agree well with previous data obtained by the capillary fluctuation method. Using the temperature dependence of {\gamma}, we estimate the nucleation free energy barrier in a wide temperature range from the classical nucleation theory. The obtained data agree very well with the results from the brute-force molecular dynamics simulations.

We study crystal nucleation of the Weeks-Chandler-Andersen (WCA) model, using the recently introduced Persistent Embryo Method (PEM). The method provides detailed characterization of pre-critical, critical and post-critical nuclei, as well as nucleation rates that compare favorably with those obtained using other methods (umbrella sampling, forward flux sampling or seeding). We further map our results to a hard sphere model allowing to compare with other existing predictions. Implications for experiments are also discussed., Comment: 27 pages, 11 figures

Measurements of temperature-dependent resistance and magnetization under hydrostatic pressures up to 2.13 GPa are reported for single-crystalline, superconducting BaBi$_3$. A temperature - pressure phase diagram is determined and the results suggest three different superconducting phases ${\alpha}$, $\beta$, and $\gamma$ in the studied pressure range. We further show that occurrence of the three superconducting phases is intuitively linked to phase transitions at higher temperature which are likely first order in nature. $T_p$, which separates phase ${\alpha}$ from $\beta$ and $\gamma$, is associated with an abrupt resistance change as pressure is increased from 0.27 GPa to 0.33 GPa. Above 0.33 GPa, an "S-shape" anomaly in the temperature-dependent resistance curve, $T_\text S$, is observed and associated with the transition between the $\beta$ and $\gamma$ phases. Further increasing of pressure above 1.05 GPa suppresses this transition and BaBi$_3$ stays in $\gamma$ phase over the whole investigated temperature range. These high-temperature anomalies are likely related to structural degrees of freedom. With the ${\alpha}$ phase being the ambient-pressure tetragonal structure ($P4/mmm$), our first-principle calculations suggest the $\beta$ phase to be cubic structure ($Pm-3m$) and the $\gamma$ phase to be a distorted tetragonal structure where the Bi atoms are moved out of the face-centered position. Finally, an analysis of the evolution of the superconducting upper critical field with pressure further confirms these transitions in the superconducting state and suggests a possible change of band structure or a Lifshitz transition near 1.54 GPa in $\gamma$ phase. Given the large atomic numbers of both Ba and Bi, our results establish BaBi$_3$ as a good candidate for the study of the interplay of structure with superconductivity in the presence of strong spin-orbit coupling., Comment: 11 pages, 10 figures

Superconductors with the A15 structure are prototypical type-II s-wave superconductors which have generated considerable interest in early superconducting material history. However, the topological nature of the electronic structure remains unnoticed so far. Here, using first-principles band structure calculations based on density functional theory, we show that the A15 superconductors (Ta$_3$Sb, Ta$_3$Sn, and Ta$_3$Pb) have nontrivial band topology in the bulk electronic band structures, leading to the formation of topological surface states near the Fermi energy. Due to the bulk superconductivity, the proximity effect in the topological surface states would induce topological superconductivity even without heterostructure of a topological insulator and an s-wave superconductor. Our results indicate that the A15 superconductors are promising candidates for the realization of the topological superconductivity and the Majorana fermion., Comment: 11 pages, 4 figures