Problem running NEB calculation

Queries about input and output files, running specific calculations, etc.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
hadasabir
Newbie
Newbie
Posts: 2
Joined: Thu Aug 09, 2012 9:09 am
License Nr.: 1227
Location: Technion City, Haifa, Israel

Problem running NEB calculation

#1 Post by hadasabir » Sun Aug 12, 2012 9:32 am

Hi all
I want to run calculation for CH4 dissociation on transition metals
I use ase simulation program
after several iteration I get the following error:
Error EDDDAV: Call to ZHEGV failed. Returncode = 80 532

my incar is:

ENCUT = 400
PREC = med
SYMPREC = 1e-6
POTIM = 1
IBRION = 3
LREAL = False
GAMMA= True

I will be glad for any help
Thanks
Hadas
Last edited by hadasabir on Sun Aug 12, 2012 9:32 am, edited 1 time in total.

kambiz
Full Member
Full Member
Posts: 105
Joined: Tue Oct 12, 2010 11:10 pm

Problem running NEB calculation

#2 Post by kambiz » Mon Aug 13, 2012 3:02 am

Hi,
Is this your complete INCAR? What about the NEB parameters?
Last edited by kambiz on Mon Aug 13, 2012 3:02 am, edited 1 time in total.

hadasabir
Newbie
Newbie
Posts: 2
Joined: Thu Aug 09, 2012 9:09 am
License Nr.: 1227
Location: Technion City, Haifa, Israel

Problem running NEB calculation

#3 Post by hadasabir » Mon Aug 13, 2012 8:46 am

Hi I am using the ase simulation program
my complete goes something like this:
from ase import *
from ase.optimize.fire import FIRE

initial = io.read('CH4.traj')
final = io.read('CH3_H.traj')

images = [initial]
images += [initial.copy() for i in range(3)]
images += [final]
neb = neb.NEB(images)

mask = [atom.symbol == 'Ni' for atom in initial]
constraint = constraints.FixAtoms(mask=mask)

neb.interpolate()

calcs = []
for i in range(3):
calc =calculators.vasp.Vasp(prec = 'Med', lreal = False,kpts=(3,3,1),gamma=True, symprec=1e-6, ispin=2, encut=400)
calcs.append(calc)

for i,image in enumerate(images[1:4]):
image.set_calculator(calcs)
image.set_constraint(constraint)

optimizer= FIRE(neb, trajectory='neb.traj')
optimizer.run(fmax=0.1)

*****************
Thanks
Last edited by hadasabir on Mon Aug 13, 2012 8:46 am, edited 1 time in total.

kambiz
Full Member
Full Member
Posts: 105
Joined: Tue Oct 12, 2010 11:10 pm

Problem running NEB calculation

#4 Post by kambiz » Tue Aug 14, 2012 12:00 pm

I have never used this program. I use Vasp TST Tools.
Have you tried to run your NEB in VASP without using any other packages?
Last edited by kambiz on Tue Aug 14, 2012 12:00 pm, edited 1 time in total.

Post Reply