First-principles theory of short-range order, electronic excitations, and spin polarization in Ni-V and Pd-V alloys

Imperfections are not 0 K: free energy of point defects in crystals

Defects determine many important properties and applications of materials, ranging from doping in semiconductors, to conductivity in mixed ionic-electronic conductors used in batteries, to active sites in catalysts. The theoretical description of defect formation in crystals has evolved substantially over the past century. Advances in supercomputing hardware, and the integration of new computational techniques such as machine learning, provide an opportunity to model longer length and time-scales than previously possible. In this Tutorial Review, we cover the description of free energies for defect formation at finite temperatures, including configurational (structural, electronic, spin) and vibrational terms. We discuss challenges in accounting for metastable defect configurations, progress such as machine learning force fields and thermodynamic integration to directly access entropic contributions, and bottlenecks in going beyond the dilute limit of defect formation. Such developments are necessary to support a new era of accurate defect predictions in computational materials chemistry.

Unlocking the hidden chemical space in cubic-phase garnet solid electrolyte for efficient quasi-all-solid-state lithium batteries

Garnet-type Li₇La₃Zr₂O₁₂ (LLZO) solid electrolytes (SE) demonstrates appealing ionic conductivity properties for all-solid-state lithium metal battery applications. However, LLZO (electro)chemical stability in contact with the lithium metal electrode is not satisfactory for developing practical batteries. To circumvent this issue, we report the preparation of various doped cubic-phase LLZO SEs without vacancy formation (i.e., Li = 7.0 such as Li₇La₃Zr_0.5Hf_0.5Sc_0.5Nb_0.5O₁₂ and Li₇La₃Zr_0.4Hf_0.4Sn_0.4Sc_0.4Ta_0.4O₁₂). The entropy-driven synthetic approach allows access to hidden chemical space in cubic-phase garnet and enables lower solid-state synthesis temperature as the cubic-phase nucleation decreases from 750 to 400 °C. We demonstrate that the SEs with Li = 7.0 show better reduction stability against lithium metal compared to SE with low lithium contents and identical atomic species (i.e., Li = 6.6 such as Li_6.6La₃Zr_0.4Hf_0.4Sn_0.4Sc_0.2Ta_0.6O₁₂). Moreover, when a Li₇La₃Zr_0.4Hf_0.4Sn_0.4Sc_0.4Ta_0.4O₁₂ pellet is tested at 60 °C in coin cell configuration with a Li metal negative electrode, a LiNi_1/3Co_1/3Mn_1/3O₂-based positive electrode and an ionic liquid-based electrolyte at the cathode|SE interface, discharge capacity retention of about 92% is delivered after 700 cycles at 0.8 mA/cm² and 60 °C.

New evaluation of the thermodynamics stability for bcc-Fe

The thermodynamic properties for bcc-Fe were predicted by combination of the first-principles calculations, the quasiharmonic approximation, the CALPHAD method and the Weiss molecular field theory. The hybrid method considers the effects of the lattice vibration, electron, intrinsic magnetism and external magnetic fields on the thermodynamic properties at finite temperature. Combined with experimental data, the calculated heat capacity without external magnetic fields was used to verify the validity of the hybrid method. Close to the Fermi level the high electronic density of states leads to a significant electronic contribution to free energy. Near the Curie temperature lattice vibrations dominant the Gibbs free energy. The order of the other three excitation contributions to Gibbs free energy from high to low is: intrinsic magnetism > electron > external magnetic fields. The investigation suggests that all the excitation contributions to Gibbs free energy are not negligible which provides a correct direction for tuning the thermodynamic properties for Fe-based alloy.

A family of superconducting boron crystals made of stacked bilayer borophenes

Monolayer borophenes tend to be easily oxidized, while thicker borophenes have stronger antioxidation properties. Herein, we proposed four novel metallic boron crystals by stacking the experimentally synthesized borophenes, and one of the crystals has been reported in our previous experiments. Bilayer units tend to act as blocks for crystals as determined by bonding analyses. Their kinetic, thermodynamic and mechanical stabilities are confirmed by our calculated phonon spectra, molecular dynamics and elastic constants. Our proposed allotropes are more stable than the boron α-Ga phase below 1000 K at ambient pressure. Some of them become more stable than the α-rh or γ-B₂₈ phases at appropriate external pressure. More importantly, our calculations show that three of the proposed crystals are phonon-mediated superconductors with critical temperatures of about 5-10 K, higher than those of most superconducting elemental solids, in contrast to typical boron crystals with significant band gaps. Our study indicates a novel preparation method for metallic and superconducting boron crystals dispensing with high pressure.

Temperature-driven phase transition of Ti2CN from first-principles calculations

First-principles evolutionary simulations are used to predict the stable compound of Ti2CN. Body-centered tetragonal I41/amd-Ti2CN is found to be more energetically favorable than other Ti2CN compounds at 0 K. The...

First-principles study of order-disorder transitions in multicomponent solid-solution alloys

