Dear community,<div><br /></div><div>I'm interested in measuring the electric field produced by a charge distribution at a given reference coordinate. To that end, I'm printing the EFIELD_CUBE files (after a single-point energy calculation with xTB within the QM/MM framework) and then I use a script to interpolate the value of the property in the given reference point. However, I have some questions:<br />- The grid spacing that results from STRIDE 1 might be too big to make an accurate interpolation, thus I'm wondering if it makes sense to use the &MGRID/CUTOFF keyword within the xTB implementation to change the number of grid points (I get the following warning in CP2K 8.2 "pool cannot reuse pw of another grid")?</div><div>- The atomic coordinates of the cube file are reordered and have been translated relative to the input coordinates, which makes it hard to extract the value of the property at the reference point. Is there any known way to avoid or correct this?</div><div>- I am unsure what are the units used in the printing of the cube files, altough I suspect it to be atomic units. Is this correct? </div><div><br /></div><div>Best regards,<br />Alexandre</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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/31bec450-8c32-46ad-a15f-d9a380f0c2a8n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/31bec450-8c32-46ad-a15f-d9a380f0c2a8n%40googlegroups.com</a>.<br />