<div dir="ltr"><div><div>Dear colleagues, </div><div><br></div><div>I'm trying to obtain thermochemistry data, like enthalpy, entropy and Gibbs energies correction, same analysis as gaussian does on the Freq calculation.</div><div>I tried with VIBRATIONAL_ANALYSIS using THERMOCHEMISTRY, however I'm not able to obtain the entropic contribution and the final Gibbs energy.</div><div><br></div><div>Is some information missing on the input?</div><div><br></div><div>Please find attached below the input:</div><br><div>&GLOBAL</div><div>  PROJECT $SLURM_JOB_NAME</div><div>  RUN_TYPE VIBRATIONAL_ANALYSIS</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &DFT</div><div>    &PRINT</div><div>     &MOMENTS</div><div>     &END</div><div>    &END</div><div>  WFN_RESTART_FILE_NAME $SLURM_JOB_NAME-RESTART.wfn</div><div>    BASIS_SET_FILE_NAME $HOME/scr/cp2k/data/BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME $HOME/scr/cp2k/data/GTH_POTENTIALS  </div><div>    CHARGE 0</div><div>    MULTIPLICITY 1</div><div>    &MGRID</div><div>       CUTOFF [Ry] 500</div><div>    &END</div><div>    &QS</div><div>       METHOD GPW </div><div>       EPS_DEFAULT 1.0E-10 </div><div>    &END</div><div>    &POISSON</div><div>       PERIODIC NONE</div><div>       POISSON_SOLVER MULTIPOLE</div><div>    &END</div><div>    &SCF                              </div><div>    SCF_GUESS RESTART</div><div>      MAX_SCF 30</div><div>      EPS_SCF 1.0E-5</div><div>      !!CHOLESKY OFF</div><div>       &OT</div><div>        PRECONDITIONER FULL_ALL</div><div>        ENERGY_GAP 2E-03</div><div>        MINIMIZER DIIS</div><div>       &END OT</div><div>      &OUTER_SCF </div><div>        MAX_SCF 30</div><div>        EPS_SCF 1.0E-5 ! must match the above</div><div>      &END</div><div>    &END SCF</div><div>    &XC</div><div>     &XC_FUNCTIONAL </div><div>         &PBE</div><div>!!         PARAMETRIZATION PBESOL</div><div>         &END</div><div>     &END XC_FUNCTIONAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL </div><div>      ABC [angstrom] 15. 15. 15.</div><div>      ALPHA_BETA_GAMMA 90. 90. 90.</div><div>      SYMMETRY NONE</div><div>      PERIODIC NONE</div><div>    &END CELL</div><div>  &TOPOLOGY</div><div>    COORD_FILE_NAME $SLURM_JOB_NAME.xyz</div><div>    COORD_FILE_FORMAT XYZ</div><div>  &END</div><div>   </div><div>    &KIND   H</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND   C</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND   N</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND   O</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND   P</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND  Co</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND   I</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>    &KIND  Ni</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE</div><div>    &END KIND</div><div><br></div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&VIBRATIONAL_ANALYSIS</div><div>  FULLY_PERIODIC</div><div>  INTENSITIES</div><div>  THERMOCHEMISTRY</div><div>  TC_TEMPERATURE 298.15</div><div>  TC_PRESSURE <font color="#000000" face="monospace"><span style="text-transform: uppercase;">1.01325000E+005</span></font></div><div>  NPROC_REP 12</div><div>!!   &MODE_SELECTIVE</div><div>!!     ATOMS 1..46</div><div>!!    &INVOLVED_ATOMS</div><div>!!      INVOLVED_ATOMS 1..46</div><div>!!    &END</div><div>!!   INITIAL_GUESS ATOMIC</div><div>!!  EPS_NORM 1.0E-04</div><div>!!   EPS_MAX_VAL 1.0E-05</div><div>!!  &END</div><div>&END    </div><div><br></div></div><div>And the output information:</div><div><br></div><div><div> VIB|                        NORMAL MODES - THERMOCHEMICAL DATA</div><div> VIB|</div><div> VIB|              Symmetry number:                                                 1</div><div> VIB|              Temperature [K]:                                            298.15</div><div> VIB|              Pressure [Pa]:                                           101325.00</div><div><br></div><div><br></div><div> VIB|              Electronic energy (U) [kJ/mol]:                  -1023152.20665924</div><div> VIB|              Zero-point correction [kJ/mol]:                       956.45314459</div><div> VIB|              Entropy [kJ/(mol K)]:                                             Nan</div><div> VIB|              Enthalpy correction (H-U) [kJ/mol]:                  1009.07755543</div><div> VIB|              Gibbs energy correction [kJ/mol]:                         Nan</div><div> VIB|              Heat capacity [kJ/(mol*K)]:                             0.38550655</div></div><div><br></div><div>Thank you.</div><div>Gerard.</div></div>