[Pw_forum] NEB Calculation with Ti and Graphene
Ettore Baldini Neto
neto.baldini at gmail.com
Fri Dec 17 23:15:58 CET 2010
Dear users,
I'm playing with a NEB calculation involving one Ti atom in a graphene monolayer.
The calculation crashes due to diagonalization problems according to the message below.
*******************************************
task # 8
from cdiaghg : error # 434
diagonalization (ZHEGV*) failed
********************************************
In what follows I copied my input file (without the atomic positions in order to make it small for the readers)
Any hint? Am I doing something very wrong?
I'm not using for now polarization for Ti for now in order to make things simpler. I'm planning to include it soon.
Also, I'm using a large ecut but I figured out that for graphene and for this pseudo 35-40Ry might work fine.
title = 'Graphene mono layer'
calculation = "neb" ,
restart_mode = 'from_scratch' ,
outdir='/home/usuarios/baldini/espresso-4.2.1/tmp/' ,
pseudo_dir = '/home/usuarios/baldini/espresso-4.2.1/pseudo' ,
etot_conv_thr = 1.0E-3 ,
forc_conv_thr = 1.0D-4 ,
nstep=10,
prefix='1MLG_Ti_NEB',
/
&SYSTEM
ibrav = 0.,
ntyp = 2 ,
nat = 33 ,
ecutwfc = 50.0 ,
ecutrho = 400.0,
occupations = 'smearing' ,
degauss = 0.02,
smearing = 'mv' ,
/
&ELECTRONS
conv_thr = 1.0d-7 ,
mixing_beta = 0.4 ,
/
&IONS
num_of_images = 4,
opt_scheme = "broyden",
path_thr = 0.03D0,
/
ATOMIC_SPECIES
C 12.0107 C.pbe-van_ak.UPF
Ti 47.8670 Ti.pbe-sp-van_ak.UPF
CELL_PARAMETERS
16.12131954 9.383047808 0.0000000
0.00000000 18.658162370 0.0000000
0.0000000 0.000000000 56.6893424
ATOMIC_POSITIONS angstrom
first_image
...
last image
.......
K_POINTS {automatic}
6.0 6.0 1.0 1.0 1.0 1.0
Thanks a lot and best regards
************************************************************************
Dr. Ettore Baldini-Neto,
Researcher at Wernher von Braun Center of Advanced Research
Campinas, Brazil
***********************************************************************
More information about the users
mailing list