Clamping enables enhanced electromechanical responses in antiferroelectric thin films.

Thin-film materials with large electromechanical responses are fundamental enablers of next-generation micro-/nano-electromechanical applications. Conventional electromechanical materials (for example, ferroelectrics and relaxors), however, exhibit severely degraded responses when scaled down to submicrometre-thick films due to substrate constraints (clamping). This limitation is overcome, and substantial electromechanical responses in antiferroelectric thin films are achieved through an unconventional coupling of the field-induced antiferroelectric-to-ferroelectric phase transition and the substrate constraints. A detilting of the oxygen octahedra and lattice-volume expansion in all dimensions are observed commensurate with the phase transition using operando electron microscopy, such that the in-plane clamping further enhances the out-of-plane expansion, as rationalized using first-principles calculations. In turn, a non-traditional thickness scaling is realized wherein an electromechanical strain (1.7%) is produced from a model antiferroelectric PbZrO3 film that is just 100 nm thick. The high performance and understanding of the mechanism provide a promising pathway to develop high-performance micro-/nano-electromechanical systems.

Ferroelectricity and multiferroicity in anti-Ruddlesden-Popper structures.

Combining ferroelectricity with other properties such as visible light absorption or long-range magnetic order requires the discovery of new families of ferroelectric materials. Here, through the analysis of a high-throughput database of phonon band structures, we identify a structural family of anti-Ruddlesden-Popper phases [Formula: see text]O (A=Ca, Sr, Ba, Eu, X=Sb, P, As, Bi) showing ferroelectric and antiferroelectric behaviors. The discovered ferroelectrics belong to the new class of hyperferroelectrics that polarize even under open-circuit boundary conditions. The polar distortion involves the movement of O anions against apical A cations and is driven by geometric effects resulting from internal chemical strains. Within this structural family, we show that [Formula: see text]O combines coupled ferromagnetic and ferroelectric order at the same atomic site, a very rare occurrence in materials physics.

LOBSTER: Local orbital projections, atomic charges, and chemical-bonding analysis from projector-augmented-wave-based density-functional theory.

We present an update on recently developed methodology and functionality in the computer program Local Orbital Basis Suite Toward Electronic-Structure Reconstruction (LOBSTER) for chemical-bonding analysis in periodic systems. LOBSTER is based on an analytic projection from projector-augmented wave (PAW) density-functional theory (DFT) computations (Maintz et al., J. Comput. Chem. 2013, 34, 2557), reconstructing chemical information in terms of local, auxiliary atomic orbitals and thereby opening the output of PAW-based DFT codes to chemical interpretation. We demonstrate how LOBSTER has been improved by taking into account time-reversal symmetry, thereby speeding up the DFT and LOBSTER calculations by a factor of 2. Over the recent years, the functionalities have also been continually expanded, including accurate projected densities of states (DOSs), crystal orbital Hamilton population (COHP) analysis, atomic and orbital charges, gross populations, and the recently introduced k-dependent COHP. The software is offered free-of-charge for non-commercial research.

Electron-Phonon beyond Fröhlich: Dynamical Quadrupoles in Polar and Covalent Solids.

We include the treatment of quadrupolar fields beyond the Fröhlich interaction in the first-principles electron-phonon vertex in semiconductors. Such quadrupolar fields induce long-range interactions that have to be taken into account for accurate physical results. We apply our formalism to Si (nonpolar), GaAs, and GaP (polar) and demonstrate that electron mobilities show large errors if dynamical quadrupoles are not properly treated.

Combining phonon accuracy with high transferability in Gaussian approximation potential models.

Machine learning driven interatomic potentials, including Gaussian approximation potential (GAP) models, are emerging tools for atomistic simulations. Here, we address the methodological question of how one can fit GAP models that accurately predict vibrational properties in specific regions of configuration space while retaining flexibility and transferability to others. We use an adaptive regularization of the GAP fit that scales with the absolute force magnitude on any given atom, thereby exploring the Bayesian interpretation of GAP regularization as an "expected error" and its impact on the prediction of physical properties for a material of interest. The approach enables excellent predictions of phonon modes (to within 0.1 THz-0.2 THz) for structurally diverse silicon allotropes, and it can be coupled with existing fitting databases for high transferability across different regions of configuration space, which we demonstrate for liquid and amorphous silicon. These findings and workflows are expected to be useful for GAP-driven materials modeling more generally.

ABINIT: Overview and focus on selected capabilities.

abinit is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe-Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the "temperature-dependent effective potential" approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which abinit relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The abinit DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. abinit has strong links with many other software projects that are briefly mentioned.

The Limited Predictive Power of the Pauling Rules.

