Engineering topics
Wolverton, C.
Publications and source records attributed to Wolverton, C..
CeTi 2 O 6 —A Promising Oxide for Solar Thermochemical Hydrogen Production
A large entropy of reduction is crucial in achieving materials capable of high-efficiency solar thermochemical hydrogen (STCH) production through two-step thermochemical water splitting cycles. We have recently demonstrated that the onsite electronic entropy of reduction attains an extreme value of 4.26 kB at 1500 K in Ce 4+ → Ce 3+ redox reactions, which explains the high performance and uniqueness of CeO 2 as an archetypal STCH material. However, ceria requires high temperatures (T > 1500 °C) to achieve a reasonable reduction extent because of its large reduction enthalpy, which is a major obstacle in practical applications. Therefore, new materials with a large entropy of reduction and lower reduction enthalpy are required. Here, we perform a systematic screening to search for Ce 4+ -based oxides which possess thermodynamics superior to CeO 2 for STCH production. We first search the Inorganic Crystal Structure Database (ICSD) and literature for Ce 4+ -based oxides and subsequently use density functional theory to compute their reduction enthalpies (i.e., oxygen vacancy formation energies). Here, we find that CeTi 2 O 6 with the brannerite structure is the most promising candidate for STCH because it possesses three essential characteristics of an STCH material: (i) a smaller reduction enthalpy compared to ceria yet large enough to split water, (ii) a high thermal stability, as reported experimentally, and (iii) a large entropy of reduction associated with Ce 4+ → Ce 3+ redox. Our proposed design strategy suggests that further exploration of Ce 4+ oxides for STCH production is warranted.
High thermoelectric performance in BaAgYTe 3 via low lattice thermal conductivity induced by bonding heterogeneity
Solid-state technology based on thermoelectric (TE) materials enables the conversion of heat into electricity offering an environmentally friendly solution to energy conservation. Here, we use density functional theory calculations to show that BaAgYTe 3 , a layered semiconductor, exhibits low lattice thermal conductivity (κ l ) and a high thermoelectric figure of merit. Our calculations reveal that the presence of bonding inhomogeneity, resulting from the coexisting rigid and fluctuating sublattices, favorably helps in both the electronic and phonon transports. We show that low κ l in this compound mainly originates from (a) the small group velocities of the acoustic modes, (b) quasilocalized low-frequency optical phonons that give rise to multiple scattering channels, and (c) strong lattice anharmonicity. While the calculations of the atomic displacement parameters and bonding analysis reveal relatively weaker bonding of Ag atoms and establish the heterogeneity in the chemical bonding, the strong anharmonicity is manifested in the large mode Gruneisen parameters. Thus, our work provides a theoretical prediction that warrants experimental verification and should encourage further exploration of potential TE materials in the same crystal family.
Bonding Hierarchy Gives Rise to High Thermoelectric Performance in Layered Zintl Compound BaAu 2 P 4
The search for new thermoelectric materials has gained rapid progress in recent years as thermoelectric technology offers the potential for environmentally friendly and sustainable energy conversion methods from waste heat to electricity. In this work, we use first-principles calculations based on density functional theory to predict high thermoelectric performance in BaAu 2 P 4 , a layered Zintl compound with a small band gap. BaAu 2 P 4 exhibits crystallographic heterogeneity in which rigid [Au 2 P 4 ] 2– units are separated by layers of Ba 2+ cations, which are bonded relatively weakly to the lattice through electrostatic interactions. The phosphorus atoms are covalently bonded to each other and form infinite chains within the crystal. While the phosphorus chains facilitate large electrical conductivity, the presence of multiple bands near the Fermi level gives rise to an enhanced Seebeck coefficient. On the other hand, the loosely bound Ba along with Au strongly scatter the heat carrying acoustic phonons, significantly reducing the lattice thermal conduction along the stacking direction. As a consequence of this bonding hierarchy (i.e., coexisting rigid and fluctuating sublattices), BaAu 2 P 4 exhibits a large power factor and low lattice thermal conductivity, which results in a high thermoelectric figure of merit ($zT$). Thus, our findings should encourage the exploration of new thermoelectric materials in the family of layered compounds with small band gaps and crystallographic heterogeneity.
Relationship between pair and higher-order correlations in solid solutions and other Ising systems
Abstract not provided
Computational investigation of inverse Heusler compounds for spintronics applications
In this work, first-principles calculations of the electronic structure, magnetism, and structural stability of inverse Heusler compounds with the chemical formula X 2 YZ are presented and discussed with a goal of identifying compounds of interest for spintronics. Compounds for which the number of electrons per atom for Y exceed that for X and for which X and Y are each one of the 3d elements, Sc-Zn, and Z is one of the group IIIA-VA elements: Al, Ga, In, Si, Ge, Sn, P, As, or Sb were considered. The formation energy per atom of each compound was calculated. By comparing our calculated formation energies to those calculated for phases in the inorganic crystal structure database of observed phases, we estimate that inverse Heuslers with formation energies within 0.052 eV/atom of the calculated convex hull are reasonably likely to be synthesizable in equilibrium. The observed trends in the formation energy and relative structural stability as the X, Y, and Z elements vary are described. In addition to the Slater-Pauling gap after 12 states per formula unit in one of the spin channels, inverse Heusler phases often have gaps after 9 states or 14 states. We describe the origin and occurrence of these gaps. We identify 14 inverse Heusler semiconductors, 51 half-metals, and 50 near-half-metals with negative formation energy. Additonally, our calculations predict 4 half-metals and 6 near-half-metals to lie close to the respective convex hull of stable phases, and thus may be experimentally realized under suitable synthesis conditions, resulting in potential candidates for future spintronics applications.
Energetics of native defects, solute partitioning, and interfacial energy of Q precipitate in Al-Cu-Mg-Si alloys
The compound Al 3 Cu 2 Mg 9 Si 7 , is known as the Q-phase and forms as a thermodynamically stable precipitate during aging in the quaternary Al-Cu-Mg-Si system. We perform atomic-scale density functional theory (DFT) calculations of defect properties, solute partitioning, and interfacial stability of the Al 3 Cu 2 Mg 9 Si 7 (Q) precipitate. We find: (i) simple native point defect (i.e., vacancies and anti-sites) thermodynamics can partially explain the experimentally observed off-stoichiometry, such as the observed variation of compositions, Al 3+δ Cu 2 Mg 9-δ Si 7 (Mg-deficient and Al-rich) in experiment. (ii) Calculated solute-partitioning energies of common solutes allow us to define general rules for site-preference in the Q-phase in terms of electronic structure and atomic radius. To validate our DFT predictions, we perform atom-probe tomography (APT) experiments for six-different elements (Zn, Ni, Mn, Ti, V, and Zr). The results show that the partitioning behavior of solutes Ni, Zn, and Mn are consistent with DFT predictions, but the transition elements (Ti, V, and Zr), which are anomalously slow diffusers in Al, partition to the Q-phase in constrast to DFT partitioning energies. (iii) For the low energy interface ($11\bar{2}$) Q //(510) Al observed in needle shaped Q-precipitate, we survey various terminations and orientations and derive a low-energy interfacial structure. Here we find this low-energy interfacial model has Cu atoms nearest to the interface, which is in agreement with previous literature on Cu interfacial segregation at the Q'//α-Al interface. The computed interfacial energy (0.52 J/m 2 ) and the corresponding structure will be useful input to future multi-scale modeling of microstructural evolution.
First-principles study of crystal structure and stability of T 1 precipitates in Al-Li-Cu alloys
Aluminum-lithium-copper alloys have a low density, high elastic modulus and high specific strength.Due to this combination of properties, alloys strengthened with the ternary (Al-Li-Cu) T 1 phase have attracted a great deal of interest especially in aerospace applications. Determining the atomic structuralinformation of the precipitate is a fundamental step in developing a basis for advanced alloy design; however, even though many experimental studies have addressed the T 1 crystal structure, it remains the subject of some controversy. Here, we use density functional theory (DFT) calculations to investigate the structure and composition of the T 1 phase by comparing the energetic stability of five previously-proposed models of the crystal structure of T 1 . The DFT formation energy of these proposed T 1 crystal structures was calculated using a special quasi-random structure (SQS) approach to describe a disordered Al-Cu sub-lattice. In conflict with the experimental phase diagram, DFT calculations of all five proposed models result in an energetically unstable T 1 phase. Here, we search for a new, lower-energy structure of T 1 using a cluster expansion approach, and find a new structural model with DFT energy that is stable (at T = 0 K), i.e., on the calculated convex hull of the Al-Li-Cu ternary system. However, this new predicted phase does not have a tie-line with Al, but the formation energy of the phase is very close to the energy required to make a tie-line with Al (Δ E = 0.013 eV/atom), which could be affected by finite temperature entropic effects (i.e., vibrational entropic stabilization).
Electrochemically Stable Coating Materials for Li, Na, and Mg Metal Anodes in Durable High Energy Batteries
Anodes made of Li, Na, or Mg metal present a rare opportunity to double the energy density of rechargeable batteries. However, these metals are highly reactive with many electrolytes and yield electronically conductive phases that allow continued electrochemical reduction of the electrolyte. This reactivity degrades cell performance over time and poses a safety risk. Surface coatings on metal anodes can limit reactivity with electrolytes and improve durability. In this paper, we screen the Open Quantum Materials Database (OQMD) to identify coatings that exhibit chemical equilibrium with the anode metals and are electronic insulators. We rank the coatings according to their electronic bandgap. We identify 92 coatings for Li anodes, 118 for Na anodes, and 97 for Mg anodes. Only two compounds that are commonly studied as Li solid electrolytes pass our screens: Li 3 N and Li 7 La 3 Hf 2 O 12 . Many of the coatings that we identify are new to the battery literature. We suggest further study of these coatings to validate their performance in cells.
First-principles/Phase-field modeling of θ′ precipitation in Al-Cu alloys
Not Available
Predicting the morphologies of γ' precipitates in cobalt-based superalloys
Cobalt-based alloys with γ/γ' microstructures have the potential to become the next generation of superalloys, but alloy compositions and processing steps must be optimized to improve coarsening, creep, and rafting behavior. While these behaviors are different than in nickel-based superalloys, alloy development can be accelerated by understanding the thermodynamic factors influencing microstructure evolution. In this work, we develop a phase field model informed by first-principles density functional theory and experimental data to predict the equilibrium shapes of Co-Al-W γ' precipitates. Three-dimensional simulations of single and multiple precipitates are performed to understand the effect of elastic and interfacial energy on coarsened and rafted microstructures; the elastic energy is dependent on the elastic stiffnesses, misfit strain, precipitate size, applied stress, and precipitate spatial distribution. We observe characteristic microstructures dependent on the type of applied stress that have the same γ' morphology and orientation seen in experiments, indicating that the elastic stresses arising from coherent γ/γ' interfaces are important for morphological evolution during creep. Here, the results also indicate that the narrow γ channels between γ' precipitates are energetically favored, and provide an explanation for the experimentally observed directional coarsening that occurs without any applied stress.
Phosphorus allotropes: Stability of black versus red phosphorus re-examined by means of the van der Waals inclusive density functional method
We evaluate the energetic stabilities of white, red, and black allotropes of phosphorus using density functional theory (DFT) and hybrid functional methods, van der Waals (vdW) corrections (DFT+vdW and hybrid+vdW), vdW density functionals, and random phase approximation (RPA). In this work, we find that stability of black phosphorus over red-V (i.e., the violet form) is not ubiquitous among these methods, and the calculated enthalpies for the reaction phosphorus (red-V)→ phosphorus (black) are scattered between –20 and 40 meV/atom. With local density and generalized gradient approximations, and hybrid functionals, mean absolute errors (MAEs) in densities of P allotropes relative to experiments are found to be around 10%–25%, whereas with vdW-inclusive methods, MAEs in densities drop below ~5%. While the inconsistency among the density functional methods could not shed light on the stability puzzle of black versus red phosphorus, comparison of their accuracy in predicting densities and the supplementary RPA results on relative stabilities indicate that opposite to the common belief, black and red phosphorus are almost degenerate, or the red-V (violet) form of phosphorus might even be the ground state.
Diffusion coefficients of transition metals in fcc cobalt
Using first-principles density functional theory (DFT), we calculate the diffusivities of 32 different solute elements—all transition metals, together with Al and Si—in fcc cobalt within the formalism of the five-frequency model. For self-diffusion in fcc cobalt, we compare the accuracy of various approximations to the exchange-correlation energy functional of DFT in estimating the activation energy, and find that only the Perdew-Burke-Ernzerhof (PBE) approximation agrees well with experimental reports and all other functionals largely overestimate it. Our calculations also show that an accurate estimation of the self-diffusion coefficient requires explicit calculation of the effective jump frequency and vacancy formation entropy via phonons. Using accurate self-diffusion data and scaling all solute-related attempt frequencies with respect to the attempt frequency for self-diffusion using a simple relation involving the atomic mass and melting temperature of the solute yields solute diffusivities in excellent agreement with experiments, where such data is available. Additionally, we find that large solutes spontaneously relax toward the nearest neighbor vacancy to relieve the misfit strain, and the extent of this relaxation correlates negatively with the migration energy. Thus, in general, larger solutes have lower migration energies and diffuse faster than smaller solutes in fcc cobalt. However, extremely large solutes, e.g., group III elements Sc, Y, Lu, tend to be trapped in an energy valley located halfway toward the vacancy, and monovacancy mediated diffusion may no longer be valid in such cases. Finally, for all the solutes considered, we systematically tabulate the diffusion-related quantities calculated—diffusion prefactors, migration and activation energies—constructing an extensive and accurate first-principles database for solute diffusion in fcc cobalt.
Computational investigation of half-Heusler compounds for spintronics applications
We present first-principles density functional calculations of the electronic structure, magnetism, and structural stability of 378 XYZ half-Heusler compounds (with X= Cr, Mn, Fe, Co, Ni, Ru, Rh; Y=Ti, V, Cr, Mn, Fe, Ni; Z= Al, Ga, In, Si, Ge, Sn, P, As, Sb). We find that a “Slater-Pauling gap” in the density of states (i.e., a gap or pseudogap after nine states in the three atom primitive cell) in at least one spin channel is a common feature in half-Heusler compounds. We find that the presence of such a gap at the Fermi energy in one or both spin channels contributes significantly to the stability of a half-Heusler compound. We calculate the formation energy of each compound and systematically investigate its stability against all other phases in the open quantum materials database (OQMD). Here, we represent the thermodynamic phase stability of each compound as its distance from the convex hull of stable phases in the respective chemical space and show that the hull distance of a compound is a good measure of the likelihood of its experimental synthesis. We find low formation energies and mostly correspondingly low hull distances for compounds with X= Co, Rh, or N, Y= Ti or V, and Z= P, As, Sb, or Si. We identify 26 18-electron semiconductors, 45 half-metals, and 34 near half-metals with negative formation energy that follow the Slater-Pauling rule of three electrons per atom. Our calculations predict several new, as-yet unknown, thermodynamically stable phases, which merit further experimental exploration—RuVAs, CoVGe, FeVAs in the half-Heusler structure, and NiScAs, RuVP, RhTiP in the orthorhombic MgSrSi-type structure. Further, two interesting zero-moment half-metals, CrMnAs and MnCrAs, are calculated to have negative formation energy. In addition, our calculations predict a number of hitherto unreported semiconducting (e.g., CoVSn and RhVGe), half-metallic (e.g., RhVSb), and near half-metallic (e.g., CoFeSb and CoVP) half-Heusler compounds to lie close to the respective convex hull of stable phases, and thus may be experimentally realized under suitable synthesis conditions, resulting in potential candidates for various semiconducting and spintronics applications.
Oxygen evolution from olivine Mn 1– x M x PO 4 ( M =Fe,Ni,Al,Mg) delithiated cathode materials
Olivine LiMnPO 4 is a promising cathode material for Li-ion batteries. One drawback of this material is the propensity of its delithiated phase, MnPO 4 , to evolve oxygen gas above approximately 200 °C. During thermal runaway of cells, this oxygen gas can burn the electrolyte and other cell components and thereby jeopardize safety. Partial substitution of Mn with M = Fe, Ni, Al, or Mg has been used to improve the lithium intercalation kinetics of Li x MnPO 4 ; however, the effect of these substitutions on oxygen evolution is not fully documented. In this paper, we calculate phase diagrams and oxygen evolution diagrams for these Mn 1–x M x PO 4 delithiated cathode materials. To generate the phase diagrams, we use subregular solid-solution models and fit the energetic parameters of these models to density functional theory calculations of special quasirandom structures. The resulting thermodynamic models describe the effect of mixing on the initial temperature of oxygen evolution and on the cumulative amount of oxygen evolution at elevated temperatures. We find that addition of Fe increases the initial temperature and decreases the cumulative amount of oxygen evolution. Mn 0.5 Fe 0.5 PO 4 exhibits an initial temperature 50 °C higher than MnPO 4 and releases 70% less oxygen gas at 300 °C. Al is insoluble in MnPO 4 , so addition of Al has no affect on the initial temperature. However, Al addition does slightly decrease the amount of oxygen evolution due to an inactive AlPO 4 component. Mg and Ni both decrease the initial temperature of oxygen evolution, and therefore may worsen the safety of MnPO 4 .
High-pressure discovery of β-NiBi
We present the discovery of a new high-pressure phase in the Ni–Bi system, β-NiBi.