Showing 1–39 of 39 results for author: Persson, K A

Overcoming systematic softening in universal machine learning interatomic potentials by fine-tuning

Authors: Bowen Deng, Yunyeong Choi, Peichen Zhong, Janosh Riebesell, Shashwat Anand, Zhuohan Li, KyuJung Jun, Kristin A. Persson, Gerbrand Ceder

Abstract: Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have seen the emergence of universal MLIPs (uMLIPs) that are pre-trained on diverse materials datasets, providing opportunities for both ready-to-use universal force fields and robust foundations for downstream machine learning refinements. However, their performance in extrap… ▽ More Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have seen the emergence of universal MLIPs (uMLIPs) that are pre-trained on diverse materials datasets, providing opportunities for both ready-to-use universal force fields and robust foundations for downstream machine learning refinements. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force under-prediction in a series of atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, phonon vibration modes, ion migration barriers, and general high-energy states. We find that the PES softening behavior originates from a systematic underprediction error of the PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. We demonstrate that the PES softening issue can be effectively rectified by fine-tuning with a single additional data point. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. This result rationalizes the data-efficient fine-tuning performance boost commonly observed with foundational MLIPs. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs. △ Less

Submitted 11 May, 2024; originally announced May 2024.

Developments and applications of the OPTIMADE API for materials discovery, design, and data exchange

Authors: Matthew L. Evans, Johan Bergsma, Andrius Merkys, Casper W. Andersen, Oskar B. Andersson, Daniel Beltrán, Evgeny Blokhin, Tara M. Boland, Rubén Castañeda Balderas, Kamal Choudhary, Alberto Díaz Díaz, Rodrigo Domínguez García, Hagen Eckert, Kristjan Eimre, María Elena Fuentes Montero, Adam M. Krajewski, Jens Jørgen Mortensen, José Manuel Nápoles Duarte, Jacob Pietryga, Ji Qi, Felipe de Jesús Trejo Carrillo, Antanas Vaitkus, Jusong Yu, Adam Zettel, Pedro Baptista de Castro , et al. (34 additional authors not shown)

Abstract: The Open Databases Integration for Materials Design (OPTIMADE) application programming interface (API) empowers users with holistic access to a growing federation of databases, enhancing the accessibility and discoverability of materials and chemical data. Since the first release of the OPTIMADE specification (v1.0), the API has undergone significant development, leading to the upcoming v1.2 relea… ▽ More The Open Databases Integration for Materials Design (OPTIMADE) application programming interface (API) empowers users with holistic access to a growing federation of databases, enhancing the accessibility and discoverability of materials and chemical data. Since the first release of the OPTIMADE specification (v1.0), the API has undergone significant development, leading to the upcoming v1.2 release, and has underpinned multiple scientific studies. In this work, we highlight the latest features of the API format, accompanying software tools, and provide an update on the implementation of OPTIMADE in contributing materials databases. We end by providing several use cases that demonstrate the utility of the OPTIMADE API in materials research that continue to drive its ongoing development. △ Less

Submitted 5 April, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

The ab initio amorphous materials database: Empowering machine learning to decode diffusivity

Authors: Hui Zheng, Eric Sivonxay, Max Gallant, Ziyao Luo, Matthew McDermott, Patrick Huck, Kristin A. Persson

Abstract: Amorphous materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of amorphous materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the… ▽ More Amorphous materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of amorphous materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the largest computed amorphous materials database to date, generated from systematic and accurate \textit{ab initio} molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductivity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching amorphous materials provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in universal machine learning potentials, impacting design beyond that of non-crystalline materials. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 28 pages, 7 figures

A foundation model for atomistic materials chemistry

Authors: Ilyes Batatia, Philipp Benner, Yuan Chiang, Alin M. Elena, Dávid P. Kovács, Janosh Riebesell, Xavier R. Advincula, Mark Asta, Matthew Avaylon, William J. Baldwin, Fabian Berger, Noam Bernstein, Arghya Bhowmik, Samuel M. Blau, Vlad Cărare, James P. Darby, Sandip De, Flaviano Della Pia, Volker L. Deringer, Rokas Elijošius, Zakariya El-Machachi, Fabio Falcioni, Edvin Fako, Andrea C. Ferrari, Annalena Genreith-Schriever , et al. (51 additional authors not shown)

Abstract: Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferabilit… ▽ More Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry. △ Less

Submitted 1 March, 2024; v1 submitted 29 December, 2023; originally announced January 2024.

Comments: 119 pages, 63 figures, 37MB PDF

Persistent anisotropy of the spin cycloid in BiFeO3 through ferroelectric switching

Authors: Peter Meisenheimer, Guy Moore, Shiyu Zhou, Hongrui Zhang, Xiaoxi Huang, Sajid Husain, Xianzhe Chen, Lane W. Martin, Kristin A. Persson, Sinéad Griffin, Lucas Caretta, Paul Stevenson, Ramamoorthy Ramesh

Abstract: A key challenge in antiferromagnetic spintronics is the control of spin configuration on nanometer scales applicable to solid-state technologies. Bismuth ferrite (BiFeO3) is a multiferroic material that exhibits both ferroelectricity and canted antiferromagnetism at room temperature, making it a unique candidate in the development of electric-field controllable magnetic devices. The magnetic momen… ▽ More A key challenge in antiferromagnetic spintronics is the control of spin configuration on nanometer scales applicable to solid-state technologies. Bismuth ferrite (BiFeO3) is a multiferroic material that exhibits both ferroelectricity and canted antiferromagnetism at room temperature, making it a unique candidate in the development of electric-field controllable magnetic devices. The magnetic moments in BiFeO3 are arranged into a spin cycloid, resulting in unique magnetic properties which are tied to the ferroelectric order. Previous understanding of this coupling has relied on average, mesoscale measurements to infer behavior. Using nitrogen vacancy-based diamond magnetometry, we show that the spin cycloid can be deterministically controlled with an electric field. The energy landscape of the cycloid is shaped by both the ferroelectric degree of freedom and strain-induced anisotropy, restricting the magnetization changes to specific ferroelectric switching events. This study provides understanding of the antiferromagnetic texture in BiFeO3 and paves new avenues for designing magnetic textures and spintronic devices. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 28 pages, 4 main text figures, 8 supp figures

Journal ref: Nature Communications 15, 2903 (2024)

Mechanistically-guided materials chemistry: synthesis of new ternary nitrides, CaZrN$_2$ and CaHfN$_2$

Authors: Christopher L. Rom, Andrew Novick, Matthew J. McDermott, Andrey A. Yakovenko, Jessica R. Gallawa, Gia Thinh Tran, Dominic C. Asebiah, Emily N. Storck, Brennan C. McBride, Rebecca C. Miller, Amy L. Prieto, Kristin A. Persson, Eric Toberer, Vladan Stevanović, Andriy Zakutayev, James R. Neilson

