Source author record

Stefano Curtarolo

Stefano Curtarolo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

38works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2026arXiv

Terahertz volume plasmon-polariton modulation in all-dielectric hyperbolic metamaterials

The development of plasmonics and related applications in the terahertz range faces limitations due to the intrinsic high electron density of standard metals. All-dielectric systems are profitable alternatives, which allows for customized modulation of the optical response upon doping. Here we focus on plasmon-based hyperbolic metamaterials realized stacking doped III-V semiconductors that have been shown to be optically active in the terahertz spectral region. By using a multi-physics multi-scale theoretical approach, we unravel the role of doping and geometrical characteristics (e.g., thickness, composition, grating) in the modulation of high-k plasmon-polariton modes across the metamaterial.

preprint2022arXiv

Physics in the Machine: Integrating Physical Knowledge in Autonomous Phase-Mapping

Application of artificial intelligence (AI), and more specifically machine learning, to the physical sciences has expanded significantly over the past decades. In particular, science-informed AI, also known as scientific AI or inductive bias AI, has grown from a focus on data analysis to now controlling experiment design, simulation, execution and analysis in closed-loop autonomous systems. The CAMEO (closed-loop autonomous materials exploration and optimization) algorithm employs scientific AI to address two tasks: learning a material system's composition-structure relationship and identifying materials compositions with optimal functional properties. By integrating these, accelerated materials screening across compositional phase diagrams was demonstrated, resulting in the discovery of a best-in-class phase change memory material. Key to this success is the ability to guide subsequent measurements to maximize knowledge of the composition-structure relationship, or phase map. In this work we investigate the benefits of incorporating varying levels of prior physical knowledge into CAMEO's autonomous phase-mapping. This includes the use of ab-initio phase boundary data from the AFLOW repositories, which has been shown to optimize CAMEO's search when used as a prior.

preprint2022arXiv

The Microscopic Diamond Anvil Cell: Stabilization of Superhard, Superconducting Carbon Allotropes at Ambient Pressure

A metallic covalently bonded carbon allotrope is predicted via first principles calculations. It is composed of an $sp^3$ carbon framework that acts as a diamond anvil cell by constraining the distance between parallel cis-polyacetylene chains. The distance between these $sp^2$ carbon atoms renders the phase metallic, and yields two well-nested nearly parallel bands that span the Fermi level. Calculations show that this phase is a conventional superconductor, with the motions of the $sp^2$ carbons being key contributors to the electron phonon coupling. The $sp^3$ carbon atoms impart superior mechanical properties, with a predicted Vickers hardness of 48~GPa. This phase, metastable at ambient conditions, could be made via cold compression of graphite to 40~GPa. A family of multifunctional materials with tunable superconducting and mechanical properties could be derived from this phase by varying the $sp^2$ versus $sp^3$ carbon content and by doping.

preprint2021arXiv

Data-driven quest for two-dimensional non-van der Waals materials

Two-dimensional (2D) materials are frequently associated with the sheets which form bulk layered compounds bonded by van der Waals (vdW) forces. The anisotropy and weak interaction between the sheets have also been the main criteria in the computational search for new 2D systems, which predicted about 2000 exfoliable compounds. However, several representatives of a new type of non-vdW 2D systems, such as hematene or ilmenene, which have no layered 3D analogues, and which, unlike, e.g. silicene, do not need to strongly interact with the substrate to be stable, were recently manufactured. The family of non-vdW 2D materials is an attractive playground for data-driven high-throughput approaches as computational design principles are still missing. Here, we outline a new set of 8 binary and 20 ternary candidates by filtering the AFLOW-ICSD database according to the structural prototype of the first two template systems realized in experiment. All materials show a strong structural relaxation upon confinement to 2D which is essential to correctly estimate the bonding strength between facets. The oxidation state of the cations at the surface of the sheets is demonstrated to regulate the inter-facet binding energy with low oxidation states leading to weak bonding. We anticipate this descriptor to be highly useful to obtain novel 2D materials, providing clear guidelines for experiments. Our calculations further indicate that the materials exhibit a vast range of appealing electronic, optical and magnetic properties, which is expected to also make them attractive for potential applications particularly in spintronics.

