ترغب بنشر مسار تعليمي؟ اضغط هنا

DFTTK: Density Functional Theory Tool Kit for High-throughput Calculations of Thermodynamic Properties at Finite Temperatures

115   0   0.0 ( 0 )
 نشر من قبل Yi Wang
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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 explore the derivation of interatomic exchange interactions in ferromagnets within density-functional theory (DFT) and the mapping of DFT results onto a spin Hamiltonian. We delve into the problem of systems comprising atoms with strong spontaneou s moments together with atoms with weak induced moments. All moments are considered as degrees of freedom, with the strong moments thermally fluctuating only in angle and the weak moments thermally fluctuating in angle and magnitude. We argue that a quadratic dependence of the energy on the weak local moments magnitude, which is a good approximation in many cases, allows for an elimination of the weak-moment degrees of freedom from the thermodynamic expressions in favor of a renormalization of the Heisenberg interactions among the strong moments. We show that the renormalization is valid at all temperatures accounting for the thermal fluctuations and resulting in temperature-independent renormalized interactions. These are shown to be the ones directly derived from total-energy DFT calculations by constraining the strong-moment directions, as is done e.g. in spin-spiral methods. We furthermore prove that within this framework the thermodynamics of the weak-moment subsystem, and in particular all correlation functions, can be derived as polynomials of the correlation functions of the strong-moment subsystem with coefficients that depend on the spin susceptibility and that can be calculated within DFT. These conclusions are rigorous under certain physical assumptions on the measure in the magnetic phase space. We implement the scheme in the full-potential linearized augmented plane wave method using the concept of spin-spiral states, considering applicable symmetry relations and the use of the magnetic force theorem. Our analytical results are corroborated by numerical calculations employing DFT and a Monte Carlo method.
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 trad eoffs. 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.
We revise critically existing approaches to evaluation of thermodynamic potentials within the Greens function calculations at finite electronic temperatures. We focus on the entropy and show that usual technical problems related to the multivalued na ture of the complex logarithm can be overcome. This results in a simple expression for the electronic entropy, which does not require any contour integration in the complex energy plane. Properties of the developed formalism are discussed and its illustrating applications to selected model systems and to bcc iron with disordered local magnetic moments are presented as well.
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).
Density functional theory (DFT) calculations are used to investigate the electronic and magnetic structures of a two-dimensional (2D) monolayer Li$_{2}$N. It is shown that bulk Li$_{3}$N is a non-magnetic semiconductor. The non-spinpolarized DFT calc ulations show that $p$ electrons of N in 2D Li$_{2}$N form a narrow band at the Fermi energy $E_{rm{F}}$ due to a low coordination number, and the density of states at the Fermi energy ($g(E_{rm{F}}$)) is increased as compared with bulk Li$_{3}$N. The large $g(E_{rm{F}}$) shows instability towards magnetism in Stoners mean field model. The spin-polarized calculations reveal that 2D Li$_{2}$N is magnetic without intrinsic or impurity defects. The magnetic moment of 1.0,$mu_{rm{B}}$ in 2D Li$_{2}$N is mainly contributed by the $p_{z}$ electrons of N, and the band structure shows half-metallic behavior. {Dynamic instability in planar Li$_{2}$N monolayer is observed, but a buckled Li$_{2}$N monolayer is found to be dynamically stable.} The ferromagnetic (FM) and antiferromagnetic (AFM) coupling between the N atoms is also investigated to access the exchange field strength. {We found that planar (buckled) 2D Li$_{2}$N is a ferromagnetic material with Curie temperature $T_{c}$ of 161 (572) K.}
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا