<div>Dear Naina,</div><div>Are you sure that all relevant parameters (especially CUTOFF, EPS_DEFAULT, EPS_SCF) are well converged? Beware that a vibrational analysis requires much tighter thresholds than force calculations. Can you provide your input file? How do the frequencies compare between both codes?</div><div>Best,</div><div>Frederick</div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">Naina Sethi schrieb am Samstag, 23. August 2025 um 16:15:05 UTC+2:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><strong>Dear all,</strong><br>
I am working on calculating the vibrational frequencies and force constants of a water molecule. In CP2K, the force constants are reported in atomic units (Hartree/Bohr²), while Gaussian lists them in mDyne/Å. The conversion factor between Hartree/Bohr² and mDyne/Å is about 15.57. However, even after applying this factor, the results from CP2K do not match the force constant values reported by Gaussian. Could anyone clarify why this discrepancy occurs?  <br><div><br></div><div>       Frc consts  --           1.7016                   8.9360               9.7649        (gaussian)</div><div>

VIB|Frc consts (a.u.)    0.017668               0.567398         0.663392     (CP2K)</div><div><br></div><div>Thanks and Regards</div><div>Naina</div></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 visit <a href="https://groups.google.com/d/msgid/cp2k/3ad31466-622b-4328-aae5-6ba3484d25b1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/3ad31466-622b-4328-aae5-6ba3484d25b1n%40googlegroups.com</a>.<br />