CeTi2O6-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 Ce4+ → Ce3+ redox reactions, which explains the high performance and uniqueness of CeO2 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 Ce4+-based oxides which possess thermodynamics superior to CeO2 for STCH production. We first search the Inorganic Crystal Structure Database (ICSD) and literature for Ce4+-based oxides and subsequently use density functional theory to compute their reduction enthalpies (i.e., oxygen vacancy formation energies). We find that CeTi2O6 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 Ce4+ → Ce3+ redox. Our proposed design strategy suggests that further exploration of Ce4+ oxides for STCH production is warranted.

High-pressure discovery of ß-NiBi.

Herein, we present the discovery of a new high-pressure phase in the Ni-Bi system, ß-NiBi, which crystallizes in the TlI structure type. The powerful technique of in situ high-pressure and high-temperature powder X-ray diffraction enabled observation of the formation of ß-NiBi and its reversible reconversion to the ambient pressure phase, α-NiBi.

High-throughput computational design of cathode coatings for Li-ion batteries.

Cathode degradation is a key factor that limits the lifetime of Li-ion batteries. To identify functional coatings that can suppress this degradation, we present a high-throughput density functional theory based framework which consists of reaction models that describe thermodynamic and electrochemical stabilities, and acid-scavenging capabilities of materials. Screening more than 130,000 oxygen-bearing materials, we suggest physical and hydrofluoric-acid barrier coatings such as WO3, LiAl5O8 and ZrP2O7 and hydrofluoric-acid scavengers such as Sc2O3, Li2CaGeO4, LiBO2, Li3NbO4, Mg3(BO3)2 and Li2MgSiO4. Using a design strategy to find the thermodynamically optimal coatings for a cathode, we further present optimal hydrofluoric-acid scavengers such as Li2SrSiO4, Li2CaSiO4 and CaIn2O4 for the layered LiCoO2, and Li2GeO3, Li4NiTeO6 and Li2MnO3 for the spinel LiMn2O4 cathodes. These coating materials have the potential to prolong the cycle-life of Li-ion batteries and surpass the performance of common coatings based on conventional materials such as Al2O3, ZnO, MgO or ZrO2.

Controlling the intercalation chemistry to design high-performance dual-salt hybrid rechargeable batteries.

We have conducted extensive theoretical and experimental investigations to unravel the origin of the electrochemical properties of hybrid Mg(2+)/Li(+) rechargeable batteries at the atomistic and macroscopic levels. By revealing the thermodynamics of Mg(2+) and Li(+) co-insertion into the Mo6S8 cathode host using density functional theory calculations, we show that there is a threshold Li(+) activity for the pristine Mo6S8 cathode to prefer lithiation instead of magnesiation. By precisely controlling the insertion chemistry using a dual-salt electrolyte, we have enabled ultrafast discharge of our battery by achieving 93.6% capacity retention at 20 C and 87.5% at 30 C, respectively, at room temperature.

High-throughput screening of high-capacity electrodes for hybrid Li-ion-Li-O2 cells.

Recent experiments have shown that lithium and oxygen can be electrochemically removed from Li5FeO4 (5Li2O·Fe2O3) and re-accommodated during discharge, creating the possibility of its use as a high-capacity electrode in a hybrid Li-ion/Li-O2 electrochemical cell. Taking this novel chemistry as a model, we use density functional theory (DFT) within a high-throughput framework to screen for analogous reactions in other materials. We search for candidate materials possessing high capacity, voltages compatible with existing electrolytes, and reasonable electrical conductivity. We identify several promising candidate materials that may operate by a similar reaction mechanism and are worthy of investigation, such as Li6MnO4, Li6CoO4, Li4MoO5 and Li8IrO6. This work paves the way for accelerated exploration of this intriguing new battery chemistry.

Origin of the high performance in GeTe-based thermoelectric materials upon Bi2Te3 doping.

