forum@abinit.org
Subject: The ABINIT Users Mailing List ( CLOSED )
List archive
- From: "Matteo Giantomassi" <Matteo.Giantomassi@uclouvain.be>
- To: forum@abinit.org
- Subject: Re: [abinit-forum] GW using HGH-PP?
- Date: Fri, 16 Jan 2009 10:38:10 +0100 (CET)
- Domainkey-signature: a=rsa-sha1; c=nofws; d=uclouvain.be; h=message-id: in-reply-to:references:date:subject:from:to:mime-version: content-type:content-transfer-encoding; q=dns; s=selucl; b=RBvVY Jdfgfw0+sqZey2wZsKXM900tpzZkmpjZWOMxlUu3PWpOcX4OcH6zpayW86b10Tl0 io+N2cTvk/0G/NmMUaaJWMy5Nx4r0ro9M0fJ0aKPkBi4ZHmhFjUpBxYzRXSanJSQ smGo6jR79AzjWxSQYtfXz2ydKmr5ukJX3FL3lA=
- Importance: Normal
> Dear Seungmi,
>
> I don't exactly if HGH pseudos are working in the GW code.
HGH pseudos contain additional projectors for the spin-orbit coupling term.
Since GW is not yet able to deal with non-collinear systems
HGH pseudos, at present, can only be used to generate a KSS file in the
collinear case employing kssform ==3
Please note however that, to the best of my knowledge, nobody has ever
used HGH to perform
GW calculations. I ran some tests and indeed the code stops complaining
about the number of
projectors found in the pseudopotential file.
Making the GW part work with HGH will require some change in the code as
well as
additional testing. For the time being I suggest use "standard"
norm-conserving pseudos
or the files available in the Goedecker-Teter-Hutter (GTH) table.
I agree with Fabien on the need of adding semicore states during the
generation of the
pseudo. Please keep in mind, however, that an accurate
pseudo with semicore states would require more than one projector per
angular state. If you don't find the elements you need in the GTH table,
you might use the opium generator to have more than one projector per a
particular l
(FHI doesn't support multiprojector pseudopotentials)
By the way:
As already explained in the forum, in case of pseudos with more than one
projector per angular state,
the treatment of the non analytic behavior of the inverse dielectric
matrix is not
coded. In this case one has to use inclvkb==0 during the screening run.
Hope it helps,
Best Regards,
Matteo Giantomassi
> But something I had repeated many times in this forum: when performing a
> GW calculation, one should always choose a pseudopotential, which
> considers the semicore states as valence. Here for indium, I think the
> 4s4p4d should be treated as valence, else the exact-exchange operator
> could be off by several eV.
> Eventually, I advise you to create your own pseudo...
>
> Cheers,
>
> Fabien
- [abinit-forum] GW using HGH-PP?, seungmi.lee, 01/14/2009
- Re: [abinit-forum] GW using HGH-PP?, Fabien Bruneval, 01/14/2009
- Re: [abinit-forum] GW using HGH-PP?, Matteo Giantomassi, 01/16/2009
- Re: [abinit-forum] GW using HGH-PP?, Philippe Blaise, 01/16/2009
- Re: [abinit-forum] GW using HGH-PP?, Matteo Giantomassi, 01/16/2009
- Re: [abinit-forum] GW using HGH-PP?, Philippe Blaise, 01/16/2009
- Re: [abinit-forum] GW using HGH-PP?, Matteo Giantomassi, 01/16/2009
- Re: [abinit-forum] GW using HGH-PP?, Fabien Bruneval, 01/14/2009
Archive powered by MHonArc 2.6.15.