Method of explicit solving KS-equation can be chosen by variable "ks_solver" in INPUT file.
When "basis_type = pw", ks_solver
can be cg
, bpcg
or dav
. The bpcg
method only supports K-point parallelism currently. The default setting cg
is recommended, which is band-by-band conjugate gradient diagonalization method. There is a large probability that the use of setting of dav
, which is block Davidson diagonalization method, can be tried to improve performance.
When "basis_type = lcao", ks_solver
can be genelpa
or scalapack_gvx
. The default setting genelpa
is recommended, which is based on ELPA (EIGENVALUE SOLVERS FOR PETAFLOP APPLICATIONS) (https://elpa.mpcdf.mpg.de/) and the kernel is auto choosed by GENELPA(https://github.com/pplab/GenELPA), usually faster than the setting of "scalapack_gvx", which is based on ScaLAPACK(Scalable Linear Algebra PACKage)
We support stochastic DFT calculation (SDFT) or mixed stochastic-deterministic DFT (MDFT) with plane-wave basis [Phys. Rev. B 106, 125132 (2022)]. Different from traditional KSDFT with the explicit diagonalization method, SDFT and MDFT calculate physical quantities with trace of the corresponding operators. The advantages of SDFT and MDFT compared to the traditional KSDFT are the ability to simulate larger sizes and higher temperatures. In our package, SDFT and MDFT can be used by setting the esolver_type
parameter to sdft
for SCF calculations or MD calculations. To start with, you can refer to two examples and an explanation of the input variables.
When we have a hamiltonian, the electronic density can be calculated with:
where the Fermi-Dirac function smearing_method
, the parameter smearing_sigma
is equal the temperature nche_sto
represents the order of the expansion.
For physical quantities represented by operator
while MDFT calculates the trace as:
where
Here the number of KS orbitals nbands
while the number of stochastic orbitals nbands_sto
.
Besides, although SDFT does not diagonalize the hamiltonian, it can also caluclate DOS and electronic conductivities with parameters out_dos
and cal_cond
separately.