[Pw_forum] Generation of H and O GIPAW PP

jérôme cuny jerome.cuny at hotmail.fr
Tue Aug 28 22:14:42 CEST 2012


Dear Davide,

I have just tried the PPs you sent me and they works perfectly. I also understand what I was doing wrong in my initial input.
So thanks a lot for the help.

However, in my initial e-mail I forget to send the output corresponding to the bad H PP, where a core contribution appears. I send it to you now (as well as the corresponding PP generation input). Maybe that will help a future pwscf user for not doing the same wrong thing as me.

Best Regards,

Jérôme Cuny



Date: Tue, 28 Aug 2012 15:12:35 +0200
From: davide.ceresoli at istm.cnr.it
To: pw_forum at pwscf.org
Subject: Re: [Pw_forum] Generation of H and O GIPAW PP

Dear Jérôme,
     the two output files you attached do not show large differences,
The core contribution is only for oxygen, in both cases.
 
What can lead to wrong/inaccurate results is the linear dependence
between GIPAW projectors (in your "not-good" calculation). Moreover,
I strongly advise to generate always two projectors per angular
momentum in the &test section.
 
I've never tried to generate hydrogen with 'p' locality. Maybe someone
else can comment. In any case, I'm attaching original the input files
for the hydrogen and oxygen pseudopotentials. You can try to add to them.
 
Best wishes,
     Davide
 
 
On 08/28/2012 02:37 PM, jérôme cuny wrote:
> Dear quantum-espresso users,
>
> I would like to share with you a problem I have for which I hope you could help me.
>
> I am trying to generate H and O norm-conserving PP to calculate the NMR
> properties of water molecules using the GIPAW module.
> I am using the 4.3.2 version of both quantum-espresso and gipaw (I join to this
> e-mail the make.sys file I used for the compilation).
>
> My problem is that, even if the generation of the PPs seems fine, their use by
> the gipaw module leads to bad results.
> I consider as good results the results I obtain using the H and O PP provided in
> the example folder of the gipaw module (see file water-nmr_1.log).
>
> Using the input file H.pbe-tm_good.in for the PP generation of H, the result I
> obtain for one water molecule in a box is not too bad (see water-nmr_1_good.log)
>
> However, my problem appear when I want to use as local part of the PP, the p
> orbital (see H.pbe-tm_bad.in and water-nmr_1_bad.log). Strangely, when I do
> this, I always obtain a contribution to the total chemical shit coming from some
> core orbitals! The same problem occur when I generate the O PP (in that case I
> use as input, the data I found in a presentation on GIPAW, see fig
> O_PP_Slide.tiff) as the core contribution become really different from the one
> of the reference calculation.
>
> So, does somebody see what I do wrong in my PP generation? Or, of course, is
> there any other source of error that could lead to this problem?
>
> I have checked the previous postings on the subject (GIPAW+PP) and I have not
> find one that answered my question. If I missed the good one, I am sorry.
>
> Best Regards,
>
> Jérôme Cuny
> ETH-Zurich/USI-Campus
>
> --------------------------------------------------------------------
> Jérôme Cuny
> Department of Chemistry and Applied Biochemistry
> ETH-Zurich
> USI-Campus, Via Giuseppe Buffi 13
> Computational Sciences/Parrinello
> 6900 Lugano, Switzerland
> Tel : +41 (0) 58666-4802       Fax : +41 (0) 58666-4817
> --------------------------------------------------------------------
>
>
 

_______________________________________________
Pw_forum mailing list
Pw_forum at pwscf.org
http://www.democritos.it/mailman/listinfo/pw_forum 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20120828/baef3677/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: water-nmr_1_bad.log
Type: application/octet-stream
Size: 7232 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20120828/baef3677/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: H.pbe-tm_bad.in
Type: application/octet-stream
Size: 476 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20120828/baef3677/attachment-0001.obj>


More information about the users mailing list