Stefano Curtarolo

The high-throughput highway to computational materials design

High-throughput computational materials design is an emerging area of materials science. By combining advanced thermodynamic and electronic-structure methods with intelligent data mining and database construction, and exploiting the power of current supercomputer architectures, scientists generate, manage and analyse enormous data repositories for the discovery of novel materials. In this Review we provide a current snapshot of this rapidly evolving field, and highlight the challenges and opportunities that lie ahead.

AFLOW: An Automatic Framework for High-throughput Materials Discovery

Abstract Recent advances in computational materials science present novel opportunities for structure discovery and optimization, including uncovering of unsuspected compounds and metastable structures, electronic structure, surface, and nano-particle properties. The practical realization of these opportunities requires systematic generation and classification of the relevant computational data by high-throughput methods. In this paper we present A flow (Automatic Flow), a software framework for high-throughput calculation of crystal structure properties of alloys, intermetallics and inorganic compounds. The A flow software is available for the scientific community on the website of the materials research consortium, aflowlib.org. Its geometric and electronic structure analysis and manipulation tools are additionally available for online operation at the same website. The combination of automatic methods and user online interfaces provide a powerful tool for efficient quantum computational materials discovery and characterization.

Predicting Crystal Structures with Data Mining of Quantum Calculations

Predicting and characterizing the crystal structure of materials is a key problem in materials research and development. It is typically addressed with highly accurate quantum mechanical computations on a small set of candidate structures, or with empirical rules that have been extracted from a large amount of experimental information, but have limited predictive power. In this Letter, we transfer the concept of heuristic rule extraction to a large library of ab initio calculated information, and we demonstrate that this can be developed into a tool for crystal structure prediction.

Uptake of gases in bundles of carbon nanotubes

Model calculations are presented that predict whether or not an arbitrary gas experiences significant absorption within carbon nanotubes and/or bundles of nanotubes. The potentials used in these calculations assume a conventional form, based on a sum of two-body interactions with individual carbon atoms; the latter employ energy and distance parameters that are derived from empirical combining rules. The results confirm intuitive expectation that small atoms and molecules are absorbed within both the interstitial channels and the tubes, while large atoms and molecules are absorbed almost exclusively within the tubes.

Charting the complete elastic properties of inorganic crystalline compounds

The elastic constant tensor of an inorganic compound provides a complete description of the response of the material to external stresses in the elastic limit. It thus provides fundamental insight into the nature of the bonding in the material, and it is known to correlate with many mechanical properties. Despite the importance of the elastic constant tensor, it has been measured for a very small fraction of all known inorganic compounds, a situation that limits the ability of materials scientists to develop new materials with targeted mechanical responses. To address this deficiency, we present here the largest database of calculated elastic properties for inorganic compounds to date. The database currently contains full elastic information for 1,181 inorganic compounds, and this number is growing steadily. The methods used to develop the database are described, as are results of tests that establish the accuracy of the data. In addition, we document the database format and describe the different ways it can be accessed and analyzed in efforts related to materials discovery and design.

A search model for topological insulators with high-throughput robustness descriptors

Topological insulators (TI) are becoming one of the most studied classes of novel materials because of their great potential for applications ranging from spintronics to quantum computers. To fully integrate TI materials in electronic devices, high-quality epitaxial single-crystalline phases with sufficiently large bulk bandgaps are necessary. Current efforts have relied mostly on costly and time-consuming trial-and-error procedures. Here we show that by defining a reliable and accessible descriptor , which represents the topological robustness or feasibility of the candidate, and by searching the quantum materials repository aflowlib.org, we have automatically discovered 28 TIs (some of them already known) in five different symmetry families. These include peculiar ternary halides, Cs{Sn,Pb,Ge}{Cl,Br,I}(3), which could have been hardly anticipated without high-throughput means. Our search model, by relying on the significance of repositories in materials development, opens new avenues for the discovery of more TIs in different and unexplored classes of systems.

Convergence of multi-valley bands as the electronic origin of high thermoelectric performance in CoSb3 skutterudites.

