<div>Dear Naina,</div><div>Considering this information (https://gaussian.com/vib/), both codes should agree. If you say that the frequencies (or wavenumbers) are matching (with what error??), the force constants should match as well because they are directly related to each other (see the document). Considering your input, I have the following suggestions</div><div>- The cutoff could be too low. Secondary, double-check EPS_DEFAULT and EPS_SCF.</div><div>- If you actually have a gas phase system: Change the periodicity to NONE (in the cell definition and in the POISSON section) and remove the FULLY_PERIODIC keyword in the VIBRATIONAL_ANALYSIS section. Add the TOPOLOGY section with the CENTER_COORDINATES section (see https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/TOPOLOGY/CENTER_COORDINATES.html ).</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 17:17: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;"><div>The frequencies and the reduced masses are very similar in both Gaussian and CP2K.(probably because both of them report in the same units).</div><div><br></div><div>I've included for you the CP2k input for frequency cal.</div><div><br></div>&GLOBAL<br>  PROJECT H2O<br>  RUN_TYPE VIBRATIONAL_ANALYSIS<br>  PRINT_LEVEL MEDIUM<br>&END GLOBAL<br>&FORCE_EVAL<br>  METHOD QS<br>  &DFT<br>    BASIS_SET_FILE_NAME BASIS_MOLOPT<br>    POTENTIAL_FILE_NAME GTH_POTENTIALS<br>    CHARGE 0<br>    MULTIPLICITY 1<br>&PRINT<br>   &MOMENTS<br>   &END<br> &END    <br>&MGRID<br>      CUTOFF 300<br>      NGRIDS 5<br>      REL_CUTOFF 70<br>    &END MGRID<br>    &QS<br>      EPS_DEFAULT 1.0E-12<br>      WF_INTERPOLATION ASPC<br>    &END QS<br>   &SCF<br>      SCF_GUESS ATOMIC<br> EPS_SCF 1.0E-7<br> MAX_SCF 1000<br> CHOLESKY INVERSE<br> ADDED_MOS 50<br>      &DIAGONALIZATION<br>       ALGORITHM STANDARD<br>      &END DIAGONALIZATION<br>      &MIXING<br>       METHOD BROYDEN_MIXING<br>      ALPHA 0.4<br>      NBROYDEN 8<br>      &END MIXING<br>    &END SCF<br>    &XC<br>      &XC_FUNCTIONAL<br>        &PBE<br>        &END PBE<br>      &END XC_FUNCTIONAL<br>      &vdW_POTENTIAL<br>        DISPERSION_FUNCTIONAL PAIR_POTENTIAL<br>         &PAIR_POTENTIAL<br>            PARAMETER_FILE_NAME dftd3.dat<br>                TYPE DFTD3<br>                REFERENCE_FUNCTIONAL PBE<br>                R_CUTOFF 15.0<br>        &END PAIR_POTENTIAL<br>     &END vdW_POTENTIAL<br>    &END XC<br>    &POISSON<br>      PERIODIC xyz<br>POISSON_SOLVER WAVELET<br>&END POISSON<br>  &END DFT<br>  &SUBSYS<br>    &CELL<br>ABC 20.0000 20.0000 20.0000<br>ALPHA_BETA_GAMMA 90.0 90.0 90.0<br>      PERIODIC xyz<br>    &END CELL    <br>&COORD<br>O    -0.762288846626      0.000000000001     -1.462238785752<br>H    -0.762288846634      0.759337000000     -0.866195785751<br>H    -0.762288846634     -0.759337000000     -0.866195785753<br>&END COORD<br>    &KIND Mg<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      POTENTIAL GTH-PBE-q10<br>    &END KIND<br>    &KIND O<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br>&KIND Y<br>ELEMENT Y<br>BASIS_SET DZVP-MOLOPT-SR-GTH<br>POTENTIAL GTH-PBE-q11<br>&END KIND<br>&KIND C<br>ELEMENT C<br>BASIS_SET DZVP-MOLOPT-GTH<br>POTENTIAL GTH-PBE-q4<br>&END KIND<br>&KIND H<br>ELEMENT H<br>BASIS_SET DZVP-MOLOPT-GTH<br>POTENTIAL GTH-PBE-q1<br>&END KIND   <br>&END SUBSYS<br>&END FORCE_EVAL<br>&VIBRATIONAL_ANALYSIS<br>FULLY_PERIODIC  <br>INTENSITIES<br>        NPROC_REP 4<br>&END VIBRATIONAL_ANALYSIS<br>&MOTION<br>  &CONSTRAINT<br>    &FIXED_ATOMS<br>    &END FIXED_ATOMS<br>  &END CONSTRAINT<br>&GEO_OPT<br>  OPTIMIZER LBFGS<br>TYPE MINIMIZATION<br>MAX_DR 4.5E-4<br>MAX_ITER 500<br>&END GEO_OPT<br>  &END MOTION<br><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Saturday, August 23, 2025 at 7:45:05 PM UTC+5:30 Naina Sethi 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"><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></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/cf8ea646-d957-4fb7-bc54-2b946cd46654n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/cf8ea646-d957-4fb7-bc54-2b946cd46654n%40googlegroups.com</a>.<br />