The Pauling rules have been used for decades to rationalise the crystal structures of ionic compounds. Despite their importance, there has been no statistical assessment of the performances of these five empirical rules so far. Here, we rigorously and automatically test all five Pauling rules for a large data set of around 5000 known oxides. We discuss each Pauling rule separately, stressing their limits and range of application in terms of chemistries and structures. We conclude that only 13 % of the oxides simultaneously satisfy the last four rules, indicating a much lower predictive power than expected.

Sr3CrN3: A New Electride with a Partially Filled d-Shell Transition Metal.

Electrides are ionic crystals in which the electrons prefer to occupy free space, serving as anions. Because the electrons prefer to be in the pockets, channels, or layers to the atomic orbitals around the nuclei, it has been challenging to find electrides with partially filled d-shell transition metals, since an unoccupied d-shell provides an energetically favorable location for the electrons to occupy. We recently predicted the existence of electrides with partially filled d-shells using high-throughput computational screening. Here, we provide experimental support using X-ray absorption spectroscopy and X-ray and neutron diffraction to show that Sr3CrN3 is indeed an electride despite its partial d-shell configuration. Our findings indicate that Sr3CrN3 is the first known electride with a partially filled d-shell transition metal, in agreement with theory, which significantly broadens the criteria for the search for new electride materials.

Measurements of Oxygen Electroadsorption Energies and Oxygen Evolution Reaction on RuO2(110): A Discussion of the Sabatier Principle and Its Role in Electrocatalysis.

We report the hydroxide (OHad) and oxide (Oad) experimental electroadsorption free energies, their dependences on pH, and their correlations to the oxygen evolution reaction (OER) electrocatalysis on RuO2(110) surface. The Sabatier principle predicts that catalyst is most active when the intermediate stabilization is moderate, not too strong such that the bound intermediate disrupts the subsequent catalytic cycle, nor too weak such that the surface is ineffective. For decades, researchers have used this concept to rationalize the activity trend of many OER electrocatalysts including RuO2, which is among the state-of-the-art OER catalysts. In this article, we report an experimental assessment of the Sabatier principle by comparing the oxygen electroadsorption energy to the OER electrocatalysis for the first time on RuO2. We find that the OHad and Oad electroadsorption energies on RuO2(110) depend on pH and obey the scaling relation. However, we did not observe a direct correlation between the OHad and Oad electroadsorption energies and the OER activity in the comparative analysis that includes both RuO2(110) and IrO2(110). Our result raises a question of whether the Sabatier principle can describe highly active electrocatalysts, where the kinetic aspects may influence the electrocatalysis more strongly than the electroadsorption energy, which captures only the thermodynamics of the intermediates and not yet kinetics.

Influence of Surface Adsorption on the Oxygen Evolution Reaction on IrO2(110).

A catalyst functions by stabilizing reaction intermediates, usually through surface adsorption. In the oxygen evolution reaction (OER), surface oxygen adsorption plays an indispensable role in the electrocatalysis. The relationship between the adsorption energetics and OER kinetics, however, has not yet been experimentally measured. Herein we report an experimental relationship between the adsorption of surface oxygen and the kinetics of the OER on IrO2(110) epitaxially grown on a TiO2(110) single crystal. The high quality of the IrO2 film grown using molecular-beam epitaxy affords the ability to extract the surface oxygen adsorption and its impact on the OER. By examining a series of electrolytes, we find that the adsorption energy changes linearly with pH, which we attribute to the electrified interfacial water. We support this hypothesis by showing that an electrolyte salt modification can lead to an adsorption energy shift. The dependence of the adsorption energy on pH has implications for the OER kinetics, but it is not the only factor; the dependence of the OER electrocatalysis on pH stipulates two OER mechanisms, one operating in acidic solution and another operating in alkaline solution. Our work points to the subtle adsorption-kinetics relationship in the OER and highlights the importance of the interfacial electrified interaction in electrocatalyst design.

A computational assessment of the electronic, thermoelectric, and defect properties of bournonite (CuPbSbS3) and related substitutions.

