Hello,
I ran into some issues when trying to perform DFT+DMFT calculations with the VASP version 6.5.1 and TRIQS following the tutorial given on the VASP Wiki site:
https://www.vasp.at/wiki/index.php/DFT% ... lculations.
Step 1 ("Perform a SCF DFT calculation") worked as expected.
In step 2 ("Convert the VASP output to TRIQS input"), the python command plo_converter.generate_and_output_as_text('plo.cfg', vasp_dir='./') produced at first a warning:
[WARNING]: Error reading Efermi from LOCPROJ, trying DOSCAR..."
and then, after some output, stopped with an error:
File "/p/project1/espo/klemenz1/apps/lib/python3.12/site-packages/triqs_dft_tools/converters/plovasp/elstruct.py", line 89, in __init__
assert nk_plo == self.nktot, "PLOCAR is inconsistent with IBZKPT (number of k-points)"
^^^^^^^^^^^^^^^^^^^^
AssertionError: PLOCAR is inconsistent with IBZKPT (number of k-points)
According to the TRIQS DFTTools documentation (https://triqs.github.io/dft_tools/unsta ... _vasp.html, at the bottom of the page, Remarks for VASP older than 6.5.0), both of these issues were known for VASP versions older than 6.5.0:
The interface works correctly only if the k-point symmetries are turned off during the VASP run (ISYM=-1).
The interface currently supports only collinear-magnetism calculation (this implies no spin-orbit coupling) and spin-polarized projectors have not been tested.
The converter needs the correct Fermi energy from VASP, which is read from the LOCPROJ file. However, VASP by default does not output this information.
As pointed out in the first bullet point, turning off the symmetry with ISYM=-1 solves the k-point inconsistency. This indeed worked in my case, however, as said in the beginning, I used 6.5.0 (and even tried with 6.5.1), where this problem should not show up anymore. With the default, ISYM=2, as in the tutorial example (where ISYM is ot explicitly specified), the problem is still there.
Also, as mentioned above, the warning corresponding to the Fermi level still showed up for VASP 6.5.0 and 6.5.1. In the locproj.F file, the line says
WRITE(99,'(4I6," # of spin, # of k-points, # of bands, # of proj" )') NS,NK,NB,NF
, i.e., the Fermi energy is not printed. The manual patches provided on https://triqs.github.io/dft_tools/unsta ... _vasp.html can only partly been done in the specified way, since they correspond to older VASP versions than 6.5.0, and in the newer versions, some of the function definitions are different.
Can you help me with my confusion?
Best,
Daniel