Abstract: Recent computational studies have predicted many new ternary nitrides, revealing synthetic opportunities in this underexplored phase space. However, synthesizing new ternary nitrides is difficult, in part because intermediate and product phases often have high cohesive energies that inhibit diffusion. Here, we report the synthesis of two new phases, calcium zirconium nitride (CaZrN$_2$) and calciu… ▽ More Recent computational studies have predicted many new ternary nitrides, revealing synthetic opportunities in this underexplored phase space. However, synthesizing new ternary nitrides is difficult, in part because intermediate and product phases often have high cohesive energies that inhibit diffusion. Here, we report the synthesis of two new phases, calcium zirconium nitride (CaZrN$_2$) and calcium hafnium nitride (CaHfN$_2$), by solid state metathesis reactions between Ca$_3$N$_2$ and $M$Cl$_4$ ($M$ = Zr, Hf). Although the reaction nominally proceeds to the target phases in a 1:1 ratio of the precursors via Ca$_3$N$_2$ + $M$Cl$_4$ $\rightarrow$ Ca$M$N$_2$ + 2 CaCl$_2$, reactions prepared this way result in Ca-poor materials (Ca$_xM_{2-x}$N$_2$, $x<1$). A small excess of Ca$_3$N$_2$ (ca. 20 mol\%) is needed to yield stoichiometric Ca$M$N$_2$, as confirmed by high-resolution synchrotron powder X-ray diffraction. In situ synchrotron X-ray diffraction studies reveal that nominally stoichiometric reactions produce Zr$^{3+}$ intermediates early in the reaction pathway, and the excess Ca$_3$N$_2$ is needed to reoxidize Zr$^{3+}$ intermediates back to the Zr$^{4+}$ oxidation state of CaZrN$_2$. Analysis of computationally-derived chemical potential diagrams rationalizes this synthetic approach and its contrast from the synthesis of MgZrN$_2$. These findings additionally highlight the utility of in situ diffraction studies and computational thermochemistry to provide mechanistic guidance for synthesis. △ Less

Submitted 30 October, 2023; originally announced October 2023.

SpinPSO: An agent-based optimization workflow for identifying global noncollinear magnetic ground-states from first-principles

Authors: Guy C. Moore, Matthew K. Horton, Kristin A. Persson

Abstract: We propose and implement a novel hybrid meta-heuristic optimization algorithm for the identification of non-collinear global ground-states in magnetic systems. The inputs to this optimization scheme are directly from non-collinear density functional theory (DFT), and the workflow is implemented in the atomate code framework, making it suitable to run on high-performance computing architectures. Th… ▽ More We propose and implement a novel hybrid meta-heuristic optimization algorithm for the identification of non-collinear global ground-states in magnetic systems. The inputs to this optimization scheme are directly from non-collinear density functional theory (DFT), and the workflow is implemented in the atomate code framework, making it suitable to run on high-performance computing architectures. The hybrid algorithm provides a seamless theoretical extension of particle swarm optimization (PSO) algorithms to continuous $\mathcal S^2$ spins, giving it the name SpinPSO. The hybrid nature of the algorithm stems from setting the dynamics of individual spins to be governed by physically motivated atomistic spin dynamics. Using this algorithm, we are able to achieve convergence to experimentally resolved magnetic ground-states for a set of diverse test case materials that exhibit exotic spin textures. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Realistic non-collinear ground states of solids with source-free exchange correlation functional

Authors: Guy C. Moore, Matthew K. Horton, Aaron D. Kaplan, Sinéad M. Griffin, Kristin A. Persson

Abstract: In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier t… ▽ More In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+$U$+$J$ functional. In the process, we compare the non-collinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analysis on probability currents, net torque, and constant terms draws additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and paves the way for future extensions and corrections to the SF corrected XC functional. △ Less

Submitted 24 May, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

Magnetism and magnetocaloric properties of Co$_{1-x}$Mn$_x$Cr$_2$O$_4$

Authors: Joya A. Cooley, Gregor Dairaghi, Guy C. Moore, Matthew K. Horton, Emily C. Schueller, Kristin A. Persson, Ram Seshadri

Abstract: Co$_{1-x}$Mn$_x$Cr$_2$O$_4$ crystallizes as a normal spinel in the cubic $Fd \overline{3}m$ space group, and the end members have been reported to display a region of collinear ferrimagnetism as well as a low-temperature spin-spiral state with variable coherence lengths from 3 nm to 10 nm in polycrystalline samples. Here, we present the synthesis of the entire solid solution, and data showing that… ▽ More Co$_{1-x}$Mn$_x$Cr$_2$O$_4$ crystallizes as a normal spinel in the cubic $Fd \overline{3}m$ space group, and the end members have been reported to display a region of collinear ferrimagnetism as well as a low-temperature spin-spiral state with variable coherence lengths from 3 nm to 10 nm in polycrystalline samples. Here, we present the synthesis of the entire solid solution, and data showing that the ferrimagnetic ordering temperature as well as the spin-spiral lock-in temperature are tunable with the Co/Mn ratio. The peak magnetocaloric entropy change was determined to be $ΔでるたS_M$ = -5.63 J kg$^{-1}$ K$^{-1}$ in an applied magnetic field change of $ΔでるたH$ = 0 T to 5 T for the Mn end-member at the ferrimagnetic ordering temperature. Using density functional theory (DFT), we explore the shortcomings of the magnetic deformation proxy to identify trends in $ΔでるたS_M$ across composition in this spinel system, and explore future extensions of theory to address these discrepancies. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Matbench Discovery -- A framework to evaluate machine learning crystal stability predictions

Authors: Janosh Riebesell, Rhys E. A. Goodall, Philipp Benner, Yuan Chiang, Bowen Deng, Alpha A. Lee, Anubhav Jain, Kristin A. Persson

Abstract: Matbench Discovery simulates the deployment of machine learning (ML) energy models in a high-throughput search for stable inorganic crystals. We address the disconnect between (i) thermodynamic stability and formation energy and (ii) in-domain vs out-of-distribution performance. Alongside this paper, we publish a Python package to aid with future model submissions and a growing online leaderboard… ▽ More Matbench Discovery simulates the deployment of machine learning (ML) energy models in a high-throughput search for stable inorganic crystals. We address the disconnect between (i) thermodynamic stability and formation energy and (ii) in-domain vs out-of-distribution performance. Alongside this paper, we publish a Python package to aid with future model submissions and a growing online leaderboard with further insights into trade-offs between various performance metrics. To answer the question which ML methodology performs best at materials discovery, our initial release explores a variety of models including random forests, graph neural networks (GNN), one-shot predictors, iterative Bayesian optimizers and universal interatomic potentials (UIP). Ranked best-to-worst by their test set F1 score on thermodynamic stability prediction, we find CHGNet > M3GNet > MACE > ALIGNN > MEGNet > CGCNN > CGCNN+P > Wrenformer > BOWSR > Voronoi tessellation fingerprints with random forest. The top 3 models are UIPs, the winning methodology for ML-guided materials discovery, achieving F1 scores of ~0.6 for crystal stability classification and discovery acceleration factors (DAF) of up to 5x on the first 10k most stable predictions compared to dummy selection from our test set. We also highlight a sharp disconnect between commonly used global regression metrics and more task-relevant classification metrics. Accurate regressors are susceptible to unexpectedly high false-positive rates if those accurate predictions lie close to the decision boundary at 0 eV/atom above the convex hull where most materials are. Our results highlight the need to focus on classification metrics that actually correlate with improved stability hit rate. △ Less

