[SOLVED] BLYP pseudopotentials

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

Moderators: ebousquet, bxu

Locked
ilya
Posts: 3
Joined: Fri Feb 04, 2011 7:26 am

[SOLVED] BLYP pseudopotentials

Post by ilya » Fri Feb 04, 2011 7:51 am

Dear All,
Could you help me to find (or get some package to construct) BLYP spin-orbit pseudopotentials, e.g., for Tl, Pb, Bi, Sb, Te, Se? Thanks in advance.
Last edited by ilya on Sat Feb 05, 2011 8:41 pm, edited 1 time in total.

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

Re: BLYP pseudopotentials

Post by jzwanzig » Fri Feb 04, 2011 12:12 pm

You will need abinit built with plug-in libxc support, and also atompaw built with libxc support. Then you can construct PAW datasets for your elements using B86 or B88 exchange and LYP correlation (that gives BLYP altogether). Then when you run calculations you can use the pawspnorb option to include at least on-site spin-orbit interaction in the PAW calculations.
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

ilya
Posts: 3
Joined: Fri Feb 04, 2011 7:26 am

Re: BLYP pseudopotentials

Post by ilya » Fri Feb 04, 2011 1:04 pm

Thanks a lot. I will try to realize such a scheme.

Mutta
Posts: 69
Joined: Fri Feb 19, 2010 11:33 am

Re: [SOLVED] BLYP pseudopotentials

Post by Mutta » Thu Nov 17, 2011 6:17 am

Dear sir
After I generated with atompaw which construct PAW datasets for using B86 or B88 exchange and LYP correlation but abinit not support
and Can I used ixc=18 for internally for GGA BLYP pseudopotentials from M. Krack, see Theor. Chem. Acc. 114, 145 (2005), available from the CP2K repository - use the LibXC instead, with ixc=-106131
see my input files below
# =======================================
# CuBO2: computation of the optimal lattice parameter
# =======================================

optcell 2

#Optimization of the lattice parameters
ionmov 2
ntime 1000
dilatmx 1.1
ecutsm 0.5

##Definition of the unit cell
acell 3*5.7456 angstrom
brvltt 7
spgaxor 2
spgroup 166
angdeg 3*28.62
ixc -106131

#Definition of the atoms
ntypat 3
znucl 29 5 8
natom 4
typat 1 2 3 3
xred
0.0 0.0 0.0
0.5 0.5 0.5
1/9 1/9 1/9
8/9 8/9 8/9

#Definition of the k-point grid
ecut 30
#pawecutdg 35
kptopt 1
ngkpt 8 8 8
nshiftk 1
shiftk 0.5 0.5 0.5

#Definition of the SCF procedure
nstep 1000
toldfe 1.0d-6

thank
Abdulmutta

Locked