Hi Tone,<br><br>Thanks for correcting me............<br>Is the linear interpolation in the Modules/path_base.f90 file used by the code only for finding the initial guess of the path (if not specified in the input). Can you please refer me to the specific subroutine where the interpolation is done in order to generate the data for the *.int file?<br>
<br>Prasenjit.<br><br><div class="gmail_quote">2008/6/12 Tone Kokalj <<a href="mailto:tone.kokalj@ijs.si">tone.kokalj@ijs.si</a>>:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="Ih2E3d">On Thu, 2008-06-12 at 17:19 +0200, Prasenjit Ghosh wrote:<br>
> As far as I know, you don't have an option to choose the type of<br>
> interpolation to obtain the data in the *.int file. The code uses a<br>
> linear interpolation technique to obtain the above mentioned data.<br>
<br>
</div>This is not true. As far as I know, the code uses the interpolation as<br>
suggested in the original NEB paper. This interpolation is kind of cubic<br>
interpolation, which uses the information also from the forces (=slope)<br>
provided by the code. According to my experience it's rather good (see<br>
the appendix in G. Henkelman, J. Chem. Phys. 113, 9978 (2000)).<br>
<br>
Regards, Tone<br>
<font color="#888888">--<br>
Tone Kokalj <<a href="mailto:tone.kokalj@ijs.si">tone.kokalj@ijs.si</a>><br>
J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia (tel:<br>
+386-1-477-3523 // fax:+386-1-477-3822)<br>
</font><div><div></div><div class="Wj3C7c"><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>PRASENJIT GHOSH,<br>POST-DOC,<br>ROOM NO: 265, MAIN BUILDING,<br>CM SECTION, ICTP,<br>STRADA COSTERIA 11,<br>TRIESTE, 34104,<br>ITALY<br>PHONE: +39 040 2240 369 (O)<br>
+39 3807528672 (M)