-
MatterSim: A Deep Learning Atomistic Model Across Elements, Temperatures and Pressures
Authors:
Han Yang,
Chenxi Hu,
Yichi Zhou,
Xixian Liu,
Yu Shi,
Jielan Li,
Guanzhi Li,
Zekun Chen,
Shuizhou Chen,
Claudio Zeni,
Matthew Horton,
Robert Pinsler,
Andrew Fowler,
Daniel Zügner,
Tian Xie,
Jake Smith,
Lixin Sun,
Qian Wang,
Lingyu Kong,
Chang Liu,
Hongxia Hao,
Ziheng Lu
Abstract:
Accurate and fast prediction of materials properties is central to the digital transformation of materials design. However, the vast design space and diverse operating conditions pose significant challenges for accurately modeling arbitrary material candidates and forecasting their properties. We present MatterSim, a deep learning model actively learned from large-scale first-principles computatio…
▽ More
Accurate and fast prediction of materials properties is central to the digital transformation of materials design. However, the vast design space and diverse operating conditions pose significant challenges for accurately modeling arbitrary material candidates and forecasting their properties. We present MatterSim, a deep learning model actively learned from large-scale first-principles computations, for efficient atomistic simulations at first-principles level and accurate prediction of broad material properties across the periodic table, spanning temperatures from 0 to 5000 K and pressures up to 1000 GPa. Out-of-the-box, the model serves as a machine learning force field, and shows remarkable capabilities not only in predicting ground-state material structures and energetics, but also in simulating their behavior under realistic temperatures and pressures, signifying an up to ten-fold enhancement in precision compared to the prior best-in-class. This enables MatterSim to compute materials' lattice dynamics, mechanical and thermodynamic properties, and beyond, to an accuracy comparable with first-principles methods. Specifically, MatterSim predicts Gibbs free energies for a wide range of inorganic solids with near-first-principles accuracy and achieves a 15 meV/atom resolution for temperatures up to 1000K compared with experiments. This opens an opportunity to predict experimental phase diagrams of materials at minimal computational cost. Moreover, MatterSim also serves as a platform for continuous learning and customization by integrating domain-specific data. The model can be fine-tuned for atomistic simulations at a desired level of theory or for direct structure-to-property predictions, achieving high data efficiency with a reduction in data requirements by up to 97%.
△ Less
Submitted 10 May, 2024; v1 submitted 8 May, 2024;
originally announced May 2024.
-
MatterGen: a generative model for inorganic materials design
Authors:
Claudio Zeni,
Robert Pinsler,
Daniel Zügner,
Andrew Fowler,
Matthew Horton,
Xiang Fu,
Sasha Shysheya,
Jonathan Crabbé,
Lixin Sun,
Jake Smith,
Bichlien Nguyen,
Hannes Schulz,
Sarah Lewis,
Chin-Wei Huang,
Ziheng Lu,
Yichi Zhou,
Han Yang,
Hongxia Hao,
Jielan Li,
Ryota Tomioka,
Tian Xie
Abstract:
The design of functional materials with desired properties is essential in driving technological advances in areas like energy storage, catalysis, and carbon capture. Generative models provide a new paradigm for materials design by directly generating entirely novel materials given desired property constraints. Despite recent progress, current generative models have low success rate in proposing s…
▽ More
The design of functional materials with desired properties is essential in driving technological advances in areas like energy storage, catalysis, and carbon capture. Generative models provide a new paradigm for materials design by directly generating entirely novel materials given desired property constraints. Despite recent progress, current generative models have low success rate in proposing stable crystals, or can only satisfy a very limited set of property constraints. Here, we present MatterGen, a model that generates stable, diverse inorganic materials across the periodic table and can further be fine-tuned to steer the generation towards a broad range of property constraints. To enable this, we introduce a new diffusion-based generative process that produces crystalline structures by gradually refining atom types, coordinates, and the periodic lattice. We further introduce adapter modules to enable fine-tuning towards any given property constraints with a labeled dataset. Compared to prior generative models, structures produced by MatterGen are more than twice as likely to be novel and stable, and more than 15 times closer to the local energy minimum. After fine-tuning, MatterGen successfully generates stable, novel materials with desired chemistry, symmetry, as well as mechanical, electronic and magnetic properties. Finally, we demonstrate multi-property materials design capabilities by proposing structures that have both high magnetic density and a chemical composition with low supply-chain risk. We believe that the quality of generated materials and the breadth of MatterGen's capabilities represent a major advancement towards creating a universal generative model for materials design.
△ Less
Submitted 29 January, 2024; v1 submitted 6 December, 2023;
originally announced December 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.
-
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.
-
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.
-
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.
-
Vapor-liquid-solid growth of highly-mismatched semiconductor nanowires with high-fidelity van der Waals layer stacking
Authors:
Edy Cardona,
Matthew K. Horton,
Daniel Paulo-Wach,
Anthony C. Salazar,
Andre Palacios Duran,
James Chavez,
Shaul Aloni,
Junqiao Wu,
Oscar D. Dubon
Abstract:
Nanobelts, nanoribbons and other quasi-one-dimensional nanostructures formed from layered, so-called, van der Waals semiconductors have garnered much attention due to their high-performance, tunable optoelectronic properties. For layered alloys made from the gallium monochalcogenides GaS, GaSe, and GaTe, near-continuous tuning of the energy bandgap across the full composition range has been achiev…
▽ More
Nanobelts, nanoribbons and other quasi-one-dimensional nanostructures formed from layered, so-called, van der Waals semiconductors have garnered much attention due to their high-performance, tunable optoelectronic properties. For layered alloys made from the gallium monochalcogenides GaS, GaSe, and GaTe, near-continuous tuning of the energy bandgap across the full composition range has been achieved in GaSe1-xSx and GaSe1-xTex alloys. Gold-catalyzed vapor-liquid-solid (VLS) growth of these alloys yields predominantly nanobelts, nanoribbons and other nanostructures for which the fast crystal growth front consists of layer edges in contact with the catalyst. We demonstrate that in the S-rich, GaS1-xTex system, unlike GaSe1-xSx and GaSe1-xTex, the Au-catalyzed VLS process yields van der Waals nanowires for which the fast growth direction is normal to the layers. The high mismatch between S and Te leads to extraordinary bowing of the GaS1-xTex alloy's energy bandgap, decreasing by at least 0.6 eV for x as small as 0.03. Calculations using density functional theory confirm the significant decrease in bandgap in S-rich GaS1-xTex. The nanowires can exceed fifty micrometers in length, consisting of tens of thousands of van der Waals-bonded layers with triangular or hexagonal cross-sections of uniform dimensions along the length of the nanowire. We propose that the low solubility of Te in GaS results in an enhancement in Te coverage around the Au catalyst-nanowire interface, confining the catalyst to the chalcogen-terminated basal plane (rather than the edges) and thereby enabling layer-by-layer, c-axis growth.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
Influential Functionals
Authors:
M. K. Horton
Abstract:
Learning about density functional approximations (DFAs), or approximations for the exchange-correlation functional, can be intimidating. Density Functional Theory is now one of the primary simulation tools for the practicing chemist or materials scientist, and its accuracy relies upon an appropriate choice of DFA. Over the past decades, there has been extensive research effort to find better DFAs,…
▽ More
Learning about density functional approximations (DFAs), or approximations for the exchange-correlation functional, can be intimidating. Density Functional Theory is now one of the primary simulation tools for the practicing chemist or materials scientist, and its accuracy relies upon an appropriate choice of DFA. Over the past decades, there has been extensive research effort to find better DFAs, and there is now a large body of literature to read through for someone learning about DFAs for the first time. In this brief report, I share an analysis that suggests which functionals and publications have been the most influential, as a potential reading list to new scientists in this area. Here, "influential" is defined as "likely to have informed the design of another functional," and not simply a measure of number of citations, or how much that functional has been used for practical applications. This analysis is not claimed to be complete.
△ Less
Submitted 2 November, 2022;
originally announced November 2022.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
New Methods for Calculating the Free Energy of Charged Defects in Solid Electrolytes
Authors:
Robert M. Horton,
Andrew J. Haslam,
Amparo Galindo,
George Jackson,
Michael W. Finnis
Abstract:
A methodology for calculating the contribution of charged defects to the configurational free energy of an ionic crystal is introduced. The temperature-independent Wang-Landau Monte Carlo technique is applied to a simple model of a solid electrolyte, consisting of charged positive and negative defects on a lattice. The electrostatic energy is computed on lattices with periodic boundary conditions,…
▽ More
A methodology for calculating the contribution of charged defects to the configurational free energy of an ionic crystal is introduced. The temperature-independent Wang-Landau Monte Carlo technique is applied to a simple model of a solid electrolyte, consisting of charged positive and negative defects on a lattice. The electrostatic energy is computed on lattices with periodic boundary conditions, and used to calculate the density of states and statistical-thermodynamic potentials of this system. The free energy as a function of defect concentration and temperature is accurately described by a regular solution model up to concentrations of 10% of defects, well beyond the range described by the ideal solution theory.
The approach, supplemented by short-ranged terms in the energy, is proposed as an alternative to free-energy methods that require a number of simulations to be carried out over a range of temperatures.
△ Less
Submitted 23 May, 2013;
originally announced May 2013.
-
Development of anomalous diffusion among crowding proteins
Authors:
Margaret R. Horton,
Felix Höfling,
Joachim O. Rädler,
Thomas Franosch
Abstract:
In cell membranes, proteins and lipids diffuse in a highly crowded and heterogeneous landscape, where aggregates and dense domains of proteins or lipids obstruct the path of diffusing molecules. In general, hindered motion gives rise to anomalous transport, though the nature of the onset of this behavior is still under debate and difficult to investigate experimentally. Here, we present a systemat…
▽ More
In cell membranes, proteins and lipids diffuse in a highly crowded and heterogeneous landscape, where aggregates and dense domains of proteins or lipids obstruct the path of diffusing molecules. In general, hindered motion gives rise to anomalous transport, though the nature of the onset of this behavior is still under debate and difficult to investigate experimentally. Here, we present a systematic study where proteins bound to supported lipid membranes diffuse freely in two dimensions, but are increasingly hindered by the presence of other like proteins. In our model system, the surface coverage of the protein avidin on the lipid bilayer is well controlled by varying the concentration of biotinylated lipid anchors. Using fluorescence correlation spectroscopy (FCS), we measure the time correlation function over long times and convert it to the mean-square displacement of the diffusing proteins. Our approach allows for high precision data and a clear distinction between anomalous and normal diffusion. It enables us to investigate the onset of anomalous diffusion, which takes place when the area coverage of membrane proteins increases beyond approximately 5%. This transition region exhibits pronounced spatial heterogeneities. Increasing the packing fraction further, transport becomes more and more anomalous, manifested in a decrease of the exponent of subdiffusion.
△ Less
Submitted 19 March, 2010;
originally announced March 2010.
-
Nanomechanical detection of antibiotic-mucopeptide binding in a model for superbug drug resistance
Authors:
J. W. Ndieyira,
M. Watari,
A. Donoso Barrera,
D. Zhou,
M. Vögtli,
M. Batchelor,
M. A. Cooper,
T. Strunz,
M. A. Horton,
C. Abell,
T. Rayment,
G. Aeppli,
R. A. Mckendry
Abstract:
The alarming growth of the antibiotic-resistant superbugs methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus (VRE) is driving the development of new technologies to investigate antibiotics and their modes of action. We report the label-free detection of vancomycin binding to bacterial cell wall precursor analogues (mucopeptides) on cantilever arrays, with 10…
▽ More
The alarming growth of the antibiotic-resistant superbugs methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant Enterococcus (VRE) is driving the development of new technologies to investigate antibiotics and their modes of action. We report the label-free detection of vancomycin binding to bacterial cell wall precursor analogues (mucopeptides) on cantilever arrays, with 10 nM sensitivity and at clinically relevant concentrations in blood serum. Differential measurements quantified binding constants for vancomycin-sensitive and vancomycin-resistant mucopeptide analogues. Moreover, by systematically modifying the mucopeptide density we gain new insights into the origin of surface stress. We propose that stress is a product of a local chemical binding factor and a geometrical factor describing the mechanical connectivity of regions affected by local binding in terms of a percolation process. Our findings place BioMEMS devices in a new class of percolative systems. The percolation concept will underpin the design of devices and coatings to significantly lower the drug detection limit and may also impact on our understanding of antibiotic drug action in bacteria.
△ Less
Submitted 15 December, 2008;
originally announced December 2008.
-
Reply to the "Comment on 'Neel order in doped quasi one-dimensional antiferromagnets'"
Authors:
Sebastian Eggert,
Ian Affleck,
Matthew D. P. Horton
Abstract:
In a recent comment by A.A. Zvyagin to our Letter [Phys. Rev. Lett. 89, 47202 (2002)] it was pointed out that the bosonization treatment of the antiferromagnetic Heisenberg chain is only valid in the low temperature and long length limit. We support this statement and quantify the limits of validity by explicitly showing the crossover to the high temperature and/or short length limit using a com…
▽ More
In a recent comment by A.A. Zvyagin to our Letter [Phys. Rev. Lett. 89, 47202 (2002)] it was pointed out that the bosonization treatment of the antiferromagnetic Heisenberg chain is only valid in the low temperature and long length limit. We support this statement and quantify the limits of validity by explicitly showing the crossover to the high temperature and/or short length limit using a combination of bosonization formulas and numerical results.
△ Less
Submitted 30 September, 2002;
originally announced October 2002.
-
Neel order in doped quasi one-dimensional antiferromagnets
Authors:
Sebastian Eggert,
Ian Affleck,
Matthew D. P. Horton
Abstract:
We study the Neel temperature of quasi one-dimensional S=1/2 antiferromagnets containing non-magnetic impurities. We first consider the temperature dependence of the staggered susceptibility of finite chains with open boundary conditions, which shows an interesting difference for even and odd length chains. We then use a mean field theory treatment to incorporate the three dimensional inter-chai…
▽ More
We study the Neel temperature of quasi one-dimensional S=1/2 antiferromagnets containing non-magnetic impurities. We first consider the temperature dependence of the staggered susceptibility of finite chains with open boundary conditions, which shows an interesting difference for even and odd length chains. We then use a mean field theory treatment to incorporate the three dimensional inter-chain couplings. The resulting Neel temperature shows a pronounced drop as a function of doping by up to a factor of 5.
△ Less
Submitted 29 March, 2002;
originally announced April 2002.
-
The Three-Magnon Contribution to the Spin Correlation Function in Integer-Spin Antiferromagnetic Chains
Authors:
Matthew Horton,
Ian Affleck
Abstract:
The exact form factor for the O(3) non-linear sigma model is used to predict the three-magnon contribution to the spin correlation function, S(q,w), near wavevector q=pi in an integer spin, one-dimensional antiferromagnet. The three-magnon contribution is extrememly broad and extremely weak; the integrated intensity is <2% of the single-magnon contribution.
The exact form factor for the O(3) non-linear sigma model is used to predict the three-magnon contribution to the spin correlation function, S(q,w), near wavevector q=pi in an integer spin, one-dimensional antiferromagnet. The three-magnon contribution is extrememly broad and extremely weak; the integrated intensity is <2% of the single-magnon contribution.
△ Less
Submitted 27 July, 1999;
originally announced July 1999.