Filled skutterudites R(x)Co4Sb12 are excellent n-type thermoelectric materials owing to their high electronic mobility and high effective mass, combined with low thermal conductivity associated with the addition of filler atoms into the void site. The favourable electronic band structure in n-type CoSb3 is typically attributed to threefold degeneracy at the conduction band minimum accompanied by linear band behaviour at higher carrier concentrations, which is thought to be related to the increase in effective mass as the doping level increases. Using combined experimental and computational studies, we show instead that a secondary conduction band with 12 conducting carrier pockets (which converges with the primary band at high temperatures) is responsible for the extraordinary thermoelectric performance of n-type CoSb3 skutterudites. A theoretical explanation is also provided as to why the linear (or Kane-type) band feature is not beneficial for thermoelectrics.

All The Catalytic Active Sites of MoS2 for Hydrogen Evolution

MoS2 presents a promising low-cost catalyst for the hydrogen evolution reaction (HER), but the understanding about its active sites has remained limited. Here we present an unambiguous study of the catalytic activities of all possible reaction sites of MoS2, including edge sites, sulfur vacancies, and grain boundaries. We demonstrate that, in addition to the well-known catalytically active edge sites, sulfur vacancies provide another major active site for the HER, while the catalytic activity of grain boundaries is much weaker. The intrinsic turnover frequencies (Tafel slopes) of the edge sites, sulfur vacancies, and grain boundaries are estimated to be 7.5 s-1 (65-75 mV/dec), 3.2 s-1 (65-85 mV/dec), and 0.1 s-1 (120-160 mV/dec), respectively. We also demonstrate that the catalytic activity of sulfur vacancies strongly depends on the density of the vacancies and the local crystalline structure in proximity to the vacancies. Unlike edge sites, whose catalytic activity linearly depends on the length, sulfur vacancies show optimal catalytic activities when the vacancy density is in the range of 7-10%, and the number of sulfur vacancies in high crystalline quality MoS2 is higher than that in low crystalline quality MoS2, which may be related with the proximity of different local crystalline structures to the vacancies.

High-Throughput Combinatorial Database of Electronic Band Structures for Inorganic Scintillator Materials

For the purpose of creating a database of electronic structures of all the known inorganic compounds, we have developed a computational framework based on high-throughput ab initio calculations (AFLOW) and an online repository (www.aflowlib.org). In this article, we report the first step of this task: the calculation of band structures for 7439 compounds intended for the research of scintillator materials for γ-ray radiation detection. Data-mining is performed to select the candidates from 193,456 compounds compiled in the Inorganic Crystal Structure Database. Light yield and scintillation nonproportionality are predicted based on semiempirical band gaps and effective masses. We present a list of materials, potentially bright and proportional, and focus on those exhibiting small effective masses and effective mass ratios.

Finding Unprecedentedly Low-Thermal-Conductivity Half-Heusler Semiconductors via High-Throughput Materials Modeling

The lattice thermal conductivity ({\kappa}{\omega}) is a key property for many potential applications of compounds. Discovery of materials with very low or high {\kappa}{\omega} remains an experimental challenge due to high costs and time-consuming synthesis procedures. High-throughput computational pre-screening is a valuable approach for significantly reducing the set of candidate compounds. In this article, we introduce efficient methods for reliably estimating the bulk {\kappa}{\omega} for a large number of compounds. The algorithms are based on a combination of machine-learning algorithms, physical insights, and automatic ab-initio calculations. We scanned approximately 79,000 half-Heusler entries in the AFLOWLIB.org database. Among the 450 mechanically stable ordered semiconductors identified, we find that {\kappa}{\omega} spans more than two orders of magnitude- a much larger range than that previously thought. {\kappa}{\omega} is lowest for compounds whose elements in equivalent positions have large atomic radii. We then perform a thorough screening of thermodynamical stability that allows to reduce the list to 77 systems. We can then provide a quantitative estimate of {\kappa}{\omega} for this selected range of systems. Three semiconductors having {\kappa}{\omega} < 5 W /(m K) are proposed for further experimental study.