search

mzjb / DeepH-pack

Deep neural networks for density functional theory Hamiltonian.
GNU Lesser General Public License v3.0
240 stars 51 forks source link

It can cause band structure instability when using ABACUS' high LCAO base orbitals. #37

Closed sgyang345 closed 1 year ago

sgyang345 commented 1 year ago

1x1

mzjb commented 1 year ago

Thank you for reporting this issue. The non-physical band structure you calculated may be due to numerical problems caused by the overcompleteness of the basis set, which leads to the singularity of the overlap matrix. As described in the section A.7. Nearly singular basis sets: Strange results from small-unit-cell periodic calculation with many k-points of the FHI-aim documentation, one approach to address this is to project out the eigenvectors corresponding to small eigenvalues by computing the eigendecomposition of the overlap matrix, before solving the generalized eigenvalue problem.

We have added this feature to the DeepH diagonalization program in this commit to enhance numerical stability. You can use the command-line parameter ill_project to determine whether to project out linearly dependent basis functions and ill_threshold to set the threshold for eigenvalue magnitudes.

sgyang345 commented 1 year ago

Thank you very much for your reply, I have been able to get a normal band structure.

1x1_1