[Pw_forum] Transition state calculation error

Cecil Humphrey Botchway humphrey2g4 at gmail.com
Tue May 3 14:48:02 CEST 2016


Hi There,
I am trying to run a transition state calculation using neb.x with methanol
C-OH bond breaking in the zeolite micro pore (the system is composed of 120
atoms). I have tried relaxing the whole structure with the adsorbate and
freezing only the zeolite leaving the adsorbate. I have even tried
increasing and decreasing the numbers of processor cores from 8-16. I
encounter the same error (which is below) regardless .


 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     task #         0
     from rdiaghg : error #       818
     S matrix not positive definite
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


What could be the possible cause?
This is how my input looks like

BEGIN
BEGIN_PATH_INPUT
&PATH
  restart_mode      = 'from_scratch'
  string_method     = 'neb',
  nstep_path        = 20,
  ds                = 2.D0,
  opt_scheme        = "broyden",
  num_of_images     = 7,
  k_max             = 0.3D0,
  k_min             = 0.2D0,
  CI_scheme         = "auto",
  path_thr          = 0.1D0,
/
END_PATH_INPUT
BEGIN_ENGINE_INPUT
&CONTROL
  prefix         = "TS1"
  outdir         = "./tmp/",
  pseudo_dir     = "/home/cbotchway/Cecil",
/
&SYSTEM
  ibrav                  = 0,
  nat                    = 120,
  ntyp                   = 5,
  ecutwfc                = 40.0D0,
  ecutrho                = 240.0D0,
  nspin                  = 2,
  starting_magnetization = 0.5D0,
  occupations            = "smearing",
  degauss                = 0.03D0,
/
&ELECTRONS
  conv_thr    = 1.D-6,
  mixing_beta = 0.2D0,
  electron_maxstep =1000000,
/
&IONS
/
ATOMIC_SPECIES
O    15.99940    O.pbe-van_ak.UPF
Si   28.0855    Si.pbe-n-van.UPF
Al   26.98154   Al.pbe-n-van.UPF
H     1.00794    H.pbe-van_ak.UPF
C    12.0107     C.pbe-van_ak.UPF
BEGIN_POSITIONS
FIRST_IMAGE

Attached to this post is a copy of my input file
Any help or advice would be appreciated
Thank you in advance
-- 
Cecil Huymphrey Botchway
Mphil. Material Science and Computational Chemistry
Department of Chemistry
Kwame Nkrumah University of Science and Technology
+23329002372
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20160503/2383f34a/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FERmepos2p.in
Type: application/octet-stream
Size: 42806 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20160503/2383f34a/attachment.obj>


More information about the users mailing list