<div dir="ltr">Dear Kamil,<div><br></div><div>this one worked:</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT piglet</div><div>  RUN_TYPE                   PINT</div><div>  PRINT_LEVEL                MEDIUM</div><div>&END GLOBAL</div><div>&MOTION</div><div>  &PINT</div><div>    P                        12</div><div>    NUM_STEPS                3</div><div>    DT                       0.25</div><div>    NRESPA                   2<span style="white-space:pre">           </span>#produces an intended warning.</div><div>    TEMP                     300.0</div><div>    TRANSFORMATION           NORMAL</div><div>    PROPAGATOR               RPMD</div><div>    HARM_INT                 EXACT</div><div>    PROC_PER_REPLICA         2</div><div>    &NORMALMODE</div><div>      MODEFACTOR             1.0</div><div>      Q_BEAD                -1.0</div><div>      Q_CENTROID            -1.0</div><div>    &END NORMALMODE</div><div>    &PIGLET</div><div>      SMATRIX_INIT          DIAGONAL</div><div>      MATRICES_FILE_NAME    ./PIGLET_P12_300K</div><div>    &END PIGLET</div><div>    &INIT</div><div>      RANDOMIZE_POS          F</div><div>      CENTROID_SPEED         F</div><div>      THERMOSTAT_SEED        12</div><div>    &END INIT</div><div>    &PRINT</div><div>      &CENTROID_GYR          OFF</div><div>      &END CENTROID_GYR</div><div>    &END PRINT</div><div>  &END PINT</div><div>&END MOTION</div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME    GTH_BASIS_SETS</div><div>    POTENTIAL_FILE_NAME    GTH_POTENTIALS </div><div><br></div><div>    &MGRID</div><div>      CUTOFF    100</div><div>      REL_CUTOFF 50</div><div>      ! SKIP_LOAD_BALANCE_DISTRIBUTED</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-10</div><div>      EPS_PGF_ORB 1.0E-10</div><div>    &END QS</div><div>    &SCF</div><div>      SCF_GUESS ATOMIC</div><div>      EPS_SCF 1.0E-5</div><div>      MAX_SCF 30</div><div><br></div><div>      &OUTER_SCF</div><div>        EPS_SCF  1.0E-7</div><div>        MAX_SCF  2</div><div>      &END</div><div>      &PRINT</div><div>        &RESTART ON</div><div>        &END</div><div>      &END</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL BLYP</div><div>      &END XC_FUNCTIONAL</div><div>    &END XC</div><div>  &END DFT</div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 5.0 5.0 5.0</div><div>    &END CELL</div><div>    &COORD</div><div>    O 0.01943016 0.956552982 0.7100297213</div><div>    H 0.15760040 0.330345154 1.477344990 </div><div>    H 0.90505028 1.198622227 0.3136827052 </div><div>    &END COORD</div><div>    &KIND H</div><div>      BASIS_SET  SZV-GTH</div><div>      POTENTIAL  GTH-BLYP-q1</div><div>     &END KIND</div><div>     &KIND O</div><div>      BASIS_SET  SZV-GTH</div><div>      POTENTIAL  GTH-BLYP-q6</div><div>     &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div>&EXT_RESTART</div><div>    RESTART_FILE_NAME piglet-1.restart</div><div>&END       </div><div><br></div><div>The file PIGLET_P12_300K can be generated at http://gle4md.org/index.html?page=matrix</div><div>One should select the proper temperature, number of beads and the output format (CP2K).</div><div><br></div><div>Of course, set in your own &FORCE_EVAL and &SUBYS sections.</div><div><br></div><div>Yours,</div><div><br></div><div>Vladimir</div><div><br></div><div><br></div><br>суббота, 27 января 2018 г., 14:05:13 UTC+1 пользователь Kamil написал:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">I am sorry, I would love to get into PIMD and PIGLET calculations with CP2K. Can I ask anybody to provide any examplary and workable PIGLET input? I would eally appreciate your kind help! Kamil</div></blockquote></div></div>