<div><br /></div><div>Hi I am running some energy+force calculation of a single benzene molecule using WB97X-D (benzene.inp). </div><div><br /></div><div>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. </div><div><br /></div><div>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. </div><div><br /></div><div>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 /><span style="color: rgb(204, 204, 204); font-family: Menlo, Monaco, "Courier New", monospace; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; word-spacing: 0px; white-space: pre; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial;">sum absolute difference 270 deg rotation: 1.651721858651293 eV/A
</span><span style="color: rgb(204, 204, 204); font-family: Menlo, Monaco, "Courier New", monospace; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; word-spacing: 0px; white-space: pre; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial;">sum absolute difference 180 deg rotation: 1.5320346141088304 eV/A
</span><span style="color: rgb(204, 204, 204); font-family: Menlo, Monaco, "Courier New", monospace; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; word-spacing: 0px; white-space: pre; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial;">sum absolute difference 90 deg rotation: 0.6562673988482665 eV/A</span></div><div><span style="color: rgb(204, 204, 204); font-family: Menlo, Monaco, "Courier New", monospace; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; word-spacing: 0px; white-space: pre; text-decoration-thickness: initial; text-decoration-style: initial; text-decoration-color: initial;"><br /></span></div><div>and this is when the total absolute forces themselves are ~1.55ev/A</div><div><br /></div><div>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?</div><div><br /></div><div>Thank you for any help anyone can provide</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/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/7c11e017-5a3e-4e0a-ad6d-c297781eacb6n%40googlegroups.com</a>.<br />