Logo Welcome to dfttk’s documentation!

_images/dfttk0.png

Project Goal

The goal of DFTTK is to make high-throughput first-principles calculations as simple as possible. The density functional theory (DFT) based software VASP is employed to perform first-principles calculations. In addition thermodynamic properties via the quasiharmonic approach, we proposed that any property, as long as it is dependent on the volume or stain, can be predicted using a quasi-static approach implemented by our group according to (i) the predicted property-volume/strain relationship from first-principles calculations directly and (ii) the volume/strain-temperature relationship of materials from the quasiharmonic approach.

For a given structure and elements, calculate all kinds of thermodynamic properties at finite temperature and pressure by first-principles approach based on density functional theory, including lattice vibration, thermal electron excitation, Seebeck coefficient, Lorenz number, effective charge carrier concentration etc.

Citing DFTTK

If you use DFTTK for work presented in a publication, we ask that you cite the following publication:

  1. Wang, M. Q. Liao, B. J. Bocklund, P. Gao, S. L. Shang, H. J. Kim, A. M. Beese, L. Q. Chen, and Z. K. Liu, DFTTK: Density functional theory toolkit for high-throughput lattice dynamics calculations, CALPHAD 75 (2021) 102355. doi:10.1016/j.calphad.2021.102355.

@article{Wang2021DFTTK,
         archivePrefix = {arXiv},
         arxivId = {2107.03966},
         author = {Yi Wang, Mingqing Liao, Brandon J. Bocklund, Peng Gao, Shun-Li Shang, Hojong Kim, Allison M. Beese, Long-Qing Chen, Zi-Kui Liu},
         doi = {10.1016/j.calphad.2021.102355},
         eprint = {2107.03966},
         issn = {0364-5916},
         journal = {CALPHAD},
         month = {december},
         pages = {102355},
         title = {{DFTTK: Density functional theory toolkit for high-throughput lattice dynamics calculations}},
         year = {2019}
}

High-throughput calculations

By its definition, the term of “first-principles” represents a philosophy that the prediction is to be based on a basic, fundamental proposition or assumption that cannot be deduced from any other proposition or assumption. This implies that the computational formulations are based on the most fundamental theory of quantum mechanics - Schrödinger equation or density functional theory (DFT) and the inputs to the calculations must be based on well-defined physical constants – the nuclear and electronic charges. In another word, once the atomic species of an assigned material are known, the theory should predict the energy of all possible crystalline structures, without invoking any phenomenological fitting parameters.

However, to perform DFT calculations in reality, it still needs the user to have extensive experiences on a variety of parameter choices and a lot of human handling on numerical or system exceptions. In the last decade, we have been working on solving the problem by integrating our experiences accumulated on high-throughout DFT calculations into a software package named as DFTTK (DFT based toolkits) and opened to the community (https://www.dfttk.org).

The main functions of DFTTK

  • Structure maker by protype and elemental substitution;

  • Robust 0 K equilibrium volume optimization;

  • Robust 0 K energy-volume curve optimization;

  • Quasiharmonic phonon calculation;

  • Born effective charge calculation;

  • Elastic constant calculations.

  • MongoDB database management

  • Thermodynamic calculations and figure plots

_images/dfttk1.png

DFTTK features

To perform DFT calculation using DFTTK, the user only needs to name the structure file called POSCAR by VASP, either prepared by user or produced by DFTTK by elemental substation on given prototype. DFTTK is developed on atomate from the Materials Project which is built on three open-source Python libraries. The main benefits of atomate are its flexibility and data management platform, in particular the numerical convergence control and computational exception handling. DFTTK is able to predict properties at finite temperatures by phonon or Debye model for both stoichiometric and solution phases, featured by:

  • High-throughput DFT calculation and postprocess;

  • Postprocess plenty of data stored in MongoDB with one simple command;

  • Compatible with Yphon package and phonopy;

  • Can recover data from certain fizzled calculations;

  • Can account thermal electron contribution to thermodynamic properties;

  • Can calculate thermodynamic properties at 0 K and a few tenth K;

  • Can perform doping calculations for semiconductors or thermoelectric materials under rigid band approximation;

  • Can account the effect of thermal expansion/temperature on Seebeck coefficient, Lorenz number, thermal carrier concentrations;

  • Automatic plot figures for more than 20 thermodynamic properties in the publishable resolution, including atomic volume, free energy, entropy, enthalpy, linear thermal expansion coefficient, isobartic specific heat, constant volume specific heat, lattice only specific heat, bulk modulus, Debye temperature, Seebeck coefficient, Lorenz number, absolute thermal electric force, etc.

The outline of flow chart of DFTTK

Provide robust density functional theory workflows for calculating thermodynamic properties in temperature and composition space. The work flow control of DFTTK is given by the following blue print.

_images/workflow.png

Robust relaxation scheme of DFTTK

_images/Robust.png

DFTTK makes it automatic to going through these ISIF options by the workflow shown in the figure in the right

The purpose of DFTTK robust relaxation is to find the lowerst energy structure within a given threshold of a given sysmmetry. This information is espetically useful in the case, for example

  1. Random sulution phase;

  2. High temperature phase which is unstable in low temperature;

  3. The calculation of energetics for endnumber for CALPHAD modeling.

The solution is to make the best use of ISIF control parameter in VASP, which controls:

ISIF

calc. force

relax stress tensor

relax ions

change cell shape

change cell volume

0

yes

no

yes

no

no

1

yes

trace only

yes

no

no

2

yes

yes

yes

no

no

3

yes

yes

yes

yes

yes

4

yes

yes

yes

yes

no

5

yes

yes

no

yes

no

6

yes

yes

no

yes

yes

7

yes

yes

no

no

yes

Contents:

Indices and tables

References

Wang2021

Wang et al., CALPHAD 75 (2021) 102355. doi:10.1016/j.calphad.2021.102355

Bocklund2019

Bocklund et al., MRS Communications 9(2) (2019) 1–10. doi:10.1557/mrc.2019.59

Dinsdale1991

Dinsdale, Calphad 15(4) (1991) 317-425, doi:10.1016/0364-5916(91)90030-N

Otis2016

Otis, Ph.D. Dissertation, The Pennsylvania State University (2016). https://etda.libraries.psu.edu/catalog/s1784k73d

Otis2017

Otis et al., JOM 69 (2017) doi:10.1007/s11837-017-2318-6

Shang2010

Shang, Wang, and Liu, Magnes. Technol. 2010 617-622 (2010).