No Arabic abstract
A central challenge in high throughput density functional theory (HT-DFT) calculations is selecting a combination of input parameters and post-processing techniques that can be used across all materials classes, while also managing accuracy-cost tradeoffs. To investigate the effects of these parameter choices, we consolidate three large HT-DFT databases: Automatic-FLOW (AFLOW), the Materials Project (MP), and the Open Quantum Materials Database (OQMD), and compare reported properties across each pair of databases for materials calculated using the same initial crystal structure. We find that HT-DFT formation energies and volumes are generally more reproducible than band gaps and total magnetizations; for instance, a notable fraction of records disagree on whether a material is metallic (up to 7%) or magnetic (up to 15%). The variance between calculated properties is as high as 0.105 eV/atom (median relative absolute difference, or MRAD, of 6%) for formation energy, 0.65 {AA}$^3$/atom (MRAD of 4%) for volume, 0.21 eV (MRAD of 9%) for band gap, and 0.15 $mu_{rm B}$/formula unit (MRAD of 8%) for total magnetization, comparable to the differences between DFT and experiment. We trace some of the larger discrepancies to choices involving pseudopotentials, the DFT+U formalism, and elemental reference states, and argue that further standardization of HT-DFT would be beneficial to reproducibility.
In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory Tool Kit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational software. It includes task submissions on all major operating systems and task execution on high-performance computing environments. The distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al3Ni, AlNi, AlNi3, Al3Ni4, and Al3Ni5, and the fcc solution phases treated using the special quasirandom structures at the compositions of Al3Ni, AlNi, and AlNi3.
The discoveries of intrinsically magnetic topological materials, including semimetals with a large anomalous Hall effect and axion insulators, have directed fundamental research in solid-state materials. Topological quantum chemistry has enabled the understanding of and the search for paramagnetic topological materials. Using magnetic topological indices obtained from magnetic topological quantum chemistry (MTQC), here we perform a high-throughput search for magnetic topological materials based on first-principles calculations. We use as our starting point the Magnetic Materials Database on the Bilbao Crystallographic Server, which contains more than 549 magnetic compounds with magnetic structures deduced from neutron-scattering experiments, and identify 130 enforced semimetals (for which the band crossings are implied by symmetry eigenvalues), and topological insulators. For each compound, we perform complete electronic structure calculations, which include complete topological phase diagrams using different values of the Hubbard potential. Using a custom code to find the magnetic co-representations of all bands in all magnetic space groups, we generate data to be fed into the algorithm of MTQC to determine the topology of each magnetic material. Several of these materials display previously unknown topological phases, including symmetry-indicated magnetic semimetals, three-dimensional anomalous Hall insulators and higher-order magnetic semimetals. We analyse topological trends in the materials under varying interactions: 60 per cent of the 130 topological materials have topologies sensitive to interactions, and the others have stable topologies under varying interactions. We provide a materials database for future experimental studies and open-source code for diagnosing topologies of magnetic materials.
This paper investigates some of the successes and failures of density functional theory in the study of high-pressure solid hydrogen at low temperature. We calculate the phase diagram, metallization pressure, phonon spectrum, and proton zero-point energy using three popular exchange-correlation functionals: the local density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation, and the semi-local Becke-Lee-Yang-Parr (BLYP) functional. We focus on the solid molecular P$6_3$/m, C2/c, Cmca-12, and Cmca structures in the pressure range from $100<P<500$ GPa over which phases I, II and III are observed experimentally. At the static level of theory, in which proton zero-point energy is ignored, the LDA, PBE and BLYP functionals give very different structural transition and metallization pressures, with the BLYP phase diagram in better agreement with experiment. Nevertheless, all three functionals provide qualitatively the same information about the band gaps of the four structures and the phase transitions between them. Going beyond the static level, we find that the frequencies of the vibron modes observed above 3000 cm$^{-1}$ depend strongly on the choice of exchange-correlation functional, although the low-frequency part of the phonon spectrum is little affected. The largest and smallest values of the proton zero-point energy, obtained using the BLYP and LDA functionals, respectively, differ by more than 10 meV/proton. Including the proton zero-point energy calculated from the phonon spectrum within the harmonic approximation improves the agreement of the BLYP and PBE phase diagrams with experiment. Taken as a whole, our results demonstrate the inadequacy of mean-field-like density functional calculations of solid molecular hydrogen in phases I, II and III and emphasize the need for more sophisticated methods.
High-entropy alloys (HEAs), which have been intensely studied due to their excellent mechanical properties, generally refer to alloys with multiple equimolar or nearly equimolar elements. According to this definition, Si-Ge-Sn alloys with equal or comparable concentrations of the three Group IV elements belong to the category of HEAs. As a result, the equimolar elements of Si-Ge-Sn alloys likely cause their atomic structures to exhibit the same core effects of metallic HEAs such as lattice distortion. Here we apply density functional theory (DFT) calculations to show that the SiGeSn HEA indeed exhibits a large local distortion effect. Unlike metallic HEAs, our Monte Carlo and DFT calculations show that the SiGeSn HEA exhibits no chemical short-range order due to the similar electronegativity of the constituent elements, thereby increasing the configurational entropy of the SiGeSn HEA. Hybrid density functional calculations show that the SiGeSn HEA remains semiconducting with a band gap of 0.38 eV, promising for economical and compatible mid-infrared optoelectronics applications. We then study the energetics of neutral single Si, Ge, and Sn vacancies and (expectedly) find wide distributions of vacancy formation energies, similar to those found in metallic HEAs. However, we also find anomalously small lower bounds (e.g., 0.04 eV for a Si vacancy) in the energy distributions, which arise from the bond reformation near the vacancy. Such small vacancy formation energies and their associated bond reformations retain the semiconducting behavior of the SiGeSn HEA, which may be a signature feature of a semiconducting HEA that differentiates from metallic HEAs.
The traditional paradigm for materials discovery has been recently expanded to incorporate substantial data driven research. With the intent to accelerate the development and the deployment of new technologies, the AFLOW Fleet for computational materials design automates high-throughput first principles calculations, and provides tools for data verification and dissemination for a broad community of users. AFLOW incorporates different computational modules to robustly determine thermodynamic stability, electronic band structures, vibrational dispersions, thermo-mechanical properties and more. The AFLOW data repository is publicly accessible online at aflow.org, with more than 1.7 million materials entries and a panoply of queryable computed properties. Tools to programmatically search and process the data, as well as to perform online machine learning predictions, are also available.