Nonlinear response calculations with PBE functional

Phonons, DFPT, electron-phonon, electric-field response, mechanical response…

Moderators: mverstra, joaocarloscabreu

Locked
ahoffman
Posts: 1
Joined: Thu Mar 01, 2018 11:41 am

Nonlinear response calculations with PBE functional

Post by ahoffman » Thu Mar 01, 2018 12:31 pm

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.

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Nonlinear response calculations with PBE functional

Post by ebousquet » Fri Mar 02, 2018 7:44 am

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

mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

Re: Nonlinear response calculations with PBE functional

Post by mverstra » Thu Sep 06, 2018 1:44 pm

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.
Matthieu Verstraete
University of Liege, Belgium

Locked