Universal Maximum Strength of Solid Metals and Alloys.

Interstitial electron density ρ_{o} is offered as a direct metric for maximum strength in metals, arising from universal properties derived from an electron gas. ρ_{o} sets the exchange-correlation parameter r_{s} in density-functional theory. It holds also for maximum shear strength τ_{max} in polycrystals [M. Chandross and N. Argibay, Phys. Rev. Lett. 124, 125501 (2020)PRLTAO0031-900710.1103/PhysRevLett.124.125501]. Elastic moduli and τ_{max} for polycrystalline (amorphous) metals are linear with ρ_{o} and melting T_{m} (glass-transition T_{g}) temperature. ρ_{o} or r_{s}, even with rule-of-mixture estimate, predicts relative strength for rapid, reliable selection of high-strength alloys with ductility, as confirmed for elements to steels to complex solid solutions, and validated experimentally.

Martensitic Transformation in Fe_{x}Mn_{80-x}Co_{10}Cr_{10} High-Entropy Alloy.

High-entropy alloys and even medium-entropy alloys are an intriguing class of materials in that structure and property relations can be controlled via alloying and chemical disorder over wide ranges in the composition space. Employing density-functional theory combined with the coherent-potential approximation to average over all chemical configurations, we tune free energies between face-centered-cubic and hexagonal-close-packed phases in Fe_{x}Mn_{80-x}Co_{10}Cr_{10} systems. Within Fe-Mn-based alloys, we show that the martensitic transformation and chemical short-range order directly correlate with the face-centered-cubic and hexagonal-close-packed energy difference and stacking-fault energies, which are in quantitative agreement with recent observation of two phase region (face-centered cubic and hexagonal closed pack) in a polycrystalline high-entropy alloy sample at x=40 at.%. Our predictions are further confirmed by single-crystal measurements on a x=40 at.% using transmission-electron microscopy, selective-area diffraction, and electron-backscattered-diffraction mapping. The results herein offer an understanding of transformation-induced or twinning-induced plasticity in this class of high-entropy alloys and a design guide for controlling the physics at the electronic level.

Revealing the Nature of Antiferroquadrupolar Ordering in Cerium Hexaboride: CeB_{6}.

The cerium hexaboride (CeB_{6}) f-electron compound displays a rich array of low-temperature magnetic phenomena, including a "magnetically hidden" order, identified as multipolar in origin via advanced x-ray scattering. From first-principles electronic-structure results, we find that the antiferroquadrupolar (AFQ) ordering in CeB_{6} arises from crystal-field splitting and yields a band structure in agreement with experiments. With interactions of p electrons between Ce and B_{6} being small, the electronic state of CeB_{6} is suitably described as Ce(4f^{1})^{3+}(e^{-})(B_{6})^{2-}. The AFQ state of orbital spins is caused by an exchange interaction induced through spin-orbit interaction, which also splits the J=5/2 state into a Γ_{8} ground state and a Γ_{7} excited state. Within the smallest antiferromagnetic (AFM) (111) configuration, an orbital-ordered AFQ state appears during charge self-consistency, and it supports the appearance of a "hidden" order. Hydrostatic pressure (either applied or chemically induced) stabilizes the AFM (AFQ) states over a ferromagnetic one, as observed at low temperatures.

Lattice Instability during Solid-Solid Structural Transformations under a General Applied Stress Tensor: Example of Si I→Si II with Metallization.

The density functional theory was employed to study the stress-strain behavior and elastic instabilities during the solid-solid phase transformation (PT) when subjected to a general stress tensor, as exemplified for semiconducting Si I and metallic Si II, where metallization precedes the PT, so stressed Si I can be a metal. The hydrostatic PT occurs at 76 GPa, while under uniaxial loading it is 11 GPa (3.7 GPa mean pressure), 21 times lower. The Si I→Si II PT is described by a critical value of the phase-field's modified transformation work, and the PT criterion has only two parameters given six independent stress elements. Our findings reveal novel, more practical synthesis routes for new or known high-pressure phases under predictable nonhydrostatic loading, where competition of instabilities can serve for phase selection rather than free energy minima used for equilibrium processing.

Characterizing Substrate-Surface Interactions on Alumina-Supported Metal Catalysts by Dynamic Nuclear Polarization-Enhanced Double-Resonance NMR Spectroscopy.

