<div dir="ltr">Hi All, I have the problem with the convergence when using SCAN+rVV10 functional for the very simple system with only one water molecule. The convergence just fluctuate around 0.002, and can not decrease more (see below). 1 OT DIIS 0.80E-01 54.9 0.00325654 -17.2027773693 -4.08E-07 2 OT DIIS 0.80E-01 27.5 0.00276549 -17.2027717490 5.62E-06 3 OT DIIS 0.80E-01 27.6 0.00376733 -17.2027529642 1.88E-05 4 OT DIIS 0.80E-01 27.5 0.00424108 -17.2027547590 -1.79E-06 5 OT DIIS 0.80E-01 27.6 0.00315550 -17.2027413567 1.34E-05 6 OT DIIS 0.80E-01 27.4 0.00191592 -17.2027409117 4.45E-07 7 OT SD 0.80E-01 27.5 0.00388064 -17.2027478523 -6.94E-06 8 OT DIIS 0.80E-01 27.6 0.00287957 -17.2027087049 3.91E-05 9 OT DIIS 0.80E-01 27.6 0.00156285 -17.2027379497 -2.92E-05 10 OT DIIS 0.80E-01 27.5 0.00129428 -17.2027445179 -6.57E-06 11 OT DIIS 0.80E-01 27.2 0.00161973 -17.2027517827 -7.26E-06 12 OT DIIS 0.80E-01 27.5 0.00220730 -17.2027652153 -1.34E-05 13 OT DIIS 0.80E-01 28.0 0.00127484 -17.2027571303 8.09E-06 14 OT DIIS 0.80E-01 27.5 0.00260157 -17.2027619579 -4.83E-06 15 OT DIIS 0.80E-01 27.8 0.00124040 -17.2027613316 6.26E-07 16 OT DIIS 0.80E-01 27.4 0.00161009 -17.2027627322 -1.40E-06 17 OT DIIS 0.80E-01 27.7 0.00270573 -17.2027633849 -6.53E-07 18 OT DIIS 0.80E-01 27.5 0.00143321 -17.2027609365 2.45E-06 19 OT DIIS 0.80E-01 27.6 0.00949434 -17.2027630484 -2.11E-06 20 OT DIIS 0.80E-01 27.6 0.00206430 -17.2027666772 -3.63E-06 21 OT DIIS 0.80E-01 27.3 0.00327081 -17.2027624981 4.18E-06 22 OT DIIS 0.80E-01 27.6 0.00140816 -17.2027603709 2.13E-06 23 OT DIIS 0.80E-01 27.6 0.00194449 -17.2027617267 -1.36E-06 24 OT DIIS 0.80E-01 27.5 0.00248205 -17.2027603408 1.39E-06 25 OT DIIS 0.80E-01 28.1 0.00317999 -17.2027635590 -3.22E-06 26 OT DIIS 0.80E-01 27.6 0.00376279 -17.2027604966 3.06E-06 27 OT DIIS 0.80E-01 27.4 0.00175908 -17.2027571292 3.37E-06 28 OT DIIS 0.80E-01 27.6 0.00113761 -17.2027575317 -4.02E-07 29 OT DIIS 0.80E-01 27.5 0.00631357 -17.2027592196 -1.69E-06 30 OT DIIS 0.80E-01 27.7 0.00129080 -17.2027564625 2.76E-06 31 OT DIIS 0.80E-01 27.3 0.00296235 -17.2027626012 -6.14E-06 32 OT DIIS 0.80E-01 27.6 0.00267161 -17.2027653970 -2.80E-06 33 OT DIIS 0.80E-01 27.6 0.00196123 -17.2027620008 3.40E-06 34 OT DIIS 0.80E-01 27.6 0.00344508 -17.2027683804 -6.38E-06 I tried many changes that was recommended from different posts in the forum (changing MINIMIZER, PRECONDITIONER, ENERGY_GAP, ...), but nothing works. I only have problem with SCAN+rVV10, I tried other functionals (B97M-rV which is also MGGA+rVV10, revPBE-D3, BLYP-D3, rVV10) and they all work. I don't know what was the problem. Could you please give me some suggestions? Has anyone perform the calculation with SCAN/SCAN+rVV10 successfully, if yes, can you send me one example of the input? I installed cp2k following this post: https://groups.google.com/forum/#!topic/cp2k/yJP-HRqx4Y0 Thanks, Huy Below is my input -------------------------------- &GLOBAL PROJECT WATER RUN_TYPE ENERGY IOLEVEL LOW &END GLOBAL &FORCE_EVAL METHOD Quickstep &DFT BASIS_SET_FILE_NAME GTH_BASIS_SETS POTENTIAL_FILE_NAME GTH_POTENTIALS &MGRID CUTOFF 1200 NGRIDS 5 &END &QS EPS_DEFAULT 1.0E-18 WF_INTERPOLATION ASPC &END QS &SCF SCF_GUESS RESTART MAX_SCF 40 EPS_SCF 1.0E-9 &OT MINIMIZER DIIS PRECONDITIONER FULL_ALL ENERGY_GAP 0.002 &END OT &OUTER_SCF MAX_SCF 40 EPS_SCF 1.0E-9 &END OUTER_SCF &PRINT &RESTART &EACH MD 0 &END &END RESTART &END PRINT &END SCF &POISSON PERIODIC NONE PSOLVER MT &END POISSON &XC &XC_FUNCTIONAL &LIBXC FUNCTIONAL MGGA_X_SCAN &END LIBXC &LIBXC FUNCTIONAL MGGA_C_SCAN &END LIBXC &END XC_FUNCTIONAL &vdW_POTENTIAL DISPERSION_FUNCTIONAL NON_LOCAL &NON_LOCAL TYPE RVV10 PARAMETERS 15.7 0.0093 KERNEL_FILE_NAME rVV10_kernel_table.dat &END NON_LOCAL &END vdW_POTENTIAL &XC_GRID XC_DERIV NN50_SMOOTH &END XC_GRID &END XC &END DFT &SUBSYS &CELL ABC 25.0 25.0 25.0 PERIODIC NONE &END CELL &TOPOLOGY COORD_FILE_NAME file.xyz COORD_FILE_FORMAT XYZ &CENTER_COORDINATES &END &END &KIND H BASIS_SET aug-QZV3P-GTH POTENTIAL GTH-PBE-q1 &END KIND &KIND O BASIS_SET aug-QZV3P-GTH POTENTIAL GTH-PBE-q6 &END KIND &KIND C BASIS_SET aug-QZV3P-GTH POTENTIAL GTH-PBE-q4 &END KIND &END SUBSYS &END FORCE_EVAL &MOTION &GEO_OPT MAX_ITER 2000 RMS_FORCE 1.0E-6 MAX_FORCE 3.0E-6 RMS_DR 1.0E-5 MAX_DR 3.0E-5 &END &END </div>