As a lead-free material, GeTe has drawn growing attention in thermoelectrics, and a figure of merit (ZT) close to unity was previously obtained via traditional doping/alloying, largely through hole carrier concentration tuning. In this report, we show that a remarkably high ZT of ∼1.9 can be achieved at 773 K in Ge0.87Pb0.13Te upon the introduction of 3 mol % Bi2Te3. Bismuth telluride promotes the solubility of PbTe in the GeTe matrix, thus leading to a significantly reduced thermal conductivity. At the same time, it enhances the thermopower by activating a much higher fraction of charge transport from the highly degenerate Σ valence band, as evidenced by density functional theory calculations. These mechanisms are incorporated and discussed in a three-band (L + Σ + C) model and are found to explain the experimental results well. Analysis of the detailed microstructure (including rhombohedral twin structures) in Ge0.87Pb0.13Te + 3 mol % Bi2Te3 was carried out using transmission electron microscopy and crystallographic group theory. The complex microstructure explains the reduced lattice thermal conductivity and electrical conductivity as well.

Ultralow thermal conductivity and high thermoelectric figure of merit in SnSe crystals.

The thermoelectric effect enables direct and reversible conversion between thermal and electrical energy, and provides a viable route for power generation from waste heat. The efficiency of thermoelectric materials is dictated by the dimensionless figure of merit, ZT (where Z is the figure of merit and T is absolute temperature), which governs the Carnot efficiency for heat conversion. Enhancements above the generally high threshold value of 2.5 have important implications for commercial deployment, especially for compounds free of Pb and Te. Here we report an unprecedented ZT of 2.6 ± 0.3 at 923 K, realized in SnSe single crystals measured along the b axis of the room-temperature orthorhombic unit cell. This material also shows a high ZT of 2.3 ± 0.3 along the c axis but a significantly reduced ZT of 0.8 ± 0.2 along the a axis. We attribute the remarkably high ZT along the b axis to the intrinsically ultralow lattice thermal conductivity in SnSe. The layered structure of SnSe derives from a distorted rock-salt structure, and features anomalously high Grüneisen parameters, which reflect the anharmonic and anisotropic bonding. We attribute the exceptionally low lattice thermal conductivity (0.23 ± 0.03 W m(-1) K(-1) at 973 K) in SnSe to the anharmonicity. These findings highlight alternative strategies to nanostructuring for achieving high thermoelectric performance.

Nonlocal first-principles calculations in Cu-Au and other intermetallic alloys.

Cu-Au is the prototypical alloy system used to exemplify ordering and compound formation, and it serves as a testbed for all new alloy theory methods. Yet, despite the importance of this system, conventional density functional theory (DFT) calculations with semilocal approximations have two dramatic failures in describing the energies of this system: (1) DFT formation energies of the observed Cu3Au and CuAu compounds are nearly a factor of 2 smaller in magnitude than experimental values, and (2) DFT predicts incorrect ordered ground states ground states for Au-rich compositions. Here, we show how modern extensions of DFT based on nonlocal interactions can rectify both of these failures. Our corrections shed light on improving the theoretical predictions for alloy systems to determine accurate formation energies, order-disorder critical temperatures, phase diagrams, and high-throughput computations.

An improved interatomic potential for xenon in UO2: a combined density functional theory/genetic algorithm approach.

We have created an improved xenon interatomic potential for use with existing UO2 potentials. This potential was fit to density functional theory calculations with the Hubbard U correction (DFT + U) using a genetic algorithm approach called iterative potential refinement (IPR). We examine the defect energetics of the IPR-fitted xenon interatomic potential as well as other, previously published xenon potentials. We compare these potentials to DFT + U derived energetics for a series of xenon defects in a variety of incorporation sites (large, intermediate, and small vacant sites). We find the existing xenon potentials overestimate the energy needed to add a xenon atom to a wide set of defect sites representing a range of incorporation sites, including failing to correctly rank the energetics of the small incorporation site defects (xenon in an interstitial and xenon in a uranium site neighboring uranium in an interstitial). These failures are due to problematic descriptions of Xe-O and/or Xe-U interactions of the previous xenon potentials. These failures are corrected by our newly created xenon potential: our IPR-generated potential gives good agreement with DFT + U calculations to which it was not fitted, such as xenon in an interstitial (small incorporation site) and xenon in a double Schottky defect cluster (large incorporation site). Finally, we note that IPR is very flexible and can be applied to a wide variety of potential forms and materials systems, including metals and EAM potentials.

