Automation methodologies and large-scale validation for GWGW, towards high-throughput GWGW calculations

The search for new materials, based on computational screening, relies on methods that accurately predict, in an automatic manner, total energy, atomic-scale geometries, and other fundamental characteristics of materials. Many technologically important material properties directly stem from the electronic structure of a material, but the usual workhorse for total energies, namely density-functional theory, is plagued by fundamental shortcomings and errors from approximate exchange-correlation functionals in its prediction of the electronic structure. At variance, the $GW$ method is currently the state-of-the-art {\em ab initio} approach for accurate electronic structure. It is mostly used to perturbatively correct density-functional theory results, but is however computationally demanding and also requires expert knowledge to give accurate results. Accordingly, it is not presently used in high-throughput screening: fully automatized algorithms for setting up the calculations and determining convergence are lacking. In this work we develop such a method and, as a first application, use it to validate the accuracy of $G_0W_0$ using the PBE starting point, and the Godby-Needs plasmon pole model ($G_0W_0^\textrm{GN}$@PBE), on a set of about 80 solids. The results of the automatic convergence study utilized provides valuable insights. Indeed, we find correlations between computational parameters that can be used to further improve the automatization of $GW$ calculations. Moreover, we find that $G_0W_0^\textrm{GN}$@PBE shows a correlation between the PBE and the $G_0W_0^\textrm{GN}$@PBE gaps that is much stronger than that between $GW$ and experimental gaps. However, the $G_0W_0^\textrm{GN}$@PBE gaps still describe the experimental gaps more accurately than a linear model based on the PBE gaps.Comment: 12 pages, 11 figure

Accuracy of generalized gradient approximation functionals for density functional perturbation theory calculations

We assess the validity of various exchange-correlation functionals for computing the structural, vibrational, dielectric, and thermodynamical properties of materials in the framework of density-functional perturbation theory (DFPT). We consider five generalized-gradient approximation (GGA) functionals (PBE, PBEsol, WC, AM05, and HTBS) as well as the local density approximation (LDA) functional. We investigate a wide variety of materials including a semiconductor (silicon), a metal (copper), and various insulators (SiO$_2$ $\alpha$-quartz and stishovite, ZrSiO$_4$ zircon, and MgO periclase). For the structural properties, we find that PBEsol and WC are the closest to the experiments and AM05 performs only slightly worse. All three functionals actually improve over LDA and PBE in contrast with HTBS, which is shown to fail dramatically for $\alpha$-quartz. For the vibrational and thermodynamical properties, LDA performs surprisingly very good. In the majority of the test cases, it outperforms PBE significantly and also the WC, PBEsol and AM05 functionals though by a smaller margin (and to the detriment of structural parameters). On the other hand, HTBS performs also poorly for vibrational quantities. For the dielectric properties, none of the functionals can be put forward. They all (i) fail to reproduce the electronic dielectric constant due to the well-known band gap problem and (ii) tend to overestimate the oscillator strengths (and hence the static dielectric constant)

First-principles calculations and experimental studies of: XYZ 2 thermoelectric compounds: Detailed analysis of van der Waals interactions

First-principles calculations can accelerate the search for novel high-performance thermoelectric materials. However, the prediction of the thermoelectric properties is strongly dependent on the approximations used for the calculations. Here, thermoelectric properties were calculated with different computational approximations (i.e., PBE-GGA, HSE06, spin-orbit coupling and DFT-D3) for three layered XYZ2 compounds (TmAgTe2, YAgTe2, and YCuTe2). In addition to the computations, the structural, electrical and thermal properties of these compounds were measured experimentally and compared to the computations. An enhanced prediction of the crystal structure and heat capacity was achieved with the inclusion of van der Waals interactions due to more accurate modeling of the interatomic forces. In particular, a large shift of the acoustic phonons and low-frequency optical phonons to lower frequencies was observed from the dispersion-optimized structure. From the phonon dispersion curves of these compounds, the ultralow thermal conductivity in the investigated XYZ2 compounds could be described by a recent developed minimum thermal conductivity model. For the prediction of the electrical conductivity, a temperature-dependent relaxation time was used, and it was limited by acoustic phonons. While HSE06 has only a small influence on the electrical properties due to a computed band gap energy of >0.25 eV, the inclusion of both van der Waals interactions and spin-orbit coupling leads to a more accurate band structure, resulting in better prediction of electrical properties. Furthermore, the experimental thermoelectric properties of YAgTe2, TmAg0.95Zn0.05Te2 and TmAg0.95Mg0.05Te2 were measured, showing an increase in zT of TmAg0.95Zn0.05Te2 by more than 35% (zT = 0.47 ± 0.12) compared to TmAgTe2

