[QE-users] Regarding initialization of elastic constants for NEB code

[Jana, Soumyajit]

Hi,

I am using NEB code to calculate vacancy jump activation energy in pure copper. I am not sure how to initialize elastic constants  ( k_max and k_min) for this system using Hartree atomic units ( hartree/ bohr^2). I tried a bunch of values but they all gave me weird results. I have also included my code here.

Any help will be greatly appreciated.

Thanks
Soumyajit


BEGIN
BEGIN_PATH_INPUT
&PATH
   restart_mode  = 'from_scratch'
   string_method  = 'neb'
   nstep_path     = 55
   num_of_images  = 9
   opt_scheme     = 'broyden'
   CI_scheme      = 'auto'
   path_thr       = 0.05
   k_max          = 0.3
   k_min          = 0.001
   ds             = 0.1
   first_last_opt = .false.
   minimum_image  = .true.
   use_freezing   = .true.
   use_masses     = .true.
/
END_PATH_INPUT

BEGIN_ENGINE_INPUT
&CONTROL
  prefix='copper_neb',
  pseudo_dir = './',
  outdir='./',
  verbosity='low',
/

&SYSTEM
   ibrav=0,
   A=7.66836091,
   B=7.66836091,
   C=7.66836091,
   cosAB=0.5,
   cosBC=0.5,
   cosAC=0.5,
   nat       = 26,
   ntyp      = 1,
   ecutwfc   = 50,
   ecutrho   = 250,
   occupations='smearing', smearing='marzari-vanderbilt', degauss=0.005,
/

&ELECTRONS
   conv_thr = 1d-7
   mixing_beta = 0.7,
/


CELL_PARAMETERS {alat}
   1.00   0.00   0.00
   0.00   1.00   0.00
   0.00   0.00   1.00

ATOMIC_SPECIES
   Cu  63.5  Cu_PBE_PAW_nonlinear.UPF
BEGIN_POSITIONS
FIRST_IMAGE
ATOMIC_POSITIONS {alat}
Cu   -0.3333   -0.3333   -0.3333
Cu   -0.3333   -0.3333   0.0000
Cu   -0.3333   -0.3333   0.3333
Cu   -0.3333   0.0000   -0.3333
Cu   -0.3333   0.0000   0.0000
Cu   -0.3333   0.0000   0.3333
Cu   -0.3333   0.3333   -0.3333
Cu   -0.3333   0.3333   0.0000
Cu   -0.3333   0.3333   0.3333
Cu   0.0000   -0.3333   -0.3333
Cu   0.0000   -0.3333   0.0000
Cu   0.0000   -0.3333   0.3333
Cu   0.0000   0.0000   -0.3333
#Cu   0.0000   0.0000   0.0000
Cu   0.0000   0.0000   0.3333
Cu   0.0000   0.3333   -0.3333
Cu   0.0000   0.3333   0.0000
Cu   0.0000   0.3333   0.3333
Cu   0.3333   -0.3333   -0.3333
Cu   0.3333   -0.3333   0.0000
Cu   0.3333   -0.3333   0.3333
Cu   0.3333   0.0000   -0.3333
Cu   0.3333   0.0000   0.0000
Cu   0.3333   0.0000   0.3333
Cu   0.3333   0.3333   -0.3333
Cu   0.3333   0.3333   0.0000
Cu   0.3333   0.3333   0.3333
LAST_IMAGE
ATOMIC_POSITIONS {alat}
Cu   -0.3333   -0.3333   -0.3333
Cu   -0.3333   -0.3333   0.0000
Cu   -0.3333   -0.3333   0.3333
Cu   -0.3333   0.0000   -0.3333
Cu   -0.3333   0.0000   0.0000
Cu   -0.3333   0.0000   0.3333
Cu   -0.3333   0.3333   -0.3333
Cu   -0.3333   0.3333   0.0000
Cu   -0.3333   0.3333   0.3333
Cu   0.0000   -0.3333   -0.3333
Cu   0.0000   -0.3333   0.0000
Cu   0.0000   -0.3333   0.3333
Cu   0.0000   0.0000   -0.3333
Cu   0.0000   0.0000   0.0000
#Cu   0.0000   0.0000   0.3333
Cu   0.0000   0.3333   -0.3333
Cu   0.0000   0.3333   0.0000
Cu   0.0000   0.3333   0.3333
Cu   0.3333   -0.3333   -0.3333
Cu   0.3333   -0.3333   0.0000
Cu   0.3333   -0.3333   0.3333
Cu   0.3333   0.0000   -0.3333
Cu   0.3333   0.0000   0.0000
Cu   0.3333   0.0000   0.3333
Cu   0.3333   0.3333   -0.3333
Cu   0.3333   0.3333   0.0000
Cu   0.3333   0.3333   0.3333


END_POSITIONS


K_POINTS {automatic}
   4 4 4  0 0 0

END_ENGINE_INPUT
END
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20191010/5379d6d9/attachment.html>