<div dir="ltr">Hi Leopold,<br><br>I looked at it again and I think now I was mistaken. You are right, the formulas for frequency and force constant are not consistent.<br><br>Matthias<br><br><br>On Sunday, May 25, 2014 11:33:06 PM UTC+2, Leopold Talirz wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Matthias,<div><br></div><div>thank you very much for your quick reply.</div><div>I am really not well familiar with normal mode analysis (nor its implementation), so please tell me if I am talking nonsense, but I thought the eigenvalues of the Hessian were proportional to konst instead of SQRT(konst)?</div><div><br></div><div>At least with the naive definition of a force constant that I am familiar with, i.e. you have a potential like U=1/2 * konst * x^2</div><div><br></div><div>Or, in other words, could you please direct me to a reference/book that I should read to understand how the force constants are defined in CP2K?</div><div><br></div><div>Thank you very much!</div><div><br></div><div>Best,</div><div>Leopold</div><div><br></div><div>On Sunday, May 25, 2014 12:38:17 PM UTC+2, Matthias Krack wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Leopold,<br><br>I think the code is correct, since the eigenvalues epsilon of the Hessian returned in H_eigval2 after its diagonalisation are equivalent to SQRT(konst/rmass). The frequencies, however, are proportional to SQRT(epsilon), i.e. SQRT(H_eigval2).<br><br>Matthias</div></blockquote></div></div></blockquote></div>