High thermoelectric performance via hierarchical compositionally alloyed nanostructures.

Previous efforts to enhance thermoelectric performance have primarily focused on reduction in lattice thermal conductivity caused by broad-based phonon scattering across multiple length scales. Herein, we demonstrate a design strategy which provides for simultaneous improvement of electrical and thermal properties of p-type PbSe and leads to ZT ~ 1.6 at 923 K, the highest ever reported for a tellurium-free chalcogenide. Our strategy goes beyond the recent ideas of reducing thermal conductivity by adding two key new theory-guided concepts in engineering, both electronic structure and band alignment across nanostructure-matrix interface. Utilizing density functional theory for calculations of valence band energy levels of nanoscale precipitates of CdS, CdSe, ZnS, and ZnSe, we infer favorable valence band alignments between PbSe and compositionally alloyed nanostructures of CdS1-xSex/ZnS1-xSex. Then by alloying Cd on the cation sublattice of PbSe, we tailor the electronic structure of its two valence bands (light hole L and heavy hole Σ) to move closer in energy, thereby enabling the enhancement of the Seebeck coefficients and the power factor.

Relative stability of normal vs. inverse spinel for 3d transition metal oxides as lithium intercalation cathodes.

Spinel oxides represent an important class of cathode materials for Li-ion batteries. Two major variants of the spinel crystal structure are normal and inverse. The relative stability of normal and inverse ordering at different stages of lithiation has important consequences in lithium diffusivity, voltage, capacity retention and battery life. In this paper, we investigate the relative structural stability of normal and inverse structures of the 3d transition metal oxide spinels with first-principles DFT calculations. We have considered ternary spinel oxides LixM2O4 with M = Ti, V, Cr, Mn, Fe, Co and Ni in both lithiated (x = 1) and delithiated (x = 0) conditions. We find that for all lithiated spinels, the normal structure is preferred regardless of the metal. We observe that the normal structure for all these oxides has a lower size mismatch between octahedral cations compared to the inverse structure. With delithiation, many of the oxides undergo a change in stability with vanadium in particular, showing a tendency to occupy tetrahedral sites. We find that in the delithiated oxide, only vanadium ions can access a +5 oxidation state which prefers tetrahedral coordination. We have also calculated the average voltage of lithiation for these spinels. The calculated voltages agree well with the previously measured and calculated values, wherever available. For the yet to be characterized spinels, our calculation provides voltage values which can motivate further experimental attention. Lastly, we observe that all the normal spinel oxides of the 3d transition metal series have a driving force for a transformation to the non-spinel structure upon delithiation.

A hybrid computational-experimental approach for automated crystal structure solution.

Crystal structure solution from diffraction experiments is one of the most fundamental tasks in materials science, chemistry, physics and geology. Unfortunately, numerous factors render this process labour intensive and error prone. Experimental conditions, such as high pressure or structural metastability, often complicate characterization. Furthermore, many materials of great modern interest, such as batteries and hydrogen storage media, contain light elements such as Li and H that only weakly scatter X-rays. Finally, structural refinements generally require significant human input and intuition, as they rely on good initial guesses for the target structure. To address these many challenges, we demonstrate a new hybrid approach, first-principles-assisted structure solution (FPASS), which combines experimental diffraction data, statistical symmetry information and first-principles-based algorithmic optimization to automatically solve crystal structures. We demonstrate the broad utility of FPASS to clarify four important crystal structure debates: the hydrogen storage candidates MgNH and NH(3)BH(3); Li(2)O(2), relevant to Li-air batteries; and high-pressure silane, SiH(4).

Morphology control of nanostructures: Na-doped PbTe-PbS system.

