A question of using NEB
Posted: Sun Mar 29, 2009 2:21 pm
Hi,
I ran my NEB calculations using 5 images with one image for one cluster node. I found that only one image is generating output files, such as, OSZICAR, OUTCAR and so on. I want to make sure first that I set up the calculations correctly. So here is my input file (INCAR):
SYSTEM = NEB
PREC = Normal !Low; MEDIUM; High; Normal; Accurate
ENMAX = 280 !Cutoff set manually (preferred)
ISTART = 0 !0:from scratch; 1:restart w/ const E-cutoff; 2:restart w/ const basis-set; 3:full restart
ICHARG = 2 !0:from initial wf; 1:from file CHGCAR; 2:superposition of charg den
NELM = 60 !Maximum electronic SCF steps
EDIFF = 1E-5 !Allowed error in total energy
EDIFFG = -0.01 !Break condition for ion relaxation
ISMEAR = 2 !-5:tetra w/ Blochl; -4:tetra w/o Blochl; -3:loop; -2:fixed partial occu; -1:Fermi; 1(N):MP method(metal); 0:Gaussian
SIGMA = 0.2 !Smearing parameter 0.2 recommended for MP method
IBRION = 1 !-1:Ion not moved; 0:MD; 1:quais-Newton; 2:CG; 3:damping factor supplied
NSW = 200 !Ionic steps
POTIM = 0.05 !Step parameter: 0.02 A stepwidth
ALGO = V !Normal; VeryFast; Fast; All; Damped
LVTOT = .FALSE. !Write total local potential, default: FALSE
STRING = -5 !String strength
IMAGES = 5 !Number of NEB intermediate images, excluding the endpoint.
ICHAIN = 0 !TS search option.0:NEB
LCLIMB = .TRUE. !Flag for climbing image algorithm.
Also, the KPOINTS:
K-Points
0
Monkhorst-Pack
3 3 1
0 0 0
Thanks a lot for your prompt help.
B
I ran my NEB calculations using 5 images with one image for one cluster node. I found that only one image is generating output files, such as, OSZICAR, OUTCAR and so on. I want to make sure first that I set up the calculations correctly. So here is my input file (INCAR):
SYSTEM = NEB
PREC = Normal !Low; MEDIUM; High; Normal; Accurate
ENMAX = 280 !Cutoff set manually (preferred)
ISTART = 0 !0:from scratch; 1:restart w/ const E-cutoff; 2:restart w/ const basis-set; 3:full restart
ICHARG = 2 !0:from initial wf; 1:from file CHGCAR; 2:superposition of charg den
NELM = 60 !Maximum electronic SCF steps
EDIFF = 1E-5 !Allowed error in total energy
EDIFFG = -0.01 !Break condition for ion relaxation
ISMEAR = 2 !-5:tetra w/ Blochl; -4:tetra w/o Blochl; -3:loop; -2:fixed partial occu; -1:Fermi; 1(N):MP method(metal); 0:Gaussian
SIGMA = 0.2 !Smearing parameter 0.2 recommended for MP method
IBRION = 1 !-1:Ion not moved; 0:MD; 1:quais-Newton; 2:CG; 3:damping factor supplied
NSW = 200 !Ionic steps
POTIM = 0.05 !Step parameter: 0.02 A stepwidth
ALGO = V !Normal; VeryFast; Fast; All; Damped
LVTOT = .FALSE. !Write total local potential, default: FALSE
STRING = -5 !String strength
IMAGES = 5 !Number of NEB intermediate images, excluding the endpoint.
ICHAIN = 0 !TS search option.0:NEB
LCLIMB = .TRUE. !Flag for climbing image algorithm.
Also, the KPOINTS:
K-Points
0
Monkhorst-Pack
3 3 1
0 0 0
Thanks a lot for your prompt help.
B