paw for covalent molecules?

External optional components of Abinit (BigDFT, Wannier90...) and codes using/providing data from/to Abinit (AtomPAW, ONCVPSP, EXC, DP, Yambo...)

Moderators: ebousquet, bxu

Locked
rundong
Posts: 16
Joined: Thu May 20, 2010 1:52 am

paw for covalent molecules?

Post by rundong » Mon Aug 16, 2010 2:58 am

Hi, there,

The other day I was optimizing O2 molecule with paw and the log file returned error message, since the 2 O atoms overlap was bigger than allowed. To my understanding, this issue should alway pops up, because for covalent molecules, the bond lengths will be always so small, resulting in big overlap. Therefore, paw will be NOT usable?!
Any insights will be greatly appreciate in advance.

rundong

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

Re: paw for covalent molecules?

Post by jzwanzig » Mon Aug 16, 2010 3:10 pm

No, PAW is still usable, but you may have to generate a PAW data set for oxygen with a small enough cut-off radius. We have a nice tutorial posted for how to make PAW data sets, you could start there.
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

rundong
Posts: 16
Joined: Thu May 20, 2010 1:52 am

Re: paw for covalent molecules?

Post by rundong » Mon Aug 16, 2010 4:46 pm

Prof Zwanziger,

Thank you very much for the insight. You saved my life! Although I did notice that the oxygen paw was usable, I got a lot of warnings. And yes I downloaded the the code on how to generate paw for ABINIT; however, I am not confident that I can produce as good paws as you did. I have used your copper paw with satisfactory result: the lattice constant was within 1% error compared with experimental one.
This post solved my problem.

Regards,
rundong

Locked