PAW carbon surface energy calculation diverges

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

Moderators: mverstra, joaocarloscabreu

Posts: 2
Joined: Fri Jun 11, 2021 12:23 pm

PAW carbon surface energy calculation diverges

Post by zanto » Mon Jun 14, 2021 8:39 am

Dear everyone,

I am trying to compute surface energy of carbon diamond structure (111) using PAW pseudopotential taken from (Type PAW (JTH v1.0), XC PBE). Everything is ok until I have 18 layers of C(111) or more with vacuum (does not matter how many angstroms of vacuum). Every computation is then oscillating and diverges. The same issue is also for WC(0001) in hexagonal structure. Everything works fine for both structures with NC SR (ONCVPSP v0.4.1), XC PBE pseudopotencial. Input and output files are attached.

I found this issue where was told that abinit is not ready for !! DFTP !! and GGA PAW pseudopotentials and there is a problem even with LDA when structure is 2D. Can this be the same problem?

Thank you

(2.62 KiB) Downloaded 347 times
(84.95 KiB) Downloaded 339 times