Submitted 4 February, 2024; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages, 18 figures, 4 tables

Assessing Thermodynamic Selectivity of Solid-State Reactions for the Predictive Synthesis of Inorganic Materials

Authors: Matthew J. McDermott, Brennan C. McBride, Corlyn Regier, Gia Thinh Tran, Yu Chen, Adam A. Corrao, Max C. Gallant, Gabrielle E. Kamm, Christopher J. Bartel, Karena W. Chapman, Peter G. Khalifah, Gerbrand Ceder, James R. Neilson, Kristin A. Persson

Abstract: Synthesis is a major challenge in the discovery of new inorganic materials. Currently, there is limited theoretical guidance for identifying optimal solid-state synthesis procedures. We introduce two selectivity metrics, primary and secondary competition, to assess the favorability of target/impurity phase formation in solid-state reactions. We used these metrics to analyze 3,520 solid-state react… ▽ More Synthesis is a major challenge in the discovery of new inorganic materials. Currently, there is limited theoretical guidance for identifying optimal solid-state synthesis procedures. We introduce two selectivity metrics, primary and secondary competition, to assess the favorability of target/impurity phase formation in solid-state reactions. We used these metrics to analyze 3,520 solid-state reactions in the literature, ranking existing approaches to popular target materials. Additionally, we implemented these metrics in a data-driven synthesis planning workflow and demonstrated its application in the synthesis of barium titanate (BaTiO$_3$). Using an 18-element chemical reaction network with first-principles thermodynamic data from the Materials Project, we identified 82,985 possible BaTiO$_3$ synthesis reactions and selected nine for experimental testing. Characterization of reaction pathways via synchrotron powder X-ray diffraction reveals that our selectivity metrics correlate with observed target/impurity formation. We discovered two efficient reactions using unconventional precursors (BaS/BaCl$_2$ and Na$_2$TiO$_3$) that produce BaTiO$_3$ faster and with fewer impurities than conventional methods, highlighting the importance of considering complex chemistries with additional elements during precursor selection. Our framework provides a foundation for predictive inorganic synthesis, facilitating the optimization of existing recipes and the discovery of new materials, including those not easily attainable with conventional precursors. △ Less

Submitted 27 September, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: Minor revisions (updated abstract, added discussion, added precursor purities, formatting). The first two authors contributed equally to this work. 53 pages, 7 figures

An equivariant graph neural network for the elasticity tensors of all seven crystal systems

Authors: Mingjian Wen, Matthew K. Horton, Jason M. Munro, Patrick Huck, Kristin A. Persson

Abstract: The elasticity tensor that describes the elastic response of a material to external forces is among the most fundamental properties of materials. The availability of full elasticity tensors for inorganic crystalline compounds, however, is limited due to experimental and computational challenges. Here, we report the materials tensor (MatTen) model for rapid and accurate estimation of the full fourt… ▽ More The elasticity tensor that describes the elastic response of a material to external forces is among the most fundamental properties of materials. The availability of full elasticity tensors for inorganic crystalline compounds, however, is limited due to experimental and computational challenges. Here, we report the materials tensor (MatTen) model for rapid and accurate estimation of the full fourth-rank elasticity tensors of crystals. Based on equivariant graph neural networks, MatTen satisfies the two essential requirements for elasticity tensors: independence of the frame of reference and preservation of material symmetry. Consequently, it provides a unified treatment of elasticity tensors for all seven crystal systems across diverse chemical spaces, without the need to deal with each separately.. MatTen was trained on a dataset of first-principles elasticity tensors garnered by the Materials Project over the past several years (we are releasing the data herein) and has broad applications in predicting the isotropic elastic properties of polycrystalline materials, examining the anisotropic behavior of single crystals, and discovering new materials with exceptional mechanical properties. Using MatTen, we have discovered a hundred new crystals with extremely large maximum directional Young's modulus and eleven polymorphs of elemental cubic metals with unconventional spatial orientation of Young's modulus. △ Less

Submitted 21 January, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Extracting Structured Seed-Mediated Gold Nanorod Growth Procedures from Literature with GPT-3

Authors: Nicholas Walker, John Dagdelen, Kevin Cruse, Sanghoon Lee, Samuel Gleason, Alexander Dunn, Gerbrand Ceder, A. Paul Alivisatos, Kristin A. Persson, Anubhav Jain

Abstract: Although gold nanorods have been the subject of much research, the pathways for controlling their shape and thereby their optical properties remain largely heuristically understood. Although it is apparent that the simultaneous presence of and interaction between various reagents during synthesis control these properties, computational and experimental approaches for exploring the synthesis space… ▽ More Although gold nanorods have been the subject of much research, the pathways for controlling their shape and thereby their optical properties remain largely heuristically understood. Although it is apparent that the simultaneous presence of and interaction between various reagents during synthesis control these properties, computational and experimental approaches for exploring the synthesis space can be either intractable or too time-consuming in practice. This motivates an alternative approach leveraging the wealth of synthesis information already embedded in the body of scientific literature by developing tools to extract relevant structured data in an automated, high-throughput manner. To that end, we present an approach using the powerful GPT-3 language model to extract structured multi-step seed-mediated growth procedures and outcomes for gold nanorods from unstructured scientific text. GPT-3 prompt completions are fine-tuned to predict synthesis templates in the form of JSON documents from unstructured text input with an overall accuracy of $86\%$. The performance is notable, considering the model is performing simultaneous entity recognition and relation extraction. We present a dataset of 11,644 entities extracted from 1,137 papers, resulting in 268 papers with at least one complete seed-mediated gold nanorod growth procedure and outcome for a total of 332 complete procedures. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Designing transparent conductors using forbidden optical transitions

Authors: Rachel Woods-Robinson, Yihuang Xiong, Jimmy-Xuan Shen, Nicholas Winner, Matthew K. Horton, Mark Asta, Alex M. Ganose, Geoffroy Hautier, Kristin A. Persson

Abstract: Many semiconductors present weak or forbidden transitions at their fundamental band gaps, inducing a widened region of transparency. This occurs in high-performing n-type transparent conductors (TCs) such as Sn-doped In2O3 (ITO), however thus far the presence of forbidden transitions has been neglected in searches for new p-type TCs. To address this, we first compute high-throughput absorption spe… ▽ More Many semiconductors present weak or forbidden transitions at their fundamental band gaps, inducing a widened region of transparency. This occurs in high-performing n-type transparent conductors (TCs) such as Sn-doped In2O3 (ITO), however thus far the presence of forbidden transitions has been neglected in searches for new p-type TCs. To address this, we first compute high-throughput absorption spectra across ~18,000 semiconductors, showing that over half exhibit forbidden or weak optical transitions at their band edges. Next, we demonstrate that compounds with highly localized band edge states are more likely to present forbidden transitions. Lastly, we search this set for p-type and n-type TCs with forbidden or weak transitions. Defect calculations yield unexplored TC candidates such as ambipolar BeSiP2, Zr2SN2 and KSe, p-type BAs, Au2S, and AuCl, and n-type Ba2InGaO5, GaSbO4, and KSbO3, among others. We share our data set via the MPContribs platform, and we recommend that future screenings for optical properties use metrics representative of absorption features rather than band gap alone. △ Less