Bournonite (CuPbSbS3) is an earth-abundant mineral with potential thermoelectric applications. This material has a complex crystal structure (space group Pmn21 #31) and has previously been measured to exhibit a very low thermal conductivity (κ < 1 W m-1 K-1 at T ≥ 300 K). In this study, we employ high-throughput density functional theory calculations to investigate how the properties of the bournonite crystal structure change with elemental substitutions. Specifically, we compute the stability and electronic properties of 320 structures generated via substitutions {Na-K-Cu-Ag}{Si-Ge-Sn-Pb}{N-P-As-Sb-Bi}{O-S-Se-Te} in the ABCD3 formula. We perform two types of transport calculations: the BoltzTraP model, which has been extensively tested, and a newer AMSET model that we have developed and which incorporates scattering effects. We discuss the differences in the model results, finding qualitative agreement except in the case of degenerate bands. Based on our calculations, we identify p-type CuPbSbSe3, CuSnSbSe3 and CuPbAsSe3 as potentially promising materials for further investigation. We additionally calculate the defect properties, finding that n-type behavior in bournonite and the selected materials is highly unlikely, and p-type behavior might be enhanced by employing Sb-poor synthesis conditions to prevent the formation of SbPb defects. Finally, we discuss the origins of various trends with chemical substitution, including the possible role of stereochemically active lone pair effects in stabilizing the bournonite structure and the effect of cation and anion selection on the calculated band gap.

Thinking Like a Chemist: Intuition in Thermoelectric Materials.

The coupled transport properties required to create an efficient thermoelectric material necessitates a thorough understanding of the relationship between the chemistry and physics in a solid. We approach thermoelectric material design using the chemical intuition provided by molecular orbital diagrams, tight binding theory, and a classic understanding of bond strength. Concepts such as electronegativity, band width, orbital overlap, bond energy, and bond length are used to explain trends in electronic properties such as the magnitude and temperature dependence of band gap, carrier effective mass, and band degeneracy and convergence. The lattice thermal conductivity is discussed in relation to the crystal structure and bond strength, with emphasis on the importance of bond length. We provide an overview of how symmetry and bonding strength affect electron and phonon transport in solids, and how altering these properties may be used in strategies to improve thermoelectric performance.

Low-dimensional transport and large thermoelectric power factors in bulk semiconductors by band engineering of highly directional electronic states.

Thermoelectrics are promising for addressing energy issues but their exploitation is still hampered by low efficiencies. So far, much improvement has been achieved by reducing the thermal conductivity but less by maximizing the power factor. The latter imposes apparently conflicting requirements on the band structure: a narrow energy distribution and a low effective mass. Quantum confinement in nanostructures and the introduction of resonant states were suggested as possible solutions to this paradox, but with limited success. Here, we propose an original approach to fulfill both requirements in bulk semiconductors. It exploits the highly directional character of some orbitals to engineer the band structure and produce a type of low-dimensional transport similar to that targeted in nanostructures, while retaining isotropic properties. Using first-principle calculations, the theoretical concept is demonstrated in Fe2YZ Heusler compounds, yielding power factors 4 to 5 times larger than in classical thermoelectrics at room temperature. Our findings are totally generic and rationalize the search of alternative compounds with similar behavior. Beyond thermoelectricity, these might be relevant also in the context of electronic, superconducting, or photovoltaic applications.

Relating voltage and thermal safety in Li-ion battery cathodes: a high-throughput computational study.

High voltage and high thermal safety are desirable characteristics of cathode materials, but difficult to achieve simultaneously. This work uses high-throughput density functional theory computations to evaluate the link between voltage and safety (as estimated by thermodynamic O2 release temperatures) for over 1400 cathode materials. Our study indicates that a strong inverse relationship exists between voltage and safety: just over half the variance in O2 release temperature can be explained by voltage alone. We examine the effect of polyanion group, redox couple, and ratio of oxygen to counter-cation on both voltage and safety. As expected, our data demonstrates that polyanion groups improve safety when comparing compounds with similar voltages. However, a counterintuitive result of our study is that polyanion groups produce either no benefit or reduce safety when comparing compounds with the same redox couple. Using our data set, we tabulate voltages and oxidation potentials for over 105 combinations of redox couple/anion, which can be used towards the design and rationalization of new cathode materials. Overall, only a few compounds in our study, representing limited redox couple/polyanion combinations, exhibit both high voltage and high safety. We discuss these compounds in more detail as well as the opportunities for designing safe, high-voltage cathodes.

Data mining approaches to high-throughput crystal structure and compound prediction.

Predicting unknown inorganic compounds and their crystal structure is a critical step of high-throughput computational materials design and discovery. One way to achieve efficient compound prediction is to use data mining or machine learning methods. In this chapter we present a few algorithms for data mining compound prediction and their applications to different materials discovery problems. In particular, the patterns or correlations governing phase stability for experimental or computational inorganic compound databases are statistically learned and used to build probabilistic or regression models to identify novel compounds and their crystal structures. The stability of those compound candidates is then assessed using ab initio techniques. Finally, we report a few cases where data mining driven computational predictions were experimentally confirmed through inorganic synthesis.

The Defects Genome of Janus Transition Metal Dichalcogenides.

2D Janus Transition Metal Dichalcogenides (TMDs) have attracted much interest due to their exciting quantum properties arising from their unique two-faced structure, broken-mirror symmetry, and consequent colossal polarization field within the monolayer. While efforts are made to achieve high-quality Janus monolayers, the existing methods rely on highly energetic processes that introduce unwanted grain-boundary and point defects with still unexplored effects on the material's structural and excitonic properties Through high-resolution scanning transmission electron microscopy (HRSTEM), density functional theory (DFT), and optical spectroscopy measurements; this work introduces the most encountered and energetically stable point defects. It establishes their impact on the material's optical properties. HRSTEM studies show that the most energetically stable point defects are single (VS and VSe) and double chalcogen vacancy (VS -VSe), interstitial defects (Mi), and metal impurities (MW) and establish their structural characteristics. DFT further establishes their formation energies and related localized bands within the forbidden band. Cryogenic excitonic studies on h-BN-encapsulated Janus monolayers offer a clear correlation between these structural defects and observed emission features, which closely align with the results of the theory. The overall results introduce the defect genome of Janus TMDs as an essential guideline for assessing their structural quality and device properties.

Small-pore hydridic frameworks store densely packed hydrogen.

Nanoporous materials have attracted great attention for gas storage, but achieving high volumetric storage capacity remains a challenge. Here, by using neutron powder diffraction, volumetric gas adsorption, inelastic neutron scattering and first-principles calculations, we investigate a magnesium borohydride framework that has small pores and a partially negatively charged non-flat interior for hydrogen and nitrogen uptake. Hydrogen and nitrogen occupy distinctly different adsorption sites in the pores, with very different limiting capacities of 2.33 H2 and 0.66 N2 per Mg(BH4)2. Molecular hydrogen is packed extremely densely, with about twice the density of liquid hydrogen (144 g H2 per litre of pore volume). We found a penta-dihydrogen cluster where H2 molecules in one position have rotational freedom, whereas H2 molecules in another position have a well-defined orientation and a directional interaction with the framework. This study reveals that densely packed hydrogen can be stabilized in small-pore materials at ambient pressures.

A substitutional quantum defect in WS2 discovered by high-throughput computational screening and fabricated by site-selective STM manipulation.

Point defects in two-dimensional materials are of key interest for quantum information science. However, the parameter space of possible defects is immense, making the identification of high-performance quantum defects very challenging. Here, we perform high-throughput (HT) first-principles computational screening to search for promising quantum defects within WS2, which present localized levels in the band gap that can lead to bright optical transitions in the visible or telecom regime. Our computed database spans more than 700 charged defects formed through substitution on the tungsten or sulfur site. We found that sulfur substitutions enable the most promising quantum defects. We computationally identify the neutral cobalt substitution to sulfur (Co S 0 ) and fabricate it with scanning tunneling microscopy (STM). The Co S 0 electronic structure measured by STM agrees with first principles and showcases an attractive quantum defect. Our work shows how HT computational screening and nanoscale synthesis routes can be combined to design promising quantum defects.

A map of single-phase high-entropy alloys.

High-entropy alloys have exhibited unusual materials properties. The stability of equimolar single-phase solid solution of five or more elements is supposedly rare and identifying the existence of such alloys has been challenging because of the vast chemical space of possible combinations. Herein, based on high-throughput density-functional theory calculations, we construct a chemical map of single-phase equimolar high-entropy alloys by investigating over 658,000 equimolar quinary alloys through a binary regular solid-solution model. We identify 30,201 potential single-phase equimolar alloys (5% of the possible combinations) forming mainly in body-centered cubic structures. We unveil the chemistries that are likely to form high-entropy alloys, and identify the complex interplay among mixing enthalpy, intermetallics formation, and melting point that drives the formation of these solid solutions. We demonstrate the power of our method by predicting the existence of two new high-entropy alloys, i.e. the body-centered cubic AlCoMnNiV and the face-centered cubic CoFeMnNiZn, which are successfully synthesized.

Limits to Hole Mobility and Doping in Copper Iodide.

Over one hundred years have passed since the discovery of the p-type transparent conducting material copper iodide, predating the concept of the "electron-hole" itself. Supercentenarian status notwithstanding, little is understood about the charge transport mechanisms in CuI. Herein, a variety of modeling techniques are used to investigate the charge transport properties of CuI, and limitations to the hole mobility over experimentally achievable carrier concentrations are discussed. Poor dielectric response is responsible for extensive scattering from ionized impurities at degenerately doped carrier concentrations, while phonon scattering is found to dominate at lower carrier concentrations. A phonon-limited hole mobility of 162 cm2 V-1 s-1 is predicted at room temperature. The simulated charge transport properties for CuI are compared to existing experimental data, and the implications for future device performance are discussed. In addition to charge transport calculations, the defect chemistry of CuI is investigated with hybrid functionals, revealing that reasonably localized holes from the copper vacancy are the predominant source of charge carriers. The chalcogens S and Se are investigated as extrinsic dopants, where it is found that despite relatively low defect formation energies, they are unlikely to act as efficient electron acceptors due to the strong localization of holes and subsequent deep transition levels.