preprint2020arXiv

Ultrathin SnTe films as a route towards all-in-one spintronics devices

Spin transistors based on a semiconducting channel attached to ferromagnetic electrodes suffer from fast spin decay and extremely low spin injection/detection efficiencies. Here, we propose an alternative all-in-one spin device whose operation principle relies on electric manipulation of the spin lifetime in two-dimensional (2D) SnTe, in which the sizable spin Hall effect eliminates the need for using ferromagnets. In particular, we explore the persistent spin texture (PST) intrinsically present in the ferroelectric phase which protects the spin from decoherence and supports extraordinarily long spin lifetime. Our first-principles calculations followed by symmetry arguments revealed that such a spin wave mode can be externally detuned by perpendicular electric field, leading to spin randomization and decrease in spin lifetime. We further extend our analysis to ultrathin SnTe films and confirm the emergence of PST as well as a moderate enhancement of intrinsic spin Hall conductivity. The recent room-temperature observation of the ferroelectric phase in 2D-SnTe suggests that novel all-electric spintronics devices are within reach.

preprint2016arXiv

A computational high-throughput search for new ternary superalloys

In 2006, a novel cobalt-based superalloy was discovered [1] with mechanical properties better than some conventional nickel-based superalloys. As with conventional superalloys, its high performance arises from the precipitate-hardening effect of a coherent L1$_2$ phase, which is in two-phase equilibrium with the fcc matrix. Inspired by this unexpected discovery of an L1$_2$ ternary phase, we performed a first-principles search through 2224 ternary metallic systems for analogous precipitate-hardening phases of the form $X_{3}$[$A_{0.5}, B_{0.5}$], where $X$ = Ni, Co, or Fe, and [$A,B$] = Li, Be, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn Ga, Sr, Y, Zr, Nb, Mo, Tc, Ru, Rh, Pd, Ag, Cd, In, Sn, Sb, Hf, Ta, W, Re, Os, Ir, Pt, Au, Hg, or Tl. We found 102 systems that have a smaller decomposition energy and a lower formation enthalpy than the Co$_{3}$(Al, W) superalloy. They have a stable two-phase equilibrium with the host matrix within the concentration range $0<x<1$ ($X_{3}$[$A_{x}, B_{1-x}$]) and have a relative lattice mismatch with the host matrix of less than or equal to 5%. These new candidates, narrowed from 2224 systems, suggest possible experimental exploration for identifying new superalloys. Of these 102 systems, 37 are new; they have no reported phase diagrams in standard databases. Based on cost, experimental difficulty, and toxicity, we limit these 37 to a shorter list of six promising candidates of immediate interest. Our calculations are consistent with current experimental literature where data exists.

preprint2016arXiv

Accurate $ab~initio$ tight-binding Hamiltonians: effective tools for electronic transport and optical spectroscopy from first principles

The calculations of electronic transport coefficients and optical properties require a very dense interpolation of the electronic band structure in reciprocal space that is computationally expensive and may have issues with band crossing and degeneracies. Capitalizing on a recently developed pseudo-atomic orbital projection technique, we exploit the exact tight-binding representation of the first principles electronic structure for the purposes of (1) providing an efficient strategy to explore the full band structure $E_n({\bf k})$, (2) computing the momentum operator differentiating directly the Hamiltonian, and (3) calculating the imaginary part of the dielectric function. This enables us to determine the Boltzmann transport coefficients and the optical properties within the independent particle approximation. In addition, the local nature of the tight-binding representation facilitates the calculation of the ballistic transport within the Landauer theory for systems with hundreds of atoms. In order to validate our approach we study the multi-valley band structure of CoSb$_3$ and a large core-shell nanowire using the ACBN0 functional. In CoSb$_3$ we point the many band minima contributing to the electronic transport that enhance the thermoelectric properties; for the core-shell nanowire we identify possible mechanisms for photo-current generation and justify the presence of protected transport channels in the wire.

