<div dir="ltr">Dear Matt,<br><br>This is really a serious bug. Now the reason is clear to me for translating the system in -ve direction of the field which i was trying to figure out fundamentally from last more than 2 months. After recompiling the code with this change, My system (Single Hydrogen atom in constant EF and single water molecule in Gaussian Pulse both for non-periodic box) is translationally invariant  in presence of EF. This is  really helpful. <br><br>Thank you so much <br><br>Pankaj<br><br>On Friday, 8 November 2013 19:14:32 UTC+1, Matt W  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Hi,<br></div><div><br></div><div>You make a very good point! Now for your argument to hold your system needs to be translationally invariant. There are two reasons (at least) why it is not.<br><br>(i) the EFIELD is non periodic (i.e. a saw tooth), so there is a huge field at the cell edge. You need to centre your molecule in the box.<br>(ii) your CUTOFF is not truly sufficient if you really want to see the molecule stay centred - there's a post from Flo discussing this concerning vibrational analysis from a month or so ago, but this is a minor point.</div><div><br></div><div>BUT, I think you have also identified a bug. According to me, there is an inconsistency between the way forces due to the field act on atoms cores and electrons.</div><div><br></div><div>changing line 233 (or thereabouts) of efield_utils.F to</div><div><br></div><div>force(ikind)%efield(:,iatom)=<wbr>force(ikind)%efield(:,iatom)-<wbr>field*zeff<br></div><div><br></div><div>i.e. changing the sign of the forces on the cores, seems to give good conservation of force in my quick tests (plus comments above). You'll have to recompile, of course.</div><div><br></div><div>Let us know if this seems to work for you (it is possible that the force on the electrons rather than the core should be inverted, but at the moment they are both being pushed in the same direction).</div><div><br></div><div>Matt</div><div><br></div>On Friday, 8 November 2013 09:10:51 UTC, Pankaj Mishra  wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<br><br>As per my fundamental knowledge, polar molecule can rotate in electric field, not translate. Being more precise, here i am talking about 1 molecule in non-periodic box. I didn't find any force to move in the field. <br><br>Even non-polar molecule and single atom (H and Ne) are also translating in EF in CP2K. In the case of water cluster, the center of mass of whole system is translating. I have not found such behavior in other package like GAMESS .<br><br>This is a very serious issue for me. Please try to suggest something if i am missing anything in my input.<br><br>Thank you<br><br>Pankaj <br><br>On Thursday, 7 November 2013 20:23:58 UTC+1, Matt W  wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hi,<br></div><div><br></div><div>just wondering: why would you not expect a polar molecule to move in an electric field?</div><div><br></div><div>Matt</div></div></blockquote></div></blockquote></div></blockquote></div>