Germanium and structure optimization

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
herman
Posts: 6
Joined: Sun Jan 09, 2011 10:55 pm

Germanium and structure optimization

Post by herman » Mon Jan 10, 2011 10:58 am

Hello everybody,

I'd have a question regarding structure optimization performed by the following parameters:

optcell 1
ionmov 3
ntime 10
dilatmx 1.05
ecutsm 0.5

Germanium structure optimization done using LDA/GGA gives the results 5.46 / 5.6 [Angstrom] which doesn't seem right according to known values. Does anybody have any experience with obtaining correct lattice constants for germanium (or basically any other element)? Is that mismatch the normal Abinit behaviour or there's something I'm missing?

I would kindy appreciate any suggestions.
Best regards,
Michal

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: Germanium and structure optimization

Post by jzwanzig » Mon Jan 10, 2011 3:47 pm

First, I would always use ionmov 2 even for simple systems. I've never had a problem with this choice.

Second, my understanding is that germanium is an (in)famous case where Kohn-Sham DFT fails, especially for LDA--that is, it predicts Ge to be a metal, rather than a semiconductor. What kind of pseudopotential are you using? Before you try to optimize the structure you should compute the band structure and see if you are getting it (at least qualitatively) correct. You may need to go beyond simple KS DFT and use exact exchange or something similar.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

herman
Posts: 6
Joined: Sun Jan 09, 2011 10:55 pm

Re: Germanium and structure optimization

Post by herman » Wed Jan 19, 2011 10:47 pm

jzwanzig wrote:First, I would always use ionmov 2 even for simple systems. I've never had a problem with this choice.

Second, my understanding is that germanium is an (in)famous case where Kohn-Sham DFT fails, especially for LDA--that is, it predicts Ge to be a metal, rather than a semiconductor. What kind of pseudopotential are you using? Before you try to optimize the structure you should compute the band structure and see if you are getting it (at least qualitatively) correct. You may need to go beyond simple KS DFT and use exact exchange or something similar.

Thank you very much for your suggestions (sorry for the delay in replying). I've been using Hartwigsen-Goedecker-Hutter psp for Ge. I'll try it your way.

I was also going to see the result of running the problem with the new Tran-Blaha exchange-correlation functional (TB/MBJ). I was wondering if I can find an appropriate psp file for that job, but no luck so far.

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: Germanium and structure optimization

Post by jzwanzig » Thu Jan 20, 2011 7:10 pm

The Tran-Blaha metaGGA is supported by libxc, and libxc can be linked into abinit as a plug-in. I'm pretty sure that PAW in abinit does not support metaGGAs, though--also, I know of no PAW generation package (at least not one compatible with abinit) that supports metaGGA's. The package that most of us in the abinit community use for PAW generation is atompaw, and this package does support the use of libxc but does not support any metaGGA's, only LDA and GGA's. I'm not sure about norm-conserving pseudopotentials, whether any available package compatible with abinit supports metaGGA's.

Having said all that I did play around with Ge using PAW and PBE, and found that the band gap vanishes (wrong answer, but expected), but that it can be opened up by using some portion of exact exchange (accessed in abinit through use of useexech, lexexch, and exchmix variables, and only in the PAW case). For my Ge PAW set the optimized cell constant is contracted by 1.7% compared to experiment (normally PBE leads to a small expansion, but note here I am not using PBE but a mix of PBE and exact exchange).

Joe
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

hicpalm
Posts: 44
Joined: Tue Feb 09, 2010 4:33 pm

Re: Germanium and structure optimization

Post by hicpalm » Fri Mar 11, 2011 10:03 am

Hi,
my Ge band gap with the GGA-PBE vanishs too. This is confirmed in the band plot.
however in the log file one can read:

Code: Select all

 === Gap info ===
  >>>> For spin  1
   Minimum optical gap =   0.6811 [eV], located at k-point      :  -0.0417  0.0000 -0.0417
   Fundamental gap     =   0.1904 [eV], Top of valence bands at :  -0.0417  0.0000 -0.0417
                                        Bottom of conduction at :   0.0000  0.0000 -0.4167


which I don't understand clearly.
thank you for comments.

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: Germanium and structure optimization

Post by jzwanzig » Wed Mar 16, 2011 11:32 am

Can you post your input file?
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

hicpalm
Posts: 44
Joined: Tue Feb 09, 2010 4:33 pm

Re: Germanium and structure optimization

Post by hicpalm » Wed Mar 16, 2011 2:44 pm

thank you for replying professor,
here is my input file (for scf calculation)

Code: Select all

 ixc 11 

 acell    7.6996626338E+00  7.6996626338E+00  7.6996626338E+00
 rprim    0.0000000000E+00  7.0710678119E-01  7.0710678119E-01
          7.0710678119E-01  0.0000000000E+00  7.0710678119E-01
          7.0710678119E-01  7.0710678119E-01  0.0000000000E+00

 ntypat 1
 znucl 32
 natom 2
 typat 1 1
 xred   -1.3877787808E-17 -1.3877787808E-17 -1.3877787808E-17
         2.5000000000E-01  2.5000000000E-01  2.5000000000E-01

 ecut 520 eV
 pawecutdg 50
 occopt 4
 tsmear 0.1 eV
 nband 19
 
 kptopt 1
 kptrlatt  -18  6  6   6  6  6   6  6 -18

 toldfe 2.0d-10



Actually, I think that it was just a problem of grid shift because when I give a gamma-centered grid with

Code: Select all

 shiftk 0 0 0


the band gap disappears. This means that in the shifted grid I missed the minimum of the valence and conduction bands
which explains that the code gives a band gap as the maximum of the valence band is no longer equal to the minimum of the conduction band.

is this correct ?
thank you professor.

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: Germanium and structure optimization

Post by jzwanzig » Wed Mar 16, 2011 2:48 pm

I think you are right, although for band gaps I never take them from the base output file but rather run a full band structure calculation. By the way your pawecutdg looks very high to me, is this set from a convergence study? This parameter needs to be checked for convergence but I never need to set it that high, and it chews up a LOT of memory.
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

hicpalm
Posts: 44
Joined: Tue Feb 09, 2010 4:33 pm

Re: Germanium and structure optimization

Post by hicpalm » Wed Mar 16, 2011 5:46 pm

although for band gaps I never take them from the base output file but rather run a full band structure calculation


the previous input was just the first step to obtain the self consistent density used as starting point for band structure calculations (though I missed the prtden 1)

By the way your pawecutdg looks very high to me, is this set from a convergence study? This parameter needs to be checked for convergence but I never need to set it that high, and it chews up a LOT of memory.


I need differences converged to 1 meV and I just want to avoid any new errors added up to those from ecut and k-point samplings ( I have even used values greater than this value). In fact convergence studies give a slightly smaller values.
For the memory it's really memory consuming and it was even the source of some buggy runs.

thank you professor for your instructive comments for this post an all other posts in this forum.

Locked