In this review, we will focus on the recent development of the order-disorder transition in metallic materials. The past decades have witnessed fast development in the first-principles methodologies and their applications to ordering transitions in multi-component alloys, particularly the high-entropy alloys. The driving force for the proceedings comes from (i) the advance of algorithms and increasingly cheaper hardware, and also (ii) the great passion to model alloys with increasing number of components. The review starts with a brief introduction of the history for the ordering transitions. More detailed scientific proceedings prior to the 1970s had been well summarized in Krivoglaz and Smirnov (1965 The Theory of Order-Disorder in Alloys (New York: Elsevier)) and Stoloff and Davies (1968 Prog. Mater. Sci. 13 1-84). In the second part, the methods to study the ordering transitions, primarily on the theoretic methods are introduced. These will include (i) KKR-CPA method and supercell methods for energetic calculations; and (ii) thermodynamic and statistical methods to compute the transition temperatures. The third part will focus on representative applications in alloys, including our own work and many others. This part supplies the primary information of this review to the readers. The fourth part will summarize the connections between ordering transitions and broader physical properties (e.g. the mechanical properties). In the last part, some concluding remarks and perspectives will be given.

Influence of defect distribution on the thermoelectric properties of FeNbSb based materials

Cooperative effects of a solid solution and phase separation could strongly scatter phonons and improve the performance of thermoelectric materials.

High-Throughput Survey of Ordering Configurations in MXene Alloys Across Compositions and Temperatures

2D transition metal carbides and nitrides known as MXenes are gaining increasing attention. About 20 of them have been synthesized (more predicted) and their applications in fields ranging from energy storage and electromagnetic shielding to medicine are being explored. To facilitate the search for double-transition-metal MXenes, we explore the structure-stability relationship for 8 MXene alloy systems, namely, (V_1-xMo_x)₃C₂, (Nb_1-xMo_x)₃C₂, (Ta_1-xMo_x)₃C₂, (Ti_1-xMo_x)₃C₂, (Ti_1-xNb_x)₃C₂, (Ti_1-xTa_x)₃C₂, (Ti_1-xV_x)₃C₂, and (Nb_1-xV_x)₃C_2, with 0 ≤ x ≤ 1, using high-throughput computations. Starting from density-functional theory calculated formation energies, we used the cluster expansion method to build quick-to-compute interactions, enabling us to scan through the formation energies of millions of alloying configurations. For the Mo-rich MXenes, (M1_1-xMo_x)₃C₂ (where M1: Ti, V, Nb, Ta) Mo atoms prefer to occupy the surface layers, and ordering persists to high temperatures, based on our Monte Carlo simulations. When Ti is alloyed with Nb or Ta, in the Ti-rich MXenes, Ti atoms prefer the surface layers (e.g., Ti-C-Nb-C-Ti sequence), and in the Nb- or Ta-rich MXenes, Ti occupies only one surface layer and the other two layers are Nb or Ta (e.g., Ti-C-Nb-C-Nb), exhibiting asymmetric ordering. However, alloying Ti with V results in solid solutions across all compositions. (Nb_1-xV_x)₃C₂ phase separates at lower temperatures but forms solid solutions at synthesis temperatures. Postsynthesis annealing at moderate temperatures (800 to 1000 K) increases the ordering for all the compositions. Lastly, by investigating the stability of their precursor MAX phases and surface-terminated MXenes, we discuss the synthesis possibilities of highly ordered MXenes.

Computationally-driven engineering of sublattice ordering in a hexagonal AlHfScTiZr high entropy alloy

Multi-principle element alloys have enormous potential, but their exploration suffers from the tremendously large range of configurations. In the last decade such alloys have been designed with a focus on random solid solutions. Here we apply an experimentally verified, combined thermodynamic and first-principles design strategy to reverse the traditional approach and to generate a new type of hcp Al-Hf-Sc-Ti-Zr high entropy alloy with a hitherto unique structure. A phase diagram analysis narrows down the large compositional space to a well-defined set of candidates. First-principles calculations demonstrate the energetic preference of an ordered superstructure over the competing disordered solid solutions. The chief ingredient is the Al concentration, which can be tuned to achieve a D019 ordering on the hexagonal lattice. The computationally designed D019 superstructure is experimentally confirmed by transmission electron microscopy and X-ray studies. Our scheme enables the exploration of a new class of high entropy alloys.

Structural γ-ε phase transition in Fe-Mn alloys from a CPA + DMFT approach

We present a computational scheme for total energy calculations of disordered alloys with strong electronic correlations. It employs the coherent potential approximation combined with the dynamical mean-field theory and allows one to study the structural transformations. The material-specific Hamiltonians in the Wannier function basis are obtained by density functional theory. The proposed computational scheme is applied to study the γ-ε structural transition in paramagnetic Fe-Mn alloys for Mn content from 10 to 20 at.%. The electronic correlations are found to play a crucial role in this transition. The calculated transition temperature decreases with increasing Mn content and is in good agreement with experiment. We demonstrate that in contrast to the α-γ transition in pure iron, the γ-ε transition in Fe-Mn alloys is driven by a combination of kinetic and Coulomb energies. The latter is found to be responsible for the decrease of the γ-ε transition temperature with Mn content.