preprint2016arXiv

Accurate Tight-Binding Hamiltonians for 2D and Layered Materials

We present a scheme to controllably improve the accuracy of tight-binding Hamiltonian matrices derived by projecting the solutions of plane-wave ab initio calculations on atomic orbital basis sets. By systematically increasing the completeness of the basis set of atomic orbitals, we are able to optimize the quality of the band structure interpolation over wide energy ranges including unoccupied states. This methodology is applied to the case of interlayer and image states, which appear several eV above the Fermi level in materials with large interstitial regions or surfaces such as graphite and graphene. Due to their spatial localization in the empty regions inside or outside of the system, these states have been inaccessible to traditional tight-binding models and even to ab initio calculations with atom-centered basis functions.

preprint2016arXiv

Addressing the lattice stability puzzle in the computational determination of intermetallic phase diagrams

The evaluation of phase stabilities of unstable elemental phases is a long-standing problem in the computational assessment of phase diagrams. Here we tackle this problem by explicitly calculating phase diagrams of intermetallic systems where its effect should be most conspicuous, binary systems of titanium with bcc transition metals (Mo, Nb, Ta and V). Two types of phase diagrams are constructed: one based on the lattice stabilities extracted from empirical data, and the other using the lattice stabilities computed from first principles. It is shown that the phase diagrams obtained using the empirical values contain clear contradictions with the experimental phase diagrams at the well known limits of low or high temperatures. Realistic phase diagrams, with a good agreement with the experimental observations, are achieved only when the computed lattice stability values are used. At intermediate temperatures, the computed phase diagrams resolve the controversy regarding the shape of the solvus in these systems, predicting a complex structure with a eutectoid transition and a miscibility gap between two bcc phases.

preprint2016arXiv

AFLUX: The LUX materials search API for the AFLOW data repositories

Automated computational materials science frameworks rapidly generate large quantities of materials data useful for accelerated materials design. We have extended the data oriented AFLOW-repository API (Application-Program-Interface, as described in Comput. Mater. Sci. 93, 178 (2014)) to enable programmatic access to search queries. A URI-based search API (Uniform Resource Identifier) is proposed for the construction of complex queries with the intent of allowing the remote creation and retrieval of customized data sets. It is expected that the new language AFLUX, acronym for Automatic Flow of LUX (light), will facilitate the creation of remote search operations on the AFLOW.org set of computational materials science data repositories.

preprint2016arXiv

High-Throughput Prediction of Finite-Temperature Properties using the Quasi-Harmonic Approximation

In order to calculate thermal properties in automatic fashion, the Quasi-Harmonic Approximation (QHA) has been combined with the Automatic Phonon Library (APL) and implemented within the AFLOW framework for high-throughput computational materials science. As a benchmark test to address the accuracy of the method and implementation, the specific heats, thermal expansion coefficients, Grüneisen parameters and bulk moduli have been calculated for 130 compounds. It is found that QHA-APL can reliably predict such values for several different classes of solids with root mean square relative deviation smaller than 28% with respect to experimental values. The automation, robustness, accuracy and precision of QHA-APL enable the computation of large material data sets, the implementation of repositories containing thermal properties, and finally can serve the community for data mining and machine learning studies.

preprint2016arXiv

Predicting Bulk Metallic Glass Forming Ability with the Thermodynamic Density of Competing Crystalline States

Despite two decades of studies, the formation of metallic glasses, very promising systems for industrial applications, still remains mostly unexplained. This lack of knowledge hinders the search for new systems, still performed with combinatorial trial and error. In the past, it was speculated that some sort of "confusion" during crystallization could play a key role during their formation. In this article, we propose a heuristic descriptor quantifying such confusion. It is based on the "thermodynamic density of competing crystalline states", parameterized from high-throughput ab-initio calculations. The existence of highly enthalpy-degenerate but geometrically-different phases frustrates the crystallization process and promotes glass formation. Two test beds are considered. A good and a bad glass-former, CuZr and NiZr, are experimentally characterized with high-throughput synthesis. The experimental results corroborate the capability of the heuristic descriptor in predicting glass forming ability through the compositional space. Our analysis is expected to deepen the understanding of the underlying mechanisms and to accelerate the discovery of novel metallic glasses.