The morphology of crystalline precipitates in a solid-state matrix is governed by complex but tractable energetic considerations driven largely by volume strain energy minimization and anisotropy of interfacial energies. Spherical precipitate morphologies are favored by isotropic systems, while anisotropic interfacial energies give energetic preference to certain crystallographically oriented interfaces, resulting in a faceted precipitate morphology. In conventional solid-solution precipitation, a precipitate's morphological evolution is mediated by surface anchoring of capping molecules, which dramatically alter the surface energy in an anisotropic manner, thereby providing exquisite morphology control during crystal growth. Herein, we present experimental evidence and theoretical validation for the role of a ternary element (Na) in controlling the morphology of nanoscale PbS crystals nucleating in a PbTe matrix, an important bulk thermoelectric system. The PbS nanostructures formed by phase separation from a PbI(2)-doped or undoped PbTe matrix have irregular morphologies. However, replacing the iodine dopant with Na (1-2 mol %) alters dramatically the morphology of the PbS precipitates. Segregation of Na at PbTe/PbS interfaces result in cuboidal and truncated cuboidal morphologies for PbS. Using analytical scanning/transmission electron microscopy and atom-probe tomography, we demonstrate unambiguously that Na partitions to the precipitates and segregates at the matrix/precipitate interfaces, inducing morphological anisotropy of PbS precipitates. First-principles and semiclassical calculations reveal that Na as a solute in PbTe has a higher energy than in PbS and that Na segregation at a (100) PbTe/PbS interface decreases the total energy of matrix/precipitate system, resulting in faceting of PbS precipitates. These results provide an impetus for a new strategy for controlling morphological evolution in matrix/precipitate systems, mediated by solute partitioning of ternary additions.

Physics of bandgap formation in Cu-Sb-Se based novel thermoelectrics: the role of Sb valency and Cu d levels.

In this paper we discuss the results of ab initio electronic structure calculations for Cu(3)SbSe(4) (Se4) and Cu(3)SbSe(3) (Se3), two narrow bandgap semiconductors of thermoelectric interest. We find that Sb is trivalent in both the compounds, in contrast to a simple nominal valence (ionic) picture which suggests that Sb should be 5 + in Se4. The gap formation in Se4 is quite subtle, with hybridization between Sb 5s and the neighboring Se 4s, 4p orbitals, position of Cu d states, and non-local exchange interaction, each playing significant roles. Thermopower calculations show that Se4 is a better p-type system. Our theoretical results for Se4 agree very well with recent experimental results obtained by Skoug et al (2011 Sci. Adv. Mater. 3 602).

Raising the thermoelectric performance of p-type PbS with endotaxial nanostructuring and valence-band offset engineering using CdS and ZnS.

We have investigated in detail the effect of CdS and ZnS as second phases on the thermoelectric properties of p-type PbS. We report a ZT of ~1.3 at 923 K for 2.5 at.% Na-doped p-type PbS with endotaxially nanostructured 3.0 at.% CdS. We attribute the high ZT to the combination of broad-based phonon scattering on multiple length scales to reduce (lattice) thermal conductivity and favorable charge transport through coherent interfaces between the PbS matrix and metal sulfide nanophase precipitates, which maintains the requisite high carrier conductivity and the associated power factor. Similar to large ionically bonded metal sulfides (ZnS, CaS, and SrS), the covalently bonded CdS can also effectively reduce the lattice thermal conductivity in p-type PbS. The presence of ubiquitous nanostructuring was confirmed by transmission electron microscopy. Valence and conduction band energy levels of the NaCl-type metal sulfides, MS (M = Pb, Cd, Zn, Ca, and Sr) were calculated from density functional theory to gain insight into the band alignment between PbS and the second phases in these materials. The hole transport is controlled by band offset minimization through the alignment of valence bands between the host PbS and the embedded second phases, MS (M = Cd, Zn, Ca, and Sr). The smallest valence band offset of about 0.13 eV at 0 K was found between PbS and CdS which is diminished further by thermal band broadening at elevated temperature. This allows carrier transport between the endotaxially aligned components (i.e., matrix and nanostructure), thus minimizing significant deterioration of the hole mobility and power factor. We conclude the thermoelectric performance of the PbS system and, by extension, other systems can be enhanced by means of a closely coupled phonon-blocking/electron-transmitting approach through embedding endotaxially nanostructured second phases.