Submitted 9 April, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Modernist Materials Synthesis: Finding Thermodynamic Shortcuts with Hyperdimensional Chemistry

Authors: James R Neilson, Matthew J McDermott, Kristin A Persson

Abstract: Synthesis remains a challenge for advancing materials science. A key focus of this challenge is how to enable selective synthesis, particularly as it pertains to metastable materials. This perspective addresses the question: how can ``spectator'' elements, such as those found in double ion exchange (metathesis) reactions, enable selective materials synthesis? By observing reaction pathways as they… ▽ More Synthesis remains a challenge for advancing materials science. A key focus of this challenge is how to enable selective synthesis, particularly as it pertains to metastable materials. This perspective addresses the question: how can ``spectator'' elements, such as those found in double ion exchange (metathesis) reactions, enable selective materials synthesis? By observing reaction pathways as they happen (\emph{in situ}) and calculating their energetics using modern computational thermodynamics, we observe transient, crystalline intermediates that suggest that many reactions attain a local thermodynamic equilibrium dictated by local chemical potentials far before achieving a global equilibrium set by the average composition. Using this knowledge, one can thermodynamically ``shortcut'' unfavorable intermediates by including additional elements beyond those of the desired target, providing access to a greater number of intermediates with advantageous energetics and selective phase nucleation. Ultimately, data-driven modeling that unites first-principles approaches with experimental insights will refine the accuracy of emerging predictive retrosynthetic models for complex materials synthesis. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Crystal Toolkit: A Web App Framework to Improve Usability and Accessibility of Materials Science Research Algorithms

Authors: Matthew Horton, Jimmy-Xuan Shen, Jordan Burns, Orion Cohen, François Chabbey, Alex M. Ganose, Rishabh Guha, Patrick Huck, Hamming Howard Li, Matthew McDermott, Joseph Montoya, Guy Moore, Jason Munro, Cody O'Donnell, Colin Ophus, Guido Petretto, Janosh Riebesell, Steven Wetizner, Brook Wander, Donald Winston, Ruoxi Yang, Steven Zeltmann, Anubhav Jain, Kristin A. Persson

Abstract: Crystal Toolkit is an open source tool for viewing, analyzing and transforming crystal structures, molecules and other common forms of materials science data in an interactive way. It is intended to help beginners rapidly develop web-based apps to explore their own data or to help developers make their research algorithms accessible to a broader audience of scientists who might not have any traini… ▽ More Crystal Toolkit is an open source tool for viewing, analyzing and transforming crystal structures, molecules and other common forms of materials science data in an interactive way. It is intended to help beginners rapidly develop web-based apps to explore their own data or to help developers make their research algorithms accessible to a broader audience of scientists who might not have any training in computer programming and who would benefit from graphical interfaces. Crystal Toolkit comes with a library of ready-made components that can be assembled to make complex web apps: simulation of powder and single crystalline diffraction patterns, convex hull phase diagrams, Pourbaix diagrams, electronic band structures, analysis of local chemical environments and symmetry, and more. Crystal Toolkit is now powering the Materials Project website frontend, providing user-friendly access to its database of computed materials properties. In the future, it is hoped that new visualizations might be prototyped using Crystal Toolkit to help explore new forms of data being generated by the materials science community, and that this in turn can help new materials scientists develop intuition for how their data behaves and the insights that might be found within. Crystal Toolkit will remain a work-in-progress and is open to contributions from the community. △ Less

Submitted 27 February, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

High-throughput optical absorption spectra for inorganic semiconductors

Authors: Ruo Xi Yang, Matthew K. Horton, Jason Munro, Kristin A. Persson

Abstract: An optical absorption spectrum constitutes one of the most fundamental material characteristics, with relevant applications ranging from material identification to energy harvesting and optoelectronics. However, the database of both experimental and computational spectra is currently lacking. In this study, we designed a computational workflow for the optical absorption spectrum and integrated the… ▽ More An optical absorption spectrum constitutes one of the most fundamental material characteristics, with relevant applications ranging from material identification to energy harvesting and optoelectronics. However, the database of both experimental and computational spectra is currently lacking. In this study, we designed a computational workflow for the optical absorption spectrum and integrated the simulated spectra into the Materials Project. Using density-functional theory, we computed the frequency-dependent dielectric function and the corresponding absorption coefficient for more than 1000 solid compounds of varying crystal structure and chemistry. The computed spectra show excellent agreement, as quantified by a high value of the Pearson correlation, with experimental results when applying the band gap correction from the HSE functional. The demonstrated calculated accuracy in the spectra suggests that the workflow can be applied in screening studies for materials with specific optical properties. △ Less

Submitted 16 September, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Inferring global dynamics from local structure in liquid electrolytes

Authors: Penelope K. Jones, Kara D. Fong, Kristin A. Persson, Alpha A. Lee

Abstract: Ion transport in concentrated electrolytes plays a fundamental role in electrochemical systems such as lithium ion batteries. Nonetheless, the mechanism of transport amid strong ion-ion interactions remains enigmatic. A key question is whether the dynamics of ion transport can be predicted by the local static structure alone, and if so what are the key structural motifs that determine transport. I… ▽ More Ion transport in concentrated electrolytes plays a fundamental role in electrochemical systems such as lithium ion batteries. Nonetheless, the mechanism of transport amid strong ion-ion interactions remains enigmatic. A key question is whether the dynamics of ion transport can be predicted by the local static structure alone, and if so what are the key structural motifs that determine transport. In this paper, we show that machine learning can successfully decompose global conductivity into the spatio-temporal average of local, instantaneous ionic contributions, and relate this ``local molar conductivity" field to the local ionic environment. Our machine learning model accurately predicts the molar conductivity of electrolyte systems that were not part of the training set, suggesting that the dynamics of ion transport is predictable from local static structure. Further, through analysing this machine-learned local conductivity field, we observe that fluctuations in local conductivity at high concentration are negatively correlated with total molar conductivity. Surprisingly, these fluctuations arise due to a long tail distribution of low conductivity ions, rather than distinct ion pairs, and are spatially correlated through both like- and unlike-charge interactions. More broadly, our approach shows how machine learning can aid the understanding of complex soft matter systems, by learning a function that attributes global collective properties to local, atomistic contributions. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 20 pages, 8 figures

A method to computationally screen for tunable properties of crystalline alloys

Authors: Rachel Woods-Robinson, Matthew K. Horton, Kristin A. Persson