preprint2016arXiv

Predicting the lattice thermal conductivity of solids by solving the Boltzmann transport equation: AFLOW - AAPL an automated, accurate and effcient framework

One of the most accurate approaches for calculating lattice thermal conductivity, $κ_l$, is solving the Boltzmann transport equation starting from third-order anharmonic force constants. In addition to the underlying approximations of ab-initio parameterization, two main challenges are associated with this path. High computational costs and lack of automation in the frameworks using this methodology affect the discovery rate of novel materials with ad-hoc properties. Here, we present the Automatic-Anharmonic-Phonon-Library, AAPL. It efficiently computes interatomic force constants by making effective use of crystal symmetry analysis, it solves the Boltzmann transport equation to obtain $κ_l$, and allows a fully integrated operation with minimum user intervention, a rational addition to the current high-throughput accelerated materials development framework AFLOW. We show an "experiment versus theory" study of the approach, we compare accuracy and speed with respect to other available packages, and for materials characterized by strong electron localization and correlation, we demonstrate that it is possible to improve accuracy without increasing computational requirements by combining AAPL with the pseudo-hybrid functional ACBN0.

preprint2016arXiv

Spectral descriptors for bulk metallic glasses based on the thermodynamics of competing crystalline phases

Metallic glasses have attracted considerable interest in recent years due to their unique combination of superb properties and processability. Predicting bulk metallic glass formers from known parameters remains a challenge and the search for new systems is still performed by trial and error. It has been speculated that some sort of "confusion" during crystallization of the crystalline phases competing with glass formation could play a key role. Here, we propose a heuristic descriptor quantifying confusion and demonstrate its validity by detailed experiments on two well-known glass forming alloy systems. With the insight provided by these results, we develop a robust model for predicting glass formation ability based on the spectral decomposition of geometrical and energetic features of crystalline phases calculated ab-initio in the AFLOW high throughput framework. Our findings indicate that the formation of metallic glass phases could be a much more common phenomenon than currently estimated, with more than 17% of binary alloy systems being potential glass formers. Our approach is capable of pinpointing favorable compositions, overcoming a major bottleneck hindering the discovery of new materials. Hence, it is demonstrated that smart descriptors, based solely on the energetics and structure of competing crystalline phases calculated from first-principles and available in online databases, others the sought-after key for accelerated discovery of novel metallic glasses.

preprint2016arXiv

Thermoelectric characterization of fine-grained Ti5O9 Magneli phase ceramics

Magneli phase Ti5O9 ceramics with 200-nm grain-size were fabricated by hot-pressing nanopowders of titanium and anatase TiO2 at 1223 K. The thermoelectric properties of these ceramics were investigated from room temperature to 1076 K. We show that the experimental variation of the electrical conductivity with temperature follows a small-polaron model and that the Seebeck coefficient can be explained by a temperature-dependent Heikes-Chaikin-Beni model. The thermoelectric figure-of-merit ZT of this nanoceramic material reaches 0.3 at 1076 K.

preprint2015arXiv

First principles thermodynamical modeling of the binodal and spinodal curves in lead chalcogenides

High-throughput ab-initio calculations, cluster expansion techniques and thermodynamic modeling have been synergistically combined to characterize the binodal and the spinodal decompositions features in the pseudo-binary lead chalcogenides PbSe-PbTe, PbS-PbTe, and PbS-PbSe. While our results agree with the available experimental data, our consolute temperatures substantially improve with respect to previous computational modeling. The computed phase diagrams corroborate that the formation of spinodal nanostructures causes low thermal conductivities in these alloys. The presented approach, making a rational use of online quantum repositories, can be extended to study thermodynamical and kinetic properties of materials of technological interest.

preprint2015arXiv

Improved predictions of the physical properties of Zn- and Cd-based wide band-gap semiconductors: a validation of the ACBN0 functional

