ABINIT Discussion Forums

Hi everybody,
I am calculate the polarization, elastic constants and piezoelectric properties. I have obtained the correct results and are published several years ago throuth the anaddb in the abinit-7.6.3.
But I do not obtain the results when I use the anaddb in the abinit-8.10.1.
The input file for the anaddb code is as follows:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
elaflag 3 !the flag for the elastic constant
piezoflag 3 !the flag for the piezoelectric constant
instrflag 1 !the flag for the internal strain tensor
dieflag 3
polflag 1

!the effective charge part
chneut 1 !enforce Born effective charge neutrality
# This line added when defaults were changed (v5.3) to keep the previous, old behaviour
asr 0
symdynmat 0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

The error is as follow:

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
--- !ERROR
src_file: anaddb.F90
src_line: 383
mpi_rank: 0
message: |
DDB file does not contain GS etotal
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

I will be very grateful if you give me some advice.
Thank you.
Q. Li

Dear Q. Li,
The DDB file should contain the DDB of the ground state calculation on the top of the DFPT ones to have the residual of the forces you have in your system. If you don't have it, run a single point energy GS to get the DDB and merge it with the DFPT ones (However, this problem should have been observed with older version of anaddb...).
Best wishes,
Eric