First principles simulations of the electrochemical lithiation and delithiation of faceted crystalline silicon.

Silicon is of significant interest as a next-generation anode material for lithium-ion batteries due to its extremely high capacity. The reaction of lithium with crystalline silicon is known to present a rich range of phenomena, including electrochemical solid state amorphization, crystallization at full lithiation of a Li(15)Si(4) phase, hysteresis in the first lithiation-delithiation cycle, and highly anisotropic lithiation in crystalline samples. Very little is known about these processes at an atomistic level, however. To provide fundamental insights into these issues, we develop and apply a first principles, history-dependent, lithium insertion and removal algorithm to model the process of lithiation and subsequent delithiation of crystalline Si. The simulations give a realistic atomistic picture of lithiation demonstrating, for the first time, the amorphization process and hinting at the formation of the Li(15)Si(4) phase. Voltages obtained from the simulations show that lithiation of the (110) surface is thermodynamically more favorable than lithiation of the (100) or (111) surfaces, providing an explanation for the drastic lithiation anisotropy seen in experiments on Si micro- and nanostructures. Analysis of the delithiation and relithiation processes also provides insights into the underlying physics of the lithiation-delithiation hysteresis, thus providing firm conceptual foundations for future design of improved Si-based anodes for Li ion battery applications.

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.

Performance of Cluster Expansions of Coverage-Dependent Adsorption of Atomic Oxygen on Pt(111).

A density functional theory (DFT) database of 66 Pt(111)/O formation energies is presented. We fit this database of formation energies to a range of cluster expansions (CEs) of systematically increasing size and flexibility. We find that the performance of the CE depends upon the property or properties of interest. Pair-wise CEs with up to third nearest neighbor interactions poorly predict all metrics. CEs with five to eight pairwise interactions and one to two triplet interactions predicted formation energies and most ground states accurately but predicted average and differential adsorption energies with modest errors. A larger CE captures average and differential adsorption energies as well as formation energies and ground states. The choice of figures in the CEs is also examined. Pair-wise figures and the linear, 1-1-3, triplet are necessary to obtain CEs that qualitatively reproduce the examined properties; however, other figures are more interchangeable. The electronic and strain components of the adsorbate-adsorbate interactions is studied by comparing a CE of DFT formation energies in which atoms were not allowed to relax to the CEs of the relaxed surface. On an unrelaxed Pt surface, interactions are shorter-ranged interactions and more repulsive at first nearest neighbor separation.

First-principles prediction of thermodynamically reversible hydrogen storage reactions in the Li-Mg-Ca-B-H system.

Introduction of economically viable hydrogen cars is hindered by the need to store large amounts of hydrogen. Metal borohydrides [LiBH(4), Mg(BH(4))(2), Ca(BH(4))(2)] are attractive candidates for onboard storage because they contain high densities of hydrogen by weight and by volume. Using a set of recently developed theoretical first-principles methods, we predict currently unknown crystal structures and hydrogen storage reactions in the Li-Mg-Ca-B-H system. Hydrogen release from LiBH(4) and Mg(BH(4))(2) is predicted to proceed via intermediate Li(2)B(12)H(12) and MgB(12)H(12) phases, while for Ca borohydride two competing reaction pathways (into CaB(6) and CaH(2), and into CaB(12)H(12) and CaH(2)) are found to have nearly equal free energies. We predict two new hydrogen storage reactions that are some of the most attractive among the presently known ones. They combine high gravimetric densities (8.4 and 7.7 wt % H(2)) with low enthalpies [approximately 25 kJ/(mol H(2))] and are thermodynamically reversible at low pressures due to low vibrational entropies of the product phases containing the [B(12)H(12)](2-) anion.