Surface adsorption is a crucial step in numerous processes, including heterogeneous catalysis, where the adsorption of key species is often used as a descriptor of efficiency. We present here an automated adsorption workflow for semiconductors which employs density functional theory calculations to generate adsorption data in a high-throughput manner. Starting from a bulk structure, the workflow performs an exhaustive surface search, followed by an adsorption structure construction step, which generates a minimal energy landscape to determine the optimal adsorbate-surface distance. An extensive set of energy-based, charge-based, geometric, and electronic descriptors tailored toward catalysis research are computed and saved to a personal user database. The application of the workflow to zinc telluride, a promising CO₂ reduction photocatalyst, is presented as a case study to illustrate the capabilities of this method and its potential as a material discovery tool.

Photoelectrocatalysts that use sunlight to power the CO2 reduction reaction will be crucial for carbon-neutral power and energy-efficient industrial processes. Scalable photoelectrocatalysts must satisfy a stringent set of criteria, such as stability under operating conditions, product selectivity, and efficient light absorption. Two-dimensional materials can offer high specific surface area, tunability, and potential for heterostructuring, providing a fresh landscape of candidate catalysts. From a set of promising bulk CO2 reduction photoelectrocatalysts, we screen for candidate monolayers of these materials, then study their catalytic feasibility and suitability. For stable monolayer candidates, we verify the presence of visible-light band gaps, check that band edges can support CO2 reduction, determine exciton binding energies, and compute surface reactivity. We find visible light absorption for SiAs, ZnTe, and ZnSe monolayers, and that due to a lack of binding, CO selectivity is possible. We thus identify SiAs, ZnTe, and ZnSe monolayers as targets for further investigation, expanding the chemical space for CO2 photoreduction candidates.

We present an end-to-end computational system for autonomous materials discovery. The system aims for cost-effective optimization in large, high-dimensional search spaces of materials by adopting a sequential, agent-based approach to deciding which experiments to carry out. In choosing next experiments, agents can make use of past knowledge, surrogate models, logic, thermodynamic or other physical constructs, heuristic rules, and different exploration-exploitation strategies. We show a series of examples for (i) how the discovery campaigns for finding materials satisfying a relative stability objective can be simulated to design new agents, and (ii) how those agents can be deployed in real discovery campaigns to control experiments run externally, such as the cloud-based density functional theory simulations in this work. In a sample set of 16 campaigns covering a range of binary and ternary chemistries including metal oxides, phosphides, sulfides and alloys, this autonomous platform found 383 new stable or nearly stable materials with no intervention by the researchers.

The photocatalytic conversion of the greenhouse gas CO₂ to chemical fuels such as hydrocarbons and alcohols continues to be a promising technology for renewable generation of energy. Major advancements have been made in improving the efficiencies and product selectiveness of currently known CO₂ reduction electrocatalysts, nonetheless, materials discovery is needed to enable economically viable, industrial-scale CO₂ reduction. We report here the largest CO₂ photocathode search to date, starting with 68860 candidate materials, using a rational first-principles computation-based screening strategy to evaluate synthesizability, corrosion resistance, visible-light absorption, and compatibility of the electronic structure with fuel synthesis. The results confirm the observation of the literature that few materials meet the stringent CO₂ photocathode requirements, with only 52 materials meeting all requirements. The results are well validated with respect to the literature, with 9 of these materials having been studied for CO₂ reduction, and the remaining 43 materials are discoveries from our pipeline that merit further investigation.

X-ray absorption spectroscopy (XAS) produces a wealth of information about the local structure of materials, but interpretation of spectra often relies on easily accessible trends and prior assumptions about the structure. Recently, researchers have demonstrated that machine learning models can automate this process to predict the coordinating environments of absorbing atoms from their XAS spectra. However, machine learning models are often difficult to interpret, making it challenging to determine when they are valid and whether they are consistent with physical theories. In this work, we present three main advances to the data-driven analysis of XAS spectra: we demonstrate the efficacy of random forests in solving two new property determination tasks (predicting Bader charge and mean nearest neighbor distance), we address how choices in data representation affect model interpretability and accuracy, and we show that multiscale featurization can elucidate the regions and trends in spectra that encode various local properties. The multiscale featurization transforms the spectrum into a vector of polynomial-fit features, and is contrasted with the commonly-used “pointwise” featurization that directly uses the entire spectrum as input. We find that across thousands of transition metal oxide spectra, the relative importance of features describing the curvature of the spectrum can be localized to individual energy ranges, and we can separate the importance of constant, linear, quadratic, and cubic trends, as well as the white line energy. This work has the potential to assist rigorous theoretical interpretations, expedite experimental data collection, and automate analysis of XAS spectra, thus accelerating the discovery of new functional materials.

