<div dir="ltr">Hi all,<div><br></div><div>I am trying to extract total electron density from a cluster of water molecules in a QM/MM calculation and I am using 6-31G* all electron basis to extract the total electron density cube of the water molecules in the QM region. As you can see in the attached image, the electron density shows "spikes". The description of the TOTAL_DENSITY keyword does gives the warning, which says that:</div><div><br></div><div><div>"..NOTE: The total density in real space might exhibit unphysical features like spikes due to the finite and thus truncated g vector expansion."<br></div></div><div><br></div><div>My question here is that is there any way by which I can reduce the prevalence of spikes in the total electron density by fine tuning GAPW parameters? I am currently using the following setup:</div><div><br></div><div><div>     &QS</div><div>       EPS_DEFAULT     1.0000000000000000E-10</div><div>       EPSFIT     9.9999999999999995E-07</div><div>       EPSISO     9.9999999999999998E-13</div><div>       EPSRHO0     1.0000000000000000E-08</div><div>       ALPHA0_HARD     1.0000000000000000E+01</div><div>       MAP_CONSISTENT  T</div><div>       LMAXN1  6</div><div>       LMAXN0  4</div><div>       QUADRATURE  GC_LOG</div><div>       EXTRAPOLATION  ASPC</div><div>       EXTRAPOLATION_ORDER  4</div><div>       METHOD  GAPW</div><div>     &END QS</div><div>     &MGRID</div><div>       CUTOFF     2.8000000000000000E+02</div><div>       COMMENSURATE  T</div><div>     &END MGRID</div><div>     &XC</div><div>       DENSITY_CUTOFF     1.0000000000000000E-10</div><div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div><div>       TAU_CUTOFF     1.0000000000000000E-10</div><div>       &XC_GRID</div><div>         XC_SMOOTH_RHO  NN50</div><div>         XC_DERIV  SPLINE2_SMOOTH</div><div>       &END XC_GRID</div><div>       &XC_FUNCTIONAL  NO_SHORTCUT</div><div>         &PBE  T</div><div>           PARAMETRIZATION  REVPBE</div><div>         &END PBE</div><div>       &END XC_FUNCTIONAL</div><div>       &VDW_POTENTIAL</div><div>         POTENTIAL_TYPE  PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>           R_CUTOFF     1.2500000000000009E+01</div><div>           TYPE  DFTD3</div><div>           PARAMETER_FILE_NAME ./dftd3.dat</div><div>           REFERENCE_FUNCTIONAL revPBE</div><div>           CALCULATE_C9_TERM  T</div><div>           REFERENCE_C9_TERM  T</div><div>           LONG_RANGE_CORRECTION  F</div><div>         &END PAIR_POTENTIAL</div><div>       &END VDW_POTENTIAL</div><div>     &END XC</div></div><div><br></div><div>.....</div><div><br></div><div><br></div><div><div>     &KIND H</div><div>       BASIS_SET 6-31G*</div><div>       POTENTIAL ALL</div><div>       LEBEDEV_GRID  80</div><div>       RADIAL_GRID  200</div><div>     &END KIND</div><div>     &KIND O</div><div>       BASIS_SET 6-31G*</div><div>       POTENTIAL ALL</div><div>       LEBEDEV_GRID  80</div><div>       RADIAL_GRID  200</div><div>     &END KIND<br></div></div><div><br></div><div>Also if there is else I can do either remove the "spikes" fully from the electron density cube, that would be really helpful as well.</div><div><br></div><div>Best Regards,</div><div>Mayank Dodia</div><div><br></div></div>