LCPQ/dft_tools
3
0
Fork 0
mirror of https://github.com/triqs/dft_tools synced 2024-12-23 04:43:42 +01:00
Code Issues Releases Wiki Activity

Merge pull request #199 from phibeck/w90

Browse Source 
Make sure Wannier90 Hamiltonian is Hermitian (remove diagonal complex entries)
This commit is contained in:
Alexander Hampel 2021-12-16 16:41:51 -05:00 committed by GitHub
commit cd1672cad7
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 13 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
python/triqs_dft_tools/converters/wannier90.py
View File

@ -486,6 +486,19 @@ class Wannier90Converter(ConverterTools):
else:
# make Fourier transform H(R) -> H(k) : it can be done one spin at a time
hamk = self.fourier_ham(hamr_full[isp])
# Sanity check if imaginary diagonal elements are zero, otherwise instabilties in lattice Gf!
diag_iterator = range(0,dim_corr_shells)
for ik in range(self.n_k):
if not numpy.allclose(hamk[ik][diag_iterator, diag_iterator].imag, 0, atol=1e-10):
mpi.report('ERROR: Wannier Hamiltonian has complex diagonal entries. '
+ f'First occurred at kpt {ik}:')
with numpy.printoptions(formatter={'float': '{:+.10f}'.format}):
mpi.report('\nWannier Hamiltonian diagonal, Fourier(H(r)), imaginary')
mpi.report(hamk[ik][diag_iterator, diag_iterator].imag)
mpi.MPI.COMM_WORLD.Abort(1)
# set imaginary part to zero
hamk[ik][diag_iterator, diag_iterator] = hamk[ik][diag_iterator, diag_iterator].real + 0*1j
# finally write hamk into hoppings
for ik in range(self.n_k):
hopping[ik, isp] = hamk[ik] - numpy.identity(numpy.max(n_orbitals)) * self.fermi_energy