High-throughput screening of perovskite alloys for piezoelectric performance and thermodynamic stability

We screen a large chemical space of perovskite alloys for systems with optimal properties to accommodate a morphotropic phase boundary (MPB) in their composition-temperature phase diagram, a crucial feature for high piezoelectric performance. We start from alloy end points previously identified in a high-throughput computational search. An interpolation scheme is used to estimate the relative energies between different perovskite distortions for alloy compositions with a minimum of computational effort. Suggested alloys are further screened for thermodynamic stability. The screening identifies alloy systems already known to host an MPB and suggests a few others that may be promising candidates for future experiments. Our method of investigation may be extended to other perovskite systems, e.g., (oxy-)nitrides, and provides a useful methodology for any application of high-throughput screening of isovalent alloy systems

The thermodynamic scale of inorganic crystalline metastability

The space of metastable materials offers promising new design opportunities for next-generation technological materials, such as complex oxides, semiconductors, pharmaceuticals, steels, and beyond. Although metastable phases are ubiquitous in both nature and technology, only a heuristic understanding of their underlying thermodynamics exists. We report a large-scale data-mining study of the Materials Project, a high-throughput database of density functional theory–calculated energetics of Inorganic Crystal Structure Database structures, to explicitly quantify the thermodynamic scale of metastability for 29,902 observed inorganic crystalline phases. We reveal the influence of chemistry and composition on the accessible thermodynamic range of crystalline metastability for polymorphic and phase-separating compounds, yielding new physical insights that can guide the design of novel metastable materials. We further assert that not all low-energy metastable compounds can necessarily be synthesized, and propose a principle of ‘remnant metastability’—that observable metastable crystalline phases are generally remnants of thermodynamic conditions where they were once the lowest free-energy phase.United States. Dept. of Energy. Office of Basic Energy Sciences (DE-AC02-05CH11231)United States. Dept. of Energy. Office of Basic Energy Sciences (contract UGA-0-41029-16/ER392000

High-throughput screening of perovskite alloys for piezoelectric performance and thermodynamic stability

We screen a large chemical space of perovskite alloys for systems with optimal properties to accommodate a morphotropic phase boundary (MPB) in their composition-temperature phase diagram, a crucial feature for high piezoelectric performance. We start from alloy end points previously identified in a high-throughput computational search. An interpolation scheme is used to estimate the relative energies between different perovskite distortions for alloy compositions with a minimum of computational effort. Suggested alloys are further screened for thermodynamic stability. The screening identifies alloy systems already known to host an MPB and suggests a few others that may be promising candidates for future experiments. Our method of investigation may be extended to other perovskite systems, e.g., (oxy-)nitrides, and provides a useful methodology for any application of high-throughput screening of isovalent alloy systems

Good Vibrations : The evolution of whisking in small mammals

Special Issue: Extreme Anatomy: Living Beyond the Edge. January 2020Abstract While most mammals have whiskers, some tactile specialists - mainly small, nocturnal and arboreal species - can actively move their whiskers in a symmetrical, cyclic movement called whisking. Whisking enables mammals to rapidly, tactually scan their environment in order to efficiently guide locomotion and foraging in complex habitats. The muscle architecture that enables whisking is preserved from marsupials to primates, prompting researchers to suggest that a common ancestor might have had moveable whiskers. Studying the evolution of whisker touch sensing is difficult, and we suggest that measuring an aspect of skull morphology that correlates with whisking would enable comparisons between extinct and extant mammals. We find that whisking mammals have larger infraorbital foramen (IOF) areas, which indicates larger infraorbital nerves and an increase in sensory acuity. While this relationship is quite variable and IOF area cannot be used to solely predict the presence of whisking, whisking mammals all have large IOF areas. Generally, this pattern holds true regardless of an animal's substrate preferences or activity patterns. Data from fossil mammals and ancestral character state reconstruction and tracing techniques for extant mammals suggest that whisking is not the ancestral state for therian mammals. Instead, whisking appears to have evolved independently as many as seven times across the clades Marsupialia, Afrosoricida, Eulipotyphla and Rodentia, with Xenarthra the only placental superordinal clade lacking whisking species. However, the term whisking only captures symmetrical and rhythmic movements of the whiskers, rather than all possible whisker movements, and early mammals may still have had moveable whiskers. This article is protected by copyright. All rights reserved.Peer reviewe