<div dir="ltr"><div>Hi All,</div>   <div>Recently, I want to run centroid molecular dynamics with CP2K. But there is no example. </div><div>The following input is got with the manual programming by me. But it cannot achieve the goal.</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT rVV10</div><div>  RUN_TYPE PINT</div><div>  PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div>  &PINT</div><div>    P                        16  #beads</div><div>    PROC_PER_REPLICA        12</div><div>    NUM_STEPS               10000</div><div>    TEMP                    5000.0</div><div>    DT                      0.1</div><div>    NRESPA                  8</div><div>    TRANSFORMATION          NORMAL</div><div>    PROPAGATOR              PIMD</div><div>    &NOSE</div><div>      NNOS 3</div><div>    &END NOSE</div><div>    &INIT       </div><div>      VELOCITY_SCALE        T</div><div>    &END INIT</div><div>  &END PINT</div><div>  &PRINT</div><div>    &VELOCITIES</div><div>    &END VELOCITIES</div><div>  &END PRINT</div><div>&END MOTION</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  STRESS_TENSOR ANALYTICAL </div><div>  &DFT</div><div>    MULTIPLICITY 1</div><div>    &MGRID</div><div>      CUTOFF 500</div><div>      NGRIDS 4</div><div>      REL_CUTOFF 50</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-12</div><div>      EPS_GVG 1.0E-6</div><div>      EPS_PGF_ORB 1.0E-6</div><div>    &END QS</div><div>    &SCF</div><div>      MAX_SCF 200</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-5</div><div>      ADDED_MOS 64</div><div>      CHOLESKY INVERSE</div><div>      &SMEAR ON</div><div>        METHOD FERMI_DIRAC</div><div>        ELECTRONIC_TEMPERATURE 5000</div><div>      &END SMEAR</div><div>      &DIAGONALIZATION</div><div>        ALGORITHM STANDARD</div><div>      &END DIAGONALIZATION</div><div>      &OT off</div><div>        MINIMIZER CG</div><div>        PRECONDITIONER FULL_ALL</div><div>      &END OT</div><div>      &MIXING</div><div>        METHOD BROYDEN_MIXING</div><div>        ALPHA 0.4</div><div>        BETA 0.5</div><div>        NBROYDEN 8</div><div>      &END MIXING</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL</div><div>        &LIBXC</div><div>          FUNCTIONAL XC_GGA_X_RPW86 XC_GGA_C_PBE</div><div>        &END LIBXC</div><div>      &END XC_FUNCTIONAL</div><div>      &VDW_POTENTIAL</div><div>       POTENTIAL_TYPE NON_LOCAL</div><div>        &NON_LOCAL</div><div>          TYPE RVV10</div><div>          CUTOFF   650</div><div>          PARAMETERS  6.3  0.0093</div><div>          KERNEL_FILE_NAME rVV10_kernel_table.dat</div><div>        &END NON_LOCAL</div><div>      &END VDW_POTENTIAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC  15.6 15.6 15.6</div><div>    &END CELL</div><div>    &COORD</div><div>    &END COORD</div><div>    &VELOCITY</div><div>    &END VELOCITY</div><div>    &KIND H                    </div><div>      BASIS_SET DZVP-GTH-PBE      </div><div>      POTENTIAL GTH-PBE-q1</div><div>    &END                                                                  </div><div>  &END SUBSYS                  </div><div>  &PRINT </div><div>    &STRESS_TENSOR</div><div>    &END STRESS_TENSOR</div><div>  &END PRINT</div><div>&END FORCE_EVAL                 </div><div><br></div></div></div>