Abstract: Conventionally, high-throughput computational materials searches start from an input set of bulk compounds extracted from material databases, and this set is screened for candidate materials for specific applications. In contrast, many functional materials, and especially semiconductors, are heavily engineered alloys or solid solutions of multiple compounds rather than a single bulk compound. To i… ▽ More Conventionally, high-throughput computational materials searches start from an input set of bulk compounds extracted from material databases, and this set is screened for candidate materials for specific applications. In contrast, many functional materials, and especially semiconductors, are heavily engineered alloys or solid solutions of multiple compounds rather than a single bulk compound. To improve our ability to design functional materials, in this work we propose a framework and open-source code to automatically construct possible "alloy pairs" and "alloy systems" and detect "alloy members" from a set of existing, experimental or calculated ordered compounds, without requiring any additional metadata beyond their crystal structure. We provide analysis tools to estimate stability across each alloy. As a demonstration, we apply this framework to all inorganic materials in the Materials Project database to create a new database of over 600,000 unique alloy pair entries that can then be used in materials discovery studies to search for materials with tunable properties. This new database has been incorporated into the Materials Project website and linked with corresponding material identifiers for any user to query and explore. Using an example of screening for p-type transparent conducting materials, we demonstrate how using this methodology reveals candidate material systems that might otherwise have been excluded by a traditional screening. This work lays a foundation from which materials databases can go beyond stoichiometric compounds, and approach a more realistic description of compositionally tunable materials. △ Less

Submitted 9 April, 2023; v1 submitted 21 June, 2022; originally announced June 2022.

Journal ref: Patterns (2023)

Expanding the Materials Search Space for Multivalent Cathodes

Authors: Ann Rutt, Jimmy-Xuan Shen, Matthew Horton, Jiyoon Kim, Jerry Lin, Kristin A. Persson

Abstract: Multivalent batteries are an energy storage technology with the potential to surpass lithium-ion batteries, however their performance has been limited by the low voltages and poor solid-state ionic mobility of available cathodes. A computational screening approach to identify high-performance multivalent intercalation cathodes among materials that do not contain the working ion of interest has bee… ▽ More Multivalent batteries are an energy storage technology with the potential to surpass lithium-ion batteries, however their performance has been limited by the low voltages and poor solid-state ionic mobility of available cathodes. A computational screening approach to identify high-performance multivalent intercalation cathodes among materials that do not contain the working ion of interest has been developed which greatly expands the search space that can be considered for materials discovery. This approach has been applied to magnesium cathodes as a proof of concept and four resulting candidate materials (NASICON V$_2$(PO$_4$)$_3$, birnessite NaMn$_4$O$_8$, tavorite MnPO$_4$F, and spinel MnO$_2$) are discussed in further detail. In examining the ion migration environment and associated Mg$^{2+}$ migration energy in these materials, local energy maxima are found to correspond with pathway positions where Mg$^{2+}$ passes through a plane of anion atoms. While previous works have established the influence of local coordination on multivalent ion mobility, these results suggest that considering both the type of local bonding environment as well as available free volume for the mobile ion along its migration pathway can be significant for improving solid-state mobility. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Topological Graph-based Analysis of Solid-State Ion Migration

Authors: Jimmy-Xuan Shen, Haoming Howard Li, Ann Rutt, Matthew K. Horton, Kristin A. Persson

Abstract: To accelerate the development of novel ion-conducting materials, we present a general graph-theoretic analysis framework for ion migration in any crystalline structure. The nodes of the graph represent metastable sites of the migrating ion and the edges represent discrete migration events between adjacent sites. Starting from a collection of possible metastable migration sites, the framework assig… ▽ More To accelerate the development of novel ion-conducting materials, we present a general graph-theoretic analysis framework for ion migration in any crystalline structure. The nodes of the graph represent metastable sites of the migrating ion and the edges represent discrete migration events between adjacent sites. Starting from a collection of possible metastable migration sites, the framework assigns a weight to the edges by calculating the individual migration energy barriers between those sites. Connected pathways in the periodic simulation cell corresponding to macroscopic ion migration are identified by searching for the lowest-cost cycle in the periodic migration graph. To exemplify the utility of the framework, we present the automatic analyses of Li migration in different polymorphs of VO(PO4), with the resulting identification of two distinct crystal structures with simple migration pathways demonstrating overall < 300 meV migration barriers. △ Less

Submitted 5 July, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

High-throughput determination of Hubbard U and Hund J values for transition metal oxides via linear response formalism

Authors: Guy C. Moore, Matthew K. Horton, Alexander M. Ganose, Martin Siron, Edward Linscott, David D. O'Regan, Kristin A. Persson

Abstract: DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund's J parameters, and the linear response (LR) methodology has proven to be computation… ▽ More DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund's J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 2000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an atomate workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. To demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-p U value results in a significantly improved agreement between the computed lattice parameters and experiment. △ Less

Submitted 27 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: 18 pages, 6 figures

Microscopic Theory of Magnetic Disorder-Induced Decoherence in Superconducting Nb Films

Authors: Evan Sheridan, Thomas F. Harrelson, Eric Sivonxay, Kristin A. Persson, M. Virginia P. Altoé, Irfan Siddiqi, D. Frank Ogletree, David I. Santiago, Sinéad M. Griffin

Abstract: The performance of superconducting qubits is orders of magnitude below what is expected from theoretical estimates based on the loss tangents of the constituent bulk materials. This has been attributed to the presence of uncontrolled surface oxides formed during fabrication which can introduce defects and impurities that create decoherence channels. Here, we develop an ab initio Shiba theory to in… ▽ More The performance of superconducting qubits is orders of magnitude below what is expected from theoretical estimates based on the loss tangents of the constituent bulk materials. This has been attributed to the presence of uncontrolled surface oxides formed during fabrication which can introduce defects and impurities that create decoherence channels. Here, we develop an ab initio Shiba theory to investigate the microscopic origin of magnetic-induced decoherence in niobium thin film superconductors and the formation of native oxides. Our ab initio calculations encompass the roles of structural disorder, stoichiometry, and strain on the formation of decoherence-inducing local spin moments. With parameters derived from these first-principles calculations we develop an effective quasi-classical model of magnetic-induced losses in the superconductor. We identify d-channel losses (associated with oxygen vacancies) as especially parasitic, resulting in a residual zero temperature surface impedance. This work provides a route to connecting atomic scale properties of superconducting materials and macroscopic decoherence channels affecting quantum systems. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 10 pages plus supplemental materials

Full spectrum optical constant interface to the Materials Project

Authors: J. J. Kas, F. D. Vila, C. D. Pemmaraju, M. P. Prange, K. A. Persson, R. X. Yang, J. J. Rehr