We study the physical properties of Zn$X$ ($X$=O, S, Se, Te) and Cd$X$ ($X$=O, S, Se, Te) in the zinc-blende, rock-salt, and wurtzite structures using the recently developed fully $ab$ $initio$ pseudo-hybrid Hubbard density functional ACBN0. We find that both the electronic and vibrational properties of these wide-band gap semiconductors are systematically improved over the PBE values and reproduce closely the experimental measurements. Similar accuracy is found for the structural parameters, especially the bulk modulus. ACBN0 results compare well with hybrid functional calculations at a fraction of the computational cost.

preprint2015arXiv

Modeling Disordered Materials with a High Throughput ab-initio Approach

Predicting material properties of disordered systems remains a long-standing and formidable challenge in rational materials design. To address this issue, we introduce an automated software framework capable of modeling partial occupation within disordered materials using a high-throughput (HT) first principles approach. At the heart of the approach is the construction of supercells containing a virtually equivalent stoichiometry to the disordered material. All unique supercell permutations are enumerated and material properties of each are determined via HT electronic structure calculations. In accordance with a canonical ensemble of supercell states, the framework evaluates ensemble average properties of the system as a function of temperature. As proof of concept, we examine the framework's final calculated properties of a zinc chalcogenide (ZnS$_{1-x}$Se$_x$), a wide-gap oxide semiconductor (Mg$_{x}$Zn$_{1-x}$O), and an iron alloy (Fe$_{1-x}$Cu$_{x}$) at various stoichiometries.

preprint2015arXiv

The AFLOW Standard for High-Throughput Materials Science Calculations

The Automatic-Flow ( AFLOW ) standard for the high-throughput construction of materials science electronic structure databases is described. Electronic structure calculations of solid state materials depend on a large number of parameters which must be understood by researchers, and must be reported by originators to ensure reproducibility and enable collaborative database expansion. We therefore describe standard parameter values for k-point grid density, basis set plane wave kinetic energy cut-off, exchange-correlation functionals, pseudopotentials, DFT+U parameters, and convergence criteria used in AFLOW calculations.

preprint2014arXiv

A high-throughput ab initio review of platinum-group alloy systems

We report a comprehensive study of the binary systems of the platinum group metals with the transition metals, using high-throughput first-principles calculations. These computations predict stability of new compounds in 38 binary systems where no compounds have been reported in the literature experimentally, and a few dozen of as yet unreported compounds in additional systems. Our calculations also identify stable structures at compound compositions that have been previously reported without detailed structural data and indicate that some experimentally reported compounds may actually be unstable at low temperatures. With these results we construct enhanced structure maps for the binary alloys of platinum group metals. These are much more complete, systematic and predictive than those based on empirical results alone.

preprint2014arXiv

A RESTful API for exchanging Materials Data in the AFLOWLIB.org consortium

The continued advancement of science depends on shared and reproducible data. In the field of computational materials science and rational materials design this entails the construction of large open databases of materials properties. To this end, an Application Program Interface (API) following REST principles is introduced for the AFLOWLIB.org materials data repositories consortium. AUIDs (Aflowlib Unique IDentifier) and AURLs (Aflowlib Uniform Resource locator) are assigned to the database resources according to a well-defined protocol described herein, which enables the client to access, through appropriate queries, the desired data for post-processing. This introduces a new level of openness into the AFLOWLIB repository, allowing the community to construct high-level work-flows and tools exploiting its rich data set of calculated structural, thermodynamic, and electronic properties. Furthermore, federating these tools would open the door to collaborative investigation of the data by an unprecedented extended community of users to accelerate the advancement of computational materials design and development.

preprint2014arXiv

Composition-spread Growth and the Robust Topological Surface State of Kondo insulator SmB6 Thin Films

