<div dir="ltr">Thank you Travis,<div><br></div><div>I didn't know that at all. Re-ordering indeed solved the problem.</div><div><br></div><div>Have a good day,</div><div><br></div><div>Jacob<br><br>On Tuesday, February 4, 2020 at 2:00:01 PM UTC+1, Travis wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>The atoms have to have the same ordering in <a href="http://i.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;">i.xyz</a> and <a href="http://f.xyz" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a>. You have NCH in one and NHC in the other.</div><div><br></div><div>-T</div><div><br><br>On Tuesday, February 4, 2020 at 7:37:46 AM UTC-5, Jacob Jensen wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><font face="arial, sans-serif">Dear all,</font><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I'm doing CI-NEB calculations for a set of reactions on different metal surfaces. During my calculations I encountered a very bizarre behavior of the code which I haven't been able to find a reason for it.</font></div><div><font face="arial, sans-serif">The problem is: I give starting (<a href="http://i.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;">i.xyz</a>) and end points (<a href="http://f.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a>) to the CP2K to do a simple CI-NEB calculation, but in a completely random way (only happens to some reactions) the code exchanges the last two lines of the final configuration (<a href="http://f.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a>) and starts the calculation with wrong coordinates. It has happened also for more than two lines in some structures.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">For instance in the coordinates provided below, "tail <a href="http://f.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a>" shows this:</font></div><div><font face="arial, sans-serif"><br></font></div><div><div><font face="arial, sans-serif"> Cu         2.5607522840        3.6179848259       10.2114493134</font></div><div><font face="arial, sans-serif"> Cu         5.1090632775        3.6199416787       10.2213122857</font></div><div><font face="arial, sans-serif"> Cu         7.6670721044        3.6155408617       10.2040766974</font></div><div><font face="arial, sans-serif"> Cu         0.0000430330        7.2313188455       10.2023231891</font></div><div><font face="arial, sans-serif"> Cu         2.5586830958        7.2299962241       10.2055589609</font></div><div><font face="arial, sans-serif"> Cu         5.1117660386        7.2291476249       10.2075861138</font></div><div><font face="arial, sans-serif"> Cu         7.6675778217        7.2313677337       10.2036234566</font></div><div><font face="arial, sans-serif">  N         6.0791607429        5.3490729602       14.3227064450</font></div><div><span style="background-color:rgb(255,255,0)"><font face="arial, sans-serif">  C         5.1523694228        5.0243296300       14.0410514879</font></span></div><div><span style="background-color:rgb(255,255,0)"><font face="arial, sans-serif">  H         4.1418019657        4.6566468208       13.6463880308</font></span></div></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">But once the job starts to run, CP2K exchanges the last two atoms (C and H) and starts with: </font></div><div><font face="arial, sans-serif"><br></font></div><div><div><font face="arial, sans-serif"> NEB|*************************<wbr>******************************<wbr>********************</font></div><div><font face="arial, sans-serif"> NEB| Geometry for Replica Nr. 2 in Angstrom</font></div><div><font face="arial, sans-serif"> NEB|         Cu         1.278231000    1.807692000    1.278231000</font></div></div><div><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif">.</font></div><div><div><font face="arial, sans-serif"> NEB|         Cu         5.112380330   -0.000210681   10.202471653</font></div><div><font face="arial, sans-serif"> NEB|         Cu         7.667615495   -0.000585929   10.201760866</font></div><div><font face="arial, sans-serif"> NEB|         Cu         0.000556842    3.615646368   10.203009475</font></div><div><font face="arial, sans-serif"> NEB|         Cu         2.560752284    3.617984826   10.211449313</font></div><div><font face="arial, sans-serif"> NEB|         Cu         5.109063277    3.619941679   10.221312286</font></div><div><font face="arial, sans-serif"> NEB|         Cu         7.667072104    3.615540862   10.204076697</font></div><div><font face="arial, sans-serif"> NEB|         Cu         0.000043033    7.231318845   10.202323189</font></div><div><font face="arial, sans-serif"> NEB|         Cu         2.558683096    7.229996224   10.205558961</font></div><div><font face="arial, sans-serif"> NEB|         Cu         5.111766039    7.229147625   10.207586114</font></div><div><font face="arial, sans-serif"> NEB|         Cu         7.667577822    7.231367734   10.203623457</font></div><div><font face="arial, sans-serif"> NEB|          N         6.079160743    5.349072960   14.322706445</font></div><div><span style="background-color:rgb(255,255,0)"><font face="arial, sans-serif"> NEB|          H         5.152369423    5.024329630   14.041051488</font></span></div><div><span style="background-color:rgb(255,255,0)"><font face="arial, sans-serif"> NEB|          C         4.141801966    4.656646821   13.646388031</font></span></div></div><div><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif">.</font></div><div><font face="arial, sans-serif">This happens without any error or warning and as you can imagine, the result of this calculation would be nothing but garbage. I have provided my input file, initial and final configurations together with the required basis sets and pseudopotentials to provide a possibility of simple test run (even with one core, cause the interpolation is at the very beginning of the output file).</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I would be grateful if you can take a look and let me know about the possible reason. Or if you have had the same experience with other structures that would be beneficial for me to know.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">I have tried with different versions of CP2K (e.g. 5.1 and 7.0) and the result is the same. Also, whenever I change <a href="http://i.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;">i.xyz</a> to <a href="http://f.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a> in :</font></div><div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"> &TOPOLOGY</font></div><div><font face="arial, sans-serif">                        COORD_FILE_FORMAT XYZ</font></div><div><font face="arial, sans-serif">                        COORD_FILE_NAME <a href="http://i.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;">i.xyz</a></font></div><div><font face="arial, sans-serif"> &END,</font></div></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">the same exchange of lines happens to <a href="http://i.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fi.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNG8E7ipZiX0LOe-snzV9VYr6iLytQ';return true;">i.xyz</a> and final replica (<a href="http://f.xyz" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Ff.xyz\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEgwd9kPFExbLobIdhccdRrAtDoaQ';return true;">f.xyz</a>) stays intact! This really baffles me a lot. </font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Any remark or sharing of experience is highly appreciated.</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Kind regards,</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Jacob Jensen</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"><br></font></div><div><br></div></div></blockquote></div></div></blockquote></div></div>