<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">

</head>

<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

Hi all,

<div><br>

</div>

<div>I am using QE 5.0.3.</div>

<div><br>

</div>

<div>I am trying to test if a saddle point from the neb.x calculation is a real saddle point. I am doing this by looking for the vibrational frequency: If a state is a saddle point, it has only one imaginary (negative) frequency. If a state is an energy minimum,

 it has no imaginary frequency. </div>

<div><br>

</div>

<div>My system consists of a crystal structure and a helium impurity in it. The crystal structures at energy minimum and the saddle point are almost identical, only the helium atom is displaced. To facilitate the frequency calculation, I used nat_todo = 1 to

 only allowing the helium to be displaced: </div>

<div><br>

</div>

<div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

phonons of initial state at Gamma</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

 &inputph</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

  nat_todo = 1</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

  tr2_ph=1.0d-14,</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

  prefix = "te",</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

  epsil=.true.,</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

 /</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

0.0 0.0 0.0</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

25</div>

<div><br>

</div>

<div>The last line, 25, indicates only the 25th atom (helium) is displaced during the ph.x calculation. Because there are 25 atoms total in the system, there should be 75 modes. The complete ph.x outputs are attached. The energy minimum state has the frequencies: </div>

<div><br>

</div>

<div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

     omega( 1) =      -2.080553 [THz] =     -69.399779 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

     omega( 2) =      -2.080553 [THz] =     -69.399779 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega( 3) =      -2.080553 [THz] =     -69.399779 [cm-1]</div>

</div>

<div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

<div style="margin: 0px; position: static; z-index: auto;">ЕЕЕ.</div>

<div style="margin: 0px; position: static; z-index: auto;">ЕЕЕ.</div>

</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega(73) =       4.621913 [THz] =     154.170425 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega(74) =       4.621913 [THz] =     154.170425 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega(75) =       4.621913 [THz] =     154.170425 [cm-1]</div>

</div>

<div><br class="webkit-block-placeholder">

</div>

<div>with omega(4~72) all equal to zero. For the saddle point one:</div>

<div><br>

</div>

<div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega( 1) =      -4.316979 [THz] =    -143.998928 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega( 2) =      -2.270462 [THz] =     -75.734475 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega( 3) =      -2.270462 [THz] =     -75.734475 [cm-1]</div>

</div>

<div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

<div style="margin: 0px; position: static; z-index: auto;">ЕЕЕ.</div>

<div style="margin: 0px; position: static; z-index: auto;">ЕЕЕ.</div>

</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega(73) =       2.168375 [THz] =      72.329209 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">

     omega(74) =       9.304582 [THz] =     310.367434 [cm-1]</div>

<div style="margin: 0px; font-size: 16px; font-family: 'Lucida Console'; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); position: static; z-index: auto;">

     omega(75) =       9.304582 [THz] =     310.367434 [cm-1]</div>

</div>

<div><br class="webkit-block-placeholder">

</div>

<div><br class="webkit-block-placeholder">

</div>

<div>My questions are:</div>

<div>1. Why are there 3 additional frequencies? ShouldnТt there only be 3 (omega 73~75) due to the DoF of one atom? </div>

<div>2. Why are the three additional frequencies negative? ShouldnТt there be one negative frequency for the saddle point, and 2 positive frequencies? and no negative frequency for the energy minimum?</div>

<div><br>

</div>

<div>Please let me know if I understand this incorrectly, or I am missing some settings.</div>

<div><br>

</div>

<div><br>

</div>

<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="color: rgb(0, 0, 0); font-family: Helvetica;  font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div style="font-variant: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div style="font-variant: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-weight: normal; font-style: normal; font-variant: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

--</div>

<div style="font-variant: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

Kan-Ju Lin, Postdoctoral Associate<br>

<font color="#a24630"><b>MIT</b></font><span style="color: rgb(0, 0, 0);"> Materials Science and Engineering<br>

Office 4-053 / </span><a href="mailto:kjlin@mit.edu">kjlin@mit.edu</a><span style="color: rgb(0, 0, 0);"><br>

<br>

<br>

<br>

</span></div>

</div>

</div>

</div>

</div>

</div>

</div>

</div>

<br>

</div>

</body>

</html>