Topological insulators are a class of materials with insulating bulk but protected conducting surfaces due to the combination of spin-orbit interactions and time-reversal symmetry. The surface states are topologically non-trivial and robust against non-magnetic backscattering, leading to interesting physics and potential quantum computing applications1, 2. Recently there has been a fast growing interest in samarium hexboride (SmB6), a Kondo insulator predicted to be the first example of a correlated topological insulator3, 4. Here we fabricated smooth thin films of nanocrystalline SmB6 films. Their transport behavior indeed shows that SmB6 is a bulk insulator with topological surface states. Upon decreasing the temperature, the resistivity \r{ho} of Sm0.14B0.86 (SmB6) films display significant increase below 50 K due to hybridization gap formation, and it shows a saturation behavior below 10 K. The saturated resistance of our textured films is similar to that of the single crystals, suggesting that this conduction is from the surface and robust against grain boundary scatterings. Point contact spectroscopy (PCS) of the film using a superconducting tip displays both a Kondo Fano resonance and Andreev reflection, suggesting the existence of both an insulating Kondo lattice and metallic surface states.

preprint2014arXiv

Finding the stable structures of W$_x$N$_1$$_-$$_x$ with an $ab$$-$$initio$ high-throughput approach

Using density functional theory calculations, many researchers have predicted that various tungsten-nitride compounds WN$_x$ ($x$ > 1) will be "ultra-compressible" or "superhard", $i.e.$ as hard as or harder than diamond. These compounds are predicted to have large bulk and shear moduli (> 200 GPa) and to be elastically and vibrationally stable. Compounds with such desirable properties must be energetically stable against decomposition into other compounds. This stability can only be found after the determination of the convex hull for W$_x$N$_1$$_-$$_x$ lines which connect the lowest enthalpy structures as a function of composition. The phase diagram of the W-N structure is uncertain, both experimentally and computationally. Complex van der Waals forces play a significant role in determining the structure of solid N$_2$. Here we use high-throughput calculations to map out the convex hull and other low energy structures for the W-N system. We find that the ground state of the system is the NbO structure, and that the WN$_2$ structures found by Wang $et$ $al$. are also stable when the van der Waals forces are neglected. Other proposed structures are above the convex hull of the W-N system. We show how the choice of density functional influences the shape of the curve and the structures that form the hull. In nitrogen-rich compounds, the choice of functional can dramatically change the structural parameters and mechanical behavior. Using any of the functionals, the bulk and shear moduli of the NbO phase are comparable to the WN$_x$ compounds that have been claimed to be ultra-compressible for superhard.

preprint2014arXiv

Finding unprecedentedly low-thermal-conductivity half-Heusler semiconductors via high-throughput materials modeling

The lattice thermal conductivity (κω) is a key property for many potential applications of compounds. Discovery of materials with very low or high κω 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 κω 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 κω spans more than two orders of magnitude- a much larger range than that previously thought. κω 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 κω for this selected range of systems. Three semiconductors having κω < 5 W /(m K) are proposed for further experimental study.

preprint2014arXiv

High-Throughput Computational Screening of thermal conductivity, Debye temperature and Grüneisen parameter using a quasi-harmonic Debye Model

The quasi-harmonic Debye approximation has been implemented within the AFLOW and Materials Project frameworks for high-throughput computational science (Automatic Gibbs Library, AGL), in order to calculate thermal properties such as the Debye temperature and the thermal conductivity of materials. We demonstrate that the AGL method, which is significantly cheaper computationally compared to the fully ab initio approach, can reliably predict the ordinal ranking of the thermal conductivity for several different classes of semiconductor materials. We also find that for the set of 182 materials investigated in this work the Debye temperature, calculated with the AGL, is often a better predictor of the ordinal ranking of the experimental thermal conductivities than the calculated thermal conductivity. The Debye temperature is thus a potential descriptor for high-throughput screening of the thermal properties of materials.

preprint2014arXiv

Low thermal conductivity and triaxial phononic anisotropy of SnSe

In this theoretical study, we investigate the origins of the very low thermal conductivity of tin selenide (SnSe) using ab-initio calculations. We obtained high-temperature lattice thermal conductivity values that are close to those of amorphous compounds. We also found a strong anisotropy between the three crystallographic axes: one of the in-plane directions conducts heat much more easily than the other. Our results are compatible with most of the experimental literature on SnSe, and differ markedly from the more isotropic values reported by a recent study.

