<div dir="ltr"><div dir="ltr"><div dir="ltr">    At first , i appreciate your help very much.<br>i try to use constraint of my atom positions.  fix_atom -->list 1 ..112 (112 atoms)<br>but it turn out a error  ,see the file "tdtest.out". what 's wrong ?  in cp2k 2.2 and 2.3 <br>it seems that we can't fix all atoms.<br><br>my input file is here.<br><br><br><br>by the way, why there is no  data of mullken affter the run.<br><br>   &MULLIKEN SILENT<br>        FILENAME ./mulliken.out<br>      &END MULLIKEN<br><br>thanks!<br>best regards.<br>xianjin<br><br><br><br><br><br>在 2013年12月18日星期三UTC+8下午9时48分28秒，<wbr>Samuel Andermatt写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr">You can constrain your system using parameters found in the section <a href="http://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/CONSTRAINT.html" target="_blank">http://manual.cp2k.org/trunk/<wbr>CP2K_INPUT/MOTION/CONSTRAINT.<wbr>html</a> .<br><br>For your case you could do an MD simulation with the positions of the atoms of your substance fixed (<a href="http://manual.cp2k.org/trunk/CP2K_INPUT/MOTION/CONSTRAINT/FIXED_ATOMS.html" target="_blank">http://manual.cp2k.org/trunk/<wbr>CP2K_INPUT/MOTION/CONSTRAINT/<wbr>FIXED_ATOMS.html</a>) . I assume your substance is in a solvent.<br>After your system equilibrated (temperatures no longer fluctuate) you could use the coordinates as input for a second simulation, where you remove the constraint of your atom positions.<br><br></div></div></blockquote></div></div></div>