<div>Hi <span>Jürg</span></div><div><span><br /></span></div><div><span><br /></span></div><div><span>Thanks for your reply, I understand there will be limitations to the invariance, i.e. not perfectly invariant, but if you glance at the forces they do not really seem invariant at all, as in lots of forces are in opposite directions etc, do you think I could have set up the input file incorrectly or something as I could not see any mistakes? <br /></span></div><div><br /></div><div>Thanks </div><div><br /></div><div>Harry </div><div><span> </span></div><div><span><br /></span></div><div><span><br /></span></div><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, June 17, 2025 at 8:45:15 AM UTC+1 Jürg Hutter wrote:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>the periodic Hamiltonian with plane wave basis sets (as partly used in CP2K or in PW codes) is not
<br>rotationally and translationally (only up to quantas of grid spacings) invariant.
<br>
<br>The variation in energy/forces depends on several factors, but exact invariance will only be achieved
<br>for infinite box size and infinite PW cutoff.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of Harry Richardson <<a href data-email-masked rel="nofollow">harrysgr...@gmail.com</a>>
<br>Sent: Monday, June 16, 2025 1:07 PM
<br>To: cp2k
<br>Subject: [CP2K:21545] Forces not equivariant upon molecule rotation
<br>
<br>
<br>Hi I am running some energy+force calculation of a single benzene molecule using WB97X-D (benzene.inp).
<br>
<br>In order to test the equivariance of the force calculations I have rotated the molecule by 180 degrees around the vector [0.5,0.5,0.5] via the file rotation.ipynb.
<br>
<br>I then calculate the forces for both (forces_0, forces_180) and use the file equi_checker.ipynb to rotate the forces back to the original space and compare them.
<br>
<br>When I do this the forces do not appear equivariant at all, in fact the magnitude of the difference in forces is of the order of the forces themselves i.e. converting force to ev/A:
<br>
<br>sum absolute difference 270 deg rotation: 1.651721858651293 eV/A sum absolute difference 180 deg rotation: 1.5320346141088304 eV/A sum absolute difference 90 deg rotation: 0.6562673988482665 eV/A
<br>
<br>and this is when the total absolute forces themselves are ~1.55ev/A
<br>
<br>Does anyone have any ideas or can see any issues with my input / outputs forces? I have tried tightening convergence and trying different potentials, none seemed to make much difference?
<br>
<br>Thank you for any help anyone can provide
<br>
<br>--
<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 data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%2540googlegroups.com&source=gmail&ust=1750323211073000&usg=AOvVaw0ekUxXc6GRKlLNoMIdQhAu">https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en&q=https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1750323211073000&usg=AOvVaw3F-JhWZ7fNszhh_4vlzzN5">https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></div>
<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+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/5db481e0-16ed-486e-ac3d-d814b0d208e5n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/5db481e0-16ed-486e-ac3d-d814b0d208e5n%40googlegroups.com</a>.<br />