<div dir="ltr">Dear Wei,<br><br>The first thing I see in your input is that you specify one colvar in SUBSYS, but in MOTION you call two colvars. This cannot work. As it is only one bond you specify as colvar, you should work with one colvar in the CONSTRAINT section. Then you should keep it "fixed" like this in section MOTION:<br><br>&CONSTRAINT<br>  &COLLECTIVE<br>    COLVAR 1<br>    INTERMOLECULAR<br>    TARGET 2.399706<br>    &RESTRAINT<br>      K 1.0<br>    &END RESTRAINT<br>  &END COLLECTIVE<br>&END CONSTRAINT<br><br>Please keep in mind that the distance you want it to be fixed (TARGET) is in bohrs, not angstroms. You need to experiment a little with the force constant for the restraint. Check it in your runs if it is sufficient to keep the bond length fixed within your desired accuracy. In principle, it should do the job. I have used this kind of constraint successfully for both, geometry optimizations and MD.<br><br>Best regards<br><br>Andreas<br><br><br>Am Donnerstag, 13. November 2014 12:31:24 UTC+1 schrieb Wei:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Pankaj, Ari and Marco,<br><br>I'd also like to fix a bond during a MD (DFT) run (not the absolute atomic positions). I tried to specify CP2K_INPUT / MOTION / CONSTRAINT / COLLECTIVE and CP2K_INPUT / FORCE_EVAL / SUBSYS / COLVAR by following an expample from regtest /Fist/regtest-12/water_3_dist_<wbr>2.inp.<br><br>However, the calculation always crashes with an error:<br><br> *<br> *** 12:16:30 ERRORL2 in topology_input:read_<wbr>constraints_section processor  ***<br> ***      0  err=-300  condition FAILED at line 561                         ***<br> *<br><br> ===== Routine Calling Stack =====<br><br>            3 topology_control<br>            2 quickstep_create_force_env<br>            1 CP2K<br> CP2K| condition FAILED at line 561<br> CP2K| Abnormal program termination, stopped by process number 0<br><br><br>Could you please give me some hint on this issue? I put my input file below. Thanks in advance!<br><br>Best Regards<br><br>Wei<br><br>input <br><br>&FORCE_EVAL<br>  METHOD QUICKSTEP<br>  &DFT<br>    BASIS_SET_FILE_NAME ./GST.qbs<br>    POTENTIAL_FILE_NAME ./GTH_POTENTIALS<br>    &MGRID<br>      CUTOFF 300<br>    &END MGRID<br>    &QS<br>      EPS_DEFAULT 1.0e-12<br>      EXTRAPOLATION ASPC<br>      EXTRAPOLATION_ORDER 0<br>      MAP_CONSISTENT<br>    &END QS<br>    &SCF<br>      EPS_SCF 1.0e-7<br>      MAX_SCF 40<br>      MAX_SCF_HIST 2<br>      &OT ON<br>        MINIMIZER DIIS<br>        PRECONDITIONER FULL_KINETIC<br>        ENERGY_GAP 0.002<br>        STEPSIZE 0.0875<br>      &END OT<br>      &OUTER_SCF ON<br>        MAX_SCF 800<br>        EPS_SCF 1.0e-7<br>      &END OUTER_SCF<br>      SCF_GUESS RANDOM<br>#     SCF_GUESS HISTORY_RESTART<br>    &END SCF<br>    &XC<br>      &XC_FUNCTIONAL PBE<br>      &END XC_FUNCTIONAL<br>    &END XC<br>  &END DFT<br>  &SUBSYS<br>    &CELL<br>      ABC   18.62463  18.62463  18.62463<br>    &END CELL<br>    &COORD<br> Ge 10.0 10.0 10.0  bond<br> Te 13.0 10.0 10.0  bond<br>    &END COORD<br>    &COLVAR<br>      &DISTANCE<br>        ATOMS 1 2<br>      &END DISTANCE<br>    &END COLVAR<br>    &KIND Ge<br>      BASIS_SET DZVP-MOLOPT-SR-GTH-q4<br>      POTENTIAL GTH-PBE-q4<br>    &END KIND<br>    &KIND Te<br>      BASIS_SET DZVP-MOLOPT-SR-GTH-q6<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br>  &END SUBSYS<br>&END FORCE_EVAL<br>&GLOBAL<br>  PREFERRED_FFT_LIBRARY FFTW<br>  PRINT_LEVEL LOW<br>  PROJECT GeTe.md<br>  RUN_TYPE MD<br>&END GLOBAL<br><br>&MOTION<br>  &MD<br>    ENSEMBLE LANGEVIN<br>    STEPS    20<br>    TEMPERATURE  1100<br>    TIMESTEP 2.0<br>    &LANGEVIN<br>      GAMMA 4.0E-3<br>      NOISY_GAMMA 3.5E-5<br>    &END LANGEVIN<br>  &END MD<br>  &CONSTRAINT<br>    &COLLECTIVE<br>      MOLNAME  bond<br>      COLVAR 1<br>    &END COLLECTIVE<br>    &COLLECTIVE<br>      MOLNAME  bond<br>      COLVAR 2<br>    &END COLLECTIVE<br>  &END CONSTRAINT<br>&END MOTION<div dir="ltr"><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br></div></blockquote></div></blockquote></div></div></blockquote></div>