<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"></head><body><div>I think that if there are two unstable directions at the saddle point this means there is another saddle point with a lower </div><div>Activation energy nearby. Just slide the neb downhill</div><div> stefano </div><div><br></div><div><br></div><div id="composer_signature"><div style="font-size:85%;color:#575757" dir="auto">Sent from my Samsung Galaxy smartphone.</div></div><div><br></div><div style="font-size:100%;color:#000000"><!-- originalMessage --><div>-------- Original message --------</div><div>From: Pascal Boulet <pascal.boulet@univ-amu.fr> </div><div>Date: 13/11/2017  19:45  (GMT+01:00) </div><div>To: PWSCF Forum <pw_forum@pwscf.org> </div><div>Subject: *spam*Re: [Pw_forum] Physical meaning of two imaginary phonon frequencies at NEB saddle point and Transition State Theory </div><div><br></div></div>Hello,<div><br></div><div>My guess: you have to remove one of the imaginary frequencies (most probably the smallest one) by following the corresponding eigenvector (=moving the atoms accordingly). Then, you will have to make sure that the transition state connects the reactants and products you are interested in.</div><div><br></div><div>HTH,</div><div>Pascal</div><div><br><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="orphans: auto; text-align: start; text-indent: 0px; widows: auto; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><span style="font-size: 16pt; font-family: Mistral;">Pascal Boulet</span></div><div><font face="Lucida Handwriting"><span style="font-size: 21px;">—</span></font></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><em style="color: rgb(34, 187, 234); font-family: verdana; font-size: 11px; line-height: 22px;">Professor in computational chemistry - DEPARTEMENT OF CHEMISTRY</em></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span style="font-size: 11px; font-family: verdana;">Aix-Marseille University - </span><span style="font-size: 11px; font-family: verdana;">ST JEROME - Avenue Escadrille Normandie Niemen - F-13013 Marseille - FRANCE</span></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-family: Tahoma, Verdana, Arial, Helvetica, sans-serif; font-size: 15px; line-height: 22px;"><span style="font-size: 11px; font-family: verdana;">Tél: +33(0)4 13 55 18 10 - Fax : +33(0)4 13 55 18 50</span></div><div><font face="verdana" style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><span style="font-size: 12px; line-height: 22px;">Site : </span></font><font face="verdana"><font color="#22bbea"><span style="font-size: 12px; line-height: 20px;"><a href="http://madirel.univ-amu.fr/pages_web_BOULET_PASCAL/infos">http://madirel.univ-amu.fr/pages_web_BOULET_PASCAL/infos</a></span></font><span style="font-size: 12px; line-height: 22px;"> - Email : </span></font><a href="mailto:pascal.boulet@univ-amu.fr" style="color: rgb(34, 187, 234); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-family: verdana; font-size: 12px; line-height: 20px; text-decoration: none;">pascal.boulet@univ-amu.fr</a></div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></div></div><br class="Apple-interchange-newline">
</div>
<br><div><div>Le 13 nov. 2017 à 12:44, JAY Antoine <<a href="mailto:Antoine.JAY@isae-supaero.fr">Antoine.JAY@isae-supaero.fr</a>> a écrit :</div><br class="Apple-interchange-newline"><blockquote type="cite">

<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
Dear all,<br>I have been using NEB to obtain the Delta E and phonon DOS at begining/end/saddle point to properly evaluate the vibrationnal entropy contribution in the transition state theory.<br>Commonly, one frequency in "negative" at the saddle point, the one that is "responssible" of the transition and the following of its eigen vector from the initial state give the minimum energy path.<br>However, when two phonon frequencies are "negatives" at the saddle point, that means that at least 3  energy minima are link together to the same saddle point.<br>The transition state theory cannot be applied there as the preexponential term is the ratio of the 3N-3 phonon frequencies of initial state over the 3N-4 frequencies of the saddle point, giving the dimension of a time (state life time)<br>Using 3N-5 frequencies without the two negatives ones would not have a physical meaning for the initial partition function.<br><br>Any idea on how to treat such a problem to evaluate the transition rate would be appreciated.<br><br>Best regards.<br><br>Antoine Jay
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://pwscf.org/mailman/listinfo/pw_forum</blockquote></div><br></div></body></html>