-
Approaches for Uncertainty Quantification of AI-predicted Material Properties: A Comparison
Authors:
Francesca Tavazza,
Kamal Choudhary,
Brian DeCost
Abstract:
The development of large databases of material properties, together with the availability of powerful computers, has allowed machine learning (ML) modeling to become a widely used tool for predicting material performances. While confidence intervals are commonly reported for such ML models, prediction intervals, i.e., the uncertainty on each prediction, are not as frequently available. Here, we in…
▽ More
The development of large databases of material properties, together with the availability of powerful computers, has allowed machine learning (ML) modeling to become a widely used tool for predicting material performances. While confidence intervals are commonly reported for such ML models, prediction intervals, i.e., the uncertainty on each prediction, are not as frequently available. Here, we investigate three easy-to-implement approaches to determine such individual uncertainty, comparing them across ten ML quantities spanning energetics, mechanical, electronic, optical, and spectral properties. Specifically, we focused on the Quantile approach, the direct machine learning of the prediction intervals and Ensemble methods.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
JARVIS-Leaderboard: A Large Scale Benchmark of Materials Design Methods
Authors:
Kamal Choudhary,
Daniel Wines,
Kangming Li,
Kevin F. Garrity,
Vishu Gupta,
Aldo H. Romero,
Jaron T. Krogel,
Kayahan Saritas,
Addis Fuhr,
Panchapakesan Ganesh,
Paul R. C. Kent,
Keqiang Yan,
Yuchao Lin,
Shuiwang Ji,
Ben Blaiszik,
Patrick Reiser,
Pascal Friederich,
Ankit Agrawal,
Pratyush Tiwary,
Eric Beyerle,
Peter Minch,
Trevor David Rhone,
Ichiro Takeuchi,
Robert B. Wexler,
Arun Mannodi-Kanakkithodi
, et al. (13 additional authors not shown)
Abstract:
Lack of rigorous reproducibility and validation are major hurdles for scientific development across many fields. Materials science in particular encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform…
▽ More
Lack of rigorous reproducibility and validation are major hurdles for scientific development across many fields. Materials science in particular encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform with multiple data modalities with both perfect and defect materials data is still lacking. This work introduces JARVIS-Leaderboard, an open-source and community-driven platform that facilitates benchmarking and enhances reproducibility. The platform allows users to set up benchmarks with custom tasks and enables contributions in the form of dataset, code, and meta-data submissions. We cover the following materials design categories: Artificial Intelligence (AI), Electronic Structure (ES), Force-fields (FF), Quantum Computation (QC) and Experiments (EXP). For AI, we cover several types of input data, including atomic structures, atomistic images, spectra, and text. For ES, we consider multiple ES approaches, software packages, pseudopotentials, materials, and properties, comparing results to experiment. For FF, we compare multiple approaches for material property predictions. For QC, we benchmark Hamiltonian simulations using various quantum algorithms and circuits. Finally, for experiments, we use the inter-laboratory approach to establish benchmarks. There are 1281 contributions to 274 benchmarks using 152 methods with more than 8 million data-points, and the leaderboard is continuously expanding. The JARVIS-Leaderboard is available at the website: https://pages.nist.gov/jarvis_leaderboard
△ Less
Submitted 26 March, 2024; v1 submitted 20 June, 2023;
originally announced June 2023.
-
Recent progress in the JARVIS infrastructure for next-generation data-driven materials design
Authors:
Daniel Wines,
Ramya Gurunathan,
Kevin F. Garrity,
Brian DeCost,
Adam J. Biacchi,
Francesca Tavazza,
Kamal Choudhary
Abstract:
The Joint Automated Repository for Various Integrated Simulations (JARVIS) infrastructure at the National Institute of Standards and Technology (NIST) is a large-scale collection of curated datasets and tools with more than 80000 materials and millions of properties. JARVIS uses a combination of electronic structure, artificial intelligence (AI), advanced computation and experimental methods to ac…
▽ More
The Joint Automated Repository for Various Integrated Simulations (JARVIS) infrastructure at the National Institute of Standards and Technology (NIST) is a large-scale collection of curated datasets and tools with more than 80000 materials and millions of properties. JARVIS uses a combination of electronic structure, artificial intelligence (AI), advanced computation and experimental methods to accelerate materials design. Here we report some of the new features that were recently included in the infrastructure such as: 1) doubling the number of materials in the database since its first release, 2) including more accurate electronic structure methods such as Quantum Monte Carlo, 3) including graph neural network-based materials design, 4) development of unified force-field, 5) development of a universal tight-binding model, 6) addition of computer-vision tools for advanced microscopy applications, 7) development of a natural language processing tool for text-generation and analysis, 8) debuting a large-scale benchmarking endeavor, 9) including quantum computing algorithms for solids, 10) integrating several experimental datasets and 11) staging several community engagement and outreach events. New classes of materials, properties, and workflows added to the database include superconductors, two-dimensional (2D) magnets, magnetic topological materials, metal-organic frameworks, defects, and interface systems. The rich and reliable datasets, tools, documentation, and tutorials make JARVIS a unique platform for modern materials design. JARVIS ensures openness of data and tools to enhance reproducibility and transparency and to promote a healthy and collaborative scientific environment.
△ Less
Submitted 25 October, 2023; v1 submitted 19 May, 2023;
originally announced May 2023.
-
High-Throughput DFT-Based Discovery of Next Generation Two-Dimensional (2D) Superconductors
Authors:
Daniel Wines,
Kamal Choudhary,
Adam J. Biacchi,
Kevin F. Garrity,
Francesca Tavazza
Abstract:
High-throughput density functional theory (DFT) calculations allow for a systematic search for conventional superconductors. With the recent interest in two-dimensional (2D) superconductors, we used a high-throughput workflow to screen over 1,000 2D materials in the JARVIS-DFT database and performed electron-phonon coupling calculations, using the McMillan-Allen-Dynes formula to calculate the supe…
▽ More
High-throughput density functional theory (DFT) calculations allow for a systematic search for conventional superconductors. With the recent interest in two-dimensional (2D) superconductors, we used a high-throughput workflow to screen over 1,000 2D materials in the JARVIS-DFT database and performed electron-phonon coupling calculations, using the McMillan-Allen-Dynes formula to calculate the superconducting transition temperature ($T_c$) for 165 of them. Of these 165 materials, we identify 34 dynamically stable structures with transition temperatures above 5 K, including materials such as W$_2$N$_3$, NbO$_2$, ZrBrO, TiClO, NaSn$_2$S$_4$, Mg$_2$B$_4$C$_2$ and the previously unreported Mg$_2$B$_4$N$_2$ ($T_c$ = 21.8 K). Finally, we performed experiments to determine the $T_c$ of selected layered superconductors (2H-NbSe$_2$, 2H-NbS$_2$, ZrSiS, FeSe) and discuss the measured results within the context of our DFT results. We aim that the outcome of this workflow can guide future computational and experimental studies of new and emerging 2D superconductors by providing a roadmap of high-throughput DFT data.
△ Less
Submitted 10 February, 2023; v1 submitted 9 November, 2022;
originally announced November 2022.
-
Systematic DFT+U and Quantum Monte Carlo benchmark of magnetic two-dimensional (2D) CrX$_3$ (X = I, Br, Cl, F)
Authors:
Daniel Wines,
Kamal Choudhary,
Francesca Tavazza
Abstract:
The search for two-dimensional (2D) magnetic materials has attracted a great deal of attention because of the experimental synthesis of 2D CrI$_3$, which has a measured Curie temperature of 45 K. Often times, these monolayers have a higher degree of electron correlation and require more sophisticated methods beyond density functional theory (DFT). Diffusion Monte Carlo (DMC) is a correlated electr…
▽ More
The search for two-dimensional (2D) magnetic materials has attracted a great deal of attention because of the experimental synthesis of 2D CrI$_3$, which has a measured Curie temperature of 45 K. Often times, these monolayers have a higher degree of electron correlation and require more sophisticated methods beyond density functional theory (DFT). Diffusion Monte Carlo (DMC) is a correlated electronic structure method that has been demonstrated successful for calculating the electronic and magnetic properties of a wide variety of 2D and bulk systems, since it has a weaker dependence on the Hubbard parameter (U) and density functional. In this study we designed a workflow that combines DFT+U and DMC in order to treat 2D correlated magnetic systems. We chose monolayer CrX$_3$ (X = I, Br, Cl, F), with a stronger focus on CrI$_3$ and CrBr$_3$, as a case study due to the fact that they have been experimentally realized and have a finite critical temperature. With this DFT+U and DMC workflow and the analytical method of Torelli and Olsen, we estimated an upper bound of 43.56 K for the T$_c$ of CrI$_3$ and 20.78 K for the T$_c$ of CrBr$_3$, in addition to analyzing the spin densities and magnetic properties with DMC and DFT+U. We expect that running this workflow for a well-known material class will aid in the future discovery and characterization of lesser known and more complex correlated 2D magnetic materials.
△ Less
Submitted 13 January, 2023; v1 submitted 21 September, 2022;
originally announced September 2022.
-
Unified Graph Neural Network Force-field for the Periodic Table
Authors:
Kamal Choudhary,
Brian DeCost,
Lily Major,
Keith Butler,
Jeyan Thiyagalingam,
Francesca Tavazza
Abstract:
Classical force fields (FF) based on machine learning (ML) methods show great potential for large scale simulations of materials. MLFFs have hitherto largely been designed and fitted for specific systems and are not usually transferable to chemistries beyond the specific training set. We develop a unified atomisitic line graph neural network-based FF (ALIGNN-FF) that can model both structurally an…
▽ More
Classical force fields (FF) based on machine learning (ML) methods show great potential for large scale simulations of materials. MLFFs have hitherto largely been designed and fitted for specific systems and are not usually transferable to chemistries beyond the specific training set. We develop a unified atomisitic line graph neural network-based FF (ALIGNN-FF) that can model both structurally and chemically diverse materials with any combination of 89 elements from the periodic table. To train the ALIGNN-FF model, we use the JARVIS-DFT dataset which contains around 75000 materials and 4 million energy-force entries, out of which 307113 are used in the training. We demonstrate the applicability of this method for fast optimization of atomic structures in the crystallography open database and by predicting accurate crystal structures using genetic algorithm for alloys.
△ Less
Submitted 16 September, 2022; v1 submitted 12 September, 2022;
originally announced September 2022.
-
Rapid Prediction of Phonon Structure and Properties using an Atomistic Line Graph Neural Network (ALIGNN)
Authors:
Ramya Gurunathan,
Kamal Choudhary,
Francesca Tavazza
Abstract:
The phonon density-of-states (DOS) summarizes the lattice vibrational modes supported by a structure, and gives access to rich information about the material's stability, thermodynamic constants, and thermal transport coefficients. Here, we present an atomistic line graph neural network (ALIGNN) model for the prediction of the phonon density of states and the derived thermal and thermodynamic prop…
▽ More
The phonon density-of-states (DOS) summarizes the lattice vibrational modes supported by a structure, and gives access to rich information about the material's stability, thermodynamic constants, and thermal transport coefficients. Here, we present an atomistic line graph neural network (ALIGNN) model for the prediction of the phonon density of states and the derived thermal and thermodynamic properties. The model is trained on a database of over 14,000 phonon spectra included in the JARVIS-DFT (Joint Automated Repository for Various Integrated Simulations: Density Functional Theory) database. The model predictions are shown to capture the spectral features of the phonon density-of-states, effectively categorize dynamical stability, and lead to accurate predictions of DOS-derived thermal and thermodynamic properties, including heat capacity $C_{\mathrm{V}}$, vibrational entropy $S_{\mathrm{vib}}$, and the isotopic phonon scattering rate $τ^{-1}_{\mathrm{i}}$. The DOS-mediated ALIGNN model provides superior predictions when compared to a direct deep-learning prediction of these material properties as well as predictions based on analytic simplifications of the phonon DOS, including the Debye or Born-von Karman models. Finally, the ALIGNN model is used to predict the phonon spectra and properties for about 40,000 additional materials listed in the JARVIS-DFT database, which are validated as far as possible against other open-sourced high-throughput DFT phonon databases.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Topological Phases of Mn$A_{2}$$X_{4}$ (A=Bi, Sb; X = Se, Te) under Magnetic Field
Authors:
Sugata Chowdhury,
Kevin F. Garrity,
Francesca Tavazza
Abstract:
The concept of electronic topology and the associated topological protection brings excellent opportunities for developing next-generation devices. Ideally, magnetic topological materials (MTM) should have their Dirac/Weyl points and/or associated mass gaps at the Fermi energy (EF) or be readily tunable such that they can be placed at EF via external perturbations such as electric field gating, ch…
▽ More
The concept of electronic topology and the associated topological protection brings excellent opportunities for developing next-generation devices. Ideally, magnetic topological materials (MTM) should have their Dirac/Weyl points and/or associated mass gaps at the Fermi energy (EF) or be readily tunable such that they can be placed at EF via external perturbations such as electric field gating, chemical substitutions, or doping. Three-dimensional antiferromagnetic (AFM) materials like Mn$Bi_{2}$$X_{4}$ (X=Se, Te) that have strong spin-orbit coupling (SOC) and broken time-reversal symmetry (TRS) due to magnetic ordering have been the subject of enormous interest.. In this work, using density functional theory (DFT), we have studied the electronic properties and topological phases of the first intrinsic magnetic topological insulator family Mn$Bi_{2}$$Te_{4}$ (MBT) in the presence of an external magnetic field. Our calculations reveal that the topological phase of bulk rhombohedral (R$\overline{3}$m) Mn$A_{2}$$X_{4}$ (A = Sb, Bi; X=Se, Te) depends on the spin direction and the chemistry. The antiferromagnetic (AFM) ground state of Mn$Sb_{2}$$Se_{4}$ (MSS) is a trivial insulator, whereas the AFM ground state of Mn$Bi_{2}$$Se_{4}$(MBS) is an Axion insulator. Both materials become nodal point or nodal line Weyl semimetals in the presence of a sufficiently strong external magnetic field. The AFM ground state of Mn$Sb_{2}$$Te_{4}$ (MST) is an Axion insulator. MST is a type-II Weyl semimetal with spins aligned in the $\hat{Z}$-direction, but becomes insulating with an inverted band gap for spins in-plane. Similarly, the AFM phase of MBT is an Axion insulator, but remains insulating with an inverted gap in the ferromagnetic phase. Additionally, we demonstrated the evolution of the topological phase of Mn$Bi_{2}$$Te_{4}$ (MBT) by substituting the Bi atoms with the Sb atoms.
△ Less
Submitted 1 June, 2022; v1 submitted 14 February, 2022;
originally announced February 2022.
-
Review of Theoretical and Computational Methods for 2D Materials Exhibiting Charge Density Waves
Authors:
Sugata Chowdhury,
Heather M. Hill,
Albert F. Rigosi,
Patrick M. Vora,
Angela R. Hight Walker,
Francesca Tavazza
Abstract:
Two-dimensional (2D) materials that exhibit charge density waves (CDWs) have generated many research endeavors in the hopes of employing their exotic properties for various quantum-based technologies. Early investigations surrounding CDWs were mostly focused on bulk materials. However, applications for quantum devices have since required devices to be constructed from few-layer material to fully u…
▽ More
Two-dimensional (2D) materials that exhibit charge density waves (CDWs) have generated many research endeavors in the hopes of employing their exotic properties for various quantum-based technologies. Early investigations surrounding CDWs were mostly focused on bulk materials. However, applications for quantum devices have since required devices to be constructed from few-layer material to fully utilize the material's properties. This field has greatly expanded over the decades, warranting a focus on the computational efforts surrounding CDWs in 2D materials. In this review, we will cover ground in the following relevant, theory-driven subtopics for TaS2 and TaSe2: summary of general computational techniques and methods, atomic structures, Raman modes, and effects of confinement and dimensionality. Through understanding how the computational methods have enabled incredible advancements in quantum materials, one may anticipate the ever-expanding directions available for continued pursuit as the field brings us through the 21st century.
△ Less
Submitted 17 January, 2022;
originally announced January 2022.
-
Bonding characteristics of the interfacial buffer layer in epitaxial graphene via density functional theory
Authors:
Alana Okullo,
Heather M. Hill,
Albert F. Rigosi,
Angela R. Hight Walker,
Francesca Tavazza,
Sugata Chowdhury
Abstract:
Monolayer epitaxial graphene is an appropriate candidate for a wide variety of electronic and optical applications. One advantage of growing graphene on the Si face of SiC is that it develops as a single crystal, as does the layer underneath, commonly referred to as the interfacial buffer layer. The properties of this supporting layer include a band gap, making it of interest to groups seeking to…
▽ More
Monolayer epitaxial graphene is an appropriate candidate for a wide variety of electronic and optical applications. One advantage of growing graphene on the Si face of SiC is that it develops as a single crystal, as does the layer underneath, commonly referred to as the interfacial buffer layer. The properties of this supporting layer include a band gap, making it of interest to groups seeking to build devices with on-off capabilities. In this work, using density functional theory, we have calculated the bonding characteristics of the buffer layer to the SiC substrate beneath. These calculations were used to determine a periodic length between the covalent bonds acting as anchor points in this interface. Additionally, it is evident that the formation of these anchor points depends on the lattice mismatch between the graphene layer and SiC.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
Recent Advances and Applications of Deep Learning Methods in Materials Science
Authors:
Kamal Choudhary,
Brian DeCost,
Chi Chen,
Anubhav Jain,
Francesca Tavazza,
Ryan Cohn,
Cheol WooPark,
Alok Choudhary,
Ankit Agrawal,
Simon J. L. Billinge,
Elizabeth Holm,
Shyue Ping Ong,
Chris Wolverton
Abstract:
Deep learning (DL) is one of the fastest growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. Recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular.…
▽ More
Deep learning (DL) is one of the fastest growing topics in materials data science, with rapidly emerging applications spanning atomistic, image-based, spectral, and textual data modalities. DL allows analysis of unstructured data and automated identification of features. Recent development of large materials databases has fueled the application of DL methods in atomistic prediction in particular. In contrast, advances in image and spectral data have largely leveraged synthetic data enabled by high quality forward models as well as by generative unsupervised DL methods. In this article, we present a high-level overview of deep-learning methods followed by a detailed discussion of recent developments of deep learning in atomistic simulation, materials imaging, spectral analysis, and natural language processing. For each modality we discuss applications involving both theoretical and experimental data, typical modeling approaches with their strengths and limitations, and relevant publicly available software and datasets. We conclude the review with a discussion of recent cross-cutting work related to uncertainty quantification in this field and a brief perspective on limitations, challenges, and potential growth areas for DL methods in materials science. The application of DL methods in materials science presents an exciting avenue for future materials discovery and design.
△ Less
Submitted 27 October, 2021;
originally announced October 2021.
-
Uncertainty Prediction for Machine Learning Models of Material Properties
Authors:
Francesca Tavazza,
Brian De Cost,
Kamal Choudhary
Abstract:
Uncertainty quantification in Artificial Intelligence (AI)-based predictions of material properties is of immense importance for the success and reliability of AI applications in material science. While confidence intervals are commonly reported for machine learning (ML) models, prediction intervals, i.e., the evaluation of the uncertainty on each prediction, are seldomly available. In this work w…
▽ More
Uncertainty quantification in Artificial Intelligence (AI)-based predictions of material properties is of immense importance for the success and reliability of AI applications in material science. While confidence intervals are commonly reported for machine learning (ML) models, prediction intervals, i.e., the evaluation of the uncertainty on each prediction, are seldomly available. In this work we compare 3 different approaches to obtain such individual uncertainty, testing them on 12 ML-physical properties. Specifically, we investigated using the Quantile loss function, machine learning the prediction intervals directly and using Gaussian Processes. We identify each approachs advantages and disadvantages and end up slightly favoring the modeling of the individual uncertainties directly, as it is the easiest to fit and, in most cases, minimizes over-and under-estimation of the predicted errors. All data for training and testing were taken from the publicly available JARVIS-DFT database, and the codes developed for computing the prediction intervals are available through JARVIS-Tools.
△ Less
Submitted 16 July, 2021;
originally announced July 2021.
-
High-throughput search for magnetic topological materials using spin-orbit spillage, machine-learning and experiments
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Nirmal J. Ghimire,
Naween Anand,
Francesca Tavazza
Abstract:
Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation, but very few high-quality candidate materials are known. In this work, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from 40000 three-dimensional materials in the…
▽ More
Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation, but very few high-quality candidate materials are known. In this work, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from 40000 three-dimensional materials in the JARVIS-DFT database (https://jarvis.nist.gov/jarvisdft). First, we screen materials with net magnetic moment > 0.5 μB and spin-orbit spillage > 0.25, resulting in 25 insulating and 564 metallic candidates. The spillage acts as a signature of spin-orbit induced band-inversion. Then, we carry out calculations of Wannier charge centers, Chern numbers, anomalous Hall conductivities, surface bandstructures, and Fermi-surfaces to determine interesting topological characteristics of the screened compounds. We also train machine learning models for predicting the spillage, bandgaps, and magnetic moments of new compounds, to further accelerate the screening process. We experimentally synthesize and characterize a few candidate materials to support our theoretical predictions.
△ Less
Submitted 5 February, 2021; v1 submitted 30 January, 2021;
originally announced February 2021.
-
Examining Experimental Raman Mode Behavior in Mono- and Bi-layer 2H-TaSe$_{2}$ via Density Functional Theory
Authors:
Sugata Chowdhury,
Heather M. Hill,
Albert F. Rigosi,
Andrew Briggs,
Helmuth Berger,
David B. Newell,
Angela R. Hight Walker,
Francesca Tavazza
Abstract:
Tantalum diselenide (TaSe$_{2}$) is a metallic transition metal dichalcogenide whose equilibrium structure and vibrational behavior strongly depends on temperature and thickness, including the emergence of charge density wave (CDW) states at very low T. In this work, observed modes for mono- and bi-layer are described across several spectral regions and com-pared to the bulk ones. Such modes, incl…
▽ More
Tantalum diselenide (TaSe$_{2}$) is a metallic transition metal dichalcogenide whose equilibrium structure and vibrational behavior strongly depends on temperature and thickness, including the emergence of charge density wave (CDW) states at very low T. In this work, observed modes for mono- and bi-layer are described across several spectral regions and com-pared to the bulk ones. Such modes, including an experimentally observed forbidden Raman mode and low frequency CDW modes, are then matched to corresponding density functional theory (DFT) predicted vibrations, to unveil their inner working. The excellent match between experimental and computational results justifies the presented vibrational visualizations of these modes. Additional support is provided by experimental phonons seen in Raman spectra as a function of temperature and thickness. These results highlight the importance of understanding interlayer interactions and their effects on mode behaviors.
△ Less
Submitted 5 December, 2020;
originally announced December 2020.
-
Predicting Anomalous Quantum Confinement Effect in van der Waals Materials
Authors:
Kamal Choudhary,
Francesca Tavazza
Abstract:
Materials with van der Waals-bonding are known to exhibit quantum confinement effect, in which the electronic bandgap of the three-dimensional (3D) realization of a material is lower than that of its two-dimensional (2D) counterpart. However, the possibility of an anomalous quantum confinement effect (AQCE) exists, where the bandgap trend is reversed. In this work, we computationally identify mate…
▽ More
Materials with van der Waals-bonding are known to exhibit quantum confinement effect, in which the electronic bandgap of the three-dimensional (3D) realization of a material is lower than that of its two-dimensional (2D) counterpart. However, the possibility of an anomalous quantum confinement effect (AQCE) exists, where the bandgap trend is reversed. In this work, we computationally identify materials for which such AQCE occurs. Using density functional theory (DFT), we compute ~1000 OptB88vdW (semi-local functional), ~50 HSE06 and ~50 PBE0 (hybrid functional) bandgaps for bulk and their corresponding monolayers in the JARVIS-DFT database. OptB88vdW identifies 65 AQCE materials, but the hybrid functionals only confirm such finding in 14 cases. Some of the AQCE systems identified through HSE06 and PBE0 are: hydroxides or oxide hydroxide compounds (AlOH2, Mg(OH)2, Mg2H2O3, Ni(OH)2, SrH2O3) as well as Sb-halogen-chalcogenide compounds (SbSBr, SbSeI) and alkali-chalcogenides (RbLiS and RbLiSe). A detailed electronic structure analysis, based on band-structure and projected density of states, shows AQCE is often characterized by lowering of the conduction band in the monolayer and corresponding changes in the pz electronic orbital contribution, with z being the non-periodic direction in the 2D case. We believe our computational results would spur the effort to validate the results experimentally and will have impact on bandgap engineering applications based on low-dimensional materials.
△ Less
Submitted 12 December, 2020; v1 submitted 19 November, 2020;
originally announced November 2020.
-
Uncertainty Quantification in Atomistic Modeling of Metals and its Effect on Mesoscale and Continuum Modeling A Review
Authors:
Joshua J. Gabriel,
Noah H. Paulson,
Thien C. Duong,
Francesca Tavazza,
Chandler A. Becker,
Santanu Chaudhuri,
Marius Stan
Abstract:
The design of next-generation alloys through the Integrated Computational Materials Engineering (ICME) approach relies on multi-scale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as Density Functional Theory (DFT) and Molecular Dynamics (MD) have been successful in predicting properties of never before studied compounds…
▽ More
The design of next-generation alloys through the Integrated Computational Materials Engineering (ICME) approach relies on multi-scale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as Density Functional Theory (DFT) and Molecular Dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies have emerged that mitigate this gap. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as Phase Field Method (PFM) and Calculation of Phase Diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist and Bayesian methodologies. Challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.
△ Less
Submitted 28 September, 2020;
originally announced September 2020.
-
The Joint Automated Repository for Various Integrated Simulations (JARVIS) for data-driven materials design
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Andrew C. E. Reid,
Brian DeCost,
Adam J. Biacchi,
Angela R. Hight Walker,
Zachary Trautt,
Jason Hattrick-Simpers,
A. Gilad Kusne,
Andrea Centrone,
Albert Davydov,
Jie Jiang,
Ruth Pachter,
Gowoon Cheon,
Evan Reed,
Ankit Agrawal,
Xiaofeng Qian,
Vinit Sharma,
Houlong Zhuang,
Sergei V. Kalinin,
Bobby G. Sumpter,
Ghanshyam Pilania,
Pinar Acar,
Subhasish Mandal,
Kristjan Haule
, et al. (3 additional authors not shown)
Abstract:
The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and d…
▽ More
The Joint Automated Repository for Various Integrated Simulations (JARVIS) is an integrated infrastructure to accelerate materials discovery and design using density functional theory (DFT), classical force-fields (FF), and machine learning (ML) techniques. JARVIS is motivated by the Materials Genome Initiative (MGI) principles of developing open-access databases and tools to reduce the cost and development time of materials discovery, optimization, and deployment. The major features of JARVIS are: JARVIS-DFT, JARVIS-FF, JARVIS-ML, and JARVIS-Tools. To date, JARVIS consists of 40,000 materials and 1 million calculated properties in JARVIS-DFT, 1,500 materials and 110 force-fields in JARVIS-FF, and 25 ML models for material-property predictions in JARVIS-ML, all of which are continuously expanding. JARVIS-Tools provides scripts and workflows for running and analyzing various simulations. We compare our computational data to experiments or high-fidelity computational methods wherever applicable to evaluate error/uncertainty in predictions. In addition to the existing workflows, the infrastructure can support a wide variety of other technologically important applications as part of the data-driven materials design paradigm. The JARVIS datasets and tools are publicly available at the website: https://jarvis.nist.gov .
△ Less
Submitted 11 July, 2021; v1 submitted 3 July, 2020;
originally announced July 2020.
-
Influence of Dimensionality on the Charge Density Wave Phase of 2H-TaSe$_{2}$
Authors:
Sugata Chowdhury,
Albert F. Rigosi,
Heather M. Hill,
Andrew Briggs,
David B. Newell,
Helmuth Berger,
Angela R. Hight Walker,
Francesca Tavazza
Abstract:
Metallic transition metal dichalcogenides like tantalum diselenide (TaSe$_{2}$) exhibit exciting behaviors at low temperatures, including the emergence of charge density wave (CDW) states. In this work, density functional theory (DFT) is used to investigate how structural, electronic, and Raman spectral properties of the CDW configuration change as a function of thickness. Such findings highlight…
▽ More
Metallic transition metal dichalcogenides like tantalum diselenide (TaSe$_{2}$) exhibit exciting behaviors at low temperatures, including the emergence of charge density wave (CDW) states. In this work, density functional theory (DFT) is used to investigate how structural, electronic, and Raman spectral properties of the CDW configuration change as a function of thickness. Such findings highlight the influence of dimensionality change (from 2D to 3D) and van der Waals (vdW) interactions on the system properties. The vdW effect is most strongly present in bulk TaSe$_{2}$ in the spectral range 165 cm$^{-1}$ to 215 cm$^{-1}$. The phonons seen in the experimental Raman spectra are compared with the results calculated from the DFT models as a function of temperature and layer number. The matching of data and calculations substantiates the model's description of the CDW structural formation as a function of thickness, which is shown in depth for 1L through 6L systems. These results highlight the importance of understanding interlayer interactions, which are pervasive in many quantum phenomena involving two-dimensional confinement.
△ Less
Submitted 20 August, 2021; v1 submitted 22 May, 2020;
originally announced May 2020.
-
Efficient Computational Design of 2D van der Waals Heterostructures: Band-Alignment, Lattice-Mismatch, Web-app Generation and Machine-learning
Authors:
Kamal Choudhary,
Kevin F Garrity,
Steven T. Hartman,
Ghanshyam Pilania,
Francesca Tavazza
Abstract:
We develop a computational database, web-apps, and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D)-heterostructures. Using density functional theory (DFT) based lattice-parameters and electronic band-energies for 674 non-metallic exfoliable 2D-materials, we generate 226,779 possible heterostructures. We classify these heterostructures into type-I, II and…
▽ More
We develop a computational database, web-apps, and machine-learning (ML) models to accelerate the design and discovery of two-dimensional (2D)-heterostructures. Using density functional theory (DFT) based lattice-parameters and electronic band-energies for 674 non-metallic exfoliable 2D-materials, we generate 226,779 possible heterostructures. We classify these heterostructures into type-I, II and III systems according to Andersons rule, which is based on the relative band-alignments of the non-interacting monolayers. We find that type-II is the most common and the type-III the least common heterostructure type. We subsequently analyze the chemical trends for each heterostructure type in terms of the periodic table of constituent elements. The band alignment data can be also used for identifying photocatalysts and high-work function 2D-metals for contacts. We validate our results by comparing them to experimental data as well as hybrid-functional predictions. Additionally, we carry out DFT calculations of a few selected systems (MoS2/WSe2, MoS2/h-BN, MoSe2/CrI3), to compare the band-alignment description with the predictions from Andersons rule. We develop web-apps to enable users to virtually create combinations of 2D materials and predict their properties. Additionally, we use ML tools to predict band-alignment information for 2D materials. The web-apps, tools and associated data will be distributed through JARVIS-Heterostructure website (https://jarvis.nist.gov/jarvish/). Our analysis, results and the developed web-apps can be applied to the screening and design applications, such as finding novel photocatalysts, photodetectors, and high-work function (WF) 2D-metal contacts.
△ Less
Submitted 12 September, 2022; v1 submitted 6 April, 2020;
originally announced April 2020.
-
Topological surface states of MnBi$_2$Te$_4$ at finite temperatures and at domain walls
Authors:
Kevin F. Garrity,
Sugata Chowdhury,
Francesca M. Tavazza
Abstract:
MnBi$_2$Te$_4$ has recently been the subject of intensive study, due to the prediction of axion insulator, Weyl semimetal, and quantum anomalous Hall insulator phases, depending on the structure and magnetic ordering. Experimental results have confirmed some aspects of this picture, but several experiments have seen zero-gap surfaces states at low temperature, in conflict with expectations. In thi…
▽ More
MnBi$_2$Te$_4$ has recently been the subject of intensive study, due to the prediction of axion insulator, Weyl semimetal, and quantum anomalous Hall insulator phases, depending on the structure and magnetic ordering. Experimental results have confirmed some aspects of this picture, but several experiments have seen zero-gap surfaces states at low temperature, in conflict with expectations. In this work, we develop a first-principles-based tight-binding model that allows for arbitrary control of the local spin direction and spin-orbit coupling, enabling us to accurately treat large unit-cells. Using this model, we examine the behavior of the topological surface state as a function of temperature, finding a gap closure only above the Néel temperature. In addition, we examine the effect of magnetic domains on the electronic structure, and we find that the domain wall zero-gap states extend over many unit-cells. These domain wall states can appear similar to the high temperature topological surface state when many domain sizes are averaged, potentially reconciling theoretical results with experiments.
△ Less
Submitted 24 February, 2021; v1 submitted 6 March, 2020;
originally announced March 2020.
-
Computational Search for Magnetic and Non-magnetic 2D Topological Materials using Unified Spin-orbit Spillage Screening
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Jie Jiang,
Ruth Pachter,
Francesca Tavazza
Abstract:
Two-dimensional topological materials (TMs) have a variety of properties that make them attractive for applications including spintronics and quantum computation. However, there are only a few such experimentally known materials. To help discover new 2D TMs, we develop a unified and computationally inexpensive approach to identify magnetic and non-magnetic 2D TMs, including gapped and semi-metalli…
▽ More
Two-dimensional topological materials (TMs) have a variety of properties that make them attractive for applications including spintronics and quantum computation. However, there are only a few such experimentally known materials. To help discover new 2D TMs, we develop a unified and computationally inexpensive approach to identify magnetic and non-magnetic 2D TMs, including gapped and semi-metallic topological classifications, in a high-throughput way using density functional theory-based spin-orbit spillage, Wannier-interpolation, and related techniques. We first compute the spin-orbit spillage for the ~1000 2D materials in the JARVIS-DFT dataset (https://www.ctcms.nist.gov/~knc6/JVASP.html ), resulting in 122 materials with high-spillage values. Then, we use Wannier-interpolation to carry-out Z2, Chern-number, anomalous Hall conductivity, Curie temperature, and edge state calculations to further support the predictions. We identify various topologically non-trivial classes such as quantum spin-hall insulators (QSHI), quantum anomalous-hall insulators (QAHI), and semimetals. For a few predicted materials, we run G0W0+SOC and DFT+U calculations. We find that as we introduce many-body effects, only a few materials retain non-trivial band-topology, suggesting the importance of high-level DFT methods in predicting 2D topological materials. However, as an initial step, the automated spillage screening and Wannier-approach provide useful predictions for finding new topological materials and to narrow down candidates for experimental synthesis and characterization.
△ Less
Submitted 11 February, 2020; v1 submitted 30 January, 2020;
originally announced January 2020.
-
Uncertainty Quantification for Materials Properties in Density Functional Theory with k-Point Density
Authors:
Joshua J. Gabriel,
Faical Yannick C. Congo,
Alexander Sinnott,
Kiran Mathew,
Thomas C. Allison,
Francesca Tavazza,
Richard G. Hennig
Abstract:
Many computational databases emerged over the last five years that report material properties calculated with density functional theory. The properties in these databases are commonly calculated to a precision that is set by choice of the basis set and the k-point density for the Brillouin zone integration. We determine how the precision of properties obtained from the Birch equation of state for…
▽ More
Many computational databases emerged over the last five years that report material properties calculated with density functional theory. The properties in these databases are commonly calculated to a precision that is set by choice of the basis set and the k-point density for the Brillouin zone integration. We determine how the precision of properties obtained from the Birch equation of state for 29 transition metals and aluminum in the three common structures -- fcc, bcc, and hcp -- correlate with the k-point density and the precision of the energy. We show that the precision of the equilibrium volume, bulk modulus, and the pressure derivative of the bulk modulus correlate comparably well with the k-point density and the precision of the energy, following an approximate power law. We recommend the k-point density as the convergence parameter because it is computationally efficient, easy to use as a direct input parameter, and correlates with property precision at least as well as the energy precision. We predict that common k-point density choices in high throughput DFT databases result in precision for the volume of 0.1%, the bulk modulus of 1%, and the pressure derivative of 10%.
△ Less
Submitted 6 January, 2020;
originally announced January 2020.
-
Density Functional Theory and Deep-learning to Accelerate Data Analytics in Scanning Tunneling Microscopy
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Charles Camp,
Sergei V. Kalinin,
Rama Vasudevan,
Maxim Ziatdinov,
Francesca Tavazza
Abstract:
We introduce the first systematic database of scanning tunneling microscope (STM) images obtained using density functional theory (DFT) for two-dimensional (2D) materials, calculated using the Tersoff-Hamann method. It currently contains data for 716 exfoliable 2D materials. Examples of the five possible Bravais lattice types for 2D materials and their Fourier-transforms are discussed. All the com…
▽ More
We introduce the first systematic database of scanning tunneling microscope (STM) images obtained using density functional theory (DFT) for two-dimensional (2D) materials, calculated using the Tersoff-Hamann method. It currently contains data for 716 exfoliable 2D materials. Examples of the five possible Bravais lattice types for 2D materials and their Fourier-transforms are discussed. All the computational STM images generated in this work will be made available on the JARVIS-DFT website (https://www.ctcms.nist.gov/~knc6/JVASP.html). We find excellent qualitative agreement between the computational and experimental STM images for selected materials. As a first example application of this database, we train a convolution neural network (CNN) model to identify Bravais lattices from the STM images. We believe the model can aid high-throughput experimental data analysis. These computational STM images can directly aid the identification of phases, analyzing defects and lattice-distortions in experimental STM images, as well as be incorporated in the autonomous experiment workflows.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
High-throughput Density Functional Perturbation Theory and Machine Learning Predictions of Infrared, Piezoelectric and Dielectric Responses
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Vinit Sharma,
Adam J. Biacchi,
Angela R. Hight Walker,
Francesca Tavazza
Abstract:
Many technological applications depend on the response of materials to electric fields, but available databases of such responses are limited. Here, we explore the infrared, piezoelectric and dielectric properties of inorganic materials by combining high-throughput density functional perturbation theory and machine learning approaches. We compute Γ-point phonons, infrared intensities, Born-effecti…
▽ More
Many technological applications depend on the response of materials to electric fields, but available databases of such responses are limited. Here, we explore the infrared, piezoelectric and dielectric properties of inorganic materials by combining high-throughput density functional perturbation theory and machine learning approaches. We compute Γ-point phonons, infrared intensities, Born-effective charges, piezoelectric, and dielectric tensors for 5015 non-metallic materials in the JARVIS-DFT database. We find 3230 and 1943 materials with at least one far and mid-infrared mode, respectively. We identify 577 high-piezoelectric materials, using a threshold of 0.5 C/m2. Using a threshold of 20, we find 593 potential high-dielectric materials. Importantly, we analyze the chemistry, symmetry, dimensionality, and geometry of the materials to find features that help explain variations in our datasets. Finally, we develop high-accuracy regression models for the highest infrared frequency and maximum Born-effective charges, and classification models for maximum piezoelectric and average dielectric tensors to accelerate discovery.
△ Less
Submitted 29 February, 2020; v1 submitted 2 October, 2019;
originally announced October 2019.
-
Data-driven Discovery of 3D and 2D Thermoelectric Materials
Authors:
Kamal Choudhary,
Kevin Garrity,
Francesca Tavazza
Abstract:
In this work, we first perform a systematic search for high-efficiency three-dimensional (3D) and two-dimensional (2D) thermoelectric materials by combining semiclassical transport techniques with density functional theory (DFT) calculations and then train machine-learning models on the thermoelectric data. Out of 36000 three-dimensional and 900 two-dimensional materials currently in the publicly…
▽ More
In this work, we first perform a systematic search for high-efficiency three-dimensional (3D) and two-dimensional (2D) thermoelectric materials by combining semiclassical transport techniques with density functional theory (DFT) calculations and then train machine-learning models on the thermoelectric data. Out of 36000 three-dimensional and 900 two-dimensional materials currently in the publicly available JARVIS-DFT database, we identify 2932 3D and 148 2D promising thermoelectric materials using a multi-steps screening procedure, where specific thresholds are chosen for key quantities like bandgaps, Seebeck coefficients and power factors. We compute the Seebeck coefficients for all the materials currently in the database and validate our calculations by comparing our results, for a subset of materials, to experimental and existing computational datasets. We also investigate the effect of chemical, structural, crystallographic and dimensionality trends on thermoelectric performance. We predict several classes of efficient 3D and 2D materials such as Ba(MgX)2 (X=P,As,Bi), X2YZ6 (X=K,Rb, Y=Pd,Pt, Z=Cl,Br), K2PtX2(X=S,Se), NbCu3X4 (X=S,Se,Te), Sr2XYO6 (X=Ta, Zn, Y=Ga, Mo), TaCu3X4 (X=S, Se,Te), and XYN (X=Ti, Zr, Y=Cl, Br). Finally, as high-throughput DFT is computationally expensive, we train machine learning models using gradient boosting decision trees (GBDT) and classical force-field inspired descriptors (CFID) for n-and p-type Seebeck coefficients and power factors, to quickly pre-screen materials for guiding the next set of DFT calculations. The dataset and tools are made publicly available at the websites: https://www.ctcms.nist.gov/~knc6/JVASP.html , https://www.ctcms.nist.gov/jarvisml/ and https://jarvis.nist.gov/ .
△ Less
Submitted 30 January, 2020; v1 submitted 14 June, 2019;
originally announced June 2019.
-
Accelerated Discovery of Efficient Solar-cell Materials using Quantum and Machine-learning Methods
Authors:
Kamal Choudhary,
Marnik Bercx,
Jie Jiang,
Ruth Pachter,
Dirk Lamoen,
Francesca Tavazza
Abstract:
Solar-energy plays an important role in solving serious environmental problems and meeting high-energy demand. However, the lack of suitable materials hinders further progress of this technology. Here, we present the largest inorganic solar-cell material search to date using density functional theory (DFT) and machine-learning approaches. We calculated the spectroscopic limited maximum efficiency…
▽ More
Solar-energy plays an important role in solving serious environmental problems and meeting high-energy demand. However, the lack of suitable materials hinders further progress of this technology. Here, we present the largest inorganic solar-cell material search to date using density functional theory (DFT) and machine-learning approaches. We calculated the spectroscopic limited maximum efficiency (SLME) using Tran-Blaha modified Becke-Johnson potential for 5097 non-metallic materials and identified 1997 candidates with an SLME higher than 10%, including 934 candidates with suitable convex-hull stability and effective carrier mass. Screening for 2D-layered cases, we found 58 potential materials and performed G0W0 calculations on a subset to estimate the prediction-uncertainty. As the above DFT methods are still computationally expensive, we developed a high accuracy machine learning model to pre-screen efficient materials and applied it to over a million materials. Our results provide a general framework and universal strategy for the design of high-efficiency solar cell materials. The data and tools are publicly distributed at: https://www.ctcms.nist.gov/~knc6/JVASP.html, https://www.ctcms.nist.gov/jarvisml/, https://jarvis.nist.gov/ and https://github.com/usnistgov/jarvis .
△ Less
Submitted 6 August, 2019; v1 submitted 15 March, 2019;
originally announced March 2019.
-
Prediction of Weyl semimetal, AFM topological insulator, nodal line semimetal, and Chern insulator phases in Bi2MnSe4
Authors:
Sugata Chowdhury,
Kevin F. Garrity,
Francesca Tavazza
Abstract:
Three dimensional materials with strong spin-orbit coupling and magnetic interactions represent an opportunity to realize a variety of rare and potentially useful topological phases. In this work, we use first principles calculations to show that the recently synthesized material Bi2MnSe4 displays a combination of band inversion and magnetic interactions, leading to several topological phases. Bi2…
▽ More
Three dimensional materials with strong spin-orbit coupling and magnetic interactions represent an opportunity to realize a variety of rare and potentially useful topological phases. In this work, we use first principles calculations to show that the recently synthesized material Bi2MnSe4 displays a combination of band inversion and magnetic interactions, leading to several topological phases. Bi2PbSe4, also studied, also displays band inversion and is a topological insulator. In bulk form, the ferromagnetic phase of Bi2MnSe4 is either a nodal line or Weyl semimetal, depending on the direction of the spins. When the spins are arranged in a layered antiferromagnetic configuration, the combination of time reversal plus a partial translation is a new symmetry, and the material instead becomes an antiferromagnetic topological insulator. However, the intrinsic TRS breaking at the surface of Bi2MnSe4 removes the typical Dirac cone feature, allowing the observation of the half-integer quantum anomalous Hall effect (AHC). Furthermore, we show that in thin film form, for some thicknesses, Bi2MnSe4 becomes a Chern insulator with a band gap of up to 58 meV. This combination of properties in a stoichiometric magnetic material makes Bi2MnSe4 an excellent candidate for displaying robust topological behavior.
△ Less
Submitted 5 November, 2018;
originally announced November 2018.
-
High-throughput Discovery of Topologically Non-trivial Materials using Spin-orbit Spillage
Authors:
Kamal Choudhary,
Kevin F. Garrity,
Francesca Tavazza
Abstract:
We present a novel methodology to identify topologically non-trivial materials based on band inversion induced by spin-orbit coupling (SOC) effect. Specifically, we compare the density functional theory (DFT) based wavefunctions with and without spin-orbit coupling and compute the spin-orbit-spillage as a measure of band-inversion. Due to its ease of calculation, without any need for symmetry anal…
▽ More
We present a novel methodology to identify topologically non-trivial materials based on band inversion induced by spin-orbit coupling (SOC) effect. Specifically, we compare the density functional theory (DFT) based wavefunctions with and without spin-orbit coupling and compute the spin-orbit-spillage as a measure of band-inversion. Due to its ease of calculation, without any need for symmetry analysis or dense k-point interpolation, the spillage is an excellent tool for identifying topologically non-trivial materials. Out of 30000 materials available in the JARVIS-DFT database, we applied this methodology to more than 4835 non-magnetic materials consisting of heavy atoms and low bandgaps. We found 1868 candidate materials with high-spillage (using 0.5 as a threshold). We validated our methodology by carrying out conventional Wannier-interpolation calculations for 289 candidate materials. We demonstrate that in addition to Z2 topological insulators, this screening method successfully identified many semimetals and topological crystalline insulators. Importantly, our approach is applicable to the investigation of disordered or distorted as well as magnetic materials, because it is not based on symmetry considerations. We discuss some individual example materials, as well as trends throughout our dataset, which is available at the websites: https://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov/.
△ Less
Submitted 3 June, 2019; v1 submitted 24 October, 2018;
originally announced October 2018.
-
Convergence and machine learning predictions of Monkhorst-Pack k-points and plane-wave cut-off in high-throughput DFT calculations
Authors:
Kamal Choudhary,
Francesca Tavazza
Abstract:
In this work, we developed an automatic convergence procedure for k-points and plane wave cut-off in density functional (DFT) calculations and applied it to more than 30000 materials. The computational framework for automatic convergence can take a user-defined input as a convergence criterion. For k-points, we converged energy per cell (EPC) to 0.001 eV/cell tolerance and compared the results wit…
▽ More
In this work, we developed an automatic convergence procedure for k-points and plane wave cut-off in density functional (DFT) calculations and applied it to more than 30000 materials. The computational framework for automatic convergence can take a user-defined input as a convergence criterion. For k-points, we converged energy per cell (EPC) to 0.001 eV/cell tolerance and compared the results with those obtained using an energy per atom (EPA) convergence criteria of 0.001 eV/atom. From the analysis of our results, we could relate k-point density and plane wave cut-off to material parameters such as density, the slope of bands, number of band-crossings, the maximum plane-wave cut-off used in pseudopotential generation, crystal systems, and the number of unique species in materials. We also identified some material species that would require more careful convergence than others. Moreover, we statistically investigated the dependence of k-points and cutoff on exchange-correlation functionals. We utilized all this data to train machine learning models to predict the k-point line density and plane-wave cut-off for generalized materials. This would provide users with a good starting point towards converged DFT calculations. The code used, and the converged data are available on the following websites: https://jarvis.nist.gov/, and https://github.com/usnistgov/jarvis .
△ Less
Submitted 8 February, 2019; v1 submitted 5 September, 2018;
originally announced September 2018.
-
Machine learning with force-field inspired descriptors for materials: fast screening and mapping energy landscape
Authors:
Kamal Choudhary,
Brian DeCost,
Francesca Tavazza
Abstract:
We present a complete set of chemo-structural descriptors to significantly extend the applicability of machine-learning (ML) in material screening and mapping energy landscape for multicomponent systems. These new descriptors allow differentiating between structural prototypes, which is not possible using the commonly used chemical-only descriptors. Specifically, we demonstrate that the combinatio…
▽ More
We present a complete set of chemo-structural descriptors to significantly extend the applicability of machine-learning (ML) in material screening and mapping energy landscape for multicomponent systems. These new descriptors allow differentiating between structural prototypes, which is not possible using the commonly used chemical-only descriptors. Specifically, we demonstrate that the combination of pairwise radial, nearest neighbor, bond-angle, dihedral-angle and core-charge distributions plays an important role in predicting formation energies, bandgaps, static refractive indices, magnetic properties, and modulus of elasticity for three-dimensional (3D) materials as well as exfoliation energies of two-dimensional (2D) layered materials. The training data consists of 24549 bulk and 616 monolayer materials taken from JARVIS-DFT database. We obtained very accurate ML models using gradient boosting algorithm. Then we use the trained models to discover exfoliable 2D-layered materials satisfying specific property requirements. Additionally, we integrate our formation energy ML model with a genetic algorithm for structure search to verify if the ML model reproduces the DFT convex hull. This verification establishes a more stringent evaluation metric for the ML model than what commonly used in data sciences. Our learnt model is publicly available on the JARVIS-ML website (https://www.ctcms.nist.gov/jarvisml ) property predictions of generalized materials.
△ Less
Submitted 19 July, 2018; v1 submitted 18 May, 2018;
originally announced May 2018.
-
Elastic properties of bulk and low-dimensional materials using Van der Waals density functional
Authors:
Kamal Choudhary,
Gowoon Cheon,
Evan Reed,
Francesca Tavazza
Abstract:
In this work, we present a high-throughput first-principles study of elastic properties of bulk and monolayer materials mainly using the vdW-DF-optB88 functional. We discuss the trends on the elastic response with respect to changes in dimensionality. We identify a relation between exfoliation energy and elastic constants for layered materials that can help to guide the search for vdW bonding in m…
▽ More
In this work, we present a high-throughput first-principles study of elastic properties of bulk and monolayer materials mainly using the vdW-DF-optB88 functional. We discuss the trends on the elastic response with respect to changes in dimensionality. We identify a relation between exfoliation energy and elastic constants for layered materials that can help to guide the search for vdW bonding in materials. We also predicted a few novel materials with auxetic behavior. The uncertainty in structural and elastic properties due to the inclusion of vdW interactions is discussed. We investigated 11,067 bulk and 257 monolayer materials. Lastly, we found that the trends in elastic constants for bulk and their monolayer counterparts can be very different. All the computational results are made publicly available at easy-to-use websites: https://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov/ . Our dataset can be used to identify stiff and flexible materials for industrial applications.
△ Less
Submitted 29 June, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
Computational screening of high-performance optoelectronic materials using OptB88vdW and TB-mBJ formalisms
Authors:
Kamal Choudhary,
Qin Zhang,
Andrew C. E. Reid,
Sugata Chowdhury,
Nhan Van Nguyen,
Zachary Trautt,
Marcus W. Newrock,
Faical Yannick Congo,
Francesca Tavazza
Abstract:
We perform high-throughput density functional theory (DFT) calculations for optoelectronic properties (electronic bandgap and frequency dependent dielectric function) using the OptB88vdW functional (OPT) and the Tran-Blaha modified Becke Johnson potential (MBJ). This data is distributed publicly through JARVIS-DFT database. We used this data to evaluate the differences between these two formalisms…
▽ More
We perform high-throughput density functional theory (DFT) calculations for optoelectronic properties (electronic bandgap and frequency dependent dielectric function) using the OptB88vdW functional (OPT) and the Tran-Blaha modified Becke Johnson potential (MBJ). This data is distributed publicly through JARVIS-DFT database. We used this data to evaluate the differences between these two formalisms and quantify their accuracy, comparing to experimental data whenever applicable. At present, we have 17,805 OPT and 7,358 MBJ bandgaps and dielectric functions. MBJ is found to predict better bandgaps and dielectric functions than OPT, so it can be used to improve the well-known bandgap problem of DFT in a relatively inexpensive way. The peak positions in dielectric functions obtained with OPT and MBJ are in comparable agreement with experiments. The data is available on our websites http://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov .
△ Less
Submitted 3 April, 2018;
originally announced April 2018.
-
High-throughput assessment of vacancy formation and surface energies of materials using classical force-fields
Authors:
Kamal Choudhary,
Adam J. Biacchi,
Supriyo Ghosh,
Lucas Hale,
Angela R. Hight Walker,
Francesca Tavazza
Abstract:
In this work, we present an open access database for surface and vacancy-formation energies using classical force-fields (FFs). These quantities are essential in understanding diffusion behavior, nanoparticle formation and catalytic activities. FFs are often designed for a specific application, hence, this database allows the user to understand whether a FF is suitable for investigating particular…
▽ More
In this work, we present an open access database for surface and vacancy-formation energies using classical force-fields (FFs). These quantities are essential in understanding diffusion behavior, nanoparticle formation and catalytic activities. FFs are often designed for a specific application, hence, this database allows the user to understand whether a FF is suitable for investigating particular defect and surface-related material properties. The FF results are compared to density functional theory and experimental data whenever applicable for validation. At present, we have 17,506 surface energies and 1,000 vacancy formation energies calculation in our database and the database is still growing. All the data generated, and the computational tools used, are shared publicly at the following websites https://www.ctcms.nist.gov/~knc6/periodic.html, https://jarvis.nist.gov and https://github.com/usnistgov/jarvis . Approximations used during the high-throughput calculations are clearly mentioned. Using some of the example cases, we show how our data can be used to directly compare different FFs for a material and to interpret experimental findings such as using Wulff construction for predicting equilibrium shape of nanoparticles. Similarly, the vacancy formation energies data can be useful in understanding diffusion related properties.
△ Less
Submitted 9 October, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
Evolution of Raman spectra in Mo$_{1-x}$W$_x$Te$_2$ alloys
Authors:
Sean M. Oliver,
Ryan Beams,
Sergiy Krylyuk,
Irina Kalish,
Arunima K. Singh,
Alina Bruma,
Francesca Tavazza,
Jaydeep Joshi,
Iris R. Stone,
Stephan J. Stranick,
Albert V. Davydov,
Patrick M. Vora
Abstract:
The structural polymorphism in transition metal dichalcogenides (TMDs) provides exciting opportunities for developing advanced electronics. For example, MoTe$_2$ crystallizes in the 2H semiconducting phase at ambient temperature and pressure, but transitions into the 1T$^\prime$ semimetallic phase at high temperatures. Alloying MoTe$_2$ with WTe$_2$ reduces the energy barrier between these two pha…
▽ More
The structural polymorphism in transition metal dichalcogenides (TMDs) provides exciting opportunities for developing advanced electronics. For example, MoTe$_2$ crystallizes in the 2H semiconducting phase at ambient temperature and pressure, but transitions into the 1T$^\prime$ semimetallic phase at high temperatures. Alloying MoTe$_2$ with WTe$_2$ reduces the energy barrier between these two phases, while also allowing access to the T$_d$ Weyl semimetal phase. The MoWTe$_2$ alloy system is therefore promising for developing phase change memory technology. However, achieving this goal necessitates a detailed understanding of the phase composition in the MoTe$_2$-WTe$_2$ system. We combine polarization-resolved Raman spectroscopy with X-ray diffraction (XRD) and scanning transmission electron microscopy (STEM) to study MoWTe$_2$ alloys over the full compositional range x from 0 to 1. We identify Raman and XRD signatures characteristic of the 2H, 1T$^\prime$, and T$_d$ structural phases that agree with density-functional theory (DFT) calculations, and use them to identify phase fields in the MoTe$_2$-WTe$_2$ system, including single-phase 2H, 1T$^\prime$, and T$_d$ regions, as well as a two-phase 1T$^\prime$ + T$_d$ region. Disorder arising from compositional fluctuations in MoWTe$_2$ alloys breaks inversion and translational symmetry, leading to the activation of an infrared 1T$^\prime$-MoTe$_2$ mode and the enhancement of a double-resonance Raman process in 2H-MoWTe$_2$ alloys. Compositional fluctuations limit the phonon correlation length, which we estimate by fitting the observed asymmetric Raman lineshapes with a phonon confinement model. These observations reveal the important role of disorder in MoWTe$_2$ alloys, clarify the structural phase boundaries, and provide a foundation for future explorations of phase transitions and electronic phenomena in this system.
△ Less
Submitted 31 March, 2017;
originally announced March 2017.
-
MPInterfaces: A Materials Project based Python Tool for High-Throughput Computational Screening of Interfacial Systems
Authors:
Kiran Mathew,
Arunima K. Singh,
Joshua J. Gabriel,
Kamal Choudhary,
Susan B. Sinnott,
Albert V. Davydov,
Francesca Tavazza,
Richard G. Hennig
Abstract:
A Materials Project based open-source Python tool, MPInterfaces, has been developed to automate the high-throughput computational screening and study of interfacial systems. The framework encompasses creation and manipulation of interface structures for solid/solid hetero-structures, solid/implicit solvents systems, nanoparticle/ligands systems; and the creation of simple system-agnostic workflows…
▽ More
A Materials Project based open-source Python tool, MPInterfaces, has been developed to automate the high-throughput computational screening and study of interfacial systems. The framework encompasses creation and manipulation of interface structures for solid/solid hetero-structures, solid/implicit solvents systems, nanoparticle/ligands systems; and the creation of simple system-agnostic workflows for in depth computational analysis using density-functional theory or empirical energy models. The package leverages existing open-source high-throughput tools and extends their capabilities towards the understanding of interfacial systems. We describe the various algorithms and methods implemented in the package. Using several test cases, we demonstrate how the package enables high-throughput computational screening of advanced materials, directly contributing to the Materials Genome Initiative (MGI), which aims to accelerate the discovery, development, and deployment of new materials.
△ Less
Submitted 24 February, 2016;
originally announced February 2016.
-
Crystal asymmetry induces single-atom chain formation in gold nanowires
Authors:
F. Tavazza,
A. Hasmy,
L. E. Levine,
A. M. Chaka,
L. Rincon,
M. Marquez,
C. Gonzalez
Abstract:
We performed density functional theory and tight-binding molecular dynamics calculations to investigate the formation of suspended linear atom chains when stretching gold nanowires along the [110] crystal orientation. We determined that chain formation can occur only when the crystal symmetry is broken in the early stages of the elongation process. Such crystallographic asymmetry can be induced…
▽ More
We performed density functional theory and tight-binding molecular dynamics calculations to investigate the formation of suspended linear atom chains when stretching gold nanowires along the [110] crystal orientation. We determined that chain formation can occur only when the crystal symmetry is broken in the early stages of the elongation process. Such crystallographic asymmetry can be induced by stretching the wire along a slightly off-axis tensile direction or by introducing thermal fluctuations. Our observation of the off-axis formation of these chains agrees with experimental findings.
△ Less
Submitted 12 March, 2007;
originally announced March 2007.