Page 1 of 1

Hessian matrix IBRION=6

Posted: Wed Aug 23, 2017 12:22 pm
by Tobias Kraemer
Hello,

it appears that the Hessian matrix is not printed in the OUTCAR file when IBRION=6 (i.e. symmetry) is used. The line SECOND DERIVATIVES (NOT SYMMETRIZED) is missing in the OUTCAR. I'd assume that VASP would fill the matrix based on the symmetry of the system. This is what I get from the VASP wiki: For IBRION=6, only symmetry inequivalent displacements are considered, and the remainder of the Hessian matrix is filled using symmetry consideration

Could someone please give some advice how to print the Hessian in the OUTCAR in this case?

Thanks

Tobias

Re: Hessian matrix IBRION=6

Posted: Mon Aug 28, 2017 1:30 pm
by alex
Hi Tobias,

it looks like it got lost over time. Go to the file finite_diff.F, search for 'clumsy and unreadable output', make the 7 lines active and recompile.
It might solve your issue.

Cheers,

alex

Re: Hessian matrix IBRION=6

Posted: Mon Sep 04, 2017 3:58 pm
by Tobias Kraemer
Alex,

thanks for your suggestion. This is indeed rather strange, but I will have a look into this and perhaps this will remedy the problem.
I don't understand why it normally works just fine for IBRION=5.

Tobias

Re: Hessian matrix IBRION=6

Posted: Sat Sep 09, 2017 8:46 am
by alex
Hi Tobias,

in the early days of VASP 4.6 the hessian was printed twice, one as is and then symmetrized. However at this time, symmetry was not used during the calculation. Maybe some of the developers had the opinion - hm, we don't need the symmetrized version of a symmetric hessian anyway and just removed the printing, leaving us with no hessian at all. :-)
We are all humans...

Cheers,

alex