[CP2K-user] [CP2K:17429] Issues with stress tensor calculations.

[Jürg Hutter]

Hi

as I wrote in my first answer, I strongly advice against using smoothing methods for the density or
the gradients. As far as I know, the stress dependence on the smoothing function is not implemented.
If the grid convergence is a problem, I would use other methods, e.g. NLCC or GAPW.

regards

Juerg Hutter

________________________________________
From: cp2k at googlegroups.com <cp2k at googlegroups.com> on behalf of Dina Kussainova <dinak at princeton.edu>
Sent: Tuesday, August 2, 2022 9:34 PM
To: cp2k
Subject: Re: [CP2K:17429] Issues with stress tensor calculations.

Dear Professor Juerg Hutter,

Thank you for your help. Turning off LONG_RANGE_CORRECTION helped to get similar stress tensors for both analytical and numerical methods. However, I have also noticed that XC_GRID parameters significantly affect stress tensor calculations as well as forces/energy calculations. Please see below snapshots of some test calculations with different parameters. Considering this, what is the best way to identify the correct parameters for my calculations?

Thanks,
Dina

[orig.png]
[test2.png]

[test8.png]

[test11.png]

[test12.png]

[test13.png]


On Monday, July 18, 2022 at 4:49:40 AM UTC-4 jgh wrote:
Hi

unfortunately, without the complete input it is not possible to reproduce your findings.
>From your input I would suggest:

- use a recent version of CP2K (version 9.1 or the recently released 2022.1)
- Don't use LONG_RANGE_CORRECTION in D3
- for debugging purposes it might also be good to switch off CALCULATE_C9_TERM
- Use EPS_SCF 1.E-7 or lower
- Don't use &XC_GRID
XC_DERIV SPLINE2
XC_SMOOTH_RHO NN50
&END XC_GRID
- Increase REL_CUTOFF to 60 or higher
- Use CELL_REF with 20% increased cell parameters to stabilize the finite difference calculation.

best regards

Juerg Hutter

________________________________________
From: cp... at googlegroups.com <cp... at googlegroups.com> on behalf of Dina Kussainova <di... at princeton.edu>
Sent: Friday, July 15, 2022 10:59 PM
To: cp2k
Subject: [CP2K:17324] Issues with stress tensor calculations.

Hi,

I am doing single-point calculations of Li2CO3 with PBE functional, DZVP basis sets, and GTH pseudopotentials. Please see attached input file. However, I am getting different values of stress tensors when I am choosing analytical or numerical methods. When I am doing NPT simulations with the analytical method for stress tensor calculations, looks like volume is also overestimated. Do you know what can be a reason for this? Any help would be appreciated.

Below are snapshots of the output files.

Thanks,
Dina

[analytical.png]

[numerical.png]

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns... at googlegroups.com<mailto:cp2k+uns... at googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/61ebfe8b-437f-42ad-affb-8d3b64a344b1n%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/61ebfe8b-437f-42ad-affb-8d3b64a344b1n%40googlegroups.com?utm_medium=email&utm_source=footer>.

--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com<mailto:cp2k+unsubscribe at googlegroups.com>.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/fddf3eab-9bb6-4b56-b127-9c019eb0bee3n%40googlegroups.com<https://groups.google.com/d/msgid/cp2k/fddf3eab-9bb6-4b56-b127-9c019eb0bee3n%40googlegroups.com?utm_medium=email&utm_source=footer>.

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/ZR0P278MB07592E9D1F70A5BAE1E601FC9F9C9%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test8.png
Type: image/png
Size: 25087 bytes
Desc: test8.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: orig.png
Type: image/png
Size: 25273 bytes
Desc: orig.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0007.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test13.png
Type: image/png
Size: 25983 bytes
Desc: test13.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0008.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test11.png
Type: image/png
Size: 24908 bytes
Desc: test11.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0009.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test12.png
Type: image/png
Size: 25111 bytes
Desc: test12.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0010.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test2.png
Type: image/png
Size: 24543 bytes
Desc: test2.png
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220803/873ad1c6/attachment-0011.png>