Identical inputs, different results at different cpu cluster

Problems running VASP, runtime errors, technical questions.
Please check whether the answer to your question is given in
the VASP online manual or has been discussed in this forum
previously!

Moderators: Moderator, Global Moderator

Identical inputs, different results at different cpu cluster

Postby smattsson » Wed May 09, 2018 12:07 pm

Hello everyone

I am currently calculating the bulk structure of CoF3 using PBE+U according to the Dudarev formulation. The cell is a rhombohedral one with two units of CoF3 in it. I have optimized the structure at my computer cluster #1 and wanted to do a frequency calculation at another cluster #2 for reasons of convenience. But when I use the optimized structure as the POSCAR and run a single-point SCF calculation at cluster #2, the E0 is higher: -.33626837E+02 at #2 vs. the reference calculation at #1: -.35146116E+02 - i.e. a difference of a whole 1.5 eV. By looking at the DOS, the #2 calculation falsely describes the system as a conductor. #1 gives a band gap of 2.8 eV.

I tested the two single-point calculations using identical input files INCAR, KPOINTS, POSCAR and POTCAR (by copy-pasting them). The only difference is the computer cluster they are running on and the core/node distribution. Does anyone have an idea what might be going wrong here? I am thankful for any replies

Best regards,
Stefan

Cpu details:
#1: aprun protocol, 8 nodes with 8 cores each
#2: mpirun protocol, 2 nodes with 8 cores each

Both systems run vasp 5.4.1 (vasp_std executable, i.e. non-collinear)


INCAR:
Code: Select all
SYSTEM = CoF3 PBE
NCORE = 8

ISTART = 0
LWAVE = .FALSE.

ISMEAR = 0 ! -5=linear tetrahedron (semiconductor/insulator), 0=gaussian
SIGMA = 0.10 ! temperature

PREC = Medium

!NSW = 100 ! ionic steps
IBRION = -1 ! 1=use quasi-newton optimization 2=conjugate gradient, -1=no optimization
ISIF = 3 ! 3=relax ions, cell shape and cell volume
ENMAX = 850 ! use together with ISIF
!ISYM = -1
POTIM = 0.25

EDIFF = 1.0E-5
EDIFFG = -0.01

ENCUT = 700
ICHARG = 2 ! 1=restart, 2=from atom solutions, 11=don't change density

ISPIN =  2 ! 2=spin polarized calculation
LORBIT = 11
MAGMOM = 4 -4 0 0 0 0 0 0

LAECHG = .TRUE.


LDAU = .TRUE.
LDAUTYPE = 2 ! 2=Dudarev
LDAUL = 2 -1 ! angular quantum number of orbitals 1=p, 2=d, 3=f
LDAUU = 6.3 0.0
LDAUJ = 1.0 0.0
LDAUPRINT = 2

LMAXMIX = 4 ! use with Hubbard


KPOINTS:
Code: Select all
K-Points
 0
Monkhorst Pack
 8 8 8
 0  0  0


POSCAR:
Code: Select all
CoF3                                   
   1.00000000000000     
     2.9515852336660924    0.0000000000000000    4.4095306744468550
    -1.4757926168330462    2.5561477937949211    4.4095306744468550
    -1.4757926168330462   -2.5561477937949211    4.4095306744468550
   Co   F
     2     6
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.5000000000000000
  0.8486837912520073  0.6513162087479927  0.2500000000000000
  0.2500000000000000  0.8486837912520073  0.6513162087479927
  0.6513162087479927  0.2500000000000000  0.8486837912520073
  0.1513162087479927  0.3486837912520073  0.7500000000000000
  0.7500000000000000  0.1513162087479927  0.3486837912520073
  0.3486837912520073  0.7500000000000000  0.1513162087479927
 
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00


I spare you the POTCAR but it is constructed from the standard PBE PAW's
smattsson
Newbie
Newbie
 
Posts: 4
Joined: Tue Sep 27, 2016 2:02 pm

Re: Identical inputs, different results at different cpu clu

Postby admin » Fri Aug 31, 2018 10:34 am

sorry, but we cannot deal with this issue based on the information you give
we need a FULL error report of both runs (including the OUTCAR,OSZICAR,job.out and job.err files)
admin
Administrator
Administrator
 
Posts: 2909
Joined: Tue Aug 03, 2004 8:18 am


Return to Bugreports

Who is online

Users browsing this forum: No registered users and 1 guest