Abstract: Optical constants characterize the interaction of materials with light and are important properties in material design. Here we present a Python-based Corvus workflow for simulations of full spectrum optical constants from the UV-VIS to hard x-ray wavelengths based on the real-space Green's function code FEFF10 and structural data from the Materials Project (MP). The Corvus workflow manager and it… ▽ More Optical constants characterize the interaction of materials with light and are important properties in material design. Here we present a Python-based Corvus workflow for simulations of full spectrum optical constants from the UV-VIS to hard x-ray wavelengths based on the real-space Green's function code FEFF10 and structural data from the Materials Project (MP). The Corvus workflow manager and its associated tools provide an interface to FEFF10 and the MP database. The workflow parallelizes the FEFF computations of optical constants over all absorption edges for each material in the MP database specified by a unique MP-ID. The workflow tools determine the distribution of computational resources needed for that case. Similarly, the optical constants for selected sets of materials can be computed in a single-shot. To illustrate the approach, we present results for nearly all elemental solids in the periodic table, as well as a sample compound, and compared with experimental results. As in x-ray absorption spectra, these results are interpreted in terms of an atomic-like background and fine-structure contributions. △ Less

Submitted 25 August, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

A representation-independent electronic charge density database for crystalline materials

Authors: Jimmy-Xuan Shen, Jason M. Munro, Matthew K. Horton, Patrick Huck, Shyam Dwaraknath, Kristin A. Persson

Abstract: In addition to being the core quantity in density functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With increasing utilization of machine-learning tools in materials sciences, a data-rich object lik… ▽ More In addition to being the core quantity in density functional theory, the charge density can be used in many tertiary analyses in materials sciences from bonding to assigning charge to specific atoms. The charge density is data-rich since it contains information about all the electrons in the system. With increasing utilization of machine-learning tools in materials sciences, a data-rich object like the charge density can be utilized in a wide range of applications. The database presented here provides a modern and user-friendly interface for a large and continuously updated collection of charge densities as part of the Materials Project. In addition to the charge density data, we provide the theory and code for changing the representation of the charge density which should enable more advanced machine-learning studies for the broader community. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Selectivity in yttrium manganese oxide synthesis via local chemical potentials in hyperdimensional phase space

Authors: Paul K. Todd, Matthew J. McDermott, Christopher L. Rom, Adam A. Corrao, Jonathan J. Denney, Shyam S. Dwaraknath, Peter G. Khalifah, Kristin A. Persson, James R. Neilson

Abstract: In sharp contrast to molecular synthesis, materials synthesis is generally presumed to lack selectivity. The few known methods of designing selectivity in solid-state reactions have limited scope, such as topotactic reactions or strain stabilization. This contribution describes a general approach for searching large chemical spaces to identify selective reactions. This novel approach explains the… ▽ More In sharp contrast to molecular synthesis, materials synthesis is generally presumed to lack selectivity. The few known methods of designing selectivity in solid-state reactions have limited scope, such as topotactic reactions or strain stabilization. This contribution describes a general approach for searching large chemical spaces to identify selective reactions. This novel approach explains the ability of a nominally "innocent" Na$_2$CO$_3$ precursor to enable the metathesis synthesis of single-phase Y$_2$Mn$_2$O$_7$ -- an outcome that was previously only accomplished at extreme pressures and which cannot be achieved with closely related precursors of Li$_2$CO$_3$ and K$_2$CO$_3$. By calculating the required change in chemical potential across all possible reactant-product interfaces in an expanded chemical space including Y, Mn, O, alkali metals, and halogens, using thermodynamic parameters obtained from density functional theory calculations, we identify reactions that minimize the thermodynamic competition from intermediates. In this manner, only the Na-based intermediates minimize the distance in the hyperdimensional chemical potential space to Y$_2$Mn$_2$O$_7$, thus providing selective access to a phase which was previously thought to be metastable. Experimental evidence validating this mechanism for pathway-dependent selectivity is provided by intermediates identified from in situ synchrotron-based crystallographic analysis. This approach of calculating chemical potential distances in hyperdimensional compositional spaces provides a general method for designing selective solid-state syntheses that will be useful for gaining access to metastable phases and for identifying reaction pathways that can reduce the synthesis temperature, and cost, of technological materials. △ Less

Submitted 15 August, 2021; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: 30 pages with 5 figures. The first two authors contributed equally to this work

Comparison of the Tetrahedron Method to Smearing Methods for the Electronic Density of States

Authors: M. Y. Toriyama, A. M. Ganose, M. Dylla, S. Anand, J. Park, M. K. Brod, J. Munro, K. A. Persson, A. Jain, G. J. Snyder

Abstract: The electronic density of states (DOS) highlights fundamental properties of materials that oftentimes dictate their properties, such as the band gap and Van Hove singularities. In this short note, we discuss how sharp features of the density of states can be obscured by smearing methods (such as the Gaussian and Fermi smearing methods) when calculating the DOS. While the common approach to reach a… ▽ More The electronic density of states (DOS) highlights fundamental properties of materials that oftentimes dictate their properties, such as the band gap and Van Hove singularities. In this short note, we discuss how sharp features of the density of states can be obscured by smearing methods (such as the Gaussian and Fermi smearing methods) when calculating the DOS. While the common approach to reach a "converged" density of states of a material is to increase the discrete k-point mesh density, we show that the DOS calculated by smearing methods can appear to converge but not to the correct DOS. Employing the tetrahedron method for Brillouin zone integration resolves key features of the density of states far better than smearing methods. △ Less

Submitted 4 March, 2021; originally announced March 2021.

OPTIMADE, an API for exchanging materials data

Authors: Casper W. Andersen, Rickard Armiento, Evgeny Blokhin, Gareth J. Conduit, Shyam Dwaraknath, Matthew L. Evans, Ádám Fekete, Abhijith Gopakumar, Saulius Gražulis, Andrius Merkys, Fawzi Mohamed, Corey Oses, Giovanni Pizzi, Gian-Marco Rignanese, Markus Scheidgen, Leopold Talirz, Cormac Toher, Donald Winston, Rossella Aversa, Kamal Choudhary, Pauline Colinet, Stefano Curtarolo, Davide Di Stefano, Claudia Draxl, Suleyman Er , et al. (31 additional authors not shown)

Abstract: The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API throug… ▽ More The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification. △ Less

Submitted 25 August, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 11 pages, 1 table

Journal ref: Sci Data 8, 217 (2021)

The role of disorder in the synthesis of metastable zinc zirconium nitrides

Authors: Rachel Woods-Robinson, Vladan Stevanović, Stephan Lany, Karen N. Heinselman, Matthew K. Horton, Kristin A. Persson, Andriy Zakutayev

Abstract: In materials science, it is often assumed that ground state crystal structures predicted by density functional theory are the easiest polymorphs to synthesize. Ternary nitride materials, with many possible metastable polymorphs, provide a rich materials space to study what influences thermodynamic stability and polymorph synthesizability. For example, ZnZrN2 is theoretically predicted at zero Kelv… ▽ More In materials science, it is often assumed that ground state crystal structures predicted by density functional theory are the easiest polymorphs to synthesize. Ternary nitride materials, with many possible metastable polymorphs, provide a rich materials space to study what influences thermodynamic stability and polymorph synthesizability. For example, ZnZrN2 is theoretically predicted at zero Kelvin to have an unusual layered "wurtsalt" ground state crystal structure with compelling optoelectronic properties, but it is unknown whether this structure can be realized experimentally under practical synthesis conditions. Here, we use combinatorial sputtering to synthesize hundreds of ZnxZr1-xNy thin film samples, and find metastable rocksalt-derived or boron-nitride-derived structures rather than the predicted wurtsalt structure. Using a statistical polymorph sampler approach, it is demonstrated that although rocksalt is the least stable polymorph at zero Kelvin, it becomes the most stable polymorph at high effective temperatures similar to those achieved using this sputter deposition method, and thus corroborates experimental results. Additional calculations show that this destabilization of the wurtsalt polymorph is due to configurational entropic and enthalpic effects, and that vibrational contributions are negligible. Specifically, rocksalt- and boron-nitride-derived structures become the most stable polymorphs in the presence of disorder because of higher tolerances to cation cross-substitution and off-stoichiometry than the wurtsalt structure. This understanding of the role of disorder tolerance in the synthesis of competing polymorphs can enable more accurate predictions of synthesizable crystal structures and their achievable material properties. △ Less