Discovering the ideal material for a new application involves determining its numerous properties, such as electronic, mechanical, or thermodynamic, to match those of its desired application. The rise of high-throughput computation has meant that large databases of material properties are now accessible to scientists. However, these databases contain far more information than might appear at first glance, since many relationships exist in the materials science literature to derive, or at least approximate, additional properties. propnet is a new computational framework designed to help scientists to automatically calculate additional information from their datasets. It does this by constructing a network graph of relationships between different materials properties and traversing this graph. Initially, propnet contains a catalog of over 100 property relationships but the hope is for this to expand significantly in the future, and contributions from the community are welcomed.

While the search for catalysts capable of directly converting methane to higher value commodity chemicals and liquid fuels has been active for over a century, a viable industrial process for selective methane activation has yet to be developed. Electronic structure calculations are playing an increasingly relevant role in this search, but large-scale materials screening efforts are hindered by computationally expensive transition state barrier calculations. The purpose of the present letter is twofold. First, we show that, for the wide range of catalysts that proceed via a radical intermediate, a unifying framework for predicting C-H activation barriers using a single universal descriptor can be established. Second, we combine this scaling approach with a thermodynamic analysis of active site formation to provide a map of methane activation rates. Our model successfully rationalizes the available empirical data and lays the foundation for future catalyst design strategies that transcend different catalyst classes.

The discovery and development of novel materials in the field of energy are essential to accelerate the transition to a low-carbon economy. Bringing recent technological innovations in automation, robotics and computer science together with current approaches in chemistry, materials synthesis and characterization will act as a catalyst for revolutionizing traditional research and development in both industry and academia. This Perspective provides a vision for an integrated artificial intelligence approach towards autonomous materials discovery, which, in our opinion, will emerge within the next 5 to 10 years. The approach we discuss requires the integration of the following tools, which have already seen substantial development to date: high-throughput virtual screening, automated synthesis planning, automated laboratories and machine learning algorithms. In addition to reducing the time to deployment of new materials by an order of magnitude, this integrated approach is expected to lower the cost associated with the initial discovery. Thus, the price of the final products (for example, solar panels, batteries and electric vehicles) will also decrease. This in turn will enable industries and governments to meet more ambitious targets in terms of reducing greenhouse gas emissions at a faster pace.

We introduce atomate, an open-source Python framework for computational materials science simulation, analysis, and design with an emphasis on automation and extensibility. Built on top of open source Python packages already in use by the materials community such as pymatgen, FireWorks, and custodian, atomate provides well-tested workflow templates to compute various materials properties such as electronic bandstructure, elastic properties, and piezoelectric, dielectric, and ferroelectric properties. Atomate also enables the computational characterization of materials by providing workflows that calculate X-ray absorption (XAS), Electron energy loss (EELS) and Raman spectra. One of the major features of atomate is that it provides both fully functional workflows as well as reusable components that enable one to compose complex materials science workflows that use a diverse set of computational tools. Additionally, atomate creates output databases that organize the results from individual calculations and contains a builder framework that creates summary reports for each computed material based on multiple simulations.

Electrocatalysis of the oxygen evolution reaction is central to several energy technologies including electrolyzers, solar fuel generators, and air-breathing batteries. Strong acid electrolytes are desirable for many implementations of these technologies, although the deployment of such device designs is often hampered by the lack of non-precious-metal oxygen evolution electrocatalysts, with Ir-based oxides comprising the only known catalysts that exhibit stable activity at low overpotential. During our exploration of the Mn-Sb-O system for precious-metal-free electrocatalysts, we discovered that Mn can be incorporated into the rutile oxide structure at much higher concentrations than previously known, and that these Mn-rich rutile alloys exhibit great catalytic activity with current densities exceeding 50 mA cm-2 at 0.58 V overpotential and catalysis onset at 0.3 V overpotential. While this activity does not surpass that of IrO2, Pourbaix analysis reveals that the Mn-Sb rutile oxide alloys have the same or better thermodynamic stability under operational conditions. By combining combinatorial composition, structure, and activity mapping with synchrotron X-ray absorption measurements and first-principles materials chemistry calculations, we provide a comprehensive understanding of these oxide alloys and identify the critical role of Sb in stabilizing the trivalent Mn octahedra that have been shown to be effective oxygen evolution reaction (OER) catalysts.