preprint2014arXiv

Materials Cartography: Representing and Mining Material Space Using Structural and Electronic Fingerprints

As the proliferation of high-throughput approaches in materials science is increasing the wealth of data in the field, the gap between accumulated-information and derived-knowledge widens. We address the issue of scientific discovery in materials databases by introducing novel analytical approaches based on structural and electronic materials fingerprints. The framework is employed to (i) query large databases of materials using similarity concepts, (ii) map the connectivity of the materials space (i.e., as a materials cartogram) for rapidly identifying regions with unique organizations/properties, and (iii) develop predictive Quantitative Materials Structure-Property Relation- ships (QMSPR) models for guiding materials design. In this study, we test these fingerprints by seeking target material properties. As a quantitative example, we model the critical temperatures of known superconductors. Our novel materials fingerprinting and materials cartography approaches contribute to the emerging field of materials informatics by enabling effective computational tools to analyze, visualize, model, and design new materials.

preprint2014arXiv

Pressure effects on the superconducting critical temperature of Li2B2

We present the structural, electronic and superconducting properties of Li2B2 under pressure within the framework of the density functional theory. The structural parameters, electronic band structure, phonon frequency of the E2g phonon mode and superconducting critical temperature Tc were calculated for pressures up to 20 GPa. We predicted that the superconducting critical temperature of Li2B2 is about 11 K and this decreases as pressure increases. We found that even though the lattice dynamics of the E2g phonon mode is similar to MgB2, the reduction of the σ-band density of states at Fermi level and the raising of the E2g phonon frequency with pressure were determinant to decrease λ and consequently Tc.

preprint2014arXiv

Reformulation of DFT+U as a pseudo-hybrid Hubbard density functional

The accurate prediction of the electronic properties of materials at a low computational expense is a necessary conditions for the development of effective high-throughput quantum-mechanics (HTQM) frameworks for accelerated materials discovery. HTQM infrastructures rely on the predictive capability of Density Functional Theory (DFT), the method of choice for the first principles study of materials properties. However, DFT suffers of approximations that result in a somewhat inaccurate description of the electronic band structure of semiconductors and insulators. In this article we introduce ACBN0, a pseudo-hybrid Hubbard density functional that yields an improved prediction of the band structure of insulators such as transition-metal oxides, as shown for TiO2, MnO, NiO and ZnO, with only a negligible increase in computational cost.

preprint2013arXiv

AFLOW: An automatic framework for high-throughput materials discovery

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 Aflow (Automatic Flow), a software framework for high-throughput calculation of crystal structure properties of alloys, intermetallics and inorganic compounds. The Aflow 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.

preprint2013arXiv

Effective and accurate representation of extended Bloch states on finite Hilbert spaces

We present a straightforward, noniterative projection scheme that can represent the electronic ground state of a periodic system on a finite atomic-orbital-like basis, up to a predictable number of electronic states and with controllable accuracy. By co-filtering the projections of plane-wave Bloch states with high-kinetic-energy components, the richness of the finite space and thus the number of exactly-reproduced bands can be selectively increased at a negligible computational cost, an essential requirement for the design of efficient algorithms for electronic structure simulations of realistic material systems and massive high-throughput investigations.

preprint2013arXiv

Stable ordered structures of binary technetium alloys from first principles

Technetium, element 43, is the only radioactive transition metal. It occurs naturally on earth in only trace amounts. Experimental investigation of its possible compounds is thus inherently difficult and limited. Half of the Tc-transition-metal systems (14 out of 28) are reported to be phase separating or lack experimental data. Using high-throughput first-principles calculations, we present a comprehensive investigation of the binary alloys of technetium with the transition metals. The calculations predict stable, ordered structures in nine of these 14 binary systems. They also predict additional compounds in all nine known compound-forming systems and in two of the five systems reported with disordered \c{hi} or σ phases. These results accentuate the incompleteness of the current knowledge of technetium alloys and should guide experimental studies to obtain accurate structural information on potential compounds, obviating some of the difficulties associated with such work.

