[Pw_forum] forces in full LDA+U scheme are not yet implemented
Raj kamal
rajkamalloyola at gmail.com
Sat Jan 23 10:38:58 CET 2016
dear QE experts
i trying to relax Limnsio4 structure with GGA + U. it shows
error ,,,,,,from force_hub : error # 1
forces in full LDA+U scheme are not yet implemented. please help me
to solve this problem.thanks in advance.
this is my input.
&CONTROL
calculation = 'relax' ,
outdir = '/home/kathirvel/espresso/tmp/' ,
pseudo_dir = '/home/kathirvel/espresso/pseudo/' ,
prefix = 'Limnsio4' ,
etot_conv_thr = 1.0D-3 ,
forc_conv_thr = 1.0D-2 ,
tprnfor=.TRUE.
/
&SYSTEM
ibrav = 8,
A = 6.31850,
B = 5.38650,
C = 4.97540,
cosBC = 0,
cosAC = 0,
cosAB = 0,
nat = 16,
ntyp = 4,
ecutwfc =30.0 ,
occupations='smearing',
smearing='mv',
degauss=0.01,
lda_plus_u = .TRUE.
lda_plus_u_kind=1
Hubbard_U(3)= 6
Hubbard_J(1,3)=1
/
&ELECTRONS
mixing_mode = 'local-TF' ,
mixing_beta = 0.7,
diagonalization = 'david' ,
conv_thr = 1.0e-6,
/
&IONS
/
ATOMIC_SPECIES
Li 6.941 Li.pbe-n-van.UPF
O 15.999 O.pbe-rrkjus.UPF
Mn 54.938 Mn.pbe-sp-van.UPF
Si 28.085 Si.pbe-rrkj.UPF
ATOMIC_POSITIONS angstrom
Li 1.59321 1.76467 4.84987
O 1.35450 3.71566 4.40074
Mn 3.15925 4.43040 4.87057
O 3.15925 0.82430 4.07092
Si -0.00000 4.52261 4.90236
O -0.00000 0.70752 4.35596
Li 1.56604 3.62183 2.36217
Li 4.75246 3.62183 2.36217
Li 4.72529 1.76467 4.84987
O 1.80475 1.67084 1.91304
O 4.51375 1.67084 1.91304
O 4.96400 3.71566 4.40074
Mn -0.00000 0.95610 2.38287
O -0.00000 4.56220 1.58322
Si 3.15925 0.86389 2.41466
O 3.15925 4.67898 1.86826
K_POINTS automatic
5 5 5 0 0 0
--
*Best regards,*
*Rajkamal.A.*
*Research Scholar,(SRM UNIV).*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20160123/48fa9792/attachment.html>
More information about the users
mailing list