Unveiling the effect of Ni on the formation and structure of Earth's inner core.

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 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations 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 that Ni may substantially impact the structure and formation process of the solid inner core.

Two-step nucleation of the Earth's inner core.

The Earth's inner core started forming when molten iron cooled below the melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies have found that it requires an unrealistic degree of undercooling to nucleate the stable, hexagonal, close-packed (hcp) phase of iron that is unlikely to be reached under core conditions and age. This contradiction is referred to as the inner core nucleation paradox. 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 does 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 in solving the inner core nucleation paradox. The two-step nucleation scenario of the inner core also opens an avenue for understanding the structure and anisotropy of the present inner core.

Accelerating the discovery of novel magnetic materials using machine learning-guided adaptive feedback.

Magnetic materials are essential for energy generation and information devices, and they play an important role in advanced technologies and green energy economies. Currently, the most widely used magnets contain rare earth (RE) elements. An outstanding challenge of notable scientific interest is the discovery and synthesis of novel magnetic materials without RE elements that meet the performance and cost goals for advanced electromagnetic devices. Here, we report our discovery and synthesis of an RE-free magnetic compound, Fe3CoB2, through an efficient feedback framework by integrating machine learning (ML), an adaptive genetic algorithm, first-principles calculations, and experimental synthesis. Magnetic measurements show that Fe3CoB2 exhibits a high magnetic anisotropy (K1 = 1.2 MJ/m3) and saturation magnetic polarization (Js = 1.39 T), which is suitable for RE-free permanent-magnet applications. Our ML-guided approach presents a promising paradigm for efficient materials design and discovery and can also be applied to the search for other functional materials.

Unveiling a Family of Dimerized Quantum Magnets, Conventional Antiferromagnets, and Nonmagnets in Ternary Metal Borides.

Dimerized quantum magnets are exotic crystalline materials where Bose-Einstein condensation of magnetic excitations can happen. However, known dimerized quantum magnets are limited to only a few oxides and halides. Here, we unveil 9 dimerized quantum magnets and 11 conventional antiferromagnets in ternary metal borides MTB4 (M = Sc, Y, La, Ce, Lu, Mg, Ca, and Al; T = V, Cr, Mn, Fe, Co, and Ni), where T atoms are arranged in structural dimers. Quantum magnetism in these compounds is dominated by strong antiferromagnetic (AFM) interactions between Cr (Cr and Mn for M = Mg and Ca) atoms within the dimers, with much weaker interactions between the dimers. These systems are proposed to be close to a quantum critical point between a disordered singlet spin-dimer phase, with a spin gap, and the ordered conventional Néel AFM phase. They greatly enrich the materials inventory that allows investigations of the spin-gap phase. Conventional antiferromagnetism in these compounds is dominated by ferromagnetic Mn (Fe for M = Mg and Ca) interactions within the dimers. The predicted stable and nonmagnetic (NM) YFeB4 phase is synthesized and characterized, providing a scarce candidate to study Fe dimers and Fe ladders in borides. The identified quantum, conventional, and NM systems provide a platform with abundant possibilities to tune the magnetic exchange coupling by doping and study the unconventional quantum phase transition and conventional magnetic transitions. This work opens new avenues for studying novel magnetism in borides arising from spin dimers and establishes a theoretical workflow for future searches for dimerized quantum magnets in other families of materials.

High-Throughput Screening for Boride Superconductors.

A high-throughput screening using density functional calculations is performed to search for stable boride superconductors from the existing materials database. The workflow employs the fast frozen-phonon method as the descriptor to evaluate the superconducting properties quickly. Twenty-three stable candidates were identified during the screening. The superconductivity was obtained earlier experimentally or computationally for almost all found binary compounds. Previous studies on ternary borides are very limited. Our extensive search among ternary systems confirmed superconductivity in known systems and found several new compounds. Among these discovered superconducting ternary borides, TaMo2B2 shows the highest superconducting temperature of ∼12 K. Most predicted compounds were synthesized previously; therefore, our predictions can be examined experimentally. Our work also demonstrates that the boride systems can have diverse structural motifs that lead to superconductivity.

Structure, bonding and electronic characteristics of amorphous Se.

