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
Problem running NEB calculation
Moderators: Global Moderator, Moderator
-
- Newbie
- Posts: 2
- Joined: Thu Aug 09, 2012 9:09 am
- License Nr.: 1227
- Location: Technion City, Haifa, Israel
Problem running NEB calculation
Last edited by hadasabir on Sun Aug 12, 2012 9:32 am, edited 1 time in total.
-
- Full Member
- Posts: 105
- Joined: Tue Oct 12, 2010 11:10 pm
Problem running NEB calculation
Hi,
Is this your complete INCAR? What about the NEB parameters?
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.
-
- Newbie
- Posts: 2
- Joined: Thu Aug 09, 2012 9:09 am
- License Nr.: 1227
- Location: Technion City, Haifa, Israel
Problem running NEB calculation
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
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.
-
- Full Member
- Posts: 105
- Joined: Tue Oct 12, 2010 11:10 pm
Problem running NEB calculation
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?
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.