The characterization of nanometer-scale interactions between carbon-containing substrates and alumina surfaces is of paramount importance to industrial and academic catalysis applications, but it is also very challenging. Here, we demonstrate that dynamic nuclear polarization surface-enhanced NMR spectroscopy (DNP SENS) allows the unambiguous description of the coordination geometries and conformations of the substrates at the alumina surface through high-resolution measurements of 13C-27Al distances. We apply this new technique to elucidate the molecular-level geometry of 13C-enriched methionine and natural abundance poly(vinyl alcohol) adsorbed on γ-Al2O3-supported Pd catalysts, and we support these results with element-specific X-ray absorption near-edge measurements. This work clearly demonstrates a surprising bimodal coordination of methionine at the Pd-Al2O3 interface.

Nanoalloy electrocatalysis: simulating cyclic voltammetry from configurational thermodynamics with adsorbates.

We simulate the adsorption isotherms for alloyed nanoparticles (nanoalloys) with adsorbates to determine cyclic voltammetry (CV) during electrocatalysis. The effect of alloying on nanoparticle adsorption isotherms is provided by a hybrid-ensemble Monte Carlo simulation that uses the cluster expansion method extended to non-exchangeable coupled lattices for nanoalloys with adsorbates. Exemplified here for the hydrogen evolution reaction, a 2-dimensional CV is mapped for Pd-Pt nanoalloys as a function of both electrochemical potential and the global Pt composition, and shows a highly non-linear alloying effect on CV. Detailed features in CV arise from the interplay among the H-adsorption in multiple sites that is closely correlated with alloy configurations, which are in turn affected by the H-coverage. The origins of specific features in CV curves are assigned. The method provides a more complete means to design nanoalloys for electrocatalysis.

Nudged-elastic band method with two climbing images: finding transition states in complex energy landscapes.

The nudged-elastic band (NEB) method is modified with concomitant two climbing images (C2-NEB) to find a transition state (TS) in complex energy landscapes, such as those with a serpentine minimal energy path (MEP). If a single climbing image (C1-NEB) successfully finds the TS, then C2-NEB finds it too. However, improved stability of C2-NEB makes it suitable for more complex cases, where C1-NEB misses the TS because the MEP and NEB directions near the saddle point are different. Generally, C2-NEB not only finds the TS, but also guarantees, by construction, that the climbing images approach it from the opposite sides along the MEP. In addition, C2-NEB provides an accuracy estimate from the three images: the highest-energy one and its climbing neighbors. C2-NEB is suitable for fixed-cell NEB and the generalized solid-state NEB.

Configurational thermodynamics of alloyed nanoparticles with adsorbates.

Changes in the chemical configuration of alloyed nanoparticle (NP) catalysts induced by adsorbates under working conditions, such as reversal in core-shell preference, are crucial to understand and design NP functionality. We extend the cluster expansion method to predict the configurational thermodynamics of alloyed NPs with adsorbates based on density functional theory data. Exemplified with PdRh NPs having O-coverage up to a monolayer, we fully detail the core-shell behavior across the entire range of NP composition and O-coverage with quantitative agreement to in situ experimental data. Optimally fitted cluster interactions in the heterogeneous system are the key to enable quantitative Monte Carlo simulations and design.

Lifshitz transition and chemical instabilities in Ba(1-x)K(x)Fe2As2 superconductors.

For solid-solution Ba1-xKxFe2As2 Fermi surface evolution is mapped via Bloch spectral functions calculated using density functional theory implemented in Korringa-Kohn-Rostoker multiple scattering theory with the coherent-potential approximation. Spectral functions reveal electronic dispersion, topology, orbital character, and broadening (electron-lifetime effects) due to chemical disorder. Dissolution of electron cylinders occurs near x∼0.9 with a nonuniform, topological (Lifshitz) transition, reducing the interband interactions; yet the dispersion maintains its dxz or dyz character. Formation energies indicate alloying at x=0.35, as observed, and a tendency for segregation on the K-rich (x>0.6) side, explaining the difficulty of controlling sample quality and the conflicting results between characterized electronic structures. Our results reveal Fermi surface transitions in alloyed samples that influence s± to nodal superconductivity and suggest the origin for deviations of common trends in Fe-based superconductors, such as Bud'ko-Ni-Canfield scaling.

Selective host-guest interaction between metal ions and metal-organic frameworks using dynamic nuclear polarization enhanced solid-state NMR spectroscopy.

The host-guest interaction between metal ions (Pt(2+) and Cu(2+) ) and a zirconium metal-organic framework (UiO-66-NH2 ) was explored using dynamic nuclear polarization-enhanced (15) N{(1) H} CPMAS NMR spectroscopy supported by X-ray absorption spectroscopy and density functional calculations. The combined experimental results conclude that each Pt(2+) coordinates with two NH2 groups from the MOF and two Cl(-) from the metal precursor, whereas Cu(2+) do not form chemical bonds with the NH2 groups of the MOF framework. Density functional calculations reveal that Pt(2+) prefers a square-planar structure with the four ligands and resides in the octahedral cage of the MOF in either cis or trans configurations.