Element-resolved thermodynamics of magnetocaloric LaFe(13-x)Si(x)

By combination of two independent approaches, nuclear resonant inelastic x-ray scattering and first-principles calculations in the framework of density functional theory, we demonstrate significant changes in the element-resolved vibrational density of states across the first-order transition from the ferromagnetic low temperature to the paramagnetic high temperature phase of LaFe(13-x)Si(x). These changes originate from the itinerant electron metamagnetism associated with Fe and lead to a pronounced magneto-elastic softening despite the large volume decrease at the transition. The increase in lattice entropy associated with the Fe subsystem is significant and contributes cooperatively with the magnetic and electronic entropy changes to the excellent magneto- and barocaloric properties.

Interplay between subsurface ordering, surface segregation, and adsorption on Pt-Ti(111) near-surface alloys

Using the first-principles cluster expansion (CE) method, we studied the subsurface ordering of Pt/Pt-Ti(111) surface alloys and the effect of this ordering on segregation and adsorption behavior. The clusters included in the CE are optimized by a genetic algorithm to better describe the interactions between Pt and Ti atoms in the subsurface layer. Similar to bulk Pt-Ti alloys, Pt-Ti(111) subsurface alloys show a strong ordering tendency. A series of stable ordered Pt-Ti subsurface structures are identified from the two-dimensional (2D) CE. As an indication of the connection between the 2D and the bulk ordering, the CE predicts a ground-state Pt(8)Ti structure in the (111) subsurface layer, which is the same ordering as the close-packed plane of the bulk Pt(8)Ti compound. We carried out Monte Carlo simulations (MC) using the CE Hamiltonian to study the finite temperature stability of the Pt-Ti subsurface structures. The MC results show that subsurface structures in the Pt-rich range have higher order-disorder transition temperatures than their Ti-rich subsurface counterparts. We calculate the binding energy of different adsorbates (O, S, H, and NO) on Pt-terminated and Ti-segregated surfaces of ordered PtTi and Pt(8)Ti subsurface alloys. The binding of these adsorbates is generally stronger on Ti-segregated surfaces than Pt-terminated surfaces. The adsorption-induced Ti surface segregation is determined by two factors: (i) the unfavorable energy penalty for the Ti atom to segregate to the clean surface and (ii) the favorable energy decrease from stronger adsorbate binding on the Ti-segregated surface. The two factors introduce similar magnitude in energy change for the S and NO adsorption on Ti-segregated surfaces of PtTi subsurface alloys. We predict an adsorption-induced Ti surface segregation that is dependent on the atomic configurations of the Ti-segregated surfaces resulting from the competition of the two factors.

Configurational electronic entropy and the phase diagram of mixed-valence oxides: the case of LixFePO4

We demonstrate that configurational electronic entropy, previously neglected, in ab initio thermodynamics of materials can qualitatively modify the finite-temperature phase stability of mixed-valence oxides. While transformations from low-T ordered or immiscible states are almost always driven by configurational disorder (i.e., random occupation of lattice sites by multiple species), in FePO4-LiFePO4 the formation of a solid solution is almost entirely driven by electronic rather than ionic configurational entropy. We argue that such an electronic entropic mechanism may be relevant to most other mixed-valence systems.

Cluster expansion of electronic excitations: Application to fcc Ni–Al alloys

The cluster expansion method is applied to electronic excitations and a set of effective cluster densities of states (ECDOS) is defined, analogous to effective cluster interactions (ECIs). The ECDOSs are used to generate alloy thermodynamic properties as well as the equation of state (EOS) of electronic excitations for the fcc Ni-Al systems. When parent clusters have a small size, the convergence of the expansion is not so good but the electronic density of state (DOS) is well reproduced. However, the integrals of the DOS such as the cluster expanded free energy, entropy, and internal energy associated with electronic excitations are well described at the level of the tetrahedron-octahedron cluster approximation, indicating that the ECDOS is applicable to produce electronic ECIs for cluster variation method (CVM) or Monte Carlo calculations. On the other hand, the Gruneisen parameter, calculated with first-principles methods, is no longer a constant and implies that the whole DOS profile should be considered for EOS of electronic excitations, where ECDOS adapts very well for disordered alloys and solid solutions.

Reliable first-principles alloy thermodynamics via truncated cluster expansions

In alloys cluster expansions (CE) are increasingly used to combine first-principles electronic-structure calculations and Monte Carlo methods to predict thermodynamic properties. As a basis-set expansion in terms of lattice geometrical clusters and effective cluster interactions, the CE is exact if infinite, but is tractable only if truncated. Yet until now a truncation procedure was not well defined and did not guarantee a reliable truncated CE. We present an optimal truncation procedure for CE basis sets that provides reliable thermodynamics. We then exemplify its importance in Ni3V, where the CE has failed unpredictably, and now show agreement to a range of measured values, predict new low-energy structures, and explain the cause of previous failures.