Submitted 11 March, 2022; v1 submitted 22 December, 2020; originally announced December 2020.

Efficient calculation of carrier scattering rates from first principles

Authors: Alex M. Ganose, Junsoo Park, Alireza Faghaninia, Rachel Woods-Robinson, Kristin A. Persson, Anubhav Jain

Abstract: The electronic transport behaviour of materials determines their suitability for technological applications. We develop an efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms form… ▽ More The electronic transport behaviour of materials determines their suitability for technological applications. We develop an efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 16 semiconductors and comparing the results against experimental measurements. The present work is amenable for use in high-throughput computational workflows and enables accurate screening of carrier mobilities, lifetimes, and thermoelectric power. △ Less

Submitted 22 December, 2020; v1 submitted 21 August, 2020; originally announced August 2020.

Comments: 11 pages, 4 figures (SI: 21 pages, 14 figures)

High-throughput search for magnetic and topological order in transition metal oxides

Authors: Nathan C. Frey, Matthew K. Horton, Jason M. Munro, Sinéad M. Griffin, Kristin A. Persson, Vivek B. Shenoy

Abstract: The discovery of intrinsic magnetic topological order in $\rm MnBi_2Te_4$ has invigorated the search for materials with coexisting magnetic and topological phases. These multi-order quantum materials are expected to exhibit new topological phases that can be tuned with magnetic fields, but the search for such materials is stymied by difficulties in predicting magnetic structure and stability. Here… ▽ More The discovery of intrinsic magnetic topological order in $\rm MnBi_2Te_4$ has invigorated the search for materials with coexisting magnetic and topological phases. These multi-order quantum materials are expected to exhibit new topological phases that can be tuned with magnetic fields, but the search for such materials is stymied by difficulties in predicting magnetic structure and stability. Here, we compute over 27,000 unique magnetic orderings for over 3,000 transition metal oxides in the Materials Project database to determine their magnetic ground states and estimate their effective exchange parameters and critical temperatures. We perform a high-throughput band topology analysis of centrosymmetric magnetic materials, calculate topological invariants, and identify 18 new candidate ferromagnetic topological semimetals, axion insulators, and antiferromagnetic topological insulators. To accelerate future efforts, machine learning classifiers are trained to predict both magnetic ground states and magnetic topological order without requiring first-principles calculations. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Comments: 9 pages, 6 figures

Journal ref: Science Advances 09 Dec 2020: Vol. 6, no. 50, eabd1076

Incorporating electronic information into Machine Learning potential energy surfaces via approaching the ground-state electronic energy as a function of atom-based electronic populations

Authors: Xiaowei Xie, Kristin A. Persson, David W. Small

Abstract: Machine Learning (ML) approximations to Density Functional Theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the med… ▽ More Machine Learning (ML) approximations to Density Functional Theory (DFT) potential energy surfaces (PESs) are showing great promise for reducing the computational cost of accurate molecular simulations, but at present they are not applicable to varying electronic states, and in particular, they are not well suited for molecular systems in which the local electronic structure is sensitive to the medium to long-range electronic environment. With this issue as the focal point, we present a new Machine Learning approach called bpopNN for obtaining efficient approximations to DFT PESs. The methodology is based on approaching the true DFT energy as a function of electron populations on atoms, which may be realized in practice with constrained DFT (CDFT). The new approach creates approximations to this function with deep neural networks. These approximations thereby incorporate electronic information naturally into a ML approach, and optimizing the model energy with respect to populations allows the electronic terms to self-consistently adapt to the environment, as in DFT. We confirm the effectiveness of this approach with a variety of calculations on Li$_n$H$_n$ clusters. △ Less

Submitted 4 March, 2020; originally announced March 2020.

Comments: 28 pages, 8 figures

Two-Dimensional Forms of Robust CO$_2$ Reduction Photocatalysts

Authors: Steven B. Torrisi, Arunima K. Singh, Joseph H. Montoya, Tathagata Biswas, Kristin A. Persson

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

Submitted 19 December, 2019; originally announced December 2019.

Comments: 23 Pages, 5 figures

Efficient Pourbaix diagrams of many-element compounds

Authors: Anjli Patel, Jens K. Nørskov, Kristin A. Persson, Joseph H. Montoya

Abstract: Pourbaix diagrams have long been an essential tool for determining the phase stability of solids and their associated ionic species under electrochemical conditions. In recent years, Pourbaix diagrams have been used for applications ranging from corrosion-resistance alloy design to electrocatalysis, and the data from which they are generated has been enhanced by the availability of materials data… ▽ More Pourbaix diagrams have long been an essential tool for determining the phase stability of solids and their associated ionic species under electrochemical conditions. In recent years, Pourbaix diagrams have been used for applications ranging from corrosion-resistance alloy design to electrocatalysis, and the data from which they are generated has been enhanced by the availability of materials data in various online databases. However, generation of multi-element Pourbaix diagrams has a critical bottleneck which makes 3-element systems difficult to analyze quickly, and 4 and 5 element systems intractable. In this work, we present a method for constructing Pourbaix diagrams which uses a pre-processing step to circumvent the most egregious computational bottleneck and make many-element Pourbaix diagrams computationally efficient. △ Less

Submitted 30 August, 2019; originally announced September 2019.

Grain Boundary Properties of Elemental Metals

Authors: Hui Zheng, Xiang-Guo Li, Richard Tran, Chi Chen, Matthew Horton, Donny Winston, Kristin Aslaug Persson, Shyue Ping Ong

Abstract: The structure and energy of grain boundaries (GBs) are essential for predicting the properties of polycrystalline materials. In this work, we use high-throughput density functional theory calculations workflow to construct the Grain Boundary Database (GBDB), the largest database of DFT-computed grain boundary properties to date. The database currently encompasses 327 GBs of 58 elemental metals, in… ▽ More The structure and energy of grain boundaries (GBs) are essential for predicting the properties of polycrystalline materials. In this work, we use high-throughput density functional theory calculations workflow to construct the Grain Boundary Database (GBDB), the largest database of DFT-computed grain boundary properties to date. The database currently encompasses 327 GBs of 58 elemental metals, including 10 common twist or symmetric tilt GBs for body-centered cubic (bcc) and face-centered cubic (fcc) systems and the $Σしぐま$7 [0001] twist GB for hexagonal close-packed (hcp) systems. In particular, we demonstrate a novel scaled-structural template approach for HT GB calculations, which reduces the computational cost of converging GB structures by a factor of $\sim 3-6$. The grain boundary energies and work of separation are rigorously validated against previous experimental and computational data. Using this large GB dataset, we develop an improved predictive model for the GB energy of different elements based on the cohesive energy and shear modulus. The open GBDB represent a significant step forward in the availability of first principles GB properties, which we believe would help guide the future design of polycrystalline materials. △ Less