Antiferroelectric Ceramics for Energy-Efficient Capacitors by Theory-Guided Discovery.

Antiferroelectric ceramics, via the electric-field-induced antiferroelectric (AFE)-ferroelectric (FE) phase transitions, show great promise for high-energy-density capacitors. Yet, currently, only 70-80% energy release is found during a charge-discharge cycle. Here, for PbZrO3-based oxides, geometric nonlinear theory of martensitic phase transitions is applied (first used to guide supercompatible shape-memory alloys) to predict the reversibility of the AFE-FE transition by using density-functional theory to assess AFE/FE interfacial lattice-mismatch strain that assures ultralow electric hysteresis and extended fatigue lifetime. A good correlation of mismatch strain with electric hysteresis, hence, with energy efficiency of AFE capacitors is observed. Guided by theory, high-throughput material search is conducted and AFE compositions with a near-perfect charge-discharge energy efficiency (98.2%), i.e., near-zero hysteresis are discovered. And the fatigue life of the capacitor reaches 79.5 million charge-discharge cycles, a factor of 80 enhancement over AFE ceramics with large electric hysteresis.

Intermetallic NaAu2 as a heterogeneous catalyst for low-temperature CO oxidation.

The enhanced stability and modified electronic structure of intermetallic compounds provide discovery of superior catalysts for chemical conversions with high activity, selectivity, and stability. We find that the intermetallic NaAu2 is an active catalyst for CO oxidation at low temperatures. From density functional theory calculations, a reaction mechanism is suggested to explain the observed low reaction barrier of CO oxidation by NaAu2, in which a CO molecule reacts directly with an adsorbed O2 to form an OOCO* intermediate. The presence of surface Na increases the binding energy of O2 and decreases the energy barrier of the transition states.

Noncrystalline-to-crystalline transformations in Pt nanoparticles.

We show that the noncrystalline-to-crystalline transition of supported Pt nanoparticles (NPs) in the subnanometer to nanometer size range is statistical in nature, and strongly affected by particle size, support, and adsorbates (here we use H2). Unlike in the bulk, a noncrystalline phase exists and is stable in small NPs, reflecting a general mesoscopic feature. Observations of >3000 particles by high-resolution transmission electron microscopy show a noncrystalline-to-crystalline transition zone that is nonabrupt; there is a size regime where disordered and ordered NPs coexist. The NP size at which this transition occurs is strongly dependent on both the adsorbate and the support, and this effect is general for late 5d transition metals. All results are reconciled via a statistical description of particle-support-adsorbate interactions.

A comprehensive search for stable Pt-Pd nanoalloy configurations and their use as tunable catalysts.

Using density-functional theory, we predict stable alloy configurations (ground states) for a 1 nm Pt-Pd cuboctahedral nanoparticle across the entire composition range and demonstrate their use as tunable alloy catalysts via hydrogen-adsorption studies. Unlike previous works, we use simulated annealing with a cluster expansion Hamiltonian to perform a rapid and comprehensive search that encompasses both high and low-symmetry configurations. The ground states show Pt(core)-Pd(shell) type configurations across all compositions but with specific Pd patterns. For catalysis studies at room temperatures, the ground states are more realistic structural models than the commonly assumed random alloy configurations. Using the ground states, we reveal that the hydrogen adsorption energy increases (decreases) monotonically with at. % Pt for the {111} hollow ({100} bridge) adsorption site. Such trends are useful for designing tunable Pd-Pt nanocatalysts for the hydrogen evolution reaction.

A generalized solid-state nudged elastic band method.

A generalized solid-state nudged elastic band (G-SSNEB) method is presented for determining reaction pathways of solid-solid transformations involving both atomic and unit-cell degrees of freedom. We combine atomic and cell degrees of freedom into a unified description of the crystal structure so that calculated reaction paths are insensitive to the choice of periodic cell. For the rock-salt to wurtzite transition in CdSe, we demonstrate that the method is robust for mechanisms dominated either by atomic motion or by unit-cell deformation; notably, the lowest-energy transition mechanism found by our G-SSNEB changes with cell size from a concerted transformation of the cell coordinates in small cells to a nucleation event in large cells. The method is efficient and can be applied to systems in which the force and stress tensor are calculated using density functional theory.

Effect of Pd alloying on structural, electronic and magnetic properties of L10Fe-Ni.

