<div dir="ltr"><div>This is unfortunately a nasty error that is not easy to solve. It is either due to pseudopotentials or to really bad structures. If the problem arises in the first NEB step, verify that your starting positions do not have atoms too close<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 3, 2016 at 2:48 PM, Cecil Humphrey Botchway <span dir="ltr"><<a href="mailto:humphrey2g4@gmail.com" target="_blank">humphrey2g4@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi There,<br></div>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 .<br> <br><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     task #         0<br>     from rdiaghg : error #       818<br>     S matrix not positive definite<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br></div><div><br><br>What could be the possible cause?<br></div>This is how my input looks like<br><br>BEGIN<br>BEGIN_PATH_INPUT<br>&PATH<br>  restart_mode      = 'from_scratch'<br>  string_method     = 'neb',<br>  nstep_path        = 20,<br>  ds                = 2.D0,<br>  opt_scheme        = "broyden",<br>  num_of_images     = 7,<br>  k_max             = 0.3D0,<br>  k_min             = 0.2D0,<br>  CI_scheme         = "auto",<br>  path_thr          = 0.1D0,<br>/<br>END_PATH_INPUT<br>BEGIN_ENGINE_INPUT<br>&CONTROL<br>  prefix         = "TS1"<br>  outdir         = "./tmp/",<br>  pseudo_dir     = "/home/cbotchway/Cecil",<br>/<br>&SYSTEM<br>  ibrav                  = 0,<br>  nat                    = 120,<br>  ntyp                   = 5,<br>  ecutwfc                = 40.0D0,<br>  ecutrho                = 240.0D0,<br>  nspin                  = 2,<br>  starting_magnetization = 0.5D0,<br>  occupations            = "smearing",<br>  degauss                = 0.03D0,<br>/<br>&ELECTRONS<br>  conv_thr    = 1.D-6,<br>  mixing_beta = 0.2D0,<br>  electron_maxstep =1000000,<br>/<br>&IONS<br>/<br>ATOMIC_SPECIES<br>O    15.99940    O.pbe-van_ak.UPF<br>Si   28.0855    Si.pbe-n-van.UPF<br>Al   26.98154   Al.pbe-n-van.UPF<br>H     1.00794    H.pbe-van_ak.UPF<br clear="all">C    12.0107     C.pbe-van_ak.UPF<br>BEGIN_POSITIONS<br>FIRST_IMAGE<br> <br>Attached to this post is a copy of my input file<br>Any help or advice would be appreciated<br>Thank you in advance<span class="HOEnZb"><font color="#888888"><br>-- <br><div><div dir="ltr"><div><div><div>Cecil Huymphrey Botchway<br></div>Mphil. Material Science and Computational Chemistry <br></div>Department of Chemistry<br></div>Kwame Nkrumah University of Science and Technology<br>+23329002372<br></div></div>
</font></span></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>