<div dir="ltr">Dear S. Ling,<br><br>I still don't see this problem, although the vibrational analysis was/should be performed only for the calculation of vibrational frequencies of water molecules in the system, not that of the crystal system. It is like as I would measure and obtain a vibrational spectrum of water at interface. <br>Maybe the solution is to turn off the smearing method ?<br>I am new in CP2K and would appreciate some hints how to write input code enabling the calculation of vibrational frequencies of the interfacial water.<br><br>Kind regards,<br><br>Phil<br><br>Am Mittwoch, 25. Januar 2017 14:51:21 UTC+1 schrieb S. Ling:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Hi</div><div><br></div><div>Your problem is associated with doing berry phase analysis with smearing. You can find more discussions on this topic in this post:</div><div><br></div><div><a href="https://groups.google.com/d/msg/cp2k/csKciEYJBk8/nUvUekPg-lwJ" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/msg/cp2k/csKciEYJBk8/nUvUekPg-lwJ';return true;" onclick="this.href='https://groups.google.com/d/msg/cp2k/csKciEYJBk8/nUvUekPg-lwJ';return true;">https://groups.google.com/d/<wbr>msg/cp2k/csKciEYJBk8/nUvUekPg-<wbr>lwJ</a></div><div><br></div><div>SL</div><div><br></div><div><br><br>On Wednesday, January 25, 2017 at 12:11:59 PM UTC, Phil G. wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear CP2K enthusiasts,<br><br>after the geometry optimization of the crystal slab and then again the geometry optimization of the system with water added and then run MD simulation (about 500 steps), I would like to obtain vibrational frequencies of the interfacial water with CP2K. The input code is written in the case of only calculate all water molecules (see the atom indices 241..390 in the INVOLVED_ATOMS section). The system consists of a crystal system LiNbO3 (2x2x2 unit cells with periodic boundary) and a water system (~50 water molecules). After some attempts to run vibrational analysis calculation I got error message, that the Berry phase moments were not implemented:<br><br> CP2K| UNIMPLEMENTED, Berry phase moments for non uniform MOs' occupation numbers not implemented<br> CP2K| Abnormal program termination, stopped by process number 48<br> CP2K| UNIMPLEMENTED, Berry phase moments for non uniform MOs' occupation numbers not implemented<br> CP2K| Abnormal program termination, stopped by process number 51<br> CP2K| UNIMPLEMENTED, Berry phase moments for non uniform MOs' occupation numbers not implemented<br> CP2K| Abnormal program termination, stopped by process number 65<br>[...]<br><br>I have not idea what the reason is and I am looking for a help or hints to solve the problem.<br>The input file is attached here.<br><br>Kind regards,<br><br>Phil<br></div></blockquote></div></div></blockquote></div>