Ovonic threshold switching (OTS) selectors can effectively improve the storage density and suppress the leakage current of advanced phase-change memory devices. As a prototypical OTS material, amorphous GeSe is widely investigated. But the attention paid to amorphous Se (i.e., the functional constituent in amorphous GeSe) has been very limited up to now. Here we have explored the structure, bonding and electronic characteristics of amorphous Se using ab initio molecular dynamics simulations. The results reveal that the Se atoms in amorphous Se tend to form 2-coordinated configurations, and they connect with each other to form long chains. The fraction of the vibrational density of state located in the high frequency range is relatively large, and the formation energy of the Se-Se bond is as large as 4.44 eV, hinting that the Se-Se bonds in chains possess high stability. In addition, the mid-gap state related to the OTS behavior is also found in amorphous Se despite the small proportion. Our findings enrich the knowledge of amorphous Se, which aids the applications of Se-based OTS selectors.

Prediction of superconductivity in metallic boron-carbon compounds from 0 to 100 GPa by high-throughput screening.

Boron-carbon compounds have been shown to have feasible superconductivity. In our earlier paper [Zheng et al., Phys. Rev. B, 2023, 107, 014508], 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 from periods 3, 4, and 5 under pressures ranging from 0 to 100 GPa. Using the high-throughput screening method of zone-center electron-phonon interaction, we found 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 a Tc of <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.

High-Throughput Screening of Strong Electron-Phonon Couplings in Ternary Metal Diborides.

We perform a high-throughput screening on phonon-mediated superconductivity in a 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 anticorrelation between the 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.

Machine Learning-Guided Discovery of Ternary Compounds Containing La, P, and Group 14 Elements.

We integrate a deep machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary compounds. Using La-Si-P as a prototype system, we demonstrate that ML-guided first-principles calculations can efficiently explore crystal structures and their relative energetic stabilities, thus greatly accelerate the pace of material discovery. A number of new La-Si-P ternary compounds with formation energies less than 30 meV/atom above the known ternary convex hull are discovered. Among them, the formation energies of La5SiP3 and La2SiP phases are only 2 and 10 meV/atom, respectively, above the convex hull. These two compounds are dynamically stable with no imaginary phonon modes. Moreover, by replacing Si with heavier-group 14 elements in the eight lowest-energy La-Si-P structures from our ML-guided predictions, a number of low-energy La-X-P phases (X = Ge, Sn, Pb) are predicted.

Topochemical Deintercalation of Li from Layered LiNiB: toward 2D MBene.

The pursuit of two-dimensional (2D) borides, MBenes, has proven to be challenging, not the least because of the lack of a suitable precursor prone to the deintercalation. Here, we studied room-temperature topochemical deintercalation of lithium from the layered polymorphs of the LiNiB compound with a considerable amount of Li stored in between [NiB] layers (33 at. % Li). Deintercalation of Li leads to novel metastable borides (Li∼0.5NiB) with unique crystal structures. Partial removal of Li is accomplished by exposing the parent phases to air, water, or dilute HCl under ambient conditions. Scanning transmission electron microscopy and solid-state 7Li and 11B NMR spectroscopy, combined with X-ray pair distribution function (PDF) analysis and DFT calculations, were utilized to elucidate the novel structures of Li∼0.5NiB and the mechanism of Li-deintercalation. We have shown that the deintercalation of Li proceeds via a "zip-lock" mechanism, leading to the condensation of single [NiB] layers into double or triple layers bound via covalent bonds, resulting in structural fragments with Li[NiB]2 and Li[NiB]3 compositions. The crystal structure of Li∼0.5NiB is best described as an intergrowth of the ordered single [NiB], double [NiB]2, or triple [NiB]3 layers alternating with single Li layers; this explains its structural complexity. The formation of double or triple [NiB] layers induces a change in the magnetic behavior from temperature-independent paramagnets in the parent LiNiB compounds to the spin-glassiness in the deintercalated Li∼0.5NiB counterparts. LiNiB compounds showcase the potential to access a plethora of unique materials, including 2D MBenes (NiB).

How to Look for Compounds: Predictive Screening and in†situ Studies in Na-Zn-Bi System.

