ABINIT Discussion Forums

Hi, the following may be a bug, I've uncovered it as I work on the ffield tutorial.

In a DFPT run, in order obtain the clamped and relaxed piezoelectric constants, it works if the phonons, strains, and electric fields are computed in a single dataset. However, it fails if the electric field perturbation is computed separately and then included through mrgddb. This seems like it can't be the expected response.

Any suggestions?

thanks,
Joe