<div dir="ltr"><div>Dear CP2K members.</div><div>I am simulating a liquid phase reaction, with PBE and BLYP functionals (and dispersion correction), and the reaction free energy do not agree with the experimental value. I ask colleagues about this, and they said me that the functionals maybe are not the best for the reaction, they suggested HCTH functional. I ran simulations with various types of HCTH (like as 120 and 407). I did not have problems with HCTH120, but with HCTH407 I have problems, because the program does not recognize the functional like as I put in the forces.inc. I used 5.1 CP2K version.</div><div><br></div><div>In the HCTH120 calculation, the XC_FUNCTIONAL section is:</div><div><br></div><div><span style="white-space:pre">            </span>&XC_FUNCTIONAL  HCTH120</div><div>                        &END XC_FUNCTIONAL</div><div>                        &VDW_POTENTIAL</div><div>                                POTENTIAL_TYPE   PAIR_POTENTIAL </div><div>                                &PAIR_POTENTIAL</div><div>                                        TYPE                    DFTD3</div><div>                                        CALCULATE_C9_TERM       TRUE</div><div>                                        PARAMETER_FILE_NAME     ${dftd3}</div><div>                                        REFERENCE_FUNCTIONAL    HCTH120</div><div>                                &END PAIR_POTENTIAL</div><div>                        &END VDW_POTENTIAL</div><div><br></div><div>If I change the number 120 for 407, the program throws an error for the functional (invalid value for enumeration:HCTH). I try to follow tha manual (https://manual.cp2k.org/cp2k-2_3-branch/CP2K_INPUT/ATOM/METHOD/XC/XC_FUNCTIONAL/HCTH.html) to use it, but the program throws an error too. The input is:</div><div><br></div><div>                        &XC_FUNCTIONAL  HCTH</div><div>                                &HCTH </div><div>                                        SECTION_PARAMETERS TRUE</div><div>                                        PARAMETER_SET 407</div><div>                                &END HCTH</div><div>                        &END XC_FUNCTIONAL</div><div>                        &VDW_POTENTIAL</div><div>                                POTENTIAL_TYPE   PAIR_POTENTIAL </div><div>                                &PAIR_POTENTIAL</div><div>                                        TYPE                    DFTD3</div><div>                                        CALCULATE_C9_TERM       TRUE</div><div>                                        PARAMETER_FILE_NAME     ${dftd3}</div><div>                                        REFERENCE_FUNCTIONAL    HCTH</div><div>                                &END PAIR_POTENTIAL</div><div><br></div><div>In despair, I tried strange things like as: Change "&XC_FUNCTIONAL  HCTH" to "&XC_FUNCTIONAL  HCTH120" ; Change "REFERENCE_FUNCTIONAL    HCTH" to "REFERENCE_FUNCTIONAL    HCTH407"... And the mixtures of these changes, but the error remains (found an unknown keyword SECTION_PARAMETERS in section HCTH).</div><div><br></div><div>Could anyone give an advice on how to write the input force.inc to use the HCTH407 type?</div><div><br></div><div>Thanks - Lucas Lodeiro</div></div>