preprint2010arXiv

High-throughput electronic band structure calculations: challenges and tools

The article is devoted to the discussion of the high-throughput approach to band structures calculations. We present scientific and computational challenges as well as solutions relying on the developed framework (Automatic Flow, AFLOW/ACONVASP). The key factors of the method are the standardization and the robustness of the procedures. Two scenarios are relevant: 1) independent users generating databases in their own computational systems (off-line approach) and 2) teamed users sharing computational information based on a common ground (on-line approach). Both cases are integrated in the framework: for off-line approaches, the standardization is automatic and fully integrated for the 14 Bravais lattices, the primitive and conventional unit cells, and the coordinates of the high symmetry k-path in the Brillouin zones. For on-line tasks, the framework offers an expandable web interface where the user can prepare and set up calculations following the proposed standard. Few examples of band structures are included. LSDA+U parameters (U, J) are also presented for Nd, Sm, and Eu.

preprint2010arXiv

Structure maps for hcp metals from first principles calculations

The ability to predict the existence and crystal type of ordered structures of materials from their components is a major challenge of current materials research. Empirical methods use experimental data to construct structure maps and make predictions based on clustering of simple physical parameters. Their usefulness depends on the availability of reliable data over the entire parameter space. Recent development of high throughput methods opens the possibility to enhance these empirical structure maps by {\it ab initio} calculations in regions of the parameter space where the experimental evidence is lacking or not well characterized. In this paper we construct enhanced maps for the binary alloys of hcp metals, where the experimental data leaves large regions of poorly characterized systems believed to be phase-separating. In these enhanced maps, the clusters of non-compound forming systems are much smaller than indicated by the empirical results alone.

preprint2009arXiv

Calculation of solubility in titanium alloys from first-principles

We present an approach to calculate the atomic bulk solubility in binary alloys based on the statistical-thermodynamic theory of dilute lattice gas. The model considers all the appropriate ground states of the alloy and results in a simple Arrhenius-type temperature dependence determined by a {\it "low-solubility formation enthalpy"}. This quantity, directly obtainable from first-principle calculations, is defined as the composition derivative of the compound formation enthalpy with respect to nearby ground states. We apply the framework and calculate the solubility of the A specie in A-Ti alloys (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). In addition to determining unknown low-temperature ground states for the eight alloys, we find qualitative agreements with solubility experimental results. The presented formalism, correct in the low-solubility limit, should be considered as an appropriate starting point for determining if more computationally expensive formalisms are otherwise needed.

preprint2009arXiv

Comparative Study of Nonproportionality and Electronic Band Structures Features in Scintillator Materials

The origin of nonproportionality in scintillator materials has been a long standing problem for more than four decades. In this manuscript, we show that, with the help of first principle modeling, the parameterization of the nonproportionality for several systems, with respect to their band structure curvature suggests a correlation between carrier effective mass and energy response. We attribute this correlation to the case where free electrons and holes are the major energy carriers. Excitonic scintillators do not show such a definitive trend. This model suggests a potential high-throughput approach for discovering novel proportional scintillators in the former class of materials.

preprint2008arXiv

Structures and topological transitions of hydrocarbon films on quasicrystalline surfaces

Lubricants can affect quasicrystalline coatings surfaces by modifying commensurability of the interfaces. We report results of the first computer simulation studies of physically adsorbed hydrocarbons on a quasicrystalline surface: methane, propane, and benzene on decagonal Al-Ni-Co. The grand canonical Monte Carlo method is employed, using novel Embedded Atom Method potentials generated from it ab initio calculations, and standard hydrocarbon interactions. The resulting adsorption isotherms and calculated structures show the films' evolution from submonolayer to condensation. We discover the presence and absence of the 5- to 6-fold topological transition, for benzene and methane, respectively, in agreement with a previsouly formulated phenomenological rule based on adsorbate-substrate size mismatch.