<div dir="ltr"><div>Hi Noam,</div><div><br></div>I'd start with using the PBE functional, as it is numerically a lot more stable than BLYP. I don't have much experience with BLYP/smoothing ...<div><br></div><div>Taking a water molecule in vacuum with PBE, forces converge smoothly with cutoff  - results below are for the input attached and cutoffs of 400,600,800 and 1000 Ry. I've not explicitly checked the translational/rotational invariance, but the sum of forces is getting very reasonable, which probably is a good hint of things being OK.</div><div><br></div><div>HTH,</div><div><br></div><div>Matt</div><div><br></div><div><div>mattw@mattw:~$ grep -A 6 ' ATOMIC FORCES in ' H2O-geoopt_???.out H2O-geoopt_????.out</div></div><div>H2O-geoopt_400.out: ATOMIC FORCES in [a.u.]<br></div><div><div>H2O-geoopt_400.out-</div><div>H2O-geoopt_400.out- # Atom   Kind   Element          X              Y              Z</div><div>H2O-geoopt_400.out-      1      1      O           0.00000009     0.00000003    -0.01711982</div><div>H2O-geoopt_400.out-      2      2      H          -0.00000004    -0.01005315     0.00720393</div><div>H2O-geoopt_400.out-      3      2      H          -0.00000004     0.01005312     0.00720394</div><div>H2O-geoopt_400.out- SUM OF ATOMIC FORCES           0.00000001    -0.00000000    -0.00271194     0.00271194</div><div>--</div><div>H2O-geoopt_600.out: ATOMIC FORCES in [a.u.]</div><div>H2O-geoopt_600.out-</div><div>H2O-geoopt_600.out- # Atom   Kind   Element          X              Y              Z</div><div>H2O-geoopt_600.out-      1      1      O           0.00000008     0.00000002    -0.01357588</div><div>H2O-geoopt_600.out-      2      2      H          -0.00000004    -0.01004780     0.00719010</div><div>H2O-geoopt_600.out-      3      2      H          -0.00000004     0.01004777     0.00719010</div><div>H2O-geoopt_600.out- SUM OF ATOMIC FORCES           0.00000001    -0.00000000     0.00080432     0.00080432</div><div>--</div><div>H2O-geoopt_800.out: ATOMIC FORCES in [a.u.]</div><div>H2O-geoopt_800.out-</div><div>H2O-geoopt_800.out- # Atom   Kind   Element          X              Y              Z</div><div>H2O-geoopt_800.out-      1      1      O           0.00000006     0.00000002    -0.01405749</div><div>H2O-geoopt_800.out-      2      2      H          -0.00000003    -0.01004908     0.00718143</div><div>H2O-geoopt_800.out-      3      2      H          -0.00000003     0.01004906     0.00718144</div><div>H2O-geoopt_800.out- SUM OF ATOMIC FORCES           0.00000000    -0.00000000     0.00030538     0.00030538</div><div>--</div><div>H2O-geoopt_1000.out: ATOMIC FORCES in [a.u.]</div><div>H2O-geoopt_1000.out-</div><div>H2O-geoopt_1000.out- # Atom   Kind   Element          X              Y              Z</div><div>H2O-geoopt_1000.out-      1      1      O           0.00000006     0.00000002    -0.01445572</div><div>H2O-geoopt_1000.out-      2      2      H          -0.00000003    -0.01004982     0.00718060</div><div>H2O-geoopt_1000.out-      3      2      H          -0.00000003     0.01004980     0.00718060</div><div>H2O-geoopt_1000.out- SUM OF ATOMIC FORCES           0.00000000    -0.00000000    -0.00009452     0.00009452</div></div><div><br></div><div><br></div><div><br></div></div>