We present a systematic study of the effect of Pd-alloying on phase stability, electronic structure, and elastic properties in L10Fe-Ni using density-functional theory. Being from the same group of the periodic table, Pd is the best candidate for chemical alloying. The Fe-Ni/Fe-Pd/Ni-Pd bond-length increases with increasing Pd-concentration, which weakens the hybridization between low lying energy states below Fermi-level. The reduced hybridization decreases the relative thermodynamic stability of L10Fe(Ni1-xPdx) untilx= 0.75. Beyond this concentration, the relative stability gets enhanced, which is attributed to a unique change in the lattice distortion (c/a). The elastic properties show a non-monotonous behavior as a function ofx, which is again due to a specific change-over in the uniaxial strain. We found that Pd alloying increases the local Fe moment and structural anisotropy of L10FeNi, which are important for applications such as microwave absorption, refrigeration systems, recording devices, imaging and sensors. We believe that the present study for the chemical alloying effect can provide critical insights toward the understanding of electronic-structure and elastic behavior of other technologically important materials.

Accelerating computational modeling and design of high-entropy alloys.

High-entropy alloys, with N elements and compositions {cν = 1,N} in competing crystal structures, have large design spaces for unique chemical and mechanical properties. Here, to enable computational design, we use a metaheuristic hybrid Cuckoo search (CS) to construct alloy configurational models on the fly that have targeted atomic site and pair probabilities on arbitrary crystal lattices, given by supercell random approximates (SCRAPs) with S sites. Our Hybrid CS permits efficient global solutions for large, discrete combinatorial optimization that scale linearly in a number of parallel processors, and linearly in sites S for SCRAPs. For example, a four-element, 128-site SCRAP is found in seconds-a more than 13,000-fold reduction over current strategies. Our method thus enables computational alloy design that is currently impractical. We qualify the models and showcase application to real alloys with targeted atomic short-range order. Being problem-agnostic, our Hybrid CS offers potential applications in diverse fields.

Machine learning assisted prediction of the Young's modulus of compositionally complex alloys.

We identify compositionally complex alloys (CCAs) that offer exceptional mechanical properties for elevated temperature applications by employing machine learning (ML) in conjunction with rapid synthesis and testing of alloys for validation to accelerate alloy design. The advantages of this approach are scalability, rapidity, and reasonably accurate predictions. ML tools were implemented to predict Young's modulus of refractory-based CCAs by employing different ML models. Our results, in conjunction with experimental validation, suggest that average valence electron concentration, the difference in atomic radius, a geometrical parameter λ and melting temperature of the alloys are the key features that determine the Young's modulus of CCAs and refractory-based CCAs. The Gradient Boosting model provided the best predictive capabilities (mean absolute error of 6.15 GPa) among the models studied. Our approach integrates high-quality validation data from experiments, literature data for training machine-learning models, and feature selection based on physical insights. It opens a new avenue to optimize the desired materials property for different engineering applications.

Incommensurate transition-metal dichalcogenides via mechanochemical reshuffling of binary precursors.

A new family of heterostructured transition-metal dichalcogenides (TMDCs) with incommensurate ("misfit") spatial arrangements of well-defined layers was prepared from structurally dissimilar single-phase 2H-MoS2 and 1T-HfS2 materials. The experimentally observed heterostructuring is energetically favorable over the formation of homogeneous multi-principle element dichalcogenides observed in related dichalcogenide systems of Mo, W, and Ta. The resulting three-dimensional (3D) heterostructures show semiconducting behavior with an indirect band gap around 1 eV, agreeing with values predicted from density functional theory. Results of this joint experimental and theoretical study open new avenues for generating unexplored metal-dichalcogenide heteroassemblies with incommensurate structures and tunable physical properties.

CoCrFeNi High-Entropy Alloy as an Enhanced Hydrogen Evolution Catalyst in an Acidic Solution.

High-entropy alloys (HEAs) have intriguing material properties, but their potential as catalysts has not been widely explored. Based on a concise theoretical model, we predict that the surface of a quaternary HEA of base metals, CoCrFeNi, should go from being nearly fully oxidized except for pure Ni sites when exposed to O2 to being partially oxidized in an acidic solution under cathodic bias, and that such a partially oxidized surface should be more active for the electrochemical hydrogen evolution reaction (HER) in acidic solutions than all the component metals. These predictions are confirmed by electrochemical and surface science experiments: the Ni in the HEA is found to be most resistant to oxidation, and when deployed in 0.5 M H2SO4, the HEA exhibits an overpotential of only 60 mV relative to Pt for the HER at a current density of 1 mA/cm2.