No Arabic abstract
At ambient pressure tin transforms from its ground-state semi-metal $alpha$-Sn (diamond structure) phase to the compact metallic $beta$-Sn phase at 13$^circ$C (286K). There may be a further transition to the simple hexagonal $gamma$-Sn above 450K. These relatively low transition temperatures are due to the small energy differences between the structures, $approx 20$,meV/atom between $alpha$- and $beta$-Sn. This makes tin an exceptionally sensitive test of the accuracy of density functionals and computational methods. Here we use the high-throughput Automatic-FLOW (AFLOW) method to study the energetics of tin in multiple structures using a variety of density functionals. We look at the successes and deficiencies of each functional. As no functional is completely satisfactory, we look Hubbard U corrections and show that the Coulomb interaction can be chosen to predict the correct phase transition temperature. We also discuss the necessity of testing high-throughput calculations for convergence for systems with small energy differences.
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.
Within 4 different crystal structures, 2280 ternary intermetallic configurations have been investigated via high-throughput density functional theory calculations in order to discover new semiconducting materials. The screening is restricted to intermetallics with the equimolar composition TMX, where T is a transition metal from the Ti, V, Cr columns, Sr, Ba, Y and La, M an element from the first line of transition metals and X a sp elements (Al, P, Si, Sn and Sb), i.e. to a list of 24 possible elements. Since the calculations are done combinatorically, every possible ternary composition is considered, even those not reported in the literature. All these TMX configurations are investigated in the 4 most reported structure-types: TiNiSi, MgAgAs, BeZrSi and ZrNiAl. With an excellent agreement between calculations and literature for the reported stable phases, we identify 472 possible stable compounds among which 21 are predicted as non-metallic. Among these 21 compositions, 4 could be considered as new semiconductors.
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.
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).
Kinetic energy (KE) approximations are key elements in orbital-free density functional theory. To date, the use of non-local functionals, possibly employing system dependent parameters, has been considered mandatory in order to obtain satisfactory accuracy for different solid-state systems, whereas semilocal approximations are generally regarded as unfit to this aim. Here, we show that instead properly constructed semilocal approximations, the Pauli-Gaussian (PG) KE functionals, especially at the Laplacian-level of theory, can indeed achieve similar accuracy as non-local functionals and can be accurate for both metals and semiconductors, without the need of system-dependent parameters.