<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8"> <meta name="Generator" content="Microsoft Word 15 (filtered medium)"> <style></style> </head> <body lang="EN-US" link="blue" vlink="purple"> <div class="WordSection1"> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Dear Dev<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">an EPS_SCF threshold value of 0.001 for the SCF convergence is by far too large which results in unconverged energies and forces screwing up your system. Likewise, I suggest a smaller EPS_DEFAULT value of at least 1.0E-12 or less. It is also highly recommended to use an OUTER_SCF with OT.<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">HTH<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Matthias<o:p></o:p></span></p> <p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p> <p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> cp...@googlegroups.com <...@googlegroups.com> <b>On Behalf Of </b>Dev Rana<br> <b>Sent:</b> Montag, 11. Mai 2020 17:25<br> <b>To:</b> cp2k <...@googlegroups.com><br> <b>Subject:</b> [CP2K:13280] SCCS Solvation Model Implementation - Molecule Radii Explodes<o:p></o:p></span></p> <p class="MsoNormal"><o:p> </o:p></p> <div> <p class="MsoNormal">Hello,<o:p></o:p></p> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">I'm trying to simulate a single carbon atom within an aluminum block using a solvation model for the aluminum. I can do this in a normal MD simulation, however to reduce computation time/cost I'd like to use a solvation model to reduce system size. Using the QS regtest-3 files as a starting point, I've created my own input file for my, now 5 atom (Al4C) system. I am finding that the bond lengths between my carbon and aluminum atoms go from something reasonable to over 5 A, and that GEO_OPT is still not complete. Any suggestions would be appreciated. I've attached the input, output, and trajectory files for your review. I've also paste the input below.<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">Best Regards,<o:p></o:p></p> </div> <div> <p class="MsoNormal">Dev<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <div> <p class="MsoNormal">&GLOBAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> PROJECT SolvationAl4C<o:p></o:p></p> </div> <div> <p class="MsoNormal"> RUN_TYPE geo_opt<o:p></o:p></p> </div> <div> <p class="MsoNormal"> PRINT_LEVEL medium<o:p></o:p></p> </div> <div> <p class="MsoNormal">&END GLOBAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">&FORCE_EVAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> METHOD Quickstep<o:p></o:p></p> </div> <div> <p class="MsoNormal"> # STRESS_TENSOR analytical<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &DFT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> BASIS_SET_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/BASIS_MOLOPT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> POTENTIAL_FILE_NAME /shared/centos7/cp2k/cp2k-6.1.0/data/POTENTIAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &QS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> METHOD GPW<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_DEFAULT 1.0E-10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EXTRAPOLATION ASPC<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END QS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &MGRID<o:p></o:p></p> </div> <div> <p class="MsoNormal"> NGRIDS 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> CUTOFF [Ry] 600<o:p></o:p></p> </div> <div> <p class="MsoNormal"> REL_CUTOFF [Ry] 250<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &RS_GRID<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTRIBUTION_TYPE distributed<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END RS_GRID<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END MGRID<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal"> &SCCS ON<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ALPHA [N*m^-1] 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal"> BETA [GPa] 1.3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DELTA_RHO 2.0E-5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DERIVATIVE_METHOD fft<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DIELECTRIC_CONSTANT 80<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_SCCS 1.0E-3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GAMMA [mN/m] 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_SCF 0.0<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MAX_ITER 100<o:p></o:p></p> </div> <div> <p class="MsoNormal">! METHOD Andreussi<o:p></o:p></p> </div> <div> <p class="MsoNormal"> METHOD Fattebert-Gygi<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MIXING 0.6<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &ANDREUSSI<o:p></o:p></p> </div> <div> <p class="MsoNormal"> RHO_MAX 0.001<o:p></o:p></p> </div> <div> <p class="MsoNormal"> RHO_MIN 0.0001<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END ANDREUSSI<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &FATTEBERT-GYGI<o:p></o:p></p> </div> <div> <p class="MsoNormal"> BETA 1.3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> RHO_ZERO 0.00078<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END FATTEBERT-GYGI<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal"> SCF_GUESS ATOMIC<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MAX_SCF 300<o:p></o:p></p> </div> <div> <p class="MsoNormal"> EPS_SCF 1.0E-3<o:p></o:p></p> </div> <div> <p class="MsoNormal">! ADDED_MOS 200<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &DIAGONALIZATION<o:p></o:p></p> </div> <div> <p class="MsoNormal">! ALGORITHM STANDARD<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &END DIAGONALIZATION<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &SMEAR TRUE<o:p></o:p></p> </div> <div> <p class="MsoNormal">! METHOD FERMI_DIRAC<o:p></o:p></p> </div> <div> <p class="MsoNormal">! ELECTRONIC_TEMPERATURE 300<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &END SMEAR<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &MIXING<o:p></o:p></p> </div> <div> <p class="MsoNormal">! METHOD BROYDEN_MIXING<o:p></o:p></p> </div> <div> <p class="MsoNormal">! ALPHA 0.1<o:p></o:p></p> </div> <div> <p class="MsoNormal">! BETA 1.0<o:p></o:p></p> </div> <div> <p class="MsoNormal">! NBROYDEN 8<o:p></o:p></p> </div> <div> <p class="MsoNormal">! &END MIXING<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &OT on<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MINIMIZER DIIS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> PRECONDITIONER FULL_SINGLE_INVERSE<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END OT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END SCF<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &XC<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &XC_FUNCTIONAL PBE<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END XC_FUNCTIONAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END XC<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &PDOS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> APPEND TRUE<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END PDOS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &SCCS ON<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &DENSITY_GRADIENT on<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END DENSITY_GRADIENT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &DIELECTRIC_FUNCTION on<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END DIELECTRIC_FUNCTION<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &POLARISATION_POTENTIAL on<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> APPEND<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END POLARISATION_POTENTIAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END SCCS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END DFT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> <o:p></o:p></p> </div> <div> <p class="MsoNormal"> &SUBSYS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ABC [angstrom] 10 10 10<o:p></o:p></p> </div> <div> <p class="MsoNormal"> PERIODIC XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MULTIPLE_UNIT_CELL 1 1 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &TOPOLOGY<o:p></o:p></p> </div> <div> <p class="MsoNormal"> COORD_FILE_NAME Al4C.xyz<o:p></o:p></p> </div> <div> <p class="MsoNormal"> COORD_FILE_FORMAT XYZ<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MULTIPLE_UNIT_CELL 1 1 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &KIND C<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ELEMENT C<o:p></o:p></p> </div> <div> <p class="MsoNormal"> BASIS_SET DZVP-MOLOPT-GTH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> POTENTIAL GTH-PBE-q4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END KIND<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &KIND Al<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ELEMENT Al<o:p></o:p></p> </div> <div> <p class="MsoNormal"> BASIS_SET DZVP-MOLOPT-SR-GTH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> POTENTIAL GTH-PBE-q3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END KIND<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END SUBSYS<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">&END FORCE_EVAL<o:p></o:p></p> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> <div> <p class="MsoNormal">&MOTION<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &GEO_OPT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> OPTIMIZER BFGS<o:p></o:p></p> </div> <div> <p class="MsoNormal"> MAX_ITER 100<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END GEO_OPT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &STRUCTURE_DATA<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 1 2<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 1 3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 1 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 1 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 2 3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 2 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 2 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 3 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 3 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> DISTANCE 4 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 2 3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 2 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 2 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 3 2<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 3 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 3 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 4 2<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 4 3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 4 5<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 5 2<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 5 3<o:p></o:p></p> </div> <div> <p class="MsoNormal"> ANGLE 1 5 4<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END STRUCTURE_DATA<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &TRAJECTORY<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END TRAJECTORY<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &VELOCITIES<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END VELOCITIES<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &FORCES<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END FORCES<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END CELL<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &RESTART<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> GEO_OPT 1<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END EACH<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END RESTART<o:p></o:p></p> </div> <div> <p class="MsoNormal"> &END PRINT<o:p></o:p></p> </div> <div> <p class="MsoNormal">&END MOTION<o:p></o:p></p> </div> </div> <div> <p class="MsoNormal"><o:p> </o:p></p> </div> </div> <p class="MsoNormal">-- <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+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>.<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/b2a95387-2350-4f1c-8eb9-32976ddffbd0%40googlegroups.com?utm_medium=email&utm_source=footer"> https://groups.google.com/d/msgid/cp2k/b2a95387-2350-4f1c-8eb9-32976ddffbd0%40googlegroups.com</a>.<o:p></o:p></p> </div> </body> </html>