Page 1 of 1

Nonlinear response calculations with PBE functional

Posted: Thu Mar 01, 2018 12:31 pm
by ahoffman
Dear all,

I recently started using ABINIT to calculate Raman spectra of nanoporous materials, as this can be done efficiently via the DFPT formalism.
Performing linear response calculations at the PBE-D3(BJ) level of theory (ixc=11, vdw_xc=7), I could accurately reproduce phonon and IR spectra obtained earlier with the VASP code.
However, I did not manage to calculate the Raman spectrum as the PBE functional is not supported when doing nonlinear response simulations.
I found a post on this forum, initiated a couple of years ago, where it was stated that the implementation of this feature was almost completed.
Are there plans to make this type of simulations possible in the near future?

Thank you in advance.

Re: Nonlinear response calculations with PBE functional

Posted: Fri Mar 02, 2018 7:44 am
by ebousquet
Dear ahoffman,
The GGA is still not yet ready for non-linear calculations... What could be done in the meantime is to fix the cell parameters as you get them in GGA (or experimental ones), relax the internal coordinates with LDA and compute phonons and Raman stuff. You can verify if the phonon frequencies you get like that with LDA is strongly different than the ones you got with GGA, if not then this is totally fine.
Best wishes,
Eric

Re: Nonlinear response calculations with PBE functional

Posted: Thu Sep 06, 2018 1:44 pm
by mverstra
Another trick is to do everything in GGA (structure GS phonons) up to the 3rd derivatives, and force those in LDA with and explicit "ixc" for the final DFPT calculations. There will be a warning, but many people use this approximation as the GGA 3rd derivatives (not to mention above) are a real pain.