Submitted 20 July, 2019; originally announced July 2019.

Comments: 32 pages, 8 figures, 2 tables in main manuscript; 6 figures, 5 tables, in 14 pages of Supplementary Information

Pawpyseed: Perturbation-extrapolation band shifting corrections for point defect calculations

Authors: Kyle Bystrom, Danny Broberg, Shyam Dwaraknath, Kristin A. Persson, Mark Asta

Abstract: Significant progress has been made recently in the automation and standardization of ab initio point defect calculations. However, the task of developing, implementing, and benchmarking charge corrections for density functional theory (DFT) point defect calculations is still an open challenge. Here we present a high-performance Python package called pawpyseed, which can read PAW DFT wave functions… ▽ More Significant progress has been made recently in the automation and standardization of ab initio point defect calculations. However, the task of developing, implementing, and benchmarking charge corrections for density functional theory (DFT) point defect calculations is still an open challenge. Here we present a high-performance Python package called pawpyseed, which can read PAW DFT wave functions and calculate the overlap between wavefunctions from different structures. Using pawpyseed, we implement a new band shifting correction derived from first order perturbation theory. We benchmark this method by calculating the transition levels of several point defects in silicon and comparing to experimental and hybrid functional results. The new band shifting method can shift single-particle energies to improve transition level predictions and can be automated and parallelized using pawpyseed, suggesting it could be a useful method for high-throughput point defect calculations. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Anisotropic work function of elemental crystals

Authors: Richard Tran, Xiang-Guo Li, Joseph Montoya, Donald Winston, Kristin Aslaug Persson, Shyue Ping Ong

Abstract: The work function is a fundamental electronic property of a solid that varies with the facets of a crystalline surface. It is a crucial parameter in spectroscopy as well as materials design, especially for technologies such as thermionic electron guns and Schottky barriers. In this work, we present the largest database of calculated work functions for elemental crystals to date. This database cont… ▽ More The work function is a fundamental electronic property of a solid that varies with the facets of a crystalline surface. It is a crucial parameter in spectroscopy as well as materials design, especially for technologies such as thermionic electron guns and Schottky barriers. In this work, we present the largest database of calculated work functions for elemental crystals to date. This database contains the anisotropic work functions of more than 100 polymorphs of about 72 elements and up to a maximum Miller index of two and three for non-cubic and cubic crystals, respectively. The database has been rigorously validated against previous experimental and computational data where available. We also propose a weighted work function based on the Wulff shape that can be compared to measurements from polycrystalline specimens, and show that this weighted work function can be modeled empirically using simple atomic parameters. Furthermore, for the first time, we were able to analyze simple bond breaking rules for metallic systems beyond a maximum Miller index of one, allowing for a more generalized investigation of work function anisotropy. △ Less

Submitted 20 February, 2019; originally announced February 2019.

2DMatPedia: An open computational database of two-dimensional materials from top-down and bottom-up approaches

Authors: Jun Zhou, Lei Shen, Miguel Dias Costa, Kristin A. Persson, Shyue Ping Ong, Patrick Huck, Yunhao Lu, Xiaoyang Ma, Yuan Ping Feng

Abstract: Two-dimensional (2D) materials have been a hot research topic in the last decade, due to novel fundamental physics in the reduced dimension and appealing applications. Systematic discovery of functional 2D materials has been the focus of many studies. Here, we present a large dataset of 2D materials, with more than 6,000 monolayer structures, obtained from both top-down and bottom-up discovery pro… ▽ More Two-dimensional (2D) materials have been a hot research topic in the last decade, due to novel fundamental physics in the reduced dimension and appealing applications. Systematic discovery of functional 2D materials has been the focus of many studies. Here, we present a large dataset of 2D materials, with more than 6,000 monolayer structures, obtained from both top-down and bottom-up discovery procedures. First, we screened all bulk materials in the database of Materials Project for layered structures by a topology-based algorithm, and theoretically exfoliate them into monolayers. Then, we generated new 2D materials by chemical substitution of elements in known 2D materials by others from the same group in the periodic table. The structural, electronic and energetic properties of these 2D materials are consistently calculated, to provide a starting point for further material screening, data mining, data analysis and artificial intelligence applications. We present the details of computational methodology, data record and technical validation of our publicly available data (http://www.2dmatpedia.org/). △ Less

Submitted 27 January, 2019; originally announced January 2019.

Elucidating the Structure of the Magnesium Aluminum Chloride Complex electrolyte for Magnesium-ion batteries

Authors: Pieremanuele Canepa, Saivenkataraman Jayaraman, Lei Cheng, Nav Nidhi Rajput, William D. Richards, Gopalakrishnan Sai Gautam, Larry A. Curtiss, Kristin A. Persson, Gerbrand Ceder

Abstract: We present a rigorous analysis of the Magnesium Aluminum Chloro Complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use \emph{ab initio} calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: \emph{i})… ▽ More We present a rigorous analysis of the Magnesium Aluminum Chloro Complex (MACC) in tetrahydrofuran (THF), one of the few electrolytes that can reversibly plate and strip Mg. We use \emph{ab initio} calculations and classical molecular dynamics simulations to interrogate the MACC electrolyte composition with the goal of addressing two urgent questions that have puzzled battery researchers: \emph{i}) the functional species of the electrolyte, and \emph{ii}) the complex equilibria regulating the MACC speciation after prolonged electrochemical cycling, a process termed as conditioning, and after prolonged inactivity, a process called aging. A general computational strategy to untangle the complex structure of electrolytes, ionic liquids and other liquid media is presented. The analysis of formation energies and grand-potential phase diagrams of Mg-Al-Cl-THF suggests that the MACC electrolyte bears a simple chemical structure with few simple constituents, namely the electro-active species MgCl$^+$ and AlCl$_4^-$ in equilibrium with MgCl$_2$ and AlCl$_3$. Knowledge of the stable species of the MACC electrolyte allows us to determine the most important equilibria occurring during electrochemical cycling. We observe that Al deposition is always preferred to Mg deposition, explaining why freshly synthesized MACC cannot operate and needs to undergo preparatory conditioning. Similarly, we suggest that aluminum displacement and depletion from the solution upon electrolyte resting (along with continuous MgCl$_2$ regeneration) represents one of the causes of electrolyte aging. Finally, we compute the NMR shifts from shielding tensors of selected molecules and ions providing fingerprints to guide future experimental investigations. △ Less

Submitted 8 November, 2015; originally announced November 2015.

Journal ref: Energy Environ. Sci. 8, 3718-3730 (2015)