Here, the combination of theoretical computations followed by rapid experimental screening and in situ diffraction studies is demonstrated as a powerful strategy for novel compounds discovery. When applied for the previously "empty" Na-Zn-Bi system, such an approach led to four novel phases. The compositional space of this system was rapidly screened via the hydride route method and the theoretically predicted NaZnBi (PbClF type, P4/nmm) and Na11 Zn2 Bi5 (Na11 Cd2 Sb5 type, P 1 ‾ ) phases were successfully synthesized, while other computationally generated compounds on the list were rejected. In addition, single crystal X-ray diffraction studies of NaZnBi indicate minor deviations from the stoichiometric 1 : 1 : 1 molar ratio. As a result, two isostructural (PbClF type, P4/nmm) Zn-deficient phases with similar compositions, but distinctly different unit cell parameters were discovered. The vacancies on Zn sites and unit cell expansion were rationalized from bonding analysis using electronic structure calculations on stoichiometric "NaZnBi". In-situ synchrotron powder X-ray diffraction studies shed light on complex equilibria in the Na-Zn-Bi system at elevated temperatures. In particular, the high-temperature polymorph HT-Na3 Bi (BiF3 type, Fm 3 ‾ m) was obtained as a product of Na11 Zn2 Bi5 decomposition above 611 K. HT-Na3 Bi cannot be stabilized at room temperature by quenching, and this type of structure was earlier observed in the high-pressure polymorph HP-Na3 Bi above 0.5 GPa. The aforementioned approach of predictive synthesis can be extended to other multinary systems.

Ternary Zinc Antimonides Unlocked Using Hydride Synthesis.

Three new sodium zinc antimonides Na11Zn2Sb5, Na4Zn9Sb9, and NaZn3Sb3 were synthesized utilizing sodium hydride NaH as a reactive sodium source. In comparison to the synthesis using sodium metal, salt-like NaH can be ball-milled, leading to the easy and uniform mixing of precursors in the desired stoichiometric ratios. Such comprehensive compositional control enables a fast screening of the Na-Zn-Sb system and identification of new compounds, followed by their preparation in bulk with high purity. Na11Zn2Sb5 crystallizes in the triclinic P1 space group (No. 2, Z = 2, a = 8.8739(6) Å, b = 10.6407(7) Å, c = 11.4282(8) Å, α = 103.453(2)°, ß = 96.997(2)°, γ = 107.517(2)°) and features polyanionic [Zn2Sb5]11- clusters with unusual 3-coordinated Zn atoms. Both Na4Zn9Sb9 (Z = 4, a = 28.4794(4) Å, b = 4.47189(5) Å, c = 17.2704(2) Å, ß = 98.3363(6)°) and NaZn3Sb3 (Z = 8, a = 32.1790(1) Å, b = 4.51549(1) Å, c = 9.64569(2) Å, ß = 98.4618(1)°) crystallize in the monoclinic C2/m space group (No. 12) and have complex new structure types. For both compounds, their frameworks are built from ZnSb4 distorted tetrahedra, which are linked via edge-, vertex-sharing, or both, while Na cations fill in the framework channels. Due to the complex structures, Na4Zn9Sb9 and NaZn3Sb3 compounds exhibit low thermal conductivities (0.97-1.26 W·m-1 K-1) at room temperature, positive Seebeck coefficients (19-32 µV/K) suggestive of holes as charge carriers, and semimetallic electrical resistivities (∼1.0-2.3 × 10-4 Ω·m). Na4Zn9Sb9 and NaZn3Sb3 decompose into the equiatomic NaZnSb above ∼800 K, as determined by in situ synchrotron powder X-ray diffraction. The discovery of multiple ternary compounds highlights the importance of judicious choice of the synthetic method.

Origin of short- and medium-range order in supercooled liquid Ge3Sb2Te6 using ab initio molecular dynamics simulations.

Phase-change materials such as Ge-Sb-Te compounds have attracted much attention due to their potential value in electrical data storage. In contrast to the amorphous and crystalline phases, supercooled liquids are far from being deeply understood despite their inevitable role in both amorphization and crystallization processes. To this end, we have studied the dynamics properties and structural characteristics of liquid and supercooled liquid Ge3Sb2Te6 during the fast cooling process. As the temperature decreases, chemical bonds become more homogeneous, but coordination numbers of Ge, Sb and Te atoms change very little. Meanwhile, the structural order of short-range configuration is obviously enhanced. Further studies suggest that Ge-centered, Sb-centered and Te-centered configurations change to the more ordered defective octahedrons mainly by adjusting the bond-angle relationship and bond length, rather than just by changing the coordination environment. It is the more ordered octahedrons that promote the formation of medium-range order. Our findings provide a deep insight into the origin of local structural order in supercooled liquid Ge3Sb2Te6, which is of great importance for the comprehensive understanding of amorphization and crystallization processes.

Stabilizing the crystal structures of NaFePO4 with Li substitutions.

Due to the high cost and insufficient resources of lithium, alternative sodium-ion batteries have been widely investigated for large-scale applications. NaFePO4 has the highest theoretical capacity of 154 mA h g-1 among the iron-based phosphates, which makes it an attractive cathode material for Na-ion batteries. Experimentally, LiFePO4 has been highly successful as a cathode material in Li-ion batteries because its olivine crystal structure provides a stable framework during battery cycling. In NaFePO4, maricite replaces olivine as the most stable phase. However, the maricite phase is experimentally found to be electrochemically inactive under normal battery operating voltages (0-4.5 V). We found that partial substitutions of Na with Li stabilize the olivine structure and may be a way to improve the performance of NaFePO4 cathodes. Using the previously developed structural LiFePO4 database, we examined the low-energy crystal structures in the system when we replace Li with Na. The known maricite and olivine NaFePO4 phases are reconfirmed and an unreported phase with energy between them is identified by our calculations. Besides, the Li-doped olivine type compound LixNa1-xFePO4 with mixed alkali ions retains better energetic stability compared with the other two types of structures of the same composition, as long as the proportion of Li exceeds 0.25. The thermodynamic stability of o-type LixNa1-xFePO4 can be further improved at finite temperatures. The primary limitation of the calculations is that we mainly focus on the zero-temperature condition; however, the relative stability of the structures may vary depending on the ambient temperature.

Predicting Complex Relaxation Processes in Metallic Glass.

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_{90}Sm_{10} metallic glass would have complex relaxation behaviors: In addition to the main (α) relaxation, the glass (i) shows a pronounced secondary (ß) relaxation at cryogenic temperatures and (ii) exhibits an anomalous relaxation process (α_{2}) accompanying α relaxation. Both of the predictions are verified by experiments. Computational simulations reveal the microscopic origins of relaxation processes: while the pronounced ß relaxation is attributed to the abundance of stringlike cooperative atomic rearrangements, the anomalous α_{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.

Computationally Driven Discovery of a Family of Layered LiNiB Polymorphs.

Two novel lithium nickel boride polymorphs, RT-LiNiB and HT-LiNiB, with layered crystal structures are reported. This family of compounds was theoretically predicted by using the adaptive genetic algorithm (AGA) and subsequently synthesized by a hydride route with LiH as the lithium source. Unique among the known ternary transition-metal borides, the LiNiB structures feature Li layers alternating with nearly planar [NiB] layers composed of Ni hexagonal rings with a B-B pair at the center. A comprehensive study using a combination of single crystal/synchrotron powder X-ray diffraction, solid-state 7 Li and 11 B NMR spectroscopy, scanning transmission electron microscopy, quantum-chemical calculations, and magnetism has shed light on the intrinsic features of these polymorphic compounds. The unique layered structures of LiNiB compounds make them ultimate precursors for exfoliation studies, thus paving a way toward two-dimensional transition-metal borides, MBenes.

Overcoming the Time Limitation in Molecular Dynamics Simulation of Crystal Nucleation: A Persistent-Embryo Approach.

The crystal nucleation from liquid in most cases is too rare to be accessed within the limited time scales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. Thus, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.

Calculation of critical nucleation rates by the persistent embryo method: application to quasi hard sphere models.

We study the 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 comparison with other existing predictions. Implications for experiments are also discussed.

Fe-Si networks and charge/discharge-induced phase transitions in Li2FeSiO4 cathode materials.

Structural phase transitions of electrode materials are responsible for poor reversibility during charge/discharge cycling in Li-ion batteries. Using previously developed structural databases, we investigate a structural landscape for LixFeSiO4 systems at x = 1. Starting with low-energy Li2FeSiO4 crystal structures, we explore the crystal structures of the material in different states of charge. The as-prepared Li2FeSiO4 materials adopt low energy structures characterized by two-dimensional (2D) Fe-Si networks. After the removal of one Li per formula unit to form LiFeSiO4, the structures with three-dimensional (3D) diamond-like Fe-Si networks become more energetically favorable without a significant impact on the charge capacity, which agrees with previous experimental and theoretical work. However, we reveal that the structure with a 3D diamond-like Fe-Si network can further transform into a new structure at x = 1. And the Li atom is hard to reinsert into these new structures. Consequently the system is prevented from returning to the Li2FeSiO4 state. We believe that the formation of this new structure plays an important role in the loss of reversible capacity of Li2FeSiO4 electrode materials.