<div dir="ltr">Are Na ions neutral? If not, you have to define the CHARGE under &DFT. </div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Aug 10, 2015 at 1:13 PM, jia <span dir="ltr"><<a href="mailto:huangji...@gmail.com" target="_blank">huangji...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><p style="text-align:center;clear:both"><a href="https://lh3.googleusercontent.com/-aDOG1a8wmP4/VcjXYKHMLHI/AAAAAAAAAMA/UeDXIFIt5h4/s1600/Screenshot%2Bfrom%2B2015-08-10.png" style="margin-left:1em;margin-right:1em" target="_blank"><img src="https://lh3.googleusercontent.com/-aDOG1a8wmP4/VcjXYKHMLHI/AAAAAAAAAMA/UeDXIFIt5h4/s320/Screenshot%2Bfrom%2B2015-08-10.png" border="0" width="320" height="244"></a></p><p style="text-align:center;clear:both">figure 1</p><p style="text-align:center;clear:both"><br></p><div><br></div><div><br></div><div>Hello,</div><div>The figure 1 is the result of DFTMD of a model system (2H2O - 2Na ) using CP2K 2.6. The trajectory of all atoms is shown here.</div><div>As you can see from the trajectory, the water molecules is moving during the 3 ps time period, but the Na atoms not. It seems that they are fixed. I have tried other systems (2H2O-2Li, 2H2O-2Cl) with similar input , no atoms is fixed there. (see figure 2)</div><div><br></div><div><br></div><p style="text-align:center;clear:both"><a href="https://lh3.googleusercontent.com/-T8CLO5k9pQ4/VcjbSO60CYI/AAAAAAAAAMM/BcEbcnVgK0E/s1600/Screenshot%2Bfrom%2B2015-08-10-fig2.png" style="margin-left:1em;margin-right:1em" target="_blank"><img src="https://lh3.googleusercontent.com/-T8CLO5k9pQ4/VcjbSO60CYI/AAAAAAAAAMM/BcEbcnVgK0E/s320/Screenshot%2Bfrom%2B2015-08-10-fig2.png" border="0" width="320" height="253"></a></p><p style="text-align:center;clear:both">figure 2</p><div><br></div><div><br></div><div>Could someone please suggest me where is wrong in the simulation? The input file is as follows. Thanks!</div><div><br></div><div>Best regards,</div><div>jia</div><div><br></div><div><br></div><div>input:</div><div><br></div><div>&FORCE_EVAL</div><div> METHOD QS</div><div> &DFT</div><div> BASIS_SET_FILE_NAME ./BASIS_MOLOPT</div><div> POTENTIAL_FILE_NAME ./GTH_POTENTIALS</div><div> &MGRID</div><div> CUTOFF 300</div><div> &END MGRID</div><div> &QS</div><div> EPS_DEFAULT 1.0E-12</div><div> WF_INTERPOLATION PS</div><div> EXTRAPOLATION_ORDER 3</div><div> &END QS</div><div> &SCF</div><div># SCF_GUESS ATOMIC</div><div> &OT ON</div><div># MINIMIZER DIIS</div><div> &END OT</div><div> SCF_GUESS RESTART</div><div> EPS_SCF 5.0E-7</div><div> MAX_SCF 300</div><div><br></div><div> &END SCF</div><div> &XC</div><div> &XC_FUNCTIONAL BLYP</div><div> &END XC_FUNCTIONAL</div><div> &vdW_POTENTIAL</div><div> DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div> &PAIR_POTENTIAL</div><div> TYPE DFTD3</div><div> CALCULATE_C9_TERM .TRUE.</div><div> REFERENCE_C9_TERM .TRUE.</div><div> LONG_RANGE_CORRECTION .TRUE.</div><div> PARAMETER_FILE_NAME ./dftd3.dat</div><div> REFERENCE_FUNCTIONAL PBE</div><div> R_CUTOFF 8.</div><div> EPS_CN 0.01</div><div> &END PAIR_POTENTIAL</div><div> &END vdW_POTENTIAL</div><div> &END XC</div><div> &END DFT</div><div> &SUBSYS</div><div> &CELL</div><div> ABC 31.00 15.6404 15.6404</div><div> &END CELL</div><div> &COORD</div><div> Na 1.2962887499096894E+01 1.2295273153310035E+01 1.4208821655264259E+01</div><div> Na 1.2291144356361610E+01 4.3563867602355567E-02 7.9483474865053152E+00</div><div> O 2.2241968510425036E+00 4.5518135434738101E+00 6.5536411797903718E+00</div><div> H 2.4867064046620353E+00 4.7403384589029143E+00 7.4756836487242460E+00</div><div> H 1.8258611320965084E+00 3.6350263378980761E+00 6.4892048508476918E+00</div><div> O 3.7704893907956558E+00 9.5503054982719924E+00 1.2806596984910234E+01</div><div> H 4.0764009642940087E+00 1.0180496826699695E+01 1.2115328686085089E+01</div><div> H 2.9239079442135907E+00 9.2048384808321888E+00 1.2405445167394735E+01</div><div> &END COORD</div><div> &KIND H</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH</div><div> POTENTIAL GTH-BLYP-q1</div><div> &END KIND</div><div><div> &KIND O</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH</div><div> POTENTIAL GTH-BLYP-q6</div><div> &END KIND</div><div> &KIND I</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH</div><div> POTENTIAL GTH-BLYP-q7</div><div> &END KIND</div><div> &KIND Li</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH</div><div> POTENTIAL GTH-BLYP-q3</div><div> &END KIND</div><div> &KIND Na</div><div> BASIS_SET DZVP-MOLOPT-SR-GTH</div><div> POTENTIAL GTH-BLYP-q9</div><div> &END KIND</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&GLOBAL</div><div> PROJECT 2w_2Na</div><div> RUN_TYPE MD</div><div> PRINT_LEVEL HIGH</div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div> &GEO_OPT</div><div> TYPE minimization</div><div> OPTIMIZER BFGS</div><div> MAX_ITER 20</div><div> &END GEO_OPT</div><div><br></div><div> &MD</div><div> &THERMOSTAT</div><div> &NOSE</div><div> LENGTH 3</div><div> YOSHIDA 3</div><div> TIMECON 1000.</div><div> MTS 2</div><div> &END NOSE</div><div> &END THERMOSTAT</div><div> ENSEMBLE NVT</div><div> STEPS 1000000</div><div> TIMESTEP 0.5</div><div> TEMPERATURE 330.0</div><div> TEMP_TOL 80</div><div> &END MD</div></div><div><br></div><div><div> &PRINT</div><div> &TRAJECTORY</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END TRAJECTORY</div><div> &VELOCITIES ON</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END VELOCITIES</div><div> &FORCES ON</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END FORCES</div><div><br></div><div> &RESTART_HISTORY</div><div> &EACH</div><div> MD 1000</div><div> &END EACH</div><div> &END RESTART_HISTORY</div><div> &RESTART</div><div> BACKUP_COPIES 1</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END RESTART</div><div> &END PRINT</div><div>&END MOTION</div></div><span class="HOEnZb"><font color="#888888"><div><br></div></font></span></div><span class="HOEnZb"><font color="#888888">
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</font></span></blockquote></div><br></div>