No k-resolved spectral function as in DMFT tutorial

Running Abinit Many-body theory
1

Dear all,
I use the version 9.4.0 of abinit. I have an issue in generating the k-resolved spectral function as in section 7 of the tutorial

https://docs.abinit.org/tutorial/dmft/#7-electronic-structure-of-srvo3-k-resolved-spectral-function

Despite I follow all the steps as described in the page, I don’t get any file

tdmft_5o_DS3_DFTDMFT_SpectralFunction_kresolved_from_realaxisself

and the calculation ends normally. Please, attached you find the log_5_DS3 and the tdmft_5.abo files for more information - I renamed them as *.in because the page was forbidding any attachment with extensions like txt, dat, log…
Thanks a lot in advance for your help!

Antonio
tdmft_5.abo.in (518 KB)
log.in (426 KB)

2

Dear Antonio,

Thanks for your message.
The error comes from a bug that has been added in the source file.
You have to comment the line “dtset%usedmft=0” in m_invars2.F90 such that you have (around line 3605):

if (dtset%nbandkss /= 0 .and. dtset%usedmft/=0) then
dtset%plowan_compute=1
!dtset%usedmft=0
else if (dtset%optdriver==3) then
dtset%plowan_compute=10
else if(dtset%optdriver==4) then

With this correction, the file tdmft_5o_DS3_DFTDMFT_SpectralFunction_kresolved_from_realaxisself is written and
the k resolved band structure can be plotted.

Next ABINIT version will include the bug correction.
Best regards
Bernard

3

Dear Bernard,
now it works. Thanks a lot for the help!

All the best

Antonio

4

Hello everyone,

I followed the instruction and commented the line dtset%usedmft=0 in m_invars2.F90 file. What is the next step if my abinit was installed through Homebrew? I cannot compile the abinit-9.4.2 file since there is no build directory for my case.
Thank you!

Best regards,
WHPeng

5

Dear Bernard,

I think the bug is not fixed yet in the version 9.6.2.

It does not show any SpectralFunction_kresolved_from_realaxisself file.

Also, I tried with “dmft_kspectral_func” (written in the doc), and it says no such variable.

Best regards,
R

6

FYI, it worked in version 9.8.3.