Hello Dr.Hutter,<div><br /></div><div>Thanks for your response. In this case, would the rest of the frequencies (3N-6) still be reliable? or they are also impacted by the fixed atoms?</div><div><br /></div><div>Thanks!</div><div><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, May 8, 2023 at 2:39:22 AM UTC-5 Jürg Hutter wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi <br>If you fix atomic coordinates, your system looses its free translations and rotations. <br>For such a system you expect 3N non-zero frequencies. <br>regards <br>JH <br> <br>________________________________________ <br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Hana <<a href data-email-masked rel="nofollow">holi...@illinois.edu</a>> <br>Sent: Friday, May 5, 2023 5:19 PM <br>To: cp2k <br>Subject: Re: [CP2K:18774] VIBRATIONAL ANALYSIS <br> <br>Hello Lucas and Krack, <br> <br>Thanks for your comments. I have attached my input file (input.txt) for your notice. My system consists of a single ion pair (<a href="http://sys.xyz" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=http://sys.xyz&source=gmail&ust=1683647703733000&usg=AOvVaw3oS9ZVGsaZvoySFZjqMn3R">sys.xyz</a>) with no periodic boundary condition and I use the GAPW method with a 400Ry cutoff (Upon convergence test, this cutoff gave me 1e-4 convergence threshold for the total energy). For the normal mode analysis I fix the atoms of one ion and just study the vibrations of the other one. Please let me know what you think. <br> <br>Thanks! <br> <br> <br>On Friday, May 5, 2023 at 2:31:52 AM UTC-5 Lucas Lodeiro wrote: <br>I put my 2 cents in the possibility that you are running your calculation using PBC, so the normal modes are not projected out to extract the translational and rotational modes. Also you probably use a "low cutoff" so the noise (energy/gradient waves due to plane waves), when the numerical derivatives are calculated through adjacent points, is important... then the frequency of your translational and rotational modes are not close to zero as expected. <br> <br>Regards - Lucas Lodeiro <br> <br>El jue, 4 may 2023 a las 12:36, Hana (<<a href data-email-masked rel="nofollow">holi...@illinois.edu</a>>) escribió: <br>Hello, <br> <br>I am trying to obtain the normal modes of my gaseous system which contains a molecule (N atoms). First I optimize the geometry of my system with a MAX_FORCE threshold of 0.0005. Second, in a further calculation, I use the following setting: <br> <br>&VIBRATIONAL_ANALYSIS <br> INTENSITIES <br> NPROC_REP 256 <br> DX 0.001 <br> &PRINT <br> &PROGRAM_RUN_INFO ON <br> &END <br> &END <br>&END <br> <br>I get 3N frequency modes all of which are nonzero. However, I expect the analysis to exclude the rotation and translation and return zeros for the first six modes. Is there any explanation for this? Here is the output (cm^-1) for the first six modes: 26.89, 30.12, 36.90, 55.93, 68.51, 90.65 <br> <br>Thanks! <br> <br>-- <br>You received this message because you are subscribed to the Google Groups "cp2k" group. <br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>. <br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%2540googlegroups.com&source=gmail&ust=1683647703733000&usg=AOvVaw17MpiIT76w8-pUi7VPldiW">https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1683647703733000&usg=AOvVaw0DATGAT9155tpfv27w9i98">https://groups.google.com/d/msgid/cp2k/0ab349dd-fe8f-4f62-bba5-c4fc2cc68a06n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. <br> <br>-- <br>You received this message because you are subscribed to the Google Groups "cp2k" group. <br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>. <br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%2540googlegroups.com&source=gmail&ust=1683647703733000&usg=AOvVaw3fhZzUVuPD7EQa89GLan4I">https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1683647703733000&usg=AOvVaw16VjVsBSQzug_d5L_6bgiR">https://groups.google.com/d/msgid/cp2k/f772710a-d0e9-4247-9f9f-58e245bf8535n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>. <br></blockquote></div> <p></p> -- <br /> You received this message because you are subscribed to the Google Groups "cp2k" group.<br /> To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br /> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b9552901-a336-4e82-bbaa-5c3ecf75c851n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/b9552901-a336-4e82-bbaa-5c3ecf75c851n%40googlegroups.com</a>.<br />