My Community

Hi,

I was running NEB calculation to find activation barrier for H in Fe system. When i run the job, both NEB & CI-NEB is giving the same result.

NEB INCAR

# Basic setup:
ISPIN = 2
LSORBIT = .FALSE.
# not using Hubbard U

# Accuracy controls:
PREC = High
# Use default energy cutoff for PREC

# Electronic loop controls:
ALGO = Normal # Normal (Davidson)
EDIFF = 0.000100
NELMIN = 2
NELM = 60
TIME = 0.400000
ISMEAR = -1 # Fermi smearing
SIGMA = 0.200000 # eV
# NBANDS - using default

# Relaxation control:
IBRION = 1 # Quasi-Newton (RMM-DIIS)
NSW = 100
POTIM = 0.500000
ISIF = 2 # Ions
EDIFFG = -0.010000
# NEB controls:
IMAGES = 7
SPRING = -5.000000

CI-NEB INCAR
ISPIN = 2
LSORBIT = .FALSE.
# not using Hubbard U

# Accuracy controls:
PREC = High
# Use default energy cutoff for PREC

# Electronic loop controls:
ALGO = Normal # Normal (Davidson)
EDIFF = 0.000100
NELMIN = 2
NELM = 60
TIME = 0.400000
ISMEAR = -1 # Fermi smearing
SIGMA = 0.200000 # eV
# NBANDS - using default

# Relaxation control:
IBRION = 1 # Quasi-Newton (RMM-DIIS)
NSW = 100
POTIM = 0.500000
ISIF = 2 # Ions
EDIFFG = -0.010000
# NEB controls:
IMAGES = 7
SPRING = -5.000000
LCLIMB = .TRUE.
# Properties:
LCHARG = .FALSE.
LWAVE = .FALSE.
LELF = .FALSE.
LVTOT = .FALSE.
LVHAR = .FALSE.

So please comment.

Regards,
Anil

It seems that you have not set the required parameters in your INCAR,

IOPT = 0
IBRION = 3
POTIM = 0

